同余方程组——中国剩余定理の板子

最新推荐文章于 2021-08-27 12:22:44 发布

FromATP

最新推荐文章于 2021-08-27 12:22:44 发布

阅读量1.4k

点赞数 1

分类专栏：板子们烧脑的数论文章标签：中国剩余定理

本文链接：https://blog.csdn.net/FromATP/article/details/54231903

版权

烧脑的数论同时被 2 个专栏收录

16 篇文章 0 订阅

订阅专栏

板子们

15 篇文章 1 订阅

订阅专栏

中国剩余定理

中国剩余定理（CRT）是求解形如

x \equiv a 1 (m o d m 1)

$x \equiv a_1(mod\space m_1)$

x \equiv a 2 (m o d m 2)

$x \equiv a_2(mod\space m_2)$

x \equiv a 3 (m o d m 3)

$x \equiv a_3(mod\space m_3)$

. . . . . .

$......$

x \equiv a k (m o d m k)

$x \equiv a_k(mod\space m_k)$
的同余方程组的一种算法。一般的中国剩余定理要求

m1,m2,...,mk $m_1,m_2,...,m_k$ 两两互质，但它的扩展算法可以处理不互质的情况。

一般情况

当 $m_1,m_2,...,m_k$ 两两互质的时候，上面提到的同余方程组一定有整数解。设 $M=m_1\times m_2\times...\times m_k$ ，也就是M是所有模数的最小公倍数，那么上面的同余方程组在模M的意义下有唯一解。

那么中国剩余定理又是如何求解这个同余方程组的呢？首先我们要知道一个东西就是：如果一个整数 $N$ 满足 $N\space mod \space m=k$ ，那么对于任意整数v，有 $(v \times N)\space mod \space m \equiv (v \times k)\space mod \space m$ 。这个东西正确性比较显然啦，设 $N=r\times m+k$ 然后往里一代就可以证明了。

那么仍然回到上面的同余方程组，如果我们能够求解出这样的一些 $N_1,N_2,...,N_k$ ，满足对于每一个 $N_i$ ，它对除了 $m_i$ 以外的所有 $m$ 取模都是0，但除以 $m_i$ 余数为1。构造出这样的一堆东西有什么用呢？结合上面说的那个东西我们可以发现，令 $N=a_1N_1+a_2N_2+...+a_kN_k$ ，则N就是同余方程的一个合法整数解。因为每个 $N_i$ 累加进去不会破坏其它模数的余数而只会影响它自己的 $m_i$ 的余数，所以这样是正确的。

于是问题就变成了如何求解 $N_i$ 。仍然用 $M$ 表示所有模数的最小公倍数，将 $M \over m_i$ 表示为 $M_i$ ，那么 $N_i=M_i \times x$ 。又因为 $N_i\space mod \space m_i=1$ ，所以 $N_i=y \times m_i+1$ 。把这两种表示 $N_i$ 的方法一联立就有了 $M_i \times x-y \times m_i=1$ 。

这里就体现了我们保证模数互质的重要性。因为模数们两两互质，所以 $m_i$ 和 $M_i$ 也互质，那么就可以把上面含有 $x,y$ 的不定方程写作 $M_i*x+m_i*y=gcd(M_i,m_i)$ 。

这个式子就是经典的扩展欧几里得能够解决的式子了。因为我们要求的是 $x$ ，那么相当于求 $M_i$ 在模 $m_i$ 意义下的逆元。这个方程一定是有解的，那么我们可以对每一个 $M_i$ 和 $m_i$ 都求一下对应的 $x$ ，同时维护一个sum，每次累加上 $x \times M_i \times a_i$ ，最后就能求出一个解。

而因为我们可以注意到 $M$ 是所有模数的最小公倍数，它模所有 $m_i$ 都是0，所以求出的这个解加上 $M$ 的整数倍仍然是合法解。那么如果要求的是最小正整数解，我们只需要把一开始求出的解对 $M$ 取模就可以了。

long long CRT(long long M){
    long long sum=0,tmp,v;
    for (int i=1;i<=cnt;i++){
        tmp=M/m[i];
        v=getInv(tmp,m[i]);
        sum=(sum+tmp*a[i]*v)%M;
    }
    return sum;
}

扩展——模数不互质

如果同余方程组扩展到模数不互质的情况，这个时候首先不能用扩展欧几里得直接求逆元，并且还有可能会出现无解的情况。这个时候就不能用上面那种方法求解了，我们需要把方程两两合并化简。

设有两个方程 $x \equiv a_1(mod \space m_1)$ ， $x \equiv a_2 (mod \space m_2)$ 。把这两个方程展开写成 $x=a_1+n_1\times m_1$ ， $x=a_2+n_2\times m_2$ 。

这两个方程左边都是 $x$ ，把它们两个联立就变成了 $n_1 \times m_1+a_1=n_2 \times m_2 +a_2$ 。再移项化简一下就得到了 $n_1 \times m_1-n_2\times m_2=a_2-a_1$ 。

这又是一个形如 $ax+by=c$ 的不定方程。结合扩展欧几里得算法可以判断如果 $(a_2-a_1)\space mod \space gcd(m_1,m_2) \neq 0$ 的话这个同余方程组就是无解的。设 $gcd(m_1,m_2)=d$ ，我们可以用扩展欧几里得求出 $n_1$ 在模 $m_2 \over d$ 意义下的值 $N$ 。但这显然不一定是真正的 $n_1$ ，所以我们把 $n_1$ 写作 $y\times {m_2 \over d}+N$ 。

把 $n_1$ 代入第一个方程 $x=a_1+n_1\times m_1$ 就可以得到 $x=a_1+(y\times {m_2 \over d}+N)\times m_1$ ，把括号展开就可以得到 $x=a_1+N\times m_1+y\times {m_1 \times m_2 \over d}$ 。

这里y虽然是一个未知数，但我们可以把这个方程写作同余方程的形式，就是 $x \equiv a1+N\times m_1 (mod\space {m_1 \times m_2 \over d})$ 。那么这两个方程就合并成了一个方程，得到了新的余数和新的模数。那么我们就可以把给出的所有方程两两合并成新的方程，最后只剩下一个方程的时候当然就可以很方便地求解了。

bool merge(long long &a1,long long &m1,long long a2,long long m2){
    long long c,d,x,a3,m3;
    c=a2-a1;d=gcd(m1,m2);
    if (c%d!=0) return false;
    c=c/d;m1=m1/d;m2=m2/d;
    x=getinv(m1,m2);
    x=(x*c)%m2;
    x=x*(m1*d)+a1; 
    m3=m1*m2*d;
    a3=(x%m3+m3)%m3;
    a1=a3;m1=m3;
    return true;
}
long long CRT(){
    long long A=a[1],M=r[1];
    for (int i=2;i<=n;i++)
      if (!merge(A,M,a[i],r[i]))
        return -1;
    return (A%M+M)%M;
}