数学模板-BSGS&&EXBSGS

最新推荐文章于 2021-07-21 16:24:07 发布

YJSchaf

最新推荐文章于 2021-07-21 16:24:07 发布

阅读量2.3k

点赞数

分类专栏：数学模板文章标签：模板数学 BSGS

本文链接：https://blog.csdn.net/YJSchaf/article/details/68942904

版权

数学同时被 2 个专栏收录

11 篇文章 0 订阅

订阅专栏

模板

6 篇文章 0 订阅

订阅专栏

BSGS

BSGS(baby step giant step)就是求同余方程

a x \equiv b (mod p)

$a^x≡b(\mod p)$
我们设

x=im−j $x = im - j$ ,将其带入式子可以得到

aim−j≡b(modp) $a^{im - j}≡b(\mod p)$
即

aim≡b∗aj(modp) $a^{im}≡b*a^{j}(\mod p)$
我们可以计算

b∗ajmodp $b*a^{j}\mod p$ (i=0 to m)的值(快速幂) 将其放入map中

map[b∗ajmodp]=j $map[b*a^{j}\mod p]=j$
之后再枚举a的im次mod p的值，如果发现再map有对应解

ans=m∗i−map[aimmodp] $ans=m*i-map[a^{im}\mod p]$
可以证明，当

m=p‾‾√ $m= \sqrt p$ 最快。
为什么

x=im−j $x = im - j$ ，即x的最大值不超过p.
有公式即

ak%p−1≡ak(modp) $a^{k\%p-1}≡a^{k}(\mod p)$ (当a与p互质),通过小费马证得。

k $k%(p-1)$ 可以看成

k−m(p−1) $k-m(p-1)$ 于是上面的式子变成了

ak(a(p−1))m≡ak(modp) ${{a^k}\over {(a^(p-1))^m}}≡a^k(\mod p)$
由于小费马，

ap−1≡1(modp) $a^{p-1}≡1(\mod p)$ ,可以得到

ak1m≡ak(modp) ${a^{k}\over 1^m}≡a^k(\mod p)$
这样就确保了x比p小。

LL BSGS(LL y,LL z,LL p)
{
    map<LL,LL> ma;
    LL m=sqrt(p),tmp=0;
    ma.clear();
    if(y%p==0&&z==0) return 1;
    if(y%p==0&&z!=0) return -1;
    for(int i=0;i<=m;i++)
    {
        if(!i) {tmp=z%p;ma[tmp]=i;continue;}
        tmp=(tmp*y)%p;
        ma[tmp]=i;
    }
    tmp=1;LL t=power(y,m,p);
    for(int i=1;i*i<=p;i++)
    {
        tmp=(tmp*t)%p;
        if(ma[tmp])
        {
            LL ans=i*m-ma[tmp];
            return ans;
        }
    }
    return -1;
}