社区讨论

板子……但萌新不会……

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

正在加载回复...