中国剩余定理

最新推荐文章于 2021-07-29 23:07:25 发布

Edward-Tsui

最新推荐文章于 2021-07-29 23:07:25 发布

阅读量151

点赞数

分类专栏： C++ 模板算法

本文链接：https://blog.csdn.net/edward00324258/article/details/102906675

版权

C++ 同时被 3 个专栏收录

21 篇文章 0 订阅

订阅专栏

算法

18 篇文章 0 订阅

订阅专栏

模板

6 篇文章 0 订阅

订阅专栏

描述

用于解决 $x\mod a_1 = c_1, x \mod a_2 = c_2 ...$ 求 $x$ 的问题。

证明

从假设可知，对任何 $\in \{1, 2, \cdots , n\}$ ，由于$\forall j \in {1, 2, \cdots , n}, ; j\neq i, ; ; \operatorname{gcd}(m_i, m_j) = 1 $，所以$ \operatorname{gcd}(m_i, M_i) = 1.$ 这说明存在整数 $t_i$ 使得 $t_i M_i \equiv 1 \pmod {m_i}.$

对于乘积 $a_i t_i M_i$ 可知：
:$a_i t_i M_i \equiv a_i \cdot 1 = a_i \pmod {m_i}, $
😒\forall j \in {1, 2, \cdots , n}, ; j\neq i, ; ; a_j t_j M_j \equiv 0 \pmod {m_i}. $
所以 $a_1 t_1 M_1 + a_2 t_2 M_2 + \cdots + a_n t_n M_n$ 满足：
😒\forall i \in {1, 2, \cdots , n}, ; ; x = a_i t_i M_i + \sum_{j \neq i} a_j t_j M_j \equiv a_i + \sum_{j \neq i} 0 = a_i \pmod {m_i}. $
这说明 $x$ 就是方程组的一个解。

代码

int CRT(int a[],int m[],int n)
{
    int M = 1;
    int ans = 0;
    for(int i=1; i<=n; i++)
        M *= m[i];
    for(int i=1; i<=n; i++)
    {
        int x, y;
        int Mi = M / m[i];
        extend_Euclid(Mi, m[i], x, y);
        ans = (ans + Mi * x * a[i]) % M;
    }
    if(ans < 0) ans += M;
    return ans;
}