POJ 3070 矩陣快速乘冪

/*很水的矩陣題，沒什麼說的... */#include<iostream>#include<cstdio>#include<cmath>using namespace std;struct matrix{     long long a[3][3];}A,B;matrix mult2(matrix A,matrix B){  /// 矩陣乘法     matrix C;    for(int i=1;i<=2;i++){   /// 第一個矩陣的行         for(int j=1;j<=2;j++){    /// 第二個矩陣的列             C.a[i][j]=0;            for(int k=1;k<=2;k++){  /// 第二個矩陣的行 或者 第一個矩陣的列                 C.a[i][j] = (C.a[i][j] + A.a[i][k]*B.a[k][j])%10000;            }        }    }    return C;}matrix mult(matrix A,long long n){  ///  二分快速冪     if(n==1) return A;    A = mult(A,n/2);    if(n%2==0)  return  A = mult2(A,A);    return  A = mult2(mult2(A,A),B);}int main(){    long long  n;    while(scanf("%lld",&n)!=EOF){        if(!n) { printf("0\n"); continue; }        if(n == -1) break;        A.a[1][1]=1; A.a[1][2]=1;        A.a[2][1]=1; A.a[2][2]=0;        B = A;        A = mult(A,n);         printf("%lld\n",A.a[1][2]%10000);    }}