poj3070 求斐波那契数列第n项 ——矩阵快速幂

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

用矩阵快速幂加速递推。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,p=;
struct Matrix{
    int a[][];
    Matrix operator * (const Matrix &y) const
    {
        Matrix x;
        memset(x.a,,sizeof x.a);
        for(int i=;i<=;i++)
            for(int k=;k<=;k++)
                for(int j=;j<=;j++)
                    x.a[i][j]+=y.a[i][k]*a[k][j];
        return x;
    }
}m,ans;
void mod()
{
    m.a[][]%=p;m.a[][]%=p;//
    m.a[][]%=p;m.a[][]%=p;
    ans.a[][]%=p;ans.a[][]%=p;
}
int main()
{
    while(scanf("%d",&n)==)
    {
        if(n==-)return ;
        memset(m.a,,sizeof m.a);
        memset(ans.a,,sizeof ans.a);
        m.a[][]=;m.a[][]=;
        m.a[][]=;m.a[][]=;
        ans.a[][]=;ans.a[][]=;
        if(!n){printf("0\n");continue;}
        if(n==){printf("1\n");continue;}
        if(n==){printf("1\n");continue;}
        n--;//!
        while(n)
        {
            if(n&)ans=ans*m;
            m=m*m;
            mod();
            n=(n>>);
        }
        printf("%d\n",ans.a[][]%p);
    }
    return ;
}