Next Article in Journal
A MODIS-Based Retrieval Model of Suspended Particulate Matter Concentration for the Two Largest Freshwater Lakes in China
Next Article in Special Issue
ePedigree Traceability System for the Agricultural Food Supply Chain to Ensure Consumer Health
Previous Article in Journal
The Role of Policy Makers and Institutions in the Energy Sector: The Case of Energy Infrastructure Governance in Nigeria
Previous Article in Special Issue
Regional Port Productivity in APEC
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Location-Routing Problem with Simultaneous Home Delivery and Customer’s Pickup for City Distribution of Online Shopping Purchases

1
College of Mechanical Engineering, Chongqing University, Chongqing 400030, China
2
Chongqing Key Laboratory of Logistics, Chongqing University, Chongqing 400030, China
3
State Key Laboratory of Mechanical Transmission, Chongqing University, Chongqing 400030, China
*
Author to whom correspondence should be addressed.
Sustainability 2016, 8(8), 828; https://0-doi-org.brum.beds.ac.uk/10.3390/su8080828
Submission received: 21 June 2016 / Revised: 7 August 2016 / Accepted: 18 August 2016 / Published: 22 August 2016
(This article belongs to the Special Issue Sustainability in Supply Chain Management)

Abstract

:
With the increasing interest in online shopping, the Last Mile delivery is regarded as one of the most expensive and pollutive—and yet the least efficient—stages of the e-commerce supply chain. To address this challenge, a novel location-routing problem with simultaneous home delivery and customer’s pickup is proposed. This problem aims to build a more effective Last Mile distribution system by providing two kinds of service options when delivering packages to customers. To solve this specific problem, a hybrid evolution search algorithm by combining genetic algorithm (GA) and local search (LS) is presented. In this approach, a diverse population generation algorithm along with a two-phase solution initialization heuristic is first proposed to give high quality initial population. Then, advantaged solution representation, individual evaluation, crossover and mutation operations are designed to enhance the evolution and search efficiency. Computational experiments based on a large family of instances are conducted, and the results obtained indicate the validity of the proposed model and method.

1. Introduction

The rapid development of e-commerce has generated growing interest in online shopping. Taking China as an example, according to China’s E-commerce Report 2015 [1], in 2015, the e-commerce transactions totaled 2.5 trillion US dollars and online shopping transactions accounted for 80% of the total e-commerce transactions with an increasing rate of 27.2% compared to 2014. The sharply increased online shopping market leads to an explosive growth of city distribution demand. The annual package delivery quantity reached 20 billion in 2015, leading China to rank first in the world. The Last Mile delivery works as the final stretch that delivers consignments to the recipients [2] and plays an indispensable role in promoting e-commerce development. A survey showed that 85% of buyers who received their orders on time would purchase online again compared with those 33% who received delayed orders [3]. Different to traditional city distribution, the customers engaging in online shopping are spatially distributed and have personalized delivery demands, rendering the Last Mile delivery the most expensive and pollutive, as well as the least efficient, stages of the whole e-commerce supply chain, accounting for 13% up to 75% of the total supply chain costs [4]. Thus the Last Mile delivery problem has recently garnered a great deal of attention from the government and all the stakeholders of e-commerce supply chain.
In the Last Mile delivery, the most widespread delivery mode is home delivery (HD), by which couriers send packages directly to a customer´s home or office. HD service provides a face-to-face opportunity with customers, but its low operational efficiency is likely to make it undesirable to handle mass orders, thereby increasing operation costs [5]. Subsequently, another kind of service named customer’s pickup (CP) has been widely accepted by both couriers and customers. More precisely, CP service allows customers to retrieve their packages at a facility close to their home or work (we refer to it as pickup point). By providing such an option that allows packages to be delivered to a facility with a large storage area, efficiency will be improved while maintaining personalized service. Customers can pick up their packages at their convenience, thus improving their satisfaction. As the key issue for the Last Mile delivery is to efficiently design the Last Mile delivery system [6]. To better deal with the Last Mile delivery problem we are facing, in this paper, we proposed a novel multi-depot location-routing problem with simultaneous home delivery and customer’s pickup (LRPSHC). In this problem, each customer can be serviced either by HD or CP; the purpose is to decide on the location of pickup points and customers’ service options, as well as schedule the vehicles in order to minimize the total operation costs. For solving this problem, a hybrid evolution search algorithm by combining genetic algorithm (GA) and local search (LS) is designed. This research allows us to operate the Last Mile delivery in an economical and sustainable way.
The remainder of this paper is structured as follows. The next section provides a comprehensive review of the recent research on Last Mile delivery and the variants of the location-routing problem, as well as the solving methods. Section 3 details the description and formulation of the proposed LRPSHC. In Section 4, the proposed hybrid evolution search algorithm is presented. Section 5 displays the results of the computational experiments and the discussion. Conclusions are given in Section 6. Finally, future research directions are given in Section 7.

2. Literature Review

Realizing that the Last Mile delivery is one of the most expensive stages of the entire e-logistics chain, a cost simulation tool was developed by Gevaers et al. [7] to simulate Last Mile delivery costs. According to their research, changes within Last Mile characteristics significantly influence cost. To more clearly show how certain factors affect the Last Mile delivery cost, customer density and delivery window size were considered by Boyer et al. [8]. Simulation results showed that greater customer density and longer delivery windows benefit the delivery efficiency enormously. Based on three delivery modes—HD, reception box and collection-and-delivery points—Wang et al. [5] undertook a quantitative study of the competitiveness of the three modes through analyzing the delivery cost structure and operation efficiency in different scenarios, which helps identify the most suitable mode for different customer distribution densities. Hayel et al. [9] proposed a queuing model to describe the Last Mile delivery system with HD and CP options. Considering the monetary and congestion effect of the two options, a game-theoretical approach was designed for determining the optimum option a consumer would make. Aized and Srai [4] provided a conceptual approach for modeling the Last Mile delivery system based on hierarchy to support the routing planning.
Another key issue related to our research is the location-routing problem (LRP). Several recent surveys are available which help provide a comprehensive overview of the studies about LRP, such as the work of Lopes et al. [10], and most recently that of Prodhon and Prins [11] and Drexl and Schneider [12]. Among these studies, LRP variants were demonstrated, including open LRP [13], LRP with time windows [14,15], LRP with split demands [16], LRP with subcontracting options [17], LRP with simultaneous delivery and pickup [18,19,20,21,22], LRP with depots connected by ring [23], and more complex variants with multi-period planning [24,25] and inventory management [26,27,28]. Another key variant of the LRP is two-echelon LRP [20,22,29,30,31].
When it comes to the solving approaches, both exact algorithms and heuristics were proposed to solve LRP and its variants. Interested readers are encouraged to examine the research of Lopes et al. [10], in which a review of the solving methods was provided. For the recent and relevant approaches, one can find ant colony optimization (ACO) [32], genetic algorithm(GA) [33,34], simulated annealing (SA) [13,35], simulated annealing with variable neighborhood search (SA-VNS) [18], local search (LS) [20], greedy randomized adaptive search with evolutionary local search (GRASP+ELS) [36], granular tabu search (GTS) [37], granular tabu search with variable neighborhood search (GVTNS) [38], multi-start iterated local search with tabu list and path relinking (MS-ILS) [31] and multi-start simulated annealing heuristic (MSA) [21]. The most related methods are the metaheuristics hybrid with GA and LS. Prins et al. [39] presented a memetic algorithm with population management (MAPM). The main idea behind MAPM is that LS is used to improve the solution generated by GA. A similar method was described by Derbel et al. [40]; iterative local search [41] is adopted to improve the solution generated by the GA. Marinakis and Marinaki [42] formulated a bi-level programming for LRP, and a bi-level GA was designed to solve the problem accordingly. In this method, GA is used to locate the facilities in the first level and the VRP in the second level is solved by expanding neighborhood search (ENS). Most recently, Lopes et al. [43] proposed a simple and effective hybrid GA with LS procedure that works as a mutation operator.
By reviewing previous studies, we find that scant research has been conducted on the LRP with simultaneous HD and CP, and thus comes the main contribution of this study: we propose a novel location-routing problem with simultaneous home delivery and customer’s pickup for the Last Mile delivery. Since customers have two available services to choose from—HD to be served by vehicles directly or CP to pick up their packages from the closest pickup points—delivery efficiency will be improved and customers’ personalized demands addressed.

