Integrated Optimization of Route and Frequency for Rail Transit Feeder Buses under the Influence of Shared Motorcycles

: In this paper, we develop a multi-objective integrated optimization method for feeder buses of rail transit based on realistic considerations. We propose a bus stop selection method that considers the influence of shared motorcycles, which can score the importance of alternative bus stops and select those with the highest scores as objectives. The objective of the model in this paper is to minimize both the travel costs of passengers and the operating costs of the bus company. This is achieved by optimizing feeder bus routes, the frequency of departures, and interchange discounts to enhance the connectivity between feeder buses and rail transit. In addition, to ensure the feasibility of generated routes in the real road network, a genetic algorithm encoded with priority is used to solve this model. We use the Xingyao Road subway station in Kunming as an example, and the results show that the optimization method is effective.


Introduction
In some medium-sized cities in China, shared motorcycles have replaced shared bicycles as the dominant form of shared transportation, having reached a mature stage of development.Furthermore, China's bus system is also well-developed, thanks to continuous efforts to prioritize public transport.However, due to tremendous operational pressures, bus companies are now focusing more on feeder buses for rail transit systems, which offer lower operating costs and better returns.Medium-sized cities such as Xi'an, Kunming, and Changsha have built and maturely operated rail transit systems, although the operational mileage is relatively limited, and the coverage area is confined.Consequently, there is a reliance on feeder buses or shared motorcycles for last-mile connectivity.Nevertheless, due to the lack of regulation on shared motorcycles, feeder buses face significant challenges in competitiveness, resulting in less-than-optimal passenger flow.
The decline in passenger flow results in a decrease in revenue for bus companies.To maintain a balance between revenues and expenses, bus companies often reduce operating costs by reducing the frequency of departures [1].However, for certain groups, such as children, the elderly, passengers with large pieces of luggage, and people who cannot ride bicycles, feeder buses are vital.Therefore, it is necessary to optimize feeder buses to better meet their needs.
Bus route optimization is a complex problem that involves multiple objectives, such as transit network design and frequency setting problems (TNDFSPs) [2,3].However, common TNDFSPs mainly treat bus systems as separate systems and neglect the temporal discrepancies between the arrival times of the two systems.This oversight can lead to extended interchange durations.In fact, feeder buses are not as closely coordinated with rail transport as shared motorcycles.Consequently, it is imperative to strengthen the connection between the two networks.
network by developing a model to minimize the travel costs of passengers and operation costs of the company, and investigated the effect of the number of bus lines on the overall efficiency of the bus network using an enumeration method.Mishra et al. [11] optimized the headway and bus stop spacing for low-demand bus routes using a multi-objective evolutionary algorithm, NSGA-II, to minimize both operator and user costs.The solution shows that the optimal values of headway and bus stop spacing are underestimated if the optimization is done based on the assumptions that are typical for high-demand routes.Zhang et al. [12] established a multi-model electric bus scheduling model to collaboratively optimize vehicle travel plans and charging schemes by considering time-sharing tariffs and orderly charging strategies.The results can effectively reduce operating costs and charging costs and improve bus service levels.Several studies have shown that multiple short-distance bus lines have lower system costs than a single long-distance bus line.At the same time, pure electric vehicles, which have been widely commissioned in recent years, are also more suitable for operating on short-distance bus routes due to range limitations [13].Therefore, a common research trend among scholars and bus operators involves gradually abandoning long-distance bus lines and developing short-distance bus lines [14][15][16].However, there is a lack of research on feeder buses, which are typically short-distance bus lines.These feeder buses typically have a high turnover rate and rarely face issues related to insufficient capacity.Furthermore, the literature has given limited attention to the interchange discount between feeder buses and rail transit, which is a crucial factor for passengers and deserves further investigation.
Once the bus optimization model is built, it needs to be solved using various algorithms.The bus optimization problem is a complex vehicle routing problem (VRP), and machine learning algorithms and heuristic algorithms are generally used in the related literature.Suh and Jeong [17] restructured urban bus routes and planned their schedules using the Naïve Bayes classification method; the data can be cleaned according to the characteristics of each bus route.The developed route improvement prediction classifier can provide accurate route optimization advice.Noor et al. [18] used artificial neural network (ANN) and support vector machine (SVM) algorithms to predict the travel time, fuel consumption, and harmful emissions of different scenarios of university shuttle buses, and it was found that the ANN model predicted better than SVM.The robust performance of machine learning algorithms significantly enhances the accuracy of optimal solutions in optimization problems.However, these algorithms necessitate deep data mining and demand high-quality, voluminous datasets.Conversely, heuristic algorithms exhibit lower data requirements and can identify feasible solutions within acceptable margins, demonstrating considerable robustness.Numerous heuristic algorithms exist for optimization, each with distinct advantages and limitations.The current research trend involves integrating various algorithms to mitigate their individual shortcomings, thereby developing more efficient combinatorial algorithms.Ahern et al. [19] used a multi-objective simulated annealing algorithm to solve two variables (the bus route and departure frequency).The algorithm is divided into three phases-the first two phases are designed to solve the objectives of passengers and operators, and the last phase is a multi-objective search for a better compromise that satisfies both objectives in some regions, addressing the issue of the algorithm potentially falling into a local optimum.Zhong et al. [20] proposed an improved particle swarm algorithm to optimize the bus rapid transit route, which addresses the issue of the traditional particle swarm algorithm being prone to local optimums.The improvements include generating two populations in the initialization stage to enhance the optimization capability and avoid a single particle from falling into the local optimum.And combining the crossover operation in the genetic algorithm to exchange the positions of two selected particles to generate a new population.Aktaş et al. [21] proposed a variable neighborhood search (VNS) algorithm that finds high-quality scheduling solutions to optimize the performance of a single bus route during peak hours in a reasonable amount of time.The algorithm also allows for re-optimization of service types and departure times during operation based on real-time demand.Among the heuristic algorithms, the genetic algorithm is the most widely used [3,4,15,22].Compared with the above algorithms, such as the particle swarm algorithm and simulated annealing algorithm, the genetic algorithm has distinct advantages in diversity preservation and global search ability, and it is not easy to be trapped in the local optimum.Moreover, by selecting appropriate encoding methods based on variable types, utilizing different operators can significantly improve both the speed and effectiveness of genetic algorithm optimization.
A few other papers have analyzed the competitive relationship between different transportation modes and feeder buses.Yang [23] studied the influence of 16 variables on the choice of using a microcirculation bus by using a multinomial logistic model.And he constructed a bilayer optimization model, where the lower layer was a model for calculating the probability of choosing a microcirculation bus to travel, thus reflecting the influence of shared bicycles.Liu et al. [24] dynamically calculated the actual travel demands of feeder buses based on the users' choice behavior between shared bicycles and feeder buses, to understand the impact of shared bicycles on feeder bus operations, build a model to optimize the feeder bus route design and vehicle allocation, and solve it with a Lagrangian relaxation algorithm.Wei et al. [25] qualitatively and quantitatively analyzed the competition and cooperation between rail transit and surrounding bus lines based on geospatial considerations.A bus route optimization model was constructed based on the co-cooperation coefficient, which provides a new optimization method for bus route adjustments under the influence of the metro.In the literature, most studies on the competitive relationship between different modes of transportation begin by analyzing choice probabilities.However, fewer studies propose specific and effective methods to optimize feeder bus services to further enhance their advantages and competitiveness.
The remaining sections of this paper are organized as follows: Section 2 proposes a method of scoring the importance of bus stops, considering the influence of shared motorcycles to provide a reference for the stop selection of feeder buses.It also presents a multi-objective integrated optimization model with the objective of achieving the highest total score, the lowest travel cost, and operation costs.The model is also solved using a genetic algorithm with priority coding.In Section 3, the route, frequency of departures, and interchange discount of the feeder bus are optimized by using the Xingyao Road station in Kunming as an example.Section 4 presents a summary of the method proposed in this paper.

