POJ 3070 Fibonacci（矩阵快速幂）

题目：http://poj.org/problem?id=3070

代码：

#include<stdio.h>
#include<string.h>

using namespace std;

const int MOD=10000;

struct point
{
    int m[2][2];
}ans,base;

point multi(point a,point b)
{
    point temp;
    for(int i=0;i<2;i++)
    {
        for(int j=0;j<2;j++)
        {
            temp.m[i][j]=0;
            for(int k=0;k<2;k++)
            {
                temp.m[i][j]=(temp.m[i][j]+a.m[i][k]*b.m[k][j])%MOD;
            }
        }
    }
    return temp;
}

int fast(int n)
{
    base.m[0][0]=base.m[0][1]=base.m[1][0]=1;
    base.m[1][1]=0;
    ans.m[0][0]=ans.m[1][1]=1;
    ans.m[0][1]=ans.m[1][0]=0;
    while(n)
    {
        if(n%2==1)
        {
            ans=multi(ans,base);
        }
        base=multi(base,base);
        n=n/2;
    }
    return ans.m[0][1];
}
int main()
{
    int n;
    while(~scanf("%d",&n)&&n!=-1)
    {
        printf("%d\n",fast(n));
    }
}

分析：

大致理解了，打不出来。55555。