用欧几里德扩展求解不定方程的整数解:
上回说到辗转相除法也就是欧几里德算法:链接
https://blog.csdn.net/qq_39652552/article/details/84073880
其中有个重要的性质 gcd(a,b) = gcd( b,a%b)
我们利用其性质做一下扩展,就可以用它做不定方程的求解。
引理:存在x,y 使得 ax + by = gcd(a,b)
利用贝祖定理证明:
如果a = 0 = b 结论显然成立。
如果b = 0(a = 0同理) 因为 gcd(a,b) = a 所以 ax + by = a 一定存在一个解,x = 1, y = 0
当a!= 0 && b!= 0
设 ax1 + by1 = gcd(a,b) 因为 gcd(a,b) = gcd(b,a%b)
所以ax1 + by1 = gcd(b,a%b)又因为 bx2 + (a%b)y2 = gcd(b,a%b)
所以ax1 + by1 = bx2 + (a%b)y2 a%b可以表示为 (a - k * b)( k = a /b(计算机里的/) )
则 ax1 + by1 = bx2 + (a - k * b)y2 将等式右边展开得到 ax1 + by1 = bx2 + ay2 - kby2
于是ax1 + by1 = ay2 + b(x2 - ky2)
所以 x1 = y2 y1 = (x2 - ky2),因此,每一层的x和y可以由下一层推出来
当b = 0时是最后一组解,递归下去算得最后一层的解,再层层递归回来就能算出来第一层的x和y
代码如下:
int qiu( int a, int b, int & x, int & y)
{
int t, r;
if( b == 0) //递归基
{
x = 1;
y = 0;
return a;
}
else
{
r = qiu( b, a%b , x , y );
t = x;
x = y;
y = t - a/b * y;
}
return r;
}
拿a = 24 b = 16 做一下示范吧
首先第一层是 24 * x1 + 16 *y1 = gcd(24,16)
二: 16 * x2 + 8 * y1 =gcd(16,8)
三: 8 * x3 = gcd (8,0)
由第三层得到x3 = 1 , y=0 然后由这两个下面这两个式子向上推,算出x1和y1
x1 = y2
y1 = (x2 - ky2)【注意,k是等于a/b哦】