如何在CPLEX中基于执行后子集添加增量约束？

核心思路 #

因为inner_beams是第一次求解后得到的子集，无法在初始模型的Subject to块中直接使用约束。正确做法是：先求解基础模型得到inner_beams，再通过主脚本（Main block）动态添加对应约束，最后重新求解更新后的模型。

正确代码实现 #

1. 基础模型定义（先求解得到inner_beams）

// 定义集合、变量、目标函数等基础内容
{int} beams = ...; // 你的beams集合定义
dvar int beam_firstchannel[beams];
dvar int beam_nomusedchannel[beams];

// 基础目标函数（示例，根据你的实际需求调整）
minimize ...;

// 基础约束（除了const02之外的其他约束）
Subject to {
    // 你的其他基础约束
    ...
}

// 定义求解后生成inner_beams的逻辑（示例，根据你的实际计算方式调整）
{int} inner_beams = {i in beams | ...}; // 这里是求解后得到子集的逻辑

2. 主脚本（Main block）动态添加约束并重新求解

Main {
    // 第一步：求解基础模型，得到inner_beams
    cplex.solve();
    
    if (cplex.getStatus() == CPLEX.OPTIMAL) {
        // 第二步：基于求解得到的inner_beams动态添加约束
        for (var i in thisOplModel.inner_beams) {
            // 创建新的约束表达式：beam_firstchannel[i] == beam_firstchannel[i-1] + beam_nomusedchannel[i-1]
            var expr = thisOplModel.beam_firstchannel[i] == 
                       thisOplModel.beam_firstchannel[i-1] + thisOplModel.beam_nomusedchannel[i-1];
            // 将约束添加到模型中
            thisOplModel.add(expr, "const02_" + i);
        }
        
        // 第三步：求解添加约束后的模型
        cplex.solve();
        
        // 执行后处理
        thisOplModel.postProcess();
    } else {
        writeln("基础模型求解失败，状态：", cplex.getStatus());
    }
}

关键修正点 #

• 移除初始Subject to中的const02约束：因为初始模型求解时inner_beams还未生成，无法生效。
• 动态创建约束：通过thisOplModel.add()方法，在得到inner_beams后逐个添加对应的约束，避免索引错误。
• 分两次求解：先得到inner_beams，再添加约束重新求解，确保约束依赖的子集已存在。

常见问题说明 #

• 如果inner_beams是通过其他逻辑（比如后处理计算得到，而非模型自带的集合推导），需要确保在第一次求解后正确赋值给thisOplModel.inner_beams，再执行约束添加步骤。
• 约束命名用const02_加索引，避免重复命名冲突。

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