Problem Statement
Feeder buses are types of regular buses.They have at least one stop near rail transit stations, linking rail transit with residential areas, office buildings, hospitals, and other facilities to meet passengers' transfer needs.Additionally, they possess all the characteristics of regular buses, such as operating on fixed routes and stops, having a flat fare, and adhering to a specified timetable.
There are several points of interest (POIs) in the study area, including residential communities, schools, commercial areas, and office buildings.Each POI has one gate or more gates for resident access, collectively denoted as set R. During a trip, a commuter starts at a particular gate and walks to the nearest parking point or bus stop, thus creating some correspondence.The set of parking points is denoted as P, forming a one-to-one relationship with the gate set R. Likewise, the set of bus stops is represented by S, and it forms a one-to-many relationship with the gate set R. This is attributed to the fact that several gates may share the same nearest bus stop.Knowing the correspondence between these three sets allows for the convenient calculation of ridership at each bus stop.Figure 1a illustrates the relationship between the three sets in the abstract, while Figure 1b provides a more intuitive demonstration of this correspondence based on actual neighborhoods.
optimization solution aims to minimize the travel costs and operating costs.In the bus stops mentioned above are treated as standard bus stops providing basi to passengers, excluding large bus stops with complex functions and situatio there are multiple platforms at a single stop.To facilitate the study, we made certain assumptions about this problem.consider that feeder buses operate over a short mileage, experience fewer delay sections, and are predictable to passengers.Thus, all buses are assumed to tra average speed [4,26], without considering any possible delays caused by interse nals or traffic congestion [27].Second, it is assumed that passengers will only trav stops with the shortest walking distance to wait for buses [26].Third, the arriva passengers at the bus stop is assumed to follow a random distribution in th Fourth, we only consider the interchange stop as the destination of passenger Thus, the OD (origin-destination) pairs between the remaining stops are not co Fifth, we assume that passengers arrive at the rail transit station and take the nea without waiting for an extra trip.Finally, this study was conducted during peak ing hours, as the passenger flow is high and concentrated during these times.By no adjustments were made to the departure frequencies and transfer discounts d peak periods.This paper defines an undirected graph, where S denotes the set of vertices, i.e., alternative stops for feeder buses, and X denotes the set of edges, i.e., alternative routes for feeder buses.Based on the selected stops that are more suitable for feeder buses, the final optimization solution aims to minimize the travel costs and operating costs.In addition, the bus stops mentioned above are treated as standard bus stops providing basic services to passengers, excluding large bus stops with complex functions and situations where there are multiple platforms at a single stop.

