The Last Mile Delivery Problem
The Last Mile Delivery Problem
X
min Cij xij (1)
(i,j)∈A
X
s.t. yi = 1 k ∈ {0, 1, ..., K},
i∈Ck
(2)
X X
xij = xji = yi ∀i ∈ V,
(i,j)∈δ + (i) (j,i)∈δ − (i)
(3)
X X
E i y i ≤ tk ≤ L i yi k ∈ {0, 1, ..., K},
i∈Ck i∈Ck
(4)
X X
tk − tℓ + Tij xij ≤ L u yu − Ev yv − (Li − Ej )xij ∀(i, j) ∈ A, i ∈ Ck , j ∈ Cℓ ,
u∈Ck v∈Cℓ
(5)
yi ∈ {0, 1} ∀i ∈ V,
(6)
xij ∈ {0, 1} ∀(i, j) ∈ A.
(7)
where δ + (i) = {(i, j) ∈ A, j 6= i} and δ − (i) = {(j, i) ∈ A, j 6= i}. The objective func-
tion (1) minimizes the overall costs. Constraints (2) ensure that exactly one vertex from
each cluster is visited. Constraints (3) are flow conservation constraints. Constraints (4)
ensure that a vertex is visited during its time-window. Constraints (5) ensure that the ar-
rival and traveling times are consistent, meanwhile eliminating subtours. Constraints (6)
and (7) are variable definitions.
When time-windows are not considered, the GTSPTW reduces to the well-known
Generalized Travelling Salesman Problem (Fischetti et al. [2]). For the multi-vehicle case,
the problem is named Generalized Vehicle Routing Problem with TW (Moccia et al. [4]).
The special case where time-windows on clients do not overlap has been recently considered
by Reyes et al. [3]. The problem is called the Vehicle Routing Problem (VRP) with
Roaming Deliveries and models the case when only trunk deliveries are considered.
2 A branch-and-cut scheme
We propose several valid inequalities for GTSPTW. Some of them are derived from valid
inequalities for the asymmetric traveling salesman problem with time-windows, others are
specific. The main constraints are the followings :
P
• Feasible path inequality. xij + h∈Sij yh ≥ yi + yj − 1 i ∈ Ck , j ∈ C ℓ ,
where Sij = {h ∈ V \ (Ck ∪ Cℓ )|Ei + Tih ≤ Lh , Eh + Thj ≤ Lj , Ei + Tih + Thj ≤ Lj }
for i ∈ Ck and j ∈ Cℓ .
Two vertices are visited either directly or there exists a connection between them.
We develop a branch-and-cut algorithm for the GTSPTW. We include at the root node
of the branch-and-bound tree all polynomial sets of inequalities while GSEC inequalities
and clique inequalities are separated in the course of the algorithm. An initial solution is
provided to the algorithm thanks to a heuristic which combines the generation of sequences
of clusters and the solution of shortest path problems with resources constraints.
The algorithm is implemented in C++ using Cplex 12.6 and the Concert technology.
Preliminary results are obtained on instances generated from GTSP instances. When
creating an instance, we guarantee that a feasible solution exists. We set the CPU time
limit to 1 hour. Instances with up to 22 clusters and 107 vertices are solved solved to
optimality. The average solution time for instances with less than 20 clusters is 19.4
seconds while it is equal to 441.6 seconds for instances with 20 and 22 clusters.
References
[1] M. Joerss, J. Schrder, F. Neuhaus, C. Klink and F. Mann, “Parcel Delivery the Future
of Last Mile”, Technical Report, McKinsey& Company, 2016.
[2] M. Fischetti, J.J. Salazar González and P. Toth, “A branch-and-cut algorithm for
the symmetric generalized traveling salesman problem”, Operations Research 45(3),
378–394, 1997.
[3] D. Reyes, M. Savelsbergh and A. Toriello, “Vehicle routing with roaming delivery
locations”, Transportation Research Part C: Emerging Technologies 80, 71–91, 2017.
[4] L. Moccia, J.F. Cordeau and G. Laporte, “An incremental tabu search heuristic for the
generalized vehicle routing problem with time windows”, Journal of the Operational
Research Society 63, 232–244, 2012.