题解 | #矩阵幂#

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