一种O2O外卖平台智能调度与路线优化方法

摘要

本发明公开了一种O2O外卖平台智能调度与路线优化方法，包括以下步骤：对外卖订单数据进行预处理；快速生成外卖订单智能指派模型；确定配送路线，估算配送成本。本发明对外卖订单数据进行预处理时考虑了配送时间、延迟成本等因素，减少了指派、路线选取环节的计算量。本发明对外卖订单进行智能快速指派提高外卖订单的指派效率。本发明利用外卖配送需先去商家取货后去客户送货的思想快速选取配送路线。本发明利用外卖订单预处理完成对动态时间、配送员位置与容量约束及配送与延迟单位成本等参数的初始化，对待指派订单按照其与已指派订单的距离智能选择配送员，利用外卖配送必先到商家取货、后向客户配送的思想为每个配送员快速选择配送路线。

说明书

技术领域

本发明属于物流调度与优化技术，特别是一种O2O外卖平台订单智能调度与路线优化方法。

背景技术

当前有很多O2O外卖平台，如百度外卖、饿了么、美团等，而外卖调度与配送作为餐饮O2O的重要支撑环节，其优劣直接影响互联网餐饮外卖的进一步扩张与发展。然而，由于区域性限制、时效性要求高、利润点相对较低等诸多原因，外卖配送物流成本居高不下，外卖调度与配送优化是外卖平台最大的难题。

目前专门针对O2O外卖平台外卖调度与配送路线的公开文献较少，多数仍然采用传统物流配送优化模型及其变形，如设施选址-分配问题(LAP)、路径优化问题(VRP)和同时送取货VRP问题(VRPSPD)。LAP重点是确定配送中心最佳位置，不允许路线巡回访问；VRP虽允许路线巡回访问，但要求开始与结束位置都必须为同一配送中心，且二者都为先从配送中心取货，后向客户配送的过程；VRPSPD为同时送取货的车辆路径问题，所有配送货物都要从配送中心装车送至客户，所有从客户装车的货物都要送至配送中心。但这些方法应用于O2O外卖平台订单调度与配送路线选取时效果都不理想：1)由于外卖配送一般会划分合理的区域范围且距离不会太远，同时为了避免影响餐品质量与口味，一般要求出餐后30分钟送达客户，有严格的配送时间约束；2)传统物流配送优化问题都有固定的配送中心(一般中心较少)且车辆都从中心出发最终回到中心的闭环结构，而外卖配送路线优化问题，必须保证先去订单的商家取餐、后配送到客户，商家、客户不能同等对待，订单的商家、客户位置为动态需求且调度时配送员位置也不固定；3)传统物流配送优化问题调度之后转化为旅行商问题(TSP)，寻找最优路线属于NP问题，多采用启发式算法求解，算法效率普遍较低，而O2O外卖平台订单非常多且比较集中，需要快速指派配送员并给出推荐路线。这些因素共同影响，使得O2O外卖平台外卖调度与配送优化问题的处理更为复杂。

发明内容

为解决现有技术存在的上述问题，本发明要提出一种同时考虑配送时间约束、商家与客户位置依据订单需求动态变化且快速高效的O2O外卖平台智能调度与路线优化方法。

为了实现上述目的，本发明的技术方案如下：一种O2O外卖平台智能调度与路线优化方法，包括以下步骤：

A、对外卖订单数据进行预处理

A1、依据当前未配送外卖订单记录，令K为已指派尚未完成配送的外卖订单数，K’为未指派配送员的外卖订单数，则Ls_k或Ls_k’分别为第k或k’个订单的商家位置，Lc_k或Lc_k’分别为第k或k’个订单的客户位置，由于位置主要用于计算距离，而外卖配送的区域范围非常小，用直线距离取代地理空间距离对结果不产生影响；k＝1,2,…,K，k’＝1,2,…,K’；

A2、令该配送区域内配送员数为N，x_ki＝1或x_k’i＝1分别表示订单k或k’指派给配送员i，x_ki＝0或x_k’i＝0分别表示订单k或k’未指派给配送员i；x_i(s,t)表示配送员i是否需要走边(s,t)，即从s位置直接沿着边(s,t)去t位置，其中是为1，否为0，边长记为l(s,t)；订单k生效时刻为T_k0,订单k的客户希望最迟送达时刻为LT_k，配送员到订单k的商家取餐的时刻为T_ks，配送员把订单k送达客户的时刻为T_kt；s,t∈{(K+K’)个订单的商家或客户位置}，i＝1,2,…,N；

