如何解决使用 Python cvxpy 和 Python_MIP() 处理不可行的解决方案CBC 求解器
在使用 cvxpy 和 Python MIP (CBC) 作为求解器的 MIP 工具中,我正在寻找一种方法来生成用于记录目的的不可行解决方案的细节。有谁知道这是否存在?我读到 CPLEX 具有这种能力(功能)来指出特定的约束违规。
解决方法
在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()
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。