社区讨论
板子……但萌新不会……
P3390【模板】矩阵快速幂参与者 6已保存回复 26
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 26 条
- 当前快照
- 1 份
- 快照标识符
- @lo2ci3wr
- 此快照首次捕获于
- 2023/10/23 11:35 2 年前
- 此快照最后确认于
- 2023/11/03 11:44 2 年前
这是一个矩阵快速幂板子……
为什么会变成这样呢(0pts)
CPP#include<bits/stdc++.h>
using namespace std;
const int Maxn=102;
const int mod=1e9+7;
int n,k;
struct matrix{
long long m[Maxn][Maxn];
}a,b;
matrix operator * (const matrix& a,const matrix& b){
matrix c;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++) c.m[i][k]=(c.m[i][k]+a.m[i][j]*b.m[j][k]%mod)%mod;
return c;
}
int main(){
scanf("%d%d",&n,&k);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
scanf("%lld",&a.m[i][j]);
for(int i=1;i<=n;i++) b.m[i][i]=1;
for(;k;k>>=1){
if(k&1) b=b*a;
a=a*a;
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++)
printf("%lld ",b.m[i][j]);
printf("\n");
}
return 0;
}
回复
共 26 条回复,欢迎继续交流。
正在加载回复...