A Modified Column Generation to Solve the Heterogeneous Fixed Fleet Open Vehicle Routing Problem

In the heterogeneous fixed fleet open vehicle routing problem (HFFOVRP), several different types of vehicles can be used to service the customers. The types of vehicles are different in terms of capacity, fixed cost, and variable cost. In this problem, the vehicles are not required to return to the depot after completing a service and the number of vehicles of each type is fixed and limited. Since this problem belongs to NP-hard problems, in this paper a compound heuristic algorithm called SISEC which includes sweep algorithm, insert, swap, and 2-opt moves, modified elite ant system (EAS), and column generation (CG) is applied to solve the HFFOVRP. We report computational results on 22 problems and solve each problem by using our SISEC. The results which were compared to the results of exact algorithms and the classic CG confirm that the proposed algorithm produces high quality solutions within an acceptable computation time.


Introduction
The open vehicle routing problem (OVRP) is a variant of vehicle routing problem (VRP) which has many applications in industry and service providing businesses [1].The description of the OVRP appeared in the literature over 30 years ago but it has attracted the attention of scientists and researchers in recent years.In the OVRP, the objective is to minimize the total traveled distance (or time) by the number of vehicles so that each customer is visited only once by a single vehicle.In general, the number of vehicles used is not imposed but it is a decision variable.Furthermore, the demand of each customer must be completely fulfilled by a single vehicle and the total demand served by each vehicle must not exceed vehicle capacity; that is, split deliveries are not allowed.Moreover, in one variant of the problem, the travel time of each vehicle should not exceed an upper limit.It should be noted that the OVRP has a unique character in that it has an open path form.This means that the vehicles are not required to return to the depot after completing a service and each vehicle ends its trip at a customer [2].This problem is defined as follows.
Let  = (, ) be a complete undirected connected graph with  = {0, 1, . . ., } as the set of vertexes and the set of arcs  = {(, ): 0 ≤ ,  ≤ } (if the graph is not complete, we can instead replace each arc with an arc which has an infinite size).Vertex 0 denotes the depot, and each vertex  ∈ {1, . . ., } is a customer with a nonnegative demand 0 <   < .The cost of travel from vertex  to vertex  is denoted by   (if  = ,   = 0 for 0 ≤  ≤ ), and it is assumed that costs are symmetric (i.e.,   =   ).The available fleet consists of  vehicles located at the depot.Furthermore, each customer must be serviced by a single vehicle and no vehicle may serve a set of customers whose total demand exceeds its capacity.Each vehicle route must start at the depot and end at the last customer it serves.The objective is to define a set of vehicle routes which minimizes the total costs [3].
In the standard version of the VRP and OVRP, a fleet of homogeneous vehicles is used and exists at a single depot [4].A sequence of customers must be visited by each vehicle so that all customers are serviced exactly once and the total traveled distance by all the vehicles is minimized.In contrast to other customers, some customers may have such large demands that they can be serviced only by large capacity vehicles.Therefore, several different vehicle types are applied in the heterogeneous version of VRP [5].There are   vehicles (unlimited or limited numbers) of type  which can be used to serve the customers.When the number of vehicles of type  is unlimited, the best composition of the fleet is determined.This problem involves two kinds of decisions including (i) the choice of a mix of vehicles among the available vehicle types, (ii) the routing for each selected vehicle.
In the fixed fleet version of the heterogeneous VRP (HFFVRP), the values of   are fixed and the number of vehicles for type  is limited.If the fleet composition is known in advance, we decide how to make the best use of heterogeneous vehicles.The HFFVRP can be defined on an undirected complete graph with one depot and  customers indexed from 1 to .A fleet of identical heterogeneous vehicles with limited capacity   for each  is based at the depot.Furthermore, the vehicles do not necessarily have the same capacity, vehicle fixed cost, and unit variable cost in the HFFVRP.In other words, a capacity   , a fixed cost   , and a variable cost   are associated with each type of vehicle  so that   is cost per unit of distance corresponding to each vehicle type .Hence,    =   ×   represents the cost of the travel from customers  to  with a vehicle of type .The HFFVRP is a further generalization with a limited availability   for each vehicle type .This problem deals with finding the minimum total transportation cost including the fixed and variable cost for a fleet of vehicles so that the following constraints are taken into account: (i) Each vehicle starts at the depot and returns to the depot.
(ii) The capacity of a vehicle cannot exceed the capacity of the corresponding vehicle type.
(iii) The used number of vehicles of type  cannot exceed   .
(iv) The demand of each customer is satisfied by exactly one vehicle in only one visit.

