【校内模拟】游戏达人（BSGS求指标）（类欧几里得）

传送门

---

TM什么毒瘤省选模拟出类欧几里得。

题解：

由于本身对类欧几里得不是很熟，考场上也没有往类欧几里得方向去想。

结果你TM告诉我最后一步就是类欧几里得

首先特判无解和全零情况。

$C^x\equiv D^y(\mathrm{m}\mathrm{o}\mathrm{d}p)$看上去并不是很好操作，注意到$p$是质数，也就是说$p$有原根，那么我们直接求出$C,D$的指标$c,d$，则现在问题就是$cx\equiv dy(\mathrm{m}\mathrm{o}\mathrm{d}p-1)$

首先，如果要解是很好解的，但是现在问题不是解这个方程，而是最小化$Ax+By$。。。考场上没有什么时间去推，也就直接宣告凉凉。

下来滚去orz了题解的推导。

接下来是一系列非常牛逼的等价转换，设现在关于x,y的等价关系为$cx\equiv dy(\mathrm{m}\mathrm{o}\mathrm{d}p)$

首先同除以三者gcd，现在$gcd(c,d,p)=1$

设$g=gcd(p,c)$，显然$gcd(d,g)=1$，则$g$一定是$y$的因数，我们可以将$B$乘上$g$来进行等价变形。

得到的方程形式没有变，设上述变量表示现在的方程，现在设$q=gcd(p,d)$，同理$q$一定是$x$的因数，则将$A$乘上$q$来等价。

现在得到方程$cx\equiv dy(\mathrm{m}\mathrm{o}\mathrm{d}p)$，最小化$Ax+By$。

由于我们能够保证现在$gcd(c,p)=1$，即$c$在$\%p$下有逆元。方程可以直接转化为$x\equiv dy(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

考虑到$x,y$都是正整数，我们需要最小化$Ax+By$，且$A,B$均为正数，显然我们确定了$y$的时候$x=dy\%p=dy-p\lfloor \frac{dy}{p}\rfloor$是$x$的最优选择，并且由于$\gcd (d,y)=1$，在$y!=p$的时候$x$一定非0，对于$x=y=p$的情况可以直接扔掉，显然不够优秀。

把求值的式子转化一下得到$(Ad+B)y-Ap\lfloor \frac{dy}{p}\rfloor ,y\in [1,p-1]$

现在要求的就是 min ⁡ { a y + b ⌊ d y p ⌋ ∣ y ∈ [ l , r ] } \min\{ay+b\lfloor\frac{dy}{p}\rfloor\mid y\in[l,r]\} min{ay+b⌊pdy​⌋∣y∈[l,r]}

神仙类欧几里得。

设 f ( l , r , a , b , c , d , e ) = min ⁡ { a y + b ⌊ c y + d e ⌋ ∣ y ∈ [ l , r ] } f(l,r,a,b,c,d,e)=\min\{ay+b\lfloor\frac{cy+d}{e}\rfloor\mid y\in[l,r]\} f(l,r,a,b,c,d,e)=min{ay+b⌊ecy+d​⌋∣y∈[l,r]}，其中$c,d\ge 0,e>0$。

然后是xjb推导，设$p(y)=\lfloor \frac{cy+d}{e}\rfloor ,L=p(l),R=p(r)$

然后大力分类讨论

1. $L=R$根据$A$的正负性，答案在两个端点取到，直接算一下就行了。
2. $c\ge e$，类欧几里得常规操作，$f(l,r,a,b,c,d,e)=f(l,r,a+b\lfloor \frac{c}{e}\rfloor ,b,c\%e,d,e)$
3. $c<e,a\ge 0$，类欧几里得常规操作，此时显然有$L<R$，由于$c<e$，对于每个$k\in [L,R]$，存在集合$Y_k$，使得$\forall y\in Y_k,p(y)=k$，现在要最小化$ay+bk$，由于$a\ge 0$，所以现在显然考虑选择最小的$y$，即$k\le p(y)<k$的最小$y$。解不等式得到$$\frac{ke-d}{c}\le y<\frac{ke+e-d}{c}$$
   转化成形式上一样的下取整，即$$\lfloor \frac{ke-d-1}{c}\rfloor <y\le \lfloor \frac{ke+e-d-1}{c}\rfloor$$
   此时由于$c<e$，总能够保证$y$在该范围内有正整数取值，且最小的$y$为$\lfloor \frac{ke-d+c-1}{c}\rfloor$。发现只有$k=L$的情况$y$有可能取不到该最小值，直接计算即可。
   发现这个形式和上面的一模一样，问题转化为$f(L+1,R,b,a,e,-d+c-1,c)$
4. $c<e,a<0$，还是类欧几里得常规操作，$y$尽量取最大值，不等式和上面一样，特判右端点，问题转化为$f(L,R-1,b,a,e,e-d-1,c)$

然后就是找一下原根，然后BSGS求指标。

注意询问有点多，BSGS需要魔改块长。

然后我平生第一次知道了sqrt算负数会自动取unsigned值

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::cout;
using std::cerr;

int p,gr;

inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=(ll)a*a%p)(b&1)&&(res=(ll)res*a%p);
	return res;
}