Notation Explanation
To facilitate the study, we made certain assumptions about this problem.First, we consider that feeder buses operate over a short mileage, experience fewer delays at intersections, and are predictable to passengers.Thus, all buses are assumed to travel at an average speed [4,26], without considering any possible delays caused by intersection signals or traffic congestion [27].Second, it is assumed that passengers will only travel to bus stops with the shortest walking distance to wait for buses [26].Third, the arrival time of passengers at the bus stop is assumed to follow a random distribution in this model.Fourth, we only consider the interchange stop as the destination of passengers [26][27][28].Thus, the OD (origin-destination) pairs between the remaining stops are not considered.Fifth, we assume that passengers arrive at the rail transit station and take the nearest train without waiting for an extra trip.Finally, this study was conducted during peak commuting hours, as the passenger flow is high and concentrated during these times.By default, no adjustments were made to the departure frequencies and transfer discounts during off-peak periods.The walking distance from gate r to the nearest bus stop meters (m)

D ijl
The travel distance of feeder bus line l on the shortest path between stops i and j meters (m) The average walking speed of commuters kilometers per hour (km/h) The average travel speed of feeder buses kilometers per hour (km/h) The number of stops that feeder bus line l has passed when it travels to stop i pieces (pcs) The number of stops for the feeder bus line l pieces (pcs) t s The average stopping time of feeder buses seconds (s) q l The rated capacity per bus on the feeder bus line l person T M n The arrival time of the nth train of rail transit -T f kl The arrival time of the kth bus of feeder bus line l - The average transfer time for commuters from the feeder bus stop to the metro stop seconds (s) The fare of one ride on feeder buses CNY C u The average travel time cost for commuters CNY C o The average operating cost per kilometer for bus companies CNY Q il The ridership of feeder bus line l at stop i person D bw r The walking distance from gate r to the nearest parking point meters (m) The riding distance from parking point x to the metro stop meters (m) The average riding speed of commuters kilometers per hour (km/h) T b The average return time for a shared motorcycle seconds (s) The fare of one ride on a shared motorcycle CNY S i The combined score of stop i -

K min
The minimum allowed frequency of departure of feeder buses times per hour The maximum allowed frequency of departure of feeder buses times per hour The minimum allowed length of the feeder bus route kilometers (km) The maximum allowed length of the feeder bus route kilometers (km) n min The minimum allowed number of stops for feeder buses pieces (pcs) The maximum allowed number of stops for feeder buses pieces (pcs)

Variables Description Units
X ijl 1 if feeder bus line l passes through stops i and j, which are adjacent to each other, and 0 otherwise - The frequency of departures of feeder bus line l times per hour

Disc
The discount rate for the feeder bus and rail transit interchange -

