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;

}