NOI2011 随机数生成器

 题目分析：矩阵快速幂模板题。我构造了一个2*2和1*2的矩阵...直接模拟矩阵乘即可。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
#include<cmath>
#include<cctype>
#include<cassert>
#include<climits>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define RepD(i,n) for(int i=n;i>=0;i--)
#define MEM(a) memset(a,0,sizeof(a))
#define MEMI(a) memset(a,127,sizeof(a))
#define MEMi(a) memset(a,128,sizeof(a))
#define INF (2139062143)
#define phiF (1000000006)
#define MAXN (1000000+10)
typedef long long LL;

LL mod,a,c,x0,n,g,A[3][3],B[3],C[3],D[3][3];

LL calc(LL x,LL y){
	long long z=0;
	for (;y;y>>=1){
		if (y&1) (z+=x)%=mod;
		(x+=x)%=mod;
	}
	return z;
}

int main(){
	//freopen("in.txt","r",stdin);
	scanf("%lld%lld%lld%lld%lld%lld",&mod,&a,&c,&x0,&n,&g);
	A[1][1]=a;A[1][2]=1;
	A[2][2]=1;
	B[1]=x0;
	B[2]=c;
	
	while (n){
		if (n%2) {
			C[1]=(calc(A[1][1],B[1])+calc(A[1][2],B[2]))%mod;
			B[1]=C[1];
		}
		n/=2;
		D[1][1]=(calc(A[1][1],A[1][1])+calc(A[1][2],A[2][1]))%mod;
		D[1][2]=(calc(A[1][1],A[1][2])+calc(A[1][2],A[2][2]))%mod;
		D[2][1]=(calc(A[2][1],A[1][1])+calc(A[2][2],A[2][1]))%mod;
		D[2][2]=(calc(A[2][1],A[1][2])+calc(A[2][2],A[2][2]))%mod;
		A[1][1]=D[1][1];A[1][2]=D[1][2];
		A[2][1]=D[2][1];A[2][2]=D[2][2];
	}
	printf("%lld",B[1]%g);
}