[51nod1113]矩阵快速幂

解题关键：模板题，方便以后熟悉

 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
struct mat{
    ll m[200][200];
};
ll mod=1e9+7;
int n,t;
mat mul(mat &A,mat &B){
    mat C={0};
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            for(int k=0;k<n;k++){
                C.m[i][j]=(C.m[i][j]+A.m[i][k]*B.m[k][j]+mod)%mod;
            }
        }
    }
    return C;
}
mat pow(mat A){
    mat B={0};
    for(int i=0;i<n;i++){
        B.m[i][i]=1;
    } 
    while(t){
        if(t&1) B=mul(B,A);
        A=mul(A,A);
        t>>=1;
    }
    return B;
}
int main(){
    cin>>n>>t;
    mat A={0};
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            cin>>A.m[i][j];
        }
    } 
    A=pow(A);
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            cout<<A.m[i][j]<<' ';
        }
        cout<<endl;
    } 
}