辗转相除法（欧几里得算法）和扩展欧几里得算法实现及证明

辗转相除法（欧几里得算法）和扩展欧几里得算法实现及证明

今天看了好长时间的数论知识点，学完之后，过了一个假期，再回头看扩展欧几里得时已经是

懵 懵 懵

于是我上网找到了它的证明！

（叙述的可能会非常啰嗦详细）

---

辗转相除法：

定义及用途：

辗转相除法， 又名欧几里德算法，是求最大公约数的一种方法。它的具体做法是：用较小数除较大数，再用出现的余数（第一余数）去除除数，再用出现的余数（第二余数）去除第一余数，如此反复，直到最后余数是0为止。如果是求两个数的最大公约数，那么最后的除数就是这两个数的最大公约数。
——摘自百度百科

证明：

我们想要证明$a$和$b$的最大公约数和$b$和$a\ mod\ b$的最大公约数相等，即$gcd(a,b)=gcd(b,a\ mod\ b)$

可以令$a = kb+r$,则$r = a-kb$,其中$k = \left \lfloor a/b \right \rfloor$,（$\left \lfloor x \right \rfloor$表示向下取整），那么$r = a\ mod\ b$。

• 若$a$,$b$的一个公约数为$d$，则$a$，$b$都可以被$d$整除，同时$r$也可被$d$整除（这里可以自己想一下）。所以$d$是$b$和$a\ mod\ b$的公约数。

• 若$b$和$a\ mod\ b$的一个公约数为$d$，则$b$,$r$可以被$d$整除，同时$a = kb+r$,所以$d$也是$a$和$b$的公约数。

综上，$a,b$和$b,a\ mod\ b$的公约数是一样的，其最大公约数也必然相等。

实现：

辗转相除法是我们很早的时候就了解的算法，实现也有很多种方法，可以用递归或非递归的方式实现，这里就不再给出

---

扩展欧几里得算法：

定义及用途：

（就是对欧几里得算法的扩展）

扩展欧几里得算法可以用来求解线性同余方程和乘法的逆元。

线性同余方程的解法步骤这里就不再赘述，扩展欧几里得算法所能实现的就是：给定$a$，$b$，求解$x$,$u$,使之满足$ax+by = gcd(a, b) =d$。接下来给出过程及证明

证明：

已知：$ax+by=1$且$a$，$b$互质

求：$x$，$y$

$ax+by=1=gcd(a,b)$

根据朴素的辗转相除法得：$gcd(a,b) = gcd(b,a\ mod\ b)$

所以：$ax+by=gcd(b,a\ mod\ b)=bx_{1}+(a\ mod\ b)y_{1}$

又因为：$a\ mod\ b=a-\left \lfloor a/b \right \rfloor*b$

所以展开并移项得：

$ax+by=bx_{1}+(a-\left \lfloor a/b \right \rfloor*b)y_{1}$

$ax+by=bx_{1}+ay_{1}-b*\left \lfloor a/b \right \rfloor y_{1}$

$ax+by=ay_{1}+b(x_{1}-\left \lfloor a/b \right \rfloor *y1)$

所以：$x=y_{1}$,$y=x_{1}-\left \lfloor a/b \right \rfloor *y1$，这就是$x$，$y$与$x_{1}$，$y_{1}$的关系

像这样递归下去，直到$b = 0$，此时$a_{n}x_{n}+0*y_{n} = gcd(a_{n},0) = a_{n}$，那么$x_{n} = 1$,$y_{n} = 任意数$，为了方便，我们取$y = 0$.

之后进行回溯，$x=y_{1}$,$y=x_{1}-\left \lfloor a/b \right \rfloor *y1$，最后就做到了在求解$a$，$b$的最大公约数的过程中也实现了对$x$，$y$的求解，这也是这个算法叫做扩展欧几里得的原因。

实现：

网上扩展欧几里得算法的模板很多，这里给出一个：

int exgcd(int a,int b)
{
    if(b==0)//边界
    {
        x=1;
        y=0;//y可以为任意数，这里取0
        return a;
    }
    int r=exgcd(b,a%b);
    int t=x;
    x=y;//更新解的过程
    y=t-a/b*y;//更新解的过程
    return r;
}

---

其他：

这里再推荐两个博客，也是有关证明的，可以看一下

https://blog.csdn.net/lincifer/article/details/49391175

http://www.cnblogs.com/frog112111/archive/2012/08/19/2646012.html

---

结束了。

有不严谨的地方还请指正。