python解一次不定方程linear Diophantus Equation

量化情绪

于 2021-07-14 19:24:50 发布

阅读量464

点赞数

分类专栏：递归文章标签： python 算法

本文链接：https://blog.csdn.net/m0_46321344/article/details/118734974

版权

递归专栏收录该内容

1 篇文章

订阅专栏

本文通过Python实现了解一次不定方程的线性丢番图方程（linearDiophantusEquation）的递归算法。首先定义了gcd函数计算最大公约数，接着利用欧几里得算法求解a'x+b'y=1的x和y，最后通过LDE2函数结合gcd和LDE解决ax+by=c的整数解问题。示例代码中展示了如何求解23x-76y=3的解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

python解一次不定方程linear Diophantus Equation

话不多说，上代码：

def gcd(a,b): #最大公约
    r=a%b
    if r==0:
        return b
    else:
        return gcd(b,r)
def LDE(a,b):#先求a‘x+b’y=1 （a',b')=1
    if b==1:
        return 0,1
    x2, y2 = LDE(b, a % b)
    x = y2
    y = x2 - a // b * y2
    return x, y
def LDE2(a,b,c): #ax+by=c 返回x,y
    d=gcd(a,b)
    #print(d)
    if c%d!=0:
        print("无整数解")
        return
    a2=a/d;b2=b/d
    x3,y3=LDE(a2,b2)
    return x3*c/d,y3*c/d
print(LDE2(23,-76,3)) #(-99.0, -30.0)