2019牛客暑期多校训练营（第五场）B generator 1 十进制矩阵快速幂

题意：，求xn

题解：用矩阵快速幂搞，但是n太大所以得转化为十进制求。然后矩阵乘法那个地方不要mod太多。。。会被卡常

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <bits/stdc++.h>
using namespace std;
typedef long long int ll;
typedef unsigned long long ull;
#define ms(s) memset(s, 0, sizeof(s))
const int INF = 0x3f3f3f3f;
const int maxn = 1e6+5;
ll  mod;
struct mat{
	ll c[2][2];
}E;
mat multi(mat a,mat b){
	 mat ans;
	 ans.c[0][0]=ans.c[0][1]=ans.c[1][0]=ans.c[1][1]=0;
	 for(int i=0;i<2;i++){
	 	for(int j=0;j<2;j++){
	 		for(int k=0;k<2;k++){
	 			ans.c[i][j]+=a.c[i][k]*b.c[k][j];
	 			ans.c[i][j]%=mod;
	 		}
	 	}
	 }
	 return ans;
}
mat quick_pow(mat a,int b){
	mat ans=E;
	while(b){
		if(b&1)ans=multi(ans,a);
		a=multi(a,a);
		b>>=1;
	}
	return ans;
}
char s[maxn];
int main() 
{
	ll  x0,x1,a,b;
	scanf("%lld%lld%lld%lld%s%lld",&x0,&x1,&a,&b,s,&mod);
	mat tmp;
	tmp.c[1][0]=b,tmp.c[1][1]=a;
	tmp.c[0][0]=E.c[0][1]=E.c[1][0]=0;
	tmp.c[0][1]=E.c[0][0]=E.c[1][1]=1;
	mat ans=E;
	for(int i=0;s[i];i++){
		ans=multi(quick_pow(ans,10),quick_pow(tmp,s[i]-'0'));
	}
	printf("%lld\n",(x0*ans.c[0][0]%mod+x1*ans.c[0][1]%mod)%mod);
    return 0;
}