【[NOI2018]屠龍勇士】
發現好像都是化掉係數之後套上\(ExCrt\)的板子
這好像是一個真正的擴充套件擴充套件中國剩餘定理
我們要處理的方程是這樣的形式
\[c_ix\equiv b_i(mod\ a_i)\]
其中\(c\)用一個\(std::multiset\)處理就好了
好像不是普通\(excrt\)可以處理的形式啊
思考一下這個方程的本質是什麼,\(c_ix=k_ia_i+b_i\)
所以如果我們有兩個方程
\[c_1x\equiv b_1(mod\ a_1)\]
\[c_2x\equiv b_2(mod\ a_2)\]
我們需要像\(crt\)那樣合併起來
我們開始化柿子了
\[x=\frac{k_1a_1+b_1}{c_1}=\frac{k_2a_2+b_2}{c_2}\]
所以就有
\[a_2c_1k_2+b_2c_1=a_1c_2k_1+b_1c_2\]
\[a_2c_1k_2=a_1c_2k_1+b_1c_2-b_2c_1\]
根據貝祖定理,這個方程有解條件是\(gcd(a_2c_1,a_1c_2)|(b_1c_2-b_2c_1)\)
如果有解的話,我們設\(t=gcd(a_2c_1,a_1c_2)\),兩邊除以\(t\)
\[\frac{a_2c_1k_2}{t}=\frac{a_1c_2k_1}{t}+\frac{b_1c_2-b_2c_1}{t}\]
顯然我們可以寫成一個同餘式
\[\frac{a_2c_1k_2}{t}\equiv \frac{b_1c_2-b_2c_1}{t}\ (mod\ \frac{a_1c_2}{t})\]
之後設\(inv=(\frac{a_2c_1}{t},\frac{a_1c_2}{t})\),即\(\frac{a_2c_1}{t}\)在\(mod\) \(\frac{a_1c_2}{t}\)意義下的乘法逆元
兩邊乘以\(inv\)
\[k_2\equiv inv*\frac{b_1c_2-b_2c_1}{t}\ (mod\ \frac{a_1c_2}{t})\]
改寫成等式
\[k_2= inv*\frac{b_1c_2-b_2c_1}{t}+y\times\frac{a_1c_2}{t}\]
我們把\(k_2\)迴帶到\(c_2x=k_2a_2+b_2\)
\[c_2x=inv*\frac{b_1c_2-b_2c_1}{t}*a_2+y\times\frac{a_1c_2a_2}{t}+b_2\]
再改寫成同餘式
\[c_2x\equiv inv*\frac{b_1c_2-b_2c_1}{t}\%(\frac{a_1c_2}{t})*a_2+b_2(mod\ \frac{a_1c_2a_2}{t})\]
我們只需要順次合併這些方程就好了,一旦出現無解就輸出\(-1\)好了
一個坑點是\(a=1\)時解出來會是\(0\),好像和實際要求不太一樣,所以對於這樣情況直接模擬特判就好了
之後因為非常的懶沒用快速乘,用了__int128
程式碼
#include<iostream>
#include<cstring>
#include<cstdio>
#include<set>
#define re register
#define maxn 100005
#define LL __int128
#define min(a,b) ((a)<(b)?(a):(b))
#define INF 9999999999999
#define max(a,b) ((a)>(b)?(a):(b))
#define set_it std::multiset<LL>::iterator
std::multiset<LL> s;
int n,m;
LL a[maxn],b[maxn],c[maxn],res[maxn];
LL gcd(LL a,LL b){ return !b?a:gcd(b,a%b);}
LL exgcd(LL a,LL b,LL &x,LL &y)
{
if(!b) return x=1,y=0,a;
LL r=exgcd(b,a%b,y,x);
y-=a/b*x;
return r;
}
inline LL pre(LL x)
{
s.insert(x);
set_it i=s.find(x);
if(i==s.begin()) return -INF;
--i;return *i;
}
inline void del(LL x)
{
set_it i=s.find(x);
s.erase(i);
}
inline int check(LL x)
{
if(s.find(x)!=s.end()) return 1;
return 0;
}
inline LL read()
{
char c=getchar();
LL x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9')
x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
inline LL inv(LL a,LL b)
{
LL x,y;
LL r=exgcd(a,b,x,y);
return (x%b+b)%b;
}
void write(LL x)
{
if(x>9) write(x/10);
putchar(x%10+48);
}
inline LL did(LL a,LL b)
{
if(a%b==0) return a/b;
return a/b+1;
}
inline void tepan()
{
LL ans=0;
for(re int i=1;i<=n;i++)
ans=max(ans,did(b[i],c[i]));
write(ans),puts("");
}
inline void solve()
{
s.clear();
LL x;
int flag=0;
n=read(),m=read();
for(re int i=1;i<=n;i++) b[i]=read();
for(re int i=1;i<=n;i++) a[i]=read(),flag|=(a[i]!=1);
for(re int i=1;i<=n;i++) res[i]=read();
for(re int i=1;i<=m;i++) x=read(),s.insert(x);
for(re int i=1;i<=n;i++)
{
if(check(b[i])) c[i]=b[i],del(b[i]);
else
{
LL now=pre(b[i]);
del(b[i]);
if(now==-INF)
{
set_it it=s.begin();
c[i]=*it;
s.erase(it);
}
else c[i]=now,del(now);
}
s.insert(res[i]);
}
if(!flag)
{
tepan();
return;
}
LL a1=a[1],b1=b[1],c1=c[1];
for(re int i=2;i<=n;i++)
{
LL a2=a[i],b2=b[i],c2=c[i];
LL r=gcd(c2*a1,c1*a2);
if((b1*c2-b2*c1)%r)
{
puts("-1");
return;
}
if((c2*a1/r)<(c1*a2/r)) std::swap(a1,a2),std::swap(b1,b2),std::swap(c1,c2);
LL P=a1*c2/r,Inv=inv(c1*a2/r,a1*c2/r);
b1=((Inv*(((b1*c2-b2*c1)/r)%P+P)%P*a2))+b2;
a1=a1*a2*c2/r;
c1=c2;
}
LL y;
if(b1%gcd(a1,c1))
{
puts("-1");
return;
}
LL r=exgcd(c1,a1,x,y);
LL t=a1/r;
x=(x*(b1/r)%t+t)%t;
write(x);
puts("");
}
int main()
{
int T;
T=read();
while(T--) solve();
return 0;
}