3. Problem Formulation

3.1. Problem Description

The proposed LRPSHC can be formally stated as follows. Letting G = ( N , A , C ) represents a weighted and undirected network containing a set of vertices N, a set of edge A and a set of connection C. Where, N = N D N P N C , N D = { 1 , , n d } stands for n d depots, N P = { n d + 1 , , n d + n p } stands for n p potential pickup points and N C = { n d + n p + 1 , , n d + n p + n c } stands for n c customers should be served. Each customer i N C has a delivery demand q i and a service time t i , t i = q i u H , where u H represents the unit demand service time. Once selected, a pickup point p N P is opened with an opening cost U p containsedcapacity B p and fixed service time u P . Set A = { ( i , j ) : i , j N } is a set of edges connecting each pair of nodes in N, with each edge ( i , j ) A associated with a routing cost c i j and a travelling time t i j which are dependent on the distance d i j between nodes i and j, i , j N . Set C = { ( p , c ) : p N P , c N C } contains the connections between pickup points and customers, a connection ( p , c ) occurs only when the distance d p c between the selected pickup point p N P and customer c N C is less than the acceptable distance D p c . For a given d N D , a capacity S d is associated with it. A fleet of identical vehicles n v with capacity Q v and maximum working time T v are available and denoted as V D . Once used, a fixed cost F v incurs.
In this research, to ensure acceptability and reliability, the following two real-world situations are considered:
(1) The preference for a customer to choose CP service is based on the acceptable distance to the closest selected pickup point, and the customer will choose CP service only when within the acceptable distance D c p as mentioned before; otherwise, HD service should be provided.
(2) Due to uncertain factors that customers face in real life, HD service may fail for the first delivery with probability p H with unit second delivery cost μ s d .
An example of the Last Mile distribution system with three depots, six candidate pickup points and 17 customers are given in Figure 1.
The aim is then to determine the following:
(1)
A set of pickup point locations;
(2)
Service options for customers;
(3)
Assignments between customers and selected pickup points;
(4)
Assignments between the selected pickup points and depots;
(5)
Vehicles´ scheduling and routing.
Without losing generality, the following assumptions are made:
(1)
As online shoppers may be clustered in the same or nearby buildings, all the customers are divided into small groups according to their locations. We then take each small group as a customer with fixed demand;
(2)
Each customer can be served by either HD or CP service;
(3)
We do not consider the operation of the second delivery caused by uncertain factors in the current scheduling, and take the second delivery cost instead.

3.2. Notions and the Proposed Model

The notions used in the formulation are summarized in Table 1.
Based on the above statements, we now formulate the model for the LRPSHC as follows:
min p N P U p y p + d N D k V D F v x d j k d + d N D k V D i N j N c i j x i j k d + d N D k V D i N j N q i p H μ s d x i j k d .
subject to
j N k V D x i j k d + p N P z i p = 1 , i N C , d N D .
d N D w p d i N C z i p , p N P .
d N D w p d = y p , p N P .
j N x i j k d = j N x j i k d , i N , k V D , d N D .
j N φ i j k d j N φ j i k d = q i , i N C , k V D , d N D .
j N φ i j k d j N φ j i k d = j N C q i z j p , i N P , k V D , d N D .
φ d i k d = c N C j N x c j k d q c + p N P c N C j N x p j k d q c z c p , i N C N P , k V D , d N D .
φ i d k d = 0 , i N P N C , k V D , d N D .
φ i j k d Q v , i , j N , k V D , d N D .
i N C q i u H x i j k d + p N P u P x p j k d + i N j N t i j x i j k d T v , k V D , d N D .
i N C j N q i x i j k d + p N P i N C q i z i p w p d S d , d N D .
i N C q i z i p B p , p N P .
d N d j N P N C k V D x d j k d n v , d N D .
φ i j k d 0 , i , j N , k V D , d N D .
x i j k d = { 0 , 1 } , i , j N , k V D , d N D .
y p = { 0 , 1 } , p N P .
z i p = { 0 , 1 } , i N C , p N P .
w p d = { 0 , 1 } , p N P , d N D .
The objective function (1) minimizes the sum of the pickup point opening cost, fixed vehicle cost, routing cost and the second delivery cost. Constraints (2) ensure that each customer is served exactly once by either HD or CP service. Constraints (3) impose that once a customer is assigned to a pickup point, this pickup point should be served by a depot. Constraints (4) make sure that once a pickup facility is selected, it should be served by a depot. Constraints (5) are the route continuity constraints that once a vehicle enters a node, it must also leave from it. Constraints (6) and (7) are the flow constraints for demand of customer and pickup point nodes respectively. Constraints (8) calculate the demand delivered by vehicles. Constraints (9) specify that the remaining demand of a vehicle is zero after serving its last customer. Constraints (10) and (11) guarantee that the capacity and working time of vehicles should be satisfied. Constraints (12) and (13) imply that vehicles should not violate the capacity constraints of the depots and pickup points respectively. Constraints (14) require that selected vehicles should not violate available vehicles. Constraints (15)–(19) define all variables.

4. Hybrid Approach

As LRP is an NP-hard optimization problem [40], the multi-pickup point and simultaneous HD and CP services included in the considered problem have further increased the complexity of the solution. In this section, a hybrid evolution search algorithm (HGALS) that combines GA and local LS is presented, which is motivated by the excellent global search ability obtainable by using GA and the ability in finding local optima by taking LS.
In this algorithm, GA-based operations are intended to guide the global evolution while LS is taken to find better local solutions both for the initial population and the new generated individuals during the evolution process. LS is also used to work as a post-optimization procedure in order to improve the best solutions in each iteration. The solution improved by LS is accepted with probability 1 when the solution is feasible; otherwise, the solution with probability pac is accepted. Additionally, a hybrid approach along with a two-phase solution generation heuristic is proposed to obtain high quality initial population; specific crossover and mutation operators are also proposed followed by a biased fitness function to evaluate the individuals. A high level overview of the presented HGALS is depicted below.
Step 1: Parameter setting: population size M0, maximum population M, crossover probability, mutation probability pm, infeasible solution accept probability pac, solution improve probability pls, terminate condition
Step 2: Taking solution generation and population initialization methods to obtain initial population P0
Step 3: Improve the population by LS
Step 4: Select individual X1 and X2 from the population following probability pc, conduct crossover operation to generate offspring Y1 and Y2, apply LS to improve Y1 and Y2 with probability pls, Y1 and Y2 are improved to be Y1′ and Y2′, accept Y1′ and Y2′ by probability 1 or pac
Step 5: Select individual X3 following probability pm to make new solution Y3 by mutation operation, apply LS to improve Y3 with probability pls, Y3 is improved to be Y3′, accept Y3′ by probability 1 or pac
Step 6: Improve the current best solution Y4 to be Y4′ by LS, accept Y4′ by probability 1 or pac
Step 7: If the population reaches M, evaluation individuals and select M0 solutions to make population P0, turn to step 3; otherwise, turn to step 8
Step 8: If the terminate condition is satisfied, turn to step 9; if not, turn to step 4
Step 9: Output the best solution.
The general structure of the presented HGALS is given in Figure 2.

