蒟蒻哪里有什么总结,只能点击%YL%

还有这位ZigZagK大佬的blog

\(\mbox{BSGS}\)

模板题:洛谷P3846 [TJOI2007]可爱的质数

给定\(a,b\)和模数\(\mbox{YL}\),求\(a^x\equiv b(\bmod\mbox{YL})\)中\(x\)的最小非负整数解。保证\(\gcd(a,\mbox{YL})=1\)。

设\(k=\lceil\sqrt{\mbox{YL}}\rceil\),令\(x=ky-c\)(\(y\in[1,k],c\in[0,k)\),故分解唯一)

于是有\(a^{ky}\equiv ba^c(\bmod\mbox{YL})\)

两边的取值都只有\(k\)种,枚举右边丢进map,再枚举左边查找一下即可得到最小的解。

复杂度\(O(\sqrt{\mbox{YL}}\log\sqrt{\mbox{YL}})\)。注意特判掉一些特殊情况。

#include<bits/stdc++.h>
#define LL long long
#define RG register
#define R RG int
using namespace std;
map<int,int>c;
inline LL qpow(RG LL b,R k,RG LL YL){//快速幂
RG LL a=1;
for(;k;k>>=1,b=b*b%YL)
if(k&1)a=a*b%YL;
return a;
}
int main(){
R a,b,YL;
cin>>YL>>a>>b;
if((b%=YL)==1)return puts("0"),0;//特判×2
if((a%=YL)==0)return puts("no solution"),0;
R k=sqrt(YL)+1,y;
RG LL p,pw=qpow(a,k,YL);
for(p=b,y=0;y<k;++y,(p*=a)%=YL)//预处理
c[(int)p]=y;
for(p=pw,y=1;y<=k;++y,(p*=pw)%=YL)//查找
if(c.count((int)p))return(cout<<(LL)k*y-c[(int)p]<<endl),0;
return puts("no solution"),0;//找不到
}

扩展\(\mbox{BSGS}\)

模板题:洛谷P4195 Spoj3105 Mod

双倍经验:洛谷SP3105 MOD - Power Modulo Inverted

可以解决\(\gcd(a,\mbox{YL})\neq1\)的情况,核心思想是试探性地约去\(\gcd\)。

同余方程\(a^x\equiv b(\bmod\mbox{YL})\),可以提一个\(a\)出来改成不定方程\(a\cdot a^{x-1}+\mbox{YL}\cdot y=b\),等于说暂时把\(a^{x-1}\)看成了未知数。

设\(g=\gcd(a,\mbox{YL})\),学完exgcd之后我们就知道了如果\(g\nmid b\)则方程组无解。否则就可以约成\(\frac a g\cdot a^{x-1}+\frac{\mbox{YL}}g\cdot y=\frac b g\)

