Multi-Depot Pickup and Delivery Problem with Resource Sharing

Resource sharing (RS) integrated into the optimization of multi-depot pickup and delivery problem (MDPDP) can greatly reduce the logistics operating cost and required transportation resources by reconfiguring the logistics network. This study formulates and solves an MDPDP with RS (MDPDPRS). First, a bi-objective mathematical programming model that minimizes the logistics cost and the number of vehicles is constructed, in which vehicles are allowed to be used multiple times by one or multiple logistics facilities. Second, a two-stage hybrid algorithm composed of a k-means clustering algorithm, a Clark-Wright (CW) algorithm, and a nondominated sorting genetic algorithm II (NSGA-II) is designed. The k-means algorithm is adopted in the first stage to reallocate customers to logistics facilities according to the Manhattan distance between them, by which the computational complexity of solving the MDPDPRS is reduced. In the second stage, CW and NSGA-II are adopted jointly to optimize the vehicle routes and find the Pareto optimal solutions. CW algorithm is used to select the initial solution, which can increase the speed of finding the optimal solution during NSGA-II. Fast nondominated sorting operator and elite strategy selection operator are utilized to maintain the diversity of solutions in NSGA-II. Third, benchmark tests are conducted to verify the performance and effectiveness of the proposed two-stage hybrid algorithm, and numerical results prove that the proposed methodology outperforms the standard NSGA-II and multi-objective particle swarm optimization algorithm. Finally, optimization results of a real-world logistics network from Chongqing confirm the applicability of the mathematical model and the designed solution algorithm. Solving the MDPDPRS provides a management tool for logistics enterprises to improve resource configuration and optimize logistics operation efficiency.


