本文介绍了使用Python cvxpy和PYTHON_MIP()(CBC解算器)处理不可行解的处理方法，对大家解决问题具有一定的参考价值，需要的朋友们下面随着小编来一起学习吧！
问题描述
在使用cvxpy和Python MIP(CBC)作为解算器的MIP工具中，我正在寻找一种方法来生成用于日志记录目的的不可行解的细节。有没有人知道这是否存在？我读到CPLEX有这个能力(功能)来指出特定的约束违规。
推荐答案
在pythondocplex中，您可以获得放松和冲突。
参见示例https://github.com/AlexFleischerParis/zoodocplex/blob/master/zoorelaxationandconflict.py
因此，在您的情况下，您可以：
- 在python docplex中转换模型
- 保存cvxpy，导出MPS文件，然后在docplex中导入MPS文件(带有mod.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()
这篇关于使用Python cvxpy和PYTHON_MIP()(CBC解算器)处理不可行解的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持编程学习网！
本站部分内容来源互联网,如果有图片或者内容侵犯您的权益请联系我们删除！