BZOJ2875: [Noi2012]随机数生成器 矩阵乘法+快速乘

→原文链接←

显然需要用矩阵乘法快速幂

矩阵为：
A 0
C 1

因为mod很大所以直接做乘法会爆掉，在做乘法运算的时候要用一下快速乘

<span style="font-size:14px;">#include<cmath>  
#include<cstdio>  
#include<cstring>  
#include<iostream>  
#include<algorithm>  
using namespace std;  
#define ll long long  
struct matrix {ll ju[3][3];}s,a;  
ll mod,b,c,x0,n,g;  
ll muti(ll p,ll q)  
{  
    ll d=0;  
    while(q)  
    {  
        if(q%2==1) d=(d+p)%mod;  
        p=(p+p)%mod;  
        q/=2;   
    }  
    return d;  
}  
matrix operator*(matrix x,matrix y)  
{  
    matrix ans;  
    for(int i=1;i<=2;i++)  
    for(int j=1;j<=2;j++)  
    ans.ju[i][j]=0;  
    for(int i=1;i<=2;i++)  
    for(int j=1;j<=2;j++)  
    for(int k=1;k<=2;k++)  
    ans.ju[i][j]=(ans.ju[i][j]+muti(x.ju[i][k],y.ju[k][j]))%mod;  
    return ans;  
}  
int main()  
{  
    scanf("%lld%lld%lld%lld%lld%lld",&mod,&b,&c,&x0,&n,&g);  
    s.ju[1][1]=x0,s.ju[1][2]=1;  
    a.ju[1][1]=b,a.ju[2][1]=c,a.ju[2][2]=1;   
    while(n)  
    {  
        if(n%2==1) s=s*a;  
        a=a*a;  
        n/=2;  
    }  
    printf("%lld",s.ju[1][1]%g);  
}  </span>