hdu1573中国剩余问题定理模数非互质的的情况

#include<iostream>
#define llint long long int
using namespace std;
const int MAXM=12;

llint exgcd(llint a, llint b, llint &x, llint &y)
{
if(b == 0){ x = 1, y = 0; return a;}
llint r = exgcd(b, a % b, x, y);
llint tmp = x;    x = y;  y = tmp - a/b * y;
return r;
}
long long int lcm(long long int a,long long int b)
{
llint x,y;
return (a*b)/(exgcd(a,b,x,y));
}
llint crr(llint mm[],llint rr[],llint n)
{
llint m1=mm[0];
llint r1=rr[0];
llint i;
bool flag;
for ( i = 0; i < n - 1; i++)
{
llint m2=mm[i+1];
llint r2=rr[i+1];
if (flag)
continue;
llint x,y;
llint d = exgcd(m1,m2,x,y);
llint c = r2 - r1;
if (c % d)
{
flag = 1;
continue;
}
llint t=m2/d;
x=(c/d*x%t+t)%t;
r1=m1*x+r1;
m1=m1*m2/d;
}
if(!flag)
{
if (r1==0&&n>1)
{
r1=mm[0];
__int64 ans=1;
llint x,y;
for (i=1;i<n;i++)
{
r1=exgcd(mm[i],r1,x,y);
}
for ( i=0;i<n;i++)
{
ans*=mm[i];
}
r1=ans/r1;
}
if (r1==0&&n==1)
r1=mm[0];
return r1;
}else
{
return -1;
}
}

int main()
{
llint T;
cin>>T;
while(T--)
{
llint N,M;
llint a[MAXM],b[MAXM];
cin>>N>>M;
llint i;
llint mt=1;
for(i=0;i<=M-1;i++)
{
cin>>a[i];
mt=lcm(mt,a[i]);
}
for(i=0;i<=M-1;i++)
{
cin>>b[i];
}

llint ans=crr(a,b,M);
//     cout<<ans<<"|"<<mt<<endl;
if(ans==-1)
{
cout<<0<<endl;
continue;
}
if(ans<N)
{
cout<<(N-ans)/mt+1<<endl;
}else
{
cout<<0<<endl;
}
//cout<<(N-ans)/mt<<endl;
}
return 0;
}