poj 2417 Discrete Logging 数论baby_step,giant_step算法
来源:互联网 发布:参获取视频长度 python 编辑:程序博客网 时间:2024/05/19 09:49
题意:
给p,b,n求最小的l使b^l==n(mod p)。
题意:
相当于在0~p-1内搜索l满足同余式,baby_step,giant_step相当于分块的二分搜索,设m=sqrt(p), 则将p分为m*m,在m块内每块进行m个数的二分搜索。
代码:
//poj 2417//sep9#include <iostream>#include <cmath>#include <algorithm>using namespace std;const int maxN=100000;typedef long long ll;struct Node{int ids;ll v;}baby[maxN];int cmp(Node x,Node y){if(x.v!=y.v)return x.v<y.v;return x.ids<y.ids;}ll pow_mod(ll a,ll b,ll mod){ll ans=1;while(b){if(b&1) ans=(ans*a)%mod;a=(a*a)%mod;b>>=1;}return ans;}int bsearch(int len,ll v){int l=0,r=len,mid;while(l<r){mid=(l+r)/2;if(baby[mid].v==v)return baby[mid].ids;else if(baby[mid].v<v)l=mid+1;elser=mid;}return -1;}int main(){ll p,b,n;while(scanf("%I64d%I64d%I64d",&p,&b,&n)==3){int m=(int)ceil(sqrt((p-1)*1.0));baby[0].ids=0,baby[0].v=1;for(int i=1;i<m;++i){baby[i].ids=i;baby[i].v=(baby[i-1].v*b)%p;}sort(baby,baby+m,cmp);int cnt=1;for(int i=1;i<m;++i)if(baby[i].v!=baby[i-1].v)baby[cnt++]=baby[i];ll bm=pow_mod(pow_mod(b,p-2,p),m,p);int ans=-1;for(int i=0;i<m;++i){int j=bsearch(cnt,n);if(j!=-1){ans=i*m+j;break;}n=(n*bm)%p;}if(ans<0)puts("no solution");else printf("%d\n",ans);}return 0;}
0 0
- poj 2417 Discrete Logging 数论baby_step,giant_step算法
- POJ 2417 Discrete Logging (baby_step,giant_step算法)
- poj 2417 Discrete Logging 求解模方程a^x=b(mod n),n为素数+模板题(baby_step giant_step)
- 【模板】baby_step giant_step poj 2417
- pku3243(baby_step, giant_step算法)
- poj_2417 (baby_step,giant_step算法)
- [poj 2417]Discrete Logging 数论 BSGS
- poj 2417 Discrete Logging(A^x=B(mod c),普通baby_step)
- poj 2417 baby_step giant_step 高次同余方程
- poj 2417 Discrete Logging(baby-step算法)
- POJ 2417/BZOJ 3239(Discrete Logging-BSGS)[Template:数论]
- Discrete Logging[POJ 2417]
- POJ 2417 Discrete Logging
- POJ 2417 Discrete Logging
- poj 2417 Discrete Logging
- POJ 2417 Discrete Logging BSGS
- [POJ 2417]Discrete Logging:BSGS
- POJ 2417 Discrete Logging 笔记
- Pthread创建线程后必须使用join或detach释放线程资源
- Jenkins不同job之间依赖关系的两种配置方式
- 析构函数
- 图像类型、MATLAB数据类型,以及类型之间的转换(二)
- 一些NGINX配置
- poj 2417 Discrete Logging 数论baby_step,giant_step算法
- 作业4.21
- 6、spss做logistic回归
- 菜鸟的安卓实习之路----android广播(3)自定义广播
- CocoaPods是什么
- 关于导入maven项目常见错误
- noj 1076 机器狗组装费用(优先队列)
- jQuery中$.get、$.post、$.getJSON和$.ajax的用法
- 终端设备开发框架