POJ-3070 Fibonacci 矩阵乘法
题目链接：http://poj.org/problem?id=3070
矩阵乘法优化，水题一枚。
1 //STATUS:C++_AC_132MS_0KB 2 #include<stdio.h> 3 #include<stdlib.h> 4 #include<string.h> 5 #include<math.h> 6 #include<iostream> 7 #include<string> 8 #include<algorithm> 9 #include<vector> 10 #include<queue> 11 #include<stack> 12 #include<map> 13 using namespace std; 14 #define LL __int64 15 #define pii pair<int,int> 16 #define Max(a,b) ((a)>(b)?(a):(b)) 17 #define Min(a,b) ((a)<(b)?(a):(b)) 18 #define mem(a,b) memset(a,b,sizeof(a)) 19 #define lson l,mid,rt<<1 20 #define rson mid+1,r,rt<<1|1 21 const int N=25,INF=0x3f3f3f3f,MOD=10000,STA=8000010; 22 const double DNF=1e13; 23 24 struct Matrix{ 25 LL ma[2][2]; 26 Matrix friend operator * (const Matrix& a,const Matrix& b){ 27 Matrix ret; 28 mem(ret.ma,0); 29 int i,j,k; 30 for(i=0;i<2;i++){ 31 for(j=0;j<2;j++){ 32 for(k=0;k<2;k++) 33 ret.ma[i][j]=(ret.ma[i][j]+a.ma[i][k]*b.ma[k][j])%MOD; 34 } 35 } 36 return ret; 37 } 38 }a,f; 39 40 int n; 41 42 void mutilpow(int k) 43 { 44 int i,j; 45 mem(f.ma,0); 46 f.ma[0][0]=f.ma[1][1]=1; 47 for(;k;k>>=1){ 48 if(k&1)f=f*a; 49 a=a*a; 50 } 51 } 52 53 int main() 54 { 55 // freopen("in.txt","r",stdin); 56 int i,j; 57 LL ans; 58 while(~scanf("%d",&n) && n!=-1) 59 { 60 a.ma[0][0]=0;a.ma[0][1]=1; 61 a.ma[1][0]=a.ma[1][1]=1; 62 if(n>=2){ 63 mutilpow(n-2); 64 ans=f.ma[1][0]+f.ma[1][1]; 65 } 66 else ans=n; 67 68 printf("%I64d\n",ans%MOD); 69 } 70 return 0; 71 }
