如何在Pulp CBC模型中检测Presolve不可行性并终止程序

使用Pulp调用CBC求解优化模型时，若输入参数错误，模型中的硬约束会导致不可行。CBC的presolve阶段会输出presolve detects 1 infeasibility，但程序仍会继续执行后续求解步骤，耗费大量时间。当前做法是等待求解完成后检测result_status，再重新运行软约束版本，效率低下。需要在presolve检测到不可行时直接终止程序，避免无意义的求解耗时。

方法1：设置CBC的abortOnInfeasible参数 #

CBC原生支持在presolve发现模型不可行时立即终止求解，只需添加命令行参数abortOnInfeasible yes。在Pulp中可以通过以下方式配置求解器：

import pulp

# 构建你的优化模型
model = pulp.LpProblem("MyOptimizationModel", pulp.LpMinimize)
# ... 添加变量、约束、目标函数 ...

# 配置CBC求解器，添加abortOnInfeasible参数
solver = pulp.PULP_CBC_CMD(options=['abortOnInfeasible yes'])
# 求解模型
model.solve(solver)

# 检查求解状态
if pulp.LpStatus[model.status] == 'Infeasible':
    print("Presolve阶段判定模型不可行，将切换为软约束版本重新求解")
    # ... 此处编写切换软约束的逻辑 ...

设置该参数后，CBC在presolve检测到不可行时会直接终止，不会进入后续迭代，大幅缩短耗时。

方法2：实时监控CBC输出日志（进阶控制） #

如果需要更精细的流程控制，可以通过subprocess直接调用CBC，实时捕获输出日志，一旦检测到presolve不可行的关键词就立即终止进程：

import subprocess
import pulp

# 生成模型的LP文件
model = pulp.LpProblem("MyModel", pulp.LpMinimize)
# ... 构建模型 ...
lp_file_path = "temp_model.lp"
model.writeLP(lp_file_path)

# 构造CBC命令行
cbc_command = [
    "cbc",
    lp_file_path,
    "abortOnInfeasible", "yes",
    "solve",
    "solution", "temp_solution.sol"
]

# 启动CBC进程并实时监控输出
process = subprocess.Popen(
    cbc_command,
    stdout=subprocess.PIPE,
    stderr=subprocess.STDOUT,
    text=True
)

infeasible_flag = False
for line in iter(process.stdout.readline, ''):
    print(line.strip())
    # 检测presolve不可行的关键词
    if "Presolve determined that the problem was infeasible" in line:
        infeasible_flag = True
        process.terminate()  # 立即终止CBC进程
        break

if infeasible_flag:
    print("模型不可行，切换至软约束版本重新求解")
    # ... 执行软约束版本的模型构建与求解 ...
else:
    # 处理可行的求解结果
    print("模型求解完成，读取结果中...")
    # ... 读取并解析solution.sol文件 ...

这种方式可以完全掌控求解流程，在检测到不可行的第一时间终止进程，避免任何不必要的计算。

配置abortOnInfeasible yes后，CBC的输出会在presolve阶段结束后立即终止，示例输出如下：

Welcome to the CBC MILP Solver 
Version: 2.10.3 
Build Date: Dec 15 2019 

command line - temp_model.lp abortOnInfeasible yes solve solution temp_solution.sol
At line 2 NAME          MODEL
At line 3 ROWS
...
Coin0008I MODEL read with 0 errors
Presolve determined that the problem was infeasible with tolerance of 1e-08
Analysis indicates model infeasible or unbounded
1 infeasibilities

此时程序会直接停止，不会进入后续的迭代步骤，耗时仅为presolve阶段的时间，通常远短于完整求解时间。

内容的提问来源于stack exchange，提问作者Bhavya Budhia