Next Article in Journal
Design Optimisation of Bi-Cruciate Retaining Total Knee Arthroplasty (TKA) Prosthesis via Taguchi Methods
Next Article in Special Issue
Stylized Model of Lévy Process in Risk Estimation
Previous Article in Journal
Fixed-Point Control of Airships Based on a Characteristic Model: A Data-Driven Approach
Previous Article in Special Issue
Metaheuristics in Business Model Development for Local Tourism Sustainability Enhancement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Vehicle Routing Problem with Simultaneous Pickup and Delivery Considering the Total Number of Collected Goods

1
School of Management, Xi’an Jiaotong University, Xi’an 710049, China
2
School of Economics and Management, Xi’an Technological University, Xi’an 710021, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(2), 311; https://doi.org/10.3390/math11020311
Submission received: 19 October 2022 / Revised: 3 January 2023 / Accepted: 5 January 2023 / Published: 6 January 2023
(This article belongs to the Special Issue Planning and Scheduling in City Logistics Optimization)

Abstract

:
As a consequence of e-commerce development, large quantities of returned goods are shipped every day. The vehicle routing problem with simultaneous delivery and pickup (VRPSDP) has become one of the most important areas of logistics management. Most related studies are aimed at minimizing travel time. However, the total number of collected goods is also very important to logistics companies. Thus, only considering the traveling time cannot reflect actual practice. To effectively optimize these operations for logistics companies, this paper introduces the vehicle routing problem with simultaneous pickup and delivery considering the total number of collected goods. Based on the principles of considering the number of collected goods, a bi-objective vehicle routing model minimizing the total travel time and maximizing the total number of collected goods simultaneously is developed. A polynomial time approximation algorithm based on the ε -constraint method is designed to address this problem, and the approximation ratio of the algorithm is analyzed. Finally, the validity and feasibility of the proposed model and algorithm are verified by test examples, and several managerial insights are derived from the sensitivity analysis.

1. Introduction

The increasingly mature e-commerce has promoted the rapid development of logistics distribution industry, and the scale of reverse logistics caused by return has also increased. According to the statistics of the National Retail Federation of the United States, the return rate of e-commerce reached 20–30%. In December 2019, UPS alone collected about 1 million return packages every day. During China’s double 11 in 2018, as of 13 November, Tmall returned 9.27 billion yuan’s commodities. The high return volume makes the logistics and distribution enterprises to provide distribution and recycling services to customers at the same time. The traditional vehicle routing problem (VRP) cannot describe the characteristics of simultaneous pick-up and delivery in the logistics and distribution process, so an increasing number of scholars have turned to its branch problem-the vehicle routing problem with simultaneous delivery and pickup (VRPSDP).
The VRPSPD refers to the situation where the vehicle starts from the distribution center to deliver goods at the designated demand point and also needs to bring the goods at the demand point back to the distribution center. How to arrange the vehicle route so that the vehicle can complete the pickup service request at the same time of delivery, or the transportation cost is the smallest when all the pickup service requests are completed. It was first proposed by Min in 1989 [1]. The optimization problem of the simultaneous pickup and delivery vehicle routing between a book center and 22 local libraries with the goal of minimizing distribution costs was studied. Its path planning has a great influence on distribution efficiency. Subsequently, this problem has attracted the attention of an increasing number of researchers. About 16.3% of the vehicle routing problem (VRP) research papers published from 2009 to 2017 belonged to the VRPSDP [2].
Almost all the studies aimed at minimizing the travel time and modeling and designing algorithms for the VRPSDP based on known or uncertain pickup demand and single or multiple vehicles. However, in practice, the total number of collected goods is also very important to logistics companies and should be considered in route optimization. In this paper, according to the actual demand and the lack of existing research, under vehicle capacity constraints and considering the travel time and the total number of collected goods at the same time, a bi-objective VRPSDP optimization model is established to minimize the travel time and maximize the total number of collected goods. An approximate algorithm is designed to solve the problem, and an example is given to verify the model and algorithm.
The rest of this study is organized as follows. Section 2 describes a literature review on the VRPSDP with known pickup demand and uncertain pickup demand. Section 3 defines the bi-objective mathematical VRPSDP model, including the description of the problem and the model. The model analysis and solution are introduced in Section 4 and Section 5 present the results of the testing example. Finally, Section 6 provides concluding statements on our work.

2. Literature Review

2.1. The VRPSPD with Known Pickup Demand

