社区讨论
灌佬多
灌水区参与者 2已保存回复 16
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 16 条
- 当前快照
- 1 份
- 快照标识符
- @lyu1jp4i
- 此快照首次捕获于
- 2024/07/20 19:23 2 年前
- 此快照最后确认于
- 2024/07/20 20:29 2 年前
本蒟蒻于近日学习高等数论,结果成功让大脑 CPU 报废。
在此,本蒟蒻祈求众灌佬能给蒟蒻提一提学习高等数论中的建议。
寄:

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 条回复,欢迎继续交流。
正在加载回复...