Model Development
The model developed in this paper is a multi-objective planning model.Firstly, the model considers the probability of choosing shared motorcycles in the study area, as well as the distinctive characteristics of bus stops.Subsequently, the model identifies feeder bus stops that can generate greater benefits for the overall system.The stops are selected by comprehensive scoring with multiple indicators, which will assist in optimizing the routes in a more purposeful way.Section 2.4 will provide a detailed description of the scoring process.The importance score of the feeder bus stops can be expressed in Equation (1), where Prior to determining the probability of selecting shared motorcycles, it is essential to first compute the travel costs associated with commuters who opt for shared motorcycles and those who choose feeder buses separately.The travel costs for feeder buses are expressed in Equation (2), and the travel costs for shared motorcycles are expressed in Equation (3). where V w denotes the average walking time from gate r to the nearest bus stop; V f denotes the travel time of feeder bus line l on the shortest path between stops i and j; n il t s is the total stopping time when feeder bus line l travels to stop i; 1 2K l is the average waiting time of passengers for feeder bus l (passenger arrivals satisfy a random distribution); T M n − T f kl − t is the waiting time for the feeder bus and rail transit interchange; D bw r V w denotes the average walking time from gate r to the corresponding parking point; D x V b denotes the average riding time from parking point x to the interchange parking point; the waiting time for the rail transit of cycling passengers can be expressed as half of the interval between two adjacent trains, i.e., ; this paper defines a return time for shared motorcycles due to insufficient capacity of the parking point, i.e., T b .During the morning rush hour, the task of returning shared motorcycles to designated parking points becomes increasingly challenging.This difficulty arises primarily due to the substantial number of riders who transfer to rail transit during this period.Consequently, the time taken to return the motorcycles is longer than during off-peak hours [29].For this paper, the scenarios of returning motorcycles can be classified into three categories:    β 1 , enough spaces left, no problem to find a parking space (scene1) β 2 , few spaces left, need time to find a parking space (scene2) β 3 , no spaces left, need to find another parking space (scene3) where β 1,2,3 indicates the proportion of the above three scenarios in the total duration of the morning rush hours, respectively, and the remaining spaces ≤ 5 represent scene 2. Correspondingly, the return time is different for each type of scenario.When there are enough spaces left, the time for commuters to return to the shared motorcycles is T b 1 ; when there are few spaces left, the time is T b 2 ; when there are no spaces left, assume the riders choose the nearest parking point to return the motorcycles, and the time is , where D b denotes the distance to the nearest parking point.Thus, the average return time for shared motorcycles is calculated by Equation (4).
In addition, the model should consider the interests of both travelers and operators by balancing the travel costs for passengers and the operating costs for the bus company.This approach aims to determine an optimal, compromised solution.The net operating costs of the bus company can be expressed as follows: The first half of this formula describes the direct operating costs of bus companies, while the second half describes the ticket revenue of bus companies.Since public transportation in China is typically considered a public good, the operating cost exceeds the profit.Thus, the above formula can also be used to describe the amount of loss incurred by bus companies.
To facilitate the optimization, this paper transforms the multi-objective planning problem into a single-objective optimization model using a weighted summation approach [30].The model can be fully represented as follows: subject to the following: Coefficients α 1 , α 2 in objective (6) denote the weights of passenger travel costs and operation costs, respectively.Constraint ( 7) is a capacity constraint, where q l is the capacity of feeder bus line l; constraint (8) is the frequency constraint, K min and K max indicate the minimum allowed frequency of departure and the maximum frequency, both obtained from an actual experience; constraint (9) ensures that the interchange discount rate cannot equal 1, and the lowest can be 0, which implies that the interchange is free; constraint (10) ensures that the feeder bus must depart from the interchange stop for rail transit; constraint (11) ensures that path (i, j) is served by, at most, one bus route [4]; constraint (12) ensures that each stop has the same number of incoming and outgoing bus routes; constraint ( 13) is a typical Miller-Tucker-Zemlin (MTZ) constraint, which aims to eliminate possible subloops in the TSP problem; constraint (14) is the path length constraint, E min and E max denote the minimum and maximum lengths of the feeder bus route, respectively, and both are obtained from actual experiences and related to the service area; constraint ( 15) is a constraint on the number of feeder bus stops, n min and n max denote the minimum and maximum stop numbers for feeder buses, respectively, and both are obtained from actual experiences and related to the service area; constraint ( 16) is a decision variable constraint when path (i, j) is selected, X ijl = 1, otherwise, 0.
The above constructs an optimization model utilizing three decision variables: route, departure frequency, and interchange discount, to address the objective of minimizing the total system cost.Consequently, this paper hypothesizes that adjusting the departure frequency and providing transfer discounts can have a positive optimization effect on the objective function, and effectively reduce the total system cost.In this paper, our objective is to provide better guidance for determining the stops at which feeder buses should stop.To achieve this, we consider the trip characteristics of each stop and propose a comprehensive rating method for the alternative stops.The model evaluates bus stops using three indicators: ridership (A 1 ), distance to the interchange stop (A 2 ), and the probability of choosing buses (A 3 ).Stops with higher scores on these indicators are more likely to be recommended for feeder buses.Ridership refers to the total number of passengers who transfer to rail transit using feeder buses and shared motorcycles, which can best reflect the importance of stops.The probability of choosing buses is determined by the travel costs of both modes of transport.This indicator can best reflect the influence of shared motorcycles on passenger preferences for buses.For instance, if the cost of a bus trip is significantly higher than that of a motorcycle trip at a certain stop, the probability of choosing feeder buses would be low, resulting in a low score.Consequently, such a stop may not be chosen for feeder buses.Additionally, since the average travel distance of buses is longer than that of motorcycles, stops that are located farther away from the interchange stop are more advantageous for feeder buses.Concerning this indicator, a higher score is assigned to stops that are situated at greater distances.
This paper examines the travel costs of feeder buses and shared motorcycles at each stop within the study area individually.Furthermore, it explores the preferences of passengers for these two travel modes.The probability of selecting each travel mode is determined by Equations ( 17) and ( 18) following the prediction generated by binomial logistic regression.
Both ridership and distance to the interchange stop can be obtained through field surveys.These three indicators have different degrees of importance in selecting stops for feeder buses and, thus, need to be assigned weights.In this paper, the study of indicators needs to consider both the preferences of practitioners and objectivity.Thus, we choose a combination of the analytic hierarchy process (AHP) and the entropy method to calculate the weights.After calculating the subjective and objective weights separately, we combine the weights of the indicators as W i .
φ i -The weights obtained from the AHP.ω i -The weights obtained from the entropy method.
Once the indicators are assigned, they can be calculated by multiplying the weights with the data for each indicator and then normalizing the values to obtain the comprehensive scores of stops.

Genetic Algorithm
The decision variables in the above model include both discrete 0-1 variables, which are used to generate routes, and continuous variables, which are used to determine the optimal frequency of departures and interchange discounts.As a result, the present model can be classified as a mixed-integer nonlinear model, which is more complex than the standard VRP and is also recognized as NP-hard.Solving NP-hard problems is challenging, and the heuristic algorithm is appropriate for addressing large-scale VRP.It can obtain good feasible solutions and has a faster convergence speed.Therefore, this paper chooses the genetic algorithm as a heuristic approach to solve the proposed model.
In some literature addressing VRP using genetic algorithms, the paths (i.e., feasible solutions) are encoded as real numbers based on the node numbers [4].However, in real road networks, some nodes may not have connections or the route between two nodes may need to pass through additional nodes.As a result, the decoding of the new chromosomes may not yield feasible routes.To address this issue, this paper proposes an indirect coding method known as the priority order [28].
In this encoding scheme, the chromosome represents a priority order rather than a specific path.The priority decreases from left to right in the chromosome.If the starting point of the path is specified, it is the first node in the path.The next node in the path is selected as the element with the highest priority among the set of points that can be directly connected to the previous node (i.e., no other nodes exist between them).This process continues until the length of the path reaches the maximum, or there are no more nodes that can be directly connected.If the starting point is not specified, the first node in the path is the first element according to the priority order.Thus, the final route must be available in the real road network.Before applying the priority order for encoding, it is important to obtain the connectivity graph of the road network.This graph represents the network's structure, where each arc indicates the existence of a path between two nodes without the need to go through another node.Figure 2 depicts the chromosome-to-path decoding process using the genetic algorithm described in this paper.The flow chart of solving the model is shown in Figure 3: In this paper, the size of the study area will determine whether one or more feeder bus routes should be planned.Once a chromosome satisfies the termination condition for decoding, the remaining unselected stops are treated as a new set, and another route is decoded starting from the first stop.This process can be repeated to obtain multiple routes.
In addition to generating routes, this paper also utilizes the genetic algorithm to determine the optimal frequency of departures and interchange discounts.These two parameters are continuous variables, and they are represented by the last two digits of the chromosome using real number coding.By optimizing these parameters along with the routes, a more convincing solution can be obtained.
The flow chart of solving the model is shown in Figure 3: The flow chart of solving the model is shown in Figure 3:

