中国剩余定理

 

 

中国剩余定理不互质情况下

POJ-2891

参考博客：http://yzmduncan.iteye.com/blog/1323599/

               https://www.cnblogs.com/linyujun/p/5199415.html

中国剩余定理不互质的情况下就是相当于把某个表达式通过另外的表达式表达出来

 假设我们现有2个表达式

X = a1 (mod m1)

X = a2 (mod m2)

等于 X = a1 + k1 * m1  = a2 + k2 * m2;

然后可以通过移项可以得到  k1 * m1 - k2 * m2 = a2 - a1  （1）

通过拓展欧几里得 我们可以知道 若此方程有整数解 必满足 (a2-a1)%gcd(m1,m2) == 0

（1）式可以变换成 k1*m1 = a2 - a1 (mod m2)

设 d = gcd(m1, m2), c = a2-a1

(m1/d) * k1 = c/d (mod m2/d)

k1 = c/d * inv(m1/d)  (mod m2/d)

设 K = c/d * inv(m1/d) 

k1 = k + y * (m2/d)

再将k1 带入 X = a1 + k1 * m1 可以得到

X =  a1 + k*m1 + y*(m1*m2/d)

即 X = a3 (mod m3)    (其中 a3 = a1 + k*m1  m3 = m1*m2/d)

然后通过这个式子对剩下的表达式进行整合 就可以求得最后的答案。

#include<cstdio>  
#define ll long long  
using namespace std;  
const int N = 100000+5;  
ll M[N], A[N];// M->mod  A -> 余数  
ll gcd(ll a, ll b) {  
    return b == 0 ? a : gcd(b, a%b);  
}  
void ex_gcd(ll a, ll b, ll &x, ll &y, ll &d)  
{  
    if(b == 0) {  
        x = 1, y = 0, d = a;  
        return ;  
    }  
    ex_gcd(b, a%b, y, x, d);  
    y -= a/b*x;  
    return ;  
}  
ll inv(ll n, ll p) {  
    ll x, y, d;  
    ex_gcd(n, p, x, y, d);  
    return (d == 1) ? (x % p + p) % p : -1;  
}  
ll china_(ll A[], ll M[], ll n)  
{  
    ll x = A[0], m = M[0];  
    for(int i = 1; i < n; i++) {  
        ll c = A[i] - x, d = gcd(M[i], m);  
        if(c % d) return -1;  
        ll k = c/d * inv(m/d, M[i]/d)%(M[i]/d);  
        x += m*k;  
        m *= M[i]/d;  
    }  
    return (x%m+m)%m;  
}  
int main()  
{  
    int n;  
    while(~scanf("%d", &n)) {  
        for(int i = 0; i < n; i++)  
            scanf("%lld%lld", &M[i], &A[i]);  
        ll ans = china_(A, M, n);  
        printf("%lld", ans);  
    }  
    return 0;  
}