Next Article in Journal
Advanced Fault Detection in Power Transformers Using Improved Wavelet Analysis and LSTM Networks Considering Current Transformer Saturation and Uncertainties
Previous Article in Journal
Information Criteria for Signal Extraction Using Singular Spectrum Analysis: White and Red Noise
Previous Article in Special Issue
A Sim-Learnheuristic for the Team Orienteering Problem: Applications to Unmanned Aerial Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

AASA: A Priori Adaptive Splitting Algorithm for the Split Delivery Vehicle Routing Problem

by
Nariman Torkzaban
1,*,†,
Anousheh Gholami
1,†,
John S. Baras
1 and
Bruce L. Golden
2,*
1
Department of Electrical and Computer Engineering, University of Maryland, College Park, MD 20742, USA
2
Robert H. Smith School of Business, University of Maryland, College Park 20742, MD, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Algorithms 2024, 17(9), 396; https://doi.org/10.3390/a17090396
Submission received: 6 June 2024 / Revised: 17 August 2024 / Accepted: 29 August 2024 / Published: 6 September 2024
(This article belongs to the Special Issue Heuristic Optimization Algorithms for Logistics)

Abstract

:
The split delivery vehicle routing problem (SDVRP) is a relaxed variant of the capacitated vehicle routing problem (CVRP) where the restriction that each customer is visited precisely once is removed. Compared with CVRP, the SDVRP allows a reduction in the total cost of the routes traveled by vehicles. The exact methods to solve the SDVRP are computationally expensive. Moreover, the complexity and difficult implementation of the state-of-the-art heuristic approaches hinder their application in real-life scenarios of the SDVRP. In this paper, we propose an easily understandable and effective approach to solve the SDVPR based on an a priori adaptive splitting algorithm (AASA) that improves the existing state of the art on a priori split strategy in terms of both solution accuracy and time complexity. In this approach, the demand of the customers is split into smaller demand values using a splitting rule in advance. Consequently, the original SDVRP instance is converted to a CVRP instance which is solved using an existing CVRP solver. While the proposed a priori splitting rule in the literature is fixed for all customers regardless of their demand and location, we suggest an adaptive splitting rule that takes into account the distance of the customers to the depot and their demand values. Our experiments show that AASA can generate solutions comparable to the state of the art, but much faster.

1. Introduction

The split delivery vehicle routing problem (SDVRP) introduced by [1] is a relaxation of the traditional capacitated vehicle routing problem (CVRP), where each customer can be visited more than once. In both CVRP and SDVRP, a number of identical vehicles having limited capacity serve a set of customers with given demands. The vehicles depart a depot and return to the depot after visiting customers. The pairwise travel costs between customers and between customers and the depot are given. The objective is to minimize the total travel cost of the vehicles. The difference between CVRP and SDVRP is that in contrast to CVRP, where each customer is required to be visited once and by only one vehicle, SDVRP allows multiple visits to a customer. Therefore, the demand of a customer may be split among multiple vehicles. The authors in [2] showed that splitting the customer demands potentially reduces the total cost by up to 50 % . Furthermore, they show that there always exists an optimal solution to the SDVRP in which there is no k-split cycle and no two routes have more than one common customer.
Apart from environmental benefits, such as reduced emissions and energy efficiency, by allowing split deliveries, SDVRP can reduce the total distance traveled and the number of vehicles required. This optimization leads to lower fuel consumption, reduced vehicle wear and tear, and decreased maintenance costs. For businesses with large delivery fleets, these savings can be substantial. Split delivery may reduce labor, inventory, and storage costs as efficient routing reduces the time drivers spend on the road, which translates to lower labor costs. Reduced travel time also minimizes overtime and improves overall productivity. Moreover, split deliveries enable better alignment with demand patterns. By avoiding the need to stockpile large quantities of inventory, businesses can reduce warehousing costs and minimize the risk of overstocking or obsolescence. These findings justify the essence of formulating and solving the SDVRP.
State-of-the-art research has applied the SDVRP to various practical scenarios. For instance, Mullaseril et al. [3] addressed the problem of distributing feed to cattle on a large ranch in Arizona, where inaccuracies in feed loading necessitate split deliveries. Their approach, modeled as a capacitated rural postman problem with time windows, demonstrated that allowing split deliveries significantly reduced the total distance traveled by the fleet in most cases. Sierksma and Tijssen [4] studied the crew exchange process between offshore gas platforms in the North Sea. They formulated this as a discrete split delivery routing problem and employed integer programming with column generation and a cluster-and-route procedure. Their experiments showed that while their column generation approach (CGRP) was more accurate, the cluster-and-route (CR) method offered faster computation times. Song et al. [5] explored the distribution of newspapers from printing plants to agents in South Korea. They implemented a two-phase solution procedure, including agent allocation and split delivery scheduling, which improved the overall delivery efficiency and reduced costs by 15% compared to manual methods. More recently, within the context of waste management, ref. [6] utilized the SDVRP framework to model and solve complex waste collection and transportation challenges, incorporating realistic collection policies and offering insights into optimizing operational costs and fleet management, where they showed how an appropriate collection policy can significantly reduce fleet operating costs by 40–60% and fleet size by approximately 50% without causing overflow.
The SDVRP is known to be an NP-hard problem [7]. Therefore, exact solutions are only applicable to the SDVRP instances of smaller scale. Such methods typically formulate the SDVRP as a mixed integer linear program (MILP) [8] and propose different variants of branching algorithms such as branch-and-cut [9,10], branch-and-price [11,12], and branch-and-price-and-cut [13,14] for various versions of the SDVRP. However, with increasing the size of the problem, the corresponding MILP models will soon become intractable for exact solutions. To tackle this issue, various traditional approximation techniques, heuristics and metaheuristics have been proposed to deal with the high time-complexity of the large SDVRP models. For instance, a max–min clustering before solving the optimization problem is proposed in [15]. Tabu search has been used in several papers [16,17,18]. In [16], the authors proposed an approach relying on simulated annealing to solve the SDVRP. Other lines of research propose methods utilizing ant colony optimization [19,20], genetic algorithms [21,22], and particle swarm optimization [23,24], to solve the SDVRP. Using column generation is another optimization technique that has been applied to the SDVRP by [8,25]. The authors in [26] applied the cutting plane method to tackle the time complexity of the SDVRP.
Recently, the authors in [27] introduced a metaheuristic approach incorporating several mathematical programming components within an iterative local search framework. We refer to this method as ILS-MIP in the rest of the paper. ILS-MIP starts from an initial set of solutions obtained using the I1 Solomon heuristic followed by a perturbation step (to escape local optima), a local search step using classical neighborhood search heuristics, and several steps based on hybrid components. The hybrid components are called only if the incumbent solution of the search is not improved after a number of steps. The hybrid components encompass a mixed integer program (MIP)-based improvement heuristic performing two operations of inserting splits or removing potential unnecessary splits. Another used hybrid component is based on converting the current SDVRP solution to a CVRP instance and using a hybrid genetic search (HGS) framework [28] specialized to solve a CVRP. An MIP model is also used in the third component of the hybrid step where a residual problem is constructed by removing the edges that are not used in any solution of the current set of solutions. The residual problem is then solved using an MIP solver. Finally, the fourth component uses the branch-and-cut (BC) framework proposed in [29]. While ILS-MIP is a powerful and effective methodology, its complexity and significant runtime for larger instances makes it challenging for addressing real-world problems.
Although the state-of-the-art algorithms for the SDVRP are known for their ability to produce high-quality solutions, they often pose challenges in terms of complexity and practical implementation. In response to these challenges, the authors in [30] introduced an efficient algorithm that decomposes the SDVRP into two sub-problems. First, a demand-splitting rule is applied to all customers. The idea of splitting the customer demands was introduced in [31]; the authors propose a heuristic solution for SDVRP with time windows based on the column generation approach. In this paper, it is assumed that the demand of each customer is served by a number of orders selected from a set of feasible orders where each order is a pre-determined discrete split of the customer demand. The authors in [30] proposed applying the a priori splitting strategy in which the demand of each customer is split into smaller demands, each representing a new customer located in the same position as the original customer. As a result, a new problem instance is generated with an increased number of customers. Second, the new problem is assumed to be an instance of the CVRP which is solved using existing powerful CVRP solvers. The transformation of the SDVRP to a CVRP not only facilitates the implementation of SDVRP but also enables leveraging the rich literature on the CVRP. The obtained solution for the resulting CVRP is then translated back to the original SDVRP.
Regarding the splitting strategy, the solutions proposed in both [30,31] are based on a static splitting rule without taking into account the specific characteristics of the problem, such as customers’ location and demand. In the method proposed by [30], a fixed splitting rule that is inspired by the US coin denominations is applied to all customer demands and the VRPH solver introduced by [32] is used to solve the resulting CVRP instance. We refer to this algorithm as VRPHAS in the rest of this paper. The numerical evaluation of VRPHAS illustrates its ability to produce acceptable sub-optimal solutions in much less time compared to the state-of-the-art heuristic algorithms. However, the adoption of the introduced splitting rule is not well justified. Instead of applying a fixed rule to all customers, a splitting rule that is adaptive, based on the specific characteristics of a customer can potentially result in an improved solution for the SDVRP. In this paper, we study the problem of defining a good splitting rule for each customer. To this end, we introduce an a priori adaptive splitting algorithm (AASA). In combination with the VRPH solver, AASA results in an improved performance in terms of the optimality gap without significantly increasing the computational complexity of VRPHAS. AASA achieves this by taking into account the distance between the customers and the depot, the value of the customers’ demand along with the vehicles’ capacity. Similar to VRPHAS, AASA can be used with any CVRP solver. Thus, it provides a method to solve the SDVRP that is easily understandable, can be implemented simply, and generates high-quality solutions very quickly.
The rest of the paper is organized as follows. In Section 2.1, we provide a formal definition of the SDVRP. Our proposed solution is explained in Section 2.2. Section 3 presents the numerical results. Finally, the conclusion is discussed in Section 4.