The VRPSDP problem was first proposed by Min in 1989 [1]. Some accurate algorithms had been designed to solve this problem, such as the branch-and-price method [3], the branch-and-cut method [4,5,6,7], and a combination of the branch-and-cut-and-price method [8,9].
Because the number of feasible solutions increases factorial with the increase of demand points, there are only a few accurate algorithms to solve the problem. The heuristic algorithm and meta-heuristic algorithm are the main algorithms to solve the simultaneous pickup and delivery vehicle routing problem. Dethloft [10] proposed a heuristic algorithm based on insertion and test instances to solve real-life problems. Kaya et al. developed a hybrid metaheuristic algorithm based on an ant colony system (ACS) and a variable neighborhood search (VNS) for VRPSDP’s solution [11]. Zachariadis solved this problem via a common local-search method for the routing aspects, integrated with a packing heuristic for developing feasible item loadings [12]. Simsir et al. established the mathematic model of VRPSDP with a single distribution center and a single type of vehicle, and constructed the artificial bee swarm algorithm to solve it [13].
Li [14] developed a metaheuristic method based on an iterated local search for the multi-depot vehicle routing problem with simultaneous deliveries and pickups (MDVRPSDP). Hs et al. [15] also discussed the VRPSDP with multiple distribution centers, established a model, and solved it. Avic [16] proposed the heterogeneous vehicle routing problem with a simultaneous pickup and delivery (HVRPSPD) model and developed a hybrid local search algorithm in which a non-monotone threshold adjusting strategy is integrated with tabu search. Zhang [17] addressed a multi-commodity many-to-many vehicle routing problem with simultaneous pickup and delivery (M-M-VRPSPD) for a fast fashion retailer in Singapore. Different from other widely studied pickup and delivery problems, the unique characteristics are collecting products from customers were encouraged to be reallocated to fulfill the demands of other customers. Wang [18] formulated and solved an optimization of multi-depot pickup and delivery problems with resource sharing.
In order to solve the VRPSDP with a time window, Wang [19] proposed a particle swarm optimization algorithm and verified its effectiveness, and the same problem was solved by and Tabu search algorithm [20].
As can be seen above, for the VRPSDP when the demand is known, scholars established the model based on the number of distribution centers, the type of vehicles, time windows et al., mostly with the goal of minimizing the transportation cost. The static optimization method was used to solve the problem. Focus on designing different algorithms to find the approximate optimal solution, trying to find a more stable and efficient algorithm [21].

2.2. The VRPSPD with Uncertain Pickup Demand

Ignoring the uncertainty in the transportation process is the main reason for the decrease in service level and the increase in cost, so some scholars have studied the situation of uncertain pickup demand. Liu et al. [22] constructed a mixed integer programming model for the VRPSDP with a single distribution center and a single type of vehicle, in which the pickup demand follows a random binomial distribution, with the goal of minimizing the expected travel time. Then, they designed a hybrid ant colony system optimization algorithm. Peng [23] established a VRPSDP model of a single distribution center and a single type of vehicle with a normal pickup distribution. Jia et al. simplified the problem to an asymmetric traveling salesman problem with a time window and solved it using a rolling time-domain scheduling algorithm [24]. Hu et al. [25] studied the VRPSDP with incompatible goods and uncertain pickup demand and handled the uncertain pickup demand through two stages: The first stage obtained the pickup demand through a probability distribution, resulting in an a priori path; and the second stage introduced the vehicle loading rate control parameters to judge and solve the service demand points. Based on pre-optimization and re-optimization strategies, Fan et al. [26] constructed a two-stage vehicle routing problem model with uncertain pickup demand. In the pre-optimization stage, the vehicle was assigned to the demand point based on a random chance constraint mechanism and vehicle loading constraint. The pre-optimization scheme was generated. In the re-optimization stage, the route of the failure point and its subsequent demand point was adjusted. Finally, a hybrid variable neighborhood search algorithm was designed to solve the problem. On this basis, the VRPSDP with random demand under the multicenter joint distribution mode of multiple distribution centers was further studied, and a two-stage path optimization model was constructed [27]. According to the characteristics of the problem, an adaptive variable neighborhood cultural gene algorithm was designed to solve the problem. With the goal of minimizing travel costs, scholars have modeled and designed algorithms to solve the situation in which single distribution centers, multi-vehicle, pick-up demand, and delivery demand were all random [28,29].
Zhou [30] studied the optimization model and algorithm of the multi-depot VRPSDP under fuzzy demand. Ma [31] assumed that the pickup demand is a fuzzy random variable, established the VRPSDP model, and designed a genetic algorithm based on fuzzy random simulation to solve it. Combined with credibility theory, Fan et al. [32] modeled the VRPSDP with fuzzy pickup demand and finally transformed it into a deterministic problem to solve.
In summary, it was found that in the VRPSDP with uncertain demands, most of the existing studies randomized or fuzzed the uncertain quantity of pickup, assuming that the quantity of pickup at the demand point obeys a certain distribution (probability distribution or fuzzy distribution). Then find out the minimum expected transportation cost and seek the optimal scheme in the average sense.
Almost all of the mentioned studies aimed at minimizing travel time and modeling and designing algorithms for the VRPSDP. Since the total number of collected goods is very important to logistics companies and should be considered in route optimization. However, no such total number of collected goods-related models has been proposed till now.

