扩展欧几里德算法解二元一次不定方程

Originum

于 2018-08-06 22:09:47 发布

阅读量2.3k

点赞数 5

分类专栏： ACM 数论扩展欧几里德理论分析文章标签：扩展欧几里德算法不定方程 ACM 数论二元一次方程

本文链接：https://blog.csdn.net/Originum/article/details/81437027

版权

扩展欧几里德算法：

已知两个不完全为 0 的非负整数 a，b，必然存在整数对 x，y ，使它们满足贝祖等式：

$ax+by = gcd(a, b) =d$

解一定存在，根据数论中的相关定理。下面给出代码：

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

解二元一次不定方程(ax + by = c)：

扩展欧几里德算法很大的一个用处就是在解不定方程。在计算时一般先是求出一对特解，再根据x与y的变化比例，构造通解。

注意我们求的所有解，都是整数解。

我们先计算 $ax+by = gcd(a, b)$ 的解，等式两边同时乘以

$c / gcd(a,b)$

原式子就被构造成

$a x(c / gcd(a,b) )+by(c / gcd(a,b) ) = c$

所以就求出了不定方程的一对特解

最低0.47元/天解锁文章