An Improved Adaptive Large Neighborhood Search Algorithm for the Heterogeneous Customized Bus Service with Multiple Pickup and Delivery Candidate Locations

In order to tackle the congestion and environmental issues, customized bus services are proposed and deployed in metropolitan areas. As emerging public transportation services, customized bus services bring passengers more convenience and accessibility. Besides, conventional customized bus services generally organize homogeneous fleet and single location selection to passengers. In this paper, to enhance the mobility and flexibility of customized buses and increase companies’ profit, we propose a new form of customized bus service with heterogeneous fleets and multiple candidate locations. First, a mixed-integer programming model (MIP) is developed to describe the customized bus problem. Comparedwith the conventionalmodel, the proposedMIP is involved in the case of one passenger withmultiple candidate pickup or delivery locations and can be solved by GUROBI on small scale, quickly and efficiently. Second, an improved adaptive large neighborhood search algorithm (ALNSip) is utilized to address the large-scale problemmore efficiently. Time slack calculationmethod is then designed to optimize vehicle timetables, which provides stable and excellent performance for searching feasible solutions. In addition, we propose two inserting operators to deal with the problem with multiple candidate locations and analyse its influence on the results. Finally, we test the performance of the proposed model and algorithm on the numerical experiments. And they are verified the effectiveness and implication in a small-scale case on a simplified Sioux waterfall network and a large-scale problem in Beijing, China.(e result shows that ALNSip outperforms other algorithms in searching for more satisfying solutions with higher efficiency. However, the GUROBI solver can obtain the solution to small-scale problems within a shorter time than ALNSip. Furthermore, it can be suggested that the heterogeneous fleets service with multiple candidate locations is helpful to facilitate collaboration among vehicles and optimize pickup and delivery routes in consequence.


Introduction
e increase in the number of private cars has imposed an impact on potential traffic jams and environmental issues, and in the meanwhile, the increasing number of people has grown interest in enjoying personalized travel services. erefore, many transportation service providers are gradually focusing their attention on seeking new services to attract more public transport users and increase passengers' accessibility. In recent years, internet-based customized buses service has become more popular in the market due to its higher accessibility and quicker response, especially in congested metropolitan areas. Since Kirby and Bhatt analysed and discussed in detail about subscribing to the bus service in 1974 [1], the customized bus system (CB) has been continuously developed into a demand-responsive transportation service. Nowadays, some transportation companies have constructed mobile application platforms for customized bus services to collect detailed travel information [2,3]. With the spatial and temporal characteristics of travel requests taken into considerable consideration, the customized bus service platform can aggregate similar requests and design bus routes. To the best of our knowledge, some customized bus services have been deployed in physical road networks of metropolitan areas such as Beijing.
According to different characteristics of travel requests, customized bus services can be usually divided into four major categories: "Customized commuter bus," "Customized school bus," "Customized business bus," and "Customized community bus," and other kinds of customized feeder or shuttle bus systems [2]. is paper focuses on the customized commuter bus, which offers commute shuttle service and takes customers from the residential community to the workplace. is is also the most important CB system. It is noted that there are some unique service features about commuters. ey usually book relatively fixed demands and change them in different weeks or months. erefore, the customized bus services will be updated according to new requests in each optimization time (e.g., one month and one week) [4]. e detailed operating procedure of the customized bus service is presented in Figure 1. Each passenger needs to input his/her request through the network and platform to complete the service-design process.
is paper mainly focuses on the process to plan and deploy new routes based on fail-to-match requests in the requests pool. If requests cannot be matched with the existing routes, these requests will be put into the requests pool. When the number of booked people meets the minimum load requirement, the platform will plan new routes deployed in the physical road network based on the requests in the requests pool.
ere are several notable features of customized bus problems. As can be seen from the above operating procedure, compared to the traditional bus routing problem, customized bus problems pay more attention to individual passenger request information. Bus route and specific timetables are updated and designed by considering more about companies' profit and various characteristics of passenger demands, including origindestination pairs and time windows. In the field of operations research, the problem regarding passengers' pickup and delivery locations with time windows can be classified as a special vehicle planning problem [5]. With the consideration of all constraints, it is worthwhile but challenging to design reasonable routes and specific timetables. Furthermore, it is also of great significance to enrich the personalized public transportation services and provide highly reliable and accessible trip services for commuters living in congested metropolitan areas. e remainder of this paper is organized as follows. e next section describes an overview of earlier research. Section 3 describes the MIP model. In Section 4, we describe ALNS ip in detail. Section 5 reports the test of the MIP and ALNS ip on numerical experiments. Finally, Section 6 summarizes the research content and future research directions of the paper.

Literature Review
Customized bus service can be regarded as a request-oriented public transportation service based on its unique service features. In earlier research, on-demand transportation was described as a set of vehicle routing problems with pickup and deliveries (VRPPD). Nagy et al. [6] divided the previous studies about the problem into three main categories: simultaneous pickups and deliveries, mixed pickups and deliveries, and delivery-first, pickup-second. In subsequent research, VRPPD is gradually evolved into different versions due to the specific problems and traffic conditions. When the pickup and delivery locations with time windows are determined, VRPPD can be seen as pickup and delivery problems with time windows (PDPTW). In general, the exact algorithm is a major way to solve multivehicle pickup and delivery problems in most studies [7][8][9]. Ropke et al. [10] proposed a new branch-and-cut-and-price algorithm to tackle the problem, which dynamically introduced valid inequalities and heuristic algorithm to enhance algorithm performance. Wang et al. [11][12][13][14] constructed related optimization models and proposed algorithms based on state-space-time network to effectively solve some pickup and delivery problems.
With the development of mobile technology, passengers can book vehicles and pay for trips through smartphones or platforms. In the operations research, this type of transportation service can be classified as dial-a-ride problems (DARP). Compared to VRPPD, DARP offers more unique characteristics, including depots, time window, vehicles capacity, ride time, and route duration. Stein [15] proposed the first DARP model and gave the bounds of both static and dynamic models. For static and dynamic DARP versions, Psaraftis et al. [9] developed a dynamic programming algorithm for the problem with a single vehicle. For the multiobjective problem in DARP, Jorgensen et al. [16] considered the total weight of multiple objects and proposed a Granular Tabu Search algorithm to solve the static DARP. In addition, Garaix et al. [17] and Schilde et al. [18] considered optimizing multiobjective problems by priority. And Paquette et al. [19] used the Pareto boundary to solve it. e original application of DARP was a nonprofit dial-a-ride system service (DAR) for the elderly and the disabled. Relevant constraints include ride and waiting time, pickup or delivery with time windows, vehicles capacity, and fleet systems [20][21][22][23]. However, due to market demand and emerging communication technology, DARP has been applied in the field of public transportation. Based on the on-demand transportation system developed in rural France, Garaix et al. [17] designed a freeriding system that maximized customer occupancy and used the column generation way to deal with it. e problem where not all passengers are collected in advance is classified as dynamic on-demand transportation services [5]. Pillac et al. [24] introduced the concept of dynamic problem and reviewed the solution methods about it. Alonso-Mora et al. [25] paid attention to the dynamic ride-sharing problem with high-capacity and proposed a pairwise shareability graph to compute feasible trips. Studies about customized bus service, one type of demandresponsive transport service, are also gradually increasing. Liu et al. [2] introduced the development and application of customized bus services in China. To solve customized bus servicedesign problems, Tong et al. [4] proposed a multicommodity optimization model and provided a Lagrangian decompositionbased solution algorithm to solve the routing problems.
Furthermore, in terms of solving route planning problems, the neighborhood search algorithm always has a high performance [26][27][28][29][30][31]. To improve the performance of neighborhood search algorithm, Ropke and Pisinger [32] proposed an adaptive large neighborhood search algorithm for the first time (ALNS). In the subsequent studies, it is found that ALNS has an outstanding performance in solving large-scale route planning problems with NP-hard, and it has been applied to different fields. To stack with the pollution routing problem (PRP), Demir et al. [33] improved the ALNS algorithm and updated some mechanisms, applicable to the problem. Polat et al. [34] proposed a mixed-integer mathematical optimization model and design a perturbation-based variable neighborhood search heuristic to solve the vehicle routing problem with simultaneous pickup and delivery with time limit (VRPSPDTL). To study the share-a-ride problem, Li et al. [35] proposed an adaptive large neighborhood search method with time slack. Yu et al. [36] and Yu et al. [37] developed an adaptive large neighborhood search heuristic to solve green vehicle routing problems in large-scale instances. In addition, ALNS is also widely used in the vehicle routing problem with drones, a supply chain, carsharing service, train stop planning, timetabling, and other fields [38][39][40][41]. Furthermore, there are also some hybridization and improvement from ALNS and other algorithms [42][43][44].
In the paper, to further improve passengers' accessibility and increase bus operators' profits, we consider adding two strategies (Strategy 1 and Strategy 2) into customized bus services.
(I) Strategy 1: the heterogeneous fleets with different operating cost and vehicle capacity can be dispatched to offer customized bus service (II) Strategy 2: passengers can choose multiple candidate locations as their pickup and delivery locations e existing literature still lacks relative research on the heterogeneous customized bus service problem with multiple pickup and delivery candidate locations. With the comprehensive consideration about the new service form, this paper builds a new MIP model based on the classical DARP model proposed by Cordeau [20] to address the practical and theoretical challenges. At the same time, to effectively solve large-scale problems, we introduce and improve ALNS (ALNS ip ) proposed by Li et al. [35] by designing corresponding operators (MI.1 and MI.2) to address new challenges with multiple candidate locations. In addition, given the complexity of the problem with time window, a time slack calculation method is utilized to adjust vehicle arriving time during searching feasible solutions.
In summary, this study has several contributions in theory and application: (1) constructing a new bus service with heterogeneous fleets and multiple candidate locations to improve conventional service quality and facilitate collaboration among different buses, (2) formulating relative optimization model with considering the different time windows, capacities, passenger assignment, and location selection to maximize operator's profit, (3) designing an improved heuristic algorithm based on ALNS, including a time slack calculation method and two inserting operators to address the large-scale problem with high efficiency.

Mathematical Model
e PDPTW and the DARP can be modelled using classical formulation (see, e.g., Cordeau [20] and Dumas et al. [7]). It is usually utilized to optimize the problem where each request is represented by a set of unique pickup location and delivery location. In this paper, to further tackle the challenges with multiple depots, heterogeneous fleets and multiple candidate locations, we propose a MIP model to deal with the problem, where the state node concept is utilized to describe the relationship between every passenger and their pickup and delivery locations.

Problem Statement.
We define the problem on a direct graph G � (S, Λ), where S is the set of all state nodes and Λ is the set of all arcs connecting state nodes. Each state node (p, i) in S � (p, i)|p ∈ P, i ∈ V means the passenger p is picked up and delivered at the physical location i. P, V, and K mean, respectively, the set of passenger requests, the set of physical locations, and the set of customized buses. For each passenger p 0 , (p 0 , i) ∈ S p is the set of state nodes for passenger p 0 picked up, and the same for (p 0 , i) ∈ S d is the set of state nodes for passenger p 0 delivered. Λ � (p, i), (q, j)|(p, i) ∈ S, (q, j) ∈ S indicates one directed arc from one state node (p, i) to another state node (q, j), where bus moves from location i to location j to offer service for passenger p and passenger q. Although one passenger may have multiple pickup and delivery locations, he/she must be picked up or delivered at a unique location.
To guarantee bus operators' profit, the number of passengers the bus serves must meet the minimum load constraint. As a result, if the loading capacity of bus k is Cap u k , the number of passengers on the bus must be more than or equal to Cap l k . Given a number of service features from customized buses, some unreasonable passenger requests may be rejected. e rejected passengers in the model are served by the vehicle k ∈ K V , where K V is the set of virtual vehicles. Expensive cost from virtual vehicles is the penalty cost for unserved customers. e remaining passengers are served by the buses k ∈ K L ∪ K S , where K L is the bus set with large capacity.
Pretreatment: before solving the model, it is necessary to find all potential state nodes S and construct the arcs Λ among these state nodes.

Related Definitions and Optimization Model.
To formulate DARP with multiple candidate locations into a mathematical model, all parameters used in the model are defined in Table 1.
e total income mainly refers to the tickets income. e cost includes the following three parts: the fixed cost of each bus dispatched, fuel cost on the arcs, and penalty cost for unserved passengers. (1) (q,j)∈S(p,i)∈S (q,j)∈S(p,i)∈S Set of all arcs connecting state nodes; Λ � (p, i), (q, j)|(p, i) ∈ S, (q, j) ∈ S indicates one directed state arc from one state node (p, i) to another state node (q, j), where bus moves from location i to location j to offer service for passenger p and passenger q K Set of all customized bus: State node means the passenger is picked up at the location: State node means the passenger is delivered at the location, Total fuel cost Customized bus k departs from depots or return depots y k (p,i),(q,j) e bus k moves from the state node (p, i) to the state node (q, j) Journal of Advanced Transportation e objective function (1) maximizes total profits. Constraints (2)-(5) are ticket income, total departure cost, penalty cost, and fuel cost, respectively. And constraints (3) and (4) can be represented by constraint (20).
For vehicles, constraints (6) and (7) ensure that each vehicle travels at most once, and each one leaves start depots and enters start depots. For some customers offering multiple locations as their candidate pickup or delivery locations, he/she can only be picked up or delivered at one location in the end, and the platform also will choose one location from their multiple pickup or delivery candidate locations. erefore, constraints (8) and (9) ensure that all passengers are served and each passenger is picked up or delivered at a unique location. In addition, constraint (10) indicates that one person is just picked up and delivered by the same vehicle. For network flow, constraints (11) and (12) enforce traffic flow balance. If vehicle k arrives in location i to serve one person, it will leave from location i to location j to serve another one. Constraint (13) guarantees that the total number of served passengers cannot exceed the capacity of vehicle. To achieve the profit goal, constraint (14) ensures a minimum number of passengers in the vehicle. Constraints (15) and (16) enforce each passenger to be picked up or delivered within he/her time windows. Constraint (17) is subtour elimination constraint. Finally, constraints (18) and (19) define the binary decision variables. If the vehicle k moves from the state node (p, i) to the state node (q, j), then y k (p,i),(q,j) � 1; otherwise, y k (p,i),(q,j) � 0. If one vehicle departs from start depots or returns to end depots, then

An Adaptive Large Neighborhood Search Algorithm for Customized Bus Service System
In this section, we describe ALNS ip and its application to the problem in great detail. e ALNS ip constructs the initial feasible solution using the greedy algorithm represented in Section 4.1. It uses several competing removal and insertion algorithms and chooses a reasonable set of operators by the score of each operator during the search. e removal algorithms are described in Section 4.3, and the insertion algorithms are described in Section 4.4. Furthermore, we propose a time slack calculation strategy to adjust timetables for searching more feasible solutions.

Initial Solution Construction.
A greedy insertion algorithm is used to search for a feasible initial solution. Requests in the set P are inserted one by one in a random order into routes. e entire process can be shown in the pseudocode of Algorithm 1.

ALNS ip Process.
To diversify the search, the ALNS ip selects suitable removal and insertion operators by a roulette wheel mechanism at every iteration. We define every UpI times as a cycle. e probabilities of choosing operator i in the i is the weight of the operator i in the k th cycle. For every UpI iteration, weights are recalculated using the scores obtained as formulation are the overall score and selected times of the operator i in k − 1 th cycle. If the operator i is chosen, its score is increased by λ q (q � 1, 2, 3 or 4), as shown in Table 2, according to its performance at the current iteration. For every UpI iteration, all scores are reset to zero. Finally, the entire process of ALNS ip is completed under U iterations, as shown in Algorithm 2.
Some parameters are shown in Table 2 by using the tuning strategy proposed by Ropke and Pisinger [32] and Li et al. [35] to obtain values. When just one parameter is adjusted, the rest will be fixed. Finally, the setting with the best average behaviour is chosen. e acceptance criterion adopts the simulated annealing criterion shown as Algorithm 2, where f ′ (S ′ ) is the objective value of the new solution, f best current (S) is the objective value of the current solution, and f best global (S) is the objective value of the global optimal solution. e probability that we accept new solution is ρ � e − (f′(S′)− f best current (S))/T , where T � T 0 ε is the current temperature, ε is the cooling rate, and T 0 is the initial temperature.

Removal Algorithms (RA).
In each iteration, a removal operator is selected from RA and removes the total number of passenger requests n p � u · P in the solution to break the current solution, where u is the random removal parameter. ese requests removed are added to the perturbation set RP. e section has designed five operators, the first three are chosen from Ropke and Pisinger [32], the fourth operator is adapted from Demir et al. [33], and the last one is motivated by Cordeau and Laporte [5].
R1-Random Removal: the operator randomly removes n p requests from the current solution put them into the set P. R2-Worst-Removal: the operator calculates the Δf p value of all passengers, where Δf p denotes the change in objective value incurred by removing passenger request p from current position. en, n p requests with the greatest value Δf p are removed and put into the set RP. R3-Shaw Removal: to diversify the search, the operate measures the similarity level R ij between passenger i and passenger j through the correlation function (2) Repeat (3) Construct a new route; (4) Insert one passenger into the route, suitable for all constraints; (5) Remove the passenger from P; (6) If the number of passengers in the current bus exceeds its capacity then (7) Return to step 3; (8) Else (9) Return to step 4; (10) End if (11) Until there is no new route can be constructed, remain passengers in P cannot be served and be refused; (12) Get all routes and passengers refused as initial solution S; ALGORITHM 1: Initial solution.
(1) Input initial solution S, the objective values of initial solution S as initial values of f best current (S) and f best global (S); Select a couple of removal and insertion operators using a roulette wheel mechanism to generate new feasible solution S′; Score of removal and insertion operators is increased by and the solution has not been accepted before then (8) Score of removal and insertion operators is increased by λ 2 ; f best current (S) ⟵ f′(S′); S←S′; (9) Else if f′(S′) ≥ f best current (S) and the solution has been accepted before then (10) Score of removal and insertion operators is increased by λ 3 ; f best current (S) ⟵ f′(S′); S←S′; (11) Else if random value ρ′ ≤ ρ then (12) Score of removal and insertion operators is increased by λ 4 ; f best current (S) ⟵ f′(S′); S←S′; (13) End if (14) if mod(i, UpI) � 0 then (15) Update weight of each operator; (16) End if (17) i ⟵ i + 1; (18) Get the efficient feasible solution S; ALGORITHM 2: ALNS ip algorithm procedure.  [32], the fourth one is motivated by Cordeau and Laporte [5], the fifth is inspired by Demir et al. [33], while the last two are the insertion operators designed for the problem with multiple pickup and delivery locations in the paper.
I1-Greedy Insertion: the requests in the set RP are sequentially removed from RP set and inserted into the best location. I2-Regret-2-Insertion: to avoid the myopia of local optimization from greedy algorithm, the operator calculates the regret-k value, where k � 2. And then it sequentially inserts the request with the greatest regret value. I3-Regret-3-Insertion: this operator is similar to I2 and makes k � 3. I4-One by one Insertion: the operator is similar to I1. e difference is that I4 inserts one request into the best location after removing it, rather than reinserting one by one after removing all requests. I5-TimeWindows Insertion: to tackle with the impact of the time feature from the problem, I5 gets the similarity level of all removed passengers in terms of time window, that is, the value of time window difference, |l 0 (p) − e 0 (p)| + |l 1 (p) − e 1 (p)|. And then, it sequentially inserts the request with the greatest value.
For the passenger with multiple locations, we need to combine one of I1-I5 and one of the last two operators to complete the insertion of the passenger. For this, one of the following two operators is randomly selected: MI.1-Random insertion based on multiple locations: MI.1 principle randomly selects a candidate location as the current pickup or delivery location of the passenger and then inserts the passenger using I1-I5. MI.2-One-by-one insertion based on multiple locations: MI.2 principle compares the performance of all candidate locations using I1-I5 separately and selects the candidate location with the best performance as the current pickup or delivery location of the passenger and inserts it. Figure 2, if we break and repair the solution using RA and IA operators without adjusting timetables, it is possible to fail to search more feasible insertion solutions. erefore, we design the instance to simulate the issue, where the maximum ride time of each passenger is 40 and the service time is 0. Time window, pickup, and delivery locations of two passengers are presented in Figure 2.

Time Slack Calculation Strategy. As seen in
If we just use Algorithm 3 without any ways to adjust timetables, it can be found that the ride time of passenger 1 is more than 40 and passenger 2 cannot be inserted into the route as shown in Figure 2. After using Algorithm 4 to change the timetables, passenger 2 can be inserted in the position and the insertion solution is feasible, as shown in Figure 3. erefore, it is very essential to adjust timetables to further improve the solution during the search. e time slack calculation strategy we proposed can test the feasibility of any insertion cases in time constraints. e symbols that may be used in Section 4.5 are shown in Table 3.
As shown in Figure 2, the greedy algorithm just guarantees time window constraint but cannot consider the impact of the maximum ride time. In addition, it may also cause more waiting time.
T tw i is the adjustable time bus can delay its arrival at the location i according to the timetable and time window: T rd i is the adjustable time bus can delay its arrival at the location i within the maximum ride time constraint:

MTS '
′ i means the adjustable time vehicle can delay its arrival at the location i, where its value may affect time constraints at other locations.

Stopping Criteria.
ere are two types of stopping criteria used to avoid local optimization: (i) When the current iteration exceeds the maximum iteration, the algorithm operation ends and the result is output (ii) When the temperature T � T 0 ε is less than 0.5, the algorithm operation ends and the result is output

Numerical Experiments
We first verify the accuracy of the above proposed MIP model by solving an instance with small-scale requests in the simplified Sioux waterfall network proposed by Tong et al. [3]. e experiments in the paper were performed on the Inter(R) Core(TM) i5-7300HQ CPU @ 2.50 GHz 8G RAM computer. e MIP model was solved by the commercial software GUROBI.

Small Numerical Experiments on a Simplified Sioux Waterfall Network.
To verify the proposed model, we test and compare it with the proposed algorithm by one example with 20 passenger requests on the simplified Sioux Falls network proposed by Tong et al. [4]. Some nodes meaning candidate stop locations (denoted by the blue and green dots), other nodes (denoted by the white dots), and transportation links are shown in Figure 4(a). e number on links means travel time. Two kinds of buses can be originally located at depot nodes, namely, nodes 1, 2, and 13 (denoted by the yellow dots), which have different load capacity, departure costs, and fuel costs. e maximum and minimum load capacity, Cap u k and Cap l k , of one with large capacity is 10 and 7, and the other with small capacity is 6 and 3. e departure cost c k d of is, respectively, 22 and 17, and the fuel cost per unit time c k f is 1 and 0.7. e ticket price is set to 10, and the penalty cost for each unserved passenger is half of the ticket price in this example. Note that, in this case, we assume that the end depot of the customized bus is virtual depots, of which distance to any point is zero. Table 4 lists related information of all passenger requests. Take passenger 1 as an example. Passenger 1 can be picked up at node 3 between time 11 and 15 and be delivered at node 7 between time 26 and 30. In addition, there are some passengers (passengers 4, 15, and 19) with more than one candidate pickup location in requests information. For passenger 4, he/she can select only one node as the final pickup location, namely, 4 or 5.
e commercial software GUROBI and ALNS ip are used to solve this instance. Finally, the route is generated as shown in Figure 4(b).

Optimal Solution for the Example.
Computational results of the example are found using, respectively, the commercial software GUROBI and algorithm ALNS ip to assess their respective performance. e computational results show that the total income is 62.6. Four passengers cannot be served, namely, passengers 5, 6, 9, and 16. Other passengers are assigned to two buses. e computational result is the optimal solution (Gap value is 0), as shown in Table 5. Although the optimal solution from the solver and ALNS ip are the same, shorter CPU time indicates the MIP solver has better performance in terms of the small-scale problems. Table 6 lists detailed vehicle routes and passenger assignments. Ten passengers are assigned to the bus with large capacity, and the rest six passengers to the bus with small one. eir routes are presented by red and orange lines in Figure 4(b).

Sensitive Analysis for Strategy 1. By the example on
Sioux Waterfall Network, we further analyse the impact of Strategy 1 on the result. e computational results are shown in Table 7.
As can be seen from Table 7, the heterogeneous fleet achieves better profit in both experiments. Furthermore, the (1) Calculate and remark arrival time, departure time, and waiting time of all locations on one route by Algorithm 3; (2) Begin with first location of the route and set u ⟵ 1; [u,n] MTS v ′ and obtain v, v is the location where MTS v ′ is the minimum value; (9) AT u � AT u + MTS u (10) Update arrival time, departure time, and waiting time of locations from the location u to n based on Algorithm 3;   result shows that the heterogeneous fleet serves more passengers than homogeneous fleets in the second experiment. It means more passengers can be served by a flexible solution from heterogeneous fleets. Homogeneous fleets cannot make flexible dispatch based on various passenger requests, which may increase cost or make more passengers unserviceable.

Benchmark Test with Unique Pickup or Delivery
Locations. Because there is no benchmark instance applied to customized bus problem, we choose DARP benchmark instances proposed by Cordeau and Laporte [5] as test instances to evaluate the performance of the ALNS ip proposed, which has similar properties as customized bus problems. Assuming all requests are serviced, the total distance is minimized as the objective function. Li et al. [35] tested, respectively, two solution evaluation approaches for ALNS (ALNS F is to allow only feasible solutions, ALNS I is to allow infeasible solutions) and found that although the results of ALNS I are better than ALNS F , ALNS I takes more time than ALNS F . Table 8 shows the calculation results of ALNS I , ALNS F , and ALNS ip on DARP instances. e first column in Table 8 represents the instance number, and the Gap(%) refers to the percentage deviation of the ALNS I , ALNS F , and ALNS ip compared to the best results from [5,45]. We test the ALNS ip on DARP benchmark instances and find it provides good results within 5000 iterations. e best result of ALNS ip is selected based on 10 runs. All results are shown in Table 8 and compared with the best results of ALNS I , ALNS F from Li et al. [35]. Table 8 shows that compared to ALNS F , the proposed ALNS ip can obtain better solutions in most instances where the results have on an average 84.16% decreased gap values with only 23.03% longer CPU time on average. In these instances, the ALNS ip generally gives better results but consumes longer time than ALNS F . But the exception is that for instances R2a, R10a, R8b, R9b, and R10b, ALNS ip searches better solutions but faster than ALNS F . Furthermore, ALNS ip generally gives similar (for instances, R1a, R2a, R4a, R7a, and R7b) or even worse results but consumes less time than ALNS I in these instances. Table 8 shows that the proposed ALNS ip searches solutions 60.71% higher gap values but 42.75% shorter CPU time on average than ALNS I . While the exception is that for instances R2a, R4a and R9a, ALNS ip obtains better solutions than ALNS I within a shorter time. As a result, it turned out that ALNS ip could get efficient and highquality results in a reasonable time. In addition, ALNS ip ensures the feasible solutions at each iteration. It is of great significance to solve the large-scale routing problems with multiple locations and time windows.

Benchmark Test with Multiple Pickup or Delivery
Locations. In the section, we improve DARP benchmark instances to further analyse the impact of Strategy 2 and test the performance of the ALNS ip in terms of multiple pickup or delivery locations problem. Some passengers are randomly appended pickup or delivery candidate locations, where new pickup or delivery locations are some locations closest to the passenger's original pickup or delivery location. Considering different scales, we choose R1a (smallscale problem), R2a (medium-scale problem), and R5a (large-scale problems) as tested instances to illustrate the following analysis.
(i) Analysis for Strategy 2: Table 9 compares the impacts of Strategy 2 or not. It can be seen that, compared with the strategy that does not consider multiple candidate locations, Strategy 2 provides more flexible and higher solution quality for the final solution. e results also show that compared with classic ALNS, ALNS ip can effectively solve the problem with multiple candidate locations in a reasonable time. (ii) Analysis for the number of candidate locations: to perform a sensitive analysis of the number of candidate locations, we test it from two perspectives, that is, the number of candidate locations of each passenger (2, 3, or 4 candidate locations) and the total number of people chosen to append candidate locations (5%, 10%, and 20% of passengers randomly). Overall, Strategy 2 with more than one location gives more applicable solutions than cases without multiple locations in terms of two perspectives. Table 10 shows that with the increase of candidate locations, it can obtain the better solution in small-scale problems. But for larger-scale instances, the results in Table 10 may be worse if one passenger has more candidate locations. Possibly, the reason is that multiple candidate locations cause more feasible solutions during the search, which expands the search area and need more iterations and time to search optimized solutions. erefore, given that passengers generally do not have too many candidate locations near their destination in practice, 2 or 3 candidate locations are more appropriate in practical issues.
No. Large

Number
Note that the number in Solution is customer's id.
(2) More candidate locations chosen by passengers in the platform make a breakthrough in the limitations of a single candidate location. From the passengers' perspective, it increases the probability for the passengers to be served. From the operators' perspective, it means more flexible route planning. Consequently, multiple candidate locations can be utilized to optimize pickup and delivery routes and serve as many passengers as possible to improve service quality.

Conclusion
For passengers living in congested metropolitan areas, the customized bus system offers flexible and high-accessibility bus service in the existing transportation network services. Compared with the general customized bus problem, this study constructs a new service form with heterogeneous fleets and multiple candidate locations to deal with practical issues. e two strategies are then analysed and evaluated their effects. To optimize the operator's profits based on the new service form proposed, a mixed-integer programming model is developed to describe the problem. It is not only suitable for solving small-scale problems but also simple to operate by GUROBI, a commercial solver. An improved ALNS ip is then designed to address a large-scale challenge. e initial ALNS algorithm has good performance in terms of routing problem. e time slack calculation method is employed to intelligently adjust vehicle timetables to optimize the wait and duration time. Two inserting operators are utilized to enhance the ability to search solution space, considering some passengers with multiple candidate locations.
Numerical experiments are constructed to evaluate the effectiveness of the proposed model and algorithm. e sensitive analysis on small-scale instances indicates that the mechanism with heterogeneous fleets and multiple candidate locations can serve more passengers and obtain more profits. A comparison among ALNS I and ALNS F indicates ALNS ip can be efficient in searching for better solutions within a reasonable time. In addition, ALNS ip also outperform other heuristic algorithms, like GA, PSO, and ABC, in search ability and avoiding premature. Finally, this study tests the customized bus service in real life to verify its theoretical research and practical implication.
is study aims to propose a new form of customized bus service based on heterogeneous fleets and multiple candidate locations. And the designed model and algorithm can offer robust solutions in different scales. In future research, our work can be conducted in the following three directions: (1) researching other models for passengers with random characteristics, who usually offer information in real time and hope to obtain dynamic service, (2) considering the cooperation mechanism between customized bus service and the other transportation systems, and (3) improving the solution algorithm through introducing distributed parallel computing to shorten computing time.
Data Availability e data of Sections 5.1 and 5.3 used to support the findings of this study are included within the article. e data of Section 5.2 are from [5].

Conflicts of Interest
e authors declare that they have no conflicts of interest.