我们把约简后的方程看作一个新的不定方程\(a'\cdot a^{x-1}+\mbox{YL}'\cdot y=b'\),如果\(\gcd(a,\mbox{YL}')\)还是不为\(1\)的话,我们就再从幂里面拿出一个\(a\)和\(\mbox{YL}'\)约。如此循环。

终于在约了\(cnt\)次后,我们得到了一个\(\gcd=1\)的方程了,变回来就是\(a'\cdot a^{x-cnt}\equiv b'(\bmod\mbox{YL}')\)

接下来的解法已经属于常规\(\mbox{BSGS}\)了,不过因为一些细节上的区别还是再写一遍。

仍然设\(k=\lceil\sqrt{\mbox{YL}'}\rceil\),令\(x=ky-c\)(\(y\in[1,k],c\in[0,k)\))

于是有\(a'\cdot a^{ky}\equiv ba^{cnt}\cdot a^c(\bmod\mbox{YL}')\),其实只是恒等号两边多乘了点东西而已。

约分多出来的复杂度是\(O(\log^2\mbox{YL})\),因为最多被约\(\log\)次,它的阶小于根号就忽略掉了。又要注意一些特判。

#include<bits/stdc++.h>
#define LL long long
#define RG register
#define R RG int
using namespace std;
map<int,int>c;
inline LL qpow(RG LL b,R k,RG LL YL){
RG LL a=1;
for(;k;k>>=1,b=b*b%YL)
if(k&1)a=a*b%YL;
return a;
}
int gcd(R x,R y){
return y?gcd(y,x%y):x;
}
int main(){
R a,b;RG LL YL;
cin>>a>>YL>>b;
while(a||b||YL){
R x=-1,a1=1,cnt=0,g,k,y,p,pw;
a%=YL;
if((b%=YL)==1){x=0;goto L;}//这个特判上面提过
while((g=gcd(a,YL))!=1){
if(b%g)goto L;//无解充要条件
++cnt;b/=g;YL/=g;
a1=a/g*(LL)a1%YL;
if(a1==b){x=cnt;goto L;}
}//注意这个特判,此时等价于a^{x-cnt}≡1(mod YL)
c.clear();//多组数据注意清空
k=sqrt(YL)+1;//常规BSGS过程
pw=qpow(a,k,YL);
p=b*qpow(a,cnt,YL)%YL;
for(y=0;y<k;++y,p=(LL)p*a%YL)
c[p]=y;
p=(LL)a1*pw%YL;//多乘一个a1
for(y=1;y<=k;++y,p=(LL)p*pw%YL)
if(c.count(p)){x=k*y-c[p];break;}
L:~x?printf("%d\n",x):puts("No Solution");
cin>>a>>YL>>b;
}
return 0;
}

BSGS及扩展BSGS总结(BSGS,map)的更多相关文章

  1. BSGS与扩展BSGS

    BSGS \(BSGS\)算法又称大步小步\((Baby-Step-Giant-Step)\)算法 \(BSGS\)算法主要用于解以下同余方程 \[A^x\equiv B(mod\ p)\]其中\(( ...

  2. BSGS和扩展BSGS

    BSGS: 求合法的\(x\)使得\(a ^ x \quad mod \quad p = b\) 先暴力预处理出\(a^0,a^1,a^2.....a^{\sqrt{p}}\) 然后把这些都存在map ...

  3. BSGS及扩展BSGS算法及例题

    \(BSGS(baby-step-giant-step)\)算法是用来解高次同余方程的最小非负整数解的算法,即形如这个的方程: \(a^x\equiv b(mod\ p)\) 其中\(p\)为质数(其 ...

  4. Codeforces 1106F Lunar New Year and a Recursive Sequence (数学、线性代数、线性递推、数论、BSGS、扩展欧几里得算法)

    哎呀大水题..我写了一个多小时..好没救啊.. 数论板子X合一? 注意: 本文中变量名称区分大小写. 题意: 给一个\(n\)阶递推序列\(f_k=\prod^{n}_{i=1} f_{k-i}b_i ...

  5. BSGS 和扩展

    BSGS BSGS,全称叫 BabyStepGiantStep,也就是大步小步 其实还是比较暴力的 它可以\(O(\sqrt p)\)的复杂度内解出: \[a^x\equiv n\pmod p,\gc ...

  6. BSGS及其扩展

    目录 定义 原理 朴素算法 数论分块 例题 Luogu2485 [SDOI2011]计算器 题解 代码 扩展 例题 Luogu4195 [模板]exBSGS/Spoj3105 Mod 代码 之前写了一 ...

  7. python3-知识扩展扫盲易忘-map,collections.Counter()的用法

    map() 会根据提供的函数对指定序列做映射. 第一个参数 function 以参数序列中的每一个元素调用 function 函数,返回包含每次 function 函数返回值的新列表. >> ...

  8. 数论算法 剩余系相关 学习笔记 (基础回顾,(ex)CRT,(ex)lucas,(ex)BSGS,原根与指标入门,高次剩余,Miller_Rabin+Pollard_Rho)

    注:转载本文须标明出处. 原文链接https://www.cnblogs.com/zhouzhendong/p/Number-theory.html 数论算法 剩余系相关 学习笔记 (基础回顾,(ex ...

  9. BSGS算法学习笔记

    从这里开始 离散对数和BSGS算法 扩展BSGS算法 离散对数和BSGS算法 设$x$是最小的非负整数使得$a^{x}\equiv b\ \ \ \pmod{m}$,则$x$是$b$以$a$为底的离散 ...

随机推荐

  1. Day7 Ubantu学习(一)

    Linux是多用户操作系统 Ubantu学习参考网址:https://www.cnblogs.com/resn/p/5800922.html 1.虚拟机网络类型的理解 bridged(桥接模式) : ...

  2. scrapy之日志等级

    scrapy之日志等级 在settings.py中配置如下项: LOG_LEVEL = 'ERROR' # 当LOG_LEVEL设置为ERROR时,在进行日志打印时,只是打印ERROR级别的日志 这样 ...

  3. Linux之查看系统配置命令

    1.查看CPU个数: [root@ifusion ~]# lscpu 2.查看CPU进程数.核数: [root@ifusion ~]# cat /proc/cpuinfo 3.查看内存总量: [roo ...

  4. 牛客OI周赛8-普及组

    https://ac.nowcoder.com/acm/contest/543#question A. 代码: #include <bits/stdc++.h> using namespa ...

  5. http/https与websocket的ws/wss的关系以及通过Nginx的配置

    http/https与websocket的ws/wss的关系 - 哒哒哒 - CSDN博客 https://blog.csdn.net/Garrettzxd/article/details/81674 ...

  6. Mysql之常用操作(2)

    Windows服务 -- 启动MySQL net start mysql -- 创建Windows服务 sc create mysql binPath= mysqld_bin_path(注意:等号与值 ...

  7. Git命令以及常见注意事项

    命令: git init -> 初始化一个git仓库 git clone -> 克隆一个本地库 git pull -> 拉取服务器最新代码 git fetch –p -> 强行 ...

  8. 【360图书馆】插入U盘自动攻击:BadUSB原理与实现

    插入U盘自动攻击:BadUSB原理与实现       漏洞背景 “BadUSB”是今年计算机安全领域的热门话题之一,该漏洞由Karsten Nohl和Jakob Lell共同发现,并在今年的Black ...

  9. Your branch is ahead of 'origin/master' by 2 commits.

    遇到这种问题,表示在你之前已经有2个commit而没有push到远程分支上,所以需要先git push origin **将本地分支提到远程仓库.也可以直接git reset --hard HEAD~ ...

  10. python数据结构与算法第三天【时间复杂度计算方法】

    最优时间复杂度(不可靠) 最坏时间复杂度(保证) 平均时间复杂度(平均状况) 不同语句的时间复杂度: (1)顺序语句:使用加法 (2)循环语句:使用乘法 (3)分支语句:使用坏时间复杂度 例如:如下代 ...