BSGS算法学习小记(大步小步算法)

BSGS 同时被 3 个专栏收录
2 篇文章 0 订阅
26 篇文章 0 订阅
75 篇文章 0 订阅

简介

先看一个式子 xyz(modp)z
现在只知道x和z,要求y。
大步小步算法(BSGS,Baby Steps Giant Steps)就是解决这个问题。

算法流程

暴搜的枚举范围

根据费马小定理: xp11
如果y已经枚举到了p-1了,继续枚举的话就会产生循环。
所以,在暴搜中y的枚举范围就是0……p-1。

如何优化暴搜

我们想一想可不可以用分块来解决枚举的y。
把y分成 p1 分别枚举行不行?
设m= p1 y=am+b ,这样枚举a和b就相当于分块枚举了。
那么现在就变成了 xam+bz(modp)
把a和b分别放在两边: xbxamz(modp)
我们可以发现左边的 xb 最多只有m个,完全可以预处理出来放进hash里面。

 ll m=ceil(sqrt(p-1));
    ll ni=qsm(x,m);ni=qsm(ni,p-2);//因为右边的指数是复数,所以需要逆元
    a.clear();
    a[1]=m+1;
    fo(j,1,m-1){
        k=k*x%p;
        if(!a[k])a[k]=j;
    }   

然后枚举右边的a就好了。

fo(i,0,m-1){
        ll u=a[y];
        if(u){
            if(u==m+1)printf("%lld\n",i*m);
            else printf("%lld\n",i*m+u);
            return;
        }
        y=y*ni%p;    
    }
    printf("Orz, I cannot find x!\n");

这就是 O(p) 的用空间换取时间的大步小步算法。
很优美的暴搜。

推荐题目

经典的裸题SDOI2011 计算器

由于本人是个蒟蒻

对于BSGS算法知道的,也只有这么多了。

  • 2
    点赞
  • 0
    评论
  • 1
    收藏
  • 打赏
    打赏
  • 扫一扫,分享海报

©️2022 CSDN 皮肤主题:编程工作室 设计师:CSDN官方博客 返回首页

打赏作者

Facico

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值