A3、给定配送员平均配送速度为v，单位时间配送成本为uc，超出指定配送时间的延迟单位时间成本为dc，通过当前配送员的手机GPS获取位置数据。计算当前所有K’+K个外卖订单的所有商家Ls_k、客户位置Lc_k的边长，边长矩阵L＝(l(s,t))_2K×2K，边长实际为两点s(s₁,s₂)、t(t₁,t₂)的欧氏距离，如式(1)所示：

B、快速生成外卖订单智能指派模型

B1、给定配送员已指派未配送的订单容量约束为Q；

B2、判断：如果K’+K≤N×Q，则全部指派，否则只取K’中先下订单的(N×Q-K)作为本次要指派的订单，其他(K’+K-N×Q)个订单下次调度时指派，同时将本次需要指派的订单作为K’。

B3、将当前所有未指派订单按照订单生成的先后顺序排序，订单序号i越小，表示订单生成越早，应该更优先被指派；设定当前未指派订单序号k’＝1；

B4、判断：当前是否有空闲配送员i，i＝1,2,…,N，有则将订单k’指派给配送员i，即x_k’i＝1，转到步骤B7，否则转到步骤B5；

B5、计算当前订单k’到所有已指派未完成订单的距离，如式(2)所示：

d(k’,k)＝l(Ls_k’,Ls_k)+l(Lc_k’,Lc_k),k＝1,2,…,K>

B6、将d(k’,k)按照距离由小到大的顺序排列，对应排序后的订单序号k_min＝1,2,…,K，之后为订单k’指派配送员：

B61、设定当前排序后的订单序号k_min＝1；

B62、订单k_min已指派的配送员为i_min，判断：如果配送员i_min已指派未完成数：

则转到步骤B63；否则转到步骤B64。

B63、令k_min＝k_min+1，判断：如果k_min≤K，转到步骤B62；否则转到步骤B8。

B64、将k’指定给配送员i_min，即x_k’imin＝1。

B7、令当前商品节点序号k’＝k’+1，K＝K+1；判断：如果k’≤K’，则转到步骤B4；否则转到步骤B8；

B8、所有订单指派完毕，停止指派。

C、确定配送路线，估算配送成本

C1、设定当前配送员序号i＝1，已指派未完成订单的配送成本Z＝0；

C2、判断：如果配送员序号i≤N，则转到步骤C3，否则转到步骤C5；

C3、判断：如果配送员i没有已指派尚未完成的订单，则令Z_i＝0，转到步骤C4，否则估算配送员i的配送成本Z_i：

C31、假定指派给配送员i的订单为k_i＝1,2,…,K_i，令D_i＝{配送员i要配送订单的所有商家位置}＝{Ls_ki|k_i＝1,2,…,K_i}，D_i’＝Φ，Φ为空集，位置变量l_i为配送员i的当前位置，Z_i＝0，当前时刻为T；

C32、计算D_i中距离l_i最近的位置d_min，令INT()表示上取整，T＝T+t_i，l_i＝d_min；

C33、判断：如果d_min为某订单k_i的商家位置Ls_ki，则令同时用订单k_i对应的客户位置Lc_ki取代D_i中的d_min，T_kis＝T，Z_i＝Z_i+t_i×uc；否则d_min为某订单的客户位置，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}，T_kit＝T，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc；

C34、判断：如果D_i≠Φ(空集)，则转到步骤C32；否则D_i’中元素的顺序即为配送顺序，转到步骤C4；

C4、令Z＝Z+Z_i，i＝i+1，转到步骤C2；

C5、停止计算，输出每个配送员i的最优配送路线D_i’和总的配送成本Z。

与现有技术相比，本发明具有以下有益效果：

1、本发明对外卖订单数据进行预处理时考虑了配送时间、延迟成本等因素

外卖订单数据包含商家、客户位置、订单时刻、订单期望配送时间要求等多指标因素，为了完成对配送员的智能快速指派和配送路线快速生成，本发明预处理过程中综合考虑配送员的当前位置、配送员到达商家及客户的时刻，预估配送员的平均配送速度、单位时间配送成本、超出指定配送时间的延迟单位时间成本等因素，并计算获得所有商家、客户的距离，减少指派、路线选取环节的计算量。

