矩阵幂

不是快速幂，只是想写离散作业，随便写的一个幂

#include<iostream>
using namespace std;
int main()
{
    int a[100][100];
    int b[100][100];
    int ans[100][100];
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            cin>>a[i][j];
            b[i][j]=a[i][j];
        }
    }
    m--;
    while(m--){
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                ans[i][j]=0;
                for(int k=1;k<=n;k++){
                    ans[i][j]+=a[i][k]*b[k][j];
                }
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                b[i][j]=ans[i][j];
            }
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            cout<<ans[i][j]<<" ";
        }
        cout<<endl;
    }
}