中国剩余定理

与我梦中

于 2022-10-23 12:19:25 发布

阅读量336

点赞数

分类专栏：数论文章标签：算法

本文链接：https://blog.csdn.net/qq_64637770/article/details/127459448

版权

数论专栏收录该内容

6 篇文章 0 订阅

订阅专栏

一.算法描述

情况1：mi两两互质

中国剩余定理是指解决求一次同余式组问题的方法，求一个最小的非负整数 $x$ ，满足

$x\equiv ai(mod \, \, \, mi) \: \: \: \: \: i \in [1.n]$

其中 $m_{i}$ 两两互质。

我们先回忆逆元的定义：对于任意 $a$ 满足：

$\frac{a}{b}\equiv ax(mod\: \: \: m)$

这样，我们称 $x$ 为 $b$ 的模 $m$ 的乘法逆元，记作 $b^{-1}$ .

这样，我们可以得到：

$\frac{a}{b}\equiv a*b^{-1}(mod\, \, \, m)$

$a\equiv a*b*b^{-1}(mod\, \, \, m)$

$b*b^{-1}\equiv 1(mod\, \, \, \, m)\cdot \cdot \cdot \cdot \cdot \cdot (*)$

这样，我们就可以给出 $x$ 的一个通式:

$x= \sum_{i=1}^{n}a_{i}M_{i}M_{i}^{-1}$

其中：

$M_{i}=\frac{M}{m_{i}},M=\prod_{i=1}^{n}=m_{1}m_{2}\cdot \cdot \cdot m_{n}$

下给出证明：

我们只需要证明，对于任意的 $i$ ，会满足 $x\equiv ai(mod \, \, \, mi)$ 即可。

而对于 $x$ 中的任意一项 $j$ :

1.若 $i=j$ ,那么 $a_{j}M_{j}M_{j}^{-1}=a_{i}M_{i}M_{i}^{-1}$ ,而 $M_{i}M_{i}^{-1}\equiv 1(mod\, \, \, m_{i})$ ,所以：

$a_{j}M_{j}M_{j}^{-1}\equiv a_{i}\left ( mod\, \, \, m_{i} \right )(i=j)$

2.若 $i\neq j$ ，那么对于 $a_{j}M_{j}M_{j}^{-1}$ 中，由于 $M_{j}=\frac{M}{m_{j}}$ ，这样 $M$ 一定是包括了 $m_{i}$ 的，所以：

$a_{j}M_{j}M_{j}^{-1}\equiv 0\left ( mod\, \, \, m_{i} \right )(i\neq j)$

综上，对于任意的 $i$ ， $x\equiv ai(mod \, \, \, mi)$ 。

下来，我们需要去求出 $M_{i}$ 关于 $m_{i}$ 的逆元 $M_{i}^{-1}$ ：

事实上，在前面扩展欧几里得算法中，我们已经可以求解：

$ax\equiv b(mod\, \, \, m)$

那我们现在要求 $M_{i}*M_{i}^{-1}\equiv 1(mod\, \, \, \, m_{i})$ 而已：

$M_{i}M_{i}^{-1}=km_{i}+1\rightarrow M_{i}M_{i}^{-1}+km_{i}=1(mod\, \, \, m_{i})$

但是我们要注意 $M_{i}$ 和 $m_{i}$ 是互质的，因而 $gcb(M_{i},m_{i})=1$ .

这样我们就可以用扩展欧几里得算法求出 $M_{i}^{-1}$ 。

情况二：mi不能保证两两互质

在mi不能保证两两互质的情况下，那么中国剩余定理就不能使用了，因为你不能保证每一个 $M_{i}$ 都存在逆元。

这时候我们就采用另一种方法：

首先我们考虑前面的两个式子：

$x\, \, \, mod\, \, \,m _{1}=a_{1}$

$x\, \, \, mod\, \, \,m _{2}=a_{2}$

那我们可以把这两个式子写成：

$x=k_{1}m_{1}+a_{1}$

$x=k_{2}m_{2}+a_{2}$

这样：

$k_{1}m_{1}+a_{1}=k_{2}m_{2}+a_{2}$

$k_{1}m_{1}-k_{2}m_{2}=a_{2}-a_{1}$

同样的，我们令 $d=gcb(a_{1},a_{2})$ ,由于 $a_{1},a_{2}$ 都是 $d$ 的倍数，那么上式有解等价于 $m_{2}-m_{1}$ 可以整除于 $d$ 。

这样，我们可以写出 $x$ 的一个通解：

$x=(k_{1}+k\frac{m_{2}}{d})m_{1}+a_{1}\cdot \cdot \cdot \cdot \cdot \cdot (*)$

为什么捏？这是因为

$x\, \, \, mod\, \, \, m_{1}=a_{1}$

$x\, \, \, mod\, \, \, m_{2}=(k_{1}m_{1}+a_{1})\, \, \, mod\, \, \, m_{2}=(k_{2}m_{2}+a_{2})\, \, \, mod\, \, \, m_{2}=a_{2}$

这样，我们将(*)式整理一下:

$x=m_{1}k_{1}+a_{1}+k\frac{m_{1}m_{2}}{d}=x+k\frac{m_{1}m_{2}}{d}$

这样，我们就把两个式子合并成了一个，其中 $\frac{m_{1}m_{2}}{d}=[m_{1},m_{2}]$ ,表示 $a_{1},a_{2}$ 的最小公倍数。

因而，在新的式子中，把 $x_{0}$ 看成新的式子的 $a$ ， $\frac{m_{1}m_{2}}{d}$ 看成是 $m$ ，即：

$x\equiv x_{0}(mod\, \, \, \, \frac{m_{1}m_{2}}{d})$

无解的条件即为： $(a_{2}-a_{1})%gcb(m_{1},m_{2})\neq 0$

基本步骤解决了，下面有几个问题需要讨论：

 k1=k1*(a2-a1)/d;
 k1 = (k1 % (m2/d) + m2/d) % (m2/d);
这里其实k1有可能是负数的(因为一般扩展欧几里得算法不止一个解叭)，那我觉得如果k1是负数的话，我们的目的是让m1k1对m1m2取余，这样也就是k对m2/d取余，由于k有可能是负数，一般步骤就是先对x取余，然后加上x，再对x取余。

最后在得到x的时候，同样让a对m取余，以防x是负数，加上：

 x = (a1 % m1 + m1) % m1;

二.具体代码

#include<iostream>
using namespace std;
typedef pair<long,long> II;
II exgcb(long a,long b){
    if(b==0) return {1,0};
    else{
        II tmp=exgcb(b,a%b);
        return {tmp.second,tmp.first-a/b*(tmp.second)};
    }
}
int main(){
    int n;
    cin>>n;
    long long  a1,m1;
    cin>>m1>>a1;
    bool flag=true;
    long long x=0;
    for(int i=1;i<=n-1;i++){
        long long a2,m2;
        cin>>m2>>a2;
        II tmp=exgcb(m1,m2);
        long long d=m1*tmp.first+m2*tmp.second;
        if((a2-a1)%d){//d有可能是负的
            flag=false;
            break;
        }
        else{
            long long k1=tmp.first;
            k1=k1*(a2-a1)/d;
            k1 = (k1 % (m2/d) + m2/d) % (m2/d);
            x=m1*k1+a1;
            a1=x%(m1/d*m2);
            m1=m1/d*m2;
        }
    }
    if(flag){
        x = (a1 % m1 + m1) % m1;
        cout<<x;
    }
    else cout<<-1;
}