[BSGS]大步小步算法
问题
BSGS被用于求解离散对数,即同余方程:
\]
求\(x\)的最小非负整数解。
保证\(A\perp P\)(互质)。
分析
首先,我们根据费马小定理,有
\]
则显然有
\]
即
\]
那么显然\(x<P-1\),我们就得到了一个\(O(P)\)的算法,然而太慢了。
考虑分块算法,对\(x\)每\(m\)分一块,则有
\]
移项整理
\]
那么我们枚举\(i\),就可以求出\(A^j\)。再对于\(j\in[0,m-1]\)的\(A^j\)存进哈希表/map,就可以得到\(x=im-j\)了。如果不考虑查询哈希表/map的时间,则时间复杂度为\(O(m+\frac{P}{m})\)。
那\(m\)应该取何值呢?求\(f(m)=m+\frac{P}{m}\)的驻点:
\]
即
\]
移项整理
\]
解得\(m=\sqrt{P}\)。
那么我们令\(m=\lceil\sqrt{P}\rceil\),就得到了一个\(O(\sqrt{P})\)的算法。
代码
\(-1\)为无解。
ll BSGS(ll a,ll b,ll p){
if(!a)return b?-1:1;
if(b==1)return 0;
map<ll,ll>mp;
ll m=ceil(sqrt(p)),ax=1;
for(int i=0;i<m;i++){
mp[ax]=i;
ax=ax*a%p;
}
ll am=pow(a,m,p),aj=am*pow(b,p-2,p)%p;
for(int i=1;i<=m;i++){
if(mp.count(aj))return m*i-mp[aj];
aj=aj*am%p;
}
return -1;
}
例题
#include<iostream>
#include<cstdio>
#include<cmath>
#include<map>
using namespace std;
typedef long long ll;
int t,k;
ll y,z,p;
ll pow(ll a,ll b,ll p){
ll ans=1;
while(b){
if(b&1)ans=ans*a%p;
a=a*a%p;
b>>=1;
}
return ans;
}
ll BSGS(ll a,ll b,ll p){
if(!a)return b?-1:1;
if(b==1)return 0;
map<ll,ll>mp;
ll m=ceil(sqrt(p)),ax=1;
for(int i=0;i<m;i++){
mp[ax]=i;
ax=ax*a%p;
}
ll am=pow(a,m,p),aj=am*pow(b,p-2,p)%p;
for(int i=1;i<=m;i++){
if(mp.count(aj))return m*i-mp[aj];
aj=aj*am%p;
}
return -1;
}
int main(){
scanf("%d%d",&t,&k);
while(t--){
scanf("%lld%lld%lld",&y,&z,&p);
if(k==1)printf("%lld\n",pow(y,z,p));
else if(k==2){
if(y%p==0)printf("Orz, I cannot find x!\n");
else printf("%lld\n",pow(y,p-2,p)*z%p);
}else{
ll ans=BSGS(y%p,z%p,p);
if(~ans)printf("%lld\n",ans);
else printf("Orz, I cannot find x!\n");
}
}
}
[BSGS]大步小步算法的更多相关文章
- [模板]大步小步算法——BSGS算法
大步小步算法用于解决:已知A, B, C,求X使得 A^x = B (mod C) 成立. 我们令x = im - j | m = ceil(sqrt(C)), i = [1, m], j = [0, ...
- 离散对数及其拓展 大步小步算法 BSGS
离散对数及其拓展 离散对数是在群Zp∗Z_{p}^{*}Zp∗而言的,其中ppp是素数.即在在群Zp∗Z_{p}^{*}Zp∗内,aaa是生成元,求关于xxx的方程ax=ba^x=bax=b的解, ...
- 离散对数&&大步小步算法及扩展
bsgs algorithm ax≡b(mod n) 大步小步算法,这个算法有一定的局限性,只有当gcd(a,m)=1时才可以用 原理 此处讨论n为素数的时候. ax≡b(mod n)(n为素数) 由 ...
- 【题解】Matrix BZOJ 4128 矩阵求逆 离散对数 大步小步算法
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4128 大水题一道 使用大步小步算法,把数字的运算换成矩阵的运算就好了 矩阵求逆?这么基础的线 ...
- 大步小步算法模板题, poj2417
大步小步模板 (hash稍微有一点麻烦, poj不支持C++11略坑) #include <iostream> #include <vector> #include <c ...
- BSGS算法(大步小步算法)
计算\(y^x ≡ z \ mod\ p\) 中 \(x\) 的解. 这个模板是最小化了\(x\) , 无解输出\(No \ Solution!\) map<ll,ll>data; ll ...
- BSGS-Junior·大步小步算法
本文原载于:http://www.orchidany.cf/2019/02/06/BSGS-junior/#more \(\rm{0x01}\) \(\mathcal{Preface}\) \(\rm ...
- 洛谷 - P4861 - 按钮 - 扩展大步小步算法
https://www.luogu.org/problemnew/show/P4861 把好像把一开始b==1的特判去掉就可以AC了. #include<bits/stdc++.h> us ...
- 浅谈BSGS(大步小步)及其扩展
用途: 一般用来求\(a^x\equiv b\,\,(mod\,p)\)的最小正整数解,其中gcd(a,p)=1 设\(u=\lceil sqrt(p)\rceil\),则式子可以转化为\(a^{iu ...
随机推荐
- BZOJ1563: [NOI2009]诗人小G(决策单调性 前缀和 dp)
题意 题目链接 Sol 很显然的一个dp方程 \(f_i = min(f_j + (sum_i - sum_j - 1 - L)^P)\) 其中\(sum_i = \sum_{j = 1}^i len ...
- ionic1跨域问题
在使用ionic是遇到跨域问题 我自己尝试简单的单间一个能返回数据的后台服务器,ionic的www目录我使用一个ionic的livereload 来当做一个简易的web服务器,在使用ionic ser ...
- java右移>> 无符号右移>>>
>>>是无符号右移,在高位补零 >>是带符号的右移,如果是正数则在高位补零,负数则补1 int a = -1; System.out.println(a>>1 ...
- Android - 通过真实案例学习解内存泄漏问题,最终发现Android原生Bug
作为一个Android新手小白,刚到新公司,最近的工作就是在学习解各类Bug.转型之初,面临各种新知识,会有压力,但是学习的过程是快乐的. 上周刚遇上一类bug,就是应用的内存泄漏问题.最终通过前辈的 ...
- 我的ORM框架
任何系统的基础,都可以算是各种数据的增删改查(CRUD).最早操作数据是直接在代码里写SQL语句,后来出现了各种ORM框架.C#下的ORM框架有很多,如微软自己的Entity Framework.第三 ...
- js从入门到精通到深入到就业
本篇博客是我参看人家代码做的总结,个人感觉非常非常好,简单.步步深入,不用花大量时间来学完正本js,只需要把其中的代码理解透彻,上班无压力(上班无压力是指js部分,包括查看框架源代码都有很大帮助) / ...
- PHP: 打印post数据,返回长度不为1但内容为空的问题
问题: 首先,html长这样: 我们可以看到textarea标签的name值为cten,那么我们进行查看Post是否能够正常获取到数据: 后台进行获取: 结果:我们可以正常获取到数据 接下来,我们进行 ...
- POJ-3067 Japan---树状数组逆序对变形
题目链接: https://vjudge.net/problem/POJ-3067 题目大意: 日本岛东海岸与西海岸分别有N和M个城市,现在修高速公路连接东西海岸的城市,求交点个数. 解题思路: 记每 ...
- 学习MyBatis之简单入门HelloWorld
转:https://blog.csdn.net/gaomb_1990/article/details/78299784 一.准备 Eclipse:Luna Service Release 1 (4.4 ...
- IBM带库加磁带操作
1.查询要弹出磁带的信息 可查询media日志,冻结,可用等,详情可查 查看带库空闲槽位 vmcheckxxx -rt tld -rn 0(0为带库名) 磁带详细信息: bpmedialist -m ...