社区讨论
爆零,T了,求助
P3390【模板】矩阵快速幂参与者 2已保存回复 2
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 2 条
- 当前快照
- 1 份
- 快照标识符
- @mi6yhxw4
- 此快照首次捕获于
- 2025/11/20 12:53 4 个月前
- 此快照最后确认于
- 2025/11/20 12:53 4 个月前
为什么都是T啊,哪里错了
CPP#include <bits/stdc++.h>
#define mod 1000000007
using namespace std;
long long read( ){
long long x=0,y=1;
char c=getchar( );
while(c>'9'||c<'0'){if(c=='-')y=-1;c=getchar( );}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar( );}
return x*y;
}
long long n,k;
struct edge{
long long m[101][101];
}e;
edge time(edge a,edge b){
edge r;
for(register int i=1;i<=n;i++)
for(register int j=1;j<=n;j++)
{
r.m[i][j]=0;
for(register int k=1;k<=n;i++)
r.m[i][j]=(r.m[i][j]+a.m[i][k]*b.m[k][j])%mod;
}
return r;
}
edge fastm(edge a,long long b){
edge r=a,base=a;
b--;
while(b!=0){
if(b%2==1) r=time(r,base);
base=time(base,base);
b>>=1;
}
return r;
}
int main( ){
n=read( );k=read( );
for(register int i=1;i<=n;i++)
for(register int j=1;j<=n;j++){
e.m[i][j]=read( );
}
e=fastm(e,k);
for(register int i=1;i<=n;i++){
for(register int j=1;j<=n;j++)
cout<<e.m[i][j]<<" ";
cout<<endl;
}
return 0;
}
回复
共 2 条回复,欢迎继续交流。
正在加载回复...