BZOJ3239: Discrete Logging

最新推荐文章于 2018-06-05 15:34:21 发布

forezxl

最新推荐文章于 2018-06-05 15:34:21 发布

阅读量242

点赞数 1

CC 4.0 BY-SA版权

分类专栏： BZOJ 数论---(EX)BSGS 蒟蒻zxl的Blog专栏文章标签： BZOJ BSGS

本文链接：https://blog.youkuaiyun.com/a1799342217/article/details/79373688

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章

订阅专栏

BZOJ

281 篇文章

订阅专栏

数论---(EX)BSGS

3 篇文章

订阅专栏

本文介绍了一道典型的BSGS算法题目，并给出了详细的代码实现。通过使用BSGS算法解决特定的数论问题，文章提供了一个可以快速查找离散对数的解决方案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

BSGS

题目传送门

BSGS裸题，找了hzw的当了个板子。
~~太懒没打哈希表~~
（缺氧的POJ用map会T）

代码：

#include<map>
#include<cmath>
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int m,a,b;
map <int,int> mp;
typedef long long LL;
inline LL ksm(LL a,LL b){
    LL ret=1; a%=m;
    while (b){
        if (b&1) ret=ret*a%m;
        a=a*a%m,b>>=1;
    }
    return ret;
}
int main(){
    while (scanf("%d%d%d",&m,&a,&b)!=EOF){
        a%=m,mp.clear();
        if ((!a)&&(!b)) { puts("1"); continue; }
        if (!a) { puts("no solution"); continue; }
        int c=ceil(sqrt(m)); LL x=1; bool flag=false; mp[1]=c+1;
        for (int i=1;i<c;i++){  x=x*a%m; if (!mp[x]) mp[x]=i; }
        for (LL i=0,n=1,tem=ksm(a,m-c-1);i<c;i++){
            int l=mp[b*n%m];
            if (l) {
                if (l==c+1) l=0; 
                printf("%lld\n",i*c+l);
                flag=true; break;
            }
            n=n*tem%m;
        }
        if (!flag) puts("no solution");
    }
    return 0;
}