bzoj2875: [Noi2012]随机数生成器 裸矩阵乘法

裸题，注意细节就可以了，不再赘述。

#include <iostream>#include <algorithm>#include <cstdio>#include <cstring>#include <cmath>using namespace std;typedef  long long sint;struct matrix{    sint a[2][2];    void cl()    {        memset(a,0,sizeof(a));    }};sint x0,a,c,g,m,n;sint mul(sint X,sint Y){    sint re=0;    while(Y)    {        if(Y&1) re+=X,re=(re>m?re-m:re);        X+=X,X=(X>m?X-m:X);        Y>>=1;    }    return re;}matrix operator * (matrix aa,matrix bb){    matrix cc;cc.cl();     for(int i=0;i<2;i++)    {        for(int j=0;j<2;j++)        {            for(int k=0;k<2;k++)            {                cc.a[i][j]+=mul(aa.a[i][k],bb.a[k][j]);                if(cc.a[i][j]>m) cc.a[i][j]-=m;            }        }    }    return cc;}sint quick(){    matrix res;    res.cl();    res.a[0][0]=1;res.a[1][1]=1;    matrix x,y;x.cl();y.cl();    x.a[0][0]=x0;x.a[0][1]=c;    y.a[0][0]=a;y.a[1][0]=y.a[1][1]=1;    while(n)    {        if(n&1)        {            res=res*y;        }        y=y*y;        n>>=1;    }    x=x*res;    return x.a[0][0]%g;}int main(){    //freopen("random20.in","r",stdin);    cin>>m>>a>>c>>x0>>n>>g;    a%=m;c%=m;    cout<<quick();    return 0;}