BZOJ2875: [Noi2012]随机数生成器 矩阵乘法+快速乘

显然需要用矩阵乘法快速幂
矩阵为：
A 0
C 1

因为mod很大所以直接做乘法会爆掉，在做乘法运算的时候要用一下快速乘

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
struct matrix {ll ju[3][3];}s,a;
ll mod,b,c,x0,n,g;
ll muti(ll p,ll q)
{
	ll d=0;
	while(q)
	{
		if(q%2==1) d=(d+p)%mod;
		p=(p+p)%mod;
		q/=2; 
	}
	return d;
}
matrix operator*(matrix x,matrix y)
{
	matrix ans;
	for(int i=1;i<=2;i++)
	for(int j=1;j<=2;j++)
	ans.ju[i][j]=0;
 	for(int i=1;i<=2;i++)
	for(int j=1;j<=2;j++)
	for(int k=1;k<=2;k++)
	ans.ju[i][j]=(ans.ju[i][j]+muti(x.ju[i][k],y.ju[k][j]))%mod;
	return ans;
}
int main()
{
	scanf("%lld%lld%lld%lld%lld%lld",&mod,&b,&c,&x0,&n,&g);
	s.ju[1][1]=x0,s.ju[1][2]=1;
	a.ju[1][1]=b,a.ju[2][1]=c,a.ju[2][2]=1; 
	while(n)
	{
		if(n%2==1) s=s*a;
		a=a*a;
		n/=2;
	}
	printf("%lld",s.ju[1][1]%g);
}