2. Problem Definition and the Proposed Solution

2.1. SDVRP Definition

Let G = ( V , E ) denote an undirected and weighted complete graph representing the network of customers and a depot. The vertex set V = { 0 , , n } represents the depot denoted by the vertex 0 and n customers indexed as 1 , , n . A non-negative weight c i j is associated with each edge of the graph ( i , j ) E that stands for the cost of traveling between customers i and j. Moreover, each customer i V has a positive demand, denoted by d i . A set of M vehicles is available to serve the customers’ demands. Each vehicle has a limited capacity denoted by Q. The objective of the SDVRP is to find a route for each vehicle such that the demands of all customers are satisfied with the minimum total cost. In this paper, we assume that the cost of traveling between customers i , j V is their Euclidean distance denoted by the function d i s t , i.e.,  c i j = d i s t ( i , j ) .

2.2. Proposed Heuristic Algorithm for SDVRP

In the following, we first describe the idea of the a priori splitting rule proposed by [30] and then discuss the core algorithmic idea of AASA based on two motivational observations. We then formalize AASA in detail.

A Priori Splitting Rule

The proposed splitting rule in [30] is inspired by the US coin denominations. The authors propose two splitting options 20 / 10 / 5 / 1 , and  25 / 10 / 5 / 1 . The first option replaces customer i, with the demand of d i , by m 20 i , m 10 i , m 5 i and m 1 i customers with the demand values of 0.2 Q , 0.1 Q , 0.05 Q and 0.01 Q , respectively. Hence, we have:
m 20 i = m a x { m Z + { 0 } : 0.2 Q m d i } ;
m 10 i = m a x { m Z + { 0 } : 0.1 Q m d i 0.2 Q m 20 i } ;
m 5 i = m a x { m Z + { 0 } : 0.05 Q m d i 0.2 Q m 20 i 0.1 Q m 10 i } ;
m 1 i = m a x { m Z + { 0 } : 0.01 Q m d i 0.2 Q m 20 i 0.1 Q m 10 i 0.05 Q m 5 i } .
A similar breakdown can be given for the second splitting rule option. Using the splitting rule 20 / 10 / 5 / 1 , the demand of customer i is split into m i = m 20 i + m 10 i + m 5 i + m 1 i smaller values. The customer i is then replaced by m i customers located at the same position as customer i and with the demand values resulting from (1)–(4). After applying the splitting rule to all customers, a new graph is constructed with m = i = 1 n m i vertices (customers). The new graph is considered as an instance of the CVRP and the VRPH solver is used to solve it. The resulting solution is also a solution to the original SDVRP instance. In the next section, we justify the essence of the adaptation of the splitting rule and propose the AASA algorithm.

2.3. A Priori Adaptive Splitting Algorithm (AASA)

The a priori splitting rule described by Equations (1)–(4) is a fixed rule, i.e., a single splitting rule is applied to all customers regardless of their demand and location. However, the information regarding the customers’ demand and location can be leveraged in the vehicle routing decisions towards reducing the total travel cost as well as the size of the resulting CVRP instance. We propose an adaptive splitting rule based on the location and demand information that is shown to effectively improve the solution. In the following, we first present the motivation behind the main building blocks of our proposed splitting rule. We then explain the AASA algorithm.

2.3.1. Motivation

