Learning and Adaptation for Optimization and Control of Complex Renewable Energy Systems
View this Special IssueResearch Article  Open Access
Ning Guo, Bin Qian, Rong Hu, Huai P. Jin, Feng H. Xiang, "A Hybrid Ant Colony Optimization Algorithm for MultiCompartment Vehicle Routing Problem", Complexity, vol. 2020, Article ID 8839526, 14 pages, 2020. https://doi.org/10.1155/2020/8839526
A Hybrid Ant Colony Optimization Algorithm for MultiCompartment Vehicle Routing Problem
Abstract
The multicompartment vehicle routing problem (MCVRP) has been applied in fuel or food delivery, waste collection, and livestock transportation. Ant colony optimization algorithm (ACO) has been recognized as an efficient method to solve the VRP and its variants. In this paper, an improved hybrid ant colony optimization algorithm (IHACO) is proposed to minimize the total mileage of the MCVRP. First, a probabilistic model is designed to guide the algorithm search towards highquality regions or solutions by considering both similar blocks of customers and customer permutations. Then, a heuristic rule is presented to generate initial individuals to initialize the probabilistic model, which can drive the search to the highquality regions faster. Moreover, a new local search using the geometry optimization is developed to execute exploitation from the promising regions. Finally, two types of variable neighborhood descent (VND) techniques based on the speedup search strategy and the first move strategy are devised to further enhance the local exploitation ability. Comparative numerical experiments with other algorithms and statistical analyses are carried out, and the results show that IHACO can achieve better solutions.
1. Introduction
The multicompartment vehicle routing problem (MCVRP) is a variant of the wellknown vehicle routing problem (VRP) [1]. Compared with the standard VRP, the MCVRP can simultaneously transport different types of products, which have to be placed in separate compartments. Since the customers probably have different demands for various products whose individual characteristics are so different that they cannot be mixed, each vehicle is divided into several compartments with a certain capacity for satisfying customers’ needs. There are many industrial fields in which the multicompartment vehicles are employed. The first example of the MCVRP arises in the process of supplying the fuels, where a lot of vehicles or ships with some tanks of various capacities are used to settle this problem [1–7]. Another example is the transportation of food, where different degrees of refrigeration goods are stored in different compartments in one vehicle [8–11]. In addition, the MCVRP has also been widely applied in the waste collection, where separate compartments are required to convey different kinds of garbage [12–14]. Despite the practical applications, the solution of MCVRP has been rarely studied. Thus, this paper aims to minimize the total mileage for such a problem.
Since the VRP is an NPhard problem and it can be reduced to the MCVRP, the latter issue is also NPhard. Some operational research algorithms were presented to resolve smallscale problems for the MCVRP [1, 15–20]. However, with the increasing problem size, these algorithms are usually of limited use owing to their excessive timeconsuming nature. Then, a variety of metaheuristic algorithms were proposed to obtain desirable solutions within a reasonable calculation time [21–24]. There are many intelligent algorithms for solving nonlinear problems [25–27]. The ant colony optimization algorithm (ACO) is one of the highly efficient metaheuristic or intelligent algorithms. The ACO is inspired by the intelligent behavior of ants in searching for food, during which the ants leave the socalled pheromone on the road. The ants exchange information with other ants through the pheromone. The higher the number of pheromones on a path, the more likely the ants are to go after that path, and vice versa. In this way of the pheromone positive feedback mechanism, the ants can find the shortest path to the food location.
Because the ACO has a strong global exploration and distributed computing capability, it has been successfully applied to handle different kinds of optimization problems. While solving the VRP and the complex problems derived from it, ACO has been usually adopted by many researchers. The existing results usually use only pheromone concentration matrix (hereinafter referred to as pheromone matrix) to construct a probabilistic model to generate new individuals. The pheromone matrix learns and accumulates similar blocks of customers from every outstanding individual. Here any two successive visited customers in an individual compose one block. Mavrovouniotis and Yang [28] proposed an ACO with immigrant schemes for the dynamic vehicle routing problem. In this ACO, only the pheromone matrix was utilized to get information of each block, and its initial value was set to , where is the corresponding path length of a solution constructed by a nearest neighbor heuristic. Schyns [29] presented an ACO to settle the responsive dynamic vehicle routing problem with time window. The pheromone matrix was adopted merely and initialized to , where is the number of customers and is the length of the route formed on the basis of the FIFO (first in first out) principle. Besides, a local search basic heuristic was developed to improve the algorithm performance. Wang et al. [30] designed a novel ACO called AMR to solve the vehicle routing problem. In AMR, only the pheromone matrix was used and its initial value was a small positive constant denoted by . In addition, the saving algorithm was introduced to enhance the performance of AMR. Huang et al. [31] employed an ACO combined with a local search method in the resolution of the routing problem with bridge inspection tasks. Only the pheromone matrix was applied, and the initial pheromone was set to , where is the total number of nodes and is the cost of the initial solution obtained by the first ant. Moreover, ACO has been rarely used to solve the MCVRP. Reed et al. [21] added a technique of clustering and local search of 2opt to ACO in their study. Only the pheromone matrix was utilized and initialized to , where is the customers’ number and is the total length of the nearest neighbor path. Motivated by the work of Reed et al. [21], Abdulkader et al. [22] proposed a hybridized ant colony (HAC) algorithm with local search schemes. The pheromone matrix was adopted alone, and the initial value was set to , where L is the length of a path which was generated randomly.
As it can be seen from above literatures, the probabilistic model only consists of the pheromone matrix in which only the relationship between two consecutive customers is considered. However, as mentioned earlier, the ACO generates new individuals based on the probabilistic model which stores the information of excellent individuals. Obviously, the more comprehensive the information of excellent individuals are learned, the easier it is to guide the algorithm to reach different highquality regions as soon as possible. Hence, only the information accumulation of similar blocks in the customer sequence is incomplete and limited. Besides, although pheromone matrices in different literatures above are initialized to different values, the values in the pheromone matrix of each literature are the same. This means that, at the beginning of the algorithm, there is no valuable information in the pheromone matrix. With the iteration of the algorithm, when the information in the pheromone matrix accumulates to a certain amount, it will guide new individuals towards highquality areas. Invisibly, the efficiency of the algorithm is reduced in this way. Moreover, it is also found that the introduction of local search on the basis of the effective neighborhood operation can further obtain higher performance of the algorithm. Hence, in this research, an improved hybrid ant colony optimization algorithm (IHACO) is devised to solve the MCVRP. The salient features of this IHACO are as follows:(1)A new matrix is designed to save the whole order information of customers. Together with the pheromone matrix, the probabilistic model of IHACO is formed to record the worthy information of optimal individuals. Compared with the previous work, this design can learn and retain valuable information more comprehensively.(2)Several initial individuals produced by a heuristic rule are used to initialize the probabilistic model so that some valuable information can be utilized at the beginning of the algorithm. Then, the speed of reaching highquality solution regions will be accelerated, and thus, the efficiency of the algorithm is raised.(3)Different from the majority of current local searches where 2opt is commonly used within the same route, a new local search is proposed by adopting the geometry optimization method to perform narrow but productive exploitation.(4)In the IHACO’s local search, three techniques are devised to enhance the exploitation capability from many aspects: two kinds of VND are used to increase the search depth of the algorithm, a speedup search strategy based on the problem’s properties is used to accelerate the evaluation process of individuals, and a first move strategy is used to enrich search regions.
The remainder of this paper is organized as follows: Section 2 introduces the mathematical model for the MCVRP; Section 3 describes the proposed IHACO for details; Section 4 presents computational experiments and analyses; and finally, Section 5 gives some conclusions and suggestions of future work.
2. Problem Description and Model Formulation
The MCVRP is defined on an undirected graph with a set of vertices and a set of arcs . Denote as the distance from to . Node in represents the depot, and are nodes of customers. The depot has an available homogeneous fleet. These identical vehicles start from the depot, visit the customers distributed in the service area, and go back to the depot in the end. All vehicles are equipped with multiple compartments whose number equals to the number of products. Each customer has the quantity to be carried for the product , and the compartment capacity reserved for this product does not exceed . Moreover, each customer is visited exactly once by only one vehicle.
To formulate the model, some sets, parameters, and variables adopted in the proposed problem are listed in Table 1.

