解决Python Gurobi约束求和中的变量作用域问题

你的代码报错核心原因是：n仅作为左侧quicksum的内部循环变量存在，作用域被限制在该推导式内，右侧表达式无法访问到这个n。根据你可能的约束逻辑，提供两种修复方案：

方案1：为每个i和n添加独立约束 #

如果你的实际需求是对每个i∈I和n∈N，让x[i,n]等于所有m≠n时(1-y[i,m])*z[i,m]的和，需要将n加入外层循环，使其成为约束的遍历变量：

model.addConstrs(
    (x[i, n] == grb.quicksum((1 - y[i, m]) * z[i, m] for m in N if m != n)
     for i in I for n in N), name="element_balance"
)

方案2：调整右侧逻辑以消除对n的依赖 #

如果你的需求确实是左侧对所有n∈N的x[i,n]求和，那么右侧的n属于逻辑错误，需要修改右侧表达式。例如，若右侧是对所有m∈N的(1-y[i,m])*z[i,m]求和：

model.addConstrs(
    (grb.quicksum(x[i, n] for n in N) == grb.quicksum((1 - y[i, m]) * z[i, m] for m in N)
     for i in I), name="total_balance"
)

关键说明 #

• 方案1对应“逐元素约束”，适用于需要为每个n单独建立等式的场景；
• 方案2对应“全局求和约束”，适用于左右两侧都是全局求和的场景。

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