In this section, we present the motivation behind our proposed splitting algorithm through illustrative examples. We study the impact of the splitting rule granularity on the total travel cost and the solution runtime considering customers’ demand values and their distances from the depot. Our investigations reveal the benefits of using coarser/finer splitting rules for different customers depending on their demand and location. By a coarser splitting rule, we mean one that entails larger splitting portions as opposed to a finer splitting rule where the portions are smaller. For instance, for a demand point with d = 200 , the splitting rule 200 = 128 + 64 + 8 is coarser compared to the rule 200 = 3 × 64 + 8 .
Impact of coarser rules for customers with higher demands: Using coarser rules for customers with high demands introduces a two-fold benefit: (i) it reduces the runtime due to a smaller number of resulting customers in the resulting CVRP and (ii) it often results in a lower total cost of delivery. To capture this effect, we employ an exponential pattern (such as 128/64/32/16/8/4/2) for the AASA rule. We illustrate this effect for the P 04 _ 3070 instance of the third dataset provided in [33], as depicted in Figure 1. In this instance, the vehicle capacity is 200 and the depot is located at the location ( 0 , 0 ) in Figure 1a,b. Each point in the plane corresponds to a customer. The routes of the vehicles in the obtained solution are denoted by different colors. Figure 1a,b show the solution of the SDVRP instance considering the 20 / 10 / 5 / 1 and 128 / 64 / 32 / 16 / 8 / 2 splitting rules, respectively. We observe that the splitting rule of 128 / 64 / 32 / 16 / 8 / 4 / 2 reduces the solver runtime by an order of 2.5 . The lower runtime is expected since using an exponential pattern in the splitting rule results in splitting the higher demands into a piece with a very high value and most likely fewer pieces in total. This effect can also be observed from the size of the generated CVRP instance. The  20 / 10 / 5 / 1 rule resulted in a new instance with 855 customers (nodes), while the 128 / 64 / 32 / 16 / 8 / 4 / 2 rule generated a network with 572 customers. Moreover, the coarser splitting rule of 128 / 64 / 32 / 16 / 8 / 4 / 2 resulted in a cost of 4424.64 that improves upon the cost of the rule 20 / 10 / 5 / 1 (proposed by the VRPHAS algorithm) by 1.34 % .
Impact of coarser rules for customers further away from the depot: The second motivation behind our proposed AASA rule is that using coarser rules for the demands with a larger distance from the depot and splitting closer demands into smaller pieces results in lower solver runtime and possibly a better solution. The intuition behind this observation is that as the distance between a customer and the depot increases, each visit to that customer would significantly increase the cost. Therefore, it is desirable to satisfy the demand of further customers within the fewest possible number of visits and that is only possible by applying coarser splitting rules to those nodes. Moreover, by splitting closer-to-depot customer demands into smaller pieces, we provide the opportunity for the vehicles returning to the depot to utilize their remaining capacity to fulfill smaller demands and, therefore, minimize the fraction of their unused capacity. It is important to note that all of these intuitive statements are based on the fact that after splitting demands and generating a new network with a larger number of customers, we use the VRPH solver and basically solve the generated instance by a number of heuristic and metaheuristic methods. Figure 2a–c illustrate the above observation for the S D 6 instance from [34]. The best-known solution for this instance is 830.86 . In Figure 2a,b, the results of using the 80 / 40 / 20 / 10 and 40 / 20 / 10 rules are shown. Although these two rules alone provide acceptable results with the costs of 863.43 and 881.21 , respectively, we observe that by using a combination of the two rules as shown in Figure 2c, we can achieve lower cost of 862.38 . In this experiment, the 32 customers are clustered into two categories of size 16. For the 16 customers closer to the depot, we used the splitting rule 40 / 20 / 10 , and the coarser splitting rule 80 / 40 / 20 / 10 is applied to the remaining 16 customers located farthest from the depot.

2.3.2. AASA Algorithm

The above observations highlight the significance of applying an adaptive splitting rule for different customers instead of using a fixed splitting rule for all customers, as in VRPHAS. Motivated by these observations, we incorporate the following factors into AASA:
  • Vehicle capacity: As implied by [30], it makes sense to look for the splitting rule as a function of the vehicle capacity. This not only presents a structured representation with vehicle capacity as the kernel allowing for better generalization but also leads the policy towards leaving the lowest portion of the vehicles’ capacities unused, and consequently, increases the vehicles’ utilization.
  • Customer location: We believe one potential improvement to the VRPHAS splitting rule is to personalize the rule for each customer considering its distance from the depot. We observe that if the customers that are located further from the depot are visited by one (or very few) vehicle(s), the total travel cost is lower. Otherwise, if a long distance needs to be traversed multiple times, the cost will be negatively impacted. This can be avoided by coarsening the splitting rule for further customers. In AASA, the impact of customers’ locations is addressed by partitioning the customers into a number of clusters based on their distance from the depot. We then use a different splitting rule for each cluster. A coarser splitting rule is used for clusters of customers further from the depot while finer rules are applied to the clusters closer to the depot.
  • Customer demand: Another important factor in the adaptation of the splitting rule is the demand of customers. Depending on how large the demand is, the splitting rule can be applied in full or partially to result in an appropriate number of additional demand points with appropriate demand values. In order to create the exponential pattern in the splitting rule, we propose to decompose the customers’ demands based on the different powers of prime numbers. For instance, the demand of a customer with d = 199 served by vehicles with capacity Q = 200 can be decomposed as d = 2 7 + 2 6 + 2 2 + 2 1 + 1 if p = 2 , and  d = 2 × 3 4 + 3 3 + 3 2 + 1 if p = 3 .
Leveraging the above factors, the AASA solution framework is shown in Figure 3. First, AASA employs a clustering algorithm to partition the customers into multiple levels. In our experiments, we used a simple distance-based clustering as follows:
l a b e l ( v ) = , if 1 L d i s t m a x < d i s t ( v , 0 ) L d i s t m a x , { 1 , , L }
where L is the number of levels, and  d i s t m a x is the maximum customer–depot distance among all customers, i.e, d i s t m a x = m a x v d i s t ( v , 0 ) . As a result of this clustering rule, first, the customers are separated with L uniformly spaced rings around the depot. Second, for the customers at each level, the same splitting rule is applied. After all the demands are split, a CVRP instance is generated. Finally, the VRPH solver is used to solve the resulting CVRP instance.
The pseudo-code for AASA is given in Algorithm 1. The algorithm starts by taking as input the SDVRP instance as a graph G = ( V , E ) , the vehicle capacity Q, the vector of customer–depot distances denoted by d i s t where d i s t ( v ) = d i s t ( v , 0 ) , and the vector of customer demands D . Next, it computes d, the greatest common divisor ( g c d ) of the demands vector and the vehicle capacity, which will be used as a parameter in the adaptive splitting rule. The vehicles may carry goods only in quantities that are multiples of d units. This is to control the granularity of the splitting rule. Next, the demand vector is scaled down by d and averaged over the customers to find μ that is used in finding s m a x determining the largest quantity of goods that is demanded by the customers in the equivalent CVRP model.
Algorithm 1 A priori Adaptive Splitting Algorithm (AASA)
Input:  G = ( V , E ) , Q , d i s t , D , p
Output:  G = ( V , E )
   1:
