首页 > 学院 > 开发设计 > 正文

[BZOJ3239][BSGS]Discrete Logging

2019-11-11 01:24:30
字体:
来源:转载
供稿:网友

题意


Ax≡B(mod P)的最小非负整数解


BSGS裸题

自己推一下:

x=im+jm=⌈q√⌉

AimAj≡B(mod P)

Aj≡B∗ine(Aim)(mod P)

已知ine(1)=1

由费马小定理得ine(Aim)=ine(A(i−1)m)∗Ap−1−m mod P

那么ine就可以处理出来了。

A0...mhash存一下

再求B∗ine(A(1...m)m)在hash表中查找。

#include <cstdio>#include <map>#include <iostream>#include <cmath>using namespace std;typedef long long ll;int p,b,n;map<int,int> Mp;inline int powf(ll x,int y,int p){ int k=1; x%=p; while(y){ if(y&1) k=k*x%p; x=x*x%p; y>>=1; } return k;}inline void solve(int y,int z,int p){ y%=p; if(!y&&!z) {puts("1");return;} if(!y) {puts("no solution");return;} ll t=ceil(sqrt(p)),k=1,ine=1; Mp.clear(); Mp[1]=0; for(int i=1;i<t;i++){ k=1ll*k*y%p; if(Mp.count(k)) continue; Mp[k]=i; } int tmp=powf(y,p-t-1,p); for(int i=0;i<t;i++){ if(Mp.count(z*ine%p)) {PRintf("%d/n",i*t+Mp[z*ine%p]);return;} ine=1ll*ine*tmp%p; } puts("no solution");}int main(){ while(~scanf("%d%d%d",&p,&b,&n)) solve(b,n,p);}
上一篇:Key Concepts (Mastering CMake)

下一篇:poj1046

发表评论 共有条评论
用户名: 密码:
验证码: 匿名发表