求解线性丢番图方程

时间:2015-12-16 23:36:11

标签: python equation linear diophantine

我正在尝试用Python编写算法来解决线性丢番图方程。

我认为我的算法是正确的,因为我已经在纸上测试了它,但是当我运行它时,它会返回奇怪的值。

我的代码:

def solve_Dioph(a,b,c):
    m1=1
    m2=0
    n1=0
    n2=1
    r1=a
    r2=b
    while r1%r2!=0:
        q=r1/r2
        aux=r1%r2
        r1=r2
        r2=aux
        aux3=n1-(n2*q)
        aux2=m1-(m2*q)
        m1=m2
        n1=n2
        m2=aux2
        n2=aux3
    return m2*c,n2*c;

它使用7个变量和3个辅助变量。用笔和纸测试后,使用以下值:

a=65 b=14 c=4

我得到了

m2=-3*4 and n2=14*4

然而,当我运行它时:

solve_Dioph(65,14,4)

它返回:

(-227/9, 16393/126)

1 个答案:

答案 0 :(得分:2)

我猜你正在使用Python 3(而不是Python 2)。在Python 3中,整数除法返回浮点结果。但是,您可以使用String fileName = Environment.getExternalStorageDiretory() + "/FolderCreatedByUser/filename.ext"; 运算符强制执行Python 2行为。尝试更改以下行:

//

要:

q=r1/r2