Next Article in Journal
Computing the Stationary Distribution of Queueing Systems with Random Resource Requirements via Fast Fourier Transform
Next Article in Special Issue
An Advanced Learning-Based Multiple Model Control Supervisor for Pumping Stations in a Smart Water Distribution System
Previous Article in Journal
On Periodic Points of the Order of Appearance in the Fibonacci Sequence
Previous Article in Special Issue
On the Resilience of Ant Algorithms. Experiment with Adapted MMAS on TSP
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Selective Vehicle Routing Problem

Department of Mathematics and Computer Science, Technical University of Cluj-Napoca, North University Center of Baia Mare, 430083 Baia Mare, Romania
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(5), 771; https://doi.org/10.3390/math8050771
Submission received: 31 March 2020 / Revised: 5 May 2020 / Accepted: 6 May 2020 / Published: 12 May 2020
(This article belongs to the Special Issue Computational Intelligence)

Abstract

:
The Generalized Vehicle Routing Problem (GVRP) is an extension of the classical Vehicle Routing Problem (VRP), in which we are looking for an optimal set of delivery or collection routes from a given depot to a number of customers divided into predefined, mutually exclusive, and exhaustive clusters, visiting exactly one customer from each cluster and fulfilling the capacity restrictions. This paper deals with a more generic version of the GVRP, introduced recently and called Selective Vehicle Routing Problem (SVRP). This problem generalizes the GVRP in the sense that the customers are divided into clusters, but they may belong to one or more clusters. The aim of this work is to describe a novel mixed integer programming based mathematical model of the SVRP. To validate the consistency of the novel mathematical model, a comparison between the proposed model and the existing models from literature is performed, on the existing benchmark instances for SVRP and on a set of additional benchmark instances used in the case of GVRP and adapted for SVRP. The proposed model showed better results against the existing models.

1. Introduction

The Vehicle Routing Problem (VRP) is one of the most well-known and famous combinatorial optimization problems, being recorded in the literature for a long time. It was in 1959 that Dantzig and Ramser [1] first considered it when they investigated a fuel distribution problem. In this context, they provided a near optimal solution based on a linear programming formulation for the most economical collection of routes of a fleet of gasoline delivery trucks between a bulk terminal and a given number of service gas stations. There are several basic aspects that justify the significance and importance of the VRP: we can mention effective solution approaches which make it practical to solve when dealing with large graphs (see the savings algorithm developed by Clarke and Wright [2]); the problem also has numerous real-life applications in industry, especially in the areas of transportation networks and supply chain management. Finding the optimal solution to VRP is NP-hard; therefore, only problems with limited size can be solved, optimally, so a lot research is developed in the area of heuristic and metaheuristic algorithms that can provide near optimal solutions within reasonable computational times for real world VRPs.
Several variations of the VRP exist in the literature: the Capacitated Vehicle Routing Problem [3], the Vehicle Routing Problem with Pickup and Delivery [4], the Vehicle Routing Problem with Time Windows [5], the Vehicle Routing Problem with Multiple Trips [6], the Open Vehicle Routing Problem [7], the Green Vehicle Routing Problem [8], the Generalized Vehicle Routing Problem [9,10,11,12,13], the Clustered Vehicle Routing Problem [14,15,16], etc.
Recently, Posada et al. [17] considered a novel variant of the VRP, called the Selective Vehicle Routing Problem (SVRP), justified by some important real applications of the problem such as response operations in humanitarian logistics in the case of disasters where relief aid should be given to any node covering a disaster area as soon as possible, design of urban transportation routes, etc. The proposed variant extends the generalized vehicle routing problem in the sense that the nodes of the graph which are divided into clusters may belong to one or more clusters. Posada et al. [17] provided a mixed integer programming formulation of the SVRP, and in an alternative way, they have rewritten some of the constraints of the proposed model using different variants of the restrictions resulting three new mathematical formulations.
Evidently, SVRP belongs to the class of generalized combinatorial optimization problems. This category of problems naturally generalizes the classical optimization problem, having the following primary features: the vertices of the underlying graph are divided into a certain number of clusters and, when considering the feasibility constraints of the initial problem, these are expressed in relation to the clusters rather than as individual vertices. Some of the most investigated generalized combinatorial optimization problems are: the generalized minimum spanning tree problem and its variants, the generalized traveling salesman problem and its variants [18,19], the generalized vehicle routing problem and its variants [20,21], the selective graph coloring problem [22,23], the clustered shortest path tree problem, etc. For further reference on the class of generalized combinatorial optimization problems, we refer to [24,25]. A closely related problem to SVRP was introduced by Pop [18] and it was called the Selective Minimum Spanning Tree Problem. In this case, the vertices were grouped within clusters such that each vertex may belong to one ore more clusters and the objective was to find a minimum cost tree spanning a subset of vertices that includes exactly one vertex from each cluster.
The present paper is organized as follows. Section 2 provides a formal definition of the selective vehicle routing problem. Section 3 presents the novel mixed integer programming based mathematical model of the investigated problem. In Section 4, we validate the consistency of the proposed model by comparing it against the existing models from the literature on two sets of instances. The concluding results and as well as further research directions are presented.

