线性丢番图方程(拓展gcd)

二元线性丢番图方程

方程形如ax + by = c(a,b,c已知整数,x,y为变量,求整数解)

几何意义上,只要这条直线上有整数坐标点就有无数解

定理:

记gcd(a,b)为d. 如果d%c != 0 无解,不然有无数解

记一个特解为x₀,y₀,则通解为x = x₀ + (b/d)n, y = y₀ - (a/d)n其中n为任意整数

如25x + 15y = 70

一个特解为x₀ = 4, y₀ = -2

通解为x = 4 +3n, y = -2 - 5n

求特解(拓展gcd)

ax + by = gcd(a,b),如果用拓展gcd可以求出gcd(a,b)的同时求出x和y

就是在gcd的过程中记录一下x和y,大概的推导过程

ax₀ + by₀ = gcd(a,b)

bx₁ + (a%b)y₁ = gcd(b, a%b)

得ax₀ + by₀ = bx₁ + (a%b)y₁

ax₀ + by₀ = bx₁ + (a - a/b * b) * y₁ (这里a/b是取整的哦)

化简得x₀ = y₁, y₀ = x₁ - (a/b)*y₁

我们得出了前后式的关系,gcd最深层的时候a是1,b是0,我们可以倒序逆推得到xy

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;
}

了解了拓展gcd之后,回到二元线性丢番图

现求ax + by = c

一:判断gcd(a,b) % c

二:用拓展gcd求特解x0,y0

三:

ax₀ + by₀ = d

两边同乘c/d得到(c/d)(ax₀ + by₀) = c

四:

ax + by = c与三对照得出特解

x₀’ = x₀ * c/d

y₀’ = y₀ * c/d

五:

通解得出x = x₀’ + (b/d)n

y = y₀’ - (a/d)n

如果要求最小正整数解的话,如果有y = x + a*t， 那么y的最小正整数y’可以这么求：

y’ = (x % t + t) % t;原因是因为t可能为负数,+t取正