这题好像只要会用set/平衡树以及裸的(Excrt)就能A啊...然而当时我虽然看出是(Excrt)却并不会...今天又学了一遍(Excrt)，趁机把这个坑给填了吧
现预处理一下，找出每条龙用哪吧剑，把所有龙都砍(tmp)刀到负血。
设之后每条龙都砍了a刀，对于第(i)条龙，剑的攻击力为(w_i)，恢复能力为(c_i)，血量为(b_i)
则根据题意，满足
[b_i-aw_i+yc_i=0
]
[aw_iequiv b_i(mod c_i)
]
将(w_i,b_i,c_i)同时除以(gcd(w_i,c_i))然后
[aequiv b_iinv(w_i)(mod c_i)
]
这道题就变成了同余方程组求解：
[egin{cases}
aequiv b_1inv(w_1)(mod c_1)\
aequiv b_2inv(w_2)(mod c_2)\
...\
aequiv 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
");
else printf("%lld
",(bi%ci+ci)%ci+tp);
}
return 0;
}