bzoj2875随机数生成器

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

矩阵乘裸题。

如果直接乘的话会爆long long，所以用加法代替乘，过程中不断取模。

加法应是快速乘，不要O(n)循环……

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
ll n,m,g,a,c,ans,x;
ll mul(ll u,ll v)
{
    ll ret=,tp=u;
    while(v)
    {
        if(v&)(ret+=tp)%=m;
        (tp+=tp)%=m;v>>=;
    }
    return ret;
}
struct Matrix{
    ll v[][];
    Matrix operator *(const Matrix &b)const
    {
        Matrix tp;memset(tp.v,,sizeof tp.v);
        for(int i=;i<=;i++)
            for(int j=;j<=;j++)
                for(int k=;k<=;k++)
                    (tp.v[i][j]+=mul(v[i][k],b.v[k][j]))%=m;
        return tp;
    }
}res,ct;
void init()
{
    scanf("%lld%lld%lld%lld%lld%lld",&m,&a,&c,&x,&n,&g);
    n--;
    res.v[][]=a;res.v[][]=;
    res.v[][]=;res.v[][]=;
    memcpy(ct.v,res.v,sizeof res.v);
}
int main()
{
    init();
    while(n)
    {
        if(n&)res=res*ct;
        ct=ct*ct;n>>=;
    }
    ans=(mul(res.v[][],x)+mul(res.v[][],c))%m%g;
    printf("%lld",ans);
    return ;
}