数学：拓展中国剩余定理

拓展中国剩余定理解决模数不互质同余方程组

如果模数互质的话，直接中国剩余定理就可以了

但是如果模数不互质又没有接触这个方法就凉凉了

推是很不好推出来的

假设我们这里有两个方程： 

x=a1∗x1+b1

x=a2∗x2+b2 

a1,a2是模数，b1,b2是余数

那么我们可以合并这两个方程：

a1∗x1+b1=a2∗x2+b2

由于x1和x2可以取负无穷到正无穷，所以符号不能约束它们，我们随便变一变形得到

a1∗x1+a2∗x2=b2−b1

然后使用拓展欧几里德算法，x和y分别是式子中的x1和x2

我们求出了一个最小正整数解x1

令k=(a1∗x1+b1)

x≡k(mod lcm(a1,a2))

一路合并下去就可以得到最终的解答了

典型例题是POJ2891,全网仅此一道？？

POJ真是交上去立刻A

#include<cstdio>
using namespace std;
const int maxn=100005;
int n;
long long a[maxn],m[maxn];
long long exgcd(long long a,long long b,long long &x,long long &y)
{
    if(b==0) {x=1;y=0;return a;}
    long long ret=exgcd(b,a%b,x,y);
    long long t=x;x=y;y=t-a/b*y;
    return ret;
}
long long crt()
{
    long long M=m[1],A=a[1];
    long long x,y;
    for(int i=2;i<=n;i++)
    {
        long long d=exgcd(M,m[i],x,y);
        if((a[i]-A)%d) return -1;  //无解
        //计算x的值 
        x*=(a[i]-A)/d;
        long long t=m[i]/d;
        x=(x%t+t)%t;
        A=M*x+A;
        M=M/d*m[i];
        A%=M; 
    }
    A=(A%M+M)%M;
    return A;
}
int main()
{
    while(scanf("%d",&n)==1)
    {
        for(int i=1;i<=n;i++)
            scanf("%lld%lld",&m[i],&a[i]);
        printf("%lld
",crt());
    }
    return 0;
}