Bzoj 3122 随机数生成器

输入含有多组数据，第一行一个正整数T，表示这个测试点内的数据组数。
接下来T行，每行有五个整数p，a，b，X1，t，表示一组数据。保证X1和t都是合法的页码。注意：P一定为质数
output

共T行，每行一个整数表示他最早读到第t页是哪一天。如果他永远不会读到第t页，输出-1。
input
3
7 1 1 3 3
7 2 2 2 0
7 2 2 2 1
output
1
3
-1
hint 0<=a<=P-1,0<=b<=P-1,2<=P<=10^9
a^(n−1)≡(Xn+b*inv(a−1))*inv(X1+b*inv(a−1))(modp) 推导出来这个式子，代码一直写不对

//#include <bits/stdc++.h>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#define X 10005
#define inF 0x3f3f3f3f
#define PI 3.141592653589793238462643383
#define IO  ios::sync_with_stdio(false),cin.tie(0), cout.tie(0);
//#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
typedef long long ll;
typedef unsigned long long Ull; //2^64
const int maxn = (int)2*1e7 + 10;
//const int MOD = 9973;//(int)1e9 + 7;
//const ll inf = 9223372036854775807;
//const int mod=76532;
int MOD;
void ex_gcd(ll a, ll b, ll &d, ll &x, ll &y) { if (!b) { x = 1; y = 0; d = a; } else { ex_gcd(b, a%b, d, y, x); y -= x * (a / b); }; }
ll gcd(ll a, ll b) { return b ? gcd(b, a%b) : a; }
ll lcm(ll a, ll b) { return b / gcd(a, b)*a; }
ll inv_exgcd(ll a, ll m) { ll d, x, y;ex_gcd(a, m, d, x, y);return d == 1 ? (x + m) % m : -1; }
ll inv1(ll b) { return b == 1 ? 1 : (MOD - MOD / b)*inv1(MOD%b) % MOD; }    //hdu1576用这个板子会爆除0错误
ll crt(int n,int *c,int *m){ll M=1,ans=0;for(int i=0;i<n;++i) M*=m[i];
for(int i=0;i<n;++i) ans=(ans+M/m[i]*c[i] %M *inv_exgcd(M/m[i],m[i]))%M; return ans;}
ll N;
const int mod = (int)1e5 + 10;
ll hs[mod],next[mod],head[mod],id[mod],top;
int read()
{int x=0;char ch=getchar();while(ch<'0'||ch>'9')ch=getchar();while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x;
}
int read()
{int x=0;char ch=getchar();while(ch<'0'||ch>'9')ch=getchar();while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x;
}
int T;
ll p,a,b,x1,t;
ll exgcd(ll a,ll b,ll &x,ll &y)
{if(b==0){x=1;y=0;return a;}ll tmp=exgcd(b,a%b,x,y);ll t=x;x=y;y=t-a/b*y;return tmp;
}
ll kpow(ll a,ll b,ll p)
{a%=p;ll ans=1;for(ll i=b;i;i>>=1,a=a*a%p)if(i&1)ans=ans*a%p;return ans;
}void insert(ll x,ll y)
{ ll k=x%mod;hs[top]=x,id[top]=y,next[top]=head[k],head[k]=top++;
}
int find(ll x)
{int k=x%mod;for(ll i=head[k];i!=-1;i=next[i])if(hs[i]==x) return id[i];
return -1;
}
ll BSGS(ll a,ll b,ll n)
{memset(head,-1,sizeof(head));top=1;if(b==1)return 0;int m=sqrt(n*1.0),j;ll x=1,p=1;for(int i=0;i<m;++i,p=p*a%n)insert(p*b%n,i);for(ll i=m; ;i+=m){if((j=find(x=x*p%n))!=-1)return i-j;if(i>n)break;}return -1;
}
ll cal1()
{ll C=(t-x1+p)%p,x,y;ll t=exgcd(b,p,x,y);if(C%t)return -1;C/=t;x=x*C%p;if(x<0)x+=p;return x+1;
}
ll cal2()
{ll c=inv_exgcd(a-1,p),A=(x1+b*c)%p;//kpow(a-1,p-2,p),A=(x1+b*c)%p,也可以用费马小定理求出逆元，ll C=(t+b*c)%p,x,y;ll t=exgcd(A,p,x,y);x=x%p+p;t=BSGS(a,x*C%p,p);if(t!=-1)return t+1;return -1;
}
ll solve()
{if(x1==t)return 1;if(a==0){if(b==t)return 2;return -1;}if(a==1)return cal1();else return cal2();
}
int main()
{T=read();while(T--){p=read();a=read();b=read();x1=read();t=read();printf("%lld\n",solve());}return 0;
}

Bzoj 3122 随机数生成器相关推荐

【BZOJ】【3671】【NOI2014】随机数生成器
贪心嗯--其实生成这个矩阵就是一个$O(n^2)$的模拟 = = 然后?字典序最小?贪心呗= =能选1就选1,然后能选2就选2-- 我们发现,对于矩阵(1,1)~(n,m),假设1的位置是(x,y) ...
luogu P3306 [SDOI2013] 随机数生成器（BSGS，数列求通项，毒瘤特判）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划发个水题的题解证明我还在() luogu P3306 [SDOI2013] 随机数生成器 Webli ...
Java中的随机数生成器：Random，ThreadLocalRandom，SecureRandom
Java中的随机数生成器:Random,ThreadLocalRandom,SecureRandom 文中的 Random即:java.util.Random, ThreadLocalRandom 即 ...
开源Math.NET基础数学类库使用(13)C#实现其他随机数生成器
原文:[原创]开源Math.NET基础数学类库使用(13)C#实现其他随机数生成器本博客所有文章分类的总目录:http://www.cnblogs.com/asxiny ...
P5147 随机数生成器 [数列]
P5147 随机数生成器数学老师看不懂系列看题目这一片代码就很晕: int work(int x) {if(x==1)return 0;else return work(rand(1,x))+1; ...
UA STAT675 统计计算I 随机数生成1 随机数生成器的一般理论
UA STAT675 统计计算I 随机数生成1 随机数生成器的一般理论 RNG的抽象表示 RNG的质量指标 RNG的统计检测在统计计算中,从某个分布中进行采样通常分为两个步骤: 生成随机数z1,z2 ...
boost::sort模块实现提供多种分布的灵活随机数生成器的测试程序
boost::sort模块实现提供多种分布的灵活随机数生成器的测试程序实现功能 C++实现代码实现功能 boost::sort模块实现提供多种分布的灵活随机数生成器的测试程序 C++实现代码 #i ...
boost::sort模块实现支持不同分布的随机数生成器的测试程序
boost::sort模块实现支持不同分布的随机数生成器的测试程序实现功能 C++实现代码实现功能 boost::sort模块实现支持不同分布的随机数生成器的测试程序 C++实现代码 #inclu ...
ITK:Mersenne Twister随机数生成器
ITK:Mersenne Twister随机数生成器内容提要 C++实现代码内容提要产生一个随机数 C++实现代码 #include "itkMersenneTwisterRandom ...

Bzoj 3122 随机数生成器

Bzoj 3122 随机数生成器相关推荐

最新文章

热门文章