正题

题目链接:https://www.luogu.com.cn/problem/P4774

---

题目大意

$n$个龙血量为$a_i$，回复能力为$p_i$，死亡后掉落剑的攻击力$t_i$，$m$把剑，攻击力为$b_i$。

从$1$开始打，每次使用不大于当前龙血量的剑中攻击力最低的一把（没有就用攻击力最低的），造成$x\times atk$点伤害，然后当前的剑坏掉。

求一个最小的$x$使得所有龙被攻击后血量是$p_i$的倍数。

$1\le n,m\le 10^5$，满足$p=1$或者$a_i\le p_i$，所有$p_i$的公倍数不超过$10^{12}$

---

解题思路

额，先用$set$处理出每个龙用哪把剑打$c_i$，然后就是对于每条龙的条件就是
$$c_{i}x\equiv a_i(mod{p}_i),c_{i}x\ge a_i$$
后面那个条件可以去掉，我们先求出满足所有$c_{i}x\ge a_i$的最小$x$，后面再调整。

然后前面那个东西可以用$EXCRT$搞了，假设我们上一次求到的答案为$ans$，目前$p_i$的公倍数是$M$，那么现在的通解就是$ans+Mx$。我们需要求出一个$x$满足
$$c_i(ans+Mx)\equiv a_i(mod{p}_i)$$
$$\Rightarrow Mx+p_{i}y=a_i-c_i\times ans$$
然后就可以扩欧合并了。

---

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#include<set>
#define ll __int128
using namespace std;
ll read(){ll x=0,f=1; char c=getchar();while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();return x*f;
}
void print(ll x)
{if(x>9)print(x/10);putchar(x%10+48);return;}
const ll N=1e5+10;
multiset<ll> s;
ll n,m,T,a[N],p[N],t[N],c[N];
ll exgcd(ll a,ll b,ll &x,ll &y){if(!b){x=1;y=0;return a;}ll d=exgcd(b,a%b,x,y);ll z=x;x=y;y=z-a/b*y;return d;
}
void work(){n=read();m=read();s.clear();ll mx=0;for(ll i=1;i<=n;i++)a[i]=read();for(ll i=1;i<=n;i++)p[i]=read();for(ll i=1;i<=n;i++)t[i]=read();for(ll i=1;i<=m;i++){ll x=read();s.insert(x);}for(ll i=1;i<=n;i++){multiset<ll>::iterator it;if(a[i]<*s.begin())it=s.begin();else it=--s.upper_bound(a[i]);c[i]=*it;s.erase(it);s.insert(t[i]);mx=max(mx,(a[i]-1)/c[i]+1);}ll M=1,x,y,ans=0;for(ll i=1;i<=n;i++){ll d=exgcd(M*c[i],p[i],x,y);ll w=a[i]-c[i]*ans,v=p[i]/d;if(w%d){puts("-1");return;}x=w/d*x%v;ans=ans+x*M;M=M*v;ans=(ans%M+M)%M;}if(ans<mx)ans+=(mx-ans+M-1)/M*M;print(ans);putchar('\n');
}
signed main()
{scanf("%lld",&T);while(T--){work();}return 0;
}