多商品流点弧模型 python+gurobi

优化模型多商品流，python+gurobi实现，带两个例子

# -*- coding: utf-8 -*-
“““
Created on Thu Apr 25 10:26:08 2019

@author: Demons
“““

from readdata_siouxfalls import readin
from gurobipy import *
import datetime
from out_inner_nodearc import out

N_list A_list A_cost A_cap OD_list OD = readin（）


m = Model（）

y = m.addVars（OD_list A_list lb = 0 ub = GRB.INFINITY name = ‘y‘）

m.setobjective（quicksum（[y[stij]*A_cost[（ij）] for i in N_list for j in N_list for s in N_list for t in N_list if （（ij） in A_list and （st） in OD_list）]）GRB.MINIMIZE）

fin=m.addConstrs（（quicksum（[y[stsj] for j in N_list if （sj） in A_list]）-quicksum（[y[stjs] for j in N_list if （js） in A_list]） == OD[（st）] for s in N_list for t in N_list if （st） in OD_list）‘flow_o‘）

fout=m.addConstrs（（quicksum（[y[sttj] for j in N_list if （tj） in A_list]）-quicksum（[y[stjt] for j in N_list if （jt） in A_list]） == -OD[（st）] for s in N_list for t in N_list if （st） in OD_list）‘flow_d‘）

fm=m.addConstrs（（quicksum（[y[stij] for j in N_list if （ij） in A_list]）-quicksum（[y[stji] for j in N_list if （ji） in A_list]） == 0 for i in N_list for s in N_list for t in N_list if （（st） in OD_list and i != s and i != t））‘flow_mid‘）

fc=m.addConstrs（（quicksum（[y[stij] for s in N_list for t in N_list if （st） in OD]） <= A_cap[（ij）] for i in N_list for j in N_list if （ij） in A_list）‘cap‘）


m.optimize（）





if m.status == GRB.OPTIMAL:
    flow_assign arc_load arc_rito = out（OD_listN_listA_listA_capy）



                

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件         91  2019-04-25 14:55  basic_nodearc\arc_five.csv

     文件        336  2019-04-25 09:59  basic_nodearc\arc_siouxfalls.csv

     文件      20598  2019-04-25 21:47  basic_nodearc\five.lp

     文件      51749  2019-04-26 16:31  basic_nodearc\gurobi.log

     文件       1523  2019-04-26 16:31  basic_nodearc\inner_nodearc.py

     文件        942  2019-04-26 10:49  basic_nodearc\inner_nodearc_Dual.py

     文件         15  2019-04-25 14:52  basic_nodearc\nodes_five.csv

     文件         85  2019-04-25 09:43  basic_nodearc\nodes_siouxfalls.csv

     文件        929  2019-04-25 15:55  basic_nodearc\out_inner_nodearc.py

     文件       1510  2019-04-26 14:23  basic_nodearc\readdata_five.py

     文件       1463  2019-04-26 14:23  basic_nodearc\readdata_siouxfalls.py

     文件        140  2019-04-25 15:01  basic_nodearc\st_OD_five.csv

     文件       4562  2019-04-25 10:21  basic_nodearc\st_OD_siouxfalls.csv

     文件      19855  2019-04-25 14:41  basic_nodearc\st_OD_siouxfalls.xlsx

     文件       1003  2019-04-26 14:23  basic_nodearc\__pycache__\out_inner_nodearc.cpython-36.pyc

     文件       1319  2019-04-26 14:23  basic_nodearc\__pycache__\readdata_siouxfalls.cpython-36.pyc

     目录          0  2019-04-26 14:23  basic_nodearc\__pycache__

     目录          0  2019-04-26 14:23  basic_nodearc

----------- ---------  ---------- -----  ----

               106120                    18