struct Map{
	static cs int magic=4898597;
	int key[magic],val[magic],waste;
	Map(){memset(key,-1,sizeof key);}
	int locate(int k)cs{
		int h=k%magic;
		while(key[h]!=-1&&key[h]!=k)h=(h+1)%magic;
		return h;
	}
	int &operator[](int k){
		int h=locate(k);
		if(key[h]==-1)key[h]=k;
		return val[h];
	}
	cs int &operator[](int k)cs{
		int h=locate(k);
		return key[h]==-1?waste:val[h];
	}
	bool find(int k)cs{return key[locate(k)]==k;}
}ma;

inline void get_gr(){
	int phi=p-1;
	std::vector<int> pr;
	for(int re i=2;(ll)i*i<=phi;++i){
		if(phi%i==0){
			pr.push_back(i);
			while(phi%i==0)phi/=i;
		}
	}
	if(phi>1)pr.push_back(phi);
	for(gr=2;;++gr){
		bool flag=true;
		for(int p:pr)if(power(gr,(::p-1)/p)==1){flag=false;break;}
		if(flag)break;
	}
}

int step,ibase;
inline void init(int T){
	get_gr();//cout<<"gr : "<<gr<<"\n";
	step=ceil(sqrt((ll)T*p));
	int now=1;
	for(int re i=0;i<step;++i){
		ma[now]=i;
		now=(ll)now*gr%p;
	}
	ibase=power(now,p-2);
}

inline int get_ind(int v){
	for(int re i=0;;++i,v=(ll)v*ibase%p){
		if(ma.find(v))
		return i*step+ma[v];
	}
	assert(0);
	return -1;
}

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

inline ll f(ll l,ll r,ll a,ll b,ll c,ll d,ll e){
	if(l>r)return 1ll<<60;
	ll L=(c*l+d)/e,R=(c*r+d)/e;
	if(L==R)return std::min(a*l,a*r)+L*b;
	if(c>=e)return f(l,r,a+b*(c/e),b,c%e,d,e);
	if(a>=0)return std::min(a*l+b*L,f(L+1,R,b,a,e,c-d-1,c));
	return std::min(a*r+b*R,f(L,R-1,b,a,e,e-d-1,c));
}

inline void solve(){
	ll a,b,c,d;
	scanf("%lld%lld%lld%lld",&a,&b,&c,&d);
	c%=p,d%=p;
	if((c==0)^(d==0)){cout<<"Damn it\n";return ;}
	if(c==0){cout<<a+b<<"\n";return ;}
	c=get_ind(c),d=get_ind(d);int p=::p-1;
	int g=gcd(gcd(c,d),p);
	c/=g,d/=g,p/=g;
	g=gcd(c,p);c/=g,p/=g,b*=g;
	d=(ll)d*inv(c,p)%p;
	g=gcd(d,p);d/=g,p/=g,a*=g;
	cout<<std::min((a+b)*p,f(1,p-1,a*d+b,-a*p,d,0,p))<<"\n";
}

signed main(){
#ifdef zxyoi
	freopen("game.in","r",stdin);
#endif
	int T;scanf("%d%d",&p,&T);init(T);
	while(T--)solve();
	return 0;
}