中国剩余定理

转自这里

 代码：

void exgcd(int a, int b, int &d, int &x, int &y)
{
    if(!b) { d = a; x = 1; y = 0; }//d为最大公约数 
    else { exgcd(b, a%b, d, y, x); y -= x*(a/b); }
}

int CRT(int a[],int m[],int n)  
{  
    int M = 1;  
    int ans = 0;  
    for(int i=1; i<=n; i++)  
        M *= m[i];  
    for(int i=1; i<=n; i++)  
    {  
        int x, y, d;  
        int Mi = M / m[i];  
        exgcd(Mi, m[i], d, x, y);
        ans = (ans + Mi * x * a[i]) % M;  
    }  
    if(ans < 0) ans += M;  
    return ans;  
}