使用 gurobipy 編寫約束

Question

我分享了 model 的參數、變量和符號：

我很難用 gurobipy 編寫方程式 7，它是 model 的約束之一。 我寫的代碼塊如下：

mdl2.addConstrs(T[i, j, k] >=  quicksum(p[l]*y[i, l, s] + s[l]*x[i, l, s] for l in N for s in ???)- d[j] - 100000*(1 - x[i, j, k])   
                for i in M
                for j in N
                for k in N) #7

你能幫我解決這個問題嗎？ 會很受歡迎的。 如果需要，我也可以分享我寫的所有 model 代碼。

Answer 1

您應該擅長s in range(k) - 總和取決於外循環中的索引k 。

為了使它更容易閱讀和理解，您可能需要像這樣切換 for 循環：

for i in M:
    for j in N:
        for k in N:
            mdl2.addConstr(
                T[i, j, k]
                >= quicksum(
                    p[l] * y[i, l, s] + s[l] * x[i, l, s] for l in N for s in range(k)
                )
                - d[j]
                - 100000 * (1 - x[i, j, k])
            )

這假設M是range(m)並且N是range(n) 。 請注意，我使用addConstr是因為在三重循環的每次迭代中只有一個約束。 在任何情況下，您都應該寫出問題的 LP 文件並檢查您正在建模的內容是否實際上符合數學公式。