POJ3070 矩阵快速幂模板

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

矩阵快速幂模板。mod写到乘法的定义部分就行了。

别忘了 I ( ) 和 i n i t ( ) 要传引用！

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,p=;
struct Matrix{
    int a[][];
    Matrix operator *(const Matrix &b)const
    {
        Matrix c;
        c.a[][]=c.a[][]=c.a[][]=c.a[][]=;
        for(int i=;i<=;i++)
            for(int j=;j<=;j++)
                for(int k=;k<=;k++)
                    (c.a[i][j]+=a[i][k]*b.a[k][j])%=p;
//                        printf("i=%d j=%d c=%d\n",i,j,b.a[i][j]);
        return c;
    }
};
Matrix r,u;
void I(Matrix &k)
{
    k.a[][]=k.a[][]=;
    k.a[][]=k.a[][]=;
}
void init(Matrix &k)
{
    k.a[][]=k.a[][]=k.a[][]=;
    k.a[][]=;
}
int main()
{
    while()
    {
        scanf("%d",&n);if(n==-)return ;
        if(!n){printf("0\n");continue;}
        n--;I(r);init(u);
        while(n)
        {
            if(n&)r=r*u;
            n>>=;u=u*u;
        }
        printf("%d\n",(r.a[][]+r.a[][])%p);
    }
}