矩阵详

#include<bits/stdc++.h>
using namespace std;//long long
int n;
long long m;
long long a[105][105];
long long ans[105][105];
long long MOD=1e9+7;
void f(){//ans=ans*a
	long long c[105][105]={0};
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			for(int k=1;k<=n;k++){
				c[i][j]=(c[i][j]+ans[i][k]*a[k][j])%MOD;//因为是矩阵所以乘法是这样
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			ans[i][j]=c[i][j];//将答案赋值给ans
		}
	}
}
void fl(){//a=a*a
	long long c[105][105]={0};
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			for(int k=1;k<=n;k++){
				c[i][j]=(c[i][j]+a[i][k]*a[k][j])%MOD;
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			a[i][j]=c[i][j];
		}
	}
}
long long read(){   
    long long f=1;
    long long x=0;
    char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-') f=-1;
        c=getchar();
    }
    while('0'<=c&&c<='9'){
        x=x*10+c-'0';
        c=getchar();
    }
    return x*f;
}
int main(){
	n=read();
	m=read();
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			a[i][j]=read();
		}
	}
	for(int i=1;i<=n;i++) ans[i][i]=1;//为什么这样初始化：你想想这样任何矩阵*ans都=这个矩阵；（就比如快速幂里把ans初始化为1）
	while(m>0){//这何不是快速幂框架乎？？？？？
		if(m%2){
			f();//ans=ans*a;
		}
		fl();//a=a*a
		m>>=1;
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			cout<<ans[i][j]%MOD<<" ";
		}
		cout<<endl;
	}
	return 0;
}