学习笔记——BSGS

Young_20220202

已于 2023-10-27 20:45:51 修改

阅读量212

点赞数

文章标签：学习笔记哈希算法

于 2023-09-29 13:09:58 首次发布

本文链接：https://blog.csdn.net/Young_20220202/article/details/133416274

版权

众所周知，北上广深是中国非常一线的城市，北京是首都，地处……

正片开始！

一、BSGS基础算法

实现目标： $A^x\equiv B(\mod P),(\gcd(P,A)=1)$ 求最小的 $x$
很明显，如果暴力枚举，时间是 $O (P)$ 的，只要题目数据范围大，就死定了。愿意的人欢迎尝试（无100警告）
于是，考虑莫队
分块
BSGS
为什么我要提前两个？
因为BSGS本质就是分块！
简单讲解一下，就是将本来 $P$ 种情况，平均分成了 $\sqrt p$ 份，对每份内进行预处理
不直观？好，那就用直观的式子吧！
$\texttt{令}x=k\sqrt p-t$
$\texttt{即}A^{k\sqrt p}\equiv A^tB (\mod p)$

于是，找个哈希表维护一下后面的即可

IO::pin>>y>>z>>p;
gp_hash_table<int,int> ht;
int f,s;
bool flg;
ht.clear();
f=ceil(sqrt(p));
s=1;
flg=1;
for(int i=1; i<=f; i++)
    s=1ll*s*y%p,ht[1ll*z*s%p]=i;
for(int j=1,k=s; j<=f; j++) {
	if(ht[k]&&flg) {
    	wt(((1ll*j*f-ht[k])%p+p)%p,'\n');
    	flg=0;
	}
			
	k=(1ll*s*k)%p;
	}
	if(flg)
		puts("Orz, I cannot find x!");

附赠模板代码

#pragma GCC optimize(1,"inline","Ofast")
#pragma GCC optimize(2,"inline","Ofast")
#pragma GCC optimize(3,"inline","Ofast")
#include<bits/stdc++.h>
#include<bits/extc++.h>
using namespace std;
using namespace __gnu_cxx;
using namespace __gnu_pbds;
namespace IO {
	class input {
	private:
		bool isdigit(char c) {
			return ('0'<=c&&c<='9');
		}
	public:
		input operator>>(int &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(short &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(bool &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(long &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(long long &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(__int128 &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(unsigned int &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(unsigned short &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(unsigned long &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(unsigned long long &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(unsigned __int128 &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=(x<<1)+(x<<3)+(c^48),c=getchar();
			if(!y)
				x=-x;
			return *this;
		}
		input operator>>(double &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=x*10+(c^48),c=getchar();
			if(!y)
				x=-x;
			if(!isdigit(c))
				if(c!='.')
					return *this;
			double z=1;
			while(isdigit(c))
				z/=10.,x=x+z*(c^48),getchar();
			return *this;
		}
		input operator>>(long double &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=x*10+(c^48),c=getchar();
			if(!y)
				x=-x;
			if(!isdigit(c))
				if(c!='.')
					return *this;
			double z=1;
			while(isdigit(c))
				z/=10.,x=x+z*(c^48),c=getchar();
			return *this;
		}
		input operator>>(float &x) {
			x=0;
			bool y=1;
			char c=getchar();
			while(!isdigit(c))
				y&=(c!='-'),c=getchar();
			while(isdigit(c))
				x=x*10+(c^48),c=getchar();
			if(!y)
				x=-x;
			if(!isdigit(c))
				if(c!='.')
					return *this;
			double z=1;
			while(isdigit(c))
				z/=10.,x=x+z*(c^48),c=getchar();
			return *this;
		}
		input operator>>(std::string &x) {
			char c=getchar();
			x.clear();
			while(!(c!=' '&&c!='\n'&&c!='	'&&c!=EOF&&c))
				c=getchar();
			while(c!=' '&&c!='\n'&&c!='	'&&c!=EOF&&c) {
				x.push_back(c);
				c=getchar();
			}
			return *this;
		}
		input operator>>(char *x) {
			char c=getchar();
			int cnt=0;
			while(!(c!=' '&&c!='\n'&&c!='	'&&c!=EOF&&c))
				c=getchar();
			while(c!=' '&&c!='\n'&&c!='	'&&c!=EOF&&c) {
				x[cnt++]=c;
				c=getchar();
			}
			return *this;
		}
		input operator>>(char x) {
			x=getchar();
			return *this;
		}
	} pin;
};
inline void wt(char ch) {
	putchar(ch);
}
template<class T>
inline void wt(T x) {
	static char ch[40];
	int p=0;
	if(x<0)
		putchar('-'),x=-x;
	do
		ch[++p]=(x%10)^48,x/=10;
	while(x);
	while(p)
		putchar(ch[p--]);
}
template<class T,class... U>
inline void wt(T x,U ...t) {
	wt(x),
	wt(t...);
}
#define int long long
#define rep(i,a,b) for(int i=a,i##end=b;i<=i##end;i++)
#define frep(i,a,b) for(int i=a,i##end=b;i>=i##end;i--)
#define lqrep(i,a,v) for(int i=hd[a],v=e[i].v;i>=i##end;i=e[i].nxt,v=e[i].v)
const int N=1e5+7;
main() {
}
//目前快速输出pout还未搞定哦

到此就结束了吗？

当然不是！！！！！！！！！！

如果没有 $\gcd(P,A)=1$ 的话，前面的一切都成了一纸空文
那该如何？
如果不需要的旅客们可以摆烂了，以下是扩展板

二、BSGS基础算法

不妨设 $\gcd(P,A)=d$
$A^x\equiv B(\mod P)->\\(A'\times d)^x\equiv B'\times d(\mod P)->\\(A'\times d)^{x-1}\equiv B'*A'^{-1}(\mod P)$
接着按上文求解即可

Young_20220202

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
学习笔记——BSGS

众所周知，北上广深是中国非常一线的城市，北京是首都，地处……正片开始！一、BSGS基础算法实现目标：Ax≡B(mod P),(gcd⁡(P,A)=1)A^x\equiv B(\mod P),(\gcd(P,A)=1)Ax≡B(modP),(gcd(P,A)=1)求最小的xxx很明显，如果暴力枚举，时间是O(P)O(P)O(P)的，只要题目数据范围大，就死定了。愿意的人欢迎尝试（无100警告）于是，考虑 ~ 莫队     &n
复制链接

扫一扫