数论--中国剩余定理

最新推荐文章于 2023-11-03 10:37:00 发布

exmy

最新推荐文章于 2023-11-03 10:37:00 发布

阅读量608

点赞数

分类专栏：数学文章标签： math 数论

本文链接：https://blog.csdn.net/xushao_Movens/article/details/51151354

版权

数学专栏收录该内容

19 篇文章 0 订阅

订阅专栏

中国剩余定理（Chinese Remainder Theorem: CRT）：求解一次同余式组。
给出以下一元线性同余方程组：

(S) : ⎧ ⎩ ⎨ x \equiv a 1 (m o d m 1) x \equiv a 2 (m o d m 2) x \equiv a 3 (m o d m 3)

$(S): \begin{cases} x \equiv a_1 (mod \text{ } m_1) \\ x \equiv a_2 (mod \text{ } m_2) \\ x \equiv a_3 (mod \text{ } m_3) \end{cases}$
定理描述：方程组有解当且仅当

m1,m2,...,mn $m_1, m_2, ..., m_n$ 两两互质，通解按如下方式构造：
设

M=∏i=1nmi $M = \prod\limits_{i = 1}^n m_i$ ,

Mi=M/mi $M_i = M/m_i$
设

ti=M−1i $t_i = M_i^{-1}$ 为

Mi $M_i$ 模

mi $m_i$ 的数论倒数（逆元）:

tiMi≡1(mod mi) $t_iM_i \equiv 1 (mod \text{ } m_i)$
则解为：

x = \sum i = 1 n a i t i M i

$x = \sum\limits_{i = 1}^n a_it_iM_i$

代码实现：

#include <iostream>
#include <algorithm>
using namespace std;

// ax + by = gcd(a, b) ==> d = gcd(a, b)
int exgcd(int a, int b, int& x, int& y)
{
    if(b == 0){ x = 1; y = 0; return a; }
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

// ax = 1 (mod m) ==> x
int mod_reverse(int a, int m)
{
    int x, y;
    int d = exgcd(a, m, x, y);
    if(d == 1) return (x % m + m) % m;
    else return -1;
}

int CRTsolve(int a[], int m[], int n)
{
    int M = 1;
    int ans = 0;
    for(int i = 0; i < n; ++i) M *= m[i];
    for(int i = 0; i < n; ++i){
        int Mi = M / m[i], ti;
        ti = mod_reverse(Mi, m[i]);
        ans += ti * Mi * a[i];
    }
    return ans % M;
}

int main()
{
    int m[] = {3, 5, 7};
    int a[] = {2, 3, 2};

    cout << CRTsolve(a, m, sizeof(a) / sizeof(a[0]));

    return 0;
}