2、本发明对外卖订单进行智能快速指派提高外卖订单的指派效率

为避免出现个别配送员订单特别多而其他配送员空闲的情况，本发明为配送员设置订单容量约束，同时订单指派时优先指派空闲配送员；并依据订单生成时间的顺序指派配送员，避免出现个别订单由于与其他订单距离较远而长时间未被配送的情况。依据单一商家或客户位置进行指派，易于出现先到达客户但尚未到商家取货的情况，鉴于此，本发明为订单指派配送员时，同时考虑商家和客户的位置，只有二者距离都很近的订单才会指派给同一个配送员。

3、本发明利用外卖配送需先去商家取货后去客户送货的思想快速选取配送路线

传统路线优化问题应用于外卖配送时会受以下三个方面的制约：1)多数路线优化求解算法都需要给定配送中心，车辆从配送中心出发最终回到配送中心的闭环配送路线，而外卖配送的商家位置和数量都是动态需求，且配送路线不是闭环结构；2)多数路线优化求解算法由于是NP问题，都采用启发式算法求解，耗费时间长，而O2O外卖平台的订单数据量非常大，要求快速给出建议路线，无法完成多次访问订单数据的任务；3)传统算法即使能实现同时送取货的任务，但所有送的货都来自于固定配送中心，所有取得货都要送到固定配送中心，而外卖配送所有的外卖都需要从随机、动态的商家取货、配送给动态、随机的客户。因此，本发明利用先到商家、后到客户的思想选取路线可以消除外卖配送的条件制约，直接依据配送员到待取餐商家和已取餐客户的距离进行判断，无须固定配送中心或商家位置、无须多次访问数据、无须建立复杂模型即可快速完成配送路线选取过程。

4、本发明的方法包括三个过程：一个外卖订单数据预处理过程、一个快速生成外卖订单智能指派过程和依据已指派订单确定配送路线的过程；利用外卖订单预处理完成对动态时间、配送员位置与容量约束及配送与延迟单位成本等参数的初始化，对待指派订单按照其与已指派订单的距离智能选择配送员，利用外卖配送必先到商家取货、后向客户配送的思想为每个配送员快速选择配送路线，并估算配送成本，为调度员提供调度依据。

附图说明

本发明共有附图3张，其中：

图1是O2O外卖平台智能调度与路线优化流程图。

图2是订单初始关系结构图。

图3是订单最终配送路线图。

具体实施方式

下面参考附图，详细描述本发明提出的O2O外卖平台智能调度与路线优化方法。

如图1所示，本发明包括三个阶段：

1、对外卖订单数据进行预处理；

2、快速生成外卖订单智能指派模型；

3、确定配送路线，估算配送成本。

下面以具体例子来说明本发明的实施过程。假设现有待指派订单5个，配送员3个，已指派未完成订单3个，各订单通过GPS获取商家和客户的位置，详见表1，表1中订单的初始结构如图2所示。

表1订单情况

具体步骤如下：

A、对外卖订单数据进行预处理

A1、依据当前未配送外卖订单记录，令K＝3为已指派尚未完成配送的外卖订单数，K’＝5为未指派配送员的外卖订单数，则Ls_k(或Ls_k’)为第k(或k’)个订单的商家位置，Lc_k(或Lc_k’)为第k(或k’)个订单的客户位置，由于位置主要用于计算距离，而外卖配送的区域范围非常小，用直线距离取代地理空间距离对结果不产生影响；k＝1,2,…,K，k’＝1,2,…,K’；

A2、令该配送区域内配送员数为N＝3，x_ki＝1(或x_k’i＝1)表示订单k(或k’)指派给配送员i，为0表示未指派给配送员i；x_i(s,t)表示配送员i是否需要走边(s,t)(从s位置直接沿着边(s,t)去t位置)，其中是为1，否为0，边长记为l(s,t)；订单k生效时刻为T_k0,订单k的客户希望最迟送达时刻为LT_k，配送员到订单k的商家取餐的时刻为T_ks，配送员把订单k送达客户的时刻为T_kt；s,t∈{(K+K’)个订单的商家或客户位置}，i＝1,2,…,N；