3. Problem Description and Modeling

3.1. Problem Description and Symbolic Description

Let G = ( V , E ) be an undirected complete graph, with vertex set V = { v 0 } V 0 and edge set E = e ( v i , v j ) | v i , v j V , where vertex v 0 represents the distribution center and V 0 = v 1 , v 2 , v n represents the demand point set of G . The distribution center v 0 has l vehicles with a maximum carrying capacity of Q , and the vehicles are sufficient to complete the distribution task. t i j is the time when the vehicle passes through edge e v i , v j , and d i and p i are the delivery and pickup quantities at demand point v i , respectively. Q i j k is the loading volume of vehicle k from demand point v i to v j . Now, we need to arrange vehicles to start from the distribution center v 0 to deliver the goods to the designated demand point and determine whether to return goods to v 0 according to the situation. The delivery demand and pickup demand of each demand point are known before departure, and the demand to be pickup cannot be split. The problem is how to determine the vehicle route and whether to pick up goods after reaching the demand point while not exceeding the carrying capacity of the vehicle and completing the delivery task so as to make the travel time as small as possible and the total number of goods collected as much as possible.
In order to facilitate the follow-up discussion, the decision variables required for modeling are given first.
x i j k = 1 if vehicle k travels directly from demand point v i to v j , and 0 otherwise. y i k = 1 if vehicle k collects goods at demand point v i , and 0 otherwise.

3.2. Modeling

Using the above-defined parameters and variables, we present the formulation of a bi-objective mathematical VRPSDP model.
f 1 = max i = 1 n k = 1 l p i y i k
f 2 = min i = 0 n j = 0 n k = 1 l t i j x i j k
subject   to   i = 0 n k = 1 l x i j k = 1 , ( j = 1 , 2 , , n )
i = 0 n x i j k = i = 0 n x j i k , ( j = 1 , 2 , , n ; k = 1 , 2 , , l )
i = 1 n x i 0 k = j = 1 n x 0 j k , ( k = 1 , 2 , , l )
Q i j k Q , ( i = 0 , 1 , 2 , , n ; j = 0 , 1 , 2 , , n ; k = 1 , 2 , , l )
i = 0 n Q i j k d j + p j y j k = i = 0 n Q j i k , ( j = 1 , 2 , , n ; k = 1 , 2 , , l )
x i j k y i k , ( i = 1 , 2 , , n ; j = 1 , 2 , , n ; k = 1 , 2 , , l )
x i j k = 0 o r 1 , ( i = 0 , 1 , , n ; j = 0 , 1 , 2 , , n ; k = 1 , 2 , , l )
y i k = 0   o r   1 , ( i = 1 , 2 , , n ; k = 1 , 2 , , l )
Objective function (1) maximizes the total number of goods collected, and objective function (2) minimizes the total travel time. Constraint (3) is a degree constraint ensuring that each demand point is visited at most once, constraint (4) preserves the flow conservation, constraint (5) indicates that the vehicles start from the distribution center v 0 and finally return to it, constraint (6) is a loading restraint, constraint (7) represents the loading balance after visiting demand points v j , constraint (8) indicates that the vehicle must have a path connection when serving any demand point, and constraints (9) and (10) impose that all variables be 0–1.

4. Model Analysis and Solution

The VRPSDP model, as defined by Equations (1)–(10), is a bi-objective mixed integer optimization programming model. For the bi-objective optimization problem, the ε -constraint method proposed by Bérubé can find the complete Pareto frontier [33]. The core idea of this method is to retain only one goal and transform other goals into constraints limited by the value of ε . Therefore, according to the basic idea of the ε -constraint method, the original problem is transformed into a single-objective problem P0 as follows.
f 2 = min i = 0 n j = 0 n k = 1 l t i j x i j k
subject to constraints (3)–(10):
i = 1 n k = 1 l p i y i k ε
This paper designed an approximate solution algorithm based on the ε -constraint method. The approximate algorithm first needed to find the TSP (traveling salesman problem) Hamiltonian cycle of the original problem and then arranged the first vehicle to distribute along the TSP Hamiltonian cycle. The vehicle visited one or several demand points in turn. If the vehicle could not meet the following conditions, it returned to the distribution center v 0 and let the next vehicle continue to deliver from v 0 . Condition 1 is that the remaining loading capacity of the vehicle is sufficient to meet the demand of the next customer for picking up and delivering goods; and condition 2 is that if the vehicle visits the next demand point, the actual load after distribution does not exceed its maximum load.
According to the above algorithm idea, the specific steps of the approximate algorithm are given as follows, where U k represents the loading capacity when vehicle k departs from depot v 0 .

