PySCIPOpt实现Benders分解调用getObjVal触发阶段错误排查

问题场景 #

• 实现逻辑：基于PySCIPOpt约束处理器开发Benders分解算法，约束处理器执行时通过model.getSolVal(sol, var)获取主问题变量取值传入子问题，子问题求解完成后添加对应割平面；子问题为独立PySCIPOpt Model实例，执行流程为先调用subproblem.optimize()求解，再调用subproblem.getObjVal()获取子问题目标值。
• 异常现象：部分算例运行时触发如下警告：

Warning: getSolObjVal cannot only be called in stage SOLVING without a valid solution (current stage: 10)

• 核心疑点：已显式调用subproblem.optimize()触发求解，但调用取值接口时子问题未处于预期求解阶段，返回阶段值为10。

根因说明 #

1. SCIP模型阶段值10对应SCIP_STAGE_SOLVED，即子问题已经完成求解流程、退出了SOLVING（阶段值9）的求解循环。警告文本为SCIP输出笔误，实际语义为：未传入有效解对象调用getSolObjVal接口，仅允许在SOLVING阶段执行，因此在SOLVED阶段调用无参getObjVal()会直接触发警告。
2. 该问题通常伴随两类隐含情况：一是子问题求解后未产出有效可行解（如判定为不可行、无界，或触发时间/节点限制终止时无可行解），直接调用取值接口本身逻辑不合法；二是无参getObjVal()的设计逻辑是读取求解过程中的动态解，不支持求解完成后的稳态解读取。

修复方案 #

• 第一步：子问题调用optimize()后，优先校验求解状态，禁止默认存在可行解
  求解完成后第一时间获取子问题状态，根据状态走对应分支：

  subproblem.optimize()
  solve_status = subproblem.getStatus()
  

  仅当solve_status为optimal（最优解）或bestfeasible（存在可行解但未证明最优）时，才执行目标值读取逻辑；若状态为infeasible（不可行）、unbounded（无界），直接生成对应可行性割/极射线割，不要调用取值接口；若触发求解限制且无可行解，按需求调整求解参数或终止流程。
• 第二步：读取目标值时显式传入最优解对象，替换无参getObjVal()调用
  求解完成后先获取子问题的最优可行解，基于解对象读取目标值，该调用方式在SOLVED阶段完全合法，不会触发阶段警告：

  if solve_status in ("optimal", "bestfeasible"):
      best_sol = subproblem.getBestSol()
      sub_obj = subproblem.getSolObjVal(best_sol)
      # 子问题变量取值也统一通过该解对象读取：subproblem.getSolVal(best_sol, sub_var)
  

• 第三步：子问题参数适配Benders场景（可选优化）
  若Benders子问题为线性规划（LP），可在子问题初始化时关闭MIP相关的预处理、启发式等模块，直接调用LP求解逻辑，减少不必要的阶段跳转，求解状态更可控。注意不要在主问题回调中对子问题调用freeTransform()等会重置模型阶段的接口，避免求解状态丢失。

内容的提问来源于stack exchange，提问作者optimal-br