BSGS

离散对数问题

有一类问题形式如下:

给定同余方程 a x ≡ b ( m o d    m ) a^x\equiv b\quad(\mod m) axb(modm) ,求其最小解。

这样的方程解起来并不简单,最暴力的想法自然是一个一个试,于是对于这类“不可做”问题考虑分块的做法。下面我们将分两种情况讨论:

( a , m ) = 1 (a,m)=1 (a,m)=1

a a a m m m 互质的情况下问题就比较好处理了。

x = A m + B x=A\sqrt m+B x=Am +B

则原来方程可化为:
a A m + B ≡ b ( m o d    m )    ⟹    a B ≡ b a − A m ( m o d    m ) \begin{aligned} a^{A\sqrt m +B}&\equiv b\quad(\mod m)\\ \implies a^B&\equiv ba^{-A\sqrt m}\quad(\mod m) \end{aligned} aAm +BaBb(modm)baAm (modm)
于是我们把 a B a^B aB 存入哈希表中,然后 O ( n ) O(\sqrt n) O(n ) 找到 A A A 的取值。

Discrete Logging

#include
#include
#include
#include
using std::ceil;
using std::sqrt;
typedef long long ll;
int qpow(int a,int x,int p)
{
	int res=1;
	for(;x;x>>=1,a=(ll)a*a%p)
		if(x&1) res=(ll)res*a%p;
	return res;
}

std::map mp;
int BSGS(int a,int b,int p)
{
	mp.clear();
	int i,sp=ceil(sqrt(p));
	int aB=1;mp[aB]=0;
	int aspInv=qpow(qpow(a,sp,p),p-2,p);
	for(i=1;i

[SDOI2011] 计算器

#include
#include
#include
#include
using std::ceil;
using std::sqrt;
typedef long long ll;

int qpow(int a,int x,int p)
{
	int res=1;
	for(;x;x>>=1,a=(ll)a*a%p)
		if(x&1) res=(ll)res*a%p;
	return res;
}

std::map mp;
int BSGS(int a,int b,int p)
{
	if(!a) return b?-1:1;
	mp.clear();
	int i,sp=ceil(sqrt(p));
	int aB=1;mp[aB]=0;
	int aspInv=qpow(qpow(a,sp,p),p-2,p);
	for(i=1;i

( a , m ) ≠ 1 (a,m)\neq1 (a,m)=1

我们既然已经有了 a a a m m m 互质情况下的做法,那只需要将 a a a m m m 化为互质的即可。

对于 a x ≡ b ( m o d    m ) a^x\equiv b\quad(\mod m) axb(modm) ( a , m ) = g (a,m)=g (a,m)=g

如果 g ∤ b g\not\mid b gb 那么无解。

如果 g ∣ b g\mid b gb 那么原方程可化为 ( a ′ g ) x ≡ b ′ g ( m o d    m ′ g ) (a'g)^x\equiv b'g\quad(\mod m'g) (ag)xbg(modmg)

同余式两边同除以 g g g 得: ( a ′ ) x − 1 ≡ b ′ ( a ′ ) − 1 ( m o d    m ′ ) (a')^{x-1}\equiv b'(a')^{-1}\quad(\mod m') (a)x1b(a)1(modm)

注意每次新的 x x x 都是原来的 x − 1 x-1 x1 ,故需要记录操作次数 c n t cnt cnt

如果 b ′ ( a ′ ) − 1 ≡ 1 ( m o d    m ′ ) b'(a')^{-1}\equiv1\quad(\mod m') b(a)11(modm) m ′ = 1 m'=1 m=1 则解为 c n t cnt cnt

不断执行上述操作,直至 ( a , m ) = 1 (a,m)=1 (a,m)=1

注意特判 a ≡ 0 ( m o d    m ) a\equiv0\quad(\mod m) a0(modm) 时的情况

[SPOJ] Mod

#include
#include
#include
#include
using std::ceil;
using std::sqrt;
typedef long long ll;

int gcd(int a,int b){return b?gcd(b,a%b):a;}
void exgcd(int a,int b,int &x,int &y)
{
	if(!b) {x=1,y=0;return;}
	exgcd(b,a%b,y,x);y-=x*(a/b);
}
int inv(int a,int p){int x,y;exgcd(a,p,x,y);return (x%p+p)%p;}

int qpow(int a,int x,int p)
{
	int res=1;
	for(;x;x>>=1,a=(ll)a*a%p)
		if(x&1) res=(ll)res*a%p;
	return res;
}

std::map mp;
int BSGS(int a,int b,int p)
{
	a%=p,b%=p;
	if(!a) return b?-1:1;
	int cnt=0;
	for(int g=gcd(a,p);g>1;g=gcd(a,p))
	{
		if(b%g) return -1;
		b/=g,p/=g;b=(ll)b*inv(a/g,p)%p;
		cnt++;
		if(b==1) return cnt;
	}
	if(p==1) return cnt;

	mp.clear();
	int i,sp=ceil(sqrt(p));
	int aB=1;mp[aB]=0;
	int aspInv=inv(qpow(a,sp,p),p);
	for(i=1;i

你可能感兴趣的:(BSGS)