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

14 篇文章 0 订阅
10 篇文章 0 订阅

扩展欧几里德算法:

已知两个不完全为 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

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

X=x*c / gcd(a,b)\, ,\, Y=y*c / gcd(a,b)

这里不存在 c / gcd(a,b) 不是整数的情况,因为当 c\, \, mod\, \, gcd(a,b) \neq 0 时,方程不存在整数解,不在考虑范围内。

证明:

a=k1*gcd(a,b)\, ,\, b=k2*gcd(a,b)

ax+by=gcd(a,b)(x*k1+y*k2)=c

若 c\, \, mod\, \, gcd(a,b) \neq 0,则 x*k1+y*k2 非整数

求一组特解代码:

bool Indefinite_equation(int a, int b, int n, int& x, int& y) {    //求一组特解
    int gcd = extgcd(a, b, x, y);       
    if (n % gcd) return false;
    int k = n / gcd;
    x *= k;
    y *= k;
    return true;
}

现在我们来考虑构造通解( 特解是X、Y ):

由于已经满足 aX+bY=c , 两个未知数的关系类似于加权和固定,若一个增大,则另一个一定按比例减小,显然满足关系:

a(X+nb)+b(Y-na)=c

X 每变化 b 的整数倍,Y 就反向变化 a 的同等倍数。似乎构造出了通解表达式 X+nb 、Y-na 。

但是我们会发现一个问题,就是 a 、b 不一定是最小的步长,可能会 “跳过” 许多解。

为了求得最小步长,我们应该对 a 、b 同时除以某个整数,使得商也是整数,就求出了每次最小的变化量。很明显,应除以

gcd(a, b)

所以,通解应该是

 X+n*b/gcd(a,b)\, \, ,\, \, Y-n*a/gcd(a,b)

求特殊的解:

最大/最小解

很多情况下,我们需要求的是 x 或 y 的最小正整数解。这里以 x 的最小解为例:

有一种方案是在求出特解的基础上,进行一个循环

  • 当特解大于零时:  while (x + b/gcd(a,b) > 0) ,x每次减 b/gcd(a,b)
  • 当特解不大于零时: while (x <= 0) ,x每次加 b/gcd(a,b)

但是这个方法存在效率问题,因为 b/gcd(a,b) 可能远小于 x , 达到退出循环条件需要花很长时间。这里我采用的是取模的方法

  • 当特解大于零时: x\, \, mod\, \, b/gcd(a,b)  便是小正整数解
  • 当特解不大于零时:x\, \, mod\, \, b/gcd(a,b) +b/gcd(a,b)  便是小正整数解

x 与 y 最接近的解

以求 x < y时 , x 与 y 最接近的解为例。通过 y - x 的值来考虑。设 x , y 的步长分别是 dx , dy ,则 y - x 的最短步长为 dx + dy。

  • 当特解 X < Y 时,\left \lfloor (y-x)/(dx+dy) \right \rfloor  便是通解中 n 的值
  • 当特解 X > Y 时,   \left \lfloor (y-x)/(dx+dy) \right \rfloor-1 便是通解中 n 的值。但若 (y - x) \, mod\, (dx + dy)=0 时,不需要减 1 

附上代码:

//求 x < y 时,x, y 最接近的解
bool Indefinite_equation(int a, int b, int n, ll& x, ll& y) {      
    int gcd = extgcd(a, b, x, y);
    if (n % gcd) return false;                         //无整数解
    int k = n / gcd, dx = b / gcd, dy = a / gcd;       //dx,dy分别为对应的最小变化值
    x *= k;
    y *= k;                                            //此为某特解
    ll cnt = (y - x) / (dx + dy);
    if (cnt < 0 && (y - x) % (dx + dy)) cnt --;        //处理x>y的情况
    x += dx * cnt;
    y -= dy * cnt;
    return true;
}

 

  • 5
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值