4.1. Solution Representation

For a GA-based heuristic, chromosome encoding strategy affects the performance of genetic operations [40]. Thus, the first and most crucial step in designing GA is to define an appropriate encoding scheme. According to the characteristics of the considered problem, a mixed encoding scheme is designed to represent a solution in which each chromosome includes two sections.
The first section is composed of three tiers of gene string with a fixed length of nc. The third tier is customer tier, which stores the customers. The second tier is used to index the service options information: if the value of a gene belongs to {nd + 1, …, nd + np}, it means that the corresponding customer in the third tier is served by CP and served by the pickup point. Otherwise, if the value of the gene is 0, it means that the customer is served by HD. The first tier is for depots and indicates the assignments between the depots and pickup points as well as depots and the customers with HD service.
The second section is to deal with the vehicle routing with a permutation of nd depots, np′ pickup points selected from NP, nc′ customers with HD service and a set of dummy zeros. Routes are started and ended at the same depots to serve the selected pickup points and customers with HD service one by one, zeros are used to terminate a route and start a new one in depots.
Figure 3 provides a solution representation of the sample given in Figure 1.

4.2. Solution and Population Initialization

As the quality of the initial solutions affects the final solution to a large extent, to obtain high quality initial solutions, we designed a two-phase heuristic named “Routing planning after the pickup points selection and services determination”. The procedure of the two-phase heuristic is depicted as follows:
Phase 1: Pickup points selection and service determination.
Step 1: Randomly select a pickup point combination set N s p N S P .
Step 2: Sort all the customers according to distance to the pickup points in N s p , let n p c be the number of the customers whose closest pickup point is p, and the corresponding customer set be N p C .
Step 3: Select the pickup point p with the largest n p c .
Step 4: If within the pickup point’s capacity, let the service option of its closest customer in N p C be CP. Delete the customer in N C and N p C . If p is in N s p , delete it from N s p and put it into N C .
Step 5: If N p C = , turn to step 2; otherwise, turn to step 3.
Step 6: If N s p , turn to step 2; otherwise, turn to step 7.
Phase 2: Route planning
Step 7: Let N d c be the node set that store the customers with HD service and selected pickup points in N C and whose closest depot is depot d. Rank all the nodes in N C with the increasing distance to the depots and put them into the corresponding sets.
Step 8: For each depot d N D and the nodes in N d c , the Traveling Salesman Problem (TSP) can be applied; the efficient algorithm proposed by Lin and Kernighan [44] can be used for generating the TSP tour.
Step 9: Split the TSP tour by the method proposed by Prins [45] without violating the vehicle capacity and working time constraints.
Step 10: Output the solution.
Based on the solution generation heuristic, the diverse population generation procedure is designed as follows:
Step 1: Randomly generate a non-redundant pickup point combination set N S P with n s p combinations, letting P 0 = .
Step 2: Conduct the solution generation algorithm to generate solution S.
Step 3: If the selected pickup point in solution S already exists in P0, turn to step 1; otherwise, put S into P0.
Step 4: If the initial population size M0 is reached, turn to step 5; otherwise, turn to step 1.
Step 5: Terminate and output the initial population.
As can be seen by taking the solution and population initialization methods, solutions are different for selected pickup points and service options, each of which stands for a unique solution space, which will cover the whole solution space to the largest extent.

4.3. Fitness Function

To compare two individual factors, the evaluation function F is defined as the sum of two components which are the objective function value and penalty cost, the evaluation function is formulated as follows:
F = F 0 + α 1 Δ S D + α 2 Δ B p + α 3 Δ n v + α 4 Δ Q v + α 5 Δ T v
where, F 0 is the objective value of a solution; Δ S D , Δ B P , Δ n v , Δ Q v and Δ T v stand for the amount that are beyond the constraints of the depot capacity, pickup point capacity, available vehicles, vehicle capacity and working time respectively; α 1 , α 2 , α 3 , α 4 and α 5 are the corresponding constant parameters that reflect the degree of the penalty.

4.4. Individual Evolution and Selection

In population-based heuristics, one of the greatest advantages is search parallel in wide solution space, rendering population diversity an important factor that affects the performance. Current widely used objective function-based individual evolution may shrink the search space and lead to premature convergence. In this study, we revise the biased fitness function proposed by Vidal et al. [46]. In this method, the objective value of a solution and its contribution to population diversity are jointly considered.
For an individual I, letting its nclose closest neighbors be stored in set Nclose, the diversity contribution of individual I to the population in Nclose is defined as the average distance to the individuals in Nclose computed by Formula (21). A normalized Hamming distance based on pickup points’ status, customers’ service options and the assignments of two individuals is adopted to show the difference between solutions. The distance is computed according to Formula (22).
Δ ( I ) = 1 n c l o s e I 1 N c l o s e D H ( I , I 1 )
D H ( I , I 1 ) = γ 1 1 n c i n c ( 1 ( s c i ( I ) s c i ( I 1 ) ) ) + γ 2 1 n p i n s ( 1 ( s p i ( I ) s p i ( I 1 ) ) ) + γ 3 1 n c i n c ( 1 ( p c i ( I ) p c i ( I 1 ) ) )
where, sci(I) stands for the assignment of customers with HD service, spi(I) stands for the assignment of pickup facilities and pci(I) stands for the assignment of customers with CP service. γ1, γ2 and γ3 are the coefficients of the three factors, γ1 + γ2 + γ3 = 1.
Let r(I) be the rank of individual I in a subpopulation with respect to F and the biased fitness function is given by Formula (23).
F E = r ( I ) Δ ( I )
It can be seen that, for each subpopulation with nclose customers, when solutions are selected by elitist strategy, excellent solutions can be protected as well as the diversity achieved.

4.5. Crossover Operation

Crossover operation aims to generate offspring by mimicking mating between parents and exchanging gene information in the hope that children will obtain the best parts of their parents, as well as diversify the search. Due to the fact that the genes and length of chromosomes vary in solutions, a modified two-point crossover is introduced, which is depicted in Figure 4. This crossover operation is conducted by exchanging the middle part of the second tier of the first section of the selected two parents between the two selected crossover positions. The process is described as follows:
Step 1: Randomly select a pair of parents P1 and P2 by probability pc from the population.
Step 2: Uniformly select two positions from [1, nc], letting the smaller one be cp1 and the other cp2.
Step 3: For each gene in position [cp1, cp2], exchange the middle part of the second tier of the first part of the two chromosomes.
Step 4: Adjust the second sections of the two chromosomes according to the mapping relationship.
By taking this crossover operation, new solutions are generated by exploring new statuses of pickup points, customers’ service options as well as the routes.

4.6. Mutation Operation

