algorithm - 求解线性丢番图方程的算法是什么: ax + by = c

Question

我在这里寻找整数解决方案。我知道它从第一对解和 gcd(a,b)|c 派生出无穷多个解。但是，我们如何找到第一对解决方案呢？有解决这个问题的算法吗？

谢谢，
陈

Answer 1

请注意，并非总有解决方案。事实上，只有当 c 是 gcd(a, b) 的倍数时，才有解决方案。

也就是说，您可以使用 extended euclidean algorithm为此。

这是一个实现它的 C++ 函数，假设 c = gcd(a, b)。我更喜欢使用递归算法:

function extended_gcd(a, b)
    if a mod b = 0
        return {0, 1}
    else
        {x, y} := extended_gcd(b, a mod b)
        return {y, x-(y*(a div b))}

int ExtendedGcd(int a, int b, int &x, int &y)
{
    if (a % b == 0)
    {
        x = 0;
        y = 1;
        return b;
    }

    int newx, newy;
    int ret = ExtendedGcd(b, a % b, newx, newy);

    x = newy;
    y = newx - newy * (a / b);
    return ret;
}

现在，如果您有 c = k*gcd(a, b) 且 k > 0，则等式变为:

ax + by = k*gcd(a, b) (1)
(a / k)x + (b / k)y = gcd(a, b) (2)

因此只需找到 (2) 的解决方案，或者找到 (1) 的解决方案并将 x 和 y 乘以 k .