欧几里得算法和扩展欧几里得算法

超级详细的基础算法和数据结构合集：
https://blog.csdn.net/GD_ONE/article/details/104061907

摘要

本文主要讲解欧几里得算法和扩展欧几里得算法。

欧几里得算法

欧几里得算法就是辗转相除法，用于求两个数的最大公约数。

设$gcd(a,b)$ 表示a和b的最大公约数。
辗转相除法的核心就是$gcd(a,b)=gcd(b,a\%b)$，当$a\%b==0$时，$b$就是最大公约数。
将以上式子转化为递归代码就是欧几里得算法。

欧几里得算法代码：

public static int gcd(int a, int b){
	if(b == 0){ // 当b==0,说明上一层的a%b == 0,所以直接返回a.
		return a;
	}
	return gcd(b, a%b);
}

那么什么是扩展欧几里得算法呢？

扩展欧几里得算法

扩展欧几里得算法是对欧几里得算法的扩展，用于求解二元一次方程的整数解。
欧几里得算法是求公约数的，跟二元一次方程能扯上什么关系？
这个关系就是 裴蜀定理：

若a,b是整数,且gcd(a,b)=d，那么对于任意的整数x,y,ax+by都一定是d的倍数，特别地，一定存在整数x,y，使ax+by=d成立

简单来说就是，如果两个整数的最大公约数是d， 那么方程$ax+by=d$一定有解。这是已经被前人证明过的定理，所以我们直接拿来用就好了，不必纠结证明过程。

那么我们知道了$ax+by=d$这个式子后怎么求解$x$和$y$呢。
$x和y$有很多组解，我们可以知道一组特解是：当$a\%b==0$时，$a==d$,此时令$x=1,y=0$，得到： a * 1 + 0 = d， 显然，x = 1， y = 0,是ax + by = d的一组特解。
那么当$a\%b!=0$时怎么办？
因为$gcd(a,b)==gcd(b,a\%b)$，它们之间是存在一定联系的，所以我们不妨试试能否利用这个式子得到通解。
已知：
$a\ast x_1+b\ast y_1=gcd(a,b)$
$b\ast x_2+(a\%b)\ast y_2=gcd(b,a\%b)$
$gcd(a,b)=gcd(b,a\%b)$;

显然：
$a\ast x_1+b\ast y_1=b\ast x_2+(a\%b)\ast y_2$， 接下来想办法消掉$\%$号，统一为常见的乘法，以便合并同类项。
因为：$a\%b=a-b\ast (a/b)$
所以： $a\ast x_1+b\ast y_1=b\ast x_2+(a-b\ast (a/b))\ast y_2$
将上式展开：$a\ast x_1+b\ast y_1=b\ast x_2+a\ast y2-b\ast (a/b)\ast y2$
将上式整理：$a\ast x_1+b\ast y_1=a\ast y2+b\ast (x2-(a/b)\ast y2)$

显然： $x1=y2,y1=x2-(a/b)\ast y2$
于是我们就得到了解之间的关系，接下来通过一组特解，就可以得到其他解。

我们知道，当$a\%b==0$时，$ax+by=d$的解是$x=1,y=0$。
而求解$gcd(a,b)$时, 当$a\%b==0$时，是递归出口，之后递归函数便会一层一层返回，我们正好可以利用这个性质，求解每一层$a\ast x+b\ast y=d$的解。

代码：

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

因为这是一个递归函数，所以当执行到$b==0$时，函数开始一层层返回，最底层是$a=d,x=1,y=0$,我们可以在函数返回的过程中依次求出每一层的解，当函数执行完毕时， 我们就可以得到原式的一个特解。
何时无解：
根据定义，对于$a\ast x+b\ast y=d，如果d是gcd(a,b)$的倍数，则该式有解，所以如果$d$,则该式无解。

最小正整数解

$a\ast x+b\ast y=d$是有很多组解的，我们怎么得到x的最小整数解呢？

设$(x_1，y_1),(x_2,y_2)是a\ast x+b\ast y=d的两组解$
则有: $a\ast x_1+b\ast y_1=a\ast x_2+b\ast y_2$
整理：$a\ast (x_1-x_2)=b\ast (y_2-y_1)$

等式两边同时除以gcd(a,b)：
得到：$a/gcd(a,b)\ast (x_1-x_2)=b/gcd(a,b)\ast (y_2-y_1)$

设$a^{\prime }=a/gcd(a,b),b^{\prime }=b/gcd(a,b)$
显然：$a^{\prime }与b^{\prime }互质$
所以对于：$a^{\prime }\ast (x_1-x_2)=b^{\prime }\ast (y_2-y_1)$，$(x_1-x_2)是b^{\prime }的整倍数,(y_2-y_1)是a^{\prime }的整倍数$
就好比$3\ast 10=2\ast 15$
所以： $x_1=x_2+b^{\prime }\ast k(k=1,2,3,...)$
即$x的值是每b^{\prime }一循环，所以直接让x\%b^{\prime }$就可以得到最小的解了，但求出的解可能是负数，所以先让$x\%b$，然后再让$x$加上一个$b^{\prime }$,然后再对$b^{\prime }$取余。
即： $x=(x\%b^{\prime }+b^{\prime })\%b^{\prime }$

模板题:
扩展欧几里得算法
代码：

import java.io.*;
import java.util.*;

public class Main{
    
    static BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
    static BufferedWriter out = new BufferedWriter(new OutputStreamWriter(System.out));
    public static int Int(String s){return Integer.parseInt(s);}
    static int x, y;
    public static int exgcd(int a, int b){
        if(b == 0){
            x = 1;
            y = 0;
            return a;
        }
        int d = exgcd(b, a%b);
        int t = y;
        y = x - a/b*y;
        x = t;
        return d;
    }
    public static void main(String[] args) throws Exception{
        int n = Int(in.readLine());
        while(n --> 0){
            String[] s = in.readLine().split(" ");
            int a = Int(s[0]);
            int b = Int(s[1]);
            exgcd(a, b);
            out.write(x+" "+y+"\n");
        }
        out.flush();
    }
} 

例题：
线性同余方程

代码：

import java.io.*;
import java.util.*;

public class Main{
    
    static BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
    static BufferedWriter out = new BufferedWriter(new OutputStreamWriter(System.out));
    public static int Int(String s){return Integer.parseInt(s);}
    static int x, y;
    public static int exgcd(int a, int b){
        if(b == 0){
            x = 1;
            y = 0;
            return a;
        }
        int d = exgcd(b, a%b);
        int t = y;
        y = x - a/b*y;
        x = t;
        return d;
    }
    public static void main(String[] args) throws Exception{
        int n = Int(in.readLine());
        while(n --> 0){
            String[] s = in.readLine().split(" ");
            int a = Int(s[0]);
            int b = Int(s[1]);
            int m = Int(s[2]);
            int d = exgcd(a, m);
            int c = b/d;
            int t = m/d;
            if(b % d != 0){
                out.write("impossible\n");
            }
            else out.write(((long)x*c)%t + "\n");
        }
        out.flush();
    }
}