扩展欧几里得

扩展欧几里得
贝祖说：ax+by = gcd(a, b) 这个公式肯定有解
额(⊙﹏⊙) 贝祖说的对！
这个公式我们写作ax+by = d，(gcd(a, b) | d)
扩展欧几里得就是求这个公式的解的方法
老样子，先上代码：

#include<cstdio>
typedef long long LL;
void extend_Eulid(LL a, LL b, LL &x, LL &y, LL &d){//注意这里的x,y,d
    if (!b) {d = a, x = 1, y = 0;}//细心的同学可能发现了这里已经把辗转相除法包含了
    else{
        extend_Eulid(b, a % b, y, x, d);
        y -= x * (a / b);
    }
}
int main(){
    LL a, b, d, x, y;
    while(~scanf("%lld%lld", &a, &b)){
        extend_Eulid(a, b, x, y, d);
        printf("%lld*a + %lld*b = %lld\n", x, y, d);
    }
}

为什么能这么解呢？
大佬给出了证明：https://blog.csdn.net/destiny1507/article/details/81750874

初步想法：由于是递归的算法，如果我们知道了这一层和上一层的关系，一层一层推下去，就可以推到最开始的。类似数学上的数学归纳法。

假设当前我们在求的时a和b的最大公约数，而我们已经求出了下一个状态：b和a%b的最大公因数，并且求出了一组x1和y1使得 b*x1+(a%b)y1=gcd*
（注意在递归算法中，永远都是先得到下面一个状态的值）

这时我们可以试着去寻找这两个相邻状态的关系：

首先我们知道：a%b=a-(a/b)*b；带入：

b*x1 + (a-(a/b)*b)*y1

= bx1 + ay1 – (a/b)by1

= ay1 + b(x1 – a/by1) = gcd 发现 x = y1 , y = x1 – a/by1

这样我们就得到了每两个相邻状态的x和y的转化，就可以在求gcd的同时对x和y进行求值了

大佬说的已经很详细了，这里给的代码是把计算和x与y交换合到一块了