Time Complexity Analysis of Algorithm 1

The time complexity of using the Christofides algorithm in step 1 is o ( n 3 ) . Steps 2 and 3 need to be calculated and compared at most n times. The maximum number of cycles in steps 3–4 is n , so the time complexity is o ( n 2 ) . The maximum number of cycles in steps 2–4 is Q , so the time complexity is o ( Q n 3 ) . In summary, the time complexity of Algorithm 1 is o ( n 3 ) .
The algorithm has been designed, and the time complexity has been analyzed. Next, we analyze the approximation ratio of Algorithm 1. The approximate ratio refers to the ratio between the objective function of the approximate solution and the exact solution. The closer its value is to 1, the closer the solution obtained by the algorithm is to the optimal solution. For example, for a cost minimization problem, let c ( A ) represent the cost corresponding to Algorithm A and c ( O P T ) represent the cost corresponding to the optimal solution. If c ( A ) α c ( O P T ) is true, then α is called the approximate ratio of Algorithm A.
Algorithm 1: The approximate algorithm
STEP 1 Calculate the TSP Hamiltonian cycle through the Christofides Algorithm [34] and record the demand points on the TSP Hamiltonian cycle as   v 1 * , v 2 * , , v n * .
STEP 2   Let   k = 1 .   The   shortest   path   S P 01 *   of   v 0   to   v 1 *   is   edge   e ( v 0 , v 1 * ) .   Let   P k = e ( v 0 , v 1 * ) ,   U k = d 1 * ,   and   i = 1 .
STEP 3   Let   i = i + 1 .   If   i > n ,   proceed   to   STEP 4 .   Otherwise ,   if Q j = 1 i 1 d j * + j = 1 i 1 p j * p i * d i *   and   U k + d i * Q ,   let   U k = U k + d i * ,   add   edge   e ( v i 1 * , v i * )   to   path   P k ,   and   repeat   STEP   3 .   Otherwise ,   add   edge   e ( v i 1 * , v 0 )   to   path P k ;   let   k = k + 1 , U k = q i * , and   P k = S P 0 i * ; and return to STEP 3.
STEP 4 Add   edge   e ( v i 1 * , v 0 )   to   path   P k ,   record   the   total   driving   time   T = k = 1 l t ( P k )   and   the   total   number   of   goods   retrieved   D = k = 1 l p ( P i ) , and output the corresponding distribution path of each vehicle.
First, the transportation time corresponding to the optimal solution of the vehicle routing problem is analyzed. The following lemma is derived from the literature [35].
Lemma 1.
The lower bound of the transportation time corresponding to the optimal solution of the vehicle routing problem is max { i = 1 n 2 d 0 i q i Q , T * ( T S P ) } , where T * ( T S P ) is the travel time of the optimal TSP Hamiltonian cycle.
Lemma 1 gives the lower bound of the travel time corresponding to the optimal solution of the vehicle routing problem. Because Problem P0 has additional constraints compared with the VRP problem, its travel time must not be less than that of the general VRP problem. Therefore, the lower bound given by Lemma 1 is also the lower bound of the travel time of Problem P0. Combined with Lemma 1, the following theorem is given.
Theorem 1.
The A* approximation ratio of the algorithm for solving Problem P0 is i = 1 n Q d 0 i i = 1 n d 0 i q i + β , where d 0 i represents the travel time from depot v 0 to demand point v i , and β is the approximation ratio of the approximate algorithm for solving the TSP path.
Proof of Theorem 1.
The solution obtained by Algorithm 1 consists of part of the TSP Hamiltonian cycle and several paths connecting the distribution center v 0 and the demand point. Here, the part related to the TSP path is called Part 1, and the path connecting the distribution center v 0 and the demand point is called Part 2. For Part 1, the overall travel time will not exceed the travel time of the entire TSP path. That is, the maximum transportation cost is T ( S T P ) . The maximum value of Part 2 is: i = 1 n d 0 i + d i 0 = i = 1 n 2 d 0 i . Then, α = c ( A * ) c ( O P T ) = i 1 n 2 d 0 i + T ( T S P ) max { i = 1 n 2 d i 0 q i Q , T * ( T S P ) } = i 1 n 2 d 0 i + T ( T S P ) max { i = 1 n 2 d 0 i q i Q , T * ( T S P ) } . If β represents the approximation ratio of the approximate algorithm for solving the TSP path, then α = i = 1 n 2 d 0 i + β T * ( T S P ) max { i = 1 n 2 d 0 i q i Q , T * ( T S P ) } . Next, compare the size of i = 1 n 2 d 0 i q i Q and T * ( T S P ) . The following two situations exist. If i = 1 n 2 d 0 i q i Q T * ( T S P ) , then α = i = 1 n 2 d 0 i + β T * ( T S P ) i = 1 n 2 d 0 i q i Q = i = 1 n Q d 0 i i = 1 n d 0 i q i + β T * ( T S P ) i = 1 n 2 d 0 i q i Q i = 1 n Q d 0 i i = 1 n d 0 i q i + β T * ( T S P ) T * ( T S P ) = i = 1 n Q d 0 i i = 1 n d 0 i q i + β . If i = 1 n 2 d 0 i q i Q T * ( T S P ) , then α = i = 1 n 2 d 0 i + β T * ( T S P ) T * ( T S P ) = i = 1 n 2 d 0 i T * ( T S P ) + β T * ( T S P ) T * ( T S P ) i = 1 n 2 d 0 i i = 1 n 2 d 0 i q i Q + β = i = 1 n Q d 0 i i = 1 n d 0 i q i + β .
Combining the above two cases, Theorem 1 is proved. □

