中国剩余定理和一元同余线性方程的详细讲解

关于中国剩余定理和一元同余线性方程的知识点以前就看过，但是不太怎么理解，今天参考了些资料，终于弄透了，下面来讲解一~

(1)中国剩余定理：若m1，m2.......mr两两互质，则下列同余方程组

                                                                                           x===a1(mod m1)

                                                                                           x===a2(mod  m2)

                                                                                           .

                                                                                           .

                                                                                           x===ar(mod mr)

有模M=m1*m2*m3*.......mr的唯一解；  

令Mi=M/mi;  那么肯定gcd(Mi,mi)==1，所以 Mi*Pi===1(mod  mi )===》  ai*Mi*Pi===ai(mod mi);

所以a1*M1*P1+a2*M2*P2+...........................ar*Mr*P就是同余方程的解了；

下面附上代码：

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

（2）上面所说的是若干个互素的整数，假如若干个整数不互素该怎么办呢？我们可以用合并方程组的思想，

    参考资料：http://yzmduncan.iteye.com/blog/1323599/

 上面公式的推导过程已经很详细了，自己动手推一推；

附上代码：

int sove()
{
  cin>>a1>>r1;
  for(int i=1;i<n;i++)
  {
    cin>>a2>>r2;
    a=a1;
    b=a2;
    c=r2-r1;
    exgcd(a,b,d,x0,y0);
    if(c%d!=0)
    {
      ishave=0;
    }
    int t=b/d;
    x0=(x0*(c/d)%t+t)%t;
    r1=a1*x0+r1;
    a1=a1*(a2/d); 
  }
  if(!ishave)
  {
    r1=-1;  
  }
  return r1;
}

如果有疑问的地方，请指教~~~