Addressing the Bike Repositioning Problem in Bike Sharing System: A Two-Stage Stochastic Programming Model

In this paper, a bike repositioning problem with stochastic demand is studied.*e problem is formulated as a two-stage stochastic programming model to optimize the routing and loading/unloading decisions of the repositioning truck at each station and depot under stochastic demands. *e goal of the model is to minimize the expected total sum of the transportation costs, the expected penalty costs at all stations, and the holding cost of the depot. A simulated annealing algorithm is developed to solve the model. Numerical experiments are conducted on a set of instances from 20 to 90 stations to demonstrate the effectiveness of the solution algorithm and the accuracy of the proposed two-stage stochastic model.


Introduction
Bike sharing system (BSS) as a means of sustainable and carbon-free transportation can effectively solve the "the first and last mile" problem of urban public transportation. BSS has gained increasing popularity in recent years due to the fact that it not only reduces urban pollution emissions and traffic congestion but also is considered as an effective way to improve the health of users [1]. A BSS consists of a depot and several stations, which scatter in different streets of the city. e success of a BSS depends on the availability of bikes. Due to the rental cycling of the user, there are often unbalanced situations at some stations of the system, i.e., there are either surplus or not enough bikes at these stations. In other words, either when there are surplus bikes in a station, they will be wasted as no users will need them, or in some stations where there are not enough bikes, a user's demand cannot be met. To ensure the availability of bikes at all stations in the system, the BSS needs to employ redistribution trucks to transfer bikes from stations to stations, to balance the number of bikes according to the demand in each station. is problem is commonly defined as the bike repositioning problem (BRP) in the BSS. e main objective of the BRP is to balance the supply and demand of bikes across stations as a mean to improve user satisfaction, while in the meantime, to reduce the transportation cost of redistribution trucks as much as possible. According to the classification of Berbeglia et al. [2], the bike repositioning problem is essentially a many-tomany pickup and delivery problem, and it is first proven to be an NP-hard problem by Benchimol et al. [3].
Typically, existing work assumes the redistribution demand at each station in the BSS is deterministic input parameters [4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. at is to say, these models heavily rely on that the loading and unloading quantities of redistribution trucks at each station be given in advance. However, this strong assumption generally could not hold in real-world conditions where the demand in urban public transportation commonly contains uncertainty and randomness [22]. As pointed out by Fricker and Gast [23], the effect of the uncertainty in user demand should not be neglected when studying the performance of a BSS. e authors have also demonstrated that incorporating knowledge of future demands can lead to more accurate and robust decisions in repositioning.
e other challenge that might not have received enough attention in BRP is the consideration of the holding cost of the depot. In the BRP, bikes transferred from a bike surplus station to another bike deficient station is required to achieve the balance of supply and demand at all bike stations in the BSS to satisfy the user's demand.
is perfect balance of supply and demand after repositioning operation must first ensure the internal balance in the BSS, i.e., the number of bikes at bike surplus stations should be equal to the number of bikes at bike deficient stations. In the actual operation process, some bikes might be lost or damaged after they are put into use. is happens quite often and will likely lead to the inherent imbalance between supply and demand in the BSS. For the bikes which cannot be balanced among the stations in the BSS, they are either picked up from or delivered back to the depot at the cost of increasing holding cost of the depot. Some BRP models have simply assumed the capacity of the depot is sufficient. As such, some bikes can be picked up from the depot to meet the demand of bike deficient stations, and/or some bikes at bike surplus stations can be brought back to the depot, in order to achieve a perfect balance between supply and demand at all stations in the BSS [4][5][6][7]. Unfortunately, these existing works have not yet considered the holding cost of the depot in the BRP model.
Stochastic programming is a technology to design mathematical programming models with stochastic demand, which helps decision makers to make more accurate decisions when the effect of the randomness of demand cannot be ignored. Although the optimization results obtained from stochastic programming cannot satisfy all future scenarios, they help to achieve better-than-average performance in dealing with stochastic demand.
To address the abovementioned concerns, in this paper, we take into account the stochastic nature of the BRP by capturing and modeling the uncertainty of redistribution demand within the system using stochastic programming. In particular, we propose a two-stage stochastic programming model for the bike repositioning problem with stochastic demand (BRPSD). In the first stage, we determine the routing decision, where "here-and-now" decision must be made before the realization of redistribution demand is known. en, moving to the second stage, we further determine the loading/unloading decisions at each station and depot, where "wait-and-see" decisions are made taking future uncertainty into account. We incorporate the holding cost into the model such that the overall objective of the BRPSD is to determine the best routes of the repositioning truck and the optimal loading/ unloading quantities at each station and depot, such that the expected total sum of the transportation costs, the penalty costs at all stations, and holding cost of the depot are minimized.
In this paper, stochastic redistribution demands are modeled by a set of discrete scenarios, and a predefined probability of occurrence is given for each scenario. e main contributions of this paper are summarized as follows: (1) A general scenario-based two-stage stochastic programming model for the bike repositioning problem with stochastic demands (BRPSD) is introduced. e proposed model takes the inventory holding cost of the depot into account in contexts where the loading and unloading activities are allowed in the depot. (2) An efficient metaheuristic, namely, simulated annealing algorithm is proposed to solve the defined BRPSD. (3) e necessity of considering redistribution demand as a stochastic demand in the BRPSD is further discussed with numerical experiments. e remainder of this paper is organized as follows. We review the relevant literature in Section 2 and describe the mathematical formulation of the problem in Section 3. In section 4, we present the improved simulated annealing algorithm to solve the model. en, we present experimental results in Section 5 and conclude this paper with discussions on potential future research in Section 6.

Literature Review
To the best of our knowledge, Benchimol et al. first formalise and address the BRP in [3]. Since then, various models and solution algorithms for the BRP have been discussed over the last decades.
We briefly review studies that are related to this paper, including heuristic and metaheuristic solution methods, for solving the deterministic BRP and models for the BRP with stochastic demands.

Heuristic Solution Algorithm for the Deterministic BRP.
According to Pal and Zhang [8], the BRP can be categorized into complete and partial repositioning problems according to the rigor of a repositioning that needs to be performed. In the former problem, the constraint for achieving a perfect balance among stations is considered as a hard constraint, and in the latter problem, it is a soft constraint.
In the complete repositioning problem, the optimal inventory of each station must be met after the repositioning operation. An iterated local search algorithm was first applied by Cruz et al. [9] to solve the single-vehicle BRP and then used by Bulhões et al. [10] to solve the multivehicle version of the BRP. Later, Pal and Zhang [8] modeled a static BRP in a free-floating bike sharing system and presented a hybrid nested large neighborhood search with a variable neighborhood descent algorithm to tackle the BRP. Dell'Amico et al. [11] developed a destroy and repair metaheuristic algorithm, which makes use of effective constructive heuristic and several local search procedures together for solving the BRP. Overall, the objective of the complete repositioning problem is mainly to pursue the minimum transportation cost by determining routing decisions. It has been a challenge to apply such methods in real-world systems effectively as there are other equally important criteria that have been neglected in these models.
In the partial repositioning problem, the ideal inventory of each station is not necessarily achieved after the repositioning operation. Nonetheless, the goal of these problems is not only to pursue the minimum transportation cost by determining routing decisions but also to minimize user dissatisfaction by determining loading and unloading decisions at each station. In the literature, a 3-step math heuristic for the multivehicle BRP is proposed by Forma et al. [6]. In their approach, stations are firstly clustered by a saving heuristic. en, the routes through the clusters are constructed. Finally, the routing sequence through the stations in each cluster and the sequence through the clusters are determined. e last two steps are formulated as a mixed integer linear model solving by CPLEX. You [12] proposes a two-phase heuristic for the multivehicle BRP under a minimum service requirement over planning. In his proposed model, each stage represents a decision process. In the first phase, loading and unloading decisions for all stations for each time slot are developed by a linear programming model. en, in the second phase, routing decisions are developed by an iterative approach through two parameter sensitive mathematical models.
In addition to the phased heuristic algorithm, many heuristics and metaheuristic solution methods have been presented to solve the deterministic BRP, to name a few, tabu search [5,13], genetic algorithm [14], chemical reaction optimization algorithm [15,16], artificial bee colony algorithm [17], constraint programming [18,19], variable neighborhood search [20], iterated local search [21], large neighborhood search algorithm [7], etc. Nonetheless, the practical applications of these existing approaches remain difficult and often bring in unexpected results. is is because uncertainty about user demands commonly exists in real-world BRPs; however, it is not yet captured by such deterministic models.

e BRP with Stochastic Demands.
To the best of our knowledge, limited literature considers stochastic demands in the BRP. Two-stage stochastic programming with recourse models have been successfully used in solving bike repositioning problems as they allow the modeler to represent routing plans and loading and unloading activities together via first-and second-stage decision variables. Dell'Amico et al. [24] appear to be the pioneer paper that proposed a scenario-based two-stage stochastic model for bike repositioning problem, considering redistribution demands of stations as random variables. e authors provide five exact procedures as well as a heuristic algorithm that combines correlation-based constructive procedures with VND local search approach to solve their BRPSD. e heuristic algorithm only considers transportation costs while evaluating the local search move and only accept a feasible solution where the feasibility of the local search move is checked by an efficient strategy.
However, the BRPSD differs from our problem in the following three aspects. First, the objective of this model is to minimize the expected total travel distance of the truck and the penalty costs at all stations, while our problem minimizes the expected total travel cost of the truck, the penalty costs at all stations, and the handling cost at the depot. Second, in our problem, we measure the local search move by using not only transportation costs but also the expected total costs. ird, we develop a method to quickly determine the feasible loading and unloading quantities at all stations and depot under any given route. erefore, it is not necessary to check the feasibility of the solution in the local search.
ere are also a few recent works that study a variation of the stochastic BRP where the stochastic demand is defined by origin-destination (OD) pair rather than the redistribution of stations. For instance, Maggioni et al. [25] propose two-stage stochastic programming models to determine the optimal number of bikes to assign to each station and the optimal number of transshipped bikes from one station to another station, respectively. ey use AMPL and CPLEX to solve the model. Yan et al. [26] apply the time-space network model to determine the optimal locations of bike stations, bike fleet allocation, and bike routing. eir solution algorithm is based on a threshold-accepting-based heuristic.
From the abovementioned review of literature, we can see that the current scholars have proposed a rich heuristic algorithm to solve the deterministic BRP. Nonetheless, to the best of our knowledge, there is no work yet that effectively employs a metaheuristic algorithm to solve the BRP with stochastic demand. e study of this paper fills this theoretical gap.

Mathematical Formulation
e bike sharing system studied in this paper comprises a depot and a set of stations. Each station has a redistribution demand for each scenario, indicating the difference between the current inventory level and the optimal inventory level for each scenario. For a certain scenario, if the redistribution demand at a station is positive, the station is defined as a pickup station; if the redistribution demand at a station is negative, then the station is defined as a delivery station and should be supplied with bikes from pickup stations; if the redistribution demand at a station is equal to 0, then the station is called a balanced station. Balanced stations must also be visited.
We consider that only a single truck with a given capacity is available in the BSS. e truck collects bikes from pickup stations and transports them to delivery stations. Some bikes can be picked up from or delivered back to the depot, and each station is allowed to be visited only once. In addition, similar to Liu et al. [16], bikes are assumed to be parked anywhere in the BSS; therefore, the station capacity constraint is removed. e decisions regarding the loading and unloading quantities at all stations and the depot in each scenario are made independently, according to the corresponding demand scenario. However, the routing decisions of the truck have to take into consideration all possible demand scenarios. is means that routing decisions are the same for all scenarios. Because the routing decisions to all scenarios are the same but the loading and unloading decisions vary in each scenario; there may exist a shortage or excess of bikes Scientific Programming for each station in each scenario. In this context, the following two-stage bike repositioning problem with stochastic demands (BRPSD) naturally arises. In the first stage, the BRPSD aims to construct a route for the truck that starts from and ends at the depot where each station is visited only once. In the second stage, the BRPSD aims to determine the optimal loading and unloading quantities at all stations and the depot to ensure the feasibility of the solution after the given route in the first stage. e objective of the BRPSD is to minimize the sum of the travel costs for the truck (first-stage objective function) and the expected total penalty costs at all stations and holding costs at the depot (second-stage objective function).

Two-Stage Stochastic Programming Model.
Our objective is to determine the sequence of all stations to visit to minimize the transportation cost, as well as the loading/ unloading decisions in each scenario to minimize the total expected penalty costs and holding costs. We formulate the problem as a two-stage stochastic programming problem. Given the discrete probability distribution of the scenario occurrences, the formulation for the SPBSD addressed in this study is as follows: e objective function (1) is defined as the weighted sum of the transportation costs and the expected penalty costs at stations and holding costs at the depot over all scenarios. Constraint (2) states that the truck departs from the depot only once and returns to the depot. Constraint (3) ensures that the truck can visit a station only once. Constraint (4) ensures that if the truck visits a station, it must leave that station. Constraint (5) is the subtour elimination constraint [27]. Constraint (6) ensures that the load on the vehicle cannot be greater than the truck capacity. Constraint (7) requires that the number of bikes loaded onto or unloaded from the truck at a given station equals the difference between the truck load before and after the station visits. Constraint (8) ensures the auxiliary variable is nonnegative. Constraints (9)-(10) ensure that the excess/slack quantities of bikes at station and the number of bikes on the truck are nonnegative integers. Constraint (11) defines x ij to be a binary variable.

Solution Algorithm
When the number of scenarios equals 1, the studied BRPSD reduces to the classical deterministic bike repositioning problems, an NP-hard combinatorial optimization problem. We employ simulated annealing (SA) to solve the studied BRPSD. SA is a metaheuristic algorithm based on local search, which can avoid falling into local optimization by accepting poor solutions with less probability in the process of iteration. It has been successfully applied to a variety of deterministic combinatorial problems [28,29] and stochastic combinatorial problems [30,31].
Our proposed algorithm, called SA BRPSD , consists of two phases: constructing the initial solution and improving the initial solution by neighborhood search mechanism. e structure of this section is organized as follows. In Section 4.1, we propose the heuristic method to generate an initial solution. During initialization, two main decisions must be made: the routing decision (in Section 4.1.1) and the loading and unloading decision (i.e., inventory decision) (in Section 4.1.2). en, in Section 4.2, we describe in detail of the local search operators used in this SA BRPSD , namely, Swap, Relocate, and 2-opt. After that the main algorithm flow based on the SA BRPSD framework is presented in Section 4.3.

Construct Initial Solution.
Since the BRPSD is a twostage decision-making process, the decision-making in this paper is divided into two stages: routing decisions are determined in the first stage, and the loading and unloading decisions are determined in the second stage. e process of constructing the initial solution is as follows, which is roughly divided into two main phases.

Construct the Initial Route for all Scenarios.
Since the route in all scenarios is the same, the greedy heuristic algorithm is used to quickly construct the initial route by ignoring the stochastic redistribution demand. In this work, the nearest principle is used to construct the initial route by inserting stations one by one.

Determine Loading and Unloading Decision for Each
Scenario. For a given scenarioξ and the corresponding redistribution demand d ξ i at station i in the scenario, the Second-Stage-Opt procedure uses a simple principle to quickly generate the initial feasible loading and unloading decision and then uses a reoptimize heuristic operator to improve the initial loading and unloading decision. e details of the Second-Stage-Opt procedure are given as follows.
Step 1: determine the initial loading and unloading quantities of all stations and the depot. e loading and unloading quantities at each station are defined as the corresponding redistribution demand. Because the demand of the depot is 0, the loading and unloading quantities are also defined as 0. As a result, the initial solution of the BRPSD is generated, including twostage decisions, namely, the routing decision and the loading and unloading decision. It is obvious that the initial solution is often not feasible.
Step 2: adjust the loading and unloading quantities of some stations to achieve a feasible initial solution. If the number of bikes on the truck when it travels directly from station i to station j in scenarioξy ξ ij > Q or y ξ ij < 0, then the initial solution is not feasible. In this case, the loading and unloading quantity of station i should be adjusted.
if y ξ ij < 0, then the loading and unloading quantity of station i is adjusted to subtract cfrom the original quantity.
Step 3 (reoptimize): the basic idea of the reoptimize operator is to optimize the solution by a heuristic method that adjusts the loading and unloading quantities between stations and the depot [5,7]. In this paper, some bikes can be collected from the depot to meet the delivery demands of some stations, and also some bikes collected from some stations can be delivered to the depot at the end of the repositioning. In this procedure, there are two cases. In the first case, the truck picks up extra bikes when it starts from the depot and drops them off at some stations. In the second case, the truck picks up more bikes from some stations and drops them off when it ends at the depot.

Neighborhood Structure.
In this paper, the local search operator aims at routing decisions. ree kinds of classical operators Swap move, Relocate move, and 2-opt move are used.
ese operators are also often used in other BRP problems. Each new neighborhood solution is generated by these three moves with equal probability. Note that, after a local search operator is executed, the new solution needs to be improved by Second-Stage-Opt procedure to obtain the second-stage decision.
Swap-N 1 : randomly select two stations from the route and then exchange them Relocate-N 2 : randomly select one station, delete it from the current route, and then insert it into another position 2-Opt-N 3 : randomly select two stations, exchange the positions of the two stations, and change the orientation of the route between the two stations

4.3.
e Proposed SA BRPSD . Algorithm 1 provides the pseudocodes of the SA BRPSD . Parameter T 0 denotes the initial temperature; T E is the final temperature; α is the coefficient used to control the cooling; I iter indicates the number of generated new solutions under a certain temperature; and Max iter is the max number of the accepted new solutions under a certain temperature. In this section, we use S 0 and S best to denote the initial solution and the best solution obtained so far, respectively. e best objective value F(S best ) is set to be the objective function value of the best solution S best .
As shown in Algorithm 1, in the inner loop of the SA BRPSD algorithm, a random new neighborhood solution S new is initially generated by Local search procedure to improve the first-stage decision and then followed by Second-Stage-Opt procedure to generate the second-stage decision. e stopping criteria for the inner loop are either it has generated I iter new solutions or accepted Max iter new solution under the current temperature T. When the current temperature T decreases to T E , the SA BRPSD algorithm terminates.

Numerical Examples
In this section, numerical examples are conducted to examine the validity and efficiency of the proposed model and solving algorithm. In addition, the necessity of considering stochastic demands is examined by analyzing the expected value of perfect information (EVPI) and the value of the stochastic solution (VSS).
In this study, we consider the instances as in [24] to conduct the numerical experiments (the instances are available at http://www.or.unimore.it/site/home/onlineresources). e size of the instances ranges from 20 to 90 stations. ere are a total of 20 instances and 30 scenarios for each instance. In any instance, the first station is chosen to be the depot. To better analyze the impact of relevant parameters on the solution of the stochastic model, for a given set N 0 , station i and station j belong the set N 0, and we redefine the parameter transportation time t ij as t ij � t ij /min {t ij }. e SA BRPSD was coded in Matlab, and all computational experiments were carried out on a computer with an Intel Core i5-4590 CPU @ 3.30 GHz and 4 GB RAM. Each instance was run 10 times, the best and average of the solutions, and the average computing time of 10 runs were used to evaluate the performance of the algorithm solution.
e SA BRPSD proposed in this study relies on five parameters, namely, the initial temperature T 0 ; the final temperature T E ; the coefficient used to control the cooling α; the number of generated solutions under a certain temperature I iter ; and the number of accepted solutions under a certain temperature Max iter . After some preliminary experiments, the results from the algorithm were obtained by e CPU (in seconds) indicates the average computing time of the SA BRPSD or the computing time for Lingo, respectively. Gap Best (%) and Gap Avg (%) indicate the performance of the SABRP relative to that of Lingo based on the best and average objective values, respectively. When |N 0 | ≥ 50, Lingo is not able to obtain any feasible solutions for the problems, and thus, the results are not shown in the tables. It also generally fails to obtain any global optimal solutions within 3600 seconds even for small-size problems (with four exceptions for the results obtained under instance Washington (20)2). From Tables 1-4, we can see that SA BRPSD is significantly faster than Lingo in all instances, i.e., it has better time efficiency to solve the BRPSD compared to Lingo. For the various sizes of problems from small to large, the average computing time of the SA BRPSD varies from 66.14 s to 70.22 s.
As shown in Tables 1-4, all of the values of Gap Best and Gap Avg are negative. is implies that the SA BRPSD performs better than Lingo. e absolute values of the average Gap Best and Gap Avg values obtained by the SA BRPSD for the tests with Q � 20 are larger than those with Q � 10. As the capacity of the repositioning truck increases, the objective function value decreases. is may be explained by the fact that if the truck capacity increases, the truck can load/unload more bikes at one station, resulting in an expansion of the solution space. Overall, we can draw the conclusion that the proposed SA BRPSD can obtain better feasible optimal solutions much faster than Lingo.

e Necessity of Considering Demand as a Stochastic
Parameter. To evaluate the performance of the stochastic model we proposed in this paper, we use the two concepts of the evaluation approach: the expected value of perfect information (EVPI) and the value of the stochastic solution (VSS) [32][33][34]. EVPI measures the maximum amount a decision maker would be ready to pay in return for perfect information about the future, and VSS measures the amount of cost saving when the decision maker uses expected values of stochastic parameters instead of stochastic parameters in the model (the amount of cost due to ignoring the uncertainty). EVPI is defined as the difference between the waitand-see solution (WS) and the here-and-now solution (HN). In the EVPI, the WS solution is the expected value of the foresight. We obtain the value by averaging the results of each stochastic scenario. e HN solution is the objective value of this two-stage stochastic model. VSS is defined as the difference between the expected results of using the expected value problem solution (EEV) and the here-andnow solution. In the VSS, the value of EEV is obtained by inputting the optimal solution of the optimization problem associated with the mean value of stochastic demands into the two-stage stochastic model. e gaps Gap VSS and Gap EVPI for VSS and EVPI are calculated by VSS/HN and EVPI/WS, respectively.
Obtain the initial solution S 0. Set the parameter for the simulated annealing, T 0 , T E , α, I iter , Max iter Set S best ⟵ S 0 ; S current ⟵ S 0 ; T ⟵ T 0.
If e − Δ/T ≥ random then S current ⟵ S new , n ⟵ n + 1; end end m ⟵ m + 1; end T ⟵ α ⃞ T; end Return: S best ALGORITHM 1: SA BRPSD . 6 Scientific Programming e test instance Chicago 20(1) used for explaining the total costs of WS, HN, EEV, EVPI, and VSS is a small-size instance consisting of 20 stations. We set the truck capacity Q � 10, C h � 1, C d � 1, and C t � 1.
In this test instance, the total costs can be solved by the SA BRPSD as follows: WS � 145.31, HN � 153.40, and EEV � 161.31. en, calculate VSS � EEV − HN � 7.90, EVPI � HN − WS � 8.09, Gap VSS � 5.15%, and Gap EVPI � 5.57%. Both values of EVPI and VSS lead to the conclusion that managers must pay 5.57% more of the total cost of the stochastic solutions to obtain the perfect redistribution demand information about the future. Ignoring the uncertainty will let managers in a situation where they have a risk to pay at least 5.15% higher if they only focus on a scenario only (deterministic solutions). e results of the sensitivity analysis on Gap VSS and Gap EVPI for the truck capacity Q, the holding cost C h , the penalty cost C d , and the travel cost C t are explained in the

Scientific Programming
following. e value of the four parameters are increased by −50%, −25%, +25%, and +50%, respectively. Tables 5-8 summarize the total costs of WS, HN, EEV, EVPI and VSS, GAP EVPI , and GAP VSS obtained from all tests. Figure 1 shows the sensitivity analysis on the results for the four parameters. Under current settings, the abovementioned four parameters have a significant impact on the solutions. From Tables 5-8, consistent with theoretical results, it can be founded that WS < HN < EEV in all tests. We can see that the average GAP VSS are all greater than the average GAP EVPI . e GAP EVPI of the stochastic model in all tests is nonnegative, indicating the problem of the BRPSD is stochastic by nature. Besides, the GAP VSS in all tests is nonnegative for the BRPSD, indicating that the problem in the study is suitable for the stochastic model.
From Figure 1, we can see that the increase of C h , C d , and C t rends the increase of HN. We also see that the increase of Q rends the decrease of HN. e parameter that has the greatest impact on HN is the unit transportation cost C t . is is likely because the contribution of transportation cost to the total objective cost is greater than the penalty cost and the holding cost.
From Figure 1, we can see that, as Q increases, GAP EVPI tends to decrease, but GAP VSS increases. It is also noticed that GAP EVPI and GAP VSS tend to increase when C h and C d increase. However, the value of GAP EVPI and GAP VSS contradicts the increase of C t . e influence of the change of the abovementioned four parameters on GAP VSS is significantly greater than that of GAP EVPI . e results illustrate that the GAP EVPI increases when C h and C d increase and Q and C t decrease, which means that the decision makers are willing to pay more to get accurate redistribution demand information about the future. GAP VSS increases by increasing the values of C h , C d , and Q or decreasing the values of C t , and consequently, the decision makers can obtain more cost saving resulting from the twostage stochastic model.

Conclusions
In this paper, a general scenario-based two-stage programming model is proposed for bike repositioning problems with stochastic redistribution demands. In the proposed model, the first-stage decisions correspond to the routing decisions and the second-stage decisions correspond to the loading and unloading decisions. e model aims to find the best routes and the ideal quantities of loading and unloading at each station and depot in order to minimize the weighted sum of the expected transportation costs, penalty costs, and holding costs. en, a simulated annealing algorithm is proposed to solve the model. Number examples are performed to evaluate the proposed model and algorithm, followed by a detailed sensitivity analysis that studies how the change of several important parameters affects the performance of the proposed model and solution algorithm. e proposed model and algorithm have important theoretical and practical significance for the BSS operators, which can reduce the operation costs of the BSS and improve user satisfaction. On the contrary, in addition to bike sharing, other sharing facilities in the sharing economy, such as car sharing, also need to be relocated, so this kind of problem has a strong application background.
In the future, we will apply robust optimization model to solve the bike repositioning problem with stochastic demand in the bike sharing system [35]. We also consider environmental issues by adding the carbon emissions to the objective function of the two-stage stochastic programming model [36][37][38]. Decision Variables x ij : Binary variable that equals to 1 if the truck travels directly from station i to station j and 0 otherwise y ξ ij : e number of bikes on the truck when it travels directly from station i to station j in scenario ξ s ξ+ i : e excess quantities of bikes at station i after relocation operation in scenario ξ s ξ− i : e slack quantities of bikes at station i after relocation operation in scenario ξ q i : Auxiliary variables associated with station i used for subtour elimination constraints.

Data Availability
All data included in this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.