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

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

矩阵为：

A 0

C 1

因为mod很大所以直接做乘法会爆掉，在做乘法运算的时候要用一下快速乘#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);
}