BSGS及其扩展
之前写了一篇关于BSGS的学习笔记。因为太过老旧,就想修改一些错误,顺便添上扩展BSGS的部分。可惜博客园不能对已发布的随笔修改编辑器,索性重新发出来。旧文已删。
定义
Baby-Step-Giant-Step算法,简称BSGS算法,又称大步小步算法,用于求方程\(a^x\equiv b(\text{mod }c)\)的最小非负整数解,此过程又称离散对数。
普通BSGS算法只能解决\((a,c)=1\)的情况,对其进行扩展之后则没有这个限制。
原理
当\((a,c)=1\)时,若\(c|a\)且\(b\not=0\),则显然无解。若\(c\not|\ \ a\),有\(a^{c-1}\equiv 1(\text{mod }c)\),即\(c-1\)是一个循环节。因此如果方程有解,必定在\([0,c-2]\)中。
朴素算法
由此得到一个朴素的算法。枚举\([0,c-2]\)中每一个数,如果相等就输出答案。时间复杂度\(O(c)\)。考虑对其进行优化,引入数论分块的思想。
数论分块
设\(x=i\times m+j(0\leq j<m)\),原方程改写为\(a^{im+j}\equiv b(\text{mod }c)\)。令\(D(i)=a^{im}\),原方程进一步改写为\(D(i)\cdot a^j\equiv (\text{mod }c)\)。
暴力枚举每个可能的\(i\),用扩展欧几里得可以求出\(a^j\)。可是我们依然无法求得\(j\)。此时注意到\(j\)的所有可能取值只有\(m\)个,因此我们预先将所有的\(a^j\)存入一个哈希表,当我们确定\(a^j\)时,直接查表可以\(O(1)\)得到是否有对应的\(j\),以及对应的\(j\)是多少。暴力枚举复杂度\(O(ilog_2c)\),预处理复杂度\(O(m)\),总复杂度\(O(ilog_2c+m)\)。注意到\(x\approx i\times m\),因此当\(m\)取\(\sqrt c\)时(事实上由于\(\log_2c\)的存在不是很严谨)总复杂度最小,为\(O(\sqrt c\ \log_2c)\)。
至此,我们采用以空间换时间,最后通过分块的思想将该算法优化至根号级别。
例题
Luogu2485 [SDOI2011]计算器
题解
操作1:快速幂;操作2:扩展欧几里得;操作3:普通BSGS。
代码
#include<bits/stdc++.h>
#include<tr1/unordered_map>
#define LL long long
using namespace std;
using namespace std::tr1;
unordered_map<int,int>mp;
int t,k;
inline LL qpow(LL a,LL k,LL p){
LL ans=1;
for(;k;a=a*a%p,k>>=1){if(k&1){ans=ans*a%p;}}
return ans;
}
inline LL inv(LL a,LL p){return qpow(a,p-2,p);}
inline LL BSGS(LL a,LL b,LL p){
if(p==1&&b==0)return 0;
if(b==1)return 0;
if(a==0||b==0){
if(a==0&&b==0)return 1;
return -1;
}
LL i,y=(LL)sqrt(p),tmp,ans;
mp.clear();
for(i=0,tmp=1;i<y;i++,tmp=tmp*a%p){
mp[b*tmp%p]=i;
}
for(i=1,ans=tmp;y*(i-1)+1<=p-2;i++,ans=ans*tmp%p){
if(mp.find(ans)==mp.end())continue;
if(i*y-mp[ans]<0)continue;
return i*y-mp[ans];
}
return -1;
}
int main()
{
int i,j;
LL y,z,p,ans;
cin>>t>>k;
while(t--){
scanf("%lld%lld%lld",&y,&z,&p);
y%=p;
if(k==1){printf("%lld\n",qpow(y,z,p));}
else if(k==2){
z%=p;
if(y==0&&z){cout<<"Orz, I cannot find x!\n";continue;}
printf("%lld\n",z*inv(y,p)%p);
}
else{
z%=p;
ans=BSGS(y,z,p);
if(ans==-1){cout<<"Orz, I cannot find x!\n";continue;}
printf("%lld\n",ans);
}
}
return 0;
}
扩展
此时\(c\)不一定是质数,无法保证解的循环性。因此试图将该问题转化为普通BSGS问题。
首先给出一个判定解是否存在的方法:
当\((a,c)\not|\ \ b \quad \text{and}\quad b\not=1\)时,方程无自然数解。
因此,当\((a,c)|b\)时,不妨设\(G=(a,c)\)。
\]
\]
令\(x'=x-1,b'=\frac bG(\frac aG)^{-1},c'=\frac cG\),则有
\]
递归进行以上过程,直至\((a,c)=1\),使用普通BSGS即可。
例题
Luogu4195 【模板】exBSGS/Spoj3105 Mod
代码
#include<bits/stdc++.h>
#include<tr1/unordered_map>
#define LL long long
using namespace std;
using namespace std::tr1;
unordered_map<int,int>mp;
#define _BUFFERSIZE 65536
static size_t _pos = _BUFFERSIZE, _len; static char _buf[_BUFFERSIZE];
inline void _getc(char &c){
if(_pos==_BUFFERSIZE){
_pos=0;
_len=fread(_buf,1,_BUFFERSIZE,stdin);
}
c=_pos<_len?_buf[_pos++]:0;
}
template<typename T>
void read(T& x){
x=0;
char c;
do _getc(c); while(!isdigit(c));
do{
x=x*10+(c-'0');
_getc(c);
}while(isdigit(c));
}
void exgcd(LL a,LL b,LL &x,LL &y){
if(!b){x=1;y=0;return;}
exgcd(b,a%b,x,y);
LL tmp=x;x=y;y=tmp-a/b*y;
}
LL inv(LL a,LL p){
LL x,y;
exgcd(a,p,x,y);
return (x%p+p)%p;
}
LL gcd(LL a,LL b){return b?gcd(b,a%b):a;}
LL BSGS(LL a,LL b,LL p){
if(p==1&&b==0)return 0;
if(b==1)return 0;
if(a==0||b==0){
if(a==0&&b==0)return 1;
return -1;
}
LL i,y=(LL)sqrt(p),tmp,ans;
mp.clear();
for(i=0,tmp=1;i<y;i++,tmp=tmp*a%p){
mp[b*tmp%p]=i;
}
for(i=1,ans=tmp;y*(i-1)+1<=p-2;i++,ans=ans*tmp%p){
if(mp.find(ans)==mp.end())continue;
if(i*y-mp[ans]<0)continue;
return i*y-mp[ans];
}
return -1;
}
LL exBSGS(LL a,LL b,LL p){
a%=p;b%=p;
if(b==1)return 0;
if(a==0){if(b==0){return 1;}return -1;}
if(b==0){
LL ans=0,d;
while((d=gcd(a,p))!=1){
ans++;p/=d;
if(p==1)return ans;
}
return -1;
}
LL ans,d,cnt=0;
while((d=gcd(a,p))!=1){
if(b%d)return -1;
cnt++;
p/=d;b=b/d*inv(a/d,p)%p;
}
ans=BSGS(a,b,p);
if(ans==-1)return -1;
return cnt+ans;
}
int main()
{
int i,j;
LL a,p,b,ans;
for(;;){
read(a);read(p);read(b);
if(!(a||b||p))break;
ans=exBSGS(a,b,p);
if(ans==-1){printf("No Solution\n");}
else{printf("%lld\n",ans);}
}
return 0;
}
BSGS及其扩展的更多相关文章
- BSGS与扩展BSGS
BSGS \(BSGS\)算法又称大步小步\((Baby-Step-Giant-Step)\)算法 \(BSGS\)算法主要用于解以下同余方程 \[A^x\equiv B(mod\ p)\]其中\(( ...
- BSGS及扩展BSGS总结(BSGS,map)
蒟蒻哪里有什么总结,只能点击%YL% 还有这位ZigZagK大佬的blog \(\mbox{BSGS}\) 模板题:洛谷P3846 [TJOI2007]可爱的质数 给定\(a,b\)和模数\(\mbo ...
- BSGS和扩展BSGS
BSGS: 求合法的\(x\)使得\(a ^ x \quad mod \quad p = b\) 先暴力预处理出\(a^0,a^1,a^2.....a^{\sqrt{p}}\) 然后把这些都存在map ...
- BSGS及扩展BSGS算法及例题
\(BSGS(baby-step-giant-step)\)算法是用来解高次同余方程的最小非负整数解的算法,即形如这个的方程: \(a^x\equiv b(mod\ p)\) 其中\(p\)为质数(其 ...
- Codeforces 1106F Lunar New Year and a Recursive Sequence (数学、线性代数、线性递推、数论、BSGS、扩展欧几里得算法)
哎呀大水题..我写了一个多小时..好没救啊.. 数论板子X合一? 注意: 本文中变量名称区分大小写. 题意: 给一个\(n\)阶递推序列\(f_k=\prod^{n}_{i=1} f_{k-i}b_i ...
- BSGS 和扩展
BSGS BSGS,全称叫 BabyStepGiantStep,也就是大步小步 其实还是比较暴力的 它可以\(O(\sqrt p)\)的复杂度内解出: \[a^x\equiv n\pmod p,\gc ...
- 数论算法 剩余系相关 学习笔记 (基础回顾,(ex)CRT,(ex)lucas,(ex)BSGS,原根与指标入门,高次剩余,Miller_Rabin+Pollard_Rho)
注:转载本文须标明出处. 原文链接https://www.cnblogs.com/zhouzhendong/p/Number-theory.html 数论算法 剩余系相关 学习笔记 (基础回顾,(ex ...
- BSGS算法学习笔记
从这里开始 离散对数和BSGS算法 扩展BSGS算法 离散对数和BSGS算法 设$x$是最小的非负整数使得$a^{x}\equiv b\ \ \ \pmod{m}$,则$x$是$b$以$a$为底的离散 ...
- Discrete Logging(POJ2417 + BSGS)
题目链接:http://poj.org/problem?id=2417 题目: 题意: 求一个最小的x满足a^x==b(mod p),p为质数. 思路: BSGS板子题,推荐一篇好的BSGS和扩展BS ...
随机推荐
- Python在项目外更改项目内引用
前言 目前有一个奇葩的需求, 将某个开源项目整合进自己的项目里去调度, 还需要在每次启动这个开源项目时, 加载不同的配置文件进去, 问题是配置文件并不是一个 conf 或者是其他的什么, 而是以 .p ...
- 一文读懂 SuperEdge 边缘容器架构与原理
前言 superedge是腾讯推出的Kubernetes-native边缘计算管理框架.相比openyurt以及kubeedge,superedge除了具备Kubernetes零侵入以及边缘自治特性, ...
- Java远程下载文件到本地(http协议和ssh2协议)
Java中java.io包为我们提供了输入流和输出流,对文件的读写基本上都依赖于这些封装好的关于流的类中来实现.前段时间遇到了以下两种需求: 1.与某系统对接,每天获取最新的图片并显示在前端页面.该系 ...
- XEE - Pikachu
概述 XXE -"xml external entity injection"既"xml外部实体注入漏洞".概括一下就是"攻击者通过向服务器注入指定的 ...
- buuctf—web—Easy Calc
启动靶机,查看网页源码,发现关键字 $("#content").val() 是什么意思: 获取id为content的HTML标签元素的值,是JQuery, ("# ...
- 【易筋经】Llinux服务器初始化及常用命令大全
Llinux服务器初始化及常用命令大全 1.关闭防火墙以及内核安全机制 systemctl stop firewalld systemctl disable firewalld ##永久性关闭 set ...
- 查看内核打印信息指令dmesg
linux系统启动的时候打印的的信息非常重要,有时候需要看这些信息但是又不想重启,可以用dmesg这条指令.
- Win2008 server R2重置登录密码Administrator
1.PE方式修改密码 背景:https://www.cnblogs.com/Crazy-Liu/p/11245730.html 上述连接中的有AD域的机器系统使用哑巴式老毛桃等启动PE出现以下: 原因 ...
- jmeter-登录获取cookie后参数化,或手动添加cookie, 再进行并发测试
以下情况其实并不适用于直接登录可以获取cookie情况,直接可以登录成功,直接添加cookie管理,cookie可以直接使用用于以下请求操作. 如果登录一次后,后续许多操作,可以将cookie管理器放 ...
- java.io.IOException: Could not find resource com/xxx/xxxMapper.xml
java.io.IOException: Could not find resource com/xxx/xxxMapper.xml 报错内容: org.apache.ibatis.exception ...