POJ 3070 矩阵快速幂

//{//1,1////1,0//}/////这个矩阵自乘n次。连续自乘n次的话就没意思了，那还不如直接上Fibonacci递推公式呢。///矩阵的魅力就在于它可以上快速幂。因为矩阵乘法满足结合律么……///注意取模#include <cstdlib>#include <cstring>#include <cstdio>#include <iostream>using namespace std;int N;struct matrix///定义数组{       int a[2][2];}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(){     origin.a[0][0]=1;     origin.a[0][1]=1;     origin.a[1][0]=1;     origin.a[1][1]=0;     memset(res.a,0,sizeof(res.a));     res.a[0][0]=res.a[1][1]=1;                  //将res.a初始化为单位矩阵}void calc(int n){     while(n)///把n分解成二进制     {             if(n&1)                    res=multiply(res,origin);             n>>=1;             origin=multiply(origin,origin);///origin， A的平方，A的四次方，A的8次方。     }    printf("%d\n",res.a[0][1]%10000);}///话说这么简单的模板真的可以么？？？= =int main(){    while(scanf("%d",&N)!=EOF)    {        if(N==-1) break;            init();            calc(N);    }    return 0;}