Basic Description
The preliminary investigation revealed that the Xingyao Road Metro Station has a high commuter flow, and there are both stable operating feeder buses and a substantial number of shared motorcycles near the entrances, highlighting a clear demand for feeder services.However, the feeder bus suffers from poor passenger flow, significantly impacted by the shared motorcycles.Thus, this paper selects the area around the Xingyao Road Station as the study area.
The specific location and scope of the study area are illustrated in Figure 4.The area is situated on the B exit side of the Xingyao Road station and encompasses 28 bus stops, 126 parking points for shared motorcycles, and one feeder bus line in regular operation.Figure 5 displays the route of the feeder bus K58 and the distribution of the 28 bus stops within the area.

Basic Description
The preliminary investigation revealed that the Xingyao Road Metro Station has a high commuter flow, and there are both stable operating feeder buses and a substantial number of shared motorcycles near the entrances, highlighting a clear demand for feeder services.However, the feeder bus suffers from poor passenger flow, significantly impacted by the shared motorcycles.Thus, this paper selects the area around the Xingyao Road Station as the study area.
The specific location and scope of the study area are illustrated in Figure 4.The area is situated on the B exit side of the Xingyao Road station and encompasses 28 bus stops, 126 parking points for shared motorcycles, and one feeder bus line in regular operation.Figure 5 displays the route of the feeder bus K58 and the distribution of the 28 bus stops within the area.Based on the methodology presented in Section 2.4.2, we need to convert the real road network in the study area into a connectivity graph, with bus stops serving as network nodes.The resulting connectivity diagram is illustrated in Figure 6.The numbers in the network nodes represent the codes of the different bus stops.Based on the methodology presented in Section 2.4.2, we need to convert the real road network in the study area into a connectivity graph, with bus stops serving as network nodes.The resulting connectivity diagram is illustrated in Figure 6.The numbers in the network nodes represent the codes of the different bus stops.There are over 30 points of interest in the study area, with a total of 106 gates.A preliminary survey found that some gates had low pedestrian activity.To reduce unnecessary workload, some of these gates were excluded from the study, resulting in a simplified research object of 57 gates.Further analysis revealed that bus stops 2, 3, 12, and 21 did not have corresponding gates with the shortest walking distances.Under the assumption that the public only walks to the nearest bus stop, these stops are deemed insignificant for the study and will be excluded in subsequent analyses.

Indicator Values
Since the morning peak on weekdays represents the week with the highest concentration of travel demand, this study focuses on investigating the study area during the morning peak hours (7:00 a.m.-9:00 a.m.) [29] from 13 March to 17 March 2023.
Xingyao Road Station is one station along Kunming Metro Line 1.During the morning peak, trains depart from this station in 5-min intervals.Zijuncun (Caiyun North Road) is a bus stop that is situated near Exit B of Xingyao Road Station.This stop serves as both There are over 30 points of interest in the study area, with a total of 106 gates.A preliminary survey found that some gates had low pedestrian activity.To reduce unnecessary workload, some of these gates were excluded from the study, resulting in a simplified research object of 57 gates.Further analysis revealed that bus stops 2, 3, 12, and 21 did not have corresponding gates with the shortest walking distances.Under the assumption that the public only walks to the nearest bus stop, these stops are deemed insignificant for the study and will be excluded in subsequent analyses.

Indicator Values
Since the morning peak on weekdays represents the week with the highest concentration of travel demand, this study focuses on investigating the study area during the morning peak hours (7:00 a.m.-9:00 a.m.) [29] from 13 March to 17 March 2023.
Xingyao Road Station is one station along Kunming Metro Line 1.During the morning peak, trains depart from this station in 5-min intervals.Zijuncun (Caiyun North Road) is a bus stop that is situated near Exit B of Xingyao Road Station.This stop serves as both the starting and ending point for bus line K58.The departure interval for K58 during the morning peak is usually between 8 and 18 min, and the total travel time approximately 15 min.
Under the assumptions in Section 2, the values of some indicators and variables in this paper are presented in Table 2. Since public transportation in China is characterized by a strong public welfare nature, the objective function of this paper prioritizes reducing passengers' travel costs, so α 1 = 0.8 and α 2 = 0.2.Regarding the decision variables, since the departure interval of the subway is 5 min, this paper matches the departure time of the feeder bus with the arrival time of the subway in order to enhance the connection between the feeder bus and the subway.Thus, the departure interval of the feeder bus is a multiple of 5, and K l denotes the values of 3, 4, 6, and 12 per hour.Disc denotes the values of 0, 0.2, 0.4, 0.6, and 0.8.The bus fare in Kunming is a unified fare system, which is not related to the riding distance, so C ft = CNY 1.8/person.