Mutation operation takes place immediately after the crossover operation, which may help GA jump out of the local optima by finding uncovered solution spaces [47]. In our research, we take gene-based mutation, which is conducted on the second section of a chromosome, in the following process:
Step 1: Select a solution P3 according to probability pm.
Step 2: Randomly select a mutation position cp3 from 1 to L(P3), where L(P3) is the length of the second section of the chromosome.
Step 3: Identify the type of the selected gene; conduct the following two cases equally.
Case 1: If it is a pickup point, each of the three operators is selected randomly:
Operator 1: Close this pickup point.
Operator 2: Exchange the assigned customers with the closest pickup point if it is open.
Operator 3: Close it and open a randomly selected one which is not in the current solution.
Case 2: If the gene is a customer or zero, remove it and insert it into a randomly selected position between 1 and L(P3).
It can be seen that, in this crossover, the operation can be performed on all search spaces by changing the status of pickup points, customers’ service options as well as the routes.

4.7. Local Search

Local search aims to improve the new generated solutions during the evolution process within a few modifications. As the neighborhood defines the way to produce a new solution, the most important issue when designing a local search is to identify effective neighborhood structure [40]. For the customers’ service options and vehicle routing in our problem, we design two kinds of neighborhood structures for service option and vehicle routing respectively.
Local search aims to improve the new generated solutions during the evolution process within a few modifications. As the neighborhood defines the way to produce a new solution, the most important issue when designing a local search is to identify effective neighborhood structure [40]. For the customers’ service options and vehicle routing in our problem, we design two kinds of neighborhood structures for service option and vehicle routing respectively.
For the service option neighborhood structure, two moves are included: they are HD to CP and CP to HD move, namely N1 and N2. More precisely, for N1, we randomly select a customer served by HD and identify its closest pickup point. If the pickup point is already selected, let the customer be served by it; otherwise, open the pickup point and serve the customer. For N2, we randomly select a pickup point and release one randomly selected customer and let him or her be served by HD.
When it comes to routing neighborhood structure, we take the nine route neighborhood moves proposed by Prins [45] which can comprehensively search the routing space.
Let T ( u ) stand for the trip visiting vertex u, x and y which are the successor of u and v in their respective trips. Define the neighborhood of vertex u as the h n c closest vertices, where h [ 0 , 1 ] is a granular threshold restricting the search to nearby vertices [48]. The following moves are concluded:
N3: If u is a customer or pickup point, remove u then insert it after v;
N4: If u is a customer or pickup point and x is a customer or dummy zero, remove them then insert ( u , x ) after v;
N5: If u is a customer or pickup point and x is a customer or dummy zero, remove them then insert ( x , u ) after v;
N6: If u and v are customers or pickup points, swap u and v;
N7: If u and v are customers or pickup points, and x is a customer or dummy zero, swap ( u , x ) and v;
N8: If u and v are customers, x and y are customers of dummy zeros, swap ( u , x ) and ( v , y ) ;
N9: If T ( u ) = T ( v ) , replace ( u , x ) and ( v , y ) by ( u , v ) and ( x , y ) ;
N10: If T ( u ) = T ( v ) , replace ( u , x ) and ( v , y ) by ( u , v ) and ( x , y ) ;
N11: If T ( u ) T ( v ) , replace ( u , x ) and ( v , y ) by ( u , y ) and ( x , v ) .
N3 to N5 correspond to insertions, whereas N6 to N8 are swaps. The aforementioned six moves can be worked on both the intra-route and inter-route. N9 is an intra-route 2-opt, while N10 and N11 are inter-route 2-opt. It can be seen that, compared to the initial moves in Prins [45], zeros are also used to formulate new neighborhood solutions, which may help to explore the search space.
When applying LS procedure, the above neighborhood moves are selected randomly to generate a new solution. The procedure terminates when either an improved solution is obtained or all neighborhood moves are checked without finding a better solution. The structure of LS can be seen in Figure 2.

4.8. Computational Complexity Analysis

