BSGS模版(求解离散对数,模数p是不是质数都可以)

抄的邝斌的模版,感觉写的真心简洁,并且,这个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;
}


  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值