2. Definition of the Selective Vehicle Routing Problem

Let G = V , E be an undirected graph with the set of nodes V = { 0 , 1 , 2 , , n } , where the node 0 represents the depot, and the set of edges E , E { { i , j } | i , j V , i j } .
The set of nodes V is divided into m + 1 nonempty subsets denoted C 0 , C 1 , , C m and called clusters, where the cluster C 0 consists of only one node namely the depot, i.e., C 0 = { 0 } and the following condition holds
V = C 0 C 1 C 2 C m .
We denote by C the set of all clusters and C = 1 , 2 , , m . Dividing the nodes of the graph in this way, we do not exclude the possibility that certain nodes to belong to more than one cluster, i.e., there are C c , C d C such that C c C d .
We associate for each node i V 0 a certain demand denoted by q i , i 1 , 2 , , n and a metric coordinate denoted by ( x , y ) . Each edge ( i , j ) E has associated a certain cost denoted by d i j and representing the Euclidean distance between the nodes i and j, with i , j V . The depot has no demand, i.e., q 0 = 0 . The total demand of each cluster is obtained by summing up the demands of all the nodes belonging to that cluster and may be satisfied by any of its nodes. In addition, we assume that at the depot is based a fleet of K of identical vehicles, each with a capacity denoted by Q and K = 1 , 2 , , K .
The goal of the Selective Vehicle Routing Problem is to determine a collection of routes, one for each vehicle, with minimum cost such that each route starts and ends at the depot, each cluster should be visited exactly once, the entering and leaving node of each cluster is the same, and the sum of all the demands of any route does not exceed the capacity of the vehicle Q. Therefore, the SVRP involves the following two related decisions:
  • Select a subset of nodes S V such that S C c = 1 , for all c { 1 , , m } = C .
  • Determine a minimum cost collection of routes in the subgraph of G induced by S, satisfying the capacity constraints.
An illustrative scheme of the SVRP and two feasible collections of routes are shown in Figure 1.
In Figure 1, we present two feasible solutions of the SVRP. We can observe that the total number of visited nodes may be different: six nodes are visited in the first feasible solution and eight in the second solution. If a belonging to the three clusters C 6 , C 7 , and C 8 is visited, then all the nodes in the clusters C 6 , C 7 , and C 8 are covered.

3. A Valid Mathematical Model of Selective Vehicle Routing Problem

With the purpose of formulating the SVRP as a mixed integer linear program (MILP), we introduce the following parameters and decision variables:
  • The parameter Λ = λ i c 1 i n 1 c m indicates if a certain node i V belongs to a cluster C c , c { 1 , , m } , thus
    λ i c = 1 , if i C c 0 , if i C c .
    This means, if the node i is an element of the cluster C c , then λ i c = 1 ; otherwise, λ i c = 0 .
  • The binary decision variable x i j k is defined as follows:
    x i j k = 1 , if the vehicle k traverse the edge i , j 0 , if the edge i , j is not part of any route .
  • The integer decision variables u i j are defined as follows: u i j = 0 if the edge i , j is not included in a solution route, and may take a positive integer value less or equal to n, showing its order in the corresponding route.
