BSGS 学习笔记

问题：求$a^x\equiv b\ (mod\ p)$的最小正整数解。

这时候就要用到BSGS（拔山盖世）算法。直接进入正题：

设$x=im-n$，

则原式等于$a^{im-n}\equiv b\ (mod\ p)$。

移项，得$a^{im}\equiv a^nb(mod\ p)$。

我们把所有$a^nb$的状态存到一个map里，然后枚举$a^{im}$，如果相等则找到最小正整数解。

当$m=\sqrt p$时，算法效率最高。则$[1,m]$枚举$n$,$[1,m]$枚举$i$。

以上说的情况是$a$与$p$互质的情况。那么不互质该怎么做呢？

我们变换一下形式：$a*a^{x-1}\equiv b\ (mod\ p)$。

移项，得$a*a^{x-1}+y*p=b$。

设$g=gcd(a,p)$，由裴蜀定理得，如果$b\ mod\ p≠0$，那么此同余方程无解。

左右两边同除$g$，得到$\frac{a}{g}*a^{x-1}+y*\frac{p}{g}=\frac{b}{g}$，即$a^{x-1}\equiv \frac{b}{g}*(\frac{a}{g})^{-1}\ (mod\ \frac{p}{g})$

重复上述步骤，直到$gcd(a,p)=1$为止，然后就可以用普通的BSGS求啦。

注意要特判一下，如果$b=1$或$p=1$或者某一时刻$(\frac{a}{g})^x=b'$，那么直接返回$0$即可。

注意各种情况下的模数！！！被这个坑了好久QAQ。

代码：

#include<bits/stdc++.h>
#define int long long
using namespace std;
map<int,int> vis;
int a,p,b;
inline int gcd(int a,int b)
{
    if(!b) return a;
    return gcd(b,a%b);
}
inline void exgcd(int a,int b,int &x,int &y)
{
    if (!b) x=,y=;
    else{
        exgcd(b,a%b,x,y);
        int t=x;x=y;y=t-a/b*y;
    }
}
inline int inv(int a,int b)
{
    int x,y;
    exgcd(a,b,x,y);
    return (x%b+b)%b;
}
inline int qcal(int a,int b,int p)
{
    int res=;
    while(b)
    {
        if (b&) res=res*a%p;
        a=a*a%p;
        b>>=;
    }
    return res;
}
inline int bsgs(int a,int b,int p)
{
    vis.clear();
    b%=p;
    int m=ceil(sqrt(p));
    for (int i=;i<=m;i++) b=b*a%p,vis[b]=i;
    b=;int tmp=qcal(a,m,p);
    for (int i=;i<=m;i++)
    {
        b=b*tmp%p;
        if (vis[b]) return (i*m-vis[b]+p)%p;
    }
    return -;
}
inline int exbsgs(int a,int b,int p)
{
    if (b==||p==) return ;
    int g=gcd(a,p),k=,na=;
    while(g>)
    {
        if (b%g!=) return -;
        b/=g;p/=g;na=na*(a/g)%p;
        k++;
        if (na==b) return k;
        g=gcd(a,p);
    }
    int f=bsgs(a,b*inv(na,p)%p,p);
    if (f==-) return -;
    return f+k;
}
signed main()
{
    cin>>a>>p>>b;
    while(a||p||b)
    {
        a%=p;b%=p;
        int t=exbsgs(a,b,p);
        if (t==-) puts("No Solution");
        else printf("%lld\n",t);
        cin>>a>>p>>b;
    }
    return ;
}