同余原理与丢番图方程

同余定理

定理： gcd(a, b) = gcd(b, a%b)
证明：
设$a=kb+r,$则r=a%b
假设d是（a,b）的一个公约数，则有d|a,d|b
而$r=a-kb,$因此$d|r,$因此$d$是（b,a%b）的公约数
假设d是（b,a%b）的公约数，则有$d|b,d|r$
而$a=kb+r,$因此$d|a,$因此$d$是$（a,b）$的公约数
因此（a,b）和（b,a%b）的公约数是相同的，其最大公约数也必然相同

int gcd(a,b){return !b?a:gcd(b,a%b);}

更相减损术

若$a<b,$ $swap(a,b)$
若$a$为奇数，$b$为偶数，$GCD（a,b）=GCD(a,b/2)$
若$a$为偶数，$b$为奇数，$GCD（a,b）=GCD(a/2,b)$
若$a$为偶数，$b$为偶数，$GCD（a,b）=GCD(a/2,b/2)$
若$a$为奇数，$b$为奇数，$GCD（a,b）=GCD(a-b,b)$

丢番图方程

裴蜀定理：丢番图方程$ax+by=c$有解当且仅当 $gcd(a,b)|c$
扩展gcd算法：求丢番图方程$ax+by=c$的整数解
推导：
$a{x}_1$ + $b{y}_1$ = $gcd(a,b)$
$b{x}_2$ + $(a$%b)$y_2$ = $gcd$(b,a%b)
由$gcd(a,b)$ = $gcd$(b, a%b)得到：
$a{x}_1$ + $b{y}_1$ = $b{x}_2$ + $(a$%b)$y_2$ = $a{y}_2$ + $b(x2-(a/b)y_2$)
所以
$x_1$ = $y_2$
$y_1$ = $x_2$ - $(a/b)y_2$

递归实现方程解

void exgcd(int a, int b, int &x, int &y) {
	if(b == 0) {x = 1; y = 0; return;}
	exgcd(b, a%b, x, y);
	int t = x; x = y; y = t - a/b*y;
}

上面求得其中一组解
（$x_1$+$\frac{kb}{gcd(a,b)}$,$y_1$+$\frac{ka}{gcd(a,b)}$）

是该方程的通解

证明：代入易证

相关题目

SDOI2009 最大公约数