Introduction
With the advancement of information technology and Internet of ings, the logistics industry is playing an increasingly important role in the development of modern businesses [1,2]. However, national and local governments worldwide are focusing on the environmental impacts of logistics and the efficient use of resources [3,4]. In a logistics network, customers send out a series of requests for delivery and pickup services, and logistics service providers (LSPs) design service plans and arrange vehicles for these requests to deliver or pickup goods [5,6]. Efficient logistics service plan can improve the operation efficiency of LSPs and resource utilization [7,8]. erefore, making an effective logistics service plan with resource sharing (RS) is essential, which not only helps to reduce the operating cost for logistics facilities but also promotes the development of green logistics and provides better logistics services for consumers [9,10].
In this study, a multi-depot pickup and delivery problem with RS (MDPDPRS) combines components from three subproblems: multi-depot vehicle routing problem (MDVRP) with pickups and deliveries (MDVRPPD), MDVRP with pickups and deliveries and time windows (MDVRPPDTW), and RS [11][12][13]. However, one of the difficult challenges in solving this problem is how to handle pickup and delivery activities among multiple depots through RS [14,15]. In the traditional MDVRPPD, each vehicle performs only one type of activity in the service route, which may be delivering or picking up goods [16][17][18].
In addition, the traditional MDVRPPDTW and single-depot vehicle routing problem with pickups and deliveries mostly only consider the optimization of logistics operational costs [19,20]. erefore, MDPDPRS focuses on how to support and achieve the efficient utilization of transportation resources with RS strategy, and optimizes the logistics network.
With regard to RS, it is often jointly adopted with collaboration or cooperation between LSPs to optimize the logistics networks with multiple depots [21,22]. RS strategy supports the sharing of customer information and transport resources to improve the resource configuration among logistics facilities to optimize the logistics network [23,24]. Here, the sharing of customer information is often enabled by customer clustering, whereas the sharing of transportation resources is related to the use of shared transportation equipment [25,26].
As for MDVRPPD, it is a crucial logistics issue with extensive applications, especially in reverse logistics [27]. ree basic types of vehicle routing problems exist in reverse logistics [28]. e first type is the vehicle routing problem with mixed deliveries and pickups, which involves customers with delivery demands, pickup demands, and delivery and pickup demands [29]. Vehicle routing problem with simultaneous delivery and pickup is the second type, which requires all customers to have both delivery and pickup demands [30]. e third type is the common MDVRPPD, which includes delivery and pickup customers in the logistics network [19,[31][32][33]. In this study, the consideration of customer service time windows makes MDVRPPD realistic.
In this study, the MDPDPRS can be formulated into a biobjective mathematical model to minimize the total logistics operating cost and number of vehicles [32,34]. On the basis of the multi-depot and RS properties of MDPDPRS, a two-stage hybrid algorithm is proposed to find the Pareto optimal solution. In the first stage, a k-means algorithm is adopted to reconfigure resources through customer clustering; thus, the MDVRPPDTW is simplified for solving [35]. e second stage focuses on finding the Pareto optimal solution for the bi-objective optimization problem [23,36].
e Clarke-Wright (CW) algorithm, which is good at constructing the initial solution of vehicle routes, and the nondominated sorting genetic algorithm (NSGA-II), which is known for its capability of finding the Pareto solution, are adopted to optimize the vehicle routing in the second stage [33,37,38]. e remainder of this study is arranged as follows. Section 2 reviews the relevant literature. Section 3 elaborates the specifics of the MDPDPRS. Section 4 explains the bi-objective mathematical model for the MDPDPRS. Section 5 presents the designed methodology for solving the MDPDPRS. Section 6 analyzes the performance and application of the proposed model formulation and solution algorithm in a real-world case study compatible to the MDPDPRS. Finally, Section 7 summarizes the conclusions and discusses potential future research.

Literature Review
MDPDPRS is mainly related to the MDVRPPD, MDVRPPDTW, and RS strategy [18,36]. MDVRPPD and MDVRPPDTW are the extension problems of MDVRP and MDVRPTW with respect to the logistics service type of customer demands, responsively [39,40]. In the widely studied MDVRP and MDVRPTW, the service types of customers are either deliveries or pickups [16,33,41]. However, in a real-world logistics network with multiple depots, customers with distribution and pickup demands often exist simultaneously and the service time windows are the additional characteristics of customers; this issue is abbreviated by scholars as MDVRPPD [11,13]. erefore, MDVRPPD and MDVRPPDTW have begun to attract the attention of scholars, and the difference between the two issues is mainly whether customers' service time window feature is considered [42,43]. In the MDVRPPDTW, the optimization of vehicle routes focuses on the confirmation of customers' service time windows and the integration of vehicles' delivery and pickup activities, which are suitable with the factors that LSPs must consider if optimizing vehicle routes [44,45]. e adoption of RS into the optimization of logistics network is also a current trend [12]. us, the MDPDPRS in this study integrates vehicles' distribution and pickup arrangement with the RS strategy to optimize the MDVRPPDW.
In contrast to the MDVRP, the MDVRPPD studies customers with delivery and pickup demands [35,42]. e MDVRPPDTW is an extension of MDVRPPD, which considers the characteristics of customer service time windows [46,47]. Many scholars have studied MDVRPPD and MDVRPPDTW considering diverse aspects and proposed different mathematical models and algorithms [30,48]. In terms of models, many of the proposed mathematical models reflect the characteristics of their problem studied by different constraints, including capacity, time windows, and priority constraints [29,49]. Ropke et al. [50] established a standard three-index model based on the characteristics of customer time window and designed an accurate algorithm to solve it. Gribkovskaia et al. [19] studied the vehicle routing problem with deliveries and pickups considering the number of times a customer has been visited and proposed a mixed integer linear programming model. Chen et al. [31] established a comprehensive mathematical model to minimize the transportation costs for unpaired vehicle routing problem with deliveries and pickups in a multi-factory production network. Conversely, hybrid heuristics algorithms (e.g., genetic algorithm and adaptive large neighborhood search algorithm) and exact algorithms based on column generation are commonly designed to solve MDVRPPD and MDVRPPDTW [50][51][52].
e proposed model and solution methodology in the above studies provide abundant reference for solving the basic MDVRPPD and MDVRPPDTW. However, few studies have optimized MDVRPPDTW with RS strategy [32,34].
Customer clustering analysis is a research aspect that groups customers based on their characteristics (e.g., location and time window), and common clustering algorithms include k-means clustering, parallel clustering, and fuzzybased customer clustering [53][54][55][56]. In comparison with other customer clustering algorithms, k-means clustering is widely adopted to solve vehicle routing problems [57,58].
Xu et al. [59] proposed an enhanced ant colony algorithm based on k-means clustering to solve dynamic vehicle routing problems and achieved good optimization results. Hakim et al. [60] designed a cluster-based method to solve a vehicle routing problem with limited vehicle capacity, and their calculation results proved the effectiveness of that method. Mourelo Ferrandez et al. [61] reduced the calculation difficulty of truck-drone in tandem delivery network by k-means clustering algorithm. Wang et al. [35] proposed a hybrid heuristic algorithm based on three-dimensional kmeans clustering and improved reference point NSGA-II to solve the multi-objective optimization model. erefore, cluster analysis can simplify the difficulty of solving vehicle routing problem [62,63].
At present, the construction of multi-objective optimization model and multi-objective optimization algorithm is the research hotspot of finding the Pareto solution of vehicle routing problems, and scholars have designed different algorithms [38,64]. NSGA-II and multi-objective particle swarm optimization (MOPSO) algorithm are two common multi-objective algorithms [65,66]. NSGA-II adopts a reference point strategy to maintain population diversity [67]. Srivastava et al. [68] proposed a NSGA-II to solve the multiobjective optimization model for MDVRPTW and verified that the method is superior to the latest method of that problem through a real-world case study. Maadanpour Safari et al. [69] adopted NSGA-II, multi-objective simulated annealing (MOSA), and MOPSO to optimize the proposed three-objective mathematical function and concluded that NSGA-II was superior to MOSA according to the results of their examples. Shafiei Nikabadi et al. [70] formed a multiobjective model for the route selection of freight fleet, optimized it with NSGA-II and MOPSO, and considered that MOPSO was superior to NSGA-II. erefore, NSGA-II and MOPSO are two commonly typical multi-objective optimization algorithms that can find Pareto solutions [71].
Many scholars have adopted the cooperation and RS strategy to optimize the multicenter logistics network [11,17,24]. Zhang et al. [15] believed that in the collaborative e-commerce truck carrier, participants can share transportation resources and customer demands to maximize the total profit of the whole alliance and improve vehicle utilization. Wang et al. [14] adopted the cooperation strategy to optimize the MDVRPPDTW with minimization of the operating cost of the transportation network and the total number of vehicles. Deng et al. [21] allowed the capabilities of the logistics facility, the vehicle resources, and the customer information to be shared through RS strategies, and they proved that this strategy can improve the utilization of logistics resources. In the study of Nourinejad et al. [72], vehicles can be used multiple times to reduce fleet size by extending vehicle reservation time. Li et al. [73] adopted the resource-sharing strategy to significantly optimize the logistics network and maximized the utilization of resources. erefore, the adoption of RS strategy not only helps optimize the logistics costs but also improves the utilization of resources to protect the environment [22,74].
In summary, the existing literature has provided rich reference materials about MDVRPPD, MDVRPPDTW, and RS, including model formulations and solution algorithms. However, the existing literature related to MDPDPRS has the following limitations. (1) Few studies on MDPDPRS have considered RS, MDVRPPD, and MDVRPPDTW. (2) e fact that a vehicle can be used multiple times on a working day is insufficiently considered in the proposed mathematical models. (3) Most of the proposed solution algorithms in the existing literature only address how to solve MDVRP, MDVRPPD, and MDVRPPDTW, and RS has not been incorporated into the designed algorithms. (4) Most of the existing literature focuses on raising problems and designing algorithms but neglects testing the proposed methods with practical cases.
In consideration of the aforementioned shortcomings, the main contributions of this study to MDPDPRS are as follows: (1) Characteristics of RS, MDVRPPD, and MDVRPPDTW are comprehensively incorporated to enrich the research on MDPDPRS. (2) On the basis of RS, this study proposes and tests that vehicles can be used multiple times and that customer information can be shared to save the transportation resources of logistics networks, which are considered in the proposed mathematical model. (3) A twostage algorithm is designed to combine RS with vehicle routing optimization to optimize the logistics network. (4) Benchmark and real-world cases are utilized to test and verify the performance and applicability of the proposed model and solution algorithm in this study.

Problem Statement
RS is an effective strategy that can optimize logistics operation costs and resource allocation in a multicenter logistics network with pickups and deliveries [21,35]. In this study, the logistics network consists of multiple distribution centers (DCs), multiple pickup centers (PCs), and multiple customers. e logistics network before and after optimization with RS, which is composed of DC1, DC2, DC3, PC1, PC2, and 52 customers (marked C1, C2, . . ., C52), is shown in Figure 1. e numerical number near the line represents the time distance between two elements (including facilities and customers).
In Figure 1(a), the unreasonable arrangement of vehicle routes and the nonsharing of resources are the main reasons for the higher operating costs of the logistics network. First, staggered driving and long-distance service are the two most significant unreasonable arrangements, and they cause additional travel costs. Second, vehicles that violate customers' time windows often occur. Arriving early and arriving late are generating penalty costs. Finally, the nonsharing of resources between facilities results in the capacity of facilities and transportation resources being left unused. For example, the vehicle no longer works after returning from C2 to DC1 in the service route DC1 ⟶ C25 ⟶ C6 ⟶ C1 ⟶ C2 ⟶ DC1.
To clearly demonstrate the discounts that RS brings to logistics network optimization, six indicators are counted, namely, total travel cost (TTC), total penalty cost (TPC), total maintenance cost (TMC), total fixed cost (TFC), number of vehicles (NV), and total operating cost (TOC), which are listed in Table 1. e unit time travel cost using a vehicle can be set to $10/h, and the unit time penalty cost for waiting or late can be set to $20/h, and the maintenance cost using a vehicle in a working period can be defined as $100, and the fixed cost of a facility being used in a working period can be set to $200.
In Table 1, the gap of TOC and NV in the initial and optimized logistics network is $1880 and 7, respectively, which are highly significant. Before optimization, TTC ($1500), TPC ($440), and TMC ($1400) are the main elements that cause additional logistics costs. In the optimized logistics network, the penalty costs are avoided, and the reduction in the number of vehicles also reduces the maintenance costs with RS. In addition, the travel costs of vehicles are also greatly reduced. erefore, RS is an effective strategy for optimizing the logistics network of the MDPDPRS.

Assumption 3.
e vehicles can be used multiple times by the same or different logistics facilities.
e relevant notations and their descriptions are listed in Table 2.

Mathematical Model.
A bi-objective mathematical model that minimizes the logistics operational costs (Equation (1)) and the number of vehicles (Equation (2)) is formed to solve the MDPDPRS. In Equation (1), the logistics operation cost is composed of four parts, which are marked by the TTC (Equation (3)), TPC (Equation (4)), TMC (Equation (5)), and TFC (Equation (6)). In Equation (2), c∈C P v∈V c∈C P k∈R v Journal of Advanced Transportation Constraint (7) ensures that each customer is served once. Constraint (8) counts the shared times of vehicle v. Constraints (9)-(11) ensure that flow conservation is achieved on each customer. Constraints (12) and (13) are the flow balance constraints of the truck. Constraints (14)- (16) ensure that the loading quantity of each vehicle and each truck cannot Fuel consumption rate of truck b per km (unit: gallon/miles) P v Gasoline price (unit: dollars/gallon) P b Gasoline price (unit: dollars/gallon) x cod � 1 if the facility providing logistics service for customer c is changed from facility o to f; otherwise, x bowf � 1, if truck k departs from facility o and travels directly from facility w to f; otherwise, exceed their capabilities. Constraints (17) and (18) count the quantity of transshipment goods between logistics facilities. Constraints (19) and (20) guarantee that the total service quantity of each facility does not exceed its capacity. Constraints (21) and (22) require that the departure time and return time of each vehicle must meet the service time window of its served facility. Constraints (23) and (24) ensure that each vehicle must provide services for customers within the customers' service time window. Constraint (25) requires that the total working time of each vehicle does not exceed its maximum working time. Constraints (26)- (28) are used to eliminate the sub-tours of each vehicle and truck. e constraints of relevant binary variables are listed in Constraints (29)-(34).

Solution Methodology for MDPDPRS
MDVRPTW and MDVRPPDTW are typical NP-hard problems [12,18,41]. Multi-objective optimization algorithm and two-stage algorithm are often designed in combination to solve MDVRPPDTW [14,23,36]. Here, a two-stage algorithm with customer clustering first and then vehicle routing optimization is designed to solve MDPDPRS. is two-stage hybrid algorithm is composed of k-means, CW, and NSGA-II algorithms, and named KCW-NSGA-II. In the first stage of KCW-NSGA-II, customers and resources are reconstructed by the k-means clustering algorithm [57,58]. e main purpose of the second stage is to optimize vehicle routes and find the Pareto optimal solution. e CW algorithm is adopted to construct the initial solution for NSGA-II [33,37,64,65,67]. e designed algorithm flow is illustrated in Figure 2. Here, Gen is the current number of iterations; MaxGen is the maximum number of iterations; r is the number of iterations of the current internal re-optimization mechanism, which is between clustering and vehicle routing optimization; and MaxR is its maximum number of iterations.
In Figure 2, the two-stage characteristics of customer clustering first and vehicle routing optimization later are clearly demonstrated. First, a k-means customer clustering mechanism based on Manhattan distance is designed. e clustering results are checked, updated, and saved after finishing the reallocation of all customers. Second, the CW algorithm is adopted to design the initial population and initial feasible solution, which accelerates the speed and possibility of NSGA-II algorithm to find the Pareto optimal solution. ird, the elite strategy and genetic operation of NSGA-II are used to iteratively optimize the generated initial solution to find the Pareto optimal solution, which is mainly embodied in the change of Gen. Fourth, a regulatory re-optimization mechanism between customer clustering and vehicle routing optimization is set up to maintain gene stability during genetic operation, and this mechanism is implemented by the re-updating of r. Finally, if Gen is updated to MaxGen, then the iterative optimization of the algorithm is finished and the found Pareto optimal solution is outputted.

K-Means Clustering Algorithm.
Customer clustering is an important measure to reduce the complexity of solving MDVRPPDTW [36]. K-means algorithm is widely used to solve MDVRPTW due to its simplicity and efficiency [14,21]. e k-means clustering pseudocode based on Manhattan distance is listed in Algorithm 1.

CW Algorithm.
A common and effective way to construct the initial solution of VRP is the CW savings algorithm, which is actually a greedy heuristic algorithm [33,37,76]. e service time windows of customers and the capacity of vehicles are the main constraints to construct the initial solution [77,78]. e pseudocode of the CW algorithm designed in this study is listed in Algorithm 2.

NSGA-II.
NSGA-II is a multi-objective optimization algorithm based on GA, which searches the Pareto optimal solution for multi-objective optimization [38,64,68]. Fast nondominated sorting operator, individual crowding distance operator, and elite strategy selection operator are the three key designs of NSGA-II [33,65]. Here, we suppose that the population is P and n individuals exist, and the individual objective function value of individual i is x i .

Fast Nondominated Sorting Operator.
e key design of NSGA-II is to find the Pareto optimal solution. To enhance the possibility of finding the Pareto solution, the fast nondominated sorting operator stratifies the population P according to the quality of individual solutions [33,64,65,67]. is method is a cyclic process of grading based on population fitness. Here, the nondominated solution set and the rank value assigned to the individual are the two indicators for fast nondominated sorting. Assume that F i represents the nondominated solution set, and i r represents the rank value of individual i. en, the pseudocode of the fast nondominated sort operator is shown in Algorithm 3.

Crowding Distance and Its Comparison Principle.
e crowding distance is designed to drive the population to converge to the Pareto optimal solution and maintain the diversity of the population, which is mainly for individuals in the same nondominant layer [38,64,68]. We assume that n individuals exist in the nondominant level of S F , and the objective function value of individual i is x i . en, the crowding distance is calculated as In Equation (35), |x i− 1 − x i+1 | represents the sum of the distance between individuals i − 1 and i + 1 in each direction of the objective function. Here, the objective function values of individuals need to be sorted before calculating the crowding distance. If the rank values of individuals i and j are i r and j r , then the crowding distance is L(i) and L(j), respectively. e individual crowding comparison strategy based on crowding distance and nondominated ranking results is as follows.

Journal of Advanced Transportation
Input: e datasets, including logistics facility and customer information, such as the coordination, time windows, and demands Output: e clustering results (1) Step 1: Select k objects as the initial clustering center (2) Step 2: Calculate the Manhattan distance between each customer and each clustering center (3) Step 3: (Re-)Assign each customer to their closest clustering center (4) Step 4: If some customers need to be adjusted among the clustering results, then enter Step 3; otherwise, go to Step 5 (5) Step 5: Update the clustering centers (6) Step 6: Output the clustering results ALGORITHM 1: Procedure of k-means algorithm.  (1) If i r > j r , then individual i is the best one. If i r < j r , then individual j is the best one.
(2) If i r � j r , then the individual with the most crowded distance is the better one.

Elite Strategy Selection Operator.
To prevent the Pareto optimal solution from being lost in the iteration process, an elite strategy selection operator is designed, which selects the optimal solution by nondominated sorting and crowding distance between the parent and offspring populations. Suppose that the current iteration is t and the parent population is P t , the offspring population is Q t , and R t is composed by P t and Q t . First, a fast nondominated sorting is performed for R t , and the crowding distance is then calculated. On the basis of the crowding distance and the nondominated layer, N individuals with high-quality solution is selected to form a new population P t+1 .

Algorithm Comparison.
e standard NSGA-II, GA-PSO, and MOPSO are adopted for comparison to verify the applicability and effectiveness of the proposed KCW-NSGA-II algorithm in solving MDPDPRS [38, 64-66, 68, 79]. e benchmark dataset C-mdvrptw (consisting of instances of 20 groups) is utilized for the test, which is mainly obtained from the website. Networking and emerging optimization and their related characteristics are listed in Table 3. To meet the characteristics of the research object in this study, depots are regarded as logistics facilities, and customers are divided into two types, that is, those with distribution demands and those with pickup demands.
In Table 3, the number of customers and logistics facilities in each instance is different. e first instance comprises four facilities and 48 customers, whereas Instance 10 includes six facilities and 288 customers. In addition, the loading capacity of vehicle used in each instance is differentiated.
Relevant parameters are properly unified to mitigate their effects on algorithm performance. ese parameters are set as follows [38, 64-66, 68, 79] Table 4. e numerical results shown in Table 4 demonstrate that the proposed algorithm KCW-NSGA-II performs better than the other three algorithms. First, the average cost of the four algorithms is $2584, $2944, $2746, and $3025, respectively. By contrast, the costs and the number of vehicles optimized by KCW-NSGA-II are the most economical solution compared with the other three algorithms in each instance. Second, the value of t-test also shows that the KCW-NSGA-II is significantly different from the other three algorithms. In addition, the proposed algorithm can obtain the optimal solution quickly. erefore, the proposed algorithm KCW-NSGA-II outperforms the other three algorithms. Moreover, this algorithm can be adjusted to address problems such as VRPMDP, VRPSDP, and PVRP.

Data Source and Relevant Parameter Setting. As an inland international logistics hub and an open highland,
Chongqing is a new first-tier city in China. erefore, as our numerical experiments, the logistics network adopted from Chongqing is appropriate to verify the applicability of this study. Six logistics facilities (i.e., DC1, DC2, DC3, PC1, PC2, and PC3) and 220 customers are the main elements of this real-world logistics network. e information and characteristics of these elements are listed in Table 5, and the spatial allocation information is plotted in Figure 3.
In Table 5, the number of customers served by the six logistics facilities is 27, 36, 39, 46, 31, and 41, respectively. In Figure 3, an obvious feature is that the customer allocation of each facility is relatively dispersed. e service area edge of each facility is not a clear division. In Table 6, the initial vehicle routes for the logistics network are shown, including the specific information of each service route.
In Table 6, the total number of vehicles used in the initial logistics network is 33, and the number of vehicles used at each facility is 5, 6, 6, 6, 4, and 6, respectively. In addition, some vehicles return to their origin early, such as V9, V10, V17, and V18, indicating that these vehicle resources are underutilized. e service vehicle routes of DC1 and the 27 customers it serves in the logistics network are shown in Figure 4.
In Figure 4, the service routes of V3, V4, and V5 are relatively complex. V5 performs delivery services for customers C23, C22, C6, C16, C2, and C5. However, these customers may be closer to DC3 on the basis of the perspective of spatial distribution. erefore, optimizing this logistics network is necessary. e values of the relevant parameters used in this real-world case study are shown in Table 7 [38,64,65,68].

Optimization Results.
Clustering customers to optimize resource allocation is the first step in optimizing the logistics network.
e customer clustering results of this logistics network by k-means algorithm are shown in Figure 5.
In Figure 5, the service relationship between customers and facilities is optimized by clustering. Each customer is covered by the logistics facility that is located close to that customer. On the whole, the service area of each logistics facility has been obviously allocated. For example, C23, C22, C6, C16, and C2 are served by DC1 before clustering; however, they are also served by DC3. Statistical analysis of customers whose service relationship has changed like those five customers is the key to handle centralized transportation. e details of the amount of goods transferred among facilities are shown in Figure 6.
In Figure 6, the amount of goods shipped from DC1 to DC2 is zero; thus, the initial customers of DC1 are not allocated to DC2 with clustering. e amount of goods transported from DC2 to DC1 is 170, which indicates that the service relationship of some customers has changed from DC2 to DC1. Here, the spatial distance between DC1 and DC2 is 33.07 km. e centralized transportation service routes of the truck are listed in Table 8.
In Table 8, T1 and T2 serve three routes. At the end of each service route, the trucks should return to their origin. e activity of T1 occurs early on each workday, and that of T2 occurs late on each workday. Table 9 shows the optimized vehicle service routes.
In Table 9, the number of vehicles used jointly by the six logistics facilities is 12. Some of the vehicles are used multiple times within and between the facilities. For example, V1 performs the route DC1 ⟶ C86 ⟶ C13 , which occur in DC1 and DC3, respectively. V10 provides service for PC2 and PC3 successively. V2, V4, and V6 are shared in DC1, DC2, and DC3, respectively.
To clarify the effect of the proposed model and algorithm, the gap of the cost and the number of vehicles in the initial and optimized logistics network are counted and listed in Table 10. Here, the TOC of the logistics facility includes the TTC, TFC, TPC, and TMC. Centralized transportation is a special project generated by the sharing of customer information and transportation resources among facilities. erefore, the costs of centralized transportation should be jointly borne by all the facilities participating in the sharing. Similarly, given that Input: e datasets of customer and facility information including location, time window, demand, and other relevant parameters Output: e initial feasible solution (1) Step 1: Mark the initial distance saving values with numbers 1, 2, 3, . . ., (n-1) (2) Step 2: Assign vehicles to each customer (3) Step 3: Calculate the distance savings of new vehicle routes, which are formed by any two routes (4) Step 4: Sort the distance savings in descending order (5) Step 5: If new vehicle routes meet the time window and capacity constraints, then go directly to Step 6; otherwise set n�n+1 and return to Step 5 (6) Step 6: Generate a new vehicle route with two vehicle routes whose distance savings are maximum (7) Step 7: Update the distance savings value through the fusion of vehicle routes (8) Step 8: Generate new vehicle route with two vehicle routes whose distance savings are maximum (9) Step 9: n � 1 (10) Step 10: If new vehicle routes meet the time window and capacity constraints, then go directly to Step 11; otherwise, return to Step 7 (11) Step 11: n � n+1 (12) Step 12: Generate new vehicle route with two vehicle routes whose distance savings are maximum (13) Step 13: If there exists a route that serves only one, then go back to Step 7; otherwise enter Step 14 (14) Step 14: Output the service route of each vehicle ALGORITHM 2: Procedure of CW algorithm.
Input: Initial population P Output: e nondominated sorting results of population P (1) Step 1: Setting PP�P, rank�1, F�1 (2) Step 2: for F�1:n (3) Step 3: for i�1:n (4) Step 4: for j�1:n && j≠i (5) Step 5: Compare the solution quality of x i and x j , and determine the dominant and nondominant relationships of i and j. End (6) Step 6: If x i is superior to all x j , then the individual i is considered to be a nondominant individual, and the nondominant ranking value of i is rank End (7) Step 7: e nondominant individuals found in the above steps constitute set S F , which is regarded as the F-level nondominant layer of the population P (8) Step 8: PP�PP\S F , i r �i r +1 (9) Step 9: If PP is an empty set, then enter Step 10, otherwise continuous this cycle (10) Step 10: Stratify population P (11) Step 11: End ALGORITHM 3: Procedure of the fast nondominated sorting operator. the vehicle is shared and some vehicles are used multiple times, the maintenance of the vehicles should be co-paid by all facilities.
In Table 10, the fixed costs of each facility are stable, which are $315, $427, $533, $578, $612, and $590, respectively. Reducing the travel, penalty, and maintenance costs is the main objective of optimizing the logistics operating cost. e TOC of the initial logistics network is $25473, whereas the optimized TOC is $16614, which indicates that the logistics network is significantly improved. Transportation resources are greatly saved, as the number of vehicles before and after optimization is 33 and 14, respectively. e gap of TTC before and after optimization is shown in Figure 7, which can directly prove the optimization effect of vehicle routes.     Journal of Advanced Transportation 13 In Figure 7, the initial TTC of DC1, DC2, DC3, PC1, PC2, and PC3 is $2210, $3581, $3730, $4900, $3067, and $3506, respectively. e optimized TTC of each facility is $2047, $1738, $1134, $2064, $2025, and $2361, respectively. e difference of TTC before and after optimization is relatively obvious. Although the transportation cost generated by the centralized transportation is $1651, the logistics network is still significantly optimized on the whole.

Analysis and Discussion.
e optimization of the logistics network is divided into four cases based on the RS strategy to analyze the influence of RS on the optimization of the logistics network. In Case 1, customer information and transportation resources are privately owned by each logistics facility. In Case 2, RS is adopted by DCs (i.e., DC1, DC2, and DC3). Correspondingly, in Case 3, the members participating in RS are PCs (i.e., PC1, PC2, and PC3). In Case 4, all six facilities receive the RS strategy. Transportation resources can be used jointly by DCs and PCs in Case 4. Table 11 and Figure 8 show the numerical optimization results of the logistics network in the four scenarios.
In Table 11 and Figure 8, Case 4 outperforms Cases 2 and 3 in most aspects. First, the TOC of Cases 2 and 3 is $17038 and $16977, respectively, which are larger than    $13020 in Case 4. Second, the number of vehicles and the TMCs of Scenario 4 is 14 and $539, respectively, which indicates that Case 4 is more conducive to saving resources. Finally, the gap of TOC between Cases 4 and 1 is $8859. erefore, if the six logistics facilities adopt RS simultaneously, then the logistics costs can be optimized better and the transportation resources can be saved considerably.

Management Insights.
In this study, the RS strategy optimizes the logistics network significantly by reallocating customers, including the logistics operation costs and the number of vehicles. erefore, the management insights obtained in this study are as follows: (1) In a multi-depot, large-scale logistics network with pickups and deliveries, customer information, facility capacity, and transportation resources can be shared to amplify resource utilization by introducing RS strategies. RS is not only conducive to the operation of logistics facilities, including the use of logistics costs and transport resources, but also conducive for providing customers with more convenient logistics services. In a logistics network with vehicle sharing, vehicles are used multiple times within and between facilities to avoid idle  Gap  TOC  NV  TOC  NV  TOC  NV  TTC  TFC  TPC  TMC  TTC  TFC  TPC  TMC  DC1  2210  315  27  179  5  2047  315  0   429  12  8859  19   DC2  3581  427  55  214  6  1738  427  0  DC3  3730  533  44  214  6  1134  533  0  PC1  4900  578  65  214  6  2064  578  0  PC2  3067  612  22  143  4  2025  612  0  PC3  3506  590  33  214  6    vehicle resources. In addition, the sharing of customer information helps facilities provide logistics services to customers with higher quality. erefore, participating actively in RS is remarkably necessary for LSPs.
(2) From the perspective of local traffic management departments and environmental departments, the optimization of the logistics network through RS not only reduces transportation resources but also relieves the local traffic pressure. In addition, the efficient use of resources promotes the development of green environment. erefore, local government departments can actively support and give incentive policies to promote RS among local logistics facilities to optimize the logistics network and ease the traffic environment. Moreover, the effective sharing mechanism in the multicenter logistics network with pickups and deliveries can enhance the efficient operation of the logistics system and promote the green and sustainable development of the local intelligent logistics system. erefore, the introduction of RS can promote the construction and development of local smart and green urban logistics with the incentive of departments and the active participation of facilities.

Conclusions
is study handles the MDPDPRS, which optimizes the logistics network by sharing customer and transportation resources. e reconfiguration of resources and customers improve the operating efficiency of the logistics network. e main contributions of this study include the following aspects. First, the MDPDPRS is modeled as a bi-objective mathematical model to optimize the total logistics operating cost and number of vehicles. Second, a two-stage hybrid algorithm is designed to solve the MDPDPRS, which contains the k-means, CW, and NSGA-II algorithms. ird, the application of the proposed mathematical model and methodology are improved by the numerical result of a realworld case study and benchmark.
In view of the shared transportation resources and customers' information, a two-stage algorithm is designed, which initially clusters customers to reconfigure the resources and then optimizes the vehicle routes. On the one hand, k-means algorithm, which clusters customers based on Manhattan distance, helps reduce the difficulty of solving MDPDPRS and enables vehicle resources to be used in a centralized manner multiple times. On the other hand, the combination of CW and NSGA-II algorithms improves the global searching capability and the speed of the algorithm in finding the Pareto optimal solutions. e numerical results of a real-word case study, which is obtained in Chongqing, China, is discussed and analyzed to improve the application and performance of the designed mathematical model and methodology in solving practical problems similar to the MDPDPRS. e gap of the logistics operating costs and number of vehicles before and after the optimization are $8859 and 19, respectively, which verify the effectiveness of the model and methodology proposed in this study. In addition, the algorithm comparison results of benchmarks (from C-mdvrptw datasets) verify that the proposed KCW-NSGA-II algorithm is superior to the standard NSGA-II, GA-PSO, and MOPSO. e results of numerical discussion on the four cases in which different RS strategies are adopted prove that RS is helpful to optimize logistics costs and save transportation resources.
In this study, the bi-objective mathematical model and methodology of MDPDPRS are designed, which provide references for the reconfiguration of resources and the optimization of logistics operation costs. In view of the limitations of the current study and the dynamic development of the logistics industry, further research can be considered from the following aspects. (1) e dynamic change in customer demands and customer satisfaction are the two aspects that the realistic LSPs focus on, and these aspects can be added into the study of MDPDPRS. (2) Constructing a dynamic mathematical programming model and designing an exact algorithm to find the exact solution of MDPDPRS are worthy of research. (3) Exploring the approaches to realize RS in a large logistics network and the means to reduce the effect of logistics transportation on the environment can be considered in the study of MDPDPRS. (4) Considering the cost sharing mechanism under the RS mechanism to promote the formation of collaboration and maintain its stability can enrich the study of MDPDPRS.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they do not have any conflicts of interest.