矩阵快速幂板子

#include<bits/stdc++.h>
using namespace std;
const int N=1010,mod=1e9+7;
long long a[N][N],ans[N][N],temp[N][N],n,k;
void calculate()
{
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		temp[i][j]=ans[i][j];
		ans[i][j]=0;
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	for(int k=1;k<=n;k++)
		ans[i][j]=(ans[i][j]+(a[i][k]*temp[k][j])%mod)%mod;
}
void change()
{
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		temp[i][j]=a[i][j];
		a[i][j]=0;
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	for(int k=1;k<=n;k++)
		a[i][j]=(a[i][j]+(temp[i][k]*temp[k][j])%mod)%mod;
}
void quickpow(long long x)
{
	while(x)
	{
		if(x%2) calculate();
		x/=2;
		change();
	}
}
int main()
{
	scanf("%lld%lld",&n,&k);
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		scanf("%lld",&a[i][j]);
		ans[i][j]=a[i][j];
	}
	quickpow(k-1);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		printf("%lld ",ans[i][j]);
		printf("\n");
	}
	return 0;
}