扩展欧几里得
贝祖说: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交换合到一块了