前提・実現したいこと
Hitchcock輸送問題を目的関数が2次の関数で最適化して解くプログラムを作ろうとしています。
解きたい輸送問題は,在庫がある倉庫が4つ,供給先の工場が5つあり,輸送コストが最小になるような組み合わせ求めるという問題です。
この下の図が倉庫と工場の関係を示すマトリックスになります。
このプログラムの目的関数が最小のときに得られる最適値を出したいのですが,
目的関数の定義の部分でquicksum()の部分でいくつもエラーメッセージが発生している状態で,解決方法がわからない状態です。
お詳しい方、ご教授を頂けませんでしょうか。よろしくお願いいたします。
発生している問題・エラーメッセージ
File "<ipython-input-33-e76473079579>", line 52 +(0.5*D*(quicksum(c[x,y]*Vxyi[x,y,i] for x in X for y in Y for i in q))^2, sence = gp.GRB.MINIMIZE) ^ SyntaxError: invalid syntax
###該当のソースコード
###Python3
import grobipy as gp
d= multidict({1:2, 2:7, 3:3 , 4:2, 5:4}) # 工場の需要量Demand
S= multidict({1:5, 2:3, 3:4, 4:6}) # 倉庫の部品数Supply
Y = [1,2,3,4,5]
X = [1,2,3,4]
q = [1,2,3,4,5,6,7,8,9]
Y,d = multidict({1:2, 2:7, 3:3 , 4:2, 5:4})
X,S = multidict({1:5, 2:3, 3:4, 4:6})
c = {(1,1):5, (1,2):2, (1,3):6, (1,4):3,
(2,1):1, (2,2):3, (2,3):4, (2,4):2,
(3,1):7, (3,2):6, (3,3):8, (3,4):2,
(4,1):3, (4,2):9, (4,3):1, (4,4):2,
(5,1):3, (5,2):5, (5,3):4, (5,4):4,
}
A=29 #重み関数
B=80
C=80
D=0.55
model= gp.Model(name="Hitchcock")
Vxyi={}
for i in q:
for y in Y: for x in X: Vxyi[x,y,i] = model.addVar(vtype='C', name = 'Vxyi(%s, %s, %s)' % (x ,y ,i))
model.update()
arcs = tuplelist([(x,y,i) for (x,y,i) in Vxyi])
for x in X:
model.addConstr(quicksum(Vxyi[x,y,i] for (x,y,i) in arcs.select('*',Y,q)) == S[x], name = 'Demand(%s)' % x)
for y in Y:
model.addConstr(quicksum(Vxyi[x,y,i] for (x,y,i) in arcs.select(X,'*',q))==d[y], name = 'Capacity(%s)' % y)
model.setObjective((-0.5 * A * quicksum((1-2Vxyi[x,y,i])^2 for (x,y,i) in Vxyi))
+ (0.5B*(quicksum((S[x]-quicksum(Vxyi[x,y,i] for y in Y for i in q))^2) for x in X)
+ (0.5C(quicksum(D[y]-quicksum(Vxyi[x,y,i] for x in X for i in q))^2) for y in Y)
+(0.5D(quicksum(c[x,y]*Vxyi[x,y,i] for x in X for y in Y for i in q))^2, sence = gp.GRB.MINIMIZE)
print('Gurobi Optimizerログ')
model.optimize()
print( 'Optimal value:', model.ObjVal)
print()
print('[解]')
if model.Status == gp.GRB.OPTIMAL:
Vxyi_opt=x.X val_opt = model.ObjVal print('最適値は {val_opt}'')``
補足情報(FW/ツールのバージョンなど)
Mac
Anaconda Jupyter notebook
Python3.7.0
また,モジュールとしてgurobiをインポートとしているのですが,他に最適なモジュールがあれば教えて頂きたいです。
直接的な解決案の回答でなくとも、参考になりそうなサイト、類似の事例を教えて頂けると幸いです。
お忙しいとは思いますが、よろしくお願いいたします。
情報に不足がありましたら、ご指摘お願いいたします
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2020/01/25 07:16