Letting the scale of the considered problem be n d n p n c , Imax is the maximum iteration. In the presented algorithm, the computational complexity for population initialization is O ( n s p ( n p 2 + n p n c + n d n c + n c 2 ) ) . In each iteration, the computational complexity for decoding, crossover operation, mutation, evaluation and local search are O ( 2 n s + n c ) , O ( 2 n p ( 2 n s + n c ) ) , O ( 2 n s + n c ) , O ( M ( n s n c + n s n p + n p n c ) ) and O ( n p ( 2 n s + n c ) 2 ) , respectively. The computational complexity of the whole algorithm can be calculated as:
O ( n s p , M , n d , n p , n c , I max ) = O ( n s p ( n p 2 + n p n c + n d n c + n c 2 ) ) + O ( I max ( ( M + 4 ) ( 2 n s + n c ) + 2 n p ( 2 n s + n c ) + M ( n s n c + n s n p + n p n c ) + 4 n p ( 2 n s + n c ) 2 ) O ( 4 I max n p ( 2 n s + n c ) 2 )
As can be seen from Formula (24), the presented method belongs to polynomial algorithm and can be solved in polynomial time. More specifically, the computational complexity of the algorithm increases with the iterations and is proportional to np and (2ns+nc)2.

5. Computational Experiments

In this section, we perform numerical experiments to evaluate the overall performance of the proposed model and method. We first provide the implementation details of HGALS for the parameters setting. As the proposed problem belongs to a variant of LRP, the performance of the HGALS in solving LRP is then evaluated based on benchmark instances. Subsequently, a real-world instance is taken to analyze the effects of simultaneous HD and CP on the Last Mile distribution system. Finally, the validity of the presented method in solving the proposed problem is verified based on generated instances. All algorithms are coded in C language and compiled on Visual Studio 2013, implemented on a PC with an Intel Core Duo CPU at 1.74 GHz and 4GB memory under Window XP system.

5.1. Parameters Setting

To find suitable parameters, extensive experiments are conducted based on the generated instances depicted in Section 5.4. For GA-based heuristics, the most important parameter is the population size; with bigger M0, excellent solutions will be improved with much smaller probability, and it will take longer to find better solutions. Conversely, the population cannot sample the solution space which may lead to premature convergence. During the tests, since we found that the suitable initial population is closely related to np and variable population strategy helps a lot, it was finally decided as M0 = 4np and M = 2M0. According to the decided M0 and M, pc, pm, pls and pac are selected to be 0.45, 0.25, 1.0 and 0.25 respectively. Other parameters: h = 0.2, α1 = 1000, α2 = 200, α3 = 100, α4 = α5 = 10, γ1 = γ3 = 0.3, γ2 = 0.4.

5.2. Tests Based on Benchmark Instance

To show that the presented method is suitable to solve LRP, the benchmark instances proposed by Prins et al. [49] are used, which can be found at http://prodhonc.free.fr/Instances/instances_us.htm. The designed operators for pickup points are performed on the depots both for the solution initialization and mutation. For the terminate condition, we take the maximum iterations without improving the best solution. The specific value depends on instance size: it is 200 for small-sized instances (n ≤ 50), 400 for mid-sized instances (50 < n ≤ 100) and 600 for large-sized instances (n > 100). During the tests, each instance was independently run 10 times; the results obtained by the presented heuristic are provided in Table 2. For the notions used in the table, LB stands for the best known lower bound, BKR stands for the best known result, CPU stands for the average running time in seconds and Gap/BKR stands for the gap to BKR.
From Table 2, we can see that all the instances are solved with an overall average runtime of 277.2 s. Average gap to the best-known solution is 0.47%. For all instances containing 50 or less customers, the presented method is successful at finding almost all the optimal solutions. Moreover, for the instances with 100 and 200 customers, the average gap to the best-known solution is not changed dramatically.
Moreover, considering the most recent and effective methods in the literature, the performance comparison is provided in Table 3 based on the CPU and Best Gap/BKR. The compared algorithms are GRASP in Prins et al. [49], MAPM in Prins et al. [39], LRGTS in Prins et al. [50], GRASP+ELS in Duhamel et al. [36], SALRP in Vincent et al. [35], ALNS in reference [51], MACO in [32], GRASP+ILP in Contardo et al. [52], GVTNS in Escobar et al. [38] and Hybrid GA in Lopes et al. [43].
Among those compared algorithms, GRASP+ILP is on average the most effective on this set of benchmark instances followed by ALNS, GVTNS, MACO and SALRP. Our algorithm obtains similar results with SALRP and performs much better than GRASP, MAPM, LRGTS and GRASP+ELS. As CPU times depend on many factors (such as computers and programming languages used) [43], the proposed method seems faster than SALRP, ALNS and GRASP+ILP.
The computational results show that our algorithm can effectively solve the instances and is applicable to LRP.

5.3. Real-World Instance

This section presents the test based on a real-world Last Mile distribution system of Shapingba District in Chongqing city. The network consists of three depots (numbered from 1 to 3), 47 candidate pickup points (with number from 4 to 50) and 136 customers (with numbers from 51 to 186). The locations of these nodes are shown in Figure 5. Other main parameters are given as follows: Dpc = 600 m, uH = 1.5 min, uP = 10 min, pH = 10%, μsd = 1.5. For the algorithm terminate conditions, to get a better solution, we take the maximum iterations Imax and Tmax into consideration simultaneously; three terminate conditions are considered according to the instance scales: (Imax, Tmax) = (1 × 104, 600 s), (2 × 104, 900 s), (3 × 104, 1200 s). After running the algorithm, the obtained pickup point locations and customer assignments are reported in Table 4, and the vehicle routing information is given in Table 5.
To show the advantages of the proposed model, we now report the comparison of the proposed model with the scenario of only HD service, which can be seen in Table 6.
From the comparison, it can be clearly seen that although added pickup points cost by providing CP service, the total cost is reduced by 1050 and accounts for 17.6%. When analyzing the reason, we find that the vehicles are largely reduced by providing HD and CP services simultaneously, which saved the fixed vehicle cost and routing cost accordingly. Besides, the increased customers served by CP reduced the second delivery cost. These two factors work together, reducing the cost to a large extent.
In this research, we take the routing cost cij by dij multiplying the unit routing cost. From the comparison of the total routing cost of the two scenarios, we can deduce that the travelling distance is reduced by 33%. As the carbon emission is closely related to the travelling distance, the proposed model can help to reduce the carbon emission effectively.
As CP service provides an important role in reducing the operation cost, to have more insights, we conduct a sensitive analysis to show how the acceptable distance Dpc influences the cost. When Dpc ranges from 100 to 800 with step 100, the results are reported in Table 7 and Figure 6. It is worth mentioning that the scenario when Dpc = 0 is with only HD service.
From the sensitive analysis, we can see that when Dpc ≤ 100 m, the results obtained are the same with the only-HD scenario, which indicates that when Dpc ≤ 100 m, CP service shows no advantages when compared to HD and there is no need to provide CP service. When 100 < Dpc ≤ 600 m, there comes a continuous fall of the total cost, which means that with the distance increasing, better solutions are obtained due to the provided CP service. After that, the results stay unchanged which is constrained by the closest customers and the capacity of pickup points. From a cost optimization perspective, Dpc = 600 m is the most acceptable distance when designing the Last Mile distribution system in this instance.
Thus, we can draw the conclusion that, within appropriate acceptable distance Dpc, simultaneous HD and CP are effective at reducing the cost within a wide range of Dpc, and a suitable distance exists for each Last Mile distribution system.

5.4. Algorithm Components Performance Analysis

In this section, a computational study is carried out to compare our approach with standard GA&LS (SGALS) to show the effectiveness of the designed algorithm components in solving this problem. In the SGALS, random solution generation, objective function-based individual evaluation and classical single point crossover and mutation on the second tier of the chromosome are considered.
We first generate a set of different scales of instances based on the real-world instance parameters. For the instances, the number of depots, pickup points and customers range from 2 to 12, 10 to 30 and 50 to 200, respectively. The instances are presented as “I<depot number>-<pickup point number>-<customer number>”. The instances are generated as follows: depots and pickup points are randomly decided within the area with a radius of 15 km to reflect the relationship between the customers and pickup points; customers are randomly assigned within a radius of 1 km of selected pickup points; the capacity of depots as well as the vehicles are assigned with the requirement that all customers could be served; the capacity of pickup points are random value in [100, 120, 150, 180, 200] and the opening costs are [60, 80, 100, 120, 150] accordingly; customers’ demands are random value in [15, 20, 25, 30, 35]. Other parameters are the same to the real-world instance in Section 5.3. The two algorithms are run ten times independently for each instance; the results are reported in Table 8.
From Table 8 we can observe that our approach can always obtain the best solutions and run faster than SGALS. The results clearly indicate that the proposed components which constitute the algorithm are necessary and useful to enhance the performance both in the search ability and the stability when comparing the results obtained by SGALS. This can be seen by the improvement of 1.83% in Best Gap and 4.27% (7%–2.73%) in Average Gap. We can also draw the conclusion that the LRPSHC is very complex since the SGASA can obtain the same best solution as HGALS only for instance I2-18-80.

6. Conclusions

This paper proposed a novel location-routing problem for the Last Mile distribution for online shopping, in which two kinds of service options, home delivery and customer’s pickup, are simultaneously provided. For solving this specific problem, a hybrid heuristic combining GA and LS is presented together with several modified algorithm components.
Computational tests and comparative analysis with other published approaches based on benchmark instances are first conducted to evaluate the efficiency of the algorithm; results show that the presented algorithm framework is suitable to solve the LRP. Then, a real-world instance is adopted to verify the proposed model. By comparing the HD-only service scenario, the proposed model shows its great advantages in reducing the operation cost and vehicles, which proved that the model can help to design an economical and environmental Last Mile distribution system. A sensitive analysis is also made to give more insights into the effects of the acceptable distance for a customer to choose CP service based on the cost. Results show that the cost could be optimized by a wide scope of customer’s acceptable distances and an optimal distance for a certain distribution system. Finally, based on the generated instances, the designed algorithm components are proved to be effective in improving the search ability and stability of the algorithm when comparing with standard GA&LS.

7. Future Research Directions

In addition to delivering packages from the depots to customers, collecting packages from customers is another challenge for the Last Mile operation. Future research will focus on location-routing problem with simultaneous home delivery and customer’s pickup, as well as simultaneous delivery and collection, which may provide an even more efficient Last Mile solution. As LRPSHC is a new variant of LRP, an exact algorithm would be needed to solve this problem. Furthermore, it is clear that LRPSHC can be applied in many retail industry cases and it is also much more complex than classical LRP. More advanced metaheuristics are urgently needed in order to solve large-scale real-world instances effectively.

Acknowledgments

We appreciate the anonymous referees and the editors for their remarkable comments. This study is supported by National Science and Technology Support Program of China under Grant No. 2015BAH46F01, the Chongqing City Key Science Program Project under Grant No. cstc2014yykfA40006, cstc2015yykfC60002.

Author Contributions

Lin Zhou and Xu Wang proposed this research and created its framework; Lin Ni contributed to amalgamating the research and gave many useful suggestions; Yun Lin helped to collect data and analysis; Lin Zhou completed the paper. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. China’s E-Commerce Report 2015. Available online: http://finance.ifeng.com/a/20160629/14541580_0.shtml (accessed on 29 June 2016). (In Chinese)
  2. Gevaers, R.; Van de Voorde, E.; Vanelslander, T. Characteristics of innovations in last-mile logistics—Using best practices, case studies and making the link with green and sustainable logistics. In Proceedings of the European Transport Conference, Noordwijkerhout, The Netherlands, 5–7 October 2009.
  3. Esper, T.L.; Jensen, T.D.; Turnipseed, F.L.; Burton, S. The last mile: An examination of effects of online retail delivery strategies on consumers. J. Bus. Logist. 2003, 24, 177–203. [Google Scholar] [CrossRef]
  4. Aized, T.; Srai, J.S. Hierarchical modelling of Last Mile logistic distribution system. Int. J. Adv. Manuf. Technol. 2014, 70, 1053–1061. [Google Scholar] [CrossRef]
  5. Wang, X.; Zhan, L.; Ruan, J.; Zhang, J. How to choose “Last Mile” delivery modes for E-fulfillment. Math. Probl. Eng. 2014, 2014, 417129. [Google Scholar] [CrossRef]
  6. Agatz, N.A.; Fleischmann, M.; Van Nunen, J.A. E-fulfillment and multi-channel distribution—A review. Eur. J. Oper. Res. 2008, 187, 339–356. [Google Scholar] [CrossRef]
  7. Gevaers, R.; Van de Voorde, E.; Vanelslander, T. Cost Modelling and Simulation of Last-mile Characteristics in an Innovative B2C Supply Chain Environment with Implications on Urban Areas and Cities. Procedia Soc. Behav. Sci. 2014, 125, 398–411. [Google Scholar] [CrossRef]
  8. Boyer, K.K.; Prud’homme, A.M.; Chung, W. The last mile challenge: Evaluating the effects of customer density and delivery window patterns. J. Bus. Logist. 2009, 30, 185–201. [Google Scholar] [CrossRef]
  9. Hayel, Y.; Quadri, D.; Jiménez, T.; Brotcorne, L. Decentralized optimization of last-mile delivery services with non-cooperative bounded rational customers. Ann. Oper. Res. 2016, 239, 451–469. [Google Scholar] [CrossRef]
  10. Lopes, R.B.; Ferreira, C.; Santos, B.S.; Barreto, S. A taxonomical analysis, current methods and objectives on location-routing problems. Int. Trans. Oper. Res. 2013, 20, 795–822. [Google Scholar] [CrossRef]
  11. Prodhon, C.; Prins, C. A survey of recent research on location-routing problems. Eur. J. Oper. Res. 2014, 238, 1–17. [Google Scholar] [CrossRef]
  12. Drexl, M.; Schneider, M. A survey of variants and extensions of the location-routing problem. Eur. J. Oper. Res. 2015, 241, 283–308. [Google Scholar] [CrossRef]
  13. Vincent, F.Y.; Lin, S.-Y. A simulated annealing heuristic for the open location-routing problem. Comput. Oper. Res. 2015, 62, 184–196. [Google Scholar]
  14. Ponboon, S.; Qureshi, A.G.; Taniguchi, E. Branch-and-price algorithm for the location-routing problem with time windows. Transp. Res. E Logist. Transp. Rev. 2016, 86, 1–19. [Google Scholar] [CrossRef]
  15. Govindan, K.; Jafarian, A.; Khodaverdi, R.; Devika, K. Two-echelon multiple-vehicle location–routing problem with time windows for optimization of sustainable supply chain network of perishable food. Int. J. Prod. Econ. 2014, 152, 9–28. [Google Scholar] [CrossRef]
  16. Wang, H.; Du, L.; Ma, S. Multi-objective open location-routing model with split delivery for optimized relief distribution in post-earthquake. Transp. Res. E Logist. Transp. Rev. 2014, 69, 160–179. [Google Scholar] [CrossRef]
  17. Stenger, A.; Schneider, M.; Schwind, M.; Vigo, D. Location routing for small package shippers with subcontracting options. Int. J. Prod. Econ. 2012, 140, 702–712. [Google Scholar] [CrossRef]
  18. Karaoglan, I.; Altiparmak, F.; Kara, I.; Dengiz, B. The location-routing problem with simultaneous pickup and delivery: Formulations and a heuristic approach. Omega 2012, 40, 465–477. [Google Scholar] [CrossRef]
  19. Rieck, J.; Ehrenberg, C.; Zimmermann, J. Many-to-many location-routing with inter-hub transport and multi-commodity pickup-and-delivery. Eur. J. Oper. Res. 2014, 236, 863–878. [Google Scholar] [CrossRef]
  20. Rahmani, Y.; Cherif-Khettaf, W.R.; Oulamara, A. A local search approach for the two–echelon multi-products location–routing problem with pickup and delivery. IFAC-PapersOnLine 2015, 48, 193–199. [Google Scholar] [CrossRef]
  21. Vincent, F.Y.; Lin, S.-W. Multi-start simulated annealing heuristic for the location routing problem with simultaneous pickup and delivery. Appl. Soft Comput. 2014, 24, 284–290. [Google Scholar]
  22. Rahmani, Y.; Ramdane Cherif-Khettaf, W.; Oulamara, A. The two-echelon multi-products location-routing problem with pickup and delivery: Formulation and heuristic approaches. Int. J. Prod. Res. 2015, 54, 1–21. [Google Scholar] [CrossRef]
  23. Gianessi, P.; Alfandari, L.; Létocart, L.; Wolfler Calvo, R. The multicommodity-ring location routing problem. Transp. Sci. 2015, 50, 541–558. [Google Scholar] [CrossRef]
  24. Prodhon, C.; Prins, C. A memetic algorithm with population management (MA|PM) for the periodic location-routing problem. In Hybrid Metaheuristics; Springer: Berlin, Germany, 2008; pp. 43–57. [Google Scholar]
  25. Klibi, W.; Lasalle, F.; Martel, A.; Ichoua, S. The stochastic multiperiod location transportation problem. Transp. Sci. 2010, 44, 221–237. [Google Scholar] [CrossRef]
  26. Guerrero, W.J.; Prodhon, C.; Velasco, N.; Amaya, C.A. Hybrid heuristic for the inventory location-routing problem with deterministic demand. Int. J. Prod. Econ. 2013, 146, 359–370. [Google Scholar] [CrossRef]
  27. Rath, S.; Gutjahr, W.J. A math-heuristic for the warehouse location–routing problem in disaster relief. Comput. Oper. Res. 2014, 42, 25–39. [Google Scholar] [CrossRef]
  28. Tang, J.; Ji, S.; Jiang, L. The Design of a sustainable location-routing-inventory model considering consumer environmental behavior. Sustainability 2016, 8, 211. [Google Scholar] [CrossRef]
  29. Contardo, C.; Hemmelmayr, V.; Crainic, T.G. Lower and upper bounds for the two-echelon capacitated location-routing problem. Comput. Oper. Res. 2012, 39, 3185–3199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Nguyen, V.-P.; Prins, C.; Prodhon, C. Solving the two-echelon location routing problem by a GRASP reinforced by a learning process and path relinking. Eur. J. Oper. Res. 2012, 216, 113–126. [Google Scholar] [CrossRef]
  31. Nguyen, V.-P.; Prins, C.; Prodhon, C. A multi-start iterated local search with tabu list and path relinking for the two-echelon location-routing problem. Eng. Appl. Artif. Intell. 2012, 25, 56–71. [Google Scholar] [CrossRef]
  32. Ting, C.-J.; Chen, C.-H. A multiple ant colony optimization algorithm for the capacitated location routing problem. Int. J. Prod. Econ. 2013, 141, 34–44. [Google Scholar] [CrossRef]
  33. Ardjmand, E.; Weckman, G.; Park, N.; Taherkhani, P.; Singh, M. Applying genetic algorithm to a new location and routing model of hazardous materials. Int. J. Prod. Res. 2015, 53, 916–928. [Google Scholar] [CrossRef]
  34. Ardjmand, E.; Young, W.A.; Weckman, G.R.; Bajgiran, O.S.; Aminipour, B.; Park, N. Applying genetic algorithm to a new bi-objective stochastic model for transportation, location, and allocation of hazardous materials. Expert Syst. Appl. 2016, 51, 49–58. [Google Scholar] [CrossRef]
  35. Vincent, F.Y.; Lin, S.-W.; Lee, W.; Ting, C.-J. A simulated annealing heuristic for the capacitated location routing problem. Comput. Ind. Eng. 2010, 58, 288–299. [Google Scholar]
  36. Duhamel, C.; Lacomme, P.; Prins, C.; Prodhon, C. A GRASP×ELS approach for the capacitated location-routing problem. Comput. Oper. Res. 2010, 37, 1912–1923. [Google Scholar] [CrossRef]
  37. Escobar, J.W.; Linfati, R.; Toth, P. A two-phase hybrid heuristic algorithm for the capacitated location-routing problem. Comput. Oper. Res. 2013, 40, 70–79. [Google Scholar] [CrossRef]
  38. Escobar, J.W.; Linfati, R.; Baldoquin, M.G.; Toth, P. A Granular Variable Tabu Neighborhood Search for the capacitated location-routing problem. Transp. Res. B Methodol. 2014, 67, 344–356. [Google Scholar] [CrossRef]
  39. Prins, C.; Prodhon, C.; Calvo, R.W. A memetic algorithm with population management (MA|PM) for the capacitated location-routing problem. In Evolutionary Computation in Combinatorial Optimization, Proceedings of the 6th European Conference, EvoCOP 2006, Budapest, Hungary, 10–12 April 2006; Springer: Berlin, Germany, 2006; pp. 183–194. [Google Scholar]
  40. Derbel, H.; Jarboui, B.; Hanafi, S.; Chabchoub, H. Genetic algorithm with iterated local search for solving a location-routing problem. Expert. Syst. Appl. 2012, 39, 2865–2871. [Google Scholar] [CrossRef]
  41. Montoya-Torres, J. R.; López Franco, J.; Nieto Isaza, S.; Felizzola Jiménez, H.; Herazo-Padilla, N. A literature review on the vehicle routing problem with multiple depots. Comput. Ind. Eng. 2015, 79, 115–129. [Google Scholar] [CrossRef]
  42. Marinakis, Y.; Marinaki, M. A bilevel genetic algorithm for a real life location routing problem. Int. J. Logist. Res. Appl. 2008, 11, 49–65. [Google Scholar] [CrossRef]
  43. Lopes, R.B.; Ferreira, C.; Santos, B.S. A simple and effective evolutionary algorithm for the capacitated location–routing problem. Comput. Oper. Res. 2016, 70, 155–162. [Google Scholar] [CrossRef]
  44. Lin, S.; Kernighan, B.W. An effective heuristic algorithm for the traveling-salesman problem. Oper. Res. 1973, 21, 498–516. [Google Scholar] [CrossRef]
  45. Prins, C. A simple and effective evolutionary algorithm for the vehicle routing problem. Comput. Oper. Res. 2004, 31, 1985–2002. [Google Scholar] [CrossRef]
  46. Vidal, T.; Crainic, T.G.; Gendreau, M.; Lahrichi, N.; Rei, W. A hybrid genetic algorithm for multidepot and periodic vehicle routing problems. Oper. Res. 2012, 60, 611–624. [Google Scholar] [CrossRef]
  47. Ahmadizar, F.; Zeynivand, M.; Arkat, J. Two-level vehicle routing with cross-docking in a three-echelon supply chain: A genetic algorithm approach. Appl. Math. Model. 2015, 39, 7065–7081. [Google Scholar] [CrossRef]
  48. Toth, P.; Vigo, D. The granular tabu search and its application to the vehicle-routing problem. Inf. J. Comput. 2003, 15, 333–346. [Google Scholar] [CrossRef]
  49. Prins, C.; Prodhon, C.; Calvo, R.W. Solving the capacitated location-routing problem by a GRASP complemented by a learning process and a path relinking. 4OR 2006, 4, 221–238. [Google Scholar] [CrossRef]
  50. Prins, C.; Prodhon, C.; Ruiz, A.; Soriano, P.; Wolfler Calvo, R. Solving the capacitated location-routing problem by a cooperative Lagrangean relaxation-granular tabu search heuristic. Transp. Sci. 2007, 41, 470–483. [Google Scholar] [CrossRef]
  51. Hemmelmayr, V.C.; Cordeau, J.-F.; Crainic, T.G. An adaptive large neighborhood search heuristic for two-echelon vehicle routing problems arising in city logistics. Comput. Oper. Res. 2012, 39, 3215–3228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Contardo, C.; Cordeau, J.-F.; Gendron, B. An exact algorithm based on cut-and-column generation for the capacitated location-routing problem. INFORMS J. Comput. 2013, 26, 88–102. [Google Scholar] [CrossRef]
Figure 1. The Last Mile distribution system with simultaneous HD and CP services.
Figure 1. The Last Mile distribution system with simultaneous HD and CP services.
Sustainability 08 00828 g001
Figure 2. The general structure of HGALS.
Figure 2. The general structure of HGALS.
Sustainability 08 00828 g002
Figure 3. Solution representation.
Figure 3. Solution representation.
Sustainability 08 00828 g003
Figure 4. Crossover operation.
Figure 4. Crossover operation.
Sustainability 08 00828 g004
Figure 5. Network of the instance.
Figure 5. Network of the instance.
Sustainability 08 00828 g005
Figure 6. Sensitive analysis of Dpc on cost.
Figure 6. Sensitive analysis of Dpc on cost.
Sustainability 08 00828 g006
Table 1. Notions used in the formulation.
Table 1. Notions used in the formulation.
SetsDescription
G Network, G = ( N , A , C )
N Node set, N = N D N P N C
N D Depot set, N D = { 1 , , n d }
N P Pickup point set, N P = { n d + 1 , , n d + n p }
N C Customer set, N C = { n d + n p + 1 , , n d + n p + n c }
A Arc set, A = { ( i , j ) : i , j N }
C Connection set, C = { ( p , c ) : p N P , c N C }
V D Available vehicle set
Parameters
n d Number of depots
n p Number of pickup points
n c Number of customers
d i j Distance of arc ( i , j ) A
t i j Travel time of arc ( i , j ) A
c i j Routing cost of arc ( i , j ) A
q i Demand of customer i N C
u H Unit demand service time for HD service
u P Fixed service time at a pickup point
μ s d Probability of the second delivery
B p Capacity of pickup point p N P
S d Capacity of depot d N D
U p Opening cost of pickup point p N P
Q v Vehicle capacity
F v Fixed vehicle cost
n v Available vehicle number
D p c Acceptable distance for CP service
Decision variables
x i j k d Equal to 1 if vehicle k V D departs from depot d N D that passes through arc ( i , j ) A ; 0, otherwise
y p Equal to 1 if pickup point p N P is selected to serve customers; 0, otherwise
z i p Equal to 1 if customer i N C is served by pickup point p N P ; 0, otherwise
w p d Equal to 1 if pickup point p N P is served by depot d N D ; 0, otherwise
φ i j k d Freight transported from depot d N D by vehicle k V D that passes through arc ( i , j ) A
Table 2. Results for the benchmark instances.
Table 2. Results for the benchmark instances.
InstanceLBBKRCPUAverageBest
ResultGap/BKRResultGap/BKR
120-5-154,79354,7931.854,8360.0854,7930
220-5-1b39,10439,1042.239,104039,1040
320-5-2a48,90848,9081.548,908048,9080
420-5-2b37,54237,5422.637,542037,5420
550-5-190,11190,11112.990,1300.0390,1110
650-5-1b67,34063,24217.863,242063,2420
750-5-288,29888,29824.388,6430.4088,6430.39
850-5-267,34067,34021.467,340067,3400
950-5-2bis84,05584,05520.884,1390.1084,0550
1050-5-2bbis51,82251,82216.251,9580.2651,8220
1150-5-386,20386,20319.386,4560.2986,4560.29
1250-5-3b61,83061,83018.161,830061,8300
Avg 13.2 0.10 0.06
13100-5-1275,993274,81492.1277,1350.84277,0350.81
14100-5-1b214,392213,615122.42,145,0340.42214,3130.33
15100-5-2194,598193,67191.8194,3660.36194,1240.23
16100-5-2b157,173157,095132.5157,5600.3157,0950
17100-5-3200,246200,079124.2201,8440.88201,6280.77
18100-5-3b152,586152,441107.7154,4841.34152,9920.36
Avg 111.8 0.69 0.42
19100-10-1290,429287,695113.3291,3391.27290,2430.89
20100-10-1b234,641230,989127.5235,0571.76233,5121.09
21100-10-2244,265243,59095.6249,1802.29245,1230.63
22100-10-2b203,988203,988130.9205,5110.74204,6670.33
23100-10-3253,344250,882133.5257,0592.46253,8651.19
24100-10-3b204,597204,317122.7205,4490.55205,2320.48
Avg 120.6 1.51 0.76
25200-10-1479,425475,294877.1484,7481.99479,9260.97
26200-10-1b378,773377,043922.3382,5951.47380,6130.94
27200-10-2450,468449,006823.6451,8350.63450,3120.29
28200-10-2b374,435374,280789.4380,6671.70375,6740.37
29200-10-3472,898469,433865.7476,2301.45473,8750.94
30200-10-3b364,178362,653900.5364,7620.58364,2010.43
Avg 863.1 1.30 0.66
Global Avg 277.2 0.77 0.47
Note: The best known solution values attained are marked in bold.
Table 3. Comparison results of the methods.
Table 3. Comparison results of the methods.
AlgorithmPerformanceAlgorithmPerformance
CPUBest Gap/BKRCPUBest Gap/BKR
GRASP96.53.57MACO176.40.40
MAPM76.71.35GRASP+ILP1163.00.12
LRGTS17.50.70GVTNS91.20.37
GRASP+ELS258.21.11Hybrid GA199.10.32
SALRP422.40.46HGALS277.20.47
ALNS4221.00.27
Table 4. Results of the instance—pickup point locations and customer assignments.
Table 4. Results of the instance—pickup point locations and customer assignments.
Pickup PointCustomerPickup PointCustomerPickup PointCustomer
554, 55, 56660, 61, 62, 65858, 59, 63, 64
1069, 75, 831173, 74, 76,771271,72
1585,861791, 92, 9322137, 138
23144, 14524141, 142, 14825146,147
26151, 152, 153, 15428132, 14929157, 159
31158, 160, 161, 162, 16533172, 17334167, 168, 170
35101, 102, 10337104, 105, 106, 1073897, 98
3995, 10840110, 111, 11241114, 115
44119, 12045122, 12346128, 129
47174, 17548177, 178, 179, 18049181, 812, 183
Table 5. Results of the instance—vehicle routing.
Table 5. Results of the instance—vehicle routing.
No.RoutingNo.Routing
11-11-12-10-67-68-172-156-26-25-143-140-136-22-23-2
21-5-51-53-6-182-31-29-163-164-166-169-34-171-2
31-8-57-66-87-84-82-192-47-176-186-185-184-49-48-2
41-70-135-139-24-150-28-133-134-81-1103-40-113-116-121-124-45-117-41-3
51-15-88-89-90-96-38-37-99-79-1113-118-44-99-125-109-39-94-17-3
61-80-35-127-126-130-156-131-46-100-78-1
Table 6. Comparison of the two scenarios.
Table 6. Comparison of the two scenarios.
ItemSimultaneous HD and CPOnly HDDifference
Vehicle number1126−15
Pickup point number30-+30
Fixed vehicle cost22005200−3000
Routing cost156233−77
Pickup point opening cost2370-+2370
Second delivery cost186529−343
Total cost49125962−1050
Table 7. Results of sensitive analysis of Dpc on cost.
Table 7. Results of sensitive analysis of Dpc on cost.
ItemDpc
100200300400500600700800
Vehicle number2621131111111111
Pickup point number-16363330303030
Fixed vehicle cost52004000260022002200220022002200
Routing cost233225109157152156156156
Pickup point opening cost-1120264024802400237023702370
Second delivery cost529417140188209186186186
Total cost59625762548950254961491249124912
Table 8. Comparison results of the two algorithms.
Table 8. Comparison results of the two algorithms.
InstanceSGALSHGALS
CPUAverageBestCPUAverageBest
ResultGapResultGapResultGapResultGap
I2-10-50272.91833.21.691806.20.19143.21804.50.991802.70
I2-15-60533.72473.24.562382.20.71413.82388.60.982365.30
I2-18-80600.13315.76.373117.10435.63202.72.753117.10
I2-20-100900440414.464098.66.5795.13958.62.883847.80
I2-22-120862.34471.46.814220.70.82852.44276.12.154186.30
I2-25-1509004967.86.134787.12.279004769.31.894680.80
I2-28-180118371669.217074.47.8212006710.62.286561.30
I2-30-20012006748.76.7465433.4912006426.91.656322.60
Agv807 7.0 2.73742 1.83 0
Note: The best known solution values attained are marked in bold.

Share and Cite

MDPI and ACS Style

Zhou, L.; Wang, X.; Ni, L.; Lin, Y. Location-Routing Problem with Simultaneous Home Delivery and Customer’s Pickup for City Distribution of Online Shopping Purchases. Sustainability 2016, 8, 828. https://0-doi-org.brum.beds.ac.uk/10.3390/su8080828

AMA Style

Zhou L, Wang X, Ni L, Lin Y. Location-Routing Problem with Simultaneous Home Delivery and Customer’s Pickup for City Distribution of Online Shopping Purchases. Sustainability. 2016; 8(8):828. https://0-doi-org.brum.beds.ac.uk/10.3390/su8080828

Chicago/Turabian Style

Zhou, Lin, Xu Wang, Lin Ni, and Yun Lin. 2016. "Location-Routing Problem with Simultaneous Home Delivery and Customer’s Pickup for City Distribution of Online Shopping Purchases" Sustainability 8, no. 8: 828. https://0-doi-org.brum.beds.ac.uk/10.3390/su8080828

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop