Research Article  Open Access
Ahmed Tarajo Buba, Lai Soon Lee, "Hybrid Differential EvolutionParticle Swarm Optimization Algorithm for Multiobjective Urban Transit Network Design Problem with Homogeneous Buses", Mathematical Problems in Engineering, vol. 2019, Article ID 5963240, 16 pages, 2019. https://doi.org/10.1155/2019/5963240
Hybrid Differential EvolutionParticle Swarm Optimization Algorithm for Multiobjective Urban Transit Network Design Problem with Homogeneous Buses
Abstract
This paper considers an urban transit network design problem (UTNDP) that deals with construction of an efficient set of transit routes and associated service frequencies on an existing road network. The UTNDP is an NPhard problem, characterized by a huge search space, multiobjective nature, and multiple constraints in which the evaluation of candidate route sets can be both time consuming and challenging. This paper proposes a hybrid differential evolution with particle swarm optimization (DEPSO) algorithm to solve the UTNDP, aiming to simultaneously optimize route configuration and service frequency with specific objectives in minimizing both the passengers’ and operators’ costs. Computational experiments are conducted based on the wellknown benchmark data of Mandl’s Swiss network and a large dataset of the public transport system of Rivera City, Northern Uruguay. The computational results of the proposed hybrid algorithm improve over the benchmark obtained in most of the previous studies. From the perspective of multiobjective optimization, the proposed hybrid algorithm is able to produce a diverse set of nondominated solutions, given the passengers’ and operators’ costs are conflicting objectives.
1. Introduction
The rapid growth of many cities in the world as a result of urbanization as well as concern over the environment impact requires systematic planning to address the corresponding increase in travel demand, pollution, and energy consumption. The solution to these problems is to utilize an efficient public transport system. Adequately designed transit network with appropriate vehicle frequency is capable of raising the utilization of the public transportation system, as well as reducing the overall system costs.
Urban transit network design problem (UTNDP) consists of determining route networks and schedules of an urban public transport system [1]. Ceder and Wilson [2] decomposed the UTNDP into a sequence of five activities: network design, frequency setting, time table development, bus scheduling, and driver scheduling. Because of its combinatorial complexity, multiobjective nature, use of assignment submodel, and spatial route layout, solving the UTNDP becomes more complex and time consuming [3, 4]. Most of the research efforts attempted to treat each of the activities in a sequential manner.
Solving the UTNDP requires several criteria to be considered to efficiently account for the quality of service offered to the passengers while at the same time reducing the cost of running the service. For instance, the passengers would prefer cheap, fast, and reliable services, while the operators need to consider the cost of running the services [1]. Hence, the different viewpoints of passengers and operators can result in different criteria for evaluating efficiency. Specifically, the UTNDP involves several conflicting objectives, such as the total passenger travel time and the number of vehicles required to operate the routes. Because of the complex nature of the components of the transit travel time, comprising invehicle travel time, waiting time, transfer time, and transfer penalties, to optimize the transit networks has been a challenging task [5].
With the improvement in search algorithms and computing technology, a variety of metaheuristic approaches have been utilized to tackle the UTNDP. Apart from solution methods using sole algorithms, Farahani et al. [6] noted that solution methods can also be developed for new and existing network design problems based on hybrids of different metaheuristics. In addition, the combination of two powerful search algorithms has been considered by several authors to achieve better quality results, particularly when one algorithm complements the other (see [7–8]).
Recent studies on UTNDP focus primarily on metaheuristics as a result of highquality solutions produced and computational performance. All the researchers in the UTNDP literature have attempted to calibrate their studies either on reallife or theoretical networks. Several authors employed a genetic algorithm (GA) with diverse formulations and coding schemes with the objective to optimize either the passenger or operator cost or both. In particular, Ngamchai and Lovell [10] and Fan and Machemehl [11] utilized a GA on theoretical networks, while Arbex and da Cunha [12] experimented on Mandl’s Swiss network. Furthermore, Gundaliya et al. [13] and Tom and Mohan [14] studied mediumsize network in Chennai, India, while Amiripour et al. [15] studied a network in Iran. More recently, Owais and Osman [16] experimented a GA on a transportation network in Rivera City, Northern Uruguay.
Over the last decade, other metaheuristic approaches such as particle swarm optimization (PSO), bee colony optimization (BCO), and differential evolution (DE) are also applied to solve the UTNDP. Blum and Mathew [17] developed an intelligent agent optimization system for the UTNDP. Bagherian et al. [18] formulated a mixed integer model and applied a discrete PSO to determine the optimal bus lines. Kechagiopoulos and Beligiannis [19] developed a PSO for tackling the urban transit routing problem (UTRP). Nikolić and Teodorović [20] developed an efficient algorithm based on BCO for UTNDP. Zhao et al. [21] proposed a mimetic algorithm (MA) to solve the urban transit network. Buba and Lee [22] proposed a DE for UTRP. Buba and Lee [23] proposed a DE to solve the UTNDP by simultaneously determining the transit network configuration and the associated frequency of service. These studies experimented on the benchmark Mandl’s Swiss network and have been extensively reviewed in Buba and Lee [23]. RuanoDaza et al. [24] studied a transit network design and frequencysetting problem in the context of bus rapid transit system (BRTS) using a globalbest harmony search coupled with a simulation model for discrete events. In the bilevel multiobjective approach, the external level handled the problem of selection of the best transit route configurations, while the internal level selected the best frequencies generated by the first level. The proposed approach is applied to a reallife BRTS of the city of Pereira, Columbia. Most recently, Jha et al. [25] studied the transit network design and frequencysetting problem for public buses. The authors used a multiobjective PSO with a number of search strategies to tackle the problem. The proposed algorithm has also experimented on the benchmark Mandl’s Swiss network.
Furthermore, Canca et al. [26] formulated an optimization model that addresses the integrated network design, line planning, and fleet investment. Later, in 2017, Canca et al. [27] formulated a model to solve the integrated rail rapid transit network design and line planning problem. An adaptive large neighborhood search metaheuristic is used to tackle the network design and line planning problems simultaneously. LópezRamos et al. [28] proposed an optimizationbased approach to simultaneously solve the network design and the frequencysetting phase in the context of railway rapid transit networks. A combined lexicographic goal programming technique and a line splitting algorithm is used to solve the model. GutiérrezJarpa et al. [29] investigated the rapid transit network design with modal competition in the case of a real city through a multiobjective mathematical model for a rapid transit network at a strategic level. In addition, an approach that can address modal competition for realistic size instances is also presented by the authors.
Over the decades, the hybridization approaches for UTNDP are gaining attention from the researchers. Zhao and Zeng [7] applied a hybrid SAGA algorithm to optimize the transit route network design problem. Subsequently, Zhao and Zeng [30] extended their method to determine the public transport network layout and headways. Liu et al. [8] proposed a hybrid SAGA strategy to solve the problem of bus route design and frequency setting. Szeto and Wu [9] proposed a hybrid approach to solve a transit network design problem for a suburban residential area in Hong Kong. Szeto and Jiang [31] proposed a model for solving UTNDP with the objective of minimizing the weighted sum of the number of transfers and the total travel time of users through a hybrid enhanced artificial BCO. Later, in 2014, Szeto and Jiang [32] developed a hybrid artificial BCO to solve a bilevel UTNDP where transit route design and frequency settings are determined simultaneously.
In this paper, we consider the UTNDP aiming to determine a set of transit routes with associated service frequencies simultaneously for transit networks with homogeneous buses, with the objective of minimizing the conflicting interests of passengers and operators. There are a few research efforts that have used the hybrid differential evolution with particle swarm optimization (DEPSO) to solve a range of continuous optimization problems (see [33–35]). However, to the best of our knowledge, no such research has been reported in the UTNDP literature. Hence, the main motivation is to propose a hybrid DEPSO to solve effectively and efficiently the UTNDP. The proposed hybrid algorithm is an extension from our previous study [23], which is designed to adapt the strength of both the approaches in solving the discrete problem aiming to achieve better quality results.
The remainder of the paper is organized as follows: the problem formulation is presented in Section 2, while the details of the proposed hybrid DEPSO framework for UTNDP is described in Section 3, followed by the experimental design in Section 4. Comparisons of results are presented in Section 5. Lastly, conclusions and future research are addressed in Section 6.
2. Problem Formulation
The problem considered in this article comprised two sequences of activities: (i) bus network design and (ii) frequency setting. The network design aspect consists of constructing a feasible set of transit routes on a given urban road network with predetermined stop points purposely to achieve an efficient transit network that optimizes passengers’ total travel time as well as the operator cost. The frequency setting aspect involves assigning service frequencies to transit routes so that several parameters including waiting times, flow capacities for transit routes, and the fleet size required for overall network operation are established.
The UTNDP is modeled as a multiobjective optimization problem so that an efficient set of transit routes and the associated frequencies are established that simultaneously produce a minimum total passenger and operator costs while all the requirements and constraints are being satisfied. In a reallife scenario, a number of constraints need to be included in the formulation such as bus line feasibility, bus capacity, and fleet size. The same design approach in Nikolić and Teodorović [20]; and Zhao et al. [21] is considered in this study. The problem can be defined as follows.
Let be a weighted graph representing the candidate transportation network with vertex set V (demand points) and edge set (street segment). is the set of feasible edges () that links the vertices and satisfying origindestination (OD) pairs. Each edge is linked to a pair of vertices , where is the source and is the sink of the edge. A symmetric matric of passenger travel demand and overall travel time associated with each OD pair are known. Consequently, should be connected to enable passengers to traverse the path between node and node through an undirected path in . The optimization model is as follows:where = set of transit routes in a solution to UTNDP; = invehicle time for travel from nodes to ; = waiting time for travel from nodes to ; = transfer time for travel from nodes to ; = demand for travel from nodes to ; = unserved demand for travel from nodes to ; = penalty factor of the unserved passengers; = transit route/bus line ; = fleet size on route ; = minimum length of a route for travel; = length of route for travel; = maximum length of a route for travel; = frequency of buses plying on route ; = minimum allowable frequency of buses operating on a route; = maximum flow occurring on any link of route ; = flows on the critical link of route ; = round trip time on route ; = fleet size available for operation on the route network; = number of transfers for travel from nodes to ; = maximum allowable transfer; = number of routes in the given route set; = maximum number of routes allowed in a given solution; = number of nodes in the route set (graph size).
The objective functions (equation (1)) represent generalized passenger cost, , in terms of total travel time of all passengers plus the unmet demand, whereas equation (2) is the operator cost, , (total fleet size), respectively. Constraint (3) provides the lower and upper limits of the route length in a given solution. Constraint (4) ensures that service frequency should not be lower than the minimum value. Constraint (5) limits the maximal flows on the critical link of transit routes. Constraint (6) defines the relationship between fleet sizes and frequencies of transit routes. Constraint (7) ensures that the sum of fleet size does not exceed the maximum allowable value from the operators’ point of view. Constraint (8) restricts the number of transfers to be below a given maximum value. Lastly, constraint (9) restricts the maximal number of routes from the operator’s perspective.
2.1. Decision Variables and Input Data
The decision variables considered in the proposed model can be categorized as either continuous or integer variable. The continuous variables include the travel time between nodes and ; , invehicle travel time between nodes and ; , waiting time between nodes and ; , transfer time between nodes and ; and , round trip time of route , . On the other hand, the integer variables include the frequency of the route ; , the number of vehicles to be assigned to route ; , length of route ; and , the maximum flow occurring on the route , .
The same input data of the UTNDP used by the previous approaches in the literature are utilized in this study. There comprises (i) the road network structure available for the vehicle routes including the designation of nodes and the available edges in the network, (ii) the travel times of each edge in the network indicating how long it takes for a vehicle to travel between any two nodes of the road network, (iii) travel demand by transit between every OD pair (origindestination matrix), and (iv) the design parameters that include penalty per transfer, the minimum frequencies of bus service on each route, the bus seating capacity (assumed all buses are homogeneous in terms of seating capacity), and the maximum load factor allowed on any transit routes [23].
2.2. Passenger Assignment
When the transit route networks are constructed, the passengers travel demand should be distributed on the candidate routes to evaluate the fitness of the vector/particle as explained in Section 3.3. We [23] have described in detail how to determine the flow of passengers on paths selected by the passengers for every OD pair on the basis of the passenger assignment procedure in Baaj and Mahmassani [3]. We also assume that transit users can use a maximum of two bus lines to get to their destination (i.e., at most one transfer) similar to Mauttone and Urquhart [36]; Yan et al. [37]; and Nikolić and Teodorović [20]. The passenger assignment is briefly performed as follows:(i)The passengers look for a direct route (i.e., a route without transfer) that would serve a given pair of OD nodes.(ii)If direct route(s) is/are located, then the routes are screened on the basis of invehicle travel time such that any route with invehicle travel time greater than the smallest value by 50% is avoided.(iii)The demand is shared with the candidate routes that scale through the screening process.(iv)However, if it is not possible to find the direct route(s), the passengers will then try to look for a set of one transfer routes with the smallest travel time to meet the travel demand.(v)If one transfer route(s) is/are found, then the total travel time is evaluated for each possible path involving one transfer, and a screening process is invoked to determine the candidate paths similar to the zero transfer, in such a way that, all paths between OD pair whose total travel time is greater than the smallest value produced by any path by more than 10% are also avoided.(vi)If it is not possible to find any path, then the demand is regarded as unserved.
We assign randomly the initial service frequencies of all routes in the solution taking into account that the least frequency of service allowed on any route is one vehicle per hour. After the first iteration, frequencies are updated and demand is allocated. The iteration process is repeated until the stopping criterion is met: (i) a fixed number of iterations, (ii) the system converges to a solution [17], or (iii) no significant improvement in frequencies attained. Consequently, the four basic properties associated with a given route are determined: nodes list, service frequency, turnaround trip times, and linkflows list [23].
2.3. Bus Line Characteristics and Travel Time Calculation
It is essential to calculate the values of basic quantities (frequency of service, required number of buses, bus headway, etc.) that is associated with any considered solution after performing the passenger assignment (Section 3.3). These parameters can be expressed mathematically. The same mathematical expressions as in [23] are adopted in this study. Additionally, we included the following equations (equation (10)–(13)) to calculate the values of the maximum and average route headways.
The passenger waiting time:where is the bus headway at route
The bus headway at the route :
Maximum route headway:
Average route headway:where is the total route headway.
3. Hybrid DE—PSO for the UTNDP
3.1. Representation
Each vector (in DE) or particle (in PSO) is a candidate route set from the given transit network configuration. A list data structure is used to represent the vectors (particles) as shown in Figure 1, which are the permutations of the m transit lines with m transit routes and their associated frequencies. A transit line comprises an adjacent sequence of transit nodes stored in a transit route list, , and its associated frequency, , stored in a sublist . A separator “” is used to demarcate the set of transit lines. These frequencies are randomly initialized to individual routes of every vector/particle of the population (see [23]. For example, Figure 1, shows that route 1 has a service frequency of 5 vehicles per hour, while route 2 has a service frequency of 4 vehicles per hour in the sublist, etc.
3.2. Initialization
The route construction heuristic proposed by Mumford [38] is utilized to generate the initial population of vectors/particles. In this heuristic algorithm, several parameters need to be defined: (i) the number of routes per route set to be predefined by the user (i) the minimum and maximum length of individual routes, and (iii) the population size of the route sets. In the UTNDP, at least two routes emanating from a list of nodes should constitute a complete route set. The routes are constructed one at a time with vertices added one after the other. However, we observe that the construction heuristic does not always produce feasible route sets (i.e. not all the nodes are connected). Hence, we [22] recently proposed an improved subroute reversal repair mechanism (iSRR) to deal with the infeasible route sets. In the literature, some requirement concerning each typical route of the route set is specified by the operator. For instance, each typical route must contain a minimum of two or three nodes and a predefined number of maximum nodes fora feasible route set. Note that a fixed population size is maintained during the implementation of the proposed algorithm.
3.3. Fitness Evaluation
The solution is evaluated when the demand is assigned and the frequency determined through the assignment model (see Section 2.2). The objective functions (1) and (2) for an individual vector/particle are considered as the fitness functions. Consequently, the overall passenger cost, as well as the operator cost, could be evaluated to be the fitness value. The passenger cost consisting of two terms in objective function (1) is calculated by summing the travel cost and the unmet demand for every OD pair. The required fleet for each bus line can be obtained from equation (6). Finally, summing the buses of all routes gives the total fleet.
3.4. Mutation and Crossover
Two mutation schemes are used in the proposed hybrid algorithm: (i) identical point mutation, which is a modified form of the mutation operator by Ngamchai and Lovell [10] and (ii) mutation operator proposed in Kechagiopoulos and Belligiannis [19]. The identical point mutation is utilized in the DE to create a noisy random vector, , as follows: a random node that has at least two identical nodes from the random vector (route set) is selected. Then, two routes that contain the random node are selected and the substrings preceding the random node in the routes interchange their position between the two routes to create a noisy random vector (see Figure 2). The second mutation operator is utilized in the PSO to update the personal and global best of the swarm (see Figure 3). However, in our case, if there is a cycle after the modification of the particle, then the repair is attempted using iSRR. The steps of the mutation are as follows:(i)Select two routes randomly: one from the particle to be modified and the other from the personal or global best.(ii)A common node is searched between these pairs of routes.(iii)If the common node is found, then copy substring (route parts) from the personal or global best between the selected nodes to modify the current particle.(iv)If a cycle is identified, then repair with iSRR.
In the case of no common node being found between the pair of the selected routes, a new pair of routes is selected by the algorithm, again randomly, in an attempt to locate a common node for swapping parts of routes. However, if a candidate node is located, then the process terminates and a new particle is selected from the swarm to undergo modification.
The crossover is applicable only to the DE. Specifically, the uniform route crossover (interstring crossover operation) inspired by Beasley et al. [39] is utilized. This crossover is performed as follows: each subroute in the trial vector is generated by interchanging the corresponding subroute either from the target or noisy random vector on the basis of a randomly generated binary crossover mask whose length is equal to the number of routes in the solution. When there is a “1” in the crossover mask, the subroute is copied from the target vector, but when there is a “0” in the mask, the subroute is copied from the noisy random vector. The process is repeated with the target and noisy random vector exchanged to yield the second trial vector (see Figure 4). Note that for each pair of vectors, it is required to generate a new crossover mask [40]. It is advisable to ensure that the trial vectors are feasible; otherwise, attempt to repair the infeasible trial vectors with iSRR.
3.5. Replacement
The population for the next generation is constructed through elitism selection for the DE, such that the best vector between the target vector and the trial vector is selected for the next generation. This is to maintain the average fitness of the population [41]. For the PSO, the fitness value is updated for the personal and global best, provided the fitness of the current particle is better than the fitness of its personal best and global best after the application of the mutation operator (see Figure 3).
The proposed framework of the proposed hybrid DEPSO is structured as follows: the population (for DE) or swarm (for PSO) is initialized using a construction heuristic proposed by Mumford [38]. The DE algorithm (see Algorithm 1) is run for a finite number of generations, I (e.g., 20) to initialize the swarm for the PSO. The PSO is the main algorithm to search for the optimal solution. Note that it is not appropriate to use the velocity vector of the classical PSO due to the discrete nature of the problem. Instead, a new mutation scheme is utilized to update the personal and global best (see Figure 3). The proposed algorithm will only alternate once, and both of the objective functions will start with the same initial population that has been recorded earlier. At the end of each DEPSO iteration, the nondominated solution is accumulated and isolated by discarding all the dominated solutions. The detailed framework of the proposed DEPSO is provided in Algorithm 2.


4. Experimental Design
Two sets of experiments are conducted to assess the effectiveness of the proposed hybrid algorithms. In the first experiment, the computational results obtained by hybrid DEPSO and PSODE are compared to verify that the hybrid DEPSO yields better solutions. In the second experiment, the computational results obtained from the hybrid DEPSO algorithm are compared with other approaches reported in the literature. All experiments are performed on the wellknown benchmark Mandl’s Swiss network [42]. The road network consists of connected edges, edge weights (i.e., travel times of link), and feasible locations of the bus stop, and others (see Section 4.1). It is also assumed that each route must contain at least three nodes (stations) as in Nikolić and Teodorović [20] and Arbex and da Cunha [12]. The proposed algorithms are coded in Python 2.7.6.4 on a computer with 1.60 GHz Intel Core™i54200 CPU with 4.00 GB of RAM.
4.1. Benchmark Data
Mandl’s Swiss network [42] is a realistic transportation network in Switzerland (see Figure 5). The graph comprises 15 vertices, the shortest travel time between the two farthest vertices is 33 minutes, 21 undirected edges, and overall passenger travel demand is 15570. The highest travel demand between a node pair is 880 passengers. The network is small and dense. In this matrix, the nonzero demand is 82% of the node pairs. The network information for the transit operator includes travel demand and existing road network data. The demand is often represented by a symmetric OD matrix whose elements may include the time, day of the week, boarding point, alighting point, and potential coverage of transit passengers. The road network information includes connectivity of link, lengths of link, travel times of link, and feasible locations of a bus stop.
4.2. Parameter Configuration and Performance Metrics
A population of 30 vectors and a computation time of 200 iterations are performed for the computational experiment. Five scenarios are considered to investigate the nondominated solutions: route sets consisting of 4, 6, 7, 8, and 12 bus lines, respectively. To make a fair comparison with other approaches in the literature, the standard parameter configuration and the performance metrics used by Nikolić and Teodorović [20]; Arbex and da Cunha [12]; and Zhao et al. [21] are provided in Tables 1 and 2, respectively.


5. Results and Discussions
5.1. Comparison of Hybrid DEPSO and Hybrid PSODE
Computational experiments on benchmark Mandl’s Swiss network are conducted to assess the effectiveness of the two proposed hybrid algorithms: hybrid DEPSO and hybrid PSODE. The initial population of both hybrid algorithms are obtained by the route set generation heuristic by Mumford [38], embedded with iSRR mechanism. The iSRR mechanism is used to repair the resulting infeasible candidate route sets. For hybrid DEPSO, the DE algorithm is initially run for a finite number of generations, I (e.g., 20) to initialize the swarm population for the PSO. The PSO is the main algorithm to search for the optimal solution in hybrid DEPSO. On the other hand, for hybrid PSODE, the PSO algorithm is initially run for a finite number of generations, I (e.g., 20) to initialize the vector population for the DE. DE is the main algorithm to search for the optimal solution in hybrid PSODE. Route sets consisting of 4, 6, 7, 8, and 12 bus lines with a maximum of 8 nodes in each route are used. For each case, 10 independent runs are performed for 200 generations, each with a population size of 30. The computational results obtained by both hybrid algorithms are given in Table 3. The entries for Avg. represents the average result after the 10 runs. The transit route network configuration constructed by the proposed hybrid DEPSO and hybrid PSODE is provided in Tables 4 and 5, respectively. In both Tables 4 and 5, column 3 gives the fleet required per route of the route set. The values in the bracket give the overall fleet size required for the given transit route network configuration. Their associated passenger costs are given in column 4.



As shown in Table 3, hybrid DEPSO produced better results in terms of best d_{0}, C_{o}, C_{p}, T_{inv}, and T_{tr} for the five scenarios considered. However, the best values achieved for the parameter T_{wt} by hybrid DEPSO is similar to hybrid PSODE. In addition, hybrid DEPSO achieves lower average T_{wt} and T_{tr} as compared with hybrid PSODE. Furthermore, both algorithms produced solutions in which the demand is met with a maximum of only one transfer for the five scenarios considered. It can be observed that hybrid DEPSO is more capable of finding better solutions from the 10 independent runs than hybrid PSODE. A statistical ttest at 5% significant level is also been carried out for both hybrid algorithms. From Table 3 (column 9 and column 10), it can be observed that there is a significant difference in the parameters (except d_{un}) in all cases between both hybrid algorithms. Relatively, hybrid DEPSO produced better results than hybrid PSODE in terms of the point estimates of the parameters. Consequently, we adopt the hybrid DEPSO for comparison with other results reported in the literature.
5.2. Comparative Experiments of Hybrid DEPSO
5.2.1. Mandl’s Swiss Network
In this experiment, the proposed hybrid DEPSO algorithm is performed on the Mandl’s Swiss network and compared with the studies of Mandl [42]; Baaj and Mahmassani [3]; Shih and Mahmassani [43]; Bagloe and Cedar [44]; Nikolić and Teodorović [20]; Zhao et al. [21]; and Buba and Lee [23]. These articles used the same passenger assignment method and the same transfer penalties to estimate passengers’ costs. The comparative results on Mandl’s Swiss network with route sets comprising 4, 6, 7, 8, and 12 bus lines are reported in Table 6.
 
In Table 6, the proposed hybrid DEPSO (last column) outperformed the DE in terms of d_{0} for all cases considered. In addition, the DEPSO algorithm achieves the best result in terms of d_{0} and d_{1} compared with previous studies in the case of route sets containing 4 and 12 bus lines. In the case of route sets with 8 and 12 bus lines, the hybrid algorithm outperformed the previous results in terms of C_{p}. The hybrid DEPSO achieves the minimum transfer time of all passengers in four out of the five cases as compared with the previous studies. Furthermore, the proposed hybrid DEPSO is able to produce significant improvement in C_{o}, as reflected by the required fleet evidencing the high quality of solutions obtained. Zhao et al. [21] have the best value for d_{0} for cases of 6, 7, and 8 bus lines. However, they have some of the worse results for C_{p} in the literature.
As mentioned previously, the UTNDP is a complex problem, given the two conflicting objectives that influence each other: higher frequencies would improve the passenger cost, but this will increase the operator cost. Consequently, Arbex and da Cunha [12] observed that it is not enough to report the best solution found in terms of percentage of direct trips (d_{0}) and the one that minimizes the number of buses, as these two objectives are related. Hence, there is a set of solutions known as Pareto optimal solutions for the UTNDP, which represent the compromise solutions between the conflicting objectives [45]. The computational results of the Pareto front and the node sequence of the best routes found in the Pareto front by the proposed algorithm are given in Tables 7 and 8, respectively.


In Table 7, some important performance indicators, including average route headways, invehicle travel times, and the maximum route headways are determined. Headways are useful indications of how frequent public transit service is, and therefore the waiting time of the passengers. Note that for these indicators, we do not have existing results to compare with our results. The maximum and average route headway lie between (7.71, 12.32) and (5.724, 8.830) minutes, respectively. The proposed hybrid DEPSO solution produced more than one route choice for most OD pairs, as implied by the low values for average waiting times, which are significantly less than average route headways.
For the Mandl’s Swiss network, the lower bound on the invehicle travel time is 155,790 minutes [17]. The lower bound represents the theoretical minimum possible passenger time if passengers incurred no waiting time or transfer penalties and were able to take a vehicle that followed the shortest path from their origin to their destination. If operator cost is not an issue and there is no maximum frequency constraint, then it is possible to achieve a total generalized passenger cost that is arbitrarily close to this limit. Table 7 (column 4) provides the additional passenger time above the theoretical minimum for the hybrid DEPSO. In this table, the hybrid DEPSO produced solutions that are 18.65% to 28.98% closer to the theoretical minimum passenger time.
The Pareto frontier for the various solutions achieved from Table 8 with passenger cost on the vertical axis and the number of buses required on the horizontal axis is depicted in Figure 6. Efficient frontier solutions are highlighted in dark crosses, whereas past literature results are authorlabeled. It can be observed that the passenger cost for the best compromising solutions is lower than the previously published results, evidencing the high quality of the solutions obtained by the proposed hybrid DEPSO algorithm.
5.2.2. Rivera City Network
To demonstrate the scalability of the proposed hybrid DEPSO, we applied the proposed algorithm to solve a realsize public transportation network of Rivera City (see Figure 7). This case study belongs to a mediumsize city of 65,000 citizens in Northern Uruguay. It comprises 84 nodes, 143 arcs, and 378 OD pairs. Because we do not have the optimal number of transit routes for the network under study, therefore, we considered the minimum number of routes is 40 and the maximum number of routes is 60, whereas the minimum and the maximum numbers of nodes per route are 10 and 30 base on preliminary run. In addition, all of the parameters in Table 1 are utilized in this methodology, except the maximum allowable fleet size, which to the best of our knowledge is not available.
The hybrid approach produced solutions with maximum and minimum fleet size (operator cost) of 126 and 86 buses, respectively (see Table 9). This implies a variation of 32% in the cost between the two extreme solutions; however, it leads to an increase in the passenger cost by 34%. In Figure 8, the inherent conflict between the network bus fleet size and passenger cost that makes it difficult to find a unique optimal solution is depicted. The different tradeoffs between the competing objectives are established and it is left for the decisionmaker to attain a biased solution either to the passenger or the operator. It is important to note that it is difficult (impossible) to minimize passenger and operator costs at the same time. Furthermore, the overall passenger travel demand is satisfied with at most one transfer. The maximum and average route headway lie between (2.05, 8.20) and (1.21, 6.31) minutes, respectively. The proposed hybrid DEPSO solution produced more than one route choice for most OD pairs, as implied by the low values for average waiting times, which are significantly less than average routes headways.

6. Conclusions and Future Research
In this article, the UTNDP is tackled through a hybrid DEPSO algorithm taking into consideration the interest of both passenger and operator. The problem is modeled as a multiobjective optimization problem with the aim to simultaneously optimize the route configuration and the associated service frequencies that produce a minimum total passengers’ and operators’ costs while satisfying a given set of constraints. The proposed hybrid DEPSO algorithm is designed to adapt the strength of both approaches to achieve better quality solutions. Experimental studies are conducted to evaluate the relative efficiency of hybrid DEPSO and hybrid PSODE where both algorithms are capable of determining efficient solutions. The hybrid DEPSO algorithms adopted for further comparative studies as it obtained better quality solutions than the hybrid PSODE. In addition, comparisons are made between the proposed hybrid DEPSO algorithm with other approaches from the literature using the benchmark Mandl’s Swiss network. Furthermore, the proposed hybrid DEPSO algorithm is applied to a public transport system of the Rivera City network. The computational results of the proposed hybrid DEPSO algorithm improve over the results obtained from the previous studies. As a multiobjective optimization problem, we have shown that the proposed hybrid DEPSO produces a diverse set of nondominated solutions. The results also dominated solutions from previous studies in the literature. The application of the proposed hybrid DEPSO to larger and more realistic problem instances with heterogeneous buses will be the direction of future research.
Data Availability
The Mandl’s Swiss network datasets supporting this study are from previously reported studies and datasets, which have been cited. The processed data are publicly available at https://users.cs.cf.ac.uk/C.L.Mumford/Research%20Topics/UTRP/Outline.html. The Rivera City network and the travel demand matrices can be can downloaded via this link: http://www.fing.edu.uy/∼mauttone/tndp.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
This research was supported by the Fundamental Research Grant Scheme (FRGS) (0101161867FR), Ministry of Higher Education, Malaysia.
References
 L. Fan, C. L. Mumford, and D. Evans, “A simple multiobjective optimization algorithm for the urban transit routing problem,” in Proceedings of the IEEE Congress on Evolutionary Computation, pp. 1–7, Trondheim, Norway, May 2009. View at: Publisher Site  Google Scholar
 A. Ceder and N. H. M. Wilson, “Bus network design,” Transportation Research Part B: Methodological, vol. 20, no. 4, pp. 331–344, 1986. View at: Publisher Site  Google Scholar
 M. H. Baaj and H. S. Mahmassani, “An AIbased approach for transit route system planning and design,” Journal of Advanced Transportation, vol. 25, no. 2, pp. 187–209, 1991. View at: Publisher Site  Google Scholar
 P. Chakroborty, “Genetic algorithms for optimal urban transit network design,” ComputerAided Civil and Infrastructure Engineering, vol. 18, no. 3, pp. 184–200, 2003. View at: Publisher Site  Google Scholar
 Y.J. Lee and V. R. Vuchic, “Transport network design with variable demand,” Journal of Transportation Engineering, vol. 131, no. 1, pp. 1–10, 2005. View at: Publisher Site  Google Scholar
 R. Z. Farahani, E. Miandoabchi, W. Y. Szeto, and H. Rashidi, “A review of urban transportation network design problems,” European Journal of Operational Research, vol. 229, no. 2, pp. 281–302, 2013. View at: Publisher Site  Google Scholar
 F. Zhao and X. Zeng, “Simulated annealinggenetic algorithm for transit network optimization,” Journal of Computing in Civil Engineering, vol. 20, no. 1, pp. 57–68, 2006. View at: Publisher Site  Google Scholar
 L. Liu, P. Olszewski, and P.C. Goh, “Combined simulated annealing and genetic algorithm approach to bus network design,” in Transport Systems Telematics, pp. 335–346, Springer, Berlin, Germany, 2010. View at: Google Scholar
 W. Y. Szeto and Y. Wu, “A simultaneous bus route design and frequency setting problem for Tin Shui Wai, Hong Kong,” European Journal of Operational Research, vol. 209, no. 2, pp. 141–155, 2011. View at: Publisher Site  Google Scholar
 S. Ngamchai and D. J. Lovell, “Optimal time transfer in bus transit route network design using a genetic algorithm,” Journal of Transportation Engineering, vol. 129, no. 5, pp. 510–521, 2003. View at: Publisher Site  Google Scholar
 W. Fan and R. B. Machemehl, “Optimal transit route network design problem with variable transit demand: genetic algorithm approach,” Journal of Transportation Engineering, vol. 132, no. 1, pp. 40–51, 2006. View at: Publisher Site  Google Scholar
 R. O. Arbex and C. B. da Cunha, “Efficient transit network design and frequencies setting multiobjective optimization by alternating objective genetic algorithm,” Transportation Research Part B: Methodological, vol. 81, pp. 355–376, 2015. View at: Publisher Site  Google Scholar
 P. Gundaliya, P. Shrivastava, and P. Dhingra, “Model for simultaneous routing and scheduling using genetic algorithm,” European Transport, vol. 6, no. 16, pp. 10–19, 2000. View at: Google Scholar
 V. M. Tom and S. Mohan, “Transit route network design using frequency coded genetic algorithm,” Journal of Transportation Engineering, vol. 129, no. 2, pp. 186–195, 2003. View at: Publisher Site  Google Scholar
 S. M. M. Amiripour, A. S. Mohaymany, and A. Ceder, “Optimal modification of urban bus network routes using a genetic algorithm,” Journal of Transportation Engineering, vol. 141, no. 3, Article ID 04014081, 2015. View at: Publisher Site  Google Scholar
 M. Owais and M. K. Osman, “Complete hierarchical multiobjective genetic algorithm for transit network design problem,” Expert Systems with Applications, vol. 114, pp. 143–154, 2018. View at: Publisher Site  Google Scholar
 J. J. Blum and T. V. Mathew, “Intelligent agent optimization of urban bus transit system design,” Journal of Computing in Civil Engineering, vol. 25, no. 5, pp. 357–369, 2011. View at: Publisher Site  Google Scholar
 M. Bagherian, S. Massah, and S. Kermanshahi, “A swarm based method for solving transit network design problem,” in Proceedings of the Australasian Transport Research Forum 2013, Brisbane, Australia, 2013. View at: Google Scholar
 P. N. Kechagiopoulos and G. N. Beligiannis, “Solving the urban transit routing problem using a particle swarm optimization based algorithm,” Applied Soft Computing, vol. 21, pp. 654–676, 2014. View at: Publisher Site  Google Scholar
 M. Nikolić and D. Teodorović, “A simultaneous transit network design and frequency setting: computing with bees,” Expert Systems with Applications, vol. 41, no. 16, pp. 7200–7209, 2014. View at: Publisher Site  Google Scholar
 H. Zhao, W. Xu, and R. Jiang, “The Memetic algorithm for the optimization of urban transit network,” Expert Systems with Applications, vol. 42, no. 7, pp. 3760–3773, 2015. View at: Publisher Site  Google Scholar
 A. T. Buba and L. S. Lee, “Differential evolution with improved subroute reversal repair mechanism for multiobjective urban transit routing problem,” Numerical Algebra, Control & Optimization, vol. 8, no. 3, pp. 351–376, 2018. View at: Google Scholar
 A. T. Buba and L. S. Lee, “A differential evolution for simultaneous transit network design and frequency setting problem,” Expert Systems with Applications, vol. 106, no. 208, pp. 277–289, 2018. View at: Publisher Site  Google Scholar
 E. RuanoDaza, C. Cobos, J. TorresJimenez, M. Mendoza, and A. Paz, “A multiobjective bilevel approach based on globalbest harmony search for defining optimal routes and frequencies for bus rapid transit systems,” Applied Soft Computing, vol. 67, pp. 567–583, 2018. View at: Publisher Site  Google Scholar
 S. B. Jha, J. K. Jha, and M. K. Tiwari, “A multiobjective metaheuristic approach for transit network design and frequency setting problem in a bus transit system,” Computers & Industrial Engineering, vol. 130, pp. 166–186, 2019. View at: Publisher Site  Google Scholar
 D. Canca, A. DeLosSantos, G. Laporte, and J. A. Mesa, “A general rapid network design, line planning and fleet investment integrated model,” Annals of Operations Research, vol. 246, no. 12, pp. 127–144, 2014. View at: Publisher Site  Google Scholar
 D. Canca, A. deLosSantos, G. Laporte, and J. A. Mesa, “An adaptive neighborhood search metaheuristic for the integrated railway rapid transit network design and line planning problem,” Computers & Operations Research, vol. 78, pp. 1–14, 2017. View at: Publisher Site  Google Scholar
 F. LópezRamos, E. Codina, Áa Marín, and A. Guarnaschelli, “Integrated approach to network design and frequency setting problem in railway rapid transit systems,” Computers & Operations Research, vol. 80, pp. 128–146, 2017. View at: Publisher Site  Google Scholar
 G. GutiérrezJarpa, G. Laporte, V. Marianov, and L. Moccia, “Multiobjective rapid transit network design with modal competition: the case of Concepción, Chile,” Computers & Operations Research, vol. 78, pp. 27–43, 2017. View at: Publisher Site  Google Scholar
 F. Zhao and X. Zeng, “Optimization of user and operator cost for largescale transit network,” Journal of Transportation Engineering, vol. 133, no. 4, pp. 240–251, 2007. View at: Publisher Site  Google Scholar
 W. Y. Szeto and Y. Jiang, “Hybrid artificial bee colony algorithm for transit network design,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2284, no. 1, pp. 47–56, 2012. View at: Publisher Site  Google Scholar
 W. Y. Szeto and Y. Jiang, “Transit route and frequency design: bilevel modeling and hybrid artificial bee colony algorithm approach,” Transportation Research Part B: Methodological, vol. 67, pp. 235–263, 2014. View at: Publisher Site  Google Scholar
 T. Hendtlass, A Combined Swarm Differential Evolution Algorithm for Optimization Problems, SpringerVerlag Berlin Heidelberg, Berlin, Germany, 2001.
 M. Pant, R. Thangaraj, and A. Abraham, “Depso: a new hybrid metaheuristic for solving global optimization problems,” New Mathematics and Natural Computation, vol. 7, no. 3, pp. 363–381, 2011. View at: Publisher Site  Google Scholar
 A. Sedki and D. Ouazar, “Hybrid particle swarm optimization and differential evolution for optimal design of water distribution systems,” Advanced Engineering Informatics, vol. 26, no. 3, pp. 582–591, 2012. View at: Publisher Site  Google Scholar
 A. Mauttone and M. E. Urquhart, “A multiobjective metaheuristic approach for the transit network design problem,” Public Transport, vol. 1, no. 4, pp. 253–273, 2010. View at: Publisher Site  Google Scholar
 Y. Yan, Z. Liu, Q. Meng, and Y. Jiang, “Robust optimization model of bus transit network design with stochastic travel time,” Journal of Transportation Engineering, vol. 139, no. 6, pp. 625–634, 2013. View at: Publisher Site  Google Scholar
 C. L. Mumford, “New heuristic and evolutionary operators for the multiobjective urban transit routing problem,” in Proceedings of the 2013 IEEE Congress on Evolutionary Computation, pp. 939–946, Cancun, Mexico, June 2013. View at: Publisher Site  Google Scholar
 D. Beasley, D. R. Bull, and R. R. Martin, “An overview of genetic algorithms: part 2, research topics,” University Computing, vol. 15, no. 4, pp. 170–181, 1993. View at: Google Scholar
 A. T. Buba and L. S. Lee, “Differential evolution for urban transit routing problem,” Journal of Computer and Communications, vol. 4, no. 14, pp. 11–25, 2016. View at: Publisher Site  Google Scholar
 A. P. Engelbrecht, Computational Intelligence. An Introduction, John Wiley & Sons Ltd, London, UK, 2007.
 C. E. Mandl, “Evaluation and optimization of urban public transportation networks,” European Journal of Operational Research, vol. 5, no. 6, pp. 396–404, 1980. View at: Publisher Site  Google Scholar
 M. Shih and H. S. Mahmassani, A Design Methodology for Bus Transit Networks with Coordinated Operations. No. SWUTC/94/600161, University of Texas, Center for Transportation Research, Austin, TX, USA, 1994.
 S. A. Bagloee and A. Ceder, “Transit network design for actual size networks,” Transportation Research Part B, Methodological, vol. 45, no. 10, pp. 1787–1804, 2011. View at: Google Scholar
 E. G. Talbi, Metaheuristics from Design to Implementation, John Wiley & Sons, Hoboken, NJ, USA, 2009.
Copyright
Copyright © 2019 Ahmed Tarajo Buba and Lai Soon Lee. 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.