Then, we can express the Selective Vehicle Routing Problem as the following mixed integer programming problem:
min i = 0 n j = 0 n κ = 1 k d i j · x i j k
s . t . i = 0 i j n x i j k x j i k = 0 , j V , k K
i = 0 n j = 0 i j n x i j k 1 , k K
i = 0 n j = 0 i j n κ = 1 k λ j c · x i j k = 1 , c C 0
j = 1 n k = 1 K x 0 j k = K
i = 0 n j = 0 i j n q i · x i j k Q , k K
n · k = 1 K x i j k u i j , i , j V , i j
j = 0 j i n u i j u j i = j = 0 j i n k = 1 K x i j k , i V
x i j k { 0 , 1 } , i , j V , k K
u i j 0 , i , j V .
Our objective is to minimize the total traveled distance. The restrictions in Equation (1) guarantee that if a vehicle reaches a certain node then it has to leave it. The restrictions in Equation (2) ensure that for every vehicle there is allocated at least one arc. The restrictions in Equation (3) guarantee that each cluster is visited exactly once. The restrictions in Equation (4) force using all the existing vehicles and in addition all the routes start and finish at the depot. The restrictions in Equation (5) provides the conditions imposed by the capacities of the vehicles. The constraints in Equations (6) and (7) are the sub-tour elimination constraints. Finally, the last constraints set the ranges of the decision variables.

4. Computational Results

To solve the proposed mathematical formulation of the SVRP, we used CPLEX 12.9. The computational experiments were performed on a PC with Procesor AMD Ryzen 9 12-Core, 64 GB RAM, Windows 10 operating system.
We report two sets of experiments: the first one compares our achieved results against the results obtained by Posada et al. [17] on a set of 36 SVRP instances generated by the same authors and the second one on a new set of instances containing 120 instances adapted from GVRP instances.

4.1. Problem Instances

Posada et al. [17] used a set of 12 GVRP instances with up to 63 nodes and adapted them in the case of SVRP by considering three levels of complexity. In their approach, the level of complexity is given by the maximum number of clusters to whom a node may belong and they considered three cases: nodes that may belong to two, three, and four clusters. This results in 12 × 3 = 36 SVRP instances. The instances were obtained by adding some random nodes on the predefined clusters. The proposed procedure by Posada et al. [17] has the disadvantage that it produces occasionally unfeasible solutions, which are repaired manually.
Next, we describe a new procedure for adapting the GVRP instances in the case of SVRP. In our case, the number of nodes and clusters and the distances between the nodes are identical for both GVRP and SVRP. We propose a natural algorithm that allows nodes close to each other to belong to more than one cluster. We identify these nodes based on the distance between them. We define the following parameters:
  • We considered four cases for the maximum number of clusters that a node may belong. If each node belongs to exactly one cluster, that corresponds to GVRP, and the other cases in which nodes may belong to two, three, or four clusters are SVRP.
  • The maximum percentage of nodes within a cluster that may belong to two clusters was established as 50%. This value can be increase based on chosen nodes that can be in three or four clusters.
  • The maximum percentage of nodes that may belong to two clusters was set to 30%. This value can also be increase based on chosen nodes that can be in three or four clusters.