Initialize L , d g c d ( Q , D ) , μ a v g ( D / d ) , s m a x [ l o g p μ ]
   2:
for i : 1 . . L
   3:
    s i d ( e x p p { 0 . . s m a x i + 1 } )
   4:
end for
   5:
l a b e l = c l u s t e r i n g ( V )                                                       ▹ l a b e l ( v ) { 1 , . . . , L }
   6:
repeat
   7:
    V = V { v }
   8:
    i l a b e l ( v )
   9:
    V = V s p l i t ( v , s i )
  10:
until   V =
  11:
return   G = ( V , E )
The algorithm takes into account L different levels for the distance-based clustering as provided in (5). Then, for each level i, the splitting rule s i is determined using the value s m a x . Next, the clustering of the customers is made according to Equation (5). The split function substitutes each node (customer) of the graph G with new customers according to the rule s i to gradually form the graph G . Note that the splitting rule ensures that the demands further away from the depot are split with coarser rules. Once the graph G is completed, the CVRP solver will solve the generated new problem. We experimented with different values of the prime number p and observed that p = 2 has the best performance overall in the current version of AASA.
For example, consider the instance S D 6 from SET-1 as shown in Figure 2. We will have d = g c d ( D , Q ) , i.e., d = g c d ( 60 , 90 , 10 ) = 10 . The simple distance-based clustering in Equation (5) with L = 2 results in two levels, where l = 1 covers the customers that are furthest from the depot and l = 2 covers the customers that are closest to the depot. We will have μ = 7.5 and s m a x = 3 . Therefore, for the first level, we have s 1 = { 10 , 20 , 40 , 80 } and for s 2 = { 10 , 20 , 40 } .
VRPH solver: Once the a priori splitting rule is applied to an SDVRP problem, any commercial solver can be used to solve the resulting CVRP instance. In this paper, we use the VRPH solver to solve CVRP instances. VRPH is a publicly available solver with a provable record of generating high-quality solutions. It takes as input a CVRP instance written to a file with a format similar to that of the Traveling Salesman Problem Library (TSPLIB) files and prints the solution to an output file. VRPH implements an open-source library of several local search heuristics for generating and improving feasible solutions to the CVRP instances. We refer the reader to [35] for a detailed description of the heuristics and the modular structure of the VRPH software. In this paper, we use VRPH as a standalone solver with the parameters set at their default values.

3. Performance Evaluation

In this section, we evaluate the performance of AASA through extensive numerical simulations. We benchmark the performance of AASA against three baseline strategies, using instances from datasets that are widely used in the literature for the comparison of various vehicle routing solutions. Table 1 presents the summary of the evaluation datasets used throughout this section. Further details about the instances of each dataset are provided in Section 3.1.

3.1. Benchmarking Instances

We consider 4 benchmarking sets that include 88 instances overall as presented in Table 1. SET-1 [34] contains 21 instances, each with customers uniformly distributed on the perimeter of concentric circles centered at the depot. The instances are sorted by the number of customers and the number of concentric circles. SET-2 [36] and SET-3 [33] instances use similar customer demand profiles but vary in the number of customers and vehicle capacity. These demands are randomly sampled from the range [ a Q , b Q ] with ( a , b ) chosen by six different scenarios as in Table 1. SET-2 contains 14 and SET-3 contains 42 instances, 6 of which are repeated. Therefore, we test only the 36 of them that are not redundant. The coordinates of the customers in SET-2 are randomly generated using the coordinates of eil51, eil76, and eil101 from TSPLIB, which are also used in SET-4 [37]. SET-4 includes 11 instances where no specific rule is preserved throughout for generating the customer demands.

3.2. Benchmarking Solutions

We benchmark the performance of the AASA algorithm against the following baseline methods:
  • No-splitting (CVRP): The baseline method that treats an SDVRP instance as a CVRP instance and runs the VRPH [32] solver to generate solutions.
  • VRPHAS: The method based on a priori splitting that is proposed by [30].
  • ILS-MIP: The metaheuristic method proposed by [27].
We note that when computing the optimality gap corresponding to each method for each experimental instance, we take the best solution known in the literature as the reference solution. Furthermore, we note that in ILS-MIP, the travel cost between each pair of customers is assumed to be their rounded Euclidean distance except for SET-1 instances where the exact Euclidean distance is taken as the cost. Hence, we compare the results of ILS-MIP and AASA only for SET-1. For the rest of the datasets, we compare the performance of AASA against the No-splitting and the VRPHS methods.

3.3. Metrics and Setup

We consider the optimality gap and the solver runtime as metrics for performance comparison between the above strategies. The optimality gap which indicates the deviation from the best-known solution is defined as:
g a p = o b j best-known solution best-known solution 100
where o b j is the objective value of the considered method. We run both VRPHAS and AASA algorithms on the same PC with an Intel Xeon processor at 3.2 GHz and 16 GB of main memory. It is important to note that the gap results we obtained for the VRPHAS method do not match the values presented in [30] for some instances. However, since we need to compare both the runtime and the gap results for the VRPHAS and AASA methods on similar hardware, we reflect the results of our run for VRPHAS. Moreover, we report the computation time of the ILS-MIP method presented in [27], including the values of time and time best corresponding to the time required for the proposed iterative algorithm to terminate and the time at which the best solution of ILS-MIP is achieved, respectively. A time limit (TL) is considered for the algorithm termination that is assumed to be equal to 1349 s by [27]. It is important to note that the ILS experiments are conducted on a PC with an Intel Core i7-8700 3.2 GHz processor and 32 GB of memory which is more powerful than our PC. Therefore, we can expect that ILS-MIP performs even worse in terms of computation time compared to AASA if both were run on the same hardware.

3.4. Numerical Results

3.4.1. CVRP Instance Attributes Comparison: AASA vs. VRPHAS

As illustrated in Section 2.3.1, one of the promises of our proposed splitting algorithm is generating smaller size CVRP instances that can be solved consistently faster. To show this, we present the size of the CVRP instances and the number of routes (i.e., vehicles) in the generated CVRP solution obtained by the splitting strategies of VRPHAS and AASA for the evaluation sets SET-1 to SET-4 in Table 2, Table 3, Table 4 and Table 5, respectively. It is observed that for all the benchmarking datasets, AASA reduces the size of the CVRP problem compared to the VRPHAS method. Moreover, the number of routes in the solution obtained by AASA is less than or equal to the VRPHAS solution except for a few instances. Although the reduced number of routes indicates the potential of AASA in producing better solutions compared to VRPHAS, it does not provide sufficient evidence for cost-efficiency since the actual routes taken by the vehicles determine the total cost. However, together with the solution gap results shown in the next section, we can conclude that AASA performs better than VRPHAS in terms of both solution quality and runtime and that is a result of the adaptive splitting approach we proposed.

