数论的回顾

扩展欧几里得

简单的欧几里得较为简单，这里就不过多赘述。
扩展欧几里得
已知$a$,$b$的最小公约数为gcd，对于方程$a\ast x+b\ast y=gcd$，已知其特解$x_0$,$y_0$那么就可得到通解为
$x=x_0+\frac{b}{gcd}\cdot t$
$y=y_0-\frac{a}{gcd}\cdot t$

考虑通过欧几里得算法计算$a$和$b$最小公约数达到递归边界时$a_n=gcd$,$b_n=0$那么此时$x_n=1$,$y_n=0$,$a_n\cdot x_n+b_n\cdot y_n=gcd$
再考虑欧几里得算法进行两层递归时
$a_i\cdot x_i+b_i\cdot y_i=gcd$
$b_i\cdot x_{i+1}+(a_i\%b_i)\cdot y_{i+1}=gcd$
展开第二个式子得到
$a_i\cdot y_{i+1}+b_i\cdot (x_{i+1}-a_i/b_i\cdot y_{i+1})=gcd$
联立上述式子则可得到
$x_i=y_{i+1}+1$
$y_i=b_i\cdot (x_{i+1}-a_i/b_i\cdot y_{i+1})$
通过从递归边界开始进行反推则可得到一个$x_0$,$y_0$的特解
在具体做题时判断形如$a\cdot x+b\cdot y=c$的式子是否有解的条件，即判断$c$是否能被$a$和$b$的最小公约数整除
代码如下：

ll exgcd(ll a,ll b,ll &x1,ll &y1){
    if(b==0){
        x1 = 1;
        y1 = 0;
        return a;
    }
    ll ret = exgcd(b, a%b, x1, y1);
    ll temp = x1;
    x1 = y1;
    y1 = temp - a/b*y1;
    return ret;
}

乘法逆元

对于一个$a\times b\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$将b称为a的乘法逆元，记$b$为$inv(a)$

扩展欧几里得求乘法逆元

当$gcd(a,p)=1$时对于$a\times b\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$可转换为$ab-py=1$,如此则可通过$exgcd$求出$a$的乘法逆元b来，注意最后要将$b$转换成最小正整数。

费马小定理求乘法逆元

费马小定理为如果$p$是一个质数，且整数$a$不是$p$的倍数，那么满足式子：
$a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$
简单证明一：考虑到$p$是一个质数，且整数$a$不是$p$的倍数。那么 { a , 2 a , 3 a , 4 a . . . ( p − 1 ) a } \{a,2a,3a,4a...(p-1)a\} {a,2a,3a,4a...(p−1)a}是$p$的一个完全剩余系，故$1\times 2\times 3\times ...\times (p-1)\equiv a\times 2a\times 3a\times ...\times (p-1)a(\mathrm{m}\mathrm{o}\mathrm{d}p)$
即$(p-1)!\equiv (p-1)!a^{p-1}$又因为$gcd((p-1)!,p)=1$所以等式两边可以约去$(p-1)!$，最后得到
$a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$
简单证明二：运用下述欧拉定理即可得证
所以在前提条件为$p$是一个质数，且整数$a$不是$p$的倍数时可以得到$inv(a)=a^{p-2}$

线性递推求乘法逆元

求阶乘的乘法逆元

若已知$inv(n!)$,则由式子
$inv((n-1)!)\cdot (n-1)!\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$
化为
$inv((n-1)!)\cdot (n-1)!\cdot n\cdot inv(n!)\equiv n\cdot inv(n!)(\mathrm{m}\mathrm{o}\mathrm{d}p)$
即可得到递推式
$inv((n-1)!)\equiv n\cdot inv(n!)(\mathrm{m}\mathrm{o}\mathrm{d}p)$
编写代码时可先通过扩展欧几里得或费马小定理得到$inv(n!)$再进行递推

递推的乘法逆元

对于当前的所求元素$i$，考虑其和模数$p$的关系为$p=a\cdot i+b$
那么$a\cdot i+b\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$
等式两边同时乘以$inv(i)$,$inv(b)$得到$a\cdot inv(b)+inv(i)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$
那么$inv(i)\equiv -a\cdot inv(b)(modp)$
最后得到递推式为$inv(i)\equiv -\lfloor \frac{p}{i}\rfloor \cdot inv(p\%i)$

中国剩余定理

考虑这样的一个式子
$\{\begin{array}{l}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ \cdots \\ x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$
中国剩余定理给出了该式子有解的条件以及解的构造方法
若$m_1,m_2\cdots ,m_n$互质，则该方程组有解。
可通过如下方式得到它的解
令 M = ∏ i = 1 n m i , M i = M m i ∀ i ∈ { 1 , 2 , ⋯   , n } M=\prod\limits_{i=1}^nm_i,M_i=\frac{M}{m_i} \forall i\in\{1,2,\cdots,n\} M=i=1∏n​mi​,Mi​=mi​M​∀i∈{1,2,⋯,n}
设$t_i{M}_i\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$即$t_i$是$M_i$在模$m_i$意义下的逆元，同时满足$t_i{M}_i\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{m}_{j,j\ne i})$
则可得到该式子解的通式为$x=a_1{t}_1{M}_1+a_2{t}_2{M}_2+\cdots +a_n{t}_n{M}_n+kM,k\in \mathbb{Z}$
又因为$M_i和m_i$互质，所以一定能用扩展欧几里得算法求出这样的$t_i$

扩展中国剩余定理

扩展中国剩余定理是$m_1,m_2\cdots ,m_n$不互质的情况下的求解，证明也不是很难，暂且占坑吧，写出其合并两个式子时的通式如下
设$gcd=gcd(m_i,m_{i+1})$
$x\equiv (inv(\frac{m1}{gcd},\frac{m2}{gcd})\cdot \frac{a_2-a_1}{gcd})\%\frac{m_2}{gcd}\cdot m_1+a_1(\mathrm{m}\mathrm{o}\mathrm{d}\frac{m_1\cdot m_2}{gcd})$

欧拉定理

欧拉函数

欧拉函数：对于一个正整数n，小于n且和n互质的正整数个数（包括1）的个数，记为$\varphi (n)$。
欧拉函数的通式为$\varphi (n)=n(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_n})$其中$p_1$，$p_2$…$p_n$是n的所有质因数。
$\varphi (n)=n\prod _{i=1}^n(1-\frac{1}{p_i})$

欧拉定理

若$gcd(a,m)=1,$则$a^{\varphi (m)}=1(\mathrm{m}\mathrm{o}\mathrm{d}m)$
所以采用欧拉定理进行降幂

扩展欧拉定理

$a^x\equiv a^{x\mathrm{m}\mathrm{o}\mathrm{d}\varphi (m)+\varphi (m)}(\mathrm{m}\mathrm{o}\mathrm{d}m)$