求卡车-无人机车辆路径规划问题（Truck-drone VRP）的CPLEX OPL/Python代码

问题说明 #

• 目标：求解**卡车-无人机车辆路径规划问题（Truck-drone vehicle routing problem）**的CPLEX OPL或Python实现代码
• 数学模型：参考指定模型
• 问题参数：
  • 配送中心（depot）数量：1
  • 卡车数量：10辆
  • 卡车容量：20
  • 客户数量：50个
  • 单个客户需求：30

本人已尝试的OPL代码 #

{string} k = ...;  // 车辆集合
{string} i = ...;  // 客户节点集合
{string} j = ...; // 配送中心集合

int demand[i]= ...; // 节点i的需求量
int cap[k]= ...; // 车辆k的容量
int dist[i][j]=...; // 节点i到节点j的距离

dvar float x[i][j];

minimize 
    sum(nodes in i, depot in j)
      dist[i][j];
      
subject to {
  forall (nodes in i)
      demand_must_be_met:
      sum(vehicle in k)
        x[k][i]==demand[i];
}  

.dat:

k = {k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12,k13,k14,k15};

i = {Taubenmarkt,Landstrasse,LLenaupark,LLandstrasse,Industriezeile,LentiaCity,LLeonding,LInterspar,KarlSteigerStrasse,Kornstrasse,Traunerkreuzung,PlusKaufStrasse,FMZ,Haid,Marchtrenk,EferdingPupping,Walding,Rohrbach,BadLeonfelden,LAltenbergerstrasse,Gallneukirchen,Unterweitersdorf,FMZNord,Mauthausen,Furnpark,FMZEnnspark,NaarerStrasse,Steyr,SteyrResthofstrasse,SteyrHaagerStrasse,Sierning,LHartheimeerstrasse,BadHall,Micheldorf,Vorchdorf,WelsBäckergasse,Dalistrasse,WelsOberfeldstrasse,MaxCenter,GrieskirchenSchlüsselberg,DruckereistrasseSEP,Laakirchen,Schwanenstadt,Ampflwang,RiedTumeltshamFMZ,Weberzeile,Schärding,BadIschlFMZ,BadIschlZentrum,VöcklabruchZentrum,Varena,FMZTimelkam,StGeorgen,Mondsee,CityCenter,Mattighofen
};

j = {Enns};

cap = [17,17,17,17,17,17,17,17,17,17,17,17,17,17,17];

demand = [2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2
];

dist = [23,22,22,24,12,24,19,18,22,23,23,16,20,33,41,35,68,52,28,37,42,63,10,4,3,17,22,21,21,28,11,36,63,52,43,37,38,43,55,75,64,64,90,82,81,106,111,108,77,76,80,87,110,123,117
];

内容的提问来源于stack exchange，提问作者HnM