欧几里得与扩展欧几里得

欧几里德算法

定义 $\gcd (a,b)$ 为 $a,b$ 的最大公约数

假设 $\gcd (a,b)\equiv \gcd (b,a\%b)$

设 $k=a\%b$，设 $k^{\prime }=\gcd (a,b)$

则 $a=x{k}^{\prime },b=y{k}^{\prime }$， { x , y } ∈ [ 1 , ∞ ] \{x,y\}\in[1,∞] {x,y}∈[1,∞]，$x\perp y$

$$\therefore k=a\%b=a-k^{\prime \prime }b=x{k}^{\prime }-k^{\prime \prime }y{k}^{\prime }=(x-k^{\prime \prime }y)k^{\prime }$$

$$\therefore k=(x-k^{\prime \prime }y)k^{\prime }$$

$$\because b=y{k}^{\prime }$$

假设 $(x-k^{\prime \prime }y)$ 与 $y$ 不互质

设 $y=pq,(x-k^{\prime \prime }y)=p^{\prime }q$

将 $y$ 代入 $(x-k^{\prime \prime }y)$

$$\therefore x-k^{\prime \prime }y=x-pq{k}^{\prime \prime }=p^{\prime }q$$

$$\therefore x-pq{k}^{\prime \prime }=p^{\prime }q$$

$$\therefore x=p^{\prime }q+pq{k}^{\prime \prime }$$

提公因式

$$x=(p^{\prime }+p{k}^{\prime \prime })q$$

$$\because a=x{k}^{\prime }=(p^{\prime }+p{k}^{\prime \prime })q{k}^{\prime }$$

$$\because b=y{k}^{\prime }=pq{k}^{\prime }$$

$$\therefore \gcd (a,b)=q{k}^{\prime }$$

与原命题不符， $\therefore (x-k^{\prime \prime }y)\perp y$

$$\because \gcd (b,a\%b)=\gcd (b,k)=\gcd (y{k}^{\prime },(x-k^{\prime \prime }y)k^{\prime })$$

$\because (x-k^{\prime \prime }y)\perp y$

$\therefore \gcd (b,a\%b)\equiv k^{\prime }\equiv \gcd (a,b)$

当 $a\%b=0,\gcd (a,b)=b$

所以就有了代码

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

对于斐波那契数列

$$\gcd (f(n),f(m))=f(\gcd (n,m))$$

扩展欧几里德

假设存在 { x , y } , \{x,y\}, {x,y}, 使得 $\gcd (a,b)\equiv ax+by$

当 $b=0$ 时

$\gcd (a,b)=a,x=1,y=0$

当 $b\ne 0$ 时

设 $a{x}^{\prime }+b{y}^{\prime }=\gcd (a,b)$

设 $b{x}^{\prime \prime }+(a\%b)y^{\prime \prime }=\gcd (b,a\%b)$

$$\because a\%b\equiv a-\lfloor \frac{a}{b}\rfloor \times b$$

$$\therefore a{x}^{\prime }+b{y}^{\prime }=b{x}^{\prime \prime }+(a-\lfloor \frac{a}{b}\rfloor \times b)y^{\prime \prime }$$

$$\therefore a{x}^{\prime }+b{y}^{\prime }=b{x}^{\prime \prime }+a{y}^{\prime \prime }-\lfloor \frac{a}{b}\rfloor b{y}^{\prime \prime }$$

$$\therefore a{x}^{\prime }+b{y}^{\prime }=a{y}^{\prime \prime }+b(x^{\prime \prime }-\lfloor \frac{a}{b}\rfloor y^{\prime \prime })$$

$$\therefore x^{\prime }=y^{\prime \prime },y^{\prime }=(x^{\prime \prime }-\lfloor \frac{a}{b}\rfloor y^{\prime \prime })$$

$∵ $ 当 $b=0$ 时存在解

$\therefore$ 解存在

然后我们就有代码了

inline void exgcd(int a,int b) {
    if(b) {
        exgcd(b,a%b);
        int t=x; x=y;
        y=k-a/b*y;
    } else y=(x=1)-1;
}

贝祖定理

定义： $ax+by=c$ ，满足 { x , y } ∈ [ 1 , ∞ ] \{x,y\}\in[1,∞] {x,y}∈[1,∞]，且 $\gcd (a,b)|c$

设 $s=\gcd (a,b)$，则 $s|a$ 且 $s|b$

且 $s|ax$ 且 $s|by$

设 $k^{\prime }\in [1,\infty ]$

则 $c=\gcd (a,b)\times k^{\prime }$

$\therefore \gcd (a,b)|c$ 成立