C bt
The fare of shared motorcycles in this paper refers to the price of the weekly ticket of Hello Bike, so C bt = CNY 1.2/ride.The Values of Indicators in Constraints The distance traveled from each bus stop to the Zijuncun (Caiyun North Road) stop is readily available on the map as an indicator in the stop scoring.This paper investigates passenger flow in two parts.Firstly, the passenger flow of the current stops of the feeder bus is examined, which represents the direct ridership of the stop and is obtained through a follow-up survey.Secondly, the passenger flow of motorcycles is considered, which represents the potential passenger flow of the stop, and some residents may shift their travel mode when the cost of the feeder bus decreases.We initially utilized real-time data from the applications of three shared motorcycle operators: Hello Bike, Meituan, and Qingju, to quantify the total number of individuals who ride shared motorcycles to Xingyao Road Station.By comparing the vehicle counts at each stop before and after the peak period, we are able to ascertain the number of trips made at each stop during peak hours.Finally, the total number of people is then allocated based on the proportion of trips to determine the number of people who transfer to rail transit via motorcycles at each bus stop.
According to the models for the travel costs of feeder buses and shared motorcycles, established in Section 2, as well as the index values mentioned earlier, the travel costs for each travel mode at each bus stop can be calculated.This paper then uses the costs of the two travel modes as independent variables and the probability of passengers choosing feeder buses as the dependent variable in a binary logistic regression analysis.In this regression, the dependent variable y can only take two values, 0 and 1.A value of y = 1 indicates that passengers choose feeder buses, while a value of y = 0 indicates that passengers choose shared motorcycles.The regression model obtained is shown below, with x 1 representing the travel cost of feeder buses and x 2 representing the travel cost of shared motorcycles.By substituting the corresponding values of x 1 and x 2 for each bus stop in the regression model, corresponding probabilities of choosing buses are obtained.ŷ = −0.779x 1 + 0.903x 2 + 1.682 (20) After completing the survey of bus stops, including three indicators of ridership, namely, (A 1 ), distance to interchange stop (A 2 ), and the probability of choosing buses (A 3 ), the information on alternative stops is obtained as shown in Table 3.In this paper, we utilize a hybrid methodology consisting of AHP and the entropy method to assign weights to the indicators discussed above.To ensure the scientific and rational nature of the weights derived from the AHP, a group of ten experts specializing in transportation planning have rated the indicators presented in this paper.Table 4 displays the output of the calculations, which were performed using the AHP, the x the combined subjective and objective assignment approach.After determining the weights for each indicator, we compute the overall score for each stop.By normalizing the computed results, we obtain the final scores for each stop, as shown in Table 5.

