这是一个矩阵快速幂板子……
为什么会变成这样呢(0pts)
#include<bits/stdc++.h>
using namespace std;
const int Maxn=102;
const int mod=1e9+7;
int n,k;
struct matrix{
long long m[Maxn][Maxn];
}a,b;
matrix operator * (const matrix& a,const matrix& b){
matrix c;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++) c.m[i][k]=(c.m[i][k]+a.m[i][j]*b.m[j][k]%mod)%mod;
return c;
}
int main(){
scanf("%d%d",&n,&k);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
scanf("%lld",&a.m[i][j]);
for(int i=1;i<=n;i++) b.m[i][i]=1;
for(;k;k>>=1){
if(k&1) b=b*a;
a=a*a;
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++)
printf("%lld ",b.m[i][j]);
printf("\n");
}
return 0;
}