社区讨论
只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 条回复,欢迎继续交流。
正在加载回复...