社区讨论

爆零,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 条回复,欢迎继续交流。

正在加载回复...