Then, the integer linear program of the MCVRP can be described as follows:subject to
Equation (1) is the objective function indicating the total mileage for the problem. Constraints (2) and (3) imply that one customer is visited just once by only one vehicle. Constraint (4) guarantees that each vehicle starts and ends its trip at the depot. Constraint (5) ensures that the total quantity delivered by each vehicle after serving customer is less than its capacity for this product. Constraint (6) imposes the connectivity requirement of the feasible routes. Both inequalities (5) and (6) are subtour elimination constraints. Constraint (7) describes the domains of the decision variable.
3. IHACO for the MCVR
This section describes the implementation of the proposed IHACO in detail, including the solution representation, establishment, and initialization of the probabilistic model, formation of the individuals, update of probability matrices, speedup VNDbased local search, and the complete IHACO process.
3.1. Solution Representation
Usually, the encoding rules of ACO for solving VRP are that each ant represents one vehicle, the path of each ant represents the driving path of each vehicle, the depot is represented by the number 0, and the customer is represented by a natural number [32]. In this paper, the decimal coding based on customer arrangement is still used, but the path of each ant which is denoted represents the path of all customers served by the depot, and the depot is a separator for different vehicles. For example, the path sequence of the ant is which means that vehicle 1 leaves the depot 0, visits customers 1, 5, 3, and 7 in turn, and then returns to the depot 0; vehicle 2 leaves the depot 0, visits customers 6, 2, 8, and 4 in turn, and finally returns to the depot 0.
3.2. Probabilistic Model and Its Initialization
The performance of ACO has strong connections with the probabilistic model. A better choice of the model is critical to improve the efficiency of the algorithm. In this research, two probability matrices are employed to form the probabilistic model. Besides, a few good individuals produced by a heuristic method are used to initialize these two probability matrices. This initialization method makes the algorithm converge to the highquality solutions faster.
3.2.1. Probabilistic Model
In the traditional ant colony algorithm, the probabilistic model usually refers to ants’ pheromone matrix which is obtained by calculating the pheromone concentration . This pheromone value represents the probability that the customer appears immediately after the customer . In other words, the pheromone matrix only reflects similar blocks of customers presented in the chosen individuals. But it is not enough to accumulate this kind of information alone. Counting the chance of customer appearing before or in the current position is also important. Hence, a new probability matrix based on the order of customers called sequence information matrix is established. Then, the probabilistic model in this paper contains a pheromone matrix and a sequence information matrix.
In a subpath visited by one vehicle, let be the appearance times that customer before or in position in the chosen sequences. Next, the specific probability value in the sequence information matrix is determined bywhere is the set of all available customers not already served by the current vehicle until position . Here is an example to compute . Suppose that the path taken by three ants are , , , and , and then, is given as follows:
Since the first location is fixed as the depot, the corresponding probability values of for the second position of each customer in are , , , and .
Then, the pheromone matrix is normalized by formula (10). Finally, let denote the probability of customer at position which is calculated by formula (11). Based on this value, the ant selects the next customer to visit:
3.2.2. Initialization of the Probabilistic Model Based on a Heuristic Rule
As mentioned in the introduction, normally, the probabilistic model (that is the pheromone matrix in ACO) is initialized to the same value at the beginning of the algorithm. In order to achieve good solutions more quickly, this paper adopts the initial solutions to initialize the proposed probabilistic model but not the same value.
Sort the distance between all customers and the depot from small to large, and take the top customers into set . Let denote the first customer to be visited by the th ant, the complete path passed by the th ant (also the th initial solution), and the length of the route taken by the th ant. The procedure of the probabilistic model initialization is given as follows: Step 1: set and and . Step 2: the ant starts from the depot and visits the first customer . Then, the ant selects the remaining customers for access based on the nearest neighbor principle. This process forms the th initial solution . Step 3: calculate the route length corresponding to . Step 4: set . If , then set ; otherwise, set. Step 5: set . Step 6: if , then go to step 2. Step 7: output and .
3.3. Ant Path Construction
In this study, each ant produces a complete solution for the MCVRP. The ant begins from the depot and continuously chooses customers to its tour as long as the constraint of capacity is satisfied. If there are no feasible customers to select, the ant comes back to the depot and restarts another trip.
Let denote the amounts of ants and the amounts of the first customers selected. Set . To guarantee the diversity of individuals, the first customer is picked randomly from customers. For example, the first ant randomly selects one of the customers as the first one to visit. Then, the first customer to be visited by the second ant will be randomly selected from the remaining customers. If , when , reset to zero.
Afterwards, each ant proceeds to add customers to its route based on the probability function (12). The probability of choosing customer as the next one in the th iteration is given by formula (13):
In formula (12), contains information about how often the customer appears in the current position , is the inverse of distance between customer and , and and denote the importance of and , respectively:
Set be a constant value and a random variable which is uniformly distributed in [0, 1]. If , select the next customer in accordance with the formulation of ; otherwise, select the customer who has the maximum value of [22].
3.4. Update Mechanism of Probability Matrices
The probabilistic model consists of two probability matrices and . They are updated by the best solution found at present. Let be the pheromone persistence, then equation (14) is used to update the pheromone concentration in the path. Meanwhile, equation (15) is adopted to update the probability matrix based on the order of customers:
In equation (14), is the total length corresponding to the optimal ant path in the th iteration:
In equation (15), is the probability matrix obtained by the th optimal individual . We refer to Section 3.2.1 for more details.
3.5. SpeedUp VNDBased Local Search
Generally, for the combinatorial optimization problems, e.g., MCVRP, there are many other highquality solutions near the highquality solutions found so far. The algorithm in this section performs the problemdependent local search on the global optimal solution , for the purpose of realizing the detailed search of the area near . In this paper, the operations of insert, swap, and intersection elimination based on VND, the speedup search strategy, and the first move strategy are adopted.
3.5.1. Basic and Pipe VND
The VND is a variant of the variable neighborhood search. The basic VND [33] explores the list of neighborhoods in a sequential way one by one. Once an improvement is detected in any neighborhood structure, the search goes back to the ﬁrst neighborhood from the list. Until all the neighborhoods encounter the local optimum, the search ends. However, the pipe VND [33] executes continuous search in the same neighborhood if an improvement is obtained. In this paper, both these VNDs are applied.
3.5.2. Operation of Insert and Swap
Assume an ant’s tour is represented as . The insert operator selects one customer like and inserts it before or after another customer like . Then, the tour becomes or . The swap operator selects two customers like and on this tour and swaps them. Then, .
3.5.3. Operation of Intersection Elimination
Since the objective function of the MCVRP is constructed to minimize the total length of the path, the optimal path of each vehicle should not be crossed.
Theorem 1. In the route of every vehicle, if two lines intersect, the total path length of the vehicle is larger than that of the total length without intersection.
Proof. As shown in Figure 1, there is one vehicle route denoted by . For any two given intersecting line segments, such as lines and which intersect at point , the end points of the two segments that can form the opposite side are connected like and . Next, the new path without intersection is obtained and denoted by . According to the theorem that the sum of any two sides in a triangle is greater than the third side, in , ; in , . Let be the total length of vehicle path with intersecting lines, then . Let be the total length after eliminating intersecting lines, then . As a result, .
From and defined above, it can be seen that if two paths intersect, there are four customers involved. According to the order of customers visited by the vehicle, only exchange positions of the second and third customers and then cross paths can be eliminated. In this paper, the geometric method is used to determine whether the crossing occurs.
Suppose the coordinates of customers 1, 3, 2, and 4 in Figure 1 are , , , and , and then, the linear equations of line segments and are shown in the following equations:In formulas (16) and (17), is the slope of the line where is located, and is the slope of the line where is located.
If and intersect, then formula (18) has a solution. In other words, inequality (19) should be satisfied. Meanwhile, it is necessary to satisfy the inequality (22) in order to make two line segments intersect within the line segment (this is the case of path intersection shown in Figure 1), rather than the vertex of the line segment or its extension line. However, equations (20) and (21) are used to calculate and , respectively:
3.5.4. SpeedUp Search Strategy
Let , , and represent the operations of inserting, swapping, and uncrossing the th and th elements in solution , respectively. Let denote the set of whose neighborhood is infeasible (the infeasibility here means that the vehicle capacity constraint is not satisfied) after executing , and is the set of whose neighborhood is infeasible after executing . Because the intersection elimination operation is only carried out within the customers served by each vehicle, there is no neighborhood that violates the capacity constraint. Let denote the feasible neighborhood set based on insert operation for customers in solution , the feasible neighborhood set based on swap operation for customers in solution , and the feasible neighborhood set based on intersection elimination operation for customers serviced by the th vehicle in solution . Set to be the number of customers visited by the th vehicle, then is the position of the first customer visited by the th vehicle in solution . Hence, there are as follows:
It can be seen from (23)–(25) that is a feasible neighborhood set that any customer inserts before or after any possible customer within the same vehicle or between different vehicles. is a feasible neighborhood set that any customer swaps with any possible customer within the same vehicle or between different vehicles. is a feasible neighborhood set that customers involved in the crossline can be interchanged within the th vehicle.
For the MCVRP, if customer is too far away from customer , the total route distance after the insert and swap operation will be greater than before, which makes the two operations meaningless. Hence, in order to avoid the search of invalid areas, the distances between all possible customers and customer are arranged from small to large, and just, the first points are taken out to insert and swap. Thus, the more compact neighborhood sets and can be obtained:
For the sake of speeding up the searching process and making the neighborhood search focus on highquality areas, the nature of the problem can be utilized to further decrease invalid operations and reduce the feasible neighborhood set. Let be the value of the objective function which is the total mileage of the route corresponding to solution and be the distance between element which is located in position of and element which is located in position of .
For , if customer is inserted before customer , ; if customer is inserted after customer , . For any , if , then . in this case can be directly excluded without calculating . If , then can be obtained directly from instead of starting from the first customer of to calculate. Therefore, a more compact and effective neighborhood set can be obtained on the basis of :
For, . And for any , if , then . in this case can be directly excluded without calculating . If , then can be obtained directly from instead of starting from the first customer of to calculate. Therefore, a more compact and effective neighborhood set can be obtained on the basis of :
For , we can use the method described in Section 3.5.3 to determine whether two paths intersect. If so, according to Theorem 1, we know . Then, can be obtained directly from instead of starting from the first customer of to calculate.
Let , , and be the optimal neighborhoods in searching and outputting , , and , respectively. These three optimal neighborhoods utilize the above speedup search method in the local search of IHACO. Meanwhile, the pipe VND described in Section 3.5.1 is also applied. The corresponding procedures are given as follows: Step 1: set and . Step 2: . Step 3: if , then . Step 4: if , then and go to step 2. Step 5: if , then . If , go to step 1. Step 6: set , and . Step 7: . Step 8: if , then . Step 9: if , then and go to step 7. Step 10: if , then . If , go to step 6. Step 11: set and . Step 12: . Step 13: if , then . Step 14: if , then and go to step 12. Step 15: if , then . If , go to step 11. Step 16: output .
3.5.5. First Move Strategy
In order to reach more different regions, a first move strategy is adopted in , , and . Here is an example to illustrate the strategy. While the first neighbor in which can improve is gained, terminates and outputs . Furthermore, the sequence in is replaced by for the purpose of implementing the search instantly from the promising areas discovered by . Obviously, the algorithm will search in the vicinity of varied suboptimal solutions, but not limited to the vicinity of the optimal solution, thus expanding the scope of the search area. The changed procedures with the above strategy, named , , and , are only different from , , and in step 2, step 7, and step 12. The pseudocodes of the changed steps are as follows: Step 2: . Step 7: . Step 12: .
3.5.6. Procedure of Local Search
Let , , and denote , , and with the pipe VND method. Moreover, the basic VND described in Section 3.5.1 is also adopted. Then, the procedure of the IHACO’s local search is given in the following: Step 1. Perturbation phase. Set ; For to Randomly select and ; ; ; End. Step 2: exploitation phase. Set ; Repeat ; ; If then and go to ; else ; If then and go to ; else ; end; end; until . Step 3: if , then . Step 4: output .
In the above procedure, step 1 is a perturbation phase, which is helpful to prevent the cycling search and overcome the local optima. The special devised step 2 performs exploitation from the region obtained by step 1.
3.6. IHACO
Denote as a generation and as the maximum generation. Based on the above subsections, the proposed IHACO is given in the following: Step 1: parameter initialization. Set the crucial parameter , , , , , , , and . Step 2: probabilistic model initialization. Set //Section 3.2.2 Step 3: ant path construction. Generate the population at generation by sampling and then evaluate its individuals to obtain the global best individual //Section 3.3 Step 4: local exploitation. Execute three types of local search to and update //Section 3.5 Step 5: update mechanism of probability matrices. Calculate the objective function value corresponding to the best individual . Update and based on and //Section 3.4 Step 6: set . If , go to step 3. Step 7: output the best solution and its objective value found so far.
As it can be seen from the above steps, not only does IHACO takes advantage of the new probabilistic modelbased global search to execute exploration to obtain promising regions within solution space but also it makes use of the speedup VNDbased local search with the first move strategy to perform exploitation from these promising regions. On account of both global and local search are well balanced, the IHACO is expected to acquire good results.
4. Computational Experiments and Statistical Analysis
In this section, three sets of numerical simulations are conducted with the aim of checking the effectiveness and efficiency of IHACO. The experiments implemented have been coded in Matlab 2019a and tested on an Intel 1.8 GHz PC with 12 GB RAM.
4.1. Experimental Setup
The generation method of instances is derived from [13]. The test set in this paper contains 14 benchmark instances with the capacity constraint, including vrpcn1vrpcn5 and vrpcn11vrpcn12 of sets a and b. These scales of problem range from 50 to 200 customers.
The main parameters of IHACO are the same as those in [22]: the number of ants , two exponential factors and , the pheromone concentration , and the constant value . There are other parameters needed to set: the upper bound of range in local search and the perturbation variable . For the sake of having a fair comparison, all the compared algorithms terminate when either of the following two criteria is reached: the first one is the same runtime limit, and the second one is less than or equal to the optimal value found in the literature. It should be noted that when setting the termination time of the proposed algorithm, two steps are adopted. We first find the running time and the running environment (that is, the CPU frequency) corresponding to the best solution of the four comprised algorithms in each instance. Suppose that the CPU frequency of this paper is , then the termination time of the proposed algorithm . In addition, to make the results clearer, in Tables 2 and 3, the best, the secondbest, and the thirdbest values in each row are marked out with the bold, the bold and underlined, and the italic and underlined fonts.
4.2. Comparison of IHACO and Other Algorithms
To test the performance of the algorithm in this paper, various simulations are carried out to compare the proposed IHACO with a classical ACO (ACS) [21], an effective hybridized ACO (HAC) [22], and two new stateoftheart approaches named HAVNS and HABC [24] for the MCVRP. In solving MCVRP, ACS is the earliest used ant colony optimization algorithm and is proved with high performance for twocompartment test problems. HAC that is presented through the inspiration of ACS combines the ant colony optimization algorithm with local search schemes. After testing new generated benchmark problem instances, the results manifest that HAC takes less computational time to get higher quality solutions for all given problems except one problem (vrpnc12a) compared to ACS. HAVNS and HABC are also heuristicbased. According to the experiment results in [24], for the 14 examples to be tested in this paper, HAVNS obtains 12 best solutions and updates 6 of them on the basis of HAC. Meanwhile, HABC gains 12 best solutions and updates 8 of them based on HAC.
The experimental results of the five algorithms are shown in Table 2. The first column in Table 2 represents the type of the instance, and the second column is the quantity of customers contained in every instance. The meaning of Obj. and T(s) in this table is the objective function value and runtime, respectively.

