求助 70pts!!!

#include<bits/stdc++.h>
using namespace std;
#define int long long 
int m,a,c,x0,x1,g,n;
struct node{
	int Matrix[3][3];
}res,ans,mid;//初始矩阵，最终矩阵，转移矩阵
int slowly_cheng(int x,int y,int mod){
	int tmp=0;
	while(y){
		if(y&1) tmp=(tmp+x)%mod;
		x=(x+x)%mod;
		y>>=1;
	}
	return tmp%mod;//龟速乘
}
node X(node a,node b){
	node c;
	memset(c.Matrix,0,sizeof(c.Matrix));
	for(int i=1;i<=2;i++){
		for(int j=1;j<=2;j++){
			for(int k=1;k<=2;k++){
				c.Matrix[i][j]=(c.Matrix[i][j]+slowly_cheng(a.Matrix[i][k],b.Matrix[k][j],m))%m;
			}
		}
	}//矩阵乘法
	return c;
}
node quickly_pow(node a,int p){
	while(p){
		if(p&1) res=X(res,a);
		a=X(a,a);
		p>>=1;
	}
	return res;
}
signed main(){
	cin>>m>>a>>c>>x0>>n>>g;
	x1=(a*x0%m+c%m)%m;
	res.Matrix[1][1]=x1;
	res.Matrix[2][1]=x0;
	res.Matrix[1][2]=1;
	res.Matrix[2][2]=1;
	mid.Matrix[1][1]=a;
	mid.Matrix[2][1]=c;
	mid.Matrix[1][2]=0;
	mid.Matrix[2][2]=1;//预处理初始矩阵和转移矩阵
	ans=quickly_pow(mid,n-1);
	cout<<ans.Matrix[1][1]%g;//第一列第一行是答案
	return 0;
}