A3、给定配送员平均配送速度为v＝30km/h(0.5km/m)，单位时间配送成本为uc＝0.2元/m，超出指定配送时间的延迟单位时间成本为dc＝0.2元/m，通过当前配送员的手机GPS获取位置数据。计算当前所有K’+K(即11)个外卖订单的所有商家Ls_k、客户位置Lc_k的边长，边长矩阵L＝(l(s,t))_2K×2K，边长实际为两点s(s₁,s₂)、t(t₁,t₂)的欧氏距离，如式(1)：

由定义可知l(s,t)＝l(t,s)，边长矩阵L为对称矩阵，可表示为上三角或下三角矩阵，具体见表2。

B、快速生成外卖订单智能指派模型

B1、给定配送员已指派未配送的订单容量约束为Q＝5；

B2、判断：K’+K＝8<N×Q＝3×5，则可以全部指派。

B3、将当前所有未指派订单按照订单生成的先后顺序排序，订单序号i越小，表示订单生成越早，应该更优先被指派，假定当前订单已经按照生成顺序排列；设定当前未指派订单序号k’＝1(即订单1’)；

k’＝1

B4、判断：当前是否有空闲配送员i(i＝1,2,…,N)，当前有空闲配送员3，则将订单k’指派给配送员3，即x_k’i＝x_1’3＝1，转到B7；

B7、当前商品节点序号k’＝k’+1＝2，K＝K+1＝4，K＝4对应订单为订单1’，

当前各订单状态下表所示；判断：k’＝2≤K’＝5，则转到步骤B4；

订单1231’->42’3’4’5’订单状态配送员1配送员2配送员2配送员3未指派未指派未指派未指派

k’＝2

B4、判断：当前是否有空闲配送员i(i＝1,2,…,N)，当前没有空闲配送员，则转到B5；

B5、计算当前订单k’到所有已指派未完成订单的距离，如式(2)

d(k’,k)＝l(Ls_k’,Ls_k)+l(Lc_k’,Lc_k),k＝1,2,…,K＝4>

如d(k’,1)＝l(Ls_k’,Ls₁)+l(Lc_k’,Lc₁)＝1.68+8.91＝10.59；

同样，得到d(k’,2)＝11.92，d(k’,3)＝14.23，d(k’,4)＝5.54

B6、将d(k’,k)按照距离由小到大的顺序排列，分别为：d(k’,4)＝5.54≤d(k’,1)＝10.59≤d(k’,2)＝11.92≤d(k’,3)＝14.23，对应排序后的订单序号k_min＝1,2,…,K，k_min＝1对应的订单为订单4，k_min＝2对应的订单为订单1，k_min＝3对应的订单为订单2，k_min＝4对应的订单为订单3，之后为订单k’指派配送员：

B61：设定当前排序后的订单序号k_min＝1(对应订单4)；

B62：订单k_min已指派的配送员为i_min(配送员3)，判断：如果配送员i_min已指派未完成数转到步骤B64。

B64：将k’指定给配送员i_min＝3，即x_k’3＝1(订单2’指派给了配送员3)。

B7、当前商品节点序号k’＝k’+1＝3，K＝K+1＝5，当前各订单状态下表所示；

判断：如果k’＝3≤K’＝5，则转到步骤B4；

k’＝3

B4、判断：当前是否有空闲配送员i(i＝1,2,…,N)，当前没有空闲配送员，则转到B5；

B5、计算当前订单k’到所有已指派未完成订单的距离，如式(2)

d(k’,k)＝l(Ls_k’,Ls_k)+l(Lc_k’,Lc_k),k＝1,2,…,K＝5(2)

得到d(k’,1)＝17.75，d(k’,2)＝7.13，d(k’,3)＝15.53，d(k’,4)＝13.87，d(k’,5)＝11.3

B6、将d(k’,k)按照距离由小到大的顺序排列，分别为：d(k’,2)＝7.13≤d(k’,5)＝11.3≤d(k’,4)＝13.87≤d(k’,3)＝15.53≤d(k’,1)＝17.75，对应排序后的订单序号k_min＝1,2,…,K，k_min＝1对应的订单为订单2，k_min＝2对应的订单为订单5，k_min＝3对应的订单为订单4，k_min＝4对应的订单为订单3，k_min＝5对应的订单为订单1，之后为订单k’指派配送员：

B61：设定当前排序后的订单序号k_min＝1(对应订单2)；

B62：订单k_min已指派的配送员为i_min(配送员2)，判断：如果配送员i_min已指派未完成数转到B64。

B64：将k’指定给配送员i_min＝2，即x_k’2＝1(订单3’指派给了配送员2)。

B7、当前商品节点序号k’＝k’+1＝4，K＝K+1＝6，当前各订单状态下表所示；

判断：如果k’＝4≤K’＝5，则转到步骤B4；

订单1231’->42’->53’->64’5’订单状态配送员1配送员2配送员2配送员3配送员3配送员2未指派未指派

k’＝4

B4、判断：当前是否有空闲配送员i(i＝1,2,…,N)，当前没有空闲配送员，则转到B5；

B5、计算当前订单k’到所有已指派未完成订单的距离，如式(2)

d(k’,k)＝l(Ls_k’,Ls_k)+l(Lc_k’,Lc_k),k＝1,2,…,K＝6>

得到d(k’,1)＝11.95,d(k’,2)＝9.85,d(k’,3)＝7.5,d(k’,4)＝5.75,d(k’,5)＝9.21,d(k’,6)＝10.75

B6、将d(k’,k)按照距离由小到大的顺序排列，分别为：d(k’,4)＝5.75≤d(k’,3)＝7.5≤d(k’,5)＝9.21≤d(k’,2)＝9.85≤d(k’,6)＝10.75≤d(k’,1)＝11.95，对应排序后的订单序号k_min＝1,2,…,K，k_min＝1对应的订单为订单4，k_min＝2对应的订单为订单3，k_min＝3对应的订单为订单5，k_min＝4对应的订单为订单2，k_min＝5对应的订单为订单6，k_min＝6对应的订单为订单1，之后为订单k’指派配送员：

B61：设定当前排序后的订单序号k_min＝1(对应订单4)；

B62：订单k_min已指派的配送员为i_min(配送员3)，判断：如果配送员i_min已指派未完成数转到步骤B64。

B64：将k’指定给配送员i_min＝3，即x_k’3＝1(订单4’指派给了配送员3)。

B7、当前商品节点序号k’＝k’+1＝5，K＝K+1＝7，当前各订单状态下表所示；

判断：如果k’＝5≤K’＝5，则转到步骤B4；

订单1231’->42’->53’->64’->75’订单状态配送员1配送员2配送员2配送员3配送员3配送员2配送员3未指派

k’＝5

B4、判断：当前是否有空闲配送员i(i＝1,2,…,N)，当前没有空闲配送员，则转到步骤B5；

B5、计算当前订单k’到所有已指派未完成订单的距离，如式(2)

d(k’,k)＝l(Ls_k’,Ls_k)+l(Lc_k’,Lc_k),k＝1,2,…,K＝7(2)

得到d(k’,1)＝9.14,d(k’,2)＝6.28,d(k’,3)＝9.8,d(k’,4)＝8.48,d(k’,5)＝7.56,d(k’,6)＝12.3,d(k’,7)＝11.03

B6、将d(k’,k)按照距离由小到大的顺序排列，分别为：d(k’,2)＝6.28≤d(k’,5)＝7.56≤d(k’,4)＝8.48≤d(k’,1)＝9.14≤d(k’,3)＝9.8≤d(k’,7)＝11.03≤d(k’,6)＝12.3，对应排序后的订单序号k_min＝1,2,…,K，k_min＝1对应的订单为订单2，k_min＝2对应的订单为订单5，k_min＝3对应的订单为订单4，k_min＝4对应的订单为订单1，k_min＝5对应的订单为订单3，k_min＝6对应的订单为订单7，k_min＝7对应的订单为订单6，之后为订单k’指派配送员：

B61：设定当前排序后的订单序号k_min＝1(对应订单2)；

B62：订单k_min已指派的配送员为i_min(配送员2)，判断：如果配送员i_min已指派未完成数转到步骤B64。

B64：将k’指定给配送员i_min＝2，即x_k’2＝1(订单5’指派给了配送员2)。

B7、当前商品节点序号k’＝k’+1＝6，K＝K+1＝8，当前各订单状态下表所示；

判断：如果k’＝6>K’＝5，则转到步骤B8；

订单1231’->42’->53’->64’->75’->8订单状态配送员1配送员2配送员2配送员3配送员3配送员2配送员3配送员2

B8、所有订单指派完毕，停止指派。

C、确定配送路线，估算配送成本Z

C1、设定当前配送员序号i＝1，已指派未完成订单的配送成本Z＝0；

i＝1

C2、判断：如果配送员序号i＝1≤N，则转到步骤C3

C3、判断：如果配送员i有已指派尚未完成的订单，则估算配送员i的配送成本Z_i：

C31、假定指派给配送员i的订单为k_i＝1,2,…,K_i，令D_i＝{配送员i要配送订单的所有商家位置}＝{Ls_ki|k_i＝1,2,…,K_i}，由于只有订单1指派给配送员1，所以D_i＝{订单1的商家位置Ls₁}，D_i’＝Φ(空集)，位置变量l_i为配送员i的当前位置，不失一般性，假定当前位置为(0,0)，Z_i＝0，当前时刻为T＝12:10；

C32、计算D_i中距离l_i最近的位置d_min，由于D_i中只有一个元素，即d_min＝Ls₁，令(分钟)，T＝T+t_i＝12:20，l_i＝d_min＝Ls₁；

C33、判断：d_min(即Ls₁)为订单1的商家位置Ls_ki，则令同时用订单1对应的客户位置Lc_ki取代D_i中的d_min，即D_i＝{Lc₁}，T_kis＝T＝12:20，Z_i＝Z_i+t_i×uc＝0+10×0.2＝2元；

C34、判断：如果D_i＝{Lc₁}≠Φ(空集)，则转到步骤C32；

C32、计算D_i中距离l_i(即Ls₁)最近的位置d_min，由于D_i中只有一个元素，即d_min＝Lc₁，令(分钟)，T＝T+t_i＝12:20+14＝12:34，l_i＝d_min＝Lc₁；

C33、判断：d_min为订单1的客户位置Lc₁，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝Φ，T_kit＝T＝12:34，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝2+14×0.2+max{12:34-12:40,0}×0.2＝4.8；

C34、判断：由于D_i＝Φ，则D_i’＝{Ls₁,Lc₁}中元素的顺序即为配送顺序，转到步骤C4；

C4、令Z＝Z+Z_i＝4.8，i＝i+1＝2，转到步骤C2；

i＝2

C2、判断：如果配送员序号i＝2≤N，则转到步骤C3；

C3、判断：如果配送员i有已指派尚未完成的订单，则估算配送员i的配送成本Z_i：

C31、假定指派给配送员i的订单为k_i＝1,2,…,K_i，令D_i＝{配送员i要配送订单的所有商家位置}＝{Ls_ki|k_i＝1,2,…,K_i}，由于订单2、3、6、8都指派给配送员2，所以D_i＝{Ls₂,Ls₃,Ls₆,Ls₈}，D_i’＝Φ，位置变量l_i为配送员i的当前位置，不失一般性，假定当前位置为(0,0)，Z_i＝0，当前时刻为T＝12:10；

C32、计算D_i中距离l_i最近的位置d_min＝Ls₈，令(分钟)，T＝T+t_i＝12:19，l_i＝d_min＝Ls₈；

C33、判断：d_min(即Ls₈)为订单8的商家位置Ls_ki，则令同时用订单8对应的客户位置Lc₈取代D_i中的d_min，即D_i＝{Ls₂,Ls₃,Ls₆,Lc₈}，T_kis＝T＝12:19，Z_i＝Z_i+t_i×uc＝0+9×0.2＝1.8元；

C34、判断：由于D_i＝{Ls₂,Ls₃,Ls₆,Lc₈}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Ls₈)最近的位置d_min，即d_min＝Lc₈，令(分钟)，T＝T+t_i＝12:19+5＝12:24，l_i＝d_min＝Lc₈；

C33、判断：d_min为订单8的客户位置Lc₈，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝{Ls₂,Ls₃,Ls₆}，T_kit＝T＝12:24，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝1.8+5×0.2+max{12:24-12:50,0}×0.2＝2.8；

C34、判断：由于D_i＝{Ls₂,Ls₃,Ls₆}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Lc₈)最近的位置d_min，即d_min＝Ls₂，令(分钟)，T＝T+t_i＝12:24+13＝12:37，l_i＝d_min＝Ls₂；

C33、判断：d_min(即Ls₂)为订单2的商家位置Ls_ki，则令同时用订单2对应的客户位置Lc₂取代D_i中的d_min，即D_i＝{Lc₂,Ls₃,Ls₆}，T_kis＝T＝12:37，Z_i＝Z_i+t_i×uc＝2.8+13×0.2＝5.4元；

