社区讨论

80pts(TLE)

P1397[NOI2013] 矩阵游戏参与者 1已保存回复 0

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
0 条
当前快照
1 份
快照标识符
@mlhfzf05
此快照首次捕获于
2026/02/11 11:00
4 周前
此快照最后确认于
2026/02/11 12:30
4 周前
查看原帖
CPP
#include<bits/stdc++.h>
using namespace std;
#define int long long
int a,b,c,d,mod=1000000007;
string n,m;
struct jz{
	int a[5][5];
	jz operator *(jz b){
		jz c;memset(c.a,0,sizeof(c.a));
		for(int i=1;i<=2;i++)for(int j=1;j<=2;j++)for(int k=1;k<=2;k++){
			c.a[i][j]=(c.a[i][j]+a[i][k]*b.a[k][j])%mod;
		}return c;
	}
}r1,r2,ans;
jz qp(jz base,string p){
	int n=p.size()-1;jz res;memset(res.a,0,sizeof(res.a));
	for(int i=1;i<=2;i++)res.a[i][i]=1;
	for(int i=n;i>0;i--){
		if(p[i]=='9'){jz tem=base*base;tem=tem*tem;tem=tem*tem;tem=tem*base;res=res*tem;}
		else if(p[i]=='8'){jz tem=base*base;tem=tem*tem;tem=tem*tem;res=res*tem;}
		else if(p[i]=='7'){jz tem=base*base;tem=tem*tem*base*base*base;res=res*tem;}
		else if(p[i]=='6'){jz tem=base*base;tem=tem*tem*base*base;res=res*tem;}
		else if(p[i]=='5'){jz tem=base*base;tem=tem*tem*base;res=res*tem;}
		else if(p[i]=='4'){jz tem=base*base;tem=tem*tem;res=res*tem;}
		else if(p[i]=='3'){jz tem=base*base*base;res=res*tem;}
		else if(p[i]=='2'){jz tem=base*base;res=res*tem;}
		else if(p[i]=='1'){res=res*base;}
		jz tem=base;base=base*base;base=base*base;base=base*base;base=base*tem*tem;
	}
	return res;
}
signed main(){
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	cin>>n>>m>>a>>b>>c>>d;n=' '+n;m=' '+m;
	r1.a[1][1]=c,r2.a[1][1]=a;
	r1.a[1][2]=d,r2.a[1][2]=b;
	r1.a[2][2]=r2.a[2][2]=1;
	ans.a[2][1]=ans.a[1][1]=1;
	int tem=m.size()-1;
	while(m[tem]=='0')m[tem--]='9';
	m[tem]--;
	jz yh=qp(r2,m);
	r2=r1*yh;tem=n.size()-1;
	while(n[tem]=='0')n[tem--]='9';
	n[tem]--;ans=qp(r2,n)*ans;
	ans=yh*ans;cout<<ans.a[1][1];
	return 0;
}

回复

0 条回复,欢迎继续交流。

正在加载回复...