#include <algorithm>
#include <cstdio>
#include <iostream>
using namespace std;
typedef long long ll;
const int maxn=15;
ll ex_gcd(ll a,ll b,ll &x,ll &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
ll g=ex_gcd(b,a%b,x,y),t;
t=x;
x=y;
y=t-(a/b)*y;
return g;
}
ll mod(ll x,ll mm)
{
return (x%mm+mm)%mm;
}
pair<ll,ll> Both(ll B[], ll C[], int n)
{//求解x=B[i](mod C[i]),总共n个线性方程组
ll x,y,g,b=0,c=1;
for(int i=0;i<n;i++)
{
ll aa=c,bb=B[i]-b,cc=C[i];
g=ex_gcd(aa,cc,x,y);
if(bb%g) return make_pair(-1,0);
x=bb/g*mod(x,cc);
b+=c*x;
c*=cc/g;
b=mod(b,c);
}
return make_pair(b,c);
}
int main()
{
ios::sync_with_stdio(false);
int t;
cin>>t;
while(t--)
{
ll m;
ll c[maxn],mo[maxn];
int n;
cin>>m>>n;
for(int i=0;i<n;i++)
cin>>mo[i];
for(int i=0;i<n;i++)
cin>>c[i];
pair<ll,ll> ans=Both(c,mo,n);
if(ans.first==-1||m<ans.first)
cout<<0<<endl;
else
cout<<(m-ans.first)/ans.second+(ans.first?1:0)<<endl;
}
}
hdu1573 同余方程
最新推荐文章于 2018-12-20 21:12:30 发布