POJ2891Strange Way to Express Integers

POJ2891

mi不互質時解同餘方程組
存一發模板吧QAQ
講解看看：這裏
那個爲什麼要mod (m2/d)還是想了很久。。
大概就是最後要在mod（m2*m1/d）意義下，然後帶入到x=m1∗t1+a1 中乘了m1 ，所以前面mod（m2∗m1/d） 就可以滿足條件。

【代碼】

#include <cstdio>
#include <iostream>
#include <algorithm>
#define N 2335
#define INF 0x7fffffff
using namespace std;
typedef long long ll;

ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n;
ll M,R;
ll m[N],r[N]; 

void Exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(b) {
        Exgcd(b,a%b,d,y,x);
        y-=(a/b)*x;
    }
    else x=1,y=0,d=a;
}

ll Ex_Crt()
{
    M=m[1],R=r[1];ll x,y,d;
    for(int i=2;i<=n;i++)
    {
        Exgcd(M,m[i],d,x,y);
        if((r[i]-R)%d) return -1;
        R=(r[i]-R)/d*x%(m[i]/d)*M+R;
        M*=m[i]/d;
        R%=M;
    }
    return (R+M)%M;
} 

int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=1;i<=n;i++) m[i]=read(),r[i]=read();
        printf("%lld\n",Ex_Crt());
    } 
    ll x,y,d;
    Exgcd(24,6,d,x,y);
    printf("%lld %lld\n",x,y);
    return 0;
}