社区讨论
10pts求调
P3390【模板】矩阵快速幂参与者 2已保存回复 2
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 2 条
- 当前快照
- 1 份
- 快照标识符
- @m3zudu9q
- 此快照首次捕获于
- 2024/11/27 20:08 去年
- 此快照最后确认于
- 2025/11/04 13:48 4 个月前
CPP
#include<iostream>
using namespace std;
long long n,k,ans[101][101],t[101][101],mod=1e9+7;
inline long long read(){
long long x=0,f=1;
char ch=getchar();
while(ch<'0' || ch>'9'){
f=-1;
ch=getchar();
}
while(ch>='0' && ch<='9'){
x=(x<<1)+(x<<3)+(ch^'0');
ch=getchar();
}
return x*f;
}
inline void jzc(bool f){
long long fl[101][101]={};
if(f){
for(int k=1;k<=n;k++){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) fl[i][j]=(fl[i][j]%mod+t[i][k]%mod*t[k][j]%mod)%mod;
}
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) t[i][j]=fl[i][j];
}
}else{
for(int k=1;k<=n;k++){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) fl[i][j]=(fl[i][j]%mod+ans[i][k]%mod*t[k][j]%mod)%mod;
}
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) ans[i][j]=fl[i][j];
}
}
}
int main(){
n=read();
k=read();
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) t[i][j]=read();
}
for(int i=1;i<=n;i++) ans[i][i]=1;
while(k){
if(k&1) jzc(0);
jzc(1);
k>>=1;
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) printf("%lld ",ans[i][j]%mod);
putchar('\n');
}
return 0;
}
回复
共 2 条回复,欢迎继续交流。
正在加载回复...