0

これが重複していないことを願っていますが、この問題に対する具体的な答えが見つからないようです。私はPythonにかなり慣れていないので、明らかかもしれませんが、エラーを見つけることができないようです。

問題は次のとおりです。ネットワークフローモデルを最適化するために gurobi を使用しています。しかし、2 つの変数を反復処理する必要があるため、制約の作成に問題があります (それが問題だと思います)。

最初のコードは次のとおりです。

from gurobipy import *
import optimization

def beer(P, S, W, SP, g, d, b, c, m, n):

    # Create new model
    m = Model("Beer-Flow")

    # Create variables
    x = { (i,j,p) : m.addVar(name = "x[%s,%s,%s]" % (i,j,p)) for i,j in W for p in P }
    y = { (i,j) : m.addVar(name = "y[%s,%s]" % (i,j)) for i,j in W }

    # Integrate variables
    m.update()

    # Set objective
    m.setObjective(quicksum(y[i,j] * d[i,j] * c for (i,j) in W ), GRB.MINIMIZE)

    # Add constraints
    for i in S:
        m.addConstr((quicksum(x[i2,j,p])for p in P for (i2,j) in W if i2 == i ) - (quicksum(x[i2,j,p]for p in P for (i2,j) in W if j == i)) == b[i,p])

    for i,j in W:
        m.addConstr(m * y[i,j] >= (quicksum(x[i,j,p] * w[p] for p in P)))
        m.addConstr(n * y[i,j] >= (quicksum(x[i,j,p] for p in P)))

    return m

# Fill sets and parameters
P, g = multidict({"full" : 12, "empty" : 4 })
S = {"Brauerei","München","Stuttgart","Göttingen","Bielefeld","Magdeburg" }
W, d = multidict({("Brauerei", "München") : 558, ("Brauerei", "Stuttgart") : 437, ("Brauerei","Göttingen") : 142,
                  ("Brauerei", "Bielefeld") : 45, ("Brauerei", "Magdeburg") : 290, ("München", "Brauerei") : 558,                           
                  ("München", "Stuttgart") : 232, ("München", "Göttingen") : 500, ("München", "Bielefeld") : 600,
                  ("München", "Magdeburg") : 523, ("Stuttgart", "Brauerei") : 437, ("Stuttgart", "München") : 232,
                  ("Stuttgart", "Göttingen") : 408, ("Stuttgart", "Bielefeld") : 492, ("Stuttgart", "Magdeburg") : 572,
                  ("Göttingen", "Brauerei") : 142, ("Göttingen", "München") : 500, ("Göttingen", "Stuttgart") : 408,
                  ("Göttingen", "Bielefeld") : 180, ("Göttingen", "Magdeburg") : 197, ("Bielefeld", "Brauerei") : 45,
                  ("Bielefeld", "München") : 600, ("Bielefeld", "Stuttgart") : 492, ("Bielefeld", "Göttingen") : 180,
                  ("Bielefeld", "Magdeburg") : 254, ("Magdeburg", "Brauerei") : 290, ("Magedeburg", "München") : 523,
                  ("Magdeburg", "Stuttgart") : 572, ("Magdeburg", "Göttingen") : 197, ("Magdeburg", "Bielefeld") : 254})      
SP, b = multidict({("München", "full") : 1840, ("Brauerei", "empty") : -1700, ("Stuttgart", "full") : 1400, ("Stuttgart", "empty") : -1550,
                   ("Göttingen","full") : 380, ("Göttingen", "empty") : - 400, ("Bielefeld", "full") : 840, ("Bielefeld", "empty") : -800,
                   ("Magdeburg", "full") : 600, ("Magdeburg", "empty") : -500, ("Brauerei", "full") : -5600, ("Brauerei", "empty") : 4950})
c = 1
m = 24000
n = 2000


# Create model for given sets and parameters
m = beer(P, S, W, SP, g, d, b, c, m, n)

# Run optimization and print results
m.optimize()
optimization.print_result(m)

機能していない部分は次のとおりです。 

for i in S:
            m.addConstr((quicksum(x[i2,j,p])for p in P for (i2,j) in W if i2 == i ) - (quicksum(x[i2,j,p]for p in P for (i2,j) in W if j == i)) == b[i,p])

エラーは次のとおりです: 'Var' 型は反復可能ではありません。たとえば、Java では、p と i の 2 つのループで実行できました (Python で試してみましたが、うまくいきませんでした)。しかし、この問題を解決する方法がわかりません。パイソンで。

前もって感謝します

4

1 に答える 1

0

間違った場所にある括弧のように見えます:

m.addConstr((quicksum(x[i2,j,p])for p in P for (i2,j) in W if i2 == i ) - (quicksum(x[i2,j,p]for p in P for (i2,j) in W if j == i)) == b[i,p])
                               ^ Closes the quicksum call on a single value, not the whole genexpr

対quicksum同じルックアップでの他の使用法:

m.addConstr((quicksum(x[i2,j,p])for p in P for (i2,j) in W if i2 == i ) - (quicksum(x[i2,j,p]for p in P for (i2,j) in W if j == i)) == b[i,p])
                                                                                             ^ Note no paren, so it's quicksum-ing the whole generator expression
于 2016-02-09T22:42:45.683 に答える