As reported from Table 2, the proposed IHACO finds the best solution for 9 instances of the 14 problem instances and updates 8 of them based on HAVNS and HABC. However, on the average value of all examples, IHACO performs best. Moreover, while the quantity of customers is more than 100, IHACO can hit 5 best solutions in 6 instances. Hence, as the problem size increases, IHACO is more effective. In addition, HAVNS and HABC were implemented on Intel Core™ i55200U @ 2.20 GHz 6,00 GO RAM, which was different from the IHACO operating environment. Considering the speed of CPU, IHACO is much faster than HAVNS and HABC through the conversion of on average. Furthermore, as it can be seen from Table 2, the time consumption corresponding to the 9 good solutions obtained in this study is less than that of the best solution in three comparative literatures. Considering the difference of computing environment, if the conversion is carried out, the time consumption of these 9 solutions will be shorter. So, IHACO is also more efficient.
4.3. Comparison of IHACO with Its Six Variants
For the purpose of evaluating the crucial operations of IHACO, IHACO is also compared with its six variants, whose abbreviations are described in the following:(1)IHACO_V1 : in IHACO’s global search, is not adopted, and only is used to realize the generation of new individuals.(2)IHACO_V2 : in IHACO’s global search, when the probabilistic model initialization phase is carried out, is initialized to the same value , and is initialized just by the individual (the individual corresponding to the objective function value ).(3)IHACO_V3 : in IHACO’s local search, and are not employed. In other words, all possible values for customer ( is replaced with ) will be assessed.(4)IHACO_V4 : in IHACO’s local search, the first move strategy is not employed. This means that, in the operations of insert and swap, after all possible positions have been tried, and the customer chooses the location that has the highest improvement to insert or swap.(5)IHACO_V5 : in IHACO’s local search, the pipe VND is not employed.(6)IHACO_V6 : in IHACO’s local search, the basic VND is not employed.
The statistical results are listed in Table 3. Besides, the graphical results for the compared algorithms are exhibited in Figure 2.

