#数论# 欧几里德算法 、扩展欧几里德算法 、费马小、逆元求解（ing）

欧几里德求gcd（辗转相除法）：

定理：

gcd(a, b) = gcd(b, a % b)
两个正整数a和b（a>b），它们的最大公约数等于a除以b的余数c和b之间的最大公约数

证明：

a可以表示成a = kb + r，则r = a % b；

1. 假设d是a, b的一个公约数，则有a % d = 0，b % d = 0，由于r = a - kb，因此r % d = 0，证明充分性；
2. 假设d 是(b, a % b)的公约数，则b % d = 0，r % d = 0，由于a = kb + r，因此a % d = 0，证明必要性；

实现代码：

int gcd(int a, int b) { //递归形式
    if(b == 0) return a;  //当 b == 0 时，（a，0）的最大公约数是a
    return gcd(b, a % b);
}

int gcd(int a, int b) { //非递归形式
    int tmp = a;
    while(b) {
        a = b;
        b = tmp % b;
    }
    return a;
}

复杂度：O( log(max(a, b)) )

扩展欧几里德算法：

用来在已知a, b求解一组x，y使得 a * x + b * y = Gcd(a, b)，扩展欧几里德常用在求解横线性方程及方程组中。

核心代码：

//if(a < 0) a = -a, c = -c;  
int exgcd(int a, int b, int &x, int &y) { //a必须大于0
    if(b == 0) {
        x = 1, y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= (a / b) * x;