5. Example

In this section, the performance of the solution method is evaluated. The proposed method in this paper was coded with Java. All experiments were performed on a PC with a 3.3-GHz processor and 8.0-GB RAM under Windows 7.

5.1. Testing Example

In this part, Solomon VRP test set R101 was adopted and modified to generate VRPSDP examples. The first 10 demand points in the database, where node 1, the distribution center, was selected; and the node data are shown in the following table. Without any loss of generality, the vehicle capacity can be assumed to be 100, and the speed of the vehicle can be assumed to be 1 (Table 1).
According to the algorithm designed above, we first find that the TSP path is v 1 v 7 v 6 v 9 v 8 v 11 v 2 v 10 v 4 v 5 v 3 . Second, according to steps 2–4 of the algorithm, the distribution path of each vehicle, the quantity of goods to be collected and the travel time are obtained, as shown in Table 2 below.

5.2. Sensitivity Analysis

In reality, logistics and distribution companies often have a variety of vehicles, that is, vehicles with different loading capacities, from which to choose. Here, we tested whether a change in the vehicle capacity will have an impact on the results. Continuing to use the above modified Solomon VRP test set R101 example, the vehicle capacity is set to 110, 120, and 150, respectively; and the solution results are shown in Table 3, Table 4, Table 5 respectively.
The comparison shows that as the vehicle capacity increases, the overall travel time shows a downward trend. The reason is that as the vehicle capacity increases, the constraints loosen, and more distribution routes become feasible paths from which to choose, so the overall travel time is reduced. Generally, it is better to choose a vehicle with a larger capacity than one with a smaller capacity for distribution. However, vehicles with large capacities often correspond to higher transport costs per unit, so they cannot be generalized in reality. Logistics and distribution companies need to weigh and select appropriate vehicles according to the actual situation.

5.3. Extended Testing Example

Table 1 shows an example with 10 demand points. In this part, the range of the demand points was extended to 20. The data of nodes 12–21 are shown in the following 6 (Table 6). The data of nodes 1–11 and other assumptions are the same as the 5.1 testing example.
According to the algorithm designed above, the distribution path of each vehicle, the quantity of goods to be collected and the travel time are obtained, as shown in Table 7 below.
It is shown that the algorithm is also applicable to the solution of more demand points.

6. Conclusions

In this paper, a bi-objective vehicle routing model with simultaneous pickup and delivery minimizing the total travel time and maximizing the total number of collected goods simultaneously was developed. An approximate algorithm based on the ε -constraint method was designed. The time complexity analysis showed that the algorithm was a polynomial time algorithm, and the approximation ratio of the algorithm was analyzed.
The model and algorithm were tested by two examples, and the influence of vehicles with different capacities on the travel time is tested by sensitivity analysis. The results showed that as the vehicle capacity increases, the overall travel time showed a downward trend, and a vehicle with a larger loading capacity was better than a vehicle with a smaller loading capacity.
In the future, the research of this paper can still be strengthened through the following work. First, the pickup quantity in this paper is known. In practice, sometimes a vehicle only knows the quantity of delivery goods but does not know the exact pickup quantity before departure; therefore, the pickup quantity can be set to be uncertain in the future to be closer to the actual situation. Second, we can further design an algorithm with a faster solution speed or an approximate algorithm that has a better approximation ratio.