We developed a natural Greedy algorithm that allocates the nodes of two clusters as follows: we sort in ascending order the edges depending on their costs; the edges containing the depot are eliminated. We take the edge with the lowest cost, and, if its nodes are from different clusters, then the node with the left side of the edge is assigned to the cluster to which the other node belongs. The corresponding modification is made in the lambda matrix. We keep considering edges until the initial constraints are met, i.e., the maximum percentage of nodes within a cluster that may belong to several clusters is 50% and the maximum percentage of nodes that may belong to more that one cluster is 30%.
To choose nodes that will belong to three clusters, we take edges in the same order as above. If the node from the left side is contained in two clusters, we search the next edges in the existing order by cost for an edge that has on the left side the same node as the chosen one and on the right side a node that is contained in another cluster that the left side node. The chosen node is attached to the cluster of the node identified on the right side of the edge. We repeat this step until the number of nodes that are contained by three clusters represents a quarter of the maximum nodes that can belong to more clusters.
In the same way, we choose nodes that are included in four clusters based on nodes that belong to three clusters, with the difference that the maximum number of this type of nodes one eighth of the maximum nodes that can belong to more clusters.
This type of algorithm, where the steps are repeated for each level of complexity, has been chosen to ensure that the results can be obtained by other researchers using different instances from one part and from another to ensure that we have a reasonable distribution of nodes for each level of complexity.
The new proposed instances of SVRP are accessible at: https://data.mendeley.com/datasets/k3vw6mm7tn.

4.2. Computational Results on the Possada et al. Instances

In Table 1, we display the achieved computational results performed for the instances generated by Posada et al. [17] (36 SVRP instances and 12 GVRP instances) and we compare them with the results provided by the same authors. The first column displays the mathematical model: the models MILP 1 , MILP 2 , MILP 3 , and MILP 4 were introduced by Posada et al. [17] and the last model denoted by MILP is the formulation that we provide in this paper. The next three columns provide the average gap, average time, and the number of achieved optimal solutions by solving the SVRP instances by each of the considered models. The last three columns provide the average gap, average time, and the number of achieved optimal solutions by solving the GVRP instances by each of the considered models.
Analyzing the results displayed in Table 1, one can notice that our proposed model outperforms the existing models, providing lower average gaps, the shortest running time, and the largest number of achieved optimal solutions. In the case of the GVRP instances, we can observe that, if we allowed the computer to run for more than 1800 s, we were able to obtain 11 optimal solutions out of 12 instances and the obtained average gap was 0.59%.

4.3. Computational Results on Our Proposed Instances

In Table 2, Table 3 and Table 4, we display the achieved computational results performed for the novel generated instances: 24 GVRP instances, 72 small and medium SVRP instances, and 48 larger SVRP instances. We considered 24 GVRP instances (for more information, see Bektas et al. [26]), with up to 39 nodes. These instances were adapted in the case of SVRP by considering as well four levels of complexity. The level of complexity is given by the maximum number of clusters to which a node may belong: nodes may belong to two, three, and four clusters. Resulting in this way, there were 24 × 3 = 72 small and medium SVRP instances and 16 × 3 = 48 larger SVRP instances. The level of complexity 1 corresponds to GVRP, in which case every node belongs to exactly one cluster.
In Table 2, we present the obtained computational results in the case of SVRP instances with the level of complexity 1, which coincides to GVRP instances. The first column displays the index of the instance. The second column contains the name of the instance and has the following format S V R P X n Y k Z C m V K , where X corresponds to the type of the instance used in GVRP, Y refers to the number of nodes, Z corresponds to the number of vehicles in the original CVRP instance, m is the number of clusters, and K is the number of vehicles in the GVRP instance. The next column represents the level of complexity of the adapted SVRP instance, in this case being 1 corresponds to GVRP. The next column delivers the optimal solution obtained by Bektas et al. [26] for solving the GVRP. The last two columns display our obtained results for solving the GVRP instance using the proposed mathematical model: the achieved solution and the necessary computational time in seconds.
Analyzing the results presented in Table 2, we can observe that for all 24 GVRP instances which correspond to SVRP instance with level of complexity, we were able to obtain the optimal solution of the problem. The computational time was below 3600 s in 21 out of 24 instances.
In Table 3 and Table 4, we present the obtained computational results in the case of small and medium and larger SVRP instances, respectively, with the levels of complexity 2, 3 and 4. The first column displays the index of the instance and the second column contains the name of the instance. The next columns display the level of complexity of the adapted SVRP instance, the optimal solution obtained by solving the problem with CPLEX, and the corresponding computational time in seconds.
Analyzing the results presented in Table 3 and Table 4, we can observe that, in 62 out of 72 small and medium SVRP instances and 34 out of 48 larger SVRP instances, corresponding to the three levels of complexity 2, 3, and 4, we were able to provide the optimal solution of the problem. For the other instances, the corresponding problems were unfeasible because, considering nodes that belong to more clusters, the capacity constraints of the vehicles to serve those clusters are no longer satisfied, and in the case of Table 4 in 12 out of 48 larger SVRP instances CPLEX was stoped after 3600 seconds.

