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
. 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.