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
",res.a[0][1]%10000);
}

int main()
{
    while(scanf("%d",&N)&&N!=-1)
    {
            init();
            calc(N);
    }
    return 0;
}