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

依旧没弄懂。。先会用模版再说。。

注意讨论(ans==-1)的情况

#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;

}