我以[半仙胡小桃]的题解为基础改了，但最后还是0分(代码很像，但非抄袭，只是找了个不用重载运算符的方法学习一下）(https://www.luogu.com.cn/article/2syex3s8)

#include<bits/stdc++.h>
using namespace std;

long long a[999][999],sum[999][999],c[999][999];
const int m=1e9+7;
void jzcf2(long long n){
	for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      c[i][j]=sum[i][j],sum[i][j]=0;
    
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      for(int k=1;k<=n;k++)
       sum[i][j]=(sum[i][j]+(a[i][k]*c[k][j])%m)%m;
}
void jzcf1(long long n){
	for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      c[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]+(c[i][k]*c[k][j])%m)%m;
}
void fastpow(long long n,long long k){
	while(k){
		if(k%2==1){
			jzcf2(n);
		}
		k/=2;
		jzcf1(n);
	}
}
int main(){long long n,k;
	cin>>n>>k;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			cin>>a[i][j];
			sum[i][j]=a[i][j];
		}
	}
	fastpow(n,k-1);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			cout<<sum[i][j]<<" ";
		}
		cout<<endl;
	}
	return 0;
}