5. Conclusions

In this paper, we describe a novel mathematical formulation based on mixed integer programming of a recently introduced variant of VRP, namely the selective vehicle routing problem (SVRP). The proposed model is compared against the existing models from the literature on the benchmark instances introduced by Posada et al. [17]. In addition, we provide a new set of SVRP instances adapted from the GVRP instances. Our proposed model achieved better results against the existing models from the literature.
Due to the complexity of the SVRP, in the future, we plan to solve the problem with heuristic or metaheuristic algorithms.

Author Contributions

Conceptualization, P.C.P. and C.S.; methodology, P.C.P.; software, C.S.; validation, P.C.P., C.S., and A.H.-M.; formal analysis, P.C.P., C.S., and A.H.-M.; investigation, P.C.P.; resources, P.C.P.; data curation, C.S.; writing—original draft preparation, P.C.P. and A.H.-M.; writing—review and editing, P.C.P.; visualization, P.C.P.; supervision, P.C.P.; project administration, C.S.; and funding acquisition, C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This paper was supported by the Post-Doctoral Programme proiectul “Entrepreneurial competences and excellence research in doctoral and postdoctoral programs—ANTREDOC”, project co-funded from European Social Fund.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Dantzig, G.B.; Ramser, J.H. The Truck Dispatching Problem. Manag. Sci. 1959, 6, 80–91. [Google Scholar] [CrossRef]
  2. Clarke, G.; Wright, J.W. Scheduling of Vehicles from a Central Depot to a Number of Delivery Points. Oper. Res. 1964, 12, 568–581. [Google Scholar] [CrossRef]
  3. Letchford, A.; Salazar-Gonzales, J. The capacitated vehicle routing problem: Stronger bounds in pseudo-polynomial time. Eur. J. Oper. Res. 2019, 227, 24–31. [Google Scholar] [CrossRef] [Green Version]
  4. Fan, J. The vehicle routing problem with simultaneous pickup and delivery based on customer satisfaction. Procedia Eng. 2011, 15, 5284–5289. [Google Scholar] [CrossRef] [Green Version]
  5. El-Sherbeny, N. Vehicle routing with time windows: An overview of exact, heuristic and metaheuristic methods. J. King Saud Univ.-Sci. 2010, 22, 123–131. [Google Scholar] [CrossRef] [Green Version]
  6. Cattaruzza, D.; Absi, N.; Feillet, D. Vehicle routing problems with multiple trips. 4OR 2016, 14, 223–259. [Google Scholar] [CrossRef]
  7. Li, F.; Golden, B.; Wasil, E. The open vehicle routing problem: Algorithms, large-scale test problems and computational results. Comput. Oper. Res. 2007, 34, 2918–2930. [Google Scholar] [CrossRef]
  8. Erdoğan, S.; Miller-Hooks, E. A Green Vehicle Routing Problem. Transp. Res. Part Logist. Transp. Rev. 2012, 48, 100–114. [Google Scholar] [CrossRef]
  9. Ghiani, G.; Improta, G. An efficient transformation of the generalized vehicle routing problem. Eur. J. Oper. Res. 2000, 122, 11–17. [Google Scholar] [CrossRef]
  10. Pintea, C.; Chira, C.; Dumitrescu, D.; Pop, P. Sensitive ants in solving the generalized vehicle routing problem. Int. J. Comput. 2011, 6, 734–741. [Google Scholar] [CrossRef] [Green Version]
  11. Pop, P.C.; Matei, O.; Sitar, C.P. An improved hybrid algorithm for solving the generalized vehicle routing problem. Neurocomputing 2013, 109, 76–83. [Google Scholar] [CrossRef]
  12. Pop, P.C.; Kara, I.; Horvat-Marc, A. New mathematical models of the generalized vehicle routing problem and extensions. Appl. Math. Model. 2012, 36, 97–107. [Google Scholar] [CrossRef]
  13. Horvat-Marc, A.; Fuksz, L.; Pop, P.C.; Dănciulescu, D. A decomposition-based method for solving the clustered vehicle routing problem. Log. J. IGPL 2017, 26, 83–95. [Google Scholar] [CrossRef]
  14. Defryn, C.; Sörensen, K. A fast two-level variable neighborhood search for the clustered vehicle routing problem. Comput. Oper. Res. 2017, 83, 78–94. [Google Scholar] [CrossRef] [Green Version]
  15. Exposito-Izquierdo, C.; Rossi, A.; Sevaux, M. A two-level solution approach to solve the clustered vehicle routing problem. Comput. Ind. Eng. 2016, 91, 274–289. [Google Scholar] [CrossRef]
  16. Pop, P.C.; Fuksz, L.; Horvat-Marc, A.; Sabo, C. A novel two-level optimization approach for clustered vehicle routing problem. Comput. Ind. Eng. 2018, 115, 304–318. [Google Scholar] [CrossRef]
  17. Posada, A.; Rivera, J.; Palacio, J. A Mixed-Integer Linear Programming Model for a Selective Vehicle Routing Problem. In Proceedings of the Applied Computer Sciences in Engineering (WEA 2018), Medellín, Colombia, 17–19 October 2018; Figueroa-García, J.C., Villegas, J.G., Orozco-Arroyave, J.R., Eds.; Communications in Computer and Information Science. Springer: Berlin/Heidelberg, Germany, 2018; Volume 916, pp. 108–119. [Google Scholar] [CrossRef]
  18. Pop, P. The generalized minimum spanning tree problem: An overview of formulations, solution procedures and latest advances. Eur. J. Oper. Res. 2020, 283, 1–15. [Google Scholar] [CrossRef]
  19. Pop, P.C.; Matei, O.; Sabo, C.; Petrovan, A. A two-level solution approach for solving the generalized minimum spanning tree problem. Eur. J. Oper. Res. 2018, 265, 478–487. [Google Scholar] [CrossRef]
  20. Fischetti, M.; González, J.J.S.; Toth, P. The symmetric generalized traveling salesman polytope. Networks 1995, 26, 113–123. [Google Scholar] [CrossRef]
  21. Fischetti, M.; González, J.J.S.; Toth, P. A Branch-and-Cut Algorithm for the Symmetric Generalized Traveling Salesman Problem. Oper. Res. 1997, 45, 378–394. [Google Scholar] [CrossRef]
  22. Demange, M.; Monnot, J.; Pop, P.; Ries, B. On the complexity of the selective graph coloring problem in some special classes of graphs. Theor. Comput. Sci. 2014, 540–541, 89–102. [Google Scholar] [CrossRef]
  23. Fidanova, S.; Pop, P. An improved hybrid ant-local search for the partition graph coloring problem. J. Comput. Appl. Math. 2016, 293, 55–61. [Google Scholar] [CrossRef]
  24. Feremans, C.; Labbé, M.; Laporte, G. Generalized network design problems. Eur. J. Oper. Res. 2003, 148, 1–13. [Google Scholar] [CrossRef]
  25. Pop, P. Generalized Network Design Problems, Modelling and Optimization; De Gruyter: Berlin, Germany, 2012. [Google Scholar]
  26. Bektaş, T.; Erdoğan, G.; Røpke, S. Formulations and Branch-and-Cut Algorithms for the Generalized Vehicle Routing Problem. Transp. Sci. 2011, 45, 299–316. [Google Scholar] [CrossRef]
