在优化模型中用Python创建四维数组

我正在努力在

python

中使用四维数组（然后在Gurobi中进行优化）

我的目标是：

model.setObjective((quicksum(r[i,j,k,l]*x[i,j,k,l] for i,j,k,l in XXX

受到若干限制：

使用以下变量：

# Define variables
x = {}
for i in range(N):
for j in range(M):
for k in range(O):
for l in range(P):
x[i,j,k,l]

我需要为我的r[I，j，k，l]输入一个4维数组，它看起来像这样：

SUPPLIER[i], CUSTOMER[j], COLOR[k], TRANSPORT[l],  r[i,j,k,l]

[1,           1,           1,         1,             12]
[1,           1,           1,         2,             34]
[1,           1,           2,         1,             34]
[1,           1,           2,         2,             23]
[1,           2,           1,         1,             67]
[1,           2,           1,         2,             34]
[1,           2,           2,         1,             87]
[1,           1,           2,         2,             53]
[2,           1,           1,         1,             12]
[2,           1,           1,         2,             34]
[2,           1,           2,         1,             23]
[2,           1,           2,         2,             12]
[2,           2,           1,         1,             24]
[2,           2,           1,         2,             78]
[2,           2,           2,         1,             09]
[2,           1,           2,         2,             12]

我已经用2D和3D做过了，但是4D阵列缺少一些东西

你有没有类似的例子可以分享，让我学习？ 你能帮我做一下语法吗？我想我需要使用numpy，但仍然不确定

---

谢谢

我想你只是遇到了范围的问题，因为范围（N）中的I意味着I=0,1,2…N-1。你写i范围（1，N+1），也写j，k，l