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 ...
随机推荐
- C++语言基础——01一切的开始
环境准备 集成开发环境(Integrated Development Environment,IDE) 竞赛中最常见的是Dev-C++,平时练习采用其他IDE也可,使用方法都是类似的. 编译器 推荐使 ...
- Dubbo 就是靠它崭露头角!(身为开源框架很重要的一点)
Hola,我是 yes. 经过了 RPC 核心和 Dubbo 微内核两篇文章后,今天终于要稍稍深入一波 Dubbo 了. 作为一个通用的 RPC 框架,性能是很重要的一环,而易用性和扩展性也极为重要. ...
- Goland 设置代码格式化
前言 之前一直喜欢 VsCode 的代码自动格式化和其他的一些功能 今天了解到原来 Goland 也有这些功能, 想想也对, 毕竟这么大 正文 Goland设置代码格式化 进入设置,按需选择要使用的, ...
- 图片质量评估论文 | 无监督SER-FIQ | CVPR2020
文章转自:同作者微信公主号[机器学习炼丹术].欢迎交流,共同进步. 论文名称:SER-FIQ: Unsupervised Estimation of Face Image Quality Based ...
- Spring Cloud实战 | 第十篇 :Spring Cloud + Seata 1.4.1 + Nacos1.4.0 整合实现微服务架构中逃不掉的话题分布式事务
Seata分布式事务在线体验地址:https://www.youlai.store 本篇完整源码地址:https://github.com/hxrui/youlai-mall 有想加入开源项目开发的童 ...
- 【JavaWeb】AJAX 请求
AJAX 请求 什么是 AJAX AJAX(Asynchronous JavaScript And XMl),即异步 JS 和 XML.是指一种创建交互式网页应用的网页开发技术. AJAX 是一种浏览 ...
- 【Flutter】功能型组件之导航返回拦截
前言 为了避免用户误触返回按钮而导致APP退出,在很多APP中都拦截了用户点击返回键的按钮,然后进行一些防误触判断,比如当用户在某一个时间段内点击两次时,才会认为用户是要退出(而非误触).Flutte ...
- 关于软件架构中的b/s
**B/S架构 b/s只需要一个浏览器,用户就可以通过不同的网址访问不同的服务器程序. 优点:开发,安装,部署,维护简单 缺点:对硬件要求过高,用户的体验会受到影响 首先是资源分类:**可以分为静态资 ...
- Azure 存储简介
Azure Storage Account(存储账户)包含所有Azure Storage的数据对象,包括Blob.Data Lake Gen2,File.Queue.Disk和Table等服务,该St ...
- 深入理解Redis之简单动态字符串
目录 SDS SDS与C字符串的区别 SDS获取字符串长度复杂度为O(1),C字符串为O(N) SDS杜绝了缓存区溢出 减少修改字符串时带来的内存重分配次数 二进制安全 Redis没有直接使用C语言传 ...