灌佬多
  • 板块灌水区
  • 楼主Frank2010
  • 当前回复16
  • 已保存回复16
  • 发布时间2024/7/20 19:22
  • 上次更新2024/7/20 20:29:44
查看原帖
灌佬多
857735
Frank2010楼主2024/7/20 19:22

本蒟蒻于近日学习高等数论,结果成功让大脑 CPU 报废。

在此,本蒟蒻祈求众灌佬能给蒟蒻提一提学习高等数论中的建议。

寄: WA声

#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;
}
2024/7/20 19:22
加载中...