同余初步，线性同余

对于整数而言：

模运算的性质

$amodb=a-b\lfloor \frac{a}{b}\rfloor$，对于负数也依然如此。
其中符号$\lfloor \rfloor$表示向下取整。

1. (a%m+b%m)%m=(a+b)%m
   将a、b取模m分为倍数部分和余数部分，无论如何取模，倍数部分都不做贡献。性质显然成立。
2. (a%m×b%m)%m=a×b%m
   取模运算的优先级定义为与乘除运算优先级相当，从左到右依次计算。
   将a、b取模m分为倍数部分和余数部分，无论如何取模，倍数部分都不做贡献。性质显然成立。
   或者将乘法变为若干次加法，也可以说明性质成立。

同余式的定义

对于a $mod$ m=b $mod$ m，我们称a、b在模m的意义下同余，即具有相同的余数。
记作$a\equiv b(modm)$，简写成$a\equiv b$。

以下性质均在同一模数下成立：

1. 由于模运算的性质：$A\equiv a,B\equiv b\Rightarrow A\pm B\equiv a\pm b$
2. 因而有：$A\equiv a\Rightarrow k\cdot A\equiv k\cdot a$
3. 因而有：$A\equiv a,B\equiv b\Rightarrow AB\equiv ab$
   这是由于：
   $A\equiv a\Rightarrow AB\equiv aB$
   $B\equiv b\Rightarrow aB\equiv ab$
   QED.
4. 进而有：$A\equiv a\Rightarrow A^k\equiv a^k$

显然有$a\equiv a$，这说明，同一模数，模意义下可以自由地做加减乘法，也可以移项，移项要变号。

乘法逆元

模意义下是没有除法的，那么如何将$ax\equiv b$这个式子转化成$a\equiv b\cdot y$的形式呢？

我们只需要找到$x^{-1}$，使得$x\cdot x^{-1}\equiv 1$，然后在$ax\equiv b$乘上$x^{-1}$，就会有：$ax\cdot x^{-1}\equiv b{x}^{-1}$，即$a\equiv b{x}^{-1}$

那么对于$x\cdot x^{-1}\equiv 1$，就称$x$与$x^{-1}$互为模$m$时的乘法逆元。

显然逆元在模意义下是唯一的，因为若有$x\cdot x^{-1}\equiv 1,x\cdot x^{\prime -1}\equiv 1,x^{-1}\not\equiv x^{\prime -1}$，则$x\cdot x^{-1}\equiv 1\equiv x\cdot x^{\prime -1}$，则$x^{-1}\equiv x^{\prime -1}$，矛盾。

枚举一下就会发现，$x$模$m$不一定会存在逆元。

计算乘法逆元

剩余系

剩余类（同余类）

给定一个正整数$m$，把所有整数根据模$m$的余数分为$[0,m-1]$这$m$类，每一类$C_r$代表可以被表示为$p\cdot m+r$形式的数。

完全剩余系（完系）

从每个模$m$剩余类中取出一个数字，构成一个大小为m的集合，称为模m的一个剩余系。

简化剩余系（缩系）

我们从每个余数与$m$互质的剩余类中取出一个数字，构成一个大小为$\phi (m)$的集合，称为模m的一个既约剩余系（简化剩余系）。

显然既约剩余系中的每一个数都与$m$互质，考虑$n=p\cdot m+r$，则有：
$\gcd (n,m)=\gcd (n\%m,m)=\gcd (r,m)=1$

欧拉定理

欧拉定理即为：
$a^{\phi (m)}\equiv 1(modm)(a\perp m)$

这意味着：$a^{\phi (m)-1}\cdot a\equiv 1(modm)$
也就是说：$a^{-1}\equiv a^{\phi (m)-1}(modm)$

因此欧拉定理可以用来寻找逆元。

证明一下欧拉定理：
设$a$表示模m的一个既约剩余系。$a_i$表示取到其第i个元素。