(a)
(b)
Compared with IHACO’s variants, it is clear from Table 3 and Figure 2 that IHACO performs best for all instances. The results indicate the importance of adopting all the crucial operations of IHACO simultaneously. It can also be found from Figure 2 that, as the problem scale increases, the difference of the objective values obtained by these comparison algorithms gradually increases, and then, the superiority of IHACO becomes more obvious. This once again verified a conclusion of the previous experiment; that is, IHACO performs better when the number of customers is large.
The above results can be further analysed. IHACO_V1 does not adopt the information sequence matrix, which may lead to the low accuracy of the information left by ants on the path and reduces the similarity between the sampled individuals and the optimal individuals. IHACO_V2 initializes the two matrixes that make up the probabilistic model to a unified value. As a result, the initial individuals have no information of excellent individuals to utilize, and then, the quality of the initial population is lowered. IHACO_V3 expands the range of invalid regions in operations of insert and swap. As a result, it could take a longer time to find highquality solutions than IHACO. IHACO_V4 must find the optimal location to execute the insert and swap operations. This causes that the search areas are relatively concentrated, which may miss some promising regions. Both IHACO_V5 and IHACO_V6 reduce the depth of exploitation, so it is not easy for them to achieve highquality solutions. Thus, the main operations removed from these six variants are all adopted by IHACO.
4.4. Superiority Verification of Using MultiCompartment Vehicles
Since some different products have different characteristic properties, they cannot be mixed. That is to say, the products cannot be transported simultaneously by vehicles with single compartment. However, they can be transported by vehicles with single compartment for several times or by vehicles with multiple compartments at one time.
In this section, for the two kinds of products involved in the problem instances adopted in this paper, singlecompartment vehicles and twocompartment vehicles are both utilized to gather the needs of customers within the service area. When using singlecompartment vehicles to collect goods from customers, the problem is decomposed into two subproblems. One is to seek the solution of VRP for collecting product 1, and another is for collecting product 2. The sum of the total distance for the two subproblems represents the final solution listed in column 3 of Table 4. The results obtained by using twocompartment vehicles are the same as reported in column 5 of Table 2. The increase percentage of total length listed in column 5 of Table 4 can be obtained from the following equation:where and are the total distance travelled by singlecompartment vehicles and twocompartment vehicles, respectively. Besides, in order to make the comparison more intuitive, plot with 95% Tukey’s honest significant difference (HSD) confidence intervals for employing these two types of compartments under two categories of data sets is shown in Figure 3.

