CPLEX CP模块中高效实现手术唯一执行逻辑约束的方法

你当前的约束逻辑虽然正确，但生成的约束数量过于庞大（7层循环带来指数级的约束数），会严重拖慢求解器的运行效率。针对“每个手术仅能执行一次”的需求，更高效的实现方式是使用求和约束，直接限制每个手术对应的所有可能执行组合的布尔变量之和为1。

优化后的约束代码 #

forall(o in Operations) {
  sum(s in Surgeons, r in Rooms, t in Times) A[o,s,r,t] == 1;
}

为什么这个方法更高效？ #

• 原来的写法会生成 |Operations| × |Surgeons|² × |Rooms|² × |Times|² 个约束，当集合规模稍大时，约束数会爆炸式增长；
• 优化后的写法仅生成 |Operations| 个约束，数量级大幅降低，求解器处理起来更轻松。

逻辑等价性说明 #

因为 A[o,s,r,t] 是0/1布尔变量，求和结果等于1意味着：

• 恰好有一组 (s,r,t) 使得 A[o,s,r,t] = 1（手术o在该组合下执行）；
• 所有其他 (s',r',t') 组合对应的 A[o,s',r',t'] 都为0（手术o不会在其他任何条件下执行）。
  这完全等价于你原来想要实现的约束逻辑。

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