Shank大步小步算法

最新推荐文章于 2022-11-15 21:18:43 发布

zyjhtutu

最新推荐文章于 2022-11-15 21:18:43 发布

阅读量2k

点赞数

分类专栏：常用模板文章标签： acm 数论 Shank大步小步算法 poj2417

本文链接：https://blog.csdn.net/zyjhtutu/article/details/38272299

版权

Shank大步小步算法用于解决模方程a^x=b(mod N)，尤其当N为素数时。利用欧拉定理，只需检查0到n-1的x值。通过预先计算前sqrt(N)个可能解并存储，再处理剩余部分，避免超时，通常使用map或hash。在poj2417题目中，这是典型应用。

摘要由CSDN通过智能技术生成

Shank大步小步算法用来求解模方程a^x=b(mod N)，其中N为素数。有欧拉定理可得a^(n-1) = 1(mod N)（因为n为素数）,所以只需要检查x = 0,1,2,,,,n-1是不是即可。当x超过n-1时，a^x就开始循环了。这里用大步小步算法来检查出x的解。先求出前m个可能取值(0,1,2,3,,m-1)（其中m常取sqrt(N)），将这前m项对N取余的解保存起来(包括得到该解的最小的x值)，为了判重和便于后来查找，通常用STL中map<int,int>。但有时候数据量太大，用map就会超时，这是一般用hash的策略或者用二分。对于大于m的部分可以通过对模方程进行变形，再到前面保存的结果里面去查询是否为可行解。其时间复杂度为(n^(1/2)*log(n)).

大步小步算法解a^x = b (mod N)（其中N为素数）的模板题目poj2417;

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<climits>
#include<cctype>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
#include<set>
#include<string>
#include<map>
#include<stack>
#define ll long long
#define MAX 100010
#define MOD 100007
#define eps 1e-8
#define INF INT_MAX

using namespace st