poj 3070 Fibonacci

http://poj.org/problem?id=3070

矩阵的快速幂，二分

#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 10000
using namespace std;
const int mod=10000;

int n;
struct node
{
    int a[4][4];
};

node multi(node s,node t)
{
    node c;
    for(int i=0; i<2; i++)
    {
        for(int j=0; j<2; j++)
        {
            c.a[i][j]=0;
            for(int k=0; k<2; k++)
            {
                c.a[i][j]=(c.a[i][j]+s.a[i][k]*t.a[k][j])%mod;
            }
        }
    }
    return c;
}


int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        if(n==-1) break;
        node m,x;
        m.a[0][0]=1; m.a[0][1]=0;
        m.a[1][0]=0; m.a[1][1]=1;
        x.a[0][0]=x.a[0][1]=x.a[1][0]=1;
        x.a[1][1]=0;;
        while(n)
        {
            if(n&1)
            {
                m=multi(m,x);
            }
            x=multi(x,x);
            n>>=1;
        }
        printf("%d
",m.a[0][1]);
    }
    return 0;
}