The aim of this paper is to study the open version of HFFVRP.
In residential and commercial waste collection, there are several applications of the HFFVRP [6].In this problem, vehicles commonly have different capacities and a specific vehicle is used when the actual service at a location must be completed such as a side loader with mechanical arms which picks up the trash cans.As a result, in logistic operations, decisions relating to routing heterogeneous fleets of vehicles are frequently taken into consideration [7].The difference between the open HFFVRP and the HFFVRP is the fact that a solution of the former consists of a set of Hamiltonian paths, rather than Hamiltonian cycles [8,9].On the other hand, the OVRP is not very realistic since most companies have heterogeneous fleets of vehicles.If using heterogeneous fixed fleet as constraints are assumed for the OVRP, we call the problem a heterogeneous fixed fleet open vehicle routing problem (HFFOVRP) which is a relaxation of the standard OVRP.In the HFFOVRP, the fleet is composed of several vehicle types.Unlike the VRP, the number of vehicles of each type is limited [10].Moreover, the vehicles have different capacities, fixed costs, and variable costs per unit distance.The idea of this problem is not only to consider the routing of the vehicles but also to take into account the composition of the vehicle fleet.This problem is more practical in distribution management and transportation [11].In more detail, the HFFOVRP appears in situations where the fleet is purchased: it combines tactical decisions (selecting the number of vehicles to be acquired) and operational ones (computing the trips and the vehicles assigned to them).
For the first time, the HFFOVRP was addressed in [12].In this paper, each arc(, ) is associated with a traveling distance   , a service time   , and the maximum allowed traveling distance   .They proposed a multistart adaptive memory programming metaheuristic called MAMP with modified tabu search algorithm for solving this problem.This algorithm is based on the framework of the adaptive memory programming procedure constructed by Taillard et al. [13].The major mechanism considered is to produce a high quality HFFOVRP by using the search memory and experience.The algorithmic efficiency and effectiveness are experimentally evaluated on a set of generated instances including 4, 15, and 20 solutions.The maximum number of iterations in all MAMPs is set to be 2500.Because the HFFOVRP is a new variant of the classical OVRP in that time, and there has been no related work done in the literature, the proposed algorithm is only compared with a tabu search algorithm because it is similar to the solution improvement procedure in MAMP.
Huachi et al. deal with both versions of the OVRP including homogeneous and heterogeneous fleet (HFFOVRP) in which their goals are to minimize traveled distance and the total costs [14].In contrast to the HFFOVRP, the number of vehicles must be minimized when the fleet is homogeneous.The objective of this work is to present a heuristic algorithm based on the iterated local search (ILS) which uses a variable neighborhood descent (VNS) procedure in the local search phase.In this phase, the VNS is combined with random variable neighborhood ordering (RVND).The developed algorithm was tested in benchmark instances with up to 480 customers and the results obtained showed that the algorithm is quite competitive with those found in the literature.This approach is an extension of [15] for the Heterogeneous Fleet Vehicle Routing Problem (HFVRP).
Yousefikhoshbakht et al. proposed a bone route algorithm (BRMTS) which uses the tabu search (TS) as an improved procedure to solve the HFFOVRP [9].This population-based algorithm was first proposed for solving the classic VRP by Tarantilis and Kiranoudis [16] which constructs a new solution using components of routes in the adaptive memory.In that work, they present a version of this method which directly produces a new solution from a component of the other solutions while using new diversification and intensification mechanisms.The BRMTS employs the generalized route for constructive algorithm (GEROCA) which was first presented in [17], for generating initial diversified solutions and modified TS improvement procedure [15].They have made two main modifications to the TS: (i) To improve the TS further, the amount of tabu tenure (TT) is considered variable for intensification and diversification policies.The TT is considered a minimum value for the intensification policy and a maximum value for the diversification policy.
(ii) The proposed TS comprises three types of neighborhood moves including 2-opt, 0-1, and 1-1 exchanges.These moves are distinguished in terms of exchanges performed to convert one tour into another.These moves are not equally performed in each iteration and probability of them is changed gradually during the run time.
Furthermore, some test problems were derived from the well-known Taillard's benchmark for HFFVRP and the results of BRMTS algorithm are compared with the results of several metaheuristic algorithms.The comparison showed that the proposed algorithm can provide better solutions within a comparatively shorter period of time.
The HFFOVRP plays an important role in supply chains specially in the first transportation step such as collecting agricultural products or in the final distribution phase toward customers.Besides, this problem is obviously NP-hard since it includes the OVRP as a special case.For such problems, the use of heuristics such as column generation (CG) is considered a reasonable approach in finding solutions.Although CG has been successfully applied to several combinatorial optimization problems, it has some shortcomings like its slow computing speed and local-convergence.Besides, despite the developments in modern heuristics, this algorithm has its own strengths and weaknesses.Therefore, many research studies have tried to develop hybrid algorithms, expecting to achieve the effectiveness and efficiency.In this paper, we have proposed an efficient hybrid heuristic called SISEC in order to improve both the performance of the algorithm and the quality of the solutions.The proposed algorithm took advantage of various versions of sweep algorithm with insert, swap, and 2-opt moves for producing a feasible solution.Furthermore, a tight integer programming model is presented.The linear programming relaxation is solved by the CG technique and an elite ant system (EAS) and several local search algorithms are used for solving generated Hamiltonian paths with some modifications to generate feasible columns efficiently.The computational results are provided on a new set of generated instances derived from HFFVRP's problems.The results show that the hybrid algorithm explores different parts of the solution space and the search method is not trapped at the local optimum.Besides, the proposed approach is able to find high quality solutions in a very short time compared to exact and heuristic algorithms.
This paper is structured as follows.In Section 2, in addition to clarifying the intricate history of HFFOVRP and providing bibliographic references, we review three algorithms which have been developed to solve this problem and its extensions.Components for CG which are able to tackle the HFFOVRP are described in Section 3. In Section 4, we first generate 22 new test problems which have 10-100 customers and the results of the proposed algorithm are reported for these instances.Moreover, the proposed algorithm is compared to exact and CG algorithms.Finally, some concluding remarks are given in Section 5.

