矩阵快速幂 POJ 3070
题意:算出Fn,矩阵.,就是求矩阵的n次方,先定义矩阵乘法,再快速幂就解决了。
#include<stdio.h>
#include<cstring>
#include<iostream>
#define ll long long
using namespace std;
struct matix
{
int a[2][2];
matix()//构造函数,初始化
{
for(int i=0;i<2;i++)
{
for(int j=0;j<2;j++)
{
a[i][j]=0;
}
}
}
};
matix operator *(matix a,matix b)//重载乘号
{
matix t;
for(int i=0;i<2;i++)
{
for(int j=0;j<2;j++)
{
for(int k=0;k<2;k++)
{
t.a[i][j]+=a.a[i][k]*b.a[k][j]%10000;
}
}
}
return t;
}
matix quick(matix a,int b)
{
matix ret;
ret.a[0][0]=1;
ret.a[1][1]=1;
while (b!=0)
{
i f(b&1)
{
ret=ret*a;
}
a=a*a;
b>>=1;
}
return ret;
}
int main ()
{
int t;
while (scanf("%d",&t)!=EOF)
{
if(t==-1)
{
break;
}
matix k;
k.a[0][0]=k.a[0][1]=k.a[1][0]=1;
matix temp=quick(k,t);
cout<<temp.a[1][0]%10000<<endl;
}
return 0;
}