3.4.2. Optimality Gap and Runtime Comparison: AASA against Baselines

Table 6 presents the performance of AASA against VRPHAS, ILS-MIP, and the no-splitting strategies for the instances of SET-1. First, we observe that the no-splitting method results in very low-quality solutions as the resulting optimality gap obtained by this method is very high compared to the VRPHAS, AASA, and ILS-MIP heuristics. Moreover, we observe that while the ILS-MIP strategy outperforms the AASA and VRPHAS solution approaches in terms of optimality gap for almost all instances of SET-1, it requires very large computation times. The average runtime for ILS-MIP is 992.38 s as opposed to the average runtime of 86.14 s and 30.56 s corresponding to the VRPHAS and AASA solution approaches, respectively. This huge time complexity hinders the practical application of the ILS-MIP approach in real-world problems with time sensitivity or scalability requirements. Figure 4 and Figure 5 illustrate the per-instance optimality gap and runtime of our proposed splitting rule against the fixed splitting rule introduced in VRPHAS. Our proposed method outperforms the VRPHAS in terms of time complexity for all instances, and it also improves upon the average optimality gap of VRPHAS. Moreover, while the VRPHAS method results in the 5.50 % optimality gap in the worst-case scenario ( S D 16 ), AASA achieves the maximum optimality gap of 3.81 % for the same instance. We also observe that AASA finds the best-known solution for four instances as opposed to the VRPHAS method obtaining the best-known solution for three instances of SET-1.
Table 7 shows the results of VRPHAS, AASA, and the no-splitting methodologies for the instances of SET-2. According to Table 7, the no-splitting approach results in a higher gap compared to both the VRPHAS and AASA rules. Furthermore, our proposed AASA strategy outperforms the VRPHAS method in terms of the average optimality gap, as well as runtime. Moreover, the per-instance results depicted in Figure 6 and Figure 7 indicate that the AASA strategy improves upon the optimality gap and solver runtime of the VRPHAS method for almost all instances of SET-2. It can also be observed that while the VRPHAS method exhibits an optimality gap of 2.51 % in the worst-case scenario ( S 76 D 4 ), AASA achieves a maximum optimality gap of 1.97 % ( S 51 D 6 ).
Similarly, the comparison of the results for SET-3 and SET-4 provided in Table 8 and Table 9 show that our proposed splitting algorithm again improves both average optimality gap and runtime of the VRPHAS splitting mechanism and, therefore, it can effectively generate higher-quality solutions, faster than VRPHAS. The improvement in the runtime is mainly owing to the rule granularity consideration introduced in AASA as discussed in Section 2.2, which results in smaller instances to be solved by the CVRP solver compared to the VRPHAS method. Figure 8 and Figure 9 illustrate the per-instance optimality gap and runtime of SET-3. According to Table 8 and Figure 8, although the VRPHAS method results in the maximum optimality gap of 3.27 % ( p 11 _ 1090 ), the worst-case optimality gap achieved by AASA is 2.54 % ( p 02 _ 7090 ). The optimality gap and run-time of SET-4 instances are presented in Figure 10 and Figure 11. Table 9 and Figure 10 illustrate that the AASA method improves the worst-case optimality gap of VRPHAS in SET-4, since AASA achieves the maximum optimality gap of 2.02 % ( e i l A 101 ) as opposed to the maximum gap of 2.35 % obtained by VRPHAS for the same instance.

4. Conclusions

In this paper, we propose an efficient and effective heuristic algorithm for the SDVRP based on a two-step procedure of splitting the demand of customers into smaller demands and solving the generated network assuming it to be a CVRP instance. Our approach focuses on enhancing the state-of-the-art splitting rule by considering customer locations and demand values. Our proposed adaptive splitting algorithm incorporates clustering customers based on their location and choosing a specific splitting rule for each cluster. It is shown to outperform the fixed splitting rule of the approach adopted by the state of the art, for all four benchmark datasets on average, with respect to both performance and computational effort.

Author Contributions

Conceptualization, N.T., A.G., J.S.B. and B.L.G.; Methodology, N.T., A.G. and J.S.B.; Software, N.T. and A.G.; Validation, A.G. and B.L.G.; Formal analysis, N.T., A.G. and B.L.G.; Investigation, N.T. and A.G.; Resources, N.T. and A.G.; Writing – original draft, N.T., A.G. and B.L.G.; Writing – review and editing, J.S.B. and B.L.G.; Visualization, J.S.B. and B.L.G.; Supervision, J.S.B. and B.L.G.; Project administration, J.S.B.; Funding acquisition, N.T., A.G. and J.S.B. All authors have read and agreed to the published version of the manuscript.

Funding

