社区讨论

只AC#11求调

P3390【模板】矩阵快速幂参与者 2已保存回复 5

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
5 条
当前快照
1 份
快照标识符
@mhizej91
此快照首次捕获于
2025/11/03 18:12
4 个月前
此快照最后确认于
2025/11/03 18:12
4 个月前
查看原帖
CPP
#include<bits/stdc++.h>
#define int long long
using namespace std;
const long long mod=1e9+7;
long long n,k,a[105][105][35],ans[105][105],r[105][105];
signed main(){
	cin>>n>>k;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			cin>>a[i][j][0];
	for(int i=1;i<=n;i++)
		ans[i][i]=1;
	for(int p=1;p<=32;p++){
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				for(int kk=1;kk<=n;kk++)
					a[i][j][p]=(a[i][j][p]+a[i][kk][p-1]*a[kk][j][p-1])%mod;
	}
	for(int p=32;p>=0;p--){
		if((k>>p)&1){
			for(int i=1;i<=n;i++)
				for(int j=1;j<=n;j++)
					r[i][j]=0;
			for(int i=1;i<=n;i++)
				for(int j=1;j<=n;j++)
					for(int kk=1;kk<=n;kk++)
						r[i][j]=(r[i][j]+ans[i][kk]*a[kk][j][p])%mod;
			for(int i=1;i<=n;i++)
				for(int j=1;j<=n;j++)
					ans[i][j]=r[i][j];
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)
			cout<<ans[i][j]<<" ";
		cout<<"\n";
	}
	return 0;
}

回复

5 条回复,欢迎继续交流。

正在加载回复...