本文主要是介绍HDU 1005(矩阵乘法)，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
A number sequence is defined as follows:
f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7.
Given A, B, and n, you are to calculate the value of f(n).
题意：给出一个递推公式，求第N项。
题解：开始试图找出一个能够直接算的递推公式，没找到，所以考虑矩阵乘法。根据给出的递推公式可以设想：
可以使用找规律：
#include <iostream>
#include <cstring>
using namespace std;
int f[1000001];
int main()
{
int a,b,n;
while(cin>>a>>b>>n)
{
if(a==0&&b==0&&n==0)
break;
f[1]=1;
f[2]=1;
int xh=0;
int i;
for(i=3;i<=1001;i++)
{
f[i]=(a*f[i-1]+b*f[i-2])%7;
if(f[i]==1&&f[i-1]==1)
{
break;
}
}
xh=i-2;
n=n%xh;
if(n==0)n=xh;
cout<<f[n]<<endl;
}
return 0;
}
使用矩阵相乘：
#include <iostream>
#include <cstdio>
using namespace std;
const int N=10001;
const __int64 INF = 10e18;
typedef struct m
{
int a[2][2];
}Mat;
typedef struct o
{
int sno; //器官序号
int day;
}Cat;
Cat cat[N]; //
Mat E;//unit matrix;
void InitE(Mat& e)
{
e.a[0][0]=1;
e.a[0][1]=0;
e.a[1][0]=0;
e.a[1][1]=1;
}
Mat MatMul(Mat x,Mat y)
{
int i,j,k;
Mat c;
for(i=0;i<2;++i)
for(j=0;j<2;++j)
{
c.a[i][j]=0;
for(k=0;k<2;++k)
{
c.a[i][j]+=x.a[i][k] * y.a[k][j];
}
c.a[i][j] %=7;
}
return c;
}
Mat MatPow(Mat d,int n) //matrix binary divide,calculate n power of matrix a
{
Mat t=E;
while(n>0)
{
if(n&1)
t=MatMul(t,d);
d = MatMul(d,d);
n>>=1;
}
return t;
}
int main()
{
int p,q,n;
while(cin>>p>>q>>n && p|q|n)
{
if(n<3)
cout<<1<<endl;
else
{
InitE(E); //init unit matrix
Mat m;
m.a[0][0] = p;
m.a[0][1] = q;
m.a[1][0] = 1;
m.a[1][1] = 0;
m = MatPow(m,n-2);
cout<<(m.a[0][0]+m.a[0][1])%7<<endl;
}
}
return 0;
}
这篇关于HDU 1005(矩阵乘法)的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
