bzoj2875随机数生成器——矩阵快速幂

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=2875

矩阵快速幂，把x和c分开求，最后加上即可；

为防止爆long long，要用快速乘。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll n,x,a,c,m,g;
ll mul(ll x,ll y)
{
    ll s=;
    while(y)
    {
        if(y&)(s+=x)%=m;
        x=(x+x)%m;
        y/=;
    }
    return s;
}
struct Matrix{
    ll aa[][];
    Matrix operator * (const Matrix &y) const
    {
        Matrix x;
        memset(x.aa,,sizeof x.aa);
        for(int i=;i<=;i++)
            for(int k=;k<=;k++)
                for(int j=;j<=;j++)
                    (x.aa[i][j]+=mul(aa[i][k],y.aa[k][j]))%=m;
        return x;
    }
}ans,s;
int main()
{
    scanf("%lld%lld%lld%lld%lld%lld",&m,&a,&c,&x,&n,&g);
    s.aa[][]=a;s.aa[][]=;
    s.aa[][]=c;s.aa[][]=;
    ans.aa[][]=;ans.aa[][]=;
    while(n)
    {
        if(n&)ans=s*ans;
        s=s*s;
        n/=;
    }
    printf("%lld",(mul(ans.aa[][],x)+ans.aa[][])%m%g);
    return ;
}