矩阵

矩阵的乘法：

Matrix：

乘法: A*B ，需要A的列数=B的行数

矩阵 C = A*B

则： C[i,j] = ∑A[i,k]*B[k,j]

 矩阵的乘法代码：

//矩阵的乘法 

void mul(int A[][MAXN],int B[][MAXN])//A=A*B
{
        int R[MAXN][MAXN];
        memset(R,0,sizeof(R));
        for(int k=0;k<N;k++)
            for(int i=0;i<N;i++)if(A[i][k])
                for(int j=0;j<N;j++)
                    R[i][j]+=A[i][k]*B[k][j];
        /*
            如果取余可能能加快计算 
            int R[MAXN][MAXN];
            for(int i=0;i<N;i++)
                for(int j=0;j<N;j++)
                {
                    long long sum = 0;
                    for(int k=0;k<N;k++)
                        sum+=A[i][k]*B[k][j];
                    R[i][j]=sum%MOD;
                }
        */
        for(int i=0;i<N;i++)
            for(int j=0;j<N;j++)
                A[i][j]=R[i][j];
}

矩阵快速幂模板：

#include<iostream>
#include<cstring>
using namespace std;

int n, k, M;
struct mat{
    int ma[65][65];
    mat() {memset(ma, 0, sizeof(ma));}//¹¹ÔìÁã¾ØÕó 
    mat(int a)  //¹¹Ô쵥λ¾ØÕó 
    {
        memset(ma, 0, sizeof(ma));
        for(int i=0; i<n; i++)
        ma[i][i] = 1;
    }
    void operator *= (const mat &x)//¾ØÕó³Ë·¨ 
    {
        mat tmp;
        int i, j, k;
        for( k=0; k<n; k++)
        for( i=0; i<n; i++)
        if(ma[i][k])
        {   long long tot=0; 
            for( j=0; j<n; j++)
                tot += ma[i][k]*x.ma[k][j];
                ma[i][j] = tot%M;
        }
        for(i=0; i<n; i++)
        for(j=0; j<n; j++)
        ma[i][j] = tmp.ma[i][j];
    }
};
mat pow_mod(mat a, int b)//¾ØÕó¿ìËÙÃÝ 
{
    mat tmp = mat(1);
    while(b>0)
    {
        if(b&1) tmp*=a;
        b>>=1;
        a*=a;
    }
    return tmp;
}

int main()
{
    while(scanf("%d%d%d", &n, &k, &M)==3)
    {
        mat A, B;
        for(int i=0; i<n; i++)
        for(int j=0; j<n; j++)
        scanf("%d", &A.ma[i][j]);
        printf("

");
        mat ans = pow_mod(A, k+1);
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            printf("%8d", ans.ma[i][j]);
            printf("
");
        }
    }
    return 0;
}