中国剩余定理证明

最新推荐文章于 2023-05-12 17:07:25 发布

转载最新推荐文章于 2023-05-12 17:07:25 发布 · 置顶 · 440 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：https://blog.csdn.net/niiick/article/details/80229217

文章标签：

#中国剩余定理 #孙子定理

知识点同时被 3 个专栏收录

13 篇文章

订阅专栏

定理

5 篇文章

订阅专栏

证明

4 篇文章

订阅专栏

本文深入探讨了中国剩余定理的原理与应用，详细解析了解同余方程组的方法，特别是针对两两互质的整数集合。文章通过数学证明阐述了定理的有效性，并提供了代码实现示例，包括使用扩展欧几里得算法求解同余式。

中国剩余定理

问题

解同余方程组

其中m1,m2,m3...mk 为两两互质的整数
求x的最小非负整数解

定理

令 $M= \prod \begin{matrix} k\\ i=1 \end{matrix}m_i$ 即M是所有 mi 的最小公倍数

ti为同余方程 $\frac{M}{mi}ti\equiv 1\left ( mod mi \right )$ 的最小非负整数解

则有一个解为 $x=\sum \begin{matrix} k\\ i=1 \end{matrix}ai\frac{M}{mi}ti$

通解为 $x+i*M\left ( i\in Z \right )$

特别的，最小非负整数解为（x mod M+M） mod M

证明

因为 $\frac{M}{mi}$ 是除mi之外的所有m的倍数

所以 $\forall k\neq i$ ， $a_i\frac{M}{m_i}t_i\equiv 0\left ( mod m_k \right )$

又有 $\frac{M}{m_i}t_i\equiv 1\left ( mod m_i \right )$

两边同时乘 $a_i$ 得 $a_i\frac{M}{m_i}t_i\equiv a_i\left ( mod m_i \right )$

带入 $x=\sum \begin{matrix} k \\ i=1\\ \end{matrix}a_i\frac{M}{m_i}t_i$

原方程组成立。

代码实现

$t_i$ 同余式的求解可以用扩展欧几里得

void exgcd(int a,int b,int &x,int &y)
{
    if(b==0){ x=1; y=0; return;}
    exgcd(b,a%b,x,y);
    int tp=x;
    x=y; y=tp-a/b*y;
}

int china()
{
    int ans=0,lcm=1,x,y;
    for(int i=1;i<=k;++i) lcm*=b[i];
    for(int i=1;i<=k;++i)
    {
        int tp=lcm/b[i];
        exgcd(tp,b[i],x,y);
        x=(x%b[i]+b[i])%b[i];//x要为最小非负整数解
        ans=(ans+tp*x*a[i])%lcm;
    }
    return (ans+lcm)%lcm;
}