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 ;
}

矩阵乘法

上一篇:Eclipse 打不开


下一篇:js代码大全(各种方法、属性)《转载》