PAO.Pyomo的MIBS求解器未正确满足z[0]≠z[1]约束问题

你的问题核心是使用PAO的MIBS求解器时，上层模型的z[0]≠z[1]约束未被正确执行，得到了明显违反约束的解。主要原因在于双层模型结构不完整：你定义了SubModel但未给下层模型添加任何目标或约束，MIBS作为双层优化求解器，无法正确处理这种空下层的场景，导致约束检查失效。

解决方案1：补全双层模型结构 #

给下层模型添加一个dummy目标或约束，让求解器能正确识别双层优化结构：

import pyomo.environ as pyo
from pyomo.environ import *
from pao.pyomo import *

N = 2
M=1024

model = ConcreteModel()
model.z = Var(range(N), domain=Integers)
for i in range(N):
    model.z[i].setlb(0)
    model.z[i].setub(N-1)

model.delta = Var(range(N), range(N), domain=Binary)

model.unique_z = ConstraintList()
for i in range(N):
    for j in range(i+1, N):
        model.unique_z.add(model.z[i] - model.z[j] + M*(1-model.delta[i, j]) >= 1)
        model.unique_z.add(model.z[i] - model.z[j] - M*model.delta[i, j] <= -1)

model.obj = Objective(expr=sum(model.z[i] for i in range(N)), sense=maximize)
model.L = SubModel(fixed=model.z)
# 添加dummy目标，补全下层模型
model.L.obj = Objective(expr=0)

with Solver('pao.pyomo.MIBS') as solver:
    # 启用日志查看求解过程
    results = solver.solve(model, tee=True)

print("obj =", model.obj.expr(), "z =", [model.z[i].value for i in range(N)])

解决方案2：简化约束定义 #

对于0-1变量的互斥约束，直接用线性等式替代大M约束，更简洁且避免潜在的求解器处理问题：

import pyomo.environ as pyo
from pyomo.environ import *
from pao.pyomo import *

N = 2

model = ConcreteModel()
model.z = Var(range(N), domain=Binary)  # 直接用Binary域替代整数+上下界

# 直接定义z[0]≠z[1]的约束：两者之和为1
model.unique_z = Constraint(expr=model.z[0] + model.z[1] == 1)

model.obj = Objective(expr=sum(model.z[i] for i in range(N)), sense=maximize)
model.L = SubModel(fixed=model.z)
model.L.obj = Objective(expr=0)

with Solver('pao.pyomo.MIBS') as solver:
    results = solver.solve(model, tee=True)

print("obj =", pyo.value(model.obj), "z =", [model.z[i].value for i in range(N)])

额外调试建议 #

• 检查约束满足情况：求解后手动验证约束值，确认解的可行性：

for idx, c in enumerate(model.unique_z, 1):
    body_val = pyo.value(c.body)
    print(f"约束{idx}: {body_val}，是否满足: {c.lower <= body_val <= c.upper}")

• 查看求解日志：启用tee=True打印求解过程日志，排查求解器是否识别到约束、是否存在可行解等问题。

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