【刷题】数学知识——扩展欧几里得、裴蜀定理

裴蜀定理

1 定理

裴蜀定理：

设$a,b$是不全为零的整数，则一定存在整数$x,y$, 使得$ax+by=gcd(a,b)$

换句话说，若$a,b$是不全为零的整数，则$ax+by=gcd(a,b)$一定存在整数解$x,y$。
证明过程在上面的链接中。
求出整数解的过程就是扩展欧几里得算法

2 推论

令$d=gcd(a,b)$，由于$a$和$b$分别去掉最大公因数$d$后，二者不再有公因子，即$\frac{a}{d}$和$\frac{b}{d}$是互质的
令$a^{\prime }=\frac{a}{d}$，$b^{\prime }=\frac{b}{d}$，$ax+by=d$左右同除$d$，得到如下推论：
推论1：

若$a^{\prime }$和$b^{\prime }$互质，则方程$a^{\prime }x+b^{\prime }y=1$有整数解。

此外还可证明推论2：

若方程$a^{\prime }x+b^{\prime }y=1$有整数解，则$a^{\prime }$和$b^{\prime }$互质

反证：若$a^{\prime }$和$b^{\prime }$不互质，则有$a^{\prime }=gcd(a^{\prime },b^{\prime })\ast p$，$b^{\prime }=gcd(a^{\prime },b^{\prime })\ast q$，此处$p,q$为互质的整数。
代入方程中，得到：
$gcd(a^{\prime },b^{\prime })\ast p\ast x+gcd(a^{\prime },b^{\prime })\ast q\ast y=1$
左右同除$gcd(a^{\prime },b^{\prime })$：
$p\ast x+q\ast y=\frac{1}{gcd(a^{\prime },b^{\prime })}$
由于$a^{\prime }$和$b^{\prime }$不互质，则$gcd(a^{\prime },b^{\prime })>1$，$\frac{1}{gcd(a^{\prime },b^{\prime })}$是小数，所以$x$和$y$不可能是整数，方程$a^{\prime }x+b^{\prime }y=1$没有整数解，矛盾。

结合推论1和推论2，可得到裴蜀定理推论：

$a^{\prime }$和$b^{\prime }$互质的充要条件是方程$a^{\prime }x+b^{\prime }y=1$有整数解

3 判断不定方程是否有整数解

$ax+by=z$有整数解的充要条件是$gcd(a,b)|z$

1、若$z$能整除$gcd(a,b)$，则$ax+by=z$有整数解
设$z=k\ast gcd(a,b)$
由裴蜀定理知存在整数解$x^{\prime },y^{\prime }$使得$a{x}^{\prime }+b{y}^{\prime }=gcd(a,b)$
则$ax+by=z$相当于$ak{x}^{\prime }+bk{y}^{\prime }=k\ast gcd(a,b)$
整数解是$x=k{x}^{\prime }$和$y=k{y}^{\prime }$

2、若$ax+by=z$有整数解，则$z$能整除$gcd(a,b)$
证明类似推论2
反证：若$z$不能整除$gcd(a,b)$
将$ax+by=z$两边同除$gcd(a,b)$
有$a^{\prime }x+b^{\prime }y=\frac{z}{gcd(a,b)}$，其中$a^{\prime }=\frac{a}{gcd(a,b)}$和$b^{\prime }=\frac{b}{gcd(a,b)}$是互质的整数
由于$\frac{z}{gcd(a,b)}$是小数，所以$x$和$y$不可能是整数，方程$ax+by=z$没有整数解，矛盾。

4 推广到n元

对不定方程$a_1{x}_1+a_2{x}_2+...+a_n{x}_n=z$，满足$gcd(a_1,a_2,...,a_n)|z$时，方程才有整数解

$gcd(a_1,a_2,...,a_n)=1$的充要条件是不定方程$a_1{x}_1+a_2{x}_2+...+a_n{x}_n=1$有整数解

---

扩展欧几里得

利用辗转相除法可以顺便得到x和y。
辗转相除法过程如下：

int exgcd(int a, int b) {
    if(!b) {
    	return a;
    }
    else {
		return exgcd(b, a % b);
	}
}

当b=0时，可以看出ax+by=gcd(a,b)的一组解是x=1,y=0。于是有：

    if(!b) {
    	x = 1, y = 0;
    	return a;
    }

在进行辗转时，$a^{\prime }=b$，$b^{\prime }=a\%b$，带入$a^{\prime }x+b^{\prime }y=d$中可得$bx+(a\%b)y=d$。
$a\%b$可以写成$a-\lfloor a/b\rfloor \ast b$
所以有$bx+(a-\lfloor a/b\rfloor b)y=d$。
我们想知道a和b的整数解x和y，因此要提出a和b，a和b的系数就是下一轮的x’和y’。
上式变为$ay+b(x-\lfloor a/b\rfloor y)=d$
于是有$x^{\prime }=y,y^{\prime }=(x-\lfloor a/b\rfloor y)$

    else {
		int d = exgcd(b, a % b, x, y);
		int tmp = x;
		x = y;
		y = tmp - a / b * y;
		return d;
	}

为了便于书写，可以交换递归时x和y的位置，省掉中间变量tmp

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

该题代码如下

#include <iostream>
using namespace std;
int n, a, b, x, y;

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

int main() {
	scanf("%d", &n);
	while(n -- ) {
		scanf("%d%d", &a, &b);
		exgcd(a, b, x, y);
		printf("%d %d\n", x, y);
	}
	return 0;
}