Results
According to the flowchart presented in Section 2, the parameters for the genetic algorithm used to solve the model are set as follows: The population size is set to 100 individuals, the maximum number of iterations is set to 200, the generation gap is set to 0.9, the crossover probability is set to 0.9, and the mutation probability is 0.1.
The genetic algorithm is iterated in the order shown in Figure 3. Figure 7 shows the change of the objective function value.The results show that stability is achieved after approximately 60 iterations.the rail transit interchange stop.Figure 8 provides a visual comparison of the route before and after the optimization.From the comparison, it is evident that the optimized route covers a larger area that is further away from the interchange stop and has a higher total travel volume, albeit with a slight detour.In some of the areas originally covered, it was more convenient for commuters to choose shared motorcycles, and the advantage of feeder buses was not as great, so they were optimized and no longer covered.Depending on the values of the two decision variables, l K and Disc , there are 20 different combination options for each route.The morning peak time studied in this paper is 2 h, so the value of l K is doubled.Table 6 shows the objective values of the 20 combinations under the optimal route, where the optimal results can be obtained when l K = 12 and Disc = 0.6.In addition, Figure 9a illustrates the trend of changes in objective val- ues resulting from different discount rates, while keeping the frequency of departures constant.It can be found that when l K = 24, the higher the interchange discount, the higher the overall cost.However, there reaches a point where the reduction in passenger costs can no longer compensate for the operational cost of the enterprise, and the enterprise faces difficulties in operation.;when l K = 6, the lower the interchange discount, the higher the overall cost, at which time the passenger cost is high; when l K = 6 or 8, the overall cost shows a trend of first falling and then rising, and the interchange discount does not consistently favor any specific direction.Figure 9b shows the impact of different After optimization, there is only one linear feeder bus line in the study area, with 12 stops.The route follows the sequence of 1-10-13-14-15-17-19-22-23-26-25-24, starting from the rail transit interchange stop.Figure 8 provides a visual comparison of the route before and after the optimization.From the comparison, it is evident that the optimized route covers a larger area that is further away from the interchange stop and has a higher total travel volume, albeit with a slight detour.In some of the areas originally covered, it was more convenient for commuters to choose shared motorcycles, and the advantage of feeder buses was not as great, so they were optimized and no longer covered.Systems 2024, 12, x FOR PEER REVIEW the rail transit interchange stop.Figure 8 provides a visual comparison of the rou and after the optimization.From the comparison, it is evident that the optimiz covers a larger area that is further away from the interchange stop and has a hig travel volume, albeit with a slight detour.In some of the areas originally covere more convenient for commuters to choose shared motorcycles, and the adva feeder buses was not as great, so they were optimized and no longer covered.Depending on the values of the two decision variables, l K and Disc , ther different combination options for each route.The morning peak time studied in th is 2 h, so the value of l K is doubled.Table 6 shows the objective values of the 20 nations under the optimal route, where the optimal results can be obtained wh Depending on the values of the two decision variables, K l and Disc, there are 20 different combination options for each route.The morning peak time studied in this paper is 2 h, so the value of K l is doubled.Table 6 shows the objective values of the 20 combinations under the optimal route, where the optimal results can be obtained when K l = 12 and Disc = 0.6.In addition, Figure 9a illustrates the trend of changes in objective values resulting from different discount rates, while keeping the frequency of departures constant.It can be found that when K l = 24, the higher the interchange discount, the higher the overall cost.However, there reaches a point where the reduction in passenger costs can no longer compensate for the operational cost of the enterprise, and the enterprise faces difficulties in operation.;when K l = 6, the lower the interchange discount, the higher the overall cost, at which time the passenger cost is high; when K l = 6 or 8, the overall cost shows a trend of first falling and then rising, and the interchange discount does not consistently favor any specific direction.Figure 9b shows the impact of different frequencies on the change in objective value while keeping the interchange discount constant.It can be found that when the frequency is reduced, the total cost always shows a trend of decreasing first and then increasing, and the lowest cost point occurs at = 12 or 8, which further indicates that too high or too low of a frequency is not a better choice for the whole.The whole running time of the optimized feeder bus is about 18 min, so the matching departure and arrival moments of the feeder bus are calculated based on the arrival moments of the rail transit as shown in Table 7 [27].The optimized feeder bus solution has led to an increase in travel costs at certain stops, due to the elimination of some stops and the increase in detour times.However, overall, the solution has proven to be beneficial.The introduction of more accurate departure times and departure intervals has helped to minimize commute times for most residents in the study area.Additionally, the incorporation of interchange discounts and the adjustment of stops have further reduced travel costs.As a result, the probability of com-  The optimized feeder bus solution has led to an increase in travel costs at certain stops, due to the elimination of some stops and the increase in detour times.However, overall, the solution has proven to be beneficial.The introduction of more accurate departure times and departure intervals has helped to minimize commute times for most residents in the study area.Additionally, the incorporation of interchange discounts and the adjustment of stops have further reduced travel costs.As a result, the probability of commuters choosing feeder buses has increased significantly.In terms of financial impact, the total cost of travel for all gates has decreased from CNY 687.94 to CNY 591.50, representing a notable decrease of 14.02%.Furthermore, the number of predicted bus trips has increased from 132 to 243, marking a substantial increase of 84.09%.
In order to further verify the positive effect of the interchange discount on integrated optimization, this study conducted an analysis where the decision variable of the interchange discount was removed.Figure 10 depicts the results of this analysis, presenting the iterative values of the objective function.The optimized route for the feeder bus without the interchange discount is as follows: 1-10-13-14-15-17-19-23-24-25-27-28, with a total of 11 stops, and the specific route is illustrated in Figure 11.Upon examining the objective value, it can be observed that the optimal result after removing the interchange discount is 40.3724.This value represents a 9.26% increase compared to the result obtained through integrated optimization.Moreover, the total travel cost amounts to CNY 608.80, indicating a 2.92% increase.Additionally, the number of predicted bus trips reduces to 216, marking a decrease of 11.11%.Although including the interchange discount may elevate the operating costs for the bus company, incorporating this variable proves to be beneficial for the feeder bus system when considering the overall perspective.

Conclusions
This paper proposes an integrated optimization model for determining the optimal routes and frequencies of feeder buses under the influence of shared motorcycles.The model assigns a score to each bus stop and aims to select stops with the highest total score.The objective is to minimize travel costs and operating costs of feeder buses, which reflects the game theory relationship between social and enterprise interests.The model incorporates three decision variables, allowing for simultaneous optimization of the bus route, the frequency of departures, and interchange discounts.
The consideration of the impact of shared motorcycles in this paper is mainly reflected in the scoring of bus stops.Shared motorcycles face drawbacks such as extended return times, slow speeds, and limited travel distances during high-demand periods.Considering these characteristics, this study constructs a model to calculate the travel cost of shared motorcycles and predict the selection probability of both shared motorcycles and feeder buses.This is done using a binomial logistic regression model.By scoring alternative bus stops based on three indicators-ridership, distance to the interchange stop, and the probability of choosing buses-the most suitable stops for feeder buses can be determined.
Finally, this paper addresses the above model by employing a genetic algorithm that utilizes priority order coding.This approach guarantees that the routes generated by the algorithm are feasible in the actual road network.Furthermore, it can generate different numbers of feeder bus routes according to the scale of the study area.The findings of this research are summarized as follows:

