BSGS 学习笔记

glorious_dream

于 2022-09-04 09:51:06 发布

阅读量273

点赞数 3

分类专栏：模板总结，学习笔记文章标签：学习算法数据结构

本文链接：https://blog.csdn.net/glorious_dream/article/details/126685863

版权

模板总结，学习笔记专栏收录该内容

30 篇文章 3 订阅

订阅专栏

在这里插入图片描述
用 $BSGS$ 算法可以在 $O(\sqrt n)$ 的时间内求出上述的同余式。

设 $\left\lceil\sqrt p\right\rceil$ ，若有解，则最小非负整数解一定在 $[0, p]$ 内，所以设 $\times m - j$ 。

将式子变形， $a^{i \times m - j} \equiv b (\mod p) \Rightarrow a^{i \times m} \equiv b \times a^j(\mod p)$ 。

这样我们枚举 $j$ ，范围 $[0, m]$ ，将 $\times a^j$ 存入哈希表。

然后枚举 $i$ ，范围 $[1, m]$ ，从哈希表中寻找第一个满足 $a^{i \times m} \equiv b \times a^j(\mod p)$ ，此时 $\times m - j$ 即为所求的答案。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <map>
#include <limits.h>
#define re register
#define int long long
#define drep(a,b,c) for(re int a(b) ; a>=(c) ; --a)
#define rep(a,b,c) 	for(re int a(b) ; a<=(c) ; ++a)
using namespace std;
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch == '-') f=-1 ; ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	return x*f;
}
inline void print(int x){
	if(x < 0) putchar('-'),x = -x;
	if(x >= 10) print(x / 10);
	putchar(x % 10 + '0');
}
int p,b,n,fl,ans=LONG_LONG_MAX;
map<int,int> mp;
signed main(){
	p = read(),b = read(),n = read();
	int m = sqrt(p);
	int mul = 1;
	rep(i,0,m-1) mp[n*mul%p] = i,mul=mul*b%p;
	int base = 1;
	rep(i,0,m){
		if(mp[base]) fl = 1,ans = min(ans,(i*m-mp[base]+p)%p);
		base = base*mul%p;
	}
	if(fl) printf("%lld\n",ans);
	else printf("no solution");
	return 0;
}