{矩阵相关}

基础：http://blog.csdn.net/q3498233/article/details/5786180

题集：

        http://www.cppblog.com/notonlysuccess/archive/2009/03/03/75405.aspx

        http://acm.hdu.edu.cn/forum/read.php?tid=15908

待补题。。。

Matrix Power Series

 POJ - 3233

题意：求矩阵和S = A + A² + A³ + … + A^k. 

构造矩阵

A A

0 1

然后这个矩阵自乘K次即可

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define CLR(m,a) memset(m,a,sizeof(m))
const int maxn=65;
int n,k,mod;
struct Matrix
{
    int n,m[maxn][maxn];
    void init(int sz)
    {
        CLR(m,0);
        n=sz;
    }
    Matrix(int sz=0){init(sz);}
    void setE()
    {
        for(int i=0;i<n;i++) m[i][i]=1;
    }
    Matrix operator* (const Matrix &a)
    {
        Matrix ans(n);
        for(int k=0;k<n;k++)
            for(int i=0;i<n;i++)
                for(int j=0;j<n;j++)
                    ans.m[i][j]=(ans.m[i][j]+(m[i][k]*a.m[k][j])%mod)%mod;
        return ans;
    }
};
int main()
{
    while(scanf("%d%d%d",&n,&k,&mod)!=EOF){
        int tot=n*2;
        Matrix a(tot),b(tot);
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                scanf("%d",&a.m[i][j]);
        for(int i=0;i<n;i++)
            for(int j=n;j<tot;j++)
                a.m[i][j]=a.m[i][j-n];
        for(int i=n;i<tot;i++) a.m[i][i]=1;
        b.setE();
        while(k){
            if(k&1) b=b*a;
            k>>=1;
            a=a*a;
        }
        for(int i=0;i<n;i++)
            for(int j=n;j<tot;j++)
                printf("%d%c",b.m[i][j],j==tot-1?'
':' ');
    }
}