python - 使用 Python cvxpy 和 Python_MIP() 处理不可行的解决方案（CBC 求解器）

Question

在使用 cvxpy 和 Python MIP (CBC) 作为求解器的 MIP 工具中，我正在寻找一种方法来生成用于记录目的的不可行解决方案的细节。有谁知道这是否存在？我正在阅读 CPLEX 具有这种能力（功能）来指出特定的约束违规。

Answer 1

在 python docplex 中你可以得到放松和冲突。

参见示例https://github.com/AlexFleischerParis/zoodocplex/blob/master/zoorelaxationandconflict.py

因此，在您的情况下，您可以：

• 在 python docplex 中打开你的模型
• 保留 cvxpy，导出 mps 文件，然后在 docplex 中导入 mps 文件（带有 model.import_model），然后在 docplex 中运行松弛

这里的模型：

from docplex.mp.model import Model
from docplex.mp.relaxer import Relaxer

from docplex.mp.conflict_refiner import ConflictRefiner

mdl = Model(name='buses')
nbbus40 = mdl.integer_var(name='nbBus40')
nbbus30 = mdl.integer_var(name='nbBus30')

mdl.add_constraint(nbbus40*40 + nbbus30*30 >= 300, 'kids')
mdl.add_constraint(nbbus40 + nbbus30 <= 7, 'maxTotalBuses')

mdl.minimize(nbbus40*500 + nbbus30*400)

mdl.solve()

mdl.report()
print(f"* solve status is: '{mdl.solve_details.status}'") #infeasible model

print()
print("------- starting relaxation")
print()

rx = Relaxer()
rx.relax(mdl)

print ("number_of_relaxations= " + str(rx.number_of_relaxations))
rx.print_information()

mdl.report()
print(f"* status after relaxation is: {mdl.solve_details.status}")
#print(mdl.solution)

print()
print("------ starting conflict refiner")
print()

cr=ConflictRefiner()
conflicts=cr.refine_conflict(mdl)
conflicts.display()