C34、判断：由于D_i＝{Lc₂,Ls₃,Ls₆}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Ls₂)最近的位置d_min，即d_min＝Ls₆，令(分钟)，T＝T+t_i＝12:37+4＝12:41，l_i＝d_min＝Ls₆；

C33、判断：d_min(即Ls₆)为订单6的商家位置Ls_ki，则令Lc₈，Ls₂,Ls₆}，同时用订单6对应的客户位置Lc₆取代D_i中的d_min，即D_i＝{Lc₂,Ls₃,Lc₆}，T_kis＝T＝12:41，Z_i＝Z_i+t_i×uc＝5.4+4×0.2＝6.2元；

C34、判断：由于D_i＝{Lc₂,Ls₃,Lc₆}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Ls₆)最近的位置d_min，即d_min＝Lc₆，令(分钟)，T＝T+t_i＝12:41+11＝12:52，l_i＝d_min＝Lc₆；

C33、判断：d_min为订单6的客户位置Lc₆，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝{Lc₂,Ls₃}，T_kit＝T＝12:52，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝6.2+11×0.2+max{12:52-12:50,0}×0.2＝8.8；

C34、判断：由于D_i＝{Lc₂,Ls₃}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Lc₆)最近的位置d_min，即d_min＝Lc₂，令(分钟)，T＝T+t_i＝12:52+11＝13:03，l_i＝d_min＝Lc₂；

C33、判断：d_min为订单2的客户位置Lc₂，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝{Ls₃}，T_kit＝T＝13:03，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝8.8+11×0.2+max{13:03-12:40,0}×0.2＝15.6；

C34、判断：由于D_i＝{Ls₃}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Lc₂)最近的位置d_min，由于D_i中只有一个元素，即d_min＝Ls₃，令(分钟)，T＝T+t_i＝13:21，l_i＝d_min＝Ls₃；

C33、判断：d_min(即Ls₃)为订单3的商家位置Ls_ki，则令同时用订单3对应的客户位置Lc₃取代D_i中的d_min，即D_i＝{Lc₃}，T_kis＝T＝13:21，Z_i＝Z_i+t_i×uc＝15.6+18×0.2＝19.2元；

C34、判断：如果D_i＝{Lc₃}≠Φ(空集)，则转到步骤C32；

C32、计算D_i中距离l_i(即Ls₃)最近的位置d_min，由于D_i中只有一个元素，即d_min＝Lc₃，令(分钟)，T＝T+t_i＝13:21+15＝13:36，l_i＝d_min＝Lc₃；

C33、判断：d_min为订单3的客户位置Lc₃，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝Φ，T_kit＝T＝13:36，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝19.2+15×0.2+max{13:36-12:40,0}×0.2＝33.4；

C34、判断：由于D_i＝Φ，则D_i’＝{Ls₈,Lc₈，Ls₂,Ls₆,Lc₆,Lc₂,Ls₃,Lc₃}中元素的顺序即为配送顺序，转到步骤C4；

C4、令Z＝Z+Z_i＝4.8+33.4＝38.2，i＝i+1＝3，转到步骤C2；

i＝3

C2、判断：如果配送员序号i＝3≤N，则转到步骤C3；

C3、判断：如果配送员i有已指派尚未完成的订单，则估算配送员i的配送成本Z_i：

C31、假定指派给配送员i的订单为k_i＝1,2,…,K_i，令D_i＝{配送员i要配送订单的所有商家位置}＝{Ls_ki|k_i＝1,2,…,K_i}，由于订单4、5、7都指派给配送员3，所以D_i＝{Ls₄,Ls₅,Ls₇}，D_i’＝Φ，位置变量l_i为配送员i的当前位置，不失一般性，假定当前位置为(0,0)，Z_i＝0，当前时刻为T＝12:10；

C32、计算D_i中距离l_i最近的位置d_min＝Ls₅，令(分钟)，T＝T+t_i＝12:17，l_i＝d_min＝Ls₅；

C33、判断：d_min(即Ls₅)为订单5的商家位置Ls_ki，则令同时用订单5对应的客户位置Lc₅取代D_i中的d_min，即D_i＝{Ls₄,Lc₅,Ls₇}，T_kis＝T＝12:17，Z_i＝Z_i+t_i×uc＝0+7×0.2＝1.4元；

