扩展欧几里得

扩展欧几里得是用来解决线性同余方程问题。

1. 什么是线性同余方程?

给定三个正整数a, b, d求下面方程的非零整数解。
$$ax+by=d$$

2. 如何求解？

首先，可以观察到这个方程是不一定有解的，比如下面这个方程是一定没有整数解的。
$$2x+4y=1$$
那么什么时候有整数解呢？有一个裴蜀定理：线性同余方程有解的充要条件是: $$d=gcd(a,b)$$ 且此gcd(a, b)是最小的整数d,其他的d均是gcd(a,b)的整数倍。
这个证明可以使用扩展欧几里得来做，因为扩展欧几里得提供了一种找到这个方程非零整数解的办法。

3. 扩展欧几里得过程
   $d=gcd(a,b)$
   $ax+by=d$
   $by+(a\mathrm{m}\mathrm{o}\mathrm{d}b)x=d$
   $x=x,y=y-(a/b)\ast x$

4. 代码

#include <bits/stdc++.h>
using namespace std;

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);
    y -= a / b * x;
    return d;
}