矩阵快速幂

矩阵快速幂，原理与整数的快速幂一模（mu）一样，只是在运算的时候是对矩阵进行运算

Code :

矩阵快速幂
能够进行幂运算的矩阵一定都是行与列元素个数相等的矩阵
#include<bits/stdc++.h>
using namespace std;
int n;

struct Matrix
{
    int mat[20][20];
};

Matrix Matrix_Mul(Matrix m1, Matrix m2)
{
    Matrix mm;
    memset(mm.mat, 0, sizeof(mm.mat));
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
            for(int k = 0; k < n; k++)
                mm.mat[i][j] += (m1.mat[i][k] * m2.mat[k][j]);
    return mm;
}


Matrix Matrix_Pow(Matrix m, int k)
{
    Matrix ans; // int ans = 1;
    for(int i = 0; i < n; i++)
    {
        memset(ans.mat[i], 0, sizeof(ans.mat[i]));
        ans.mat[i][i] = 1;
    }

    while(k != 0)
    {
        if(k & 1)
            ans = Matrix_Mul(ans, m); //ans *= m;
        m = Matrix_Mul(m, m); // m*= m;
        k >>= 1;
    }
    return ans;
}


int main()
{
    Matrix mm;
    scanf("%d", &n); //n行 能进行快速幂运算的都是方阵；
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
            scanf("%d", &mm.mat[i][j]);
    int k;
    scanf("%d", &k);
    mm = Matrix_Pow(mm, k);
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
            cout << mm.mat[i][j] << " ";
        cout << endl;
    }
    return 0;
}

这回就美滋滋的开始进行利用吧！！