社区讨论

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

正在加载回复...