首先证明： ∀ a i ∈ a , { a i ⋅ x ∣ x ∈ a } 也是模 m 的一个既约剩余系 \forall a_i\in a,\{a_i\cdot x|x\in a\}也是模m的一个既约剩余系 ∀ai​∈a,{ai​⋅x∣x∈a}也是模m的一个既约剩余系
先证明，集合 { a i ⋅ x ∣ x ∈ a } \{a_i\cdot x|x\in a\} {ai​⋅x∣x∈a}的模意义下的元素是两两不同的：
假设有$a_i\cdot x\equiv a_i\cdot y(modm)(x\not\equiv y)$，那么根据我们下面说到的第三种求逆元的方法的理论，我们知道，必存在$a_i$的逆元，因此$x\equiv y$，矛盾。

接着证明， ∀ X ∈ { a i ⋅ x ∣ x ∈ a } : X ⊥ m \forall X\in\{a_i\cdot x|x\in a\}:X\perp m ∀X∈{ai​⋅x∣x∈a}:X⊥m，也就是集合中的每一个元素都与m是互质的。
这很显然，由于$a_i\perp m,x\perp m$，则两者都与m无公共因子，立即得出$a_i\cdot x\perp m$。

那么我们证明了 ∀ a i ∈ a , { a i ⋅ x ∣ x ∈ a } 也是模 m 的一个既约剩余系 \forall a_i\in a,\{a_i\cdot x|x\in a\}也是模m的一个既约剩余系 ∀ai​∈a,{ai​⋅x∣x∈a}也是模m的一个既约剩余系。
定义这个集合为$S$

则有$\prod a\equiv \prod S(modm)$
$S$的每一个元素都可以提出一个$a_i$，作为既约剩余系，它共有$\phi (m)$个元素：
$\prod a\equiv a_i^{\phi (m)}\prod a(modm)$

仿照上面的证明，立即得知$\prod a\perp m$，所以根据下文所说，我们知道存在逆元，可以除掉：

$a_i^{\phi (m)}\equiv 1(modm)$

$a_i$是模m既约剩余系的任意一个元素，因此：
$\forall x\perp m:x^{\phi (m)}\equiv 1(modm)$

可以先筛出欧拉函数值，然后$x^{\phi (m)}$可以用快速幂求。
时间复杂度O(n)~O(logn)

费马小定理

费马小定理说的是：
$x^{p-1}\equiv 1(modp)(p\in \mathbb{P})$

因此费马小定理也可以用来求逆元：
$x^{p-2}\equiv x^{-1}(modp)$

费马小定理是欧拉定理的特殊情况，带入$\phi (p)=p-1$立即得到费马小定理。

费马小定理常用快速幂求逆元。

时间复杂度O(logn)

同余方程转不定方程，Exgcd求乘法逆元

注意到同余方程$x\cdot x^{-1}\equiv 1(modm)$实际上等于相当于不定方程
$x\cdot x^{-1}+m\cdot y=1$，其中$x,m$是常数，$x^{-1},y$是未知数，均为整数。

那么存在逆元相当于此不定方程有整数解，也即$\gcd (x,m)|1$，则$\gcd (x,m)=1$，即$x\perp m$时存在逆元。

因此可以用exgcd求出$x^{-1}$，即为逆元。

求解一次同余方程

含有未知数的同余式称为同余方程。

对于一元一次同余方程$ax\equiv b(modm)$，有以下两种解法：

1. 同余方程转不定方程：
   此同余方程等价于不定方程$ax+my=b$，可以用exgcd求解。
   同时我们得到了同余方程有解的充分必要条件：$\gcd (a,m)|b$。
2. 乘法逆元：
   若$a\perp m$，则存在$a$在模$m$意义下的乘法逆元，可立即得出$x\equiv b\cdot a^{-1}$。
