Python递归死循环问题-有关数学建模整数规划

zfcq

代码是解决数学建模中整数规划问题的代码

1. import math
   
2. from scipy.optimize import linprog
   
3. import sys
   
4. 
   
5. def integerPro(c, A, b, Aeq, beq, t=1.0E-12):
   
6.     res = linprog(c, A_ub = A, b_ub=b, A_eq=Aeq, b_eq=beq)
   
7.     if(type(res.x) is float):
   
8.         bestX = [sys.maxsize]*len(c)
   
9.     else:
   
10.         bestX = res.x
    
11.     bestVal = sum([x*y for x,y in zip(c,bestX)])
    
12.     if all(((x-math.floor(x))<t or (math.ceil(x)-x)<t ) for x in bestX):
    
13.   
    
14.         return (bestVal, bestX)
    
15.    
    
16.     else:
    
17.         ind = [i for i,x in enumerate(bestX) if (x-math.floor(x))>t and (math.ceil(x) - x)>t][0]
    
18.         newCon1 = [0]*len(A[0])
    
19.         newCon2 = [0]*len(A[0])
    
20.         newCon1[ind] = -1
    
21.         newCon2[ind] = 1
    
22.         newA1 = A.copy()
    
23.         newA2 = A.copy()
    
24.         newA1.append(newCon1)
    
25.         newA2.append(newCon2)
    
26.         newB1 = b.copy()
    
27.         newB2 = b.copy()
    
28.         newB1.append(-math.ceil(bestX[ind]))
    
29.         newB2.append(math.floor(bestX[ind]))
    
30.         r1 = integerPro(c, newA1, newB1, Aeq, beq)
    
31.         r2 = integerPro(c, newA2, newB2, Aeq, beq)
    
32.         if r1[0] < r2[0]:
    
33.             return r1
    
34.         else:
    
35.             return r2
    
36. 
    
37.         
    
38. if __name__ == '__main__':
    
39.     c = [3,4,1]
    
40.     A = [[-1,-6, -2],[-2,0,0]]
    
41.     b = [-5, -3]
    
42.     Aeq = [[0,0,0]]
    
43.     beq=[0]
    
44.     print(integerPro(c, A, b, Aeq, beq))
    

复制代码

运行结果

运行结果应为：
(8.0, array([2, 0., 2.]))

主要是求解决整数规划-分支定界相关问题的方法，可以发有点相关解决方法，或者直接用其他库。。

suchocolate

递归教学实验可以，生产环境尽量不用。
上面是递归超限了，所以error了，把代码改成不用递归。

学渣李某人

非要递归的话开头加上:

1. sys.setrecursionlimit(1000000)
   

复制代码

zfcq

suchocolate 发表于 2021-7-22 10:06
递归教学实验可以，生产环境尽量不用。
上面是递归超限了，所以error了，把代码改成不用递归。

那不就成了线性规划问题了吗，这不用递归的话，整数最优解难道要自己一点点试？

kogawananari

把代码改成尾递归

zfcq

学渣李某人 发表于 2021-7-22 10:13
非要递归的话开头加上:

不输出结果，直接 RESTART: Shell 了

学渣李某人

别用idle或者别用递归, 昨天有人问过这个问题

学渣李某人

你这个问题和这个很像:
https://fishc.com.cn/forum.php?mod=viewthread&tid=199343

zfcq

学渣李某人 发表于 2021-7-22 10:23
你这个问题和这个很像:
https://fishc.com.cn/forum.php?mod=viewthread&tid=199343

我这个应该是可以运行成功的，因为有确切结果，我怕我是代码逻辑问题

学渣李某人

zfcq 发表于 2021-7-22 10:28
我这个应该是可以运行成功的，因为有确切结果，我怕我是代码逻辑问题

你先试试用别的编辑器运行一下吧, 实在不行可以用vscode打断点调试

3236654291

我也不知道，不过悬赏40鱼币还是我头一次见