[NOI2012]随机数生成器

题解：

很显然是一道矩阵优化dp

然而表示我很智障地把式子一个个带入

然后就发现了为什么会有那些部分分（大概用扩欧是70吧）

注意用矩阵计算的时候要用快速乘（当然想写高精那也随便，时间无限宽裕）

代码：

#include <bits/stdc++.h>
using namespace std;
#define ll long long
ll m,aa,c,x,n,g,mo1;
struct re{
    ll jz[][];
}a;
ll ksc(ll x,ll y)
{
    if (y==) return();
    if (y==) return(x);
    ll tmp=ksc(x,y/);
    tmp=(tmp*)%mo1;
    if (y%==) tmp=(tmp+x)%mo1;
    return(tmp);
}
re XX(re x,re y)
{
    re tmp;
    memset(tmp.jz,,sizeof(tmp.jz));
    for (ll i=;i<=;i++)
      for (ll j=;j<=;j++)
        for (ll k=;k<=;k++)
          tmp.jz[i][k]+=ksc(x.jz[i][j],y.jz[j][k]),
          tmp.jz[i][k]%=mo1;
    return(tmp);
}
re fast_pow(ll x)
{
    if (x==) return(a);
    re b=fast_pow(x/);
    b=XX(b,b);
    if (x%) b=XX(b,a);
    return(b);
}
int main()
{
    cin>>mo1>>aa>>c>>x>>n>>g;
    a.jz[][]=aa; a.jz[][]=c;
    a.jz[][]=; a.jz[][]=;
    re cc=fast_pow(n);
    ll ans=((ksc(cc.jz[][],x)+cc.jz[][])%mo1)%g;
    cout<<ans;
    return ;
}