如何获取PuLP求解LP后得到的约束条件对应值？

嘿，我刚好能帮你解决这个问题！PuLP本身确实没有直接提供获取约束条件计算值的内置方法，但咱们可以手动计算出每个约束左边表达式在最优解下的实际值，操作起来很简单。

步骤1：先构建一个示例LP问题 #

首先咱们先写一个简单的线性规划例子，方便后续演示：

import pulp

# 创建最小化问题
prob = pulp.LpProblem("Sample_LP", pulp.LpMinimize)

# 定义变量
x = pulp.LpVariable('x', lowBound=0)
y = pulp.LpVariable('y', lowBound=0)

# 添加目标函数
prob += 3*x + 4*y, "Total_Cost"

# 添加约束（这里给约束命名方便后续识别）
prob += 2*x + y <= 10, "Resource_A_Limit"
prob += x + 3*y <= 12, "Resource_B_Limit"
prob += x >= 2, "Minimum_x_Requirement"

# 求解（关闭日志输出）
prob.solve(pulp.PULP_CBC_CMD(msg=False))

步骤2：计算每个约束的实际值 #

求解完成后，我们可以遍历问题中的所有约束，提取约束的表达式，然后代入变量的最优值计算结果：

# 遍历所有约束并计算其实际值
for constraint_name, constraint in prob.constraints.items():
    # 获取约束左侧的线性表达式
    constraint_expr = constraint.expr
    # 计算表达式在最优解下的数值
    computed_value = pulp.value(constraint_expr)
    # 获取约束的右侧值和约束类型
    rhs_value = constraint.rhs
    sense = constraint.sense

    # 打印结果
    print(f"约束名称: {constraint_name}")
    print(f"  约束左侧计算值: {computed_value:.2f}")
    print(f"  约束右侧值: {rhs_value:.2f}")
    print(f"  约束类型: {sense}")
    print("---")

运行这段代码后，你会看到类似这样的输出：

约束名称: Resource_A_Limit
  约束左侧计算值: 10.00
  约束右侧值: 10.00
  约束类型: <=
---
约束名称: Resource_B_Limit
  约束左侧计算值: 9.60
  约束右侧值: 12.00
  约束类型: <=
---
约束名称: Minimum_x_Requirement
  约束左侧计算值: 2.00
  约束右侧值: 2.00
  约束类型: >=

步骤3：单独获取特定约束的计算值 #

如果你只需要某一个特定约束的计算值，可以直接通过约束名称访问：

# 获取指定约束的计算值
target_constraint = prob.constraints["Resource_A_Limit"]
target_value = pulp.value(target_constraint.expr)
print(f"Resource_A_Limit的计算值: {target_value:.2f}")

原理说明 #

PuLP的每个约束对象都包含expr属性，这个属性存储了约束左侧的线性表达式（比如2*x + y）。当问题求解完成后，变量的最优值已经被赋值，此时调用pulp.value(expr)会自动代入变量的最优值，计算出表达式的实际结果。

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