bzoj 2875: [Noi2012]随机数生成器

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #define ll long long
 using namespace std;
 ll x[][],y[][],m,a,c,x0,n,g;
 ll chen1(ll a1,ll a2)
 {
     long long a3;
     if(!a2)
       return ;
     a3=chen1(a1,a2>>);
     a3=(a3+a3)%m;
     if(a2&)
       a3=(a3+a1)%m;
     return a3;
 }
 void chen(ll x1[][],ll x2[][])
 {
     long long x3[][];
     memset(x3,,sizeof(x3));
     for(int i=;i<;i++)
       for(int j=;j<;j++)
         for(int k=;k<;k++)
           x3[i][j]=(x3[i][j]+chen1(x1[i][k],x2[k][j]))%m;
     for(int i=;i<;i++)
       for(int j=;j<;j++)
         x1[i][j]=x3[i][j];
 }
 int main()
 {
     scanf("%lld%lld%lld%lld%lld%lld",&m,&a,&c,&x0,&n,&g);
     x[][]=;
     x[][]=;
     y[][]=a;
     y[][]=c;
     y[][]=;
     for(;n;)
       {
         if(n%)
           chen(x,y);
         chen(y,y);
         n>>=;
       }
     printf("%lld\n",(chen1(x[][],x0)+x[][])%m%g);
     return ;
 }

矩阵乘法