Figure 1. An illustrative example of the SVRP and two feasible solutions.
Figure 1. An illustrative example of the SVRP and two feasible solutions.
Mathematics 08 00771 g001
Table 1. Summary of the computational results on Possada et al. [17] instances.
Table 1. Summary of the computational results on Possada et al. [17] instances.
ModelSVRPGVRP
Average GapAverage Time (seconds)Number of Optimal SolutionsAverage GapAverage Time (seconds)Number of Optimal Solutions
MILP 1 21.26%1848.221933.90%2742.763
MILP 2 21.38%1663.402029.19%2505.064
MILP 3 19.51%1602.862229.83%2555.154
MILP 4 21.46%1848.881933.91%2978.263
MILP0.0004%21.05302.26%1738.957
0.59%23,063.5511
Table 2. Computational results achieved by our mathematical model in the case of SVRP instances with level of complexity 1.
Table 2. Computational results achieved by our mathematical model in the case of SVRP instances with level of complexity 1.
Instance IndexInstance’s Name ComplexityLevel of SolutionOPT [26]Our Results
SolutionTime
1.SVRP-A-n32-k5-C11-V21386386574.578
2.SVRP-A-n33-k6-C11-V21370370106.329
3.SVRP-A-n33-k6-C17-V314654655231.750
4.SVRP-A-n34-k5-C12-V2141941933,748.204
5.SVRP-A-n37-k5-C13-V21347347965.516
6.SVRP-A-n37-k5-C19-V31432432335.703
7.SVRP-A-n38-k5-C13-V213673673278.234
8.SVRP-B-n31-k5-C11-V21356356588.157
9.SVRP-B-n38-k6-C13-V2137037021,129.14
10.SVRP-B-n39-k5-C13-V212802801147.344
11.SVRP-P-n16-k8-C6-V411701700.250
12.SVRP-P-n16-k8-C8-V512392390.516
13.SVRP-P-n19-k2-C10-V211471470.422
14.SVRP-P-n19-k2-C7-V111111110.203
15.SVRP-P-n20-k2-C10-V211541540.625
16.SVRP-P-n20-k2-C7-V111171170.187
17.SVRP-P-n21-k2-C11-V211601600.813
18.SVRP-P-n21-k2-C7-V111171170.484
19.SVRP-P-n22-k2-C11-V211621621.078
20.SVRP-P-n22-k2-C8-V111111110.390
21.SVRP-P-n22-k8-C11-V5131431423.250
22.SVRP-P-n22-k8-C8-V412492491.313
23.SVRP-P-n23-k8-C12-V51312312145.422
24.SVRP-P-n23-k8-C8-V3117417430.672
Table 3. Experimental results for solving small and medium SVRP instances with levels of complexity 2, 3 and 4.
Table 3. Experimental results for solving small and medium SVRP instances with levels of complexity 2, 3 and 4.
Instance IndexInstance’s NameLevel of ComplexityOptimal SolutionTime (seconds)Level of ComplexityOptimal SolutionTime (seconds)Level of ComplexityOptimal SolutionTime (seconds)
1.SVRP-A-n32-k5-C11-V22319.930.7813315.670.9374310.140.640
2.SVRP-A-n33-k6-C11-V22285.880.6093285.880.4214285.880.500
3.SVRP-A-n33-k6-C17-V32406.054.5783393.178.9064391.662.297
4.SVRP-A-n34-k5-C12-V22319.930.7813319.200.5784313.430.500
5.SVRP-A-n37-k5-C13-V22314.811.6723316.361.5314313.061.078
6.SVRP-A-n37-k5-C19-V32433.772.6413433.774.3594428.415.437
7.SVRP-A-n38-k5-C13-V22307.451.4383305.230.9064316.500.656
8.SVRP-B-n31-k5-C11-V22334.000.4533330.710.6404335.260.375
9.SVRP-B-n38-k6-C13-V22297.48144.8133300/343.4844293.909.094
10.SVRP-B-n39-k5-C13-V22267.9335.5163267.8116.2034243.8212.609
11.SVRP-P-n16-k8-C6-V42171.480.0323171.480.0314unfeasible
12.SVRP-P-n16-k8-C8-V52241.680.0623241.680.0464249.300.015
13.SVRP-P-n19-k2-C10-V22144.480.1723133.000.1414123.710.094
14.SVRP-P-n19-k2-C7-V1295.430.141395.430.094488.070.078
15.SVRP-P-n20-k2-C10-V22125.750.1883128.650.1724134.350.110
16.SVRP-P-n20-k2-C7-V1297.950.140394.260.156497.950.140
17.SVRP-P-n21-k2-C11-V22134.050.2353134.050.1884131.800.172
18.SVRP-P-n21-k2-C7-V12106.100.140392.020.1094103.150.203
19.SVRP-P-n22-k2-C11-V22unfeasible 3unfeasible 4unfeasible
20.SVRP-P-n22-k2-C8-V1288.090.413377.520.442477.180.526
21.SVRP-P-n22-k8-C11-V52295.280.3293295.280.3284292.490.328
22.SVRP-P-n22-k8-C8-V42213.510.3593213.510.3754208.860.234
23.SVRP-P-n23-k8-C12-V52unfeasible 3unfeasible 4unfeasible
24.SVRP-P-n23-k8-C8-V32unfeasible 3unfeasible 4unfeasible
Table 4. Experimental results for solving larger SVRP with levels of complexity 2, 3 and 4.
Table 4. Experimental results for solving larger SVRP with levels of complexity 2, 3 and 4.
Instance IndexInstance’s NameLevel of ComplexityBest SolutionTime (seconds)Level of ComplexityBest SolutionTime (seconds)Level of ComplexityBest SolutionTime (seconds)
1.A-n39-k5-C13-V22283.103.0473288.018.5944282.661.5
2.A-n39-k6-C13-V22305.671.5313313.510.5944321.751.203
3A-n48-k7-C16-V32355.88110.0473366.65232.0944337.273600 *
4A-n54-k7-C18-V32378.53964.9843360.5526.4224391.42198.984
5A-n55-k9-C19-V32377.20107.3283unfeasible-4unfeasible-
6A-n61-k9-C21-V42381.453600 *3369.651182.54340.2181.86
7A-n65-k9-C22-V32409.573600 *3370.9949.2974391.823600 *
8B-n45-k5-C15-V22378.1144.4383353.132.6724348.444.281
9B-n45-k5-C23-V32453.99146.0473436.3130.4374435.863.031
10B-n50-k8-C17-V32475.57820.0313475.57781.5784457.71250.156
11B-n64-k9-C22-V42429.273600 *3432.943600 *4405.95225.515
12B-n68-k9-C34-V52569.073600 *3632.693600 *4601.063600 *
13P-n50-k10-C17-V42263.0692.8753266.01143.8444253.44201.859
14P-n55-k10-C19-V42254.8192.0783258.3585.6724260.5699.656
15P-n60-k15-C20-V52300.30677.9843300.35235.1874304.88251.922
16P-n76-k5-C26-V22299.873600 *3319.413600 *4288.403600 *
* The time values marked with 3600 * means that the algorithm was interrupted after 3600 s.

Share and Cite

MDPI and ACS Style

Sabo, C.; Pop, P.C.; Horvat-Marc, A. On the Selective Vehicle Routing Problem. Mathematics 2020, 8, 771. https://doi.org/10.3390/math8050771

AMA Style

Sabo C, Pop PC, Horvat-Marc A. On the Selective Vehicle Routing Problem. Mathematics. 2020; 8(5):771. https://doi.org/10.3390/math8050771

Chicago/Turabian Style

Sabo, Cosmin, Petrică C. Pop, and Andrei Horvat-Marc. 2020. "On the Selective Vehicle Routing Problem" Mathematics 8, no. 5: 771. https://doi.org/10.3390/math8050771

APA Style

Sabo, C., Pop, P. C., & Horvat-Marc, A. (2020). On the Selective Vehicle Routing Problem. Mathematics, 8(5), 771. https://doi.org/10.3390/math8050771

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