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

#include<bits/stdc++.h>
using namespace std;
#define rd(x) scanf("%d",&x)
#define rdd(x,y) scanf("%d%d",&x,&y)
#define rddd(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define rds(s) scanf("%s",s)
#define rep(i,n) for(int i=0;i<n;i++)
#define LL long long
const int N = 2e5+10;
const int inf=1e9;
const int MOD=1e9+7;
int n,m,k;
LL gcd(LL a,LL b){
    return b==0?a:gcd(b,a%b);
}
void add(LL&a,LL b,LL md){
    a+=b;
    if(a>=md) a-=md;
}
LL mul(LL a,LL b,LL md){
    LL ret=0;
    while(b){
        if(b&1) add(ret,a,md);
        add(a,a,md);
        b>>=1;
    }
    return ret;
}
void Exgcd(LL a,LL b,LL&x,LL& y){
    if(b==0){
        x=1;y=0;
        return;
    }
    LL x2,y2;
    Exgcd(b,a%b,x2,y2);
    x=y2;
    y=x2-a/b*y2;
}
LL rev(LL a,LL b){
    LL x,y;
    Exgcd(a,b,x,y);
    x%=b;
    x+=b;
    x%=b;
    return x;
}
int a[15];
int b[15];
bool merge(LL a1,LL n1,LL a2,LL n2,LL&a3,LL& n3){
     LL d=gcd(n1,n2);
     LL a=a2-a1;
     if(a%d) return false;
     a/=d;
     a%=(n2/d);a+=(n2/d);a%=(n2/d);
     LL k=mul(a,rev(n1/d,n2/d),n2/d);
     n3=n1*n2/d;
     add(a1,mul(n1,k,n3),n3);
    a1%=n3;
    a1+=n3;
    a1%=n3;
    a3=a1;
    return true;
}
LL china_reminder(int len,int *a,int *b){
    LL a1,n1,a2,n2;
   a1=a[0];n1=b[0];
    for(int i=1;i<m;i++){
            a2=a[i];n2=b[i];
            LL a3,n3;
            if(!merge(a1,n1,a2,n2,a3,n3)) return -1;
            a1=a3;n1=n3;
    }
    a1%=n1;a1+=n1;a1%=n1;
    return a1;
}
int main()
{
#ifndef ONLINE_JUDGE
  freopen("aaa","r",stdin);
#endif
    int T,q;
    rd(T);
    while(T--){
        rdd(n,m);
        rep(i,m) rd(b[i]);
        rep(i,m) rd(a[i]);
        LL M=1;
        rep(i,m) M=M*b[i]/gcd(M,b[i]);
        LL ans=china_reminder(m,a,b);
        if(ans==-1 || ans>n) puts("0");
        else{
            LL ret=1+(n-ans)/M;
            if(ans==0) ret--;
            printf("%I64d\n",ret);
        }
    }
    return 0;
}