poj 3070 Fibonacci (矩阵快速幂乘/模板)

题意:给你一个n,输出Fibonacci (n)%10000的结果

思路:裸矩阵快速幂乘,直接套模板

代码:

#include <cstdio>
#include <cstring>
#include <iostream> using namespace std; typedef long long ll;
const int N=,M=,P=;
const int MOD=;
struct Matrix
{
ll m[N][N];
}; Matrix A={,,
,}; Matrix I={,,
,}; Matrix multi(Matrix a,Matrix b)
{
Matrix ans;
for(int i=;i<N;i++)
{
for(int j=;j<M;j++)
{
ans.m[i][j]=;
for(int k=;k<P;k++)
{
ans.m[i][j]+=a.m[i][k]*b.m[k][j]%MOD;
}
ans.m[i][j]%=MOD;
}
}
return ans;
} Matrix power(Matrix a,int k)
{
Matrix ans=I,p=a;
while(k)
{
if(k&)
{
ans=multi(ans,p);
}
k>>=;
p=multi(p,p);
}
return ans;
} int main(int argc, char const *argv[])
{
int n;
while(scanf("%d",&n)!=-)
{
if(n==-) break;
if(n==)
{
cout<<""<<endl;
continue;
}
Matrix ans=power(A,n-);
printf("%lld\n",ans.m[][] );
}
return ;
}
上一篇:IDC:IDC 清单


下一篇:Java学习之J2EE