扩展欧几里德算法证明及例题

ax + by = gcd(a,b) = gcd(b,a%b) = b$x_1$ + (a - $\frac{a}{b}b)y_1=a{y}_1+b(x_1-\frac{a}{b}{y}_1$)

考虑gcd的调用顺序 有 x = $y_1$, y = ($x_1$ - $\frac{a}{b}{y}_1$)

当b == 0 时，则有 x = 1,y = 0;

故有代码为

int exgcd(int a,int b,int &x,int &y){
	if(b == 0){
		x = 1,y = 0;
		return a;
	}
	int d = exgcd(b,a%b,y,x);// 这里直接交换了x,y
	y -= a/b * x;
	return d;
}

性质1：
假设求出了一组合法的解。(x,y) 对于另外的一组解 (m,n)

ax + by = am + bn -> a(x-m) = b(y-n)

设g = gcd(a,b) $a_1=\frac{a}{g}$ $b_1=\frac{b}{g}$

$a_1$(x-m) = $b_1$(y-n)

因为$a_1和b_1$ 互质，故(x-m) = $k{b}_1$ 同时解除(y-n) = $-k{a}_1$

故(m,n) = ($x-k{b}_1$,$y+k{a}_1$)

性质2：

$ax+by=c$ 若 c | gcd(a,b) 则有解 反之，无解。

例题：
同余方程
有理数取余
青蛙的约会

待更