扩展欧几里得算法详解

欧几里得法（辗转相除法）

在扩展欧几里得之前先讲一下欧几里得法。
定理：设$a=qb+r$，其中$a，b，q，r$，都是正整数,且，则有：
        $gcd(a,b)=gcd(b,r)$
证明：
只需证明$a$与$b$和$b$与$r$有相同的最大公因子即可。设$d$为$a$与$b$的 最大公约数，那么$d|a$且$d|b$。又因为$r=qb-a$，且$d|a，d|b，$所以$d|r$。此时$a$与$b$的最大公约数为$d$，又因为$d|$r&且$d|b$，即$r$与$b$的最大公约数也为$d$，所以上式成立。

辗转相除法：设整数$a，b，$且$b\ne 0$，做带余除法$a=qb+r$，若$r=0$时，$gcd(a,b)=b$。否则$a^{\prime }=b，b^{\prime }=r$；做带余除法$a^{\prime }=q{b}^{\prime }+r^{\prime }$,直到$r=0$时，$gcd(a,b)=b^{\prime }$。由上面定理可知该法正确性。

拉梅定理：用辗转相除法计算两个正整数的最大公因子时，所需的除法次数不会超过两个整数较小的那个数的倍数的5倍。

算法实现：

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

有时也可以直接使用c语言的库函数__gcd()。

扩展欧几里得算法

定义：设a和b不全为0，则存在整数x，y，使得：
       $gcd(a,b)=ax+by$
证明：
记$a=r_0,b=r_1$，对$a,b$做辗转相除法，即：
  $r_i=q_{i+1}{r}_{i+1}+r_{i+2},i=0,1,...,k-2$
  $r_{k-1}=q_k{r}_k$
此时$gcd(a,b)=r_k$,把上式改写为
  $r_i=r_{i-2}-q_{i-1}{r}_{i+1},i=2,3,...,k$
一直往前代入，$r_k$即可表示为$r_0$和$r_1$的线性组合，即$a$和$b$的线性组合,得证上式。

由扩展欧几里得可知整数$a$与$b$互质的充分必要条件：存在整数$x$和$y$，使得$xa+yb=1$

扩展欧几里得的算法实现由递归实现，递归边界为$b==0，a=gcd(a,b)$，这时有$a\times 1+b\times 0=gcd(a,b)$，但此时的0，1并不是需要求出的$x，y$，只要找到每次递归$x，y$变化的关系就可以递归求出$x，y$的值。
我们假设$x_{i-1}\times a+y_{i-1}\times b=gcd(a,b)$，那么$x_i\times b+y_i\times a\%b=gcd(b,a\%b)$,$x_i$为第$i$次递归$x$的值，$x_0$即为所求的$x$。
因为$a\%b=a-(\frac{a}{b})\times b$，代入可得：
$x_i\times b+y_i\times (a-(\frac{a}{b})\times b)=gcd(b,a\%b)$
$b\times (x_i-\frac{a}{b}\times y_i)+y_i\times a=gcd(a,b)$
所以$x_{i-1}=y_i,y_{i-1}=x_i-\frac{a}{b}\times y1$

算法实现：

int exgcd(int a,int b,int&x,int&y)
{
    if(b==0)
    {
        x=1,y=0;
        return a;
    }
    else 
    {
        int r=gcdEx(b,a%b,x,y); //r=GCD(a,b)=GCD(b, a%b)
        int t=x;
        x=y;
        y=t-a/b*y ;
        return r ;
    }
}

扩展欧几里得算法的应用主要有三方面：
1.求解不定方程。
2.求解模的逆元。
3.求解同余方程。

这里简单介绍以下求解模的逆元：
如果$a\times x\equiv 1(modp)$,且$gcd(a,p)=1$，即$a$与$p$互质，则称$a$关于模$p$的乘法逆元为$x$。
设$a\times x=1(modp)$，即：
$a\times x-m\times p=1$ ($m$为整数)
其中gcd(a,p)=1，即求解$x，p$，设$b=p，y=m$，即求解
  $gcd(a,b)=ax+by$
求出$x$即为$a$关于模$p$的乘法逆元。