NOI 2018 屠龙勇士（exgcd，excrt）

题目
对于方程$ax\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}p)$
记$\gcd (a,p)=g$
则求出$ax\equiv g(\mathrm{m}\mathrm{o}\mathrm{d}p)$的一组特解$x_0$后。
若$b\mathrm{m}\mathrm{o}\mathrm{d}g\ne 0$
则无解。
否则原方程等价于$x\equiv \frac{x_{0}b}{g}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{p}{g})$

对于同余方程组$x\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}c),x\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}d)$
可以设$x=k_{1}c+a=k_{2}d+b$
则$k_{1}c-k_{2}d=b-a$
用$\mathrm{e}\mathrm{x}\mathrm{g}\mathrm{c}\mathrm{d}$求出$k_1$的一个特解$k_0$（记住中间要判断无解）
则可以合并这两个同余方程为$x\equiv k_{0}c+a(\mathrm{m}\mathrm{o}\mathrm{d}\frac{cd}{\gcd (c,d)})$

然后就可以解决此题了。
$\mathrm{A}CCode$

#include<bits/stdc++.h>
#define LL long long
#define maxn 100005
using namespace std;

int n,m;
LL a[maxn],p[maxn],b[maxn],c[maxn];
LL mul(LL a,LL b,LL p){
	LL r = a * b - (LL)((long double)a / p * b) * p;
	if(r < 0) return r + p;
	if(r >= p) return r - p;
	return r;
}
multiset<LL>st;
void exgcd(LL a,LL b,LL &x,LL &y,LL &gcd){
	if(!b) gcd=a,x=1,y=0;
	else exgcd(b,a%b,y,x,gcd),
		y -= (a/b) * x;
}

int main(){
	int T;
	scanf("%d",&T);
	for(;T--;){
		scanf("%d%d",&n,&m);
		for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
		for(int i=1;i<=n;i++) scanf("%lld",&p[i]);
		for(int i=1;i<=n;i++) scanf("%lld",&b[i]);
		st.clear();
		for(int i=1;i<=m;i++) scanf("%lld",&c[i]),st.insert(c[i]);
		LL lim = 1 , x = 0 , mod = 1;
		for(int i=1;i<=n;i++){
			set<LL>::iterator it = st.upper_bound(a[i]);
			if(it != st.begin()) it--;
			LL v = *it , sx , sy , gx;
			st.erase(it);
			exgcd(v,p[i],sx,sy,gx);
			if(a[i] % gx){
				mod = -1;
				break;
			}
			sx = mul(sx,a[i]/gx,p[i]/gx);
			gx = p[i] / gx;
			lim = max(lim , (a[i] + v - 1) / v);
			
			LL k1,k2,gd;
			exgcd(gx,mod,k1,k2,gd);
			if((x-sx) % gd){
				mod= -1;
				break;
			}
			mod = mod * (gx / gd);
			k1 = mul(k1 , (x-sx) / gd , mod);
			x = (mul(k1 , gx , mod) + sx) % mod;
			st.insert(b[i]);
		}
		if(mod == -1){
			puts("-1");
			continue;
		}
		LL K = (lim - x + mod - 1) / mod;
		x = max(x , K * mod + x);
		printf("%lld\n",x);
	}
}