Author Contributions

Conceptualization, N.W.; Writing—original draft, Q.G.; Writing—review & editing, Q.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Major projects of National Natural Science Foundation of China, grant number 72192830 & 72192834, the Key Project of National Natural Science Foundation of China, grant number 71732006, Xi’an Science and Technology Plan Soft Science Project, grant number 22RKYJ0069.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Min, H. The multiple vehicle routing problem with simultaneous delivery and pick-up points. Transp. Res. 1989, 23, 377–386. [Google Scholar] [CrossRef]
  2. Awd, H.; Elshaer, R. A Taxonomic Review of Metaheuristic Algorithms for Solving the Vehicle Routing Problem and Its Variants. Comput. Ind. Eng. 2019, 140, 106242. [Google Scholar]
  3. Dell’Amico, M.; Righini, G.; Salani, M. A branch-and-price approach to the vehicle routing problem with simultaneous distribution and collection. Transp. Sci. 2006, 40, 237–247. [Google Scholar] [CrossRef] [Green Version]
  4. Rieck, J.; Zimmermann, J. A branch-and-cut approach to the vehicle routing problem with simultaneous delivery and pick-up. In Operations Research Proceedings 2008; Fleischmann, B., Borgwardt, K.H., Klein, R., Tuma, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 301–306. [Google Scholar]
  5. Subramanian, A.; Uchoa, E.; Ochi, L.S. New lower bounds for the vehicle routing problem with simultaneous pickup and delivery. In Experimental Algorithms. SEA 2010. Lecture Notes in Computer Science; Festa, P., Ed.; Springer: Berlin/Heidelberg, Germany, 2010; Volume 6049. [Google Scholar]
  6. Subramanian, A.; Uchoa, E.; Pessoa, A.A.; Ochi, L.S. Branch-and-cut with lazy separation for the vehicle routing problem with simultaneous pickup and delivery. Oper. Res. Lett. 2011, 39, 338–341. [Google Scholar] [CrossRef]
  7. Archetti, C.; Speranza, M.G.; Boccia, M.; Sforza, A.; Sterle, C. A branch-and-cut algorithm for the inventory routing problem with pickups and deliveries. Eur. J. Oper. Res. 2020, 282, 886–895. [Google Scholar] [CrossRef]
  8. Subramanian, A.; Uchoa, E.; Pessoa, A.A.; Ochi, L.S. Branch-cut-and-price for the vehicle routing problem with simultaneous pickup and delivery. Optim. Lett. 2013, 7, 1569–1581. [Google Scholar] [CrossRef]
  9. Li, C.; Gong, L.; Luo, Z.; Lim, A. A branch-and-price- and-cut algorithm for a pickup and delivery problem in retailing. Omega 2019, 89, 71–91. [Google Scholar] [CrossRef]
  10. Dethloff, J. Vehicle routing and reverse logistics: The vehicle routing problem with simultaneous delivery and pick-up. OR Spektrum 2001, 23, 79–96. [Google Scholar] [CrossRef]
  11. Kalayci, C.B.; Kaya, C. An ant colony system empowered variable neighbourhood search algorithm for the vehicle routing problem with simultaneous pickup and delivery. Expert Syst. Appl. 2016, 66, 163–175. [Google Scholar] [CrossRef]
  12. Zachariadis, E.E.; Tarantilis, C.D.; Kiranoudis, C.T. Vehicle routing strategies for pick-up and delivery service under two dimensional loading constraints. Oper. Res. 2017, 17, 115–143. [Google Scholar] [CrossRef]
  13. Simsir, F.; Ekmekci, D. A metaheuristic solution approach to capacitied vehicle routing and network optimization. Eng. Sci. Technol. Int. J. 2019, 22, 727–735. [Google Scholar] [CrossRef]
  14. Li, J.; Pardalos, P.; Sun, H.; Pei, J.; Zhang, Y. Iterated local search embedded adaptive neighbourhood selection approach for the multi-depot vehicle routing problem with simultaneous deliveries and pickups. Expert Syst. Appl. 2015, 42, 3551–3561. [Google Scholar] [CrossRef]
  15. Hs, A.; Yc, A.; Hg, B. Collection and distribution of returned remanufactured products in a vehicle routing problem with pickup and delivery considering sustainable and green criteria. J. Clean. Prod. 2018, 172, 960–970. [Google Scholar]
  16. Avci, M.; Topaloglu, S. A hybrid metaheuristic algorithm for heterogeneous vehicle routing problem with simultaneous pickup and delivery. Expert Syst. Appl. 2016, 53, 160–171. [Google Scholar] [CrossRef]
  17. Zhang, Z.; Cheang, B.; Li, C.; Lim, A. Multi-commodity demand fulfillment via simultaneous pickup and delivery for a fast fashion retailer. Comput. Oper. Res. 2019, 103, 81–96. [Google Scholar] [CrossRef]
  18. Wang, Y.; Ran, L.; Guan, X.; Zou, Y. Multi-Depot pickup and delivery problem with resource sharing. J. Adv. Transp. 2021, 2021, 5182989. [Google Scholar] [CrossRef]
  19. Wang, C.; Mu, D.; Zhao, F.; Sutherland, J. A parallel simulated annealing method for the vehicle routing problem with simultaneous pickup-delivery and time windows. Comput. Ind. Eng. 2015, 83, 111–122. [Google Scholar] [CrossRef]
  20. Carolina, L.; Guillermo, G.; Enrique, C.; Moltedo, A.; Johnson, F.; Paredes, F. An improved particle swarm optimization algorithm for the VRP with simultaneous pickup and delivery and time windows. Lat. Am. Trans. 2018, 16, 1732–1740. [Google Scholar]
  21. Hornstra, R.P.; Silva, A.; Roodbergen, K.J.; Coelho, L.C. The vehicle routing problem with simultaneous pickup and delivery and handing costs. Comput. Oper. Res. 2019, 115, 104858. [Google Scholar] [CrossRef]
  22. Liu, Q.; Liu, W.J. Study on modeling and optimization of SVRPSDP based on ACS-RSM. In Proceedings of the IEEE Conference on Grey Systems and Intelligent Services, Nanjing, China, 15–18 September 2011; pp. 791–794. [Google Scholar]
  23. Peng, X. Research on Reverse Logistics Vehicle Routing Problem under Stochastic Recovery Demand with Time Windows. Master’s Thesis, Wuhan University of Technology, Wuhan, China, 2013. [Google Scholar]
  24. Jia, Y.J.; Gu, H.Y.; Xi, Y.G. Rolling horizon scheduling algorithm for dynamic vehicle scheduling system. J. Southeast Univ. 2005, 21, 92–96. [Google Scholar]
  25. Hu, Z.-H.; Sheu, J.-B.; Zhao, L.; Lu, C.-C. A dynamic closed-loop vehicle routing problem with uncertainty and incompatible goods. Transp. Res. Part C 2015, 55, 273–297. [Google Scholar] [CrossRef]
  26. Fan, H.; Liu, P.; Wu, J.; Lu, C.C. Hybrid genetic algorithm with variable neighborhood descent for the vehicle routing problem with simultaneous stochastic pickup and deterministic delivery. Syst. Eng.-Theory Pract. 2019, 39, 2646–2659. [Google Scholar]
  27. Fan, H.; Liu, P.; Liu, T. The multi-depot vehicle routing problem with simultaneous deterministic delivery and stochastic pickup based on joint distribution. ACTA Autom. Sin. 2021, 47, 1646–1660. [Google Scholar]
  28. Minis, I.; Tatarakis, A. Stochastic single vehicle routing problem with delivery and pick up and a defined customer sequence. Eur. J. Oper. Res. 2011, 213, 37–51. [Google Scholar] [CrossRef]
  29. Zhu, L.; Sheu, J.B. Failure-specific cooperative recourse strategy for simultaneous pickup and delivery problem with stochastic demands. Eur. J. Oper. Res. 2018, 271, 896–912. [Google Scholar] [CrossRef]
  30. Zhou, R. Research on Optimization Models and Algorithms for Vehicle Routing Problem with Simultaneous Pickup and Delivery. Ph.D. Thesis, Hefei University of Technology, Hefei, China, 2016. [Google Scholar]
  31. Ma, Y.; Li, Z.; Yan, F.; Feng, C. A hybrid priority-based genetic algorithm for simultaneous pickup and delivery problems in reverse logistics with time windows and multiple decision-makers. Soft Comput. 2019, 23, 6697–6714. [Google Scholar] [CrossRef]
  32. Fan, H.; Liu, H.; Liu, P.; Ren, X. Heterogeneous fleet vehicle routing problem with simultaneous deterministic delivery and fuzzy pickup. Control Theory Appl. 2021, 38, 15. [Google Scholar]
  33. Bérubé, J.; Gendreau, M.; Potvin, J. An exact ε-constraint method for bi-objective combinatorial optimization problems: Application to the traveling salesman problem with profits. Eur. J. Oper. Res. 2009, 194, 39–50. [Google Scholar] [CrossRef]
  34. Christofides, N. Worst Case Analysis of a New Heuristic for the Traveling Salesman Problem; Report 388 Graduate School of Industrial Administration; Carnegie Mellon University: Pittsburgh, PA, USA, 1976. [Google Scholar]
  35. Bompadre, A.; Dror, M.; Orlin, J. Improved bounds for vehicle routing solutions. Discret. Optim. 2006, 3, 299–316. [Google Scholar] [CrossRef]
