poj3070 Fibonacci 矩阵快速幂

学了线代之后 终于明白了矩阵的乘法。。

于是 第一道矩阵快速幂。。

实在是太水了。。。

这差不多是个模板了

#include <cstdlib>#include <cstring>#include <cstdio>#include <iostream>using namespace std;int N;struct matrix{       int a[3][3];}origin,res;matrix multiply(matrix x,matrix y){       matrix temp;       memset(temp.a,0,sizeof(temp.a));       for(int i=0;i<2;i++)       {               for(int j=0;j<2;j++)               {                       for(int k=0;k<2;k++)                       {                               temp.a[i][j]+=(x.a[i][k]*y.a[k][j]%10000);                               temp.a[i][j]%=10000;                       }               }       }       return temp;}void init(){     memset(res.a,0,sizeof(res.a));     res.a[0][0]=res.a[1][1]=1;     res.a[1][0]=res.a[0][1]=0;     origin.a[0][0]=origin.a[0][1]=origin.a[1][0]=1;     origin.a[1][1]=0;}void calc(int n){     while(n)     {             if(n&1)                    res=multiply(res,origin);             n>>=1;             origin=multiply(origin,origin);     }     printf("%d\n",res.a[0][1]%10000);}int main(){    while(scanf("%d",&N)&&N!=-1)    {            init();            calc(N);    }    return 0;}