The research of N. Torkzaban, A. Gholami, and J. Baras was partially supported by ONR grant N00014-17-1-2622, and by a grant from Lockheed Martin Chair in Systems Engineering.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dror, M.; Trudeau, P. Savings by split delivery routing. Transp. Sci. 1989, 23, 141–145. [Google Scholar] [CrossRef]
  2. Archetti, C.; Savelsbergh, M.; Speranza, M. Worst-case analysis for split delivery vehicle routing problems. Transp. Sci. 2006, 40, 226–234. [Google Scholar] [CrossRef]
  3. Mullaseril, P.; Dror, M.; Leung, J. Split-delivery routing heuristics in livestock feed distribution. J. Oper. Res. Soc. 1997, 48, 107–116. [Google Scholar] [CrossRef]
  4. Sierksma, G.; Tijssen, G. Routing helicopters for crew exchanges on off-shore locations. Ann. Oper. Res. 1998, 76, 261–286. [Google Scholar] [CrossRef]
  5. Song, S.; Lee, K.; Kim, G. A practical approach to solving a newspaper logistics problem using a digital map. Comput. Ind. Eng. 2002, 43, 315–330. [Google Scholar] [CrossRef]
  6. Luo, K.; Zhao, W.; Zhang, R. A multi-day waste collection and transportation problem with selective collection and split delivery. Appl. Math. Model. 2024, 126, 753–771. [Google Scholar] [CrossRef]
  7. Archetti, C.; Feillet, D.; Gendreau, M.; Speranza, M.G. Complexity of the VRP and SDVRP. Transp. Res. Part C Emerg. Technol. 2011, 19, 741–750. [Google Scholar] [CrossRef]
  8. Archetti, C.; Bianchessi, N.; Speranza, M.G. A column generation approach for the split delivery vehicle routing problem. Networks 2011, 58, 241–254. [Google Scholar] [CrossRef]
  9. Hernández-Pérez, H.; Salazar-González, J.-J. A Branch-and-cut algorithm for the split-demand one-commodity pickup-anddelivery travelling salesman problem. Eur. J. Oper. Res. 2021, 297, 467–483. [Google Scholar] [CrossRef]
  10. Wolfinger, D.; Salazar-González, J.-J. The pickup and delivery problem with split loads and transshipments: A branch-and-cut solution approach. Eur. J. Oper. Res. 2020, 289, 470–484. [Google Scholar] [CrossRef]
  11. Casazza, M.; Ceselli, A.; Calvo, R.W. A route decomposition approach for the single commodity Split Pickup and Split Delivery Vehicle Routing Problem. Eur. J. Oper. Res. 2019, 289, 897–911. [Google Scholar] [CrossRef]
  12. Li, J.L.; Qin, H.; Shen, H.X.; Tong, X.; Xu, Z. Exact algorithms for the multiple depot vehicle scheduling problem with heterogeneous vehicles, split loads and toll-by-weight scheme. Comput. Ind. Eng. 2022, 168, 108137. [Google Scholar] [CrossRef]
  13. Li, J.; Qin, H.; Baldacci, R.; Zhu, W. Branch-and-price-and-cut for the synchronized vehicle routing problem with split delivery, proportional service time and multiple time windows. Transp. Res. Part E Logist. Transp. Rev. 2020, 140, 101955. [Google Scholar] [CrossRef]
  14. Gschwind, T.; Bianchessi, N.; Irnich, S. Stabilized branch-price-and-cut for the commodity-constrained split delivery vehicle routing problem. Eur. J. Oper. Res. 2019, 278, 91–104. [Google Scholar] [CrossRef]
  15. Min, J.N.; Jin, C.; Lu, L.J. Maximum-minimum distance clustering method for split-delivery vehicle-routing problem: Case studies and performance comparisons. Adv. Prod. Eng. Manag. 2019, 14, 125–135. [Google Scholar] [CrossRef]
  16. Wang, J.; Jagannathan, A.K.; Zuo, X.; Murray, C.C. Two-layer simulated annealing and tabu search heuristics for a vehicle routing problem with cross docks and split deliveries. Comput. Ind. Eng. 2017, 112, 84–98. [Google Scholar] [CrossRef]
  17. Xia, Y.; Fu, Z. An adaptive tabu search algorithm for the open vehicle routing problem with split deliveries by order. Wirel. Pers. Commun. 2018, 103, 595–609. [Google Scholar] [CrossRef]
  18. Xia, Y.; Fu, Z.; Tsai, S.B.; Wang, J. A new TS algorithm for solving low-carbon logistics vehicle routing problem with split deliveries by backpack—From a green operation perspective. Int. J. Environ. Res. Public Health 2018, 15, 949. [Google Scholar] [CrossRef] [PubMed]
  19. Ma, X.; Liu, C. Improved Ant Colony Algorithm for the Split Delivery Vehicle Routing Problem. Appl. Sci. 2024, 14, 5090. [Google Scholar] [CrossRef]
  20. Yang, W.; Wang, D.; Pang, W.; Tan, A.-H.; Zhou, Y. Goods consumed during transit in split delivery vehicle routing problems: Modeling and solution. IEEE Access 2020, 8, 110336–110350. [Google Scholar] [CrossRef]
  21. Jiang, Y.; Bian, B.; Liu, Y. Integrated multi-item packaging and vehicle routing with split delivery problem for fresh agri-product emergency supply at large-scale epidemic disease context. J. Traffic Transp. Eng. (Engl. Ed.) 2020, 8, 196–208. [Google Scholar] [CrossRef]
  22. Fan, H.M.; Zhang, X.; Ren, X.X.; Liu, P. Optimization of multi-depot open split delivery vehicle routing problem with simultaneous delivery and pick-up. Syst. Eng.-Theory Pract. 2021, 41, 1521–1534. [Google Scholar]
  23. Shi, J.; Zhang, J.; Wang, K.; Fang, X. Particle swarm optimization for split delivery vehicle routing problem. Asia-Pac. J. Oper. Res. 2018, 35, 1840006. [Google Scholar] [CrossRef]
  24. Qing, D.S.; Deng, Q.L.; Li, J.J.; Liu, S.; Liu, X.; Zeng, S.P. Split vehicle route planning with full load demand based on particle swarm optimization. J. Control Decis. 2021, 36, 1397–1406. [Google Scholar]
  25. Jin, M.; Liu, K.; Eksioglu, B. A column generation approach for the split delivery vehicle routing problem. Oper. Res. Lett. 2008, 36, 265–270. [Google Scholar] [CrossRef]
  26. Bianchessi, N.; Drexl, M.; Irnich, S. The split delivery vehicle routing problem with time windows and customer inconvenience constraints. Transp. Sci. 2019, 53, 1067–1084. [Google Scholar] [CrossRef]
  27. Alvarez, A.; Munari, P. A Matheuristic Approach for Split Delivery Vehicle Routing Problems. 2022. Available online: http://www.optimization-online.org/DB_FILE/2022/02/8790.pdf (accessed on 31 August 2023).
  28. Vidal, T. Hybrid genetic search for the CVRP: Open-source implementation and SWAP* neighborhood. Comput. Oper. Res. 2022, 140, 105643. [Google Scholar] [CrossRef]
  29. Munari, P.; Savelsbergh, M. Compact formulations for split delivery routing problems. Transp. Sci. 2022, 56, 1022–1043. [Google Scholar] [CrossRef]
  30. Chen, P.; Golden, B.; Wang, X.; Wasil, E. A novel approach to solve the split delivery vehicle routing problem. Int. Trans. Oper. Res. 2017, 24, 27–41. [Google Scholar] [CrossRef]
  31. Salani, M.; Vacca, I. Branch and price for the vehicle routing problem with discrete split deliveries and time windows. Eur. J. Oper. Res. 2011, 213, 470–477. [Google Scholar] [CrossRef]
  32. Groer, C. VRPH. 2011. Available online: https://github.com/coin-or/VRPH (accessed on 1 August 2024).
  33. Archetti, C.; Speranza, M.G.; Savelsbergh, M.W. An optimization-based heuristic for the split delivery vehicle routing problem. Transp. Sci. 2008, 42, 22–31. [Google Scholar] [CrossRef]
  34. Chen, S.; Golden, B.; Wasil, E. The split delivery vehicle routing problem: Applications, algorithms, test problems, and computational results. Netw. Int. J. 2007, 49, 318–329. [Google Scholar] [CrossRef]
  35. Groer, C.S.; Golden, B.; Wasil, E. A Library of Local Search Heuristics for the Vehicle Routing Problem. Math. Program. Comput. 2010, 2, 79–101. [Google Scholar] [CrossRef]
  36. Belenguer, J.M.; Martinez, M.C.; Mota, E. A lower bound for the split delivery vehicle routing problem. Oper. Res. 2000, 48, 801–810. [Google Scholar] [CrossRef]
  37. Reinhelt, G. TSPLIB: A Library of Sample Instances for the TSP (and Related Problems) from Various Sources and of Various Types. 2014. Available online: http://comopt.ifi.uni-heidelberg.de/software/TSPLIB95 (accessed on 31 August 2023).
