社区讨论

编译失败?

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 条回复,欢迎继续交流。

正在加载回复...