6×6矩阵CPLEX模型正常，14×14矩阵报错的原因及解决方法

问题背景 #

需构建以二进制矩阵Z为决策变量的模型，要求Z为置换矩阵：

• 约束18：每列元素之和等于1
• 约束19：每行元素之和等于1
  6×6矩阵场景下模型运行正常，但14×14矩阵时求解器报错，疑似列约束失效。

问题代码 #

Z = np.empty((Total_T, Total_T), dtype = object)
for k in range(Total_T):
    for i in range(Total_T):
     Z[k][i] = mdl.binary_var(name='Z' + str(k+1) + str(',') + str(i+1))

#Constraint 18
#mdl.add_constraints(mdl.sum(Z[k][i] for k in range(Total_T)) == 1 for i in range(Total_T))

# Constraint 19
mdl.add_constraints(mdl.sum(Z[k][i] for i in range(Total_T)) == 1 for k in range(Total_T))

问题根源 #

1. 核心约束缺失：约束18（列和为1）被注释掉，模型仅添加了行和为1的约束，未满足置换矩阵的完整要求。
   • 6×6场景下，求解器可能碰巧找到同时满足行和、列和为1的解，因此未报错；
   • 14×14场景下，仅行和约束的解空间与模型隐含需求（或目标函数）冲突，导致求解器报错。
2. 变量命名使用逗号分隔，虽不会直接报错，但可能增加调试难度，存在潜在解析风险。

解决方法 #

1. 恢复列约束：取消约束18的注释，同时添加行和、列和约束，构建标准置换矩阵模型：

   # Constraint 18：每列元素和为1
   mdl.add_constraints(mdl.sum(Z[k][i] for k in range(Total_T)) == 1 for i in range(Total_T))
   
   # Constraint 19：每行元素和为1
   mdl.add_constraints(mdl.sum(Z[k][i] for i in range(Total_T)) == 1 for k in range(Total_T))
   

2. 优化变量命名：改用下划线分隔索引，提升可读性，避免特殊字符潜在问题：

   Z[k][i] = mdl.binary_var(name=f"Z_{k+1}_{i+1}")
   

3. 验证模型可行性：若恢复约束后仍报错，需检查目标函数或其他隐含约束是否与置换矩阵要求冲突，确保模型整体逻辑自洽。

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