POJ 3070 Fibonacci

1、本题用到矩阵快速幂，注意矩阵乘法不满足交换律。

2、将矩阵写成类会更易编程。一次AC~

#include<cstdio>
#include<cstring>
using namespace std;
const int MAXN=2;
const int MAXM=2;
struct Matrix{
    int n,m;
 int a[MAXN][MAXM];
 void clear(){
  n=m=0;
  memset(a,0,sizeof(a));
 }
 Matrix operator + (const Matrix &b) const{
     Matrix temp;
  temp.n=n;temp.m=m;
  for(int i=0;i<n;i++)
   for(int j=0;j<m;j++)
    temp.a[i][j]=a[i][j]+b.a[i][j];
  return temp;
 }
 Matrix operator - (const Matrix &b) const{
     Matrix temp;
  temp.n=n;temp.m=m;
  for(int i=0;i<n;i++)
   for(int j=0;j<m;j++)
    temp.a[i][j]=a[i][j]-b.a[i][j];
  return temp;
 }
 Matrix operator * (const Matrix &b) const{
     Matrix temp;
  temp.clear();
  temp.n=n;temp.m=b.m;
  for(int i=0;i<n;i++)
   for(int j=0;j<b.m;j++)
    for(int k=0;k<m;k++)
     temp.a[i][j]+=a[i][k]*b.a[k][j];
  return temp;
 }
 Matrix operator % (const int mod) const{
     Matrix temp;
  temp.n=n;temp.m=m;
  for(int i=0;i<n;i++)
   for(int j=0;j<m;j++)
    temp.a[i][j]=a[i][j]%mod;
  return temp;
 }
};
Matrix pow_mod(Matrix A,int i,int mod){
    if(i==0){
    Matrix u;
    u.clear();
    u.n=MAXN;u.m=MAXM;
    u.a[0][0]=u.a[1][1]=1%mod;
    return u;
 }
 Matrix temp=pow_mod(A,i>>1,mod);
 temp=temp*temp%mod;
 if(i&1) temp=temp*A%mod;
 return temp;
}
Matrix A,An,X;
int main(){
 int n;
    A.n=MAXN;A.m=MAXM;
 A.a[0][0]=0;
 A.a[0][1]=A.a[1][0]=A.a[1][1]=1;
 while(scanf("%d",&n)==1){
     if(n==-1) break;
        An=pow_mod(A,n,10000);
        Matrix X;
  X.n=MAXN;X.m=1;
  X.a[0][0]=0;X.a[1][0]=1;
  X=An*X;
  printf("%d\n",X.a[0][0]);
 }
 return 0;
}