欧几里得算法

注意：本篇文章会涉及到较多的数论的符号，如果有些符号您不认识，可以跳转到数论之基础一文查看某符号的定义

$\cdot 欧几里得算法$

$1、gcd(a,b)=gcd(b,a\%b)$

即$(a,b)$ 与$(b,a\%b)$的最大公约数相同，这是辗转相除法的基本原理；

$2、lcm(a,b)=\frac{a\times b}{gcd(a,b)}$

$lcm$是最大公倍数

为了避免 $a\times b$ 爆$longlong$，最好写成$a/gcd(a,b)\ast b$的形式；

$3、$辗转相除法的代码实现：

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

$\cdot 拓展欧几里得算法$

前置知识：裴蜀定理

$1、$一般用于求解$ax+by=gcd(a,b)$的解。

​ 当$b=0$时，$ax+by=a$，故而$x=1,y=0$.

​ 当$b\ne 0$时：
$$\begin{array}{r}\because gcd(a,b)=gcd(b,a\%b)\\ 则:b{x}^{{}^{\prime }}+(a\%b)\times y^{{}^{\prime }}=gcd(b,a\%b)\\ b{x}^{{}^{\prime }}+(a-\lfloor a/b\rfloor )\times y^{{}^{\prime }}=gcd(b,a\%b)\\ a{y}^{{}^{\prime }}+b\times (x^{{}^{\prime }}-\lfloor a/b\rfloor \times y^{{}^{\prime }})=gcd(b,a\%b)\\ \therefore x=y^{{}^{\prime }},y=x^{{}^{\prime }}-\lfloor a/b\rfloor \ast y^{{}^{\prime }}\end{array}$$

则有C++代码实现如下：

int exgcd(int a, int b, int &x, int &y)
{
    if(b == 0) {
		x = 1, y = 0;
        return a;
    }
    int x1, y1, gcd;
    gcd = exgcd(b, a%b, x1, y1);
    x = y1, y = x1 - a/b*y1;  //C++中的两个整数相除自动会下取整
    return gcd;
}

另一种证明方法及代码实现：

当$b=0$时，$ax+by=a$，故而$x=1,y=0$.

当$b\ne 0$时：
$$\begin{array}{r}\because gcd(a,b)=gcd(b,a\%b)=gcd(a\%b,b)\\ b{y}^{{}^{\prime }}+a\%b\times x^{{}^{\prime }}=gcd(a\%b,b)\\ b{y}^{{}^{\prime }}+x^{{}^{\prime }}(a-\lfloor a/b\rfloor \times b)=gcd(a\%b,b)\\ a{x}^{{}^{\prime }}+b(y^{{}^{\prime }}-\lfloor a/b\rfloor \times x^{{}^{\prime }})=gcd(a\%b,b)\\ \therefore x=x^{{}^{\prime }},y=y^{{}^{\prime }}-\lfloor a/b\rfloor \times x^{{}^{\prime }}\end{array}$$
则有C++代码实现如下：

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

我个人喜好第二种证明方式，因为比较好理解，代码也相对而言比较容易实现。

$2、$对于求解更一般的方程$ax+by=c$：

设$d=gcd(a,b)$，则其有解当且仅当$d|c.$

求解方法：

用拓展欧几里得求出$a{x}_0+b{y}_0=d$的解.

则$a(x_0\ast c/d)+b(y_0\ast c/d)=c$

故特解为 $x^{{}^{\prime }}=x_0\ast c/d,y^{{}^{\prime }}=y_0\ast c/d$.

而$通解=特解+齐次解$，齐次解即为$ax+by=0$的解

故而通解为 $x=x^{{}^{\prime }}+k\ast b/d,y=y^{{}^{\prime }}-k\ast a/d,k\in z$

若令$t=b/d,$则对于$x$的最小非负整数解为$(x^{{}^{\prime }}\%t+t)\%t$

然后将$x$代入原式，即可得到$y$的解.

$3、$应用：求解一次同余方程$ax\equiv b(\mathrm{mod}m)$

则等价于求：
$$\begin{gathered} ax=m\times (-y)+b \\ ax+my=b \end{gathered}$$
有解的条件为$gcd(a,m)|b$，然后用拓展欧几里得求解即可.

特别的，当$b=1$且$a$与$m$互质时，则所求$x$即为$a$的逆元.

原创文章，如需转载，请标明出处。