3. 技巧做法：
   首先注意到$ax\equiv b(\mathrm{mod}b)\iff ax+my=b$
   设$\gcd (a,m)=g$，若$g|b$有解，否则无解。
   因此不定方程两边同时除以$g$：
   $\frac{a}{g}\cdot x+\frac{my}{g}=\frac{b}{g}$
   这意味着同余方程也可以同除以$g$：
   $\frac{a}{g}\cdot x\equiv \frac{b}{g}(\mathrm{mod}\frac{m}{g})$
   （随便说一句，既然可以同时除，那么就可以同时乘对吧.）
   此时$\frac{a}{g}\perp \frac{m}{g}$，除过来：
   $x\equiv \frac{b}{g}\cdot {\left(\frac{a}{g}\right)}^{-1}$
   由于这种方法需要求逆元，因此实用性不足（这要求筛出欧拉函数，或要求$\frac{m}{g}$为质数，但事实上，$\frac{m}{g}$不一定为质数），但是它具有很大的启发性。

这说明单个同余方程的解满足形式：
$x\equiv x_0\left(mod\frac{m}{gcd(a,m)}\right)$
$x_0\equiv \frac{b}{g}\cdot {\left(\frac{a}{g}\right)}^{-1}\left(mod\frac{m}{gcd(a,m)}\right)$

批量求逆元

模m

计算 { 1 , 2 , . . . , m − 2 , m − 1 } \{1,2,...,m-2,m-1\} {1,2,...,m−2,m−1}模m的乘法逆元：

有$k!=k\times (k-1)\times (k-2)...\times 2\times 1$：
$(k-1){!}^{-1}\equiv k{!}^{-1}\cdot k$，$k^{-1}\equiv (k-1)!\times k{!}^{-1}$，只需要求$n{!}^{-1}$

模质数M

计算 { 1 , 2 , . . . , M − 2 , M − 1 } \{1,2,...,M-2,M-1\} {1,2,...,M−2,M−1}模质数M的乘法逆元：

其中$1^{-1}\equiv 1$，对于$x\ge 2$，有：
$\lfloor \frac{M}{x}\rfloor \cdot x+M\%x\equiv M\equiv 0(\mathrm{mod}M)$
则：
$x\equiv -\left(M\%x\right){\lfloor \frac{M}{x}\rfloor }^{-1}(\mathrm{mod}M)$
$x^{-1}\equiv {\left(-\left(M\%x\right){\lfloor \frac{M}{x}\rfloor }^{-1}\right)}^{-1}(\mathrm{mod}M)$

我们知道$(xy{)}^{-1}\equiv x^{-1}{y}^{-1}$，还知道$(-1{)}^{-1}\equiv 1$，还知道${\left(x^{-1}\right)}^{-1}\equiv x$，因此：
$x^{-1}\equiv -{\left(M\%x\right)}^{-1}\lfloor \frac{M}{x}\rfloor (\mathrm{mod}M)$

而$M\%x$的逆元早就求出来了，因此是可以做的。

inv[0]=inv[1]=1;
for(int i=2;i<=N;i++) (inv[i]=(M-1)*(M/i)%M*inv[M%i])%=M;
	

互为逆元

对于$m\perp x$，则存在逆元：
$x\cdot x^{-1}\equiv 1(\mathrm{mod}m)$

因此同余方程转不定方程：
$x\cdot x^{-1}+m\cdot y=1$

最后求出的$x^{-1}$即为$x$模$m$的逆元，那么求出的$y$是什么呢？

不难发现，$y\equiv m^{-1}(\mathrm{mod}x)$

这个现象本身就十分有趣了，且在一种中国剩余定理的方法内会用到。

事实上，这还有个高级版本：
$x\cdot x^{-1}+y\cdot y^{-1}=c$

则：
$x\cdot x^{-1}\equiv c(\mathrm{mod}y)$
$y\cdot y^{-1}\equiv c(\mathrm{mod}x)$
$x\cdot x^{-1}\equiv c(\mathrm{mod}{y}^{-1})$
$y\cdot y^{-1}\equiv c(\mathrm{mod}{x}^{-1})$

后记

于是皆大欢喜。