社区讨论

灌佬多

灌水区参与者 2已保存回复 16

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
16 条
当前快照
1 份
快照标识符
@lyu1jp4i
此快照首次捕获于
2024/07/20 19:23
2 年前
此快照最后确认于
2024/07/20 20:29
2 年前
查看原帖
本蒟蒻于近日学习高等数论,结果成功让大脑 CPU 报废。
在此,本蒟蒻祈求众灌佬能给蒟蒻提一提学习高等数论中的建议。
寄: WA声
CPP
#include<iostream>
#include<cstring>
#include<cstdio>
const int N=105,mod=1e9+5;
int n;
long long k;
struct matrix{
    int a[N][N];
    matrix operator*(const matrix& you)const{
        matrix res;
        for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){
            res.a[i][j]=0;
            for(int k=1;k<=n;k++)res.a[i][j]=(res.a[i][j]+(long long)a[i][k]*you.a[k][j])%mod;
        }
        return res;
    }
}a,ans;
matrix pow(matrix a,int k){
    matrix res,now=a;
    memset(res.a,0,sizeof res.a);
    for(int i=1;i<=n;i++)res.a[i][i]=1;
    while(k){
        if(k&1)res=res*now;
        now=now*now;
        k>>=1;
    }
    return res;
}
int main(){
    scanf("%d%lld",&n,&k);
    for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)scanf("%d",&a.a[i][j]);
    ans=pow(a,k);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)printf("%d ",ans.a[i][j]);
        printf("\n");
    }
    return 0;
}

回复

16 条回复,欢迎继续交流。

正在加载回复...