社区讨论
编译失败?
P3390【模板】矩阵快速幂参与者 23已保存回复 41
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 41 条
- 当前快照
- 1 份
- 快照标识符
- @ly3879wy
- 此快照首次捕获于
- 2024/07/02 01:00 2 年前
- 此快照最后确认于
- 2024/10/04 15:18 去年
CPP
#include<bits/stdc++.h>
#define ll long long
ll n,k,a[101][101],d[101][101],m=1e9+7;
using namespace std;
void ab(ll a[101][101],ll b[101][101]){
int c[101][101];
memset(c,0,sizeof(c));
for(int i=0;i<n;i++)
for(int k=0;k<n;k++)
for(int j=0;j<n;j++)c[i][j]=(c[i][j]+b[i][k]*a[k][j])%m;
for(int i=0;i<n;i++){
for(int j=0;j<n;j++)a[i][j]=c[i][j];
}return;
}
void qp(ll a[101][101],ll b[101][101]){
while(k){
if(k&1ll)ab(a,b);
ab(b,b);
k>>=1ll;
}return;
}
int main(){
cin>>n>>k;
for(int i=0;i<n;i++){
d[i][i]=1;
for(int j=0;j<n;j++)cin>>a[i][j];
}qp(d,a);
for(int i=0;i<n;i++){
for(int j=0;j<n;j++)cout<<d[i][j]<<' ';
cout<<endl;
}return 0;
}
回复
共 41 条回复,欢迎继续交流。
正在加载回复...