记录编号 |
197157 |
评测结果 |
AAAAAAAAAAAAAAAAAAAA |
题目名称 |
[NOI 2012]随机数生成器 |
最终得分 |
100 |
用户昵称 |
devil |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.008 s |
提交时间 |
2015-10-23 16:43:06 |
内存使用 |
0.31 MiB |
显示代码纯文本
#include <cstdlib>
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <map>
#include <stack>
#include <vector>
#include <map>
#include <queue>
#include <ctime>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef unsigned int uint;
const int inf=1061109567;
const int maxn=1010;
const int maxm=10000010;
const int mod=1000000000;
struct Matrix
{
ull num[3][3];
Matrix() {memset(num,0,sizeof(num));}
} A,I,B,ans;
ull m;
ull mul(ull a,ull b)
{
ull sum=0;
while(b)
{
if(b&1) sum=(sum+a)%m;
a=(a+a)%m;b>>=1;
}
return sum;
}
Matrix operator * (Matrix a,Matrix b)
{
Matrix c;
for(int i=1;i<=2;i++)
{
for(int j=1;j<=2;j++)
{
for(int k=1;k<=2;k++)
{
c.num[i][j]=(c.num[i][j]+mul(a.num[i][k],b.num[k][j]))%m;
}
}
}
return c;
}
Matrix power(Matrix a,ull b)
{
Matrix sum;sum=I;
while(b)
{
if(b&1) sum=sum*a;
a=a*a;b>>=1;
}
return sum;
}
int main()
{
freopen("randoma.in","r",stdin);
freopen("randoma.out","w",stdout);
//clock_t st=clock();
ull a,c,x0,n,g;
cin>>m>>a>>c>>x0>>n>>g;
A.num[1][1]=a;A.num[1][2]=0;
A.num[2][1]=c;A.num[2][2]=1;
B.num[1][1]=x0;B.num[1][2]=1;
I.num[1][1]=1;I.num[1][2]=0;
I.num[2][1]=0;I.num[2][2]=1;
ans=power(A,n);ans=B*ans;
ans.num[1][1]%=g;
cout<<ans.num[1][1]<<endl;
//clock_t ed=clock();
//printf("\nTime used : %.5lf Ms\n",double(ed-st)/CLOCKS_PER_SEC);
return 0;
}