Table 1. Testing data (from modified Solomon VRP test set R101).
Table 1. Testing data (from modified Solomon VRP test set R101).
Node Number1234567891011
X Axis3541355555152520105530
Y Axis3549174520303050436060
Number of Pickup Goods032171319255035101426
Number of Delivery Goods0411231673629403523
Table 2. Results of path optimization when the vehicle capacity is 100.
Table 2. Results of path optimization when the vehicle capacity is 100.
VehicleRouteThe Number of Collected GoodsThe Travel Time
Vehicle 1 v 1 v 7 v 6 v 1 7542
Vehicle 2 v 1 v 9 v 8 v 11 v 1 7178
Vehicle 3 v 1 v 2 v 10 v 1 4665
Vehicle 4 v 1 v 4 v 5 v 3 v 1 4986
Total 241271
Table 3. Results of path optimization when the vehicle capacity is 110.
Table 3. Results of path optimization when the vehicle capacity is 110.
VehicleRouteThe Number of Collected GoodsThe Travel Time
Vehicle 1 v 1 v 7 v 6 v 1 7542
Vehicle 2 v 1 v 9 v 8 v 11 v 1 7178
Vehicle 3 v 1 v 2 v 10 v 4 v 1 5970
Vehicle 4 v 1 v 5 v 3 v 1 3663
Total 241253
Table 4. Results of path optimization when the vehicle capacity is 120.
Table 4. Results of path optimization when the vehicle capacity is 120.
VehicleRouteThe Number of Collected GoodsThe Travel Time
Vehicle 1 v 1 v 7 v 6 v 9 v 1 8561
Vehicle 2 v 1 v 8 v 11 v 2 v 1 9366
Vehicle 3 v 1 v 10 v 4 v 5 v 3 v 1 63110
Total 241237
Table 5. Results of path optimization when the vehicle capacity is 150.
Table 5. Results of path optimization when the vehicle capacity is 150.
VehicleRouteThe Number of Collected GoodsThe Travel Time
Vehicle 1 v 1 v 7 v 6 v 9 v 8 v 1 12069
Vehicle 2 v 1 v 11 v 2 v 10 v 4 v 5 v 1 104124
Vehicle 3 v 1 v 3 v 1 1736
Total 241229
Table 6. Extended testing data (from modified Solomon VRP test set R101).
Table 6. Extended testing data (from modified Solomon VRP test set R101).
Node Number12131415161718192021
X Axis2050301530105201545
Y Axis6535251052030406065
Number of Pickup Goods2381021282117243048
Number of Delivery Goods27112414355023183432
Table 7. Results of path optimization when the vehicle capacity is 100.
Table 7. Results of path optimization when the vehicle capacity is 100.
VehicleRouteThe Number of Collected GoodsThe Travel Time
Vehicle 1 v 1 v 14 v 7 v 1 6037
Vehicle 2 v 1 v 6 v 18 v 17 v 1 6371
Vehicle 3 v 1 v 15 v 16 v 3 v 5 v 13 v 1 93112
Vehicle 4 v 1 v 4 v 10 v 21 v 1 7580
Vehicle 5 v 1 v 2 v 11 v 12 v 1 8176
Vehicle 6 v 1 v 20 v 8 v 1 6564
Vehicle 7 v 1 v 9 v 19 v 1 3453
Total 471493
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Q.; Wang, N. The Vehicle Routing Problem with Simultaneous Pickup and Delivery Considering the Total Number of Collected Goods. Mathematics 2023, 11, 311. https://doi.org/10.3390/math11020311

AMA Style

Guo Q, Wang N. The Vehicle Routing Problem with Simultaneous Pickup and Delivery Considering the Total Number of Collected Goods. Mathematics. 2023; 11(2):311. https://doi.org/10.3390/math11020311

Chicago/Turabian Style

Guo, Qinge, and Nengmin Wang. 2023. "The Vehicle Routing Problem with Simultaneous Pickup and Delivery Considering the Total Number of Collected Goods" Mathematics 11, no. 2: 311. https://doi.org/10.3390/math11020311

APA Style

Guo, Q., & Wang, N. (2023). The Vehicle Routing Problem with Simultaneous Pickup and Delivery Considering the Total Number of Collected Goods. Mathematics, 11(2), 311. https://doi.org/10.3390/math11020311

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