bsgs 算法 求高次同余 方程 a^x %p = b 的形式

利用指数循环节,枚举i,达到根号n的复杂度

ll qpow(int a,int b,int p)
{
	ll ans=1;
	while(b)
	{
		if(b&1)ans=ans*a%p;
		a=a*a%p;
		b/=2;
	}
	return ans;
}
int bsgs(int a,int b,int p)// a^x \equiv  b (mod p)
{
	map<int,int>hash;hash.clear();
	b%=p;
	int t=(int)sqrt(p)+1;
	for(int j=0;j<t;j++)
	{
		int val=(ll)b*qpow(a,j,p)%p;//b*a^j
		hash[val]=j;
	}
	a=qpow(a,t,p);//a^t
	if(a==0)return b==0?1:-1;//0 % b =0  只有b是0才有解 
	for(int i=0;i<=t;i++)
	{
		ll val=qpow(a,i,p);
		int j= hash.find(val)==hash.end() ? -1 : hash[val];
		if(j>=0&&i*t-j>=0)return i*t-j;//最非负整数解 
	}
	return -1;
}

 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值