同余的一些知识点整理

同余

定义

数论讲义上给的定义是

给定一个正整数 $m$ ，如果用 $m$ 去除两个整数 $a$ 和 $b$ 所得到的余数相同，我们就说 $a,b$ 对模数 $m$ 同余，记作 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$.

同时给出 $a,b$ 对模数 $m$ 同余的充要条件是

$m|(a-b)$

欧几里得算法

欧几里得算法，即辗转相除法，用于计算两个数的最大公约数 $\gcd$。

公式 $\gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b),b\ne 0$

同时，若 $a\mathrm{m}\mathrm{o}\mathrm{d}b=0$ , $b$ 即为最大公约数

int gcd(int a, int b){
    return b ? gcd(b, a % b) : b;
}

裴蜀定理

百度百科的定义如下

若 $a,b$ 是整数，且 $\gcd (a,b)=d$，那对任意的整数 $x,y,ax+by$ 都一定是 $d$ 的倍数，特别地，一定存在整数 $x,y$， 使得 $ax+by=\gcd (a,b)$

也就是说二元一次不定方程 $ax+by=c$ 有解的充要条件是 $\gcd (a,b)|c$

扩展欧几里得算法

根据裴蜀定理有
$$\begin{gathered} ax+by=\gcd (a,b) \\ b{x}_0+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_0=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b) \end{gathered}$$

同时根据欧几里得算法有
$$\begin{gathered} \gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b) \\ ax+by=b{x}_0+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_0 \\ ax+by=b{x}_0+(a-\lfloor \frac{a}{b}\rfloor \ast b)y_0 \\ ax+by=a{y}_0+b\ast (x_0-\lfloor \frac{a}{b}\rfloor \ast y_0) \end{gathered}$$
因此我们得到了 $x,y$ 的一组解
$$\{\begin{array}{r}x=y_0\\ y=x_0-\lfloor \frac{a}{b}\rfloor \ast y_0\end{array}$$

特别地，当 $b=0$ ，显然有 $x_0=1,y_0=0$

因此我们可以在 $\gcd$ 的同时递归更新 $x,y$

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

以上程序给出了二元一次不定方程 $ax+by=\gcd (a,b)$ 的一组解。但实际应用中一般是给你$a,b,c$，让你求
$$ax+by=c$$
的最小整数解 $x$

根据裴蜀定理，我们知道二元一次不定方程 $ax+by=c$ 有解的充要条件是 $\gcd (a,b)|c$

对于 $\gcd (a,b)\nmid c$，显然是无解。

设 $d=\gcd (a,b)$，假设我们已经求出$ax+by=d$ 的一组解 $x_0,y_0$，两边同乘 $\frac{c}{d}$，得
$$a{x}_0\ast \frac{c}{d}+b{y}_0\ast \frac{c}{d}=d\ast \frac{c}{d}=c$$
因此我们得到一组解
$$\{\begin{array}{r}{x}^{\prime }=x_0\ast \frac{c}{d}\\ y^{\prime }=y_0\ast \frac{c}{d}\end{array}$$

当然这个 $x^{\prime }$ 显然不是最小整数解，考虑将 $x^{\prime }$ 减 $k\ast \frac{b}{d}$ 。为保证等式成立，$y^{\prime }$ 应该加上 $k\ast \frac{a}{d}$

因此我们得到了方程的一般解
$$\{\begin{array}{r}x=x_0\ast \frac{c}{d}-k\ast \frac{b}{d}\\ y=y_0\ast \frac{c}{d}+k\ast \frac{a}{d}\end{array}$$
而对于的最小整数解 $x$ ，我们可以通过取模运算来得到

x *= abs(c / d);
x = (x % abs(b / d) + abs(b / d)) % abs(b / d); 

扩展欧几里得算法解同余方程

给定同余方程
$$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$$
求最小整数解 $x$

根据同余的充要条件，有
$$b|(ax-1)$$
不妨设
$$ax-1=b{y}^{\prime }$$

移项得
$$ax-b{y}^{\prime }=1$$

设
$$y=-y^{\prime }$$
因此原题就是求二元一次不定方程 $ax+by=1$ 的解。

逆元

给定正整数 $a$ 和 $m$ ，如果有 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$， 那么把这个同余方程中 $x$ 的最小整数解叫做 $a$ 模 $m$ 的逆元，其中这个 $x$ 记作 $a^{-1}$ ，其中逆元又被称为数论倒数。

求解逆元的方法即为上述扩展欧几里得算法解同余方程。

特别地，当 $m$ 为质数，根据费马小定理有 $x=a^{m-2}$

中国剩余定理（孙子定理）

给定同余方程组
$$\{\begin{array}{r}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ x\equiv a_3(\mathrm{m}\mathrm{o}\mathrm{d}{m}_3)\\ \cdots \\ x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$$
且 $m_i$ 两两互质，则方程有唯一解。

设
$$\begin{gathered} M=\prod _{i=1}^n{m}_i \\ {M}_i=\frac{M}{m_i} \\ {M}_i{t}_i\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i) \end{gathered}$$
则唯一解
$$x=(\sum _{i=1}^n{a}_i{M}_i{t}_i)\mathrm{m}\mathrm{o}\mathrm{d}M$$
说明一点，这个解是构造出来的，不是求出来的。因此不给出证明，读者可带入自行验证即可领会构造原理。