Figure 1. Impact of using splitting rules with different granularity for different demands.
Figure 1. Impact of using splitting rules with different granularity for different demands.
Algorithms 17 00396 g001
Figure 2. Impact of using splitting rules with different granularity according to the customers’ distance to the depot.
Figure 2. Impact of using splitting rules with different granularity according to the customers’ distance to the depot.
Algorithms 17 00396 g002
Figure 3. AASA framework.
Figure 3. AASA framework.
Algorithms 17 00396 g003
Figure 4. Per-instance gap of SET-1.
Figure 4. Per-instance gap of SET-1.
Algorithms 17 00396 g004
Figure 5. Per-Instance solver runtime of SET-1.
Figure 5. Per-Instance solver runtime of SET-1.
Algorithms 17 00396 g005
Figure 6. Per-instance gap of SET-2.
Figure 6. Per-instance gap of SET-2.
Algorithms 17 00396 g006
Figure 7. Per-instance solver runtime of SET-2.
Figure 7. Per-instance solver runtime of SET-2.
Algorithms 17 00396 g007
Figure 8. Per-instance gap of SET-3.
Figure 8. Per-instance gap of SET-3.
Algorithms 17 00396 g008
Figure 9. Per-instance solver runtime of SET-3.
Figure 9. Per-instance solver runtime of SET-3.
Algorithms 17 00396 g009
Figure 10. Per-instance AASA gap of SET-4.
Figure 10. Per-instance AASA gap of SET-4.
Algorithms 17 00396 g010
Figure 11. Per-instance solver runtime of SET-4.
Figure 11. Per-instance solver runtime of SET-4.
Algorithms 17 00396 g011
Table 1. Evaluation datasets.
Table 1. Evaluation datasets.
DatasetNumber of InstancesNumber of Customers (N)Vehicle Capacity (Q)Customers’ Demands
SET-121 [ 8 , 288 ] 100 { 60 , 90 }
SET-214 { 50 , 75 , 100 } 160randomly from [ a Q , b Q ]   1
SET-342 [ 50 , 199 ] [ 140 , 200 ] randomly from [ a Q , b Q ]   1
SET-411 [ 21 , 100 ] [ 112 , 8000 ] No pattern
1  ( a , b ) { ( 0.01 , 0.1 ) , ( 0.1 , 0.3 ) , ( 0.1 , 0.5 ) , ( 0.1 , 0.9 ) , ( 0.3 , 0.7 ) , ( 0.7 , 0.9 ) } , corresponding to six different cases.
Table 2. CVRP instance attributes: VRPHAS vs. AASA on SET-1.
Table 2. CVRP instance attributes: VRPHAS vs. AASA on SET-1.
InstanceCVRP SizeNumber of Routes
VRPHASAASAVRPHASAASA
SD1321887
SD264421413
SD364361313
SD496541920
SD5128842626
SD6128722928
SD71601103434
SD81921344141
SD91921263939
SD102561685556
SD113202226868
SD123202206868
SD133842728083
SD14480330103100
SD15576403123121
SD16576324128129
SD17640444138140
SD18640440141139
SD19768535167165
SD20960666210208
SD211152756259259
Table 3. CVRP instance attributes: VRPHAS vs. AASA on SET-2.
Table 3. CVRP instance attributes: VRPHAS vs. AASA on SET-2.
InstanceCVRP SizeNumber of Routes
VRPHASAASAVRPHASAASA
S51D113910933
S51D21671501010
S51D32091621717
S51D42462473232
S51D52532312828
S51D63072465050
S76D121817644
S76D22632451616
S76D33272832424
S76D43523524140
S101D127621855
S101D23453152121
S101D34283723433
S101D55004725856
Table 4. CVRP instance attributes: VRPHAS vs. AASA on SET-3.
Table 4. CVRP instance attributes: VRPHAS vs. AASA on SET-3.
InstanceCVRP SizeNumber of Routes
VRPHASAASAVRPHASAASA
p01_1101219233
p01_10301761481111
p01_10502071931717
p01_10902682043333
p01_30702571983333
p01_70903262405050
p02_11018216855
p02_10302042311717
p02_10502972832727
p02_10903823774847
p02_30703702934848
p02_70904623777575
p03_11034624466
p03_10304113332323
p03_10504713763737
p03_10906084766767
p03_30705713696666
p03_709071159110099
p04_1105163151010
p04_10306134673434
p04_10507055475555
p04_10909156769898
p04_30708556499898
p04_70901063879149149
p05_1106574751212
p05_10308046164343
p05_10509206216769
p05_10901169824123124
p05_30701119806124126
p05_709014011082199197
p11_11042223388
p11_10304893882727
p11_10505604754243
p11_10907225437878
p11_30706805007776
p11_7090849663119119
Table 5. CVRP instance attributes: VRPHAS vs. AASA on SET-4.
Table 5. CVRP instance attributes: VRPHAS vs. AASA on SET-4.
InstanceCVRP SizeNumber of Routes
VRPHASAASAVRPHASAASA
eil22694744
eil23744433
eil301105744
eil331087244
eil5115713966
eilA762282151010
eilB762542151515
eilC7623120188
eilD7621621577
eilA10134724888
eilB1012652481414
Table 6. Performance of VRPHAS, AASA, and ILS-MIP on SET-1.
Table 6. Performance of VRPHAS, AASA, and ILS-MIP on SET-1.
InstanceBest-KnownVRPHASAASAILS-MIPNo Splitting
SolutionGap (%)Time (s)CostGap (%)Time (s)GapTime Time best  (s)Gap (%)
SD1228.280.000.51228.280.001.60 0.00 25.10–0.005.13
SD2708.280.001.23708.280.002.85 0.00 47.51–0.0012.94
SD3430.400.041.25430.580.042.57 0.05 37.68–0.0011.51
SD4630.620.072.35631.040.073.57 0.07 75.45–0.0014.17
SD51389.941.223.911390.570.055.66 0.05 216.69–0.8015.10
SD6830.860.053.85831.240.055.09 0.04 202.60–92.4215.47
SD73640.000.555.473640.000.007.79 0.00 TL–0.0020.87
SD85068.280.007.645068.280.0010.43 0.00 TL–0.0023.11
SD92042.880.727.502051.060.408.97 0.07 TL–3918717.47
SD102683.730.9513.732704.880.7913.68 0.38 TL–960.4219.23
SD1113,280.000.1520.5913,480.001.5018.43 0.00 TL–0.0026.50
SD127213.620.7320.897238.880.3420.99 0.09 TL–1223.3121.99
SD1310,105.860.0525.9010,110.580.0522.05 0.05 TL–1091.8423.49
SD1410,717.530.7042.5710,804.840.8132.03 0.22 TL–967.3023.16
SD1515,094.480.9359.4315,232.600.9147.32 0.37 TL–318.4424.01
SD163379.335.5054.883508.163.8130.16 0.11 TL–0.0127.83
SD1726,493.560.2566.9826,962.641.5650.14 0.1 TL–1328.0526.82
SD1814,202.530.6972.0414,278.700.5454.75 0.46 TL–1347.0523.92
SD1919,995.691.00101.0520,197.041.0070.14 0.64 TL–1328.1724.82
SD2039,635.510.17162.9040,236.661.51104.51 0.25 TL–1328.4427.15
SD2111,271.063.41240.3611,576.962.71129.31 0.2 TL–0.0727.76
Average 0.8286.14 0.7630.56 0.15 992.38–494.1720.60
Table 7. Performance of VRPHAS, AASA, and the no-splitting methodologies on SET-2.
Table 7. Performance of VRPHAS, AASA, and the no-splitting methodologies on SET-2.
InstanceBest-KnownVRPHASAASANo Splitting
SolutionGap (%)Time (s)CostGap (%)Time (s)Gap (%)
S51D1459.500.0010.98459.500.007.920.00
S51D2709.290.7715.04713.800.6311.221.71
S51D3948.061.8820.28959.921.2513.043.19
S51D41562.011.9428.241583.881.4022.947.32
S51D51333.671.0330.681343.1090.7120.818.89
S51D62169.101.8643.812211.901.9723.3710.75
S76D1598.940.0023.18598.940.0014.550.78
S76D21087.401.8232.481104.881.6122.432.79
S76D31427.861.9346.991436.990.6428.052.52
S76D42079.762.5155.892120.161.9440.095.41
S101D1726.591.2334.03729.720.4342.311.14
S101D2 1378.432.0052.441399.671.5434.451.67
S101D31874.811.8570.641905.431.6341.054.13
S101D52791.221.6195.712844.511.9162.5411.59
Average 1.3939.94 1.0826.124.42
Table 8. Performance of VRPHAS, AASA, and the no-splitting methodologies on SET-3.
Table 8. Performance of VRPHAS, AASA, and the no-splitting methodologies on SET-3.
InstanceBest-KnownVRPHASAASANo Splitting
SolutionGap (%)Time (s)CostGap (%)Time (s)Gap (%)
p01_110459.500.009.84459.500.006.630.00
p01_1030757.171.7417.2770.931.8111.083.29
p01_10501005.752.5622.521021.181.5315.764.08
p01_10901488.581.4832.701507.141.2419.6112.75
p01_30701481.712.2332.841503.661.4817.0213.81
p01_70902156.142.2449.822203.932.2221.4111.41
p02_110617.850.6621.82618.130.0412.080.65
p02_10301109.621.6528.861125.651.4421.832.50
p02_10501502.051.0622.521508.980.4626.212.72
p02_10902298.582.0061.012334.331.5541.199.21
p02_30702219.972.2358.802269.502.2328.0212.87
p02_70903223.41.5789.763305.532.5440.5212.64
p03_110752.621.7855.12762.541.3122.081.51
p03_10301458.461.8737.261476.861.2635.672.09
p03_10501996.762.5458.402032.251.7741.763.93
p03_10903085.692.29143.323143.691.8856.5612.17
p03_30702989.32.19136.103034.281.5040.0216.50
p03_70904387.321.66186.144464.651.7684.4813.72
p04_110919.171.35116.08925.190.6533.830.84
p04_10302016.972.71152.442056.711.9755.372.97
p04_10502849.662.40196.822889.621.4072.145.09
p04_10904545.462.44310.844624.181.73112.4114.72
p04_30704334.712.87272.444407.021.6799.6519.64
p04_70906395.411.94469.626555.472.582.0715.09
p05_1101074.182.74186.661086.271.1266.51.42
p05_10302478.41.83241.442515.451.4982.722.46
p05_10503471.411.54311.943510.991.1487.564.01
p05_10905521.572.02561.545619.521.77145.0212.67
p05_30705409.762.41383.285412.850.06142.8115.57
p05_70908192.031.96884.688332.361.71255.0517.28
p11_1101031.112.1174.121044.441.2920.201.55
p11_10302881.82.8994.062923.861.4640.641.95
p11_10504219.012.29108.544268.791.1855.272.71
p11_10906854.093.27188.306994.562.0564.0613.22
p11_30706671.041.44162.366759.991.3359.0115.79
p11_709010,204.811.99231.6010,406.251.97101.0219.93
Average 1.96166.94 1.4533.388.08
Table 9. Performance of VRPHAS, AASA, and the no-splitting methodologies on SET-4.
Table 9. Performance of VRPHAS, AASA, and the no-splitting methodologies on SET-4.
InstanceBest-KnownVRPHASAASANo Splitting
SolutionGap (%)Time (s)CostGap (%)Time (s)Gap (%)
eil22375.280.003.28375.280.003.110.00
eil23568.560.005.92568.560.002.920.00
eil30497.531.506.89505.011.503.641.48
eil33826.411.366.89837.671.365.081.36
eil51524.610.0012.58524.610.0010.950.00
eilA76823.890.4324.60830.970.8619.771.61
eilB761009.041.9631.861025.241.6018.651.86
eilC76738.670.5226.78742.490.5517.480.64
eilD76684.530.8123.82689.470.7219.410.82
eilA101812.512.3555.39828.982.0223.371.92
eilB1011076.261.5730.901087.701.0623.761.53
Average 0.9620.80 0.8813.631.02
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

Torkzaban, N.; Gholami, A.; Baras, J.S.; Golden, B.L. AASA: A Priori Adaptive Splitting Algorithm for the Split Delivery Vehicle Routing Problem. Algorithms 2024, 17, 396. https://doi.org/10.3390/a17090396

AMA Style

Torkzaban N, Gholami A, Baras JS, Golden BL. AASA: A Priori Adaptive Splitting Algorithm for the Split Delivery Vehicle Routing Problem. Algorithms. 2024; 17(9):396. https://doi.org/10.3390/a17090396

Chicago/Turabian Style

Torkzaban, Nariman, Anousheh Gholami, John S. Baras, and Bruce L. Golden. 2024. "AASA: A Priori Adaptive Splitting Algorithm for the Split Delivery Vehicle Routing Problem" Algorithms 17, no. 9: 396. https://doi.org/10.3390/a17090396

APA Style

Torkzaban, N., Gholami, A., Baras, J. S., & Golden, B. L. (2024). AASA: A Priori Adaptive Splitting Algorithm for the Split Delivery Vehicle Routing Problem. Algorithms, 17(9), 396. https://doi.org/10.3390/a17090396

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