这题好像只要会用set/平衡树以及裸的\(Excrt\)就能A啊...然而当时我虽然看出是\(Excrt\)却并不会...今天又学了一遍\(Excrt\)，趁机把这个坑给填了吧
现预处理一下，找出每条龙用哪吧剑，把所有龙都砍\(tmp\)刀到负血。
设之后每条龙都砍了a刀，对于第\(i\)条龙，剑的攻击力为\(w_i\)，恢复能力为\(c_i\)，血量为\(b_i\)
则根据题意，满足
\[ b_i-aw_i+yc_i=0 \]
\[ aw_i\equiv b_i(mod\ c_i) \]
将\(w_i,b_i,c_i\)同时除以\(gcd(w_i,c_i)\)然后
\[ a\equiv b_iinv(w_i)(mod\ c_i) \]
这道题就变成了同余方程组求解：
\[ \begin{cases} a\equiv b_1inv(w_1)(mod\ c_1)\\ a\equiv b_2inv(w_2)(mod\ c_2)\\ ...\\ a\equiv b_ninv(w_n)(mod\ c_n)\\ \end{cases} \]
解出\(a\)，答案就是\(a+tmp\)
代码
#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef multiset<ll> MT;
const int maxn=1e5+100;
int n,m;
ll b[maxn],c[maxn],w[maxn],a,bi,ci,k[maxn],sxz[maxn];
bool err;
MT st;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
void exgcd(ll a,ll b,ll &x,ll &y){b?(exgcd(b,a%b,y,x),y-=a/b*x):(x=1,y=0);}
ll inv(ll a,ll b){ll x,y; exgcd(a,b,x,y); return (x%b+b)%b;}
ll e(ll a,ll b,ll p){
ll x=a*b-(ll)((long double)a*b/p+0.5)*p;
return x<0?x+p:x;
}
void excrt(ll bj,ll cj){
if(bi==0&&ci==0){
bi=bj,ci=cj;
return;
}
ll C=bj-bi,d=gcd(ci,cj),P=ci/d*cj;
if(C%d!=0){err=1; return;}
ll K=e(C/d%(cj/d),inv(ci/d,cj/d),(cj/d));
bi=(e(K,ci,P)+bi)%P,ci=P;
}
int main(){
int t;
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%lld",&b[i]);
for(int i=1;i<=n;i++) scanf("%lld",&c[i]);
for(int i=1;i<=n;i++) scanf("%lld",&w[i]);
for(int i=1;i<=n;i++)
st.clear(),bi=ci=0,err=0;
for(int i=1;i<=m;i++)
scanf("%lld",&a),st.insert(a),sxz[i]=a;
sort(sxz+1,sxz+m+1);
ll tp=0;
for(int i=1;i<=n;i++){
MT::iterator p=st.lower_bound(b[i]);
if(p==st.end()||(p!=st.begin()&&*p>b[i]))
p--;
k[i]=*p,tp=max(tp,(b[i]-1)/k[i]+1);
st.erase(p);
st.insert(w[i]);
}
for(int i=1;i<=n;i++){
b[i]-=tp*k[i],b[i]=(b[i]%c[i]+c[i])%c[i];
int lp=gcd(k[i],c[i]);
if(b[i]%lp!=0){
err=1;
break;
}
c[i]/=lp,b[i]/=lp,k[i]/=lp;
b[i]=e(b[i],inv(k[i],c[i]),c[i]);
}
for(int i=1;i<=n;i++){
excrt(b[i],c[i]);
if(err)
break;
}
if(err) printf("-1\n");
else printf("%lld\n",(bi%ci+ci)%ci+tp);
}
return 0;
}