Conclusions
This paper proposes an integrated optimization model for determining the optimal routes and frequencies of feeder buses under the influence of shared motorcycles.The model assigns a score to each bus stop and aims to select stops with the highest total score.The objective is to minimize travel costs and operating costs of feeder buses, which reflects the game theory relationship between social and enterprise interests.The model incorporates three decision variables, allowing for simultaneous optimization of the bus route, the frequency of departures, and interchange discounts.
The consideration of the impact of shared motorcycles in this paper is mainly reflected in the scoring of bus stops.Shared motorcycles face drawbacks such as extended return times, slow speeds, and limited travel distances during high-demand periods.Considering these characteristics, this study constructs a model to calculate the travel cost of shared motorcycles and predict the selection probability of both shared motorcycles and feeder buses.This is done using a binomial logistic regression model.By scoring alternative bus stops based on three indicators-ridership, distance to the interchange stop, and the probability of choosing buses-the most suitable stops for feeder buses can be determined.
Finally, this paper addresses the above model by employing a genetic algorithm that utilizes priority order coding.This approach guarantees that the routes generated by the algorithm are feasible in the actual road network.Furthermore, it can generate different numbers of feeder bus routes according to the scale of the study area.The findings of this research are summarized as follows: 1.
A well-coordinated connection between feeder buses and rail transit, along with a suitable interchange discount, has been proven to considerably reduce the expenses of bus travel and encourage residents to opt for buses as their preferred mode of transportation.In the context of this study, the cost of bus travel is projected to reduce by 14.02%, and it is anticipated that the number of individuals opting for bus travel will increase by 84.09%.This advantage is currently difficult to achieve for shared motorcycles, making the smooth transfer between feeder buses and rail transit a critical factor in providing feeder buses with a competitive advantage over shared motorcycles.2.
Increasing the frequency of departures and offering an interchange discount may lead to higher operating costs for the bus company.The benefits of these optimization measures include attracting more residents to choose feeder buses and increasing the operating revenue, thereby enabling bus companies to break free from the vicious cycle of declining ridership.

3.
The method proposed in this paper for comprehensively scoring stops based on three indicators can lead to better driving routes and stops for feeder buses under the influence of shared motorcycles.On the one hand, it reduces the coverage of feeder buses in areas closer to rail transit stations and reduces competition at stops where shared motorcycles have obvious advantages; on the other hand, it strengthens the coverage of areas with higher travel volume and further from rail transit stations, which can attract hidden passenger traffic.
In future work, we will consider refining the interchange between buses and other modes of transport, exploring the key factors that may affect the level of bus service or passengers' travel choices in the interchange.Additionally, we will evaluate the current state of bus transfer services.Beyond the traffic flow and location characteristics of bus stations, numerous other aspects warrant a thorough investigation, such as the type of bus station and its facilities, major hubs with transfer and distribution functions, and the comfort level of bus stop facilities, all of which may influence passenger preferences for buses.This paper optimizes and adjusts the routes and stops of feeder buses in the case study, but this may inadvertently reduce the overall accessibility of the study area.Therefore, in subsequent will incorporate accessibility as an objective in our models to ensure regional bus accessibility while enhancing bus efficiency.Finally, by segmenting the speeds of different stages in the bus operation process and accounting for varying degrees of delays, we aim to improve the study's relevance and validate the fault tolerance of the optimization results.

21 Figure 2 .
Figure 2. The process of decoding a chromosome into a route.

Figure 2 .
Figure 2. The process of decoding a chromosome into a route.

Figure 2 .
Figure 2. The process of decoding a chromosome into a route.

Figure 3 .
Figure 3. Flowchart for solving the model in this paper.Figure 3. Flowchart for solving the model in this paper.

Figure 3 .
Figure 3. Flowchart for solving the model in this paper.Figure 3. Flowchart for solving the model in this paper.
Systems 2024, 12, x FOR PEER REVIEW 12 of 21

Figure 4 .
Figure 4.The specific location of the study area.Figure 4. The specific location of the study area.

Figure 4 .
Figure 4.The specific location of the study area.Figure 4. The specific location of the study area.

Figure 4 .
Figure 4.The specific location of the study area.

Figure 5 .
Figure 5. Current operating bus route and stop layouts.

Figure 5 .
Figure 5. Current operating bus route and stop layouts.

21 Figure 6 .
Figure 6.Bus stop connectivity diagram in the study area.

Figure 6 .
Figure 6.Bus stop connectivity diagram in the study area.

Figure 8 .
Figure 8. Optimization results for routes and stops.

Figure 8 .
Figure 8. Optimization results for routes and stops.

Figure 8 .
Figure 8. Optimization results for routes and stops.

Figure 9 .
Figure 9. Trends in objective values due to a certain variable when the other variable is fixed.(a) When the frequency of departures is fixed.(b) When the interchange discount is fixed.

Figure 9 .
Figure 9. Trends in objective values due to a certain variable when the other variable is fixed.(a) When the frequency of departures is fixed.(b) When the interchange discount is fixed.

Systems 2024 , 21 Figure 10 .
Figure 10.Iterative process for optimal values without the interchange discount.Figure 10.Iterative process for optimal values without the interchange discount.

Figure 10 .
Figure 10.Iterative process for optimal values without the interchange discount.Figure 10.Iterative process for optimal values without the interchange discount.

Figure 10 .
Figure 10.Iterative process for optimal values without the interchange discount.

Figure 11 .
Figure 11.Optimization results for routes and stops without interchange discounts.

Figure 11 .
Figure 11.Optimization results for routes and stops without interchange discounts.

Table 1
introduces the notations involved in the problem.

Table 1 .
General subscripts, sets, input parameters, and decision variables.

Table 1
introduces the notations involved in the problem.

Table 1 .
General subscripts, sets, input parameters, and decision variables.

Table 2 .
The values of some indicators and variables in the models.

Values of Indicators and Variables in the Objective Function
[31] average hourly wage of employees in Kunming is CNY 44/h, so C u = CNY 26/h[31].C ft

Table 3 .
Information on the three indicators for alternative stops.

Table 4 .
Weights of the three indicators.

Table 5 .
Final scores for each bus stop.

Table 6 .
Objective values for different frequency and interchange discount combination schemes under the optimal route.

Table 7 .
Optimized arrival timetable for feeder buses.