第十五届吉林省赛The 15th Jilin Provincial Collegiate Programming Contest C.Random Number Generator(数学 BSGS)
Posted quinn18
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了第十五届吉林省赛The 15th Jilin Provincial Collegiate Programming Contest C.Random Number Generator(数学 BSGS)相关的知识,希望对你有一定的参考价值。
文章目录
题意
:
题解
:
BSGS
代码
:
#include <bits/stdc++.h>
#include <unordered_map>
using namespace std;
#define ll long long
#define int long long
const ll N=1e8+5;
ll ksm(ll a, ll p, ll mod)
ll res=1;
while(p)
if(p&1) res=res*a%mod;
a=a*a%mod;
p>>=1;
return res;
ll exgcd(ll a, ll b, ll &x, ll &y)
if(!b)
x=1;
y=0;
return a;
ll ans=exgcd(b, a%b, y, x);
y=y-a/b*x;
return ans;
ll inv(ll a, ll mod) //存在逆元条件:gcd(a,mod)=1
ll x, y;
ll g=exgcd(a, mod, x, y);
if(g!=1) return -1;
return (x%mod+mod)%mod;
ll bsgs(ll a,ll b,ll p)
b%=p;
if(b==1||p==1)return 0;
ll n=sqrt(p);
static unordered_map<ll,ll>Bmp;
Bmp.clear();
ll inva=inv(ksm(a,n-1,p),p)*b%p;
for(ll i=n-1;i>=0;i--)
Bmp[inva]=i;inva=inva*a%p;
ll ta=1,powa=ksm(a,n,p);
for(ll k=0;k<=p;k+=n)
if(Bmp.count(ta))return k+Bmp[ta];
ta=ta*powa%p;
return -1;
signed main()
ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
int a, b, M, x0, x;
cin>>a>>b>>M>>x0>>x;
bool q=0;
if(x==x0) q=1;
else if(a==0)
if(b==x) q=1;
else q=0;
else if(a==1)
if(b) q=1;
else q=0;
else
int aa=ksm(a-1,M-2,M)%M;
int t=((x+b*aa)%M*ksm(x0+b%M*aa%M,M-2,M)%M)%M;
if(bsgs(a,t,M)==-1) q=0;
else q=1;
if(q) cout<<"YES"<<endl;
else cout<<"NO"<<endl;
return 0;
总结
Qaq
以上是关于第十五届吉林省赛The 15th Jilin Provincial Collegiate Programming Contest C.Random Number Generator(数学 BSGS)的主要内容,如果未能解决你的问题,请参考以下文章
第十五届吉林程序设计省赛The 15th Jilin Provincial Collegiate Programming Contest I. Nim Game(树状数组,线性基)
第十五届吉林程序设计省赛The 15th Jilin Provincial Collegiate Programming Contest I. Nim Game(树状数组,线性基)