P2044 随机数生成器

QAQ
递推矩阵：|x[n],1|=|x[n-1],1|*|a,0|
　　　　　　　　　　　　　　|ｃ,1|
但是两个long long 值乘会爆类型，于是把乘的过程分解，在过程中mod m即可AC

#include <cstdio>
#include <iostream> 
#include <cstring>
#define ll long long 
using namespace std;
ll m,a,c,x0,n,g;
ll ans[2][3];
ll dx[3][3];
ll c1[3][3];
ll multi(ll y ,ll cnt ) 
{
    if (!cnt) return 0 ;
    if (cnt==1) return y % m ;
    ll rec=multi( y , cnt / 2 ) ;
    rec=(rec+rec)%m ;
    if (cnt%2) rec=(rec+y)%m ;
    return rec ;
}
void ans_cf()
{
    for(int i=1;i<=1;i++)
     for(int j=1;j<=2;j++)
      c1[i][j]=ans[i][j],ans[i][j]=0;
    for(int i=1;i<=1;i++)
     for(int j=1;j<=2;j++)
      for(int k=1;k<=2;k++)
       ans[i][j]=(ans[i][j]+multi(c1[i][k],dx[k][j]))%m;
}
void x_cf()
{
    for(int i=1;i<=2;i++)
     for(int j=1;j<=2;j++)
      c1[i][j]=dx[i][j],dx[i][j]=0;
    for(int i=1;i<=2;i++)
     for(int j=1;j<=2;j++)
      for(int k=1;k<=2;k++)
       dx[i][j]=(dx[i][j]+multi(c1[i][k],c1[k][j]))%m;  
}
void fast_pow(ll k)
{
    while(k)
    {
        if(k%2) ans_cf();
        k/=2;
        x_cf();
    }
}
int main()
{
    cin>>m>>a>>c>>x0>>n>>g;
    dx[1][1]=a,dx[2][1]=c,dx[2][2]=1;
    ans[1][1]=x0,ans[1][2]=1;
    fast_pow(n);
    printf("%lld",ans[1][1]%g);
}