BSGS模版（求解离散对数，模数p是不是质数都可以）

最新推荐文章于 2021-04-01 22:33:51 发布

zhhx2001

最新推荐文章于 2021-04-01 22:33:51 发布

阅读量976

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/zhhx2001/article/details/51852982

版权

数论专栏收录该内容

42 篇文章 0 订阅

订阅专栏

抄的邝斌的模版，感觉写的真心简洁，并且，这个hash判重的方法也很好

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;

ll b,n,p;
/***************************************/下面就是BSGS的函数集，包括所需变量
const int N=100009;
const int mod=76543;
ll hs[N],id[N],head[N],next[N],tot;
void insert(ll x,ll y)
{
	int k=x%mod;
	hs[++tot]=x;id[tot]=y;next[tot]=head[k];head[k]=tot;
}
ll find(ll x)
{
	int k=x%mod;
	for (int i=head[k];i;i=next[i])
	if (hs[i]==x) return id[i];
	return -1;
}
ll BSGS(ll a,ll b,ll n)
{
	if (b==1) return 0;
	memset(head,0,sizeof(head));
	tot=0;
	ll j,m=sqrt(n*1.0),x=1,p=1;
	for (int i=0;i<m;i++,p=p*a%n) insert(p*b%n,i);
	for (ll i=m;i<=n;i+=m)
		if ((j=find(x=x*p%n))!=-1) return i-j;
	return -1;
}
/**************************************/
int main()
{
	while(scanf("%lld%lld%lld",&p,&b,&n)!=EOF)//poj提交一定要有“！=EOF”      ！！！！！！！！！
	{
		ll ans=BSGS(b,n,p);
		if (ans==-1) printf("no solution\n");else printf("%lld\n",ans);
	} 
	return 0;
}