It can be obviously seen from Table 4 and Figure 3 that using singlecompartment vehicles brings the notable increase in the total length of each route. The reason is that although the total route cost is reduced due to the increase of vehicle capacity while transporting the same product, each customer must be served twice in order to ensure that different products are all delivered. More specifically, it can be noticed from Table 4 that the total distance raises 47.72% on average when the vehicles with the single compartment are employed to visit the customers. Moreover, it is worth mentioning that the total distance could be reduced by up to 70.22%, while twocompartment vehicles are selected. The statistical results clearly display the profit of using multicompartment vehicles.
5. Conclusion
The goal of this paper is to propose an alternative solution for the NPhard problem, which is dedicated to minimize the total travelled mileage for the multicompartment vehicle routing problem (MCVRP). A novel improved hybrid ant colony algorithm (IHACO) is proposed, aiming at promoting the global search to reach different regions of highquality solutions. A new probabilistic model which utilizes two probability matrices to learn different useful information of highquality solutions is developed. Then, a heuristic rule based on the nearest neighbor is used to initialize this probabilistic model so that the algorithm’s speed of reaching the promising regions is enhanced. Meanwhile, a new local search is constructed to execute the effective exploitation. The efficiency of this method is reflected in the fact that it does not need to judge the capacity constraint and evaluate the objective function. It only needs to judge whether the two routes cross or not, and then, it can directly carry out the exchange operation between two customers. Moreover, two types of VND combined with the speedup strategy and the first move strategy are employed to further perform the deep, narrow but rich search in the good solution space that has been found. Finally, numerical experiment results based on a benchmark dataset display the effectiveness and efficiency of IHACO. Future work is to develop an effective ACObased algorithm to deal with the multiobjective MCVRP and the fuzzy travel timebased MCVRP.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was supported by the National Natural Science Foundation of China under Grant nos. 61963022 and 51665025 and the Applied Basic Research Key Project of Yunnan Province.
References
 F. Cornillier, F. F. Boctor, G. Laporte, and J. Renaud, “An exact algorithm for the petrol station replenishment problem,” Journal of the Operational Research Society, vol. 59, no. 5, pp. 607–615, 2008. View at: Publisher Site  Google Scholar
 F. Cornillier, F. F. Boctor, G. Laporte, and J. Renaud, “A heuristic for the multiperiod petrol station replenishment problem,” European Journal of Operational Research, vol. 191, no. 2, pp. 295–305, 2008. View at: Publisher Site  Google Scholar
 F. Cornillier, G. Laporte, F. F. Boctor, and J. Renaud, “The petrol station replenishment problem with time windows,” Computers & Operations Research, vol. 36, no. 3, pp. 919–935, 2009. View at: Publisher Site  Google Scholar
 F. F. Boctor, J. Renaud, and F. Cornillier, “Trip packing in petrol stations replenishment,” Omega, vol. 39, no. 1, pp. 86–98, 2011. View at: Publisher Site  Google Scholar
 F. Cornillier, F. Boctor, and J. Renaud, “Heuristics for the multidepot petrol station replenishment problem with time windows,” European Journal of Operational Research, vol. 220, no. 2, pp. 361–369, 2012. View at: Publisher Site  Google Scholar
 E. E. HalvorsenWeare and K. Fagerholt, “Routing and scheduling in a liquefied natural gas shipping problem with inventory and berth constraints,” Annals of Operations Research, vol. 203, no. 1, pp. 167–186, 2013. View at: Publisher Site  Google Scholar
 R. Lahyani, L. C. Coelho, M. Khemakhem, G. Laporte, and F. Semet, “A multicompartment vehicle routing problem arising in the collection of olive oil in Tunisia,” Omega, vol. 51, pp. 1–10, 2015. View at: Publisher Site  Google Scholar
 E. D. Chajakis and M. Guignard, “Scheduling deliveries in vehicles with multiple compartments,” Journal of Global Optimization, vol. 26, no. 1, pp. 43–78, 2003. View at: Publisher Site  Google Scholar
 M. Ostermeier and A. Hübner, “Vehicle selection for a multicompartment vehicle routing problem,” European Journal of Operational Research, vol. 269, no. 2, pp. 682–694, 2018. View at: Publisher Site  Google Scholar
 S. Martins, M. Ostermeier, P. Amorim, A. Hübner, and B. AlmadaLobo, “Productoriented time window assignment for a multicompartment vehicle routing problem,” European Journal of Operational Research, vol. 276, no. 3, pp. 893–909, 2019. View at: Publisher Site  Google Scholar
 A. Huebner and M. Ostermeier, “A multicompartment vehicle routing problem with loading and unloading costs,” Transportation Science, vol. 53, no. 1, pp. 282–300, 2017. View at: Google Scholar
 L. Muyldermans and G. Pang, “On the benefits of cocollection: experiments with a multicompartment vehicle routing algorithm,” European Journal of Operational Research, vol. 206, no. 1, pp. 93–103, 2010. View at: Publisher Site  Google Scholar
 M. Elbek and S. Wøhlk, “A variable neighborhood search for the multiperiod collection of recyclable materials,” European Journal of Operational Research, vol. 249, no. 2, pp. 540–550, 2016. View at: Publisher Site  Google Scholar
 E. MofidNakhaee and F. Barzinpour, “A multicompartment capacitated arc routing problem with intermediate facilities for solid waste collection using hybrid adaptive large neighborhood search and whale algorithm,” Waste Management & Research, vol. 37, no. 1, pp. 38–47, 2019. View at: Publisher Site  Google Scholar
 P. Avella, M. Boccia, and A. Sforza, “Solving a fuel delivery problem by heuristic and exact approaches,” European Journal of Operational Research, vol. 152, no. 1, pp. 170–179, 2004. View at: Publisher Site  Google Scholar
 L. M. Hvattum, K. Fagerholt, and V. A. Armentano, “Tank allocation problems in maritime bulk shipping,” Computers & Operations Research, vol. 36, no. 11, pp. 3051–3060, 2009. View at: Publisher Site  Google Scholar
 M. F. Uzar and B. Catay, “Distribution planning of bulk lubricants at BP Turkey,” Omega, vol. 40, no. 6, pp. 870–881, 2012. View at: Google Scholar
 X. Li, S. C. H. Leung, and P. Tian, “A multistart adaptive memorybased tabu search algorithm for the heterogeneous fixed fleet open vehicle routing problem,” Expert Systems with Applications, vol. 39, no. 1, pp. 365–374, 2012. View at: Publisher Site  Google Scholar
 T. Henke, M. G. Speranza, and G. Wäscher, “The multicompartment vehicle routing problem with flexible compartment sizes,” European Journal of Operational Research, vol. 246, no. 3, pp. 730–743, 2015. View at: Publisher Site  Google Scholar
 T. Henke, M. G. Speranza, and G. Wäscher, “A branchandcut algorithm for the multicompartment vehicle routing problem with flexible compartment sizes,” Annals of Operations Research, vol. 275, no. 2, pp. 321–338, 2019. View at: Publisher Site  Google Scholar
 M. Reed, A. Yiannakou, and R. Evering, “An ant colony algorithm for the multicompartment vehicle routing problem,” Applied Soft Computing, vol. 15, pp. 169–176, 2014. View at: Publisher Site  Google Scholar
 M. M. S. Abdulkader, Y. Gajpal, and T. Y. Elmekkawy, “Hybridized ant colony algorithm for the multi compartment vehicle routing problem,” Applied Soft Computing, vol. 37, pp. 196–203, 2015. View at: Publisher Site  Google Scholar
 P. V. Silvestrin and M. Ritt, “An iterated tabu search for the multicompartment vehicle routing problem,” Computers & Operations Research, vol. 81, pp. 192–202, 2017. View at: Publisher Site  Google Scholar
 I. Kaabachi, H. Yahyaoui, S. Krichen, and A. Dekdouk, “Measuring and evaluating hybrid metaheuristics for solving the multicompartment vehicle routing problem,” Measurement, vol. 141, pp. 407–419, 2019. View at: Publisher Site  Google Scholar
 Z. Wan, G. Wang, and B. Sun, “A hybrid intelligent algorithm by combining particle swarm optimization with chaos searching technique for solving nonlinear bilevel programming problems,” Swarm and Evolutionary Computation, vol. 8, pp. 26–32, 2013. View at: Publisher Site  Google Scholar
 Q. Chen, H. Shi, and M. Sun, “Echo state networkbased backstepping adaptive iterative learning control for strictfeedback systems: an errortracking approach,” IEEE Transactions on Cybernetics, vol. 99, pp. 1–14, 2019. View at: Publisher Site  Google Scholar
 Z. C. Li, B. Qian, R. Hu, L. L. Chang, and J. B. Yang, “An elitist nondominated sorting hybrid algorithm for multiobjective flexible jobshop scheduling problem with sequencedependent setups,” KnowledgeBased Systems, vol. 173, no. 1, pp. 83–112, 2019. View at: Publisher Site  Google Scholar
 M. Mavrovouniotis and S. Yang, “Ant algorithms with immigrants schemes for the dynamic vehicle routing problem,” Information Sciences, vol. 294, pp. 456–477, 2015. View at: Publisher Site  Google Scholar
 M. Schyns, “An ant colony system for responsive dynamic vehicle routing,” European Journal of Operational Research, vol. 245, no. 3, pp. 704–718, 2015. View at: Publisher Site  Google Scholar
 X. Wang, T.M. Choi, H. Liu, and X. Yue, “Novel ant colony optimization methods for simplifying solution construction in vehicle routing problems,” IEEE Transactions on Intelligent Transportation Systems, vol. 17, no. 11, pp. 3132–3141, 2016. View at: Publisher Site  Google Scholar
 S.H. Huang, Y.H. Huang, C. A. Blazquez, and G. ParedesBelmar, “Application of the ant colony optimization in the resolution of the bridge inspection routing problem,” Applied Soft Computing, vol. 65, pp. 443–461, 2018. View at: Publisher Site  Google Scholar
 J. E. Belf and P. R. Mcmullen, “Ant colony optimization techniques for the vehicle routing problem,” Advanced Engineering Informatics, vol. 18, no. 1, pp. 41–48, 2004. View at: Google Scholar
 A. Duarte, N. Mladenovi, J. SánchezOro, and R. Todosijevi, “Variable neighborhood descent,” in Handbook of Heuristics, R. Martí, P. Panos, and M. Resende, Eds., Springer, Cham, Switzerland, 2016. View at: Google Scholar
Copyright
Copyright © 2020 Ning Guo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.