扩展欧几里德算法是用来在已知a, b求解一组x,y,使它们满足贝祖等式: ax + by = gcd(a, b) = d(解一定存在,根据数论中的相关定理)。扩展欧几里德常用在求解模线性方程及方程组中。
对于两个不全为0的整数a、b,必存在一组解x,y,
使得ax+by==gcd(a,b)
由定理:ax+by==gcd(a,b),若b==0,则gcd(a,0)==a。
原式变为ax+by==a --> x==1,y==0。
若x,y表示第一次递归时的值,x1,y1表示第二次递归时的值。
则gcd(a,b)==gcd(b,a%b),有ax+by==b*x1+(a%b)*y1。
b*x1+(a%b)*y1==b*x1+(a-(a/b)*b)*y1==a*y1+b*(x1-(a/b)*y1),
最终得到ax+by==a*y1+b*(x1-(a/b)*y1)
即上一深度的x等于下一深度的y1,上一深度的y等于下一深度的x1-(a/b)*y1。
对于ax+by==c,x1=x*(c/gcd(a,b)),y1=y*(c/gcd(a,b))
//如果GDC(a,b) = d,则存在x,y,使d = ax + by
//extended_euclid(a,b) = ax + by
int extended_euclid(int a, int b, int &x, int &y)
{
int d;
if (b == 0) { x = 1; y = 0; return a; }
d = extended_euclid(b, a%b, y, x);
y -= a / b*x;
return d;
}
拓展欧几里得的几个应用
用来在已知a, b求解一组x,y使得ax+by = Gcd(a, b) =d(解一定存在,根据数论中的相关定理)。扩展欧几里德常用在求解模线性方程及方程组中.