车辆段位置（配电2）（mpex20）
/***************************************************************//* *//*S A S S A M P L E L I B R A R Y*//* *//*名称：mpex20*//*标题：车辆段位置（配电2）（mpex20）*//*产品：或*//*系统：所有*//*程序：OPTMODEL*//*数据：*//* *//*支持：更新：*//*参考：*//*MISC：数学编程示例20*//*示例手册*//* *//***************************************************************/数据arc_data；输入j$12。我12美元。成本；数据线；纽卡斯尔利物浦0.5伯明翰利物浦0.5伯明翰-布莱顿0.3伦敦利物浦1.0伦敦布赖顿0.5埃克塞特利物浦0.2埃克塞特·布莱顿0.2C1利物浦1.0C1布莱顿2.0C1伯明翰1.0C2纽卡斯尔1.5C2伯明翰0.5C2伦敦1.5C3利物浦1.5C3纽卡斯尔0.5C3伯明翰0.5C3伦敦2.0C3埃克塞特0.2C4利物浦2.0C4纽卡斯尔1.5C4伯明翰1.0C4埃克塞特1.5C5伯明翰0.5C5伦敦0.5C5埃克塞特0.5C6利物浦1.0C6纽卡斯尔1.0C6伦敦1.5C6埃克塞特1.5布里斯托尔利物浦0.6布里斯托尔·布莱顿0.4北安普敦利物浦0.4北安普敦-布莱顿0.3C1布里斯托尔1.2C2布里斯托尔0.6北安普敦C2 0.4C3布里斯托尔0.5C4北安普敦0.5C5布里斯托尔0.3C5北安普敦0.6C6布里斯托尔0.8C6北安普敦0.9;数据customer_data；输入客户$需求；数据线；15000加元C2 10000C3 40000C4 35000C5 60000元20000加元;数据因子_data；输入工厂10美元。容量；数据线；利物浦150000布莱顿200000;数据仓库data；输入仓库$12。生产成本节约；数据线；纽卡斯尔70000 0 10000伯明翰500000。伦敦100000。埃克塞特40000 0 5000布里斯托尔30000120000北安普敦25000 4000 0;数据expand_depot_data；输入仓库$12。吞吐量成本；数据线；伯明翰200003000;%设max_num_depots=4；过程操作模型；设置<str，str>ARCS；num成本{ARCS}；将数据arc_data读入ARCS=[ij]cost；设置<str>FACTORIES；num容量｛FACTORIES｝；将数据factory_data读入FACTORIES=[factory]容量；设置<str>DEPOTS；num吞吐量{DEPOTS}；num open_cost{DEPOTS}；num close_savings{DEPOTS}；将数据depot_data读入DEPOTS=[仓库]吞吐量open_cost=成本关闭节省=节省；设置EXPAND_DEPOTS；num expand_throughput{expand_DEPOTS}；num expand_cost{expand_DEPOTS}；将数据expand_depot_data读取到expand_DEPOTS=[仓库]expand_throughput=吞吐量expand_cost=成本；设置<str>客户；num需求{客户}；将数据customer_data读入CUSTOMERS=[客户]需求；set NODES=FACTORIES联合DEPOTS联合CUSTOMERS；num电源{NODES}初始化0；对于工厂}中的{i，供应[i]=容量[i]；对于客户}中的{i，供应[i]=-需求[i]；var流{ARCS}>=0；con Flow_balance_con{i位于节点}中：ARCS}流[i，j]中的求和{<（i），j><=电源[i]；var IsOpen{DEPOTS}二进制；var展开{Expand_DEPOTS}二进制文件；con最大值_节点_ con：DEPOTS中的sum{i}IsOpen[i]<=&max_num_DEPOTS；con Depot_con{i in DEPOTS}：ARCS}流[i，j]中的和{<（i），j><=吞吐量[i]*IsOpen[i]+（如果i在EXPAND_DEPOTS中，则EXPAND_throughput[i]*EXPAND[i]）；Expand_DEPOTS}中的con Expand_con{i：展开[i]<=IsOpen[i]；对于DEPOTS中的{i:close_savings[i]=.}do；close_savings[i]=0；修复IsOpen[i]=1；结束；impvar固定成本=总和{DEPOTS中的仓库}（open_cost[仓库]*IsOpen[仓库]-关闭_保存[仓库]*（1-IsOpen[仓库]）+sum{depot in EXPAND_DEPOTS}EXPAND_cost[depot]*EXPAND[depote]；impvar VariableCost=ARCS}成本[i，j]*流量[i，j]中的总和{i，j>；最小总成本=固定成本+可变成本；解决；打印固定成本可变成本总成本；在ARCS中打印{i，j>：Flow[i，j].sol>0}Flow；打印IsOpen展开；退出；