C34、判断：由于D_i＝{Ls₄,Lc₅,Ls₇}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Ls₅)最近的位置d_min，即d_min＝Ls₄，令(分钟)，T＝T+t_i＝12:17+6＝12:23，l_i＝d_min＝Ls₄；

C33、判断：d_min(即Ls₄)为订单4的商家位置Ls_ki，则令同时用订单4对应的客户位置Lc₄取代D_i中的d_min，即D_i＝{Lc₄,Lc₅,Ls₇}，T_kis＝T＝12:23，Z_i＝Z_i+t_i×uc＝1.4+6×0.2＝2.6元；

C34、判断：由于D_i＝{Lc₄,Lc₅,Ls₇}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Ls₄)最近的位置d_min，即d_min＝Lc₄，令(分钟)，T＝T+t_i＝12:23+11＝12:34，l_i＝d_min＝Lc₄

C33、判断：d_min为订单4的客户位置Lc₄，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝{Lc₅,Ls₇}，T_kit＝T＝12:34，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝2.6+11×0.2+max{12:34-12:40,0}×0.2＝4.8；

C34、判断：由于D_i＝{Lc₅,Ls₇}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Lc₄)最近的位置d_min，即d_min＝Lc₅，令(分钟)，T＝T+t_i＝12:34+6＝12:40，l_i＝d_min＝Lc₅；

C33、判断：d_min为订单5的客户位置Lc₅，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝{Ls₇}，T_kit＝T＝12:40，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝4.8+6×0.2+max{12:40-12:40,0}×0.2＝6；

C34、判断：由于D_i＝{Ls₇}≠Φ(空集)，则转到步骤C32；

C32、查表2获得D_i中距离l_i(即Lc₅)最近的位置d_min，由于D_i中只有一个元素，即d_min＝Ls₇，令(分钟)，T＝T+t_i＝12:40+12＝12:52，l_i＝d_min＝Ls₇；

C33、判断：d_min(即Ls₇)为订单7的商家位置Ls_ki，则令D_i’＝D_i’∪{d_min}＝{Ls₅,Ls₄,Lc₄,Lc₅,Ls₇}，同时用订单对7应的客户位置Lc₇取代D_i中的d_min，即D_i＝{Lc₇}，T_kis＝T＝12:52，Z_i＝Z_i+t_i×uc＝6+12×0.2＝8.4元；

C34、判断：如果D_i＝{Lc₇}≠Φ(空集)，则转到步骤C32；

C32、计算D_i中距离l_i(即Ls₇)最近的位置d_min，由于D_i中只有一个元素，即d_min＝Lc₇，令(分钟)，T＝T+t_i＝12:52+8＝13:00，l_i＝d_min＝Lc₇；

C33、判断：d_min为订单7的客户位置Lc₇，则同时将d_min从D_i中移除，即D_i＝D_i-{d_min}＝Φ，T_kit＝T＝13:00，Z_i＝Z_i+t_i×uc+max{T_kit-LT_k,0}×dc＝8.4+8×0.2+max{13:00-12:50,0}×0.2＝12；

C34、判断：由于D_i＝Φ，则D_i’＝{Ls₅,Ls₄,Lc₄,Lc₅,Ls₇,Lc₇}中元素的顺序即为配送顺序，转到步骤C4；

C4、令Z＝Z+Z_i＝38.2+12＝50.2，i＝i+1＝4，转到步骤C2；

C2、判断：如果配送员序号i＝4>N，则转到步骤C5；

C5、停止计算，输出每个配送员i的最优配送路线D_i’和总的配送成本Z：每个配送员i的最优配送路线为D_i’，分别为：D₁’＝{Ls₁,Lc₁}，D₂’＝{Ls₈,Lc₈，Ls₂,Ls₆,Lc₆,Lc₂,Ls₃,Lc₃}，D₃’＝{Ls₅,Ls₄,Lc₄,Lc₅,Ls₇,Lc₇}；总的配送成本Z＝50.2元，配送路线如图3所示。

上面是本发明提出的方法的一种实施方式，但在某些步骤上，可以进行适当改变，以适应具体情况的要求。例如，在第一阶段的步骤A3利用两点位置求边长，可以根据实际条件需要适当调整为空间距离。例如，在第三阶段估算配送成本时，配送员的配送路线选取策略可以适当调整。

本发明不局限于本实施例，任何在本发明披露的技术范围内的等同构思或者改变，均列为本发明的保护范围。