矩阵快速幂求调悬关

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod=1e9+7;
int n,k;
struct mat{
	int size,s[110][110];
	mat operator *(mat m){
		mat mm={size};
		for(int i=1;i<=size;i++)
			for(int j=1;j<=size;j++)
				for(int k=1;k<=size;mm.s[i][j]+=s[i][k]*m.s[k][j]%mod,s[i][j]%=mod,k++);
		return mm;
	}
}e,mt,ans;
mat ppow(mat a,int b){
	if(!b)
		return e;
	if(b==1)
		return a;
	if(b==2)
		return a*a;
	mat mid=ppow(a,b>>1);
	if(b&1)
		return mid*mid*a;
	return mid*mid;
}
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>k;
	for(int i=1;i<=n;e.s[i][i]=1,i++);
	e.size=mt.size=n;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;cin>>mt.s[i][j],j++);
	ans=ppow(mt,k);
	for(int i=1;i<=n;cout<<'\n',i++)
		for(int j=1;j<=n;cout<<ans.s[i][j]<<' ',j++);
	return 0;
}