扩展欧几里德解不定方程

对于不定方程 $ax+by=c$，设 $k^{\prime }=\gcd (a,b)$

$\because$ 存在 { x , y } \{x,y\} {x,y} 使得 $\gcd (a,b)\equiv ax+by\equiv k^{\prime }$

当 $c\%k^{\prime }\ne 0$ 无解

当 $c\%k^{\prime }=0$

可用 $\text{exgcd}$ 算出 $\gcd (a,b)=ax+by$，即 $ax+by=k^{\prime }$ 的一组解

解得 $x=x_0,y=y_0$

$\therefore x_0\times \frac{c}{k^{\prime }},y_0\times \frac{c}{k^{\prime }}$

$\therefore x=x_0\times \frac{c}{k^{\prime }}\times k^{\prime \prime },y=y_0\times \frac{c}{k^{\prime }}\times k^{\prime \prime }$

$$\therefore \frac{ac{k}^{\prime \prime }{x}_0}{k^{\prime }}+\frac{bc{k}^{\prime \prime }{y}_0}{k^{\prime }}=c$$

$$\therefore ac{k}^{\prime \prime }{x}_0+bc{k}^{\prime \prime }{y}_0=c{k}^{\prime }$$

$$\therefore (a{x}_0+b{y}_0)c{k}^{\prime \prime }=c{k}^{\prime }$$

$$\therefore c{k}^{\prime }{k}^{\prime \prime }=c{k}^{\prime }$$

$$\therefore k^{\prime \prime }=1$$

$$\therefore x=x_0\times \frac{c}{k^{\prime }},y=y_0\times \frac{c}{k^{\prime }}$$

$$\therefore a{x}_0\times \frac{c}{k^{\prime }}+b{y}_0\times \frac{c}{k^{\prime }}=c$$

$$\therefore (a{x}_0+b{y}_0)\times \frac{c}{k^{\prime }}=c$$

$$\therefore a{x}_0+b{y}_0=k^{\prime }$$

所以此式恒等成立

扩展欧几里德解同余方程

对于同余方程 $ax\equiv c\mathrm{m}\mathrm{o}\mathrm{d}b$

$$\therefore ax-c\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}b$$

$$\therefore ax+by\equiv c$$

所以，只要求 $\text{exgcd}$ 即可

乘法逆元

通常用于求

$$\frac{a}{b}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

若

$$a\times x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

且 $a$ 与 $x$ 互质，则可称 $x$ 为 $a$ 的逆元，也称 $x$ 为 $a$ 的倒数

对于

$$\frac{a}{b}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

只需求出

$$(\frac{1}{b}(\mathrm{m}\mathrm{o}\mathrm{d}p))\times a$$

也就是求出 $b$ 的逆元再乘上 $a$

拓展欧几里德解乘法逆元

利用拓展欧几里德解同余方程 $ax\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 时 $c=1$ 的情况

也就是转换为 $ax+by\equiv 1$ 的情况

然后用 $\text{exgcd}$ 求解

快速幂解乘法逆元

通过费马小定理可知，$\forall a\in [1,\infty ]$ 与素数 $p$

有 $a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

有 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

所以有 $a^{p-1}\equiv ax(\mathrm{m}\mathrm{o}\mathrm{d}p)$

所以有 $a^{p-2}\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}p)$

只要用快速幂求出 $a^{p-2}(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的值就行了

线性乘法逆元

对于求一连串数字 $(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的逆元，上面的方法就显得慢

考虑逆元之间的递推关系

首先我们有 $1^{-1}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$

不妨设 $p=k\times i+r$，也就是 $k=\lfloor \frac{p}{i}\rfloor$

则 $k\times i+r\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$

乘上 $i^{-1}\times r^{-1}$

得

$$k\times r^{-1}+i^{-1}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

$$i^{-1}\equiv -k\times r^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

$$i^{-1}\equiv -\lfloor \frac{p}{i}\rfloor \times r^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

$$i^{-1}\equiv -\lfloor \frac{p}{i}\rfloor \times (p(\mathrm{m}\mathrm{o}\mathrm{d}i){)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

然后我们就有代码了

ans[1]=1;
for(int i=1;i<p;i++) ans[i]=(p-p/i)*ans[p%i]%p;

阶乘逆元

我们有

$$ans[i+1]=\frac{1}{(i+1)!}$$

$$ans[i+1]\times (i+1)=\frac{i+1}{(i+1)!}=\frac{1}{i}=ans[i]$$

$$ans[i]=ans[i+1]\times (i+1)$$

所以就可以从后往前求出所有逆元