Related Works
The heterogeneous open vehicle routing problem (HOVRP) can be studied in two different versions.Some researchers consider an assumption that there is an unlimited number of vehicles of each type and try to find the optimal set of vehicles to be scheduled in the problem.This is called the fleet size and mix OVRP (FSMOVRP).However, to the best of our knowledge, no work has addressed this version of VRP which is more practical in distribution management and transportation.Other researchers study the case where there is a fixed vehicle fleet and try to schedule this fleet of vehicles to the customers in an optimal way.This problem is called heterogeneous fixed fleet OVRP (HFFOVRP) as mentioned before.This HFFOVRP is more realistic than FSMOVRP and has attracted more attention in literature.In the fixed fleet version of the OVRP, the values of   are fixed.More specifically, the number of vehicles of type  is limited and the fleet composition is known in advance.We must decide how to make the best use of a fixed fleet of heterogeneous vehicles.
The overarching idea is that many linear programs like HFFOVRP are too large to consider all the variables explicitly [18].Since most of the variables will be nonbasic and assume a value of zero in the optimal solution, only a subset of variables need to be considered in theory when solving the problem.Column generation leverages this idea to generate only the variables which have the potential to improve the objective function-that is, to find variables with negative reduced cost (assuming without loss of generality that the problem is a minimization problem).So, the column generation is one of the effective algorithms in order to solve an NP-hard problems.
The HFFOVRP is defined on an undirected connected graph  = (, ), where  = {0, 1, . . ., } is a vertex set and  = {(, ): ,  ∈ } is an edge set.Vertex 0 represents the depot and each other vertex  ∈  \ {0} is a customer with a demand   .A distance   is associated with each edge (, ) ∈ .There is a fixed fleet of  different types of vehicles located at the depot and the number of vehicles of each type is   ( ∈ {1, 2, . . ., }).It should be noted that there is no guarantee that a feasible solution can be found because the number of vehicles of each type is fixed.It means that if this happens in real-life problems, other vehicles have to be hired in order to serve the customers.A capacity   , a fixed cost   , and a variable cost   are associated with each type of vehicle .In this problem, a sequence of deliveries is generated for a heterogeneous fleet of vehicles, in order to fulfill a known client requirement such that (1) each route starts at the depot and ends at an arbitrary customer, ( The HFFOVRP consists of defining a set of routes and vehicles as assigned to each route so that no more vehicles than those available are used; customers' demands are satisfied; each customer is visited exactly once; each vehicle route starts at the depot and does not finish at the depot; and the capacity of each kind of vehicle is not exceeded.The objective of the HFFOVRP is to minimize the sum of the costs of all the routes subject to the mentioned constraints.The formulation of the HFFOVRP requires the following additional notation: Now, HFFOVRP can be formulated as follows: In this formulation, the objective function (2) minimizes the total costs.Constraints (3) indicate that each customer is assigned only once to a feasible route-vehicle pair.Constraints (4) ensure that a maximum number of used vehicles of each type does not exceed the number of available vehicles of the same type.Finally, constraints (5) describe that each arc in the network has the value 1 if it is used and 0 otherwise.In this type of formulation, each column corresponds to one route-vehicle pair.As number of columns is extremely large, this problem cannot be solved directly and consequently a CG approach is used to solve this problem.

The Proposed Algorithm
Now the proposed algorithm which uses the EAS as an improved procedure is presented in this section.In our approach, we have made three main modifications to the CG: (1) To produce a high quality initial solution, the modified sweep algorithm, swap, insert, and 2-opt move are used.(2) A route-vehicle pair in open vehicle routing (, ) is feasible when route  starts at the depot and the sum of demand of customers on route  is not larger than the capacity of the vehicle type .In the proposed algorithm, each (, ) is solved by the modified elite ant system (MEAS) proposed by Yousefikhoshbakht et al. [19].(3) To improve the initial solution, the CG and several local search algorithms are applied.
In this section, first, the sweep, insert, swap, and 2opt algorithms are explained and then the MEAS and the proposed CG are described in more detail.

The Proposed Modified Sweep Algorithm.
In order to be applied, the proposed algorithm requires a high quality initial solution.Here, a modified sweep algorithm combined with insert, swap, and 2-opt moves is used to obtain initial solutions.The classic powerful sweep algorithm is relatively simple and was introduced by Gillett and Miller in 1974 [20].It perhaps is one of the most greedy construction algorithms for the traveling salesman problem (TSP) and VRP.The sweeping clockwise or counterclockwise of this algorithm starts from customer  which has not been visited, from the unvisited customer with the smallest angle to the unvisited customer with the largest angle until all customers are included.In other words, this algorithm has the simple idea of inserting the unvisited customer with the nearest angle of the last inserted customer in the route.This method can produce a set of initial feasible open routes solution for HFFOVRP in order to use as an initial solution for other algorithms.This algorithm ranks and links demand points by their polar coordinate angle.To apply the proposed modified sweep algorithm (PMSA), we assume that the depot is node 0 and each customer has an (, ) coordinate as its location.Therefore, polar coordinate angle for each customer  is calculated by formula (6).It should be noted that if () − (0) < 0 and () − (0) ≥ 0 then − < () < 0 and 0 ≤ () ≤ , respectively.Besides, if customer  and customer  have the same polar angle, we put a customer which has a shorter distance with respect to the depot earlier in the list: In the proposed algorithm, first, the polar coordinates of all  customers are calculated and ordered as (1) ≤ (2) ≤ ⋅ ⋅ ⋅ ≤ ().Then, the route assigned to the chosen vehicle  starts with the unrouted customer farthest from the depot such as .The next customer  to be inserted in the route will be the one who has not been served yet and who is the nearest polar coordinates to the customers of the route.If    is load carrying capacity of vehicle  in customer , that customer  also has to be admissible in    +   ≤   .This process is repeated until no customer is admissible in the current route.When this happens, a new vehicle is selected and the whole process is repeated until all the customers are routed.In the proposed algorithm we assume a vehicle with infinite capacity and fixed cost.If the modified sweep algorithm cannot produce a feasible solution, all the remaining customers are allocated to this artificial vehicle.
After producing a feasible or unfeasible solution, the insert and sweep algorithm is used to obtain a feasible solution with high quality.These neighborhood structures are an important key feature in the performance of most of effective algorithms because it determines the extent and the quality of the solution space explored.These moves are distinguished in terms of exchanges performed to convert one tour into another and they can be described as follows.
Besides, in the swap algorithm two nodes from different routes are exchanged [4].Consequently, if it is supposed that the initial tour consists of the set of nodes (0, . . .,  − 1, ,  + 1, . . .,  − 1, ,  + 1, . . ., 0), the improved one is constructed as (0, . . .,  − 1, ,  + 1, . . .,  − 1, ,  + 1, . . ., 0).The same procedure is conducted in the case of multiple routes.The move is allowed when it is considered favorable for the performance of the entire algorithm in terms of objective and constraints.This move is demonstrated in Figure 1(b).The most commonly encountered move is the 2-opt.In multiple routes, edges (,  + 1) and (,  + 1) belong to different routes but they form a criss-cross again.The 2-opt move is applied exactly in the same way as is the case in multiple routes.The move is allowed when it is considered favorable for the performance of the entire algorithm in terms of objectives and constraints.This move is demonstrated in Figure 1(c).

Column Generation
Approach.Two problems are in CG algorithm including the Master Problem (MP) and the Pricing Subproblem (PSP).The MP is linear programming problem whose aim is to find the minimal cost and to produce the shadow prices of the temporary optimal solution to be used in the Pricing Subproblem.Since the number of columns of MP might be large, a Restricted MP is used to initiate the computations.The goal of the Pricing Subproblem is to generate additional columns for MP.In this section, the formulations of the Master Problem and Pricing Subproblem are presented.

Master Problem.
The MP can be formulated as a set covering formulation.This formulation is a linear programming problem to which each customer is assigned at least one feasible route-vehicle pair since the arc costs satisfy the triangle inequality and each customer is visited exactly once to minimize the costs.In addition, in constraint (5), the integrality of variable    is eliminated: Creating all feasible routes is considered as an NP-hard problem.The main idea of column generation approach is to use only a small number of feasible routes in order to find the optimal solution out of a large set of possible feasible routes and additional routes are added only when needed.Therefore, at first, we considered the Restricted Master Problem (RMP) containing only the routes which have been generated by the proposed modified sweep algorithm.If the generated routes for vehicle type  are indicated by    ⊂   , then the RMP is created by replacing   with    in the MP.The RMP is solved from this solution and we are able to obtain shadow prices for each of the constraints in the RMP.This information is then utilized in the objective function of the Pricing Subproblem.Each time the Pricing Subproblem is being solved, new routes are being generated as needed and inserted in the set   = ⋃ ∈    .We solve the RMP by solver by modified elite ant system and vector  * = ( *  1 , . . .,  *   ,  * 1 , . . .,  *  ) is the values of the dual variables (shadow price) corresponding to constraints ( 8) and (9).

Pricing Subproblem. The Pricing Subproblem (PSP) is a new problem created to identify a new route-vehicle pair
with the minimum reduced cost.The PSP uses the shadow price information from RMP to generate promising a routevehicle pair that is also feasible.We formulated a Pricing Subproblem for vehicle type  (PSP()) and consequently have  problems.We consider vehicle type  and  ∈   .This open route is as follows: The reduced cost    of  ∈   is defined as where  *  0 =  *  .Therefore, the reduced cost    is equal to the sum of the fixed cost and the cost of route  on the directed subgraph   = (,   ) for vehicle type  with arc modified cost defined as follows: Now, the PSP formulation can be solved by modifying the elite ant system which is one of the most important versions of the ant colony optimization (ACO).The ACO is used to solve the traveling salesman problem (TSP) in 1991 [19].This algorithm simulates the natural ant behavior for food finding and applies it for solving the NP-hard problems for which no effective algorithm has been found yet.Studies on real ants show that, despite the fact that they do not enjoy a complete sense of vision, they can find the shortest path from the food sources to the nest by some evaporated material called pheromone [4].This material is secreted by the ants when they move from one place to another to find the shortest path.The first version of ACO family was the ant system (AS) used for solving small scale TSP instances.Although it gained good solutions for small scale problems of TSP type, the AS algorithm failed to reach an acceptable efficiency in large scale problems of TSP in comparison with the famous algorithms of that time.Therefore, much effort was made to solve the problem efficiently [4].The first modification applied on the AS algorithm was the usage of elitist strategy called EAS published in 1996.In EAS algorithm like the AS, the probability of movement from node  to node  which is not visited yet by ant "" is presented as follows: where    is the collection of nodes which have not been visited by ant "" yet,   () is the value of pheromone on the arc joining  to ,   () is the heuristic information for the ant visibility measure (e.g., defined as the reciprocal of the distance between node  and node  for the TSP), and ,  are the controlling parameters by user which determine the ratio of importance of ant's visibility measure compared to the value of pheromone released on arc(, ) [4].
This process further highlights the arcs belonging to the best route in any iteration and causes them to be updated according to the value of the best route   .Note that the above operator indicates that the less the value of   is, the more pheromone is released on the arcs.
Updating the pheromone simulates the changes in values of pheromone in any iteration and mainly it is one of the reasons that algorithms are different.Generally, two operations motivate this updating procedure in EAS algorithm [4]: (i) 1st, releasing new pheromone on the arcs, locally and globally.This operation leads to an increase in pheromone secretion on the arcs.Based on this algorithm, in addition to the local releasing of pheromone on the arcs through which the ants have passed, pheromone is released on the arcs belonging to the best route  * and they are encouraged with the constant coefficient  in the following way [4]: (ii) 2nd, the evaporation of pheromone.At the end of each iteration of the algorithm, the value of pheromone left on the arcs is decreased by the constant coefficient 0 ≤  ≤ 1.Thus, the new footprint of pheromone has an average weight between the value of the pheromone left on the arcs, and the value of new pheromone released in the arcs.Thus the formula of updating pheromone in the EAS algorithm is where Δ   () is the formula of locally updating the pheromone [4].Thus ants passing over the arc between nodes  and  release pheromone on the respective edge which is equal to the inverse of the cost of the tour (  ()) taken by ants by formula (17).In this formula,   is the collection of arcs passed over by the ant "": Therefore, the result of PSP() will be one route for vehicle type .If the best objective value of PSP() is negative, the route-vehicle pair having the minimum reduced cost    will be added to the set    in RMP.After this step, the primal model shown in Section 3 is solved by an exact algorithm and the best solution is considered .For this solution, the insert and swap algorithms are used to further improve the best solutions and then the new produced route-vehicle pairs are added to the set    in RMP.In this level, the RMP will be solved again and shadow prices are obtained for each of the constraints in RMP.For vehicle type , arc modified cost    is obtained and the process is repeated until no routes with negative reduced cost are identified.The proposed algorithm will be described as follows:

Pseudo Code of the Proposed Approach
(1) Find initial sets   of routes by sweep, insert, swap, and 2-opt algorithms for the MP.
(2) Solve the MP by using Solver CPLEX 12.3 and obtain the shadow prices of the optimal solution.
(4) For vehicle type , solve () by using the MEAS and find feasible route-vehicle pairs with negative reduced costs and add them to the sets    .( 5) Solve the primal model in Section 3 by exact algorithm and produce .Then, the insert and swap algorithms are used to further improve .
(6) Add route-vehicle pairs to the sets    .
(7) If the solution obtained from the last RMP is integer, the approach terminates; otherwise replace   in HFFOVRP formulation by    ; then solve it and the approach terminates.

Computational Results
To evaluate the validity of our proposed SISEC algorithm compared to other algorithms including exact algorithm (EA), EAS, genetic algorithm (GA), and CG for the HFFOVRP, the performance of our algorithms was tested on a set of 22 benchmark instances derived from well-known Taillard's benchmarks in this section.These algorithms can be downloaded from the website [21].The proposed algorithm was coded in AIMMS with solver GUROBI 4.5.One of the most advanced development environments for building optimization-based operations research applications and advanced planning systems is AIMMS.This commercial linear programming software could find optimal solutions for the small scale of the problems, like the HFFOVRP, and hence can be used to evaluate the accuracy of the proposed model.Based on this information, AIMMS obtained the optimal solution only for a few instances and in the other instance automatically terminated before reaching to optimal solution.Also in the large problems, this algorithm failed to obtain a feasible solution.
Besides, all the experiments were implemented on a PC with Pentium 4 at 3.5 GHZ and 8 GB RAM running Windows 7 Home Basic Operating system.First, the benchmark problems are described and then some numerical results of comparison between the proposed algorithm and some algorithms are presented on a set of HFFOVRP benchmark problems.Table 1 shows these problems are built of 10-100 nodes including the depot which all randomly located over a square with no service time.Euclidean distances are used in all the problems.Besides, they have a fixed fleet with capacity restrictions and with no route length restrictions.The proposed compound heuristic algorithm is compared to mentioned algorithms in the Table 2.The information of this table includes the number of customers (), the best solution costs of EA, EAS, GA, CG, and SISEC, and their CPU time, respectively, the best known solutions (BKSs), and gap between the best known solution (BKS) and our hybrid algorithm.The gap is a simple criterion to measure the efficiency of an algorithm and to compute the relative percentage deviation of its solution from the BKSs on specific benchmark.These values are calculated by the following formula.A zero gap indicates that best known solutions should be obtained: It should be noted that while EA, CG, and SISEC have been run until the software stops automatically, the stop condition of EAS and GA is the best know solution iterated ten times.The results of Table 2 show that the EA fails to find optimal solutions for most of the problems especially  In more detail, the results indicate that SISEC is a competitive approach compared to the EA.For instances 1, 2, 3, and 4, the gap between SISEC and BKS is zero.However, in other instances, the proposed algorithm finds better solution than the EA.Moreover, computational results of the proposed algorithm and CG show that these algorithms have a close competition for three instances including 1, 2, and 3.In addition, the proposed algorithm produces better 19 solutions than CG.In more detail, the SISEC has better solutions with more than 5% of gap in 10, 13, 20, and 21.Moreover, the proposed algorithm obtains nearly the solutions of CG for 4 instances including 5, 6, 8, and 9; that is, the gap is below 1% and for other instances the gap is between 1% and 5%.For overall, the performance of the proposed algorithm is better in reaching the suboptimal solution than the CG.
Figure 2 presents the comparison of CPU time between SISEC and CG algorithms for the 22 instances.In this figure, the horizontal axis shows the name of instances and the  In Table 3, columns 1 and 2 show the names of instances and the number of customers (), respectively.Furthermore, columns 3, 4, and 5 contain the EA, CG, and SISEC on the same CPU time which is shown in column 6.Finally, the BKSs for all instances are shown in column 7. From the comparison between exact algorithm and SISEC on five instances including 1, 2, 3, 4, and 5, it can be seen that the proposed algorithm in 4 instances has been able to find better solutions than the exact algorithm.In more detail, the exact algorithm has produced equal solutions only in instance 1.Furthermore, the CG has had a weak performance in general and has not been able to produce better solutions in all instances compared to SISEC.In other words, the CG has obtained equal solutions only on instances 1, 2, and 8. Therefore, the SISEC has been able to escape local optimum points and has caused a satisfactory improvement in the performance of the algorithm.
Figure 3 shows the comparison of our algorithm and CPLEX 12.3 (EA) and CG for five instances including 1, 2, 3, 4, and 5.In this figure, the horizontal axis shows the name of instances and the vertical axis indicates the results of three mentioned algorithms.This figure shows that proposed SISEC clearly yields better solutions than the other algorithms.

Conclusions
This paper investigates the HFFOVRP in the transportation system and presents a new combined heuristic algorithm based on CG called SISEC which has allowed us to improve the quality of several of the solutions obtained by EA, EAS, GA, and classical CG.Computational results generally have shown that the proposed algorithm yields better results compared to the exact and CG algorithms in terms of both the solution quality and the CPU time.It seems that using metaheuristic algorithms like ant colony optimization can lead to gaining better initial solutions.Furthermore, powerful metaheuristic algorithms can be used for solving this problem and other versions of HFFOVRP.Future projects will focus on working on such ideas and making them operational.
the SISEC − Value of the BKS) Value of the BKS .

Figure 2 :
Figure 2: The ratio of SISEC CPU time to CG CPU time.
difference between the open and closed versions.To transform the open version into the closed one, it suffices to set the cost to zero for traveling from any customer to the depot.However, if travel costs are symmetric, things are more subtle.Indeed, it is proved that, somewhat surprisingly, the open version turns out to be more general than the closed one in the sense that any closed HFFVRP on  customers can be transformed into an open HFFVRP on  customers, but there is no transformation in the reverse direction.
At first sight, having open routes instead of closed ones looks like a minor modification.Indeed, if travel costs are asymmetric, there is essentially no is the set of feasible routes for vehicle type ,  is the set of all of feasible routes ( ∈ ,  = ⋃ ∈   ),    is the cost of route  for vehicle type ((, ) ∈   .
= ∑ (,)∈   ), and   is the flow variables specifying the quantity of goods which vehicle  is carrying when it leaves customer  to service customer .The traveling cost of each edge (, ) ∈  by a vehicle of type  is

Table 1 :
Data for problems.

Table 2 :
Comparison results of the SISEC with some heuristic and metaheuristic algorithms.

Table 3 :
Comparison results between algorithms in the same run time.