1. Introduction
Bicycle sharing systems (BSSs) have been introduced in many cities around the world as a new public transportation system. Currently, more than 1,700 programs of BSSs are operating in the world [
1], for instance, BIXI in Montréal [
2], Velib’ in Paris [
3], docomo bike share in Tokyo [
4], citybike in Vienna [
5], Citi Bike in New York [
6], Capital bikeshare in Washington D.C. [
7], BULEbikes in Boston [
8] and Bicing in Barcelona [
9].
A BSS consists of many distributed automatic rental ports and many rental bicycles. A registered user can rent a bicycle at a port, use it for a short trip and return it to a different port. A BSS is not only a convenient service, but also it can reduce greenhouse gas and traffic congestion issues. For example, Zhang et al. (2018) [
10] have estimated that the introduction of a BSS at Shanghai saved 8358 tonnes of petrol, and reduced CO
and
emissions 25,240, and 64 tonnes, respectively in 2016.
In a BSS, almost all ports are in either an empty or full state after a while. An empty bicycle port implies that a user cannot rent a bicycle, whereas a full bicycle port implies that a user cannot return a bicycle. Thus, full or empty states lead to dissatisfaction of users. In addition, broken bicycles are left unattended. Consequently, users’ frustration is further increased. To address these issues, several capacitated vehicles need to relocate bicycles among ports so that sufficient bicycles are available at each port. If operators can execute repositioning works efficiently, users will not have experience that they cannot rent or return bicycles and BSS companies can reduce labor costs and carbon emissions of vehicles. Thus, both users and BSS companies benefit.
Recently, many studies have been conducted on repositioning bicycles of a BSS. The repositioning of BSS can be divided into two cases.
The first one is a static case. In the static case problem, the number of bicycles does not change at each port while the bicycles are restored by the vehicle, because the system assumes that the bicycles are not moving much at midnight or after operation hours. As for the static case problem, Dell’Amico et al. (2014) [
11] proposed a bike sharing rebalancing problem (BRP) and formulated BRP as a mixed integer linear programming problem, and they solved BRP using the branch-and-cut algorithm for the 22 real world instances collected from the web sites. Then, they proposed a destroy and repair algorithm for the BRP, and tested their algorithm on benchmark instances [
12]. Although these studies [
11,
12] are interesting, they did not consider a time limit constraint: vehicles must return to a depot within a given time limit. Then, it is not guaranteed that all operations will be finished within the working hour in BRP.
On the other hand, Rainer-Harbach et al. (2013) [
13] proposed a general variable neighborhood search (VNS) with an embedded variable neighborhood descent for finding efficient vehicle tours of a balancing bicycle sharing system (BBSS). In the BBSS, they considered the time limit constraint, and that a vehicle can visit a port multiple times. Then, they tested their algorithm on benchmark instances of real-world data provided by Citybike in Vienna with up to 90 ports, and the result of numerical experiments showed that their algorithm obtained good solutions in a reasonable time. In addition, Raidl et al. (2013) [
14] proposed efficiently determining optimal loading operation of VNS method for finding efficient vehicle tours of the BBSS. They tested the same instances with Rainer-Harbach et al. (2013), and they reported that their method showed better performance than the original VNS. Moreover, Papazek et al. (2013) [
15] also address the BBSS, and they proposed two GRASP (Greedy Randomized Adaptive Search Procedure) heuristic approaches and conducted numerical experiments by varying sizes of instances that ranged from 30 to 700. The numerical experiments show that combined PILOT (Preferred Iterative Look Ahead Technique) and GRASP method obtained good performance than VNS for large size instances. In the BBSS, when vehicles start and come back to a depot, the vehicles must be empty. However, if broken bicycles are found in the process of rebalancing operation, they must be repaired at the depot after the rebalancing operation, because the depot has many spare bicycles. Thus, the constraint that vehicles must be empty when vehicles start and come back to the depot is not realistic [
13,
14,
15].
The second one of the repositioning of the bicycle sharing system is a dynamic case. In the dynamic case problem, the number of bicycles is changed at each port while the bicycles are restored by the vehicle. Kloimüllner et al. (2014) [
16] extended the BBSS to a dynamic case problem, and proposed a model of a dynamic balancing bicycle sharing system (DBBSS). Then, they applied their previous approach to the DBBSS. The numerical experiments showed that their approach was also effective to DBBSS. Chiariotti et al. (2018) [
17] also addressed the dynamic case problem. Then, they tested instances with 280 ports generated by New York City’s CitiBike, and the result of numerical experiments showed that their dynamic strategy obtained better performance than static approach. In the numerical experiments for the dynamic repositioning of BSS, Kloimüllner et al. (2014) and Chiariotti et al. (2018) used Vienna and New York City’s data, respectively. However, it is not clear whether their approaches can be applied to other cities.
Table 1 summarizes the related literatures.
In our previous study, we formulated a new bicycle sharing system routing problem as “multiple-vehicle bike sharing system routing problem (mBSSRP).” The mBSSRP finds tours within the shortest travel time when multiple vehicles load and unload bicycles at each port. In the mBSSRP, we consider the time limit constraint, and vehicles can start and finish at the depot with bicycles. Although we treated the static case problems, the mBSSRP is close to a real-life situation. As a result of solving the mBSSRP using a mixed-integer linear programming (MIP) solver [
18], an optimal solution is obtained for small-size instances. However, we cannot obtain an optimal solution for large-size instances within a reasonable time frame. Thus, we proposed local search methods to obtain approximate solutions of mBSSRP within a short time [
19]. Nevertheless, owing to some strict constraints of mBSSRP, the proposed methods could not obtain feasible solutions for some instances. Therefore, we proposed a “multiple-vehicle bike sharing system routing problem with soft constraints (mBSSRP-S) ” that removes some constraints from the mBSSRP and appends violations of the mBSSRP to the objective function as penalties [
20]. In Ref. [
20], we also proposed a method that controls the execution of the Or-opt [
21] and the CROSS-exchange [
22] with the tabu search method [
23,
24,
25] for solving mBSSRP and mBSSRP-S. Numerical experiments show that solving mBSSRP-S to obtain feasible solutions of mBSSRP achieves good performance when compared to solving directly mBSSRP [
20]. However, computational costs increase because neighborhood solutions increase as a result of including infeasible solutions of mBSSRP. Thus, it is necessary to reduce the computational costs while maintaining the performance.
Regarding the reduction in computational costs, there are several discussions for the combinatorial optimization problems. For example, Hasegawa et al. (2002) [
26] and Matsuura et al. (2008) [
27] introduced a method of limiting the search space to short paths when they solved the TSP (Traveling Salesman Problem). In particular, Hasegawa et al. (2002) used only the 10% nearest neighbors of the cities. Then, Matsuura et al. (2008) used only the near neighbors of the cities. This is because most good TSP solutions have short paths. Motohashi et al. (2008) [
28] used a neighborhood list to efficiently solve the TSP. Moreover, Motohashi et al. (2009) [
29] used two different candidate lists, that is, the nearest neighbors and quadrant neighbors to solve a large scale TSP. Finally, Taillard et al. (2019) [
30] showed that the computational costs of heuristic methods can be significantly reduced by constructing a set of edges that are considered to be included in the optimal solution at first.
Regarding the BSS, Pal et al. (2016) [
31] used a candidate list based on nearest neighbors to efficiently solve the static complete rebalancing problem. Then, Dell’Amico et al. (2016) [
12] used some formulations to determine whether constraints are satisfied or not to speed up the computation of bike sharing rebalancing problem. Thus, it is important to reduce the computational costs for finding a good rebalancing tour of the BSS, because a BSS has hundreds of ports and it is not practical to wait for hours to find a rebalancing tour of the BSS. Then, the goal of this study is to propose efficient search strategies with low computational costs while maintaining the performance of the mBSSRP.
To reduce computational costs in this study, we focused on the search process and investigated the rate of the Or-opt and the CROSS-exchange executions. The execution rate of the Or-opt that inserts ports in the normal order was 66%, that of the CROSS-exchange that exchanges ports in the normal order was 20%, and the reverse orders of insertion of the Or-opt and exchange of the CROSS-exchange were rarely executed. In addition, we found that once a feasible solution of the mBSSRP is obtained, our conventional method [
20] searches for the feasible solutions successively.
From these results of the investigation, in this paper, we propose two search strategies: the first one is to reduce neighborhood solutions to obtain a feasible solution in a short time before finding a feasible solution of the mBSSRP from the mBSSRP-S, and the second one is to change the problem to be solved (mBSSRP or mBSSRP-S) after a feasible solution is obtained and to search good near-optimal solutions in a short time. As the first search strategy, we propose two search methods for reducing the number of neighborhood solutions in the Or-opt and the CROSS-exchange, and compare their performance with our previous results. The results of the numerical experiments indicate that by limiting the neighborhood solutions to the normal order insertion of the Or-opt and the normal order exchange of the CROSS-exchange, we obtained feasible solutions within a short time while maintaining the performance.
Next, as the second search strategy after obtaining a feasible solution, we proposed four methods in which the searching space is further reduced. The problem to be solved is changed from mBSSRP-S to mBSSRP to explore only feasible solutions. The results of the numerical experiments show that by reducing the number of the neighborhood solution, only searching for the normal order insertion of the Or-opt and the normal order exchange of the CROSS-exchange with hard constraints after obtaining a feasible solution can obtain short tours within a short time. From these numerical experiments, we have succeeded in reducing the computational time while maintaining the performances. Furthermore, we found that after finding a feasible solution, we could find short tours within a short time by solving the hard constraint problem instead of the soft constraint problem.
This paper is structured as follows. In
Section 2, we describe the mBSSRP and mBSSRP-S. In
Section 3, we present the heuristic method. Then,
Section 4 presents the computational costs reduction method. Numerical experiments are reported in
Section 5. Finally, conclusions and future works are provided in
Section 6.
4. Computational Cost Reduction Method
To reduce neighborhood solution and computational costs, we investigate the percentage of the Or-opt and the CROSS-exchange executions used in the tabu search part. The feasible neighborhood solutions of the mBSSRP-S include the infeasible neighborhood solutions of the mBSSRP. Therefore, when we solve the mBSSRP-S, it takes more time to determine the next solution than in the case of mBSSRP. To reduce the computational cost of solving mBSSRP-S, we need to remove unnecessary neighborhood solutions. The percentage of the insertion and exchange operations performed by the proposed method was approximately 66% for the Or-opt, which inserts ports in the normal order (
Figure 4b), and approximately 20% for the CROSS-exchange, which exchanges ports in the normal order (
Figure 5b), implying that the reverse order of insertion and exchange was rarely executed. Thus, we did not search for less frequently performed reverse order operations to reduce the computational costs. In our previous study, we have shown that once a feasible solution of the mBSSRP is obtained, the proposed method searches for the feasible solutions successively [
20]. Next, to find a good feasible solution of the mBSSRP, we examine whether we should change the problem to be solved (mBSSRP or mBSSRP-S) before and after a feasible solution is obtained.
Table 2 summarizes the search strategy before obtaining a feasible solution of the mBSSRP. We propose two search methods (1B-S and 1C-S) in this case. In
Table 2, the first search method 1A-S searches for all neighborhood solutions of the Or-opt and the CROSS-exchange. Thus, this search method is the same as the conventional search method [
20]. The second search method 1B-S only explores insertion of the normal order of the Or-opt which is the most executed (
Figure 4b). The third search method 1C-S searches for the insertion of the normal order of the Or-opt and exchange of the normal order of the CROSS-exchange (
Figure 4b and
Figure 5b). The goal of the search methods before obtaining a feasible solution is to find a feasible solution within a short time for all trials.
In addition, to further reduce the computational costs the search strategy is changed after a feasible solution of mBSSRP is obtained. In this paper, we propose four new search methods (2B-S, 2B-H, 2C-S and 2C-H).
Table 3 shows the search strategy after finding a feasible solution of the mBSSRP. In
Table 3, the first search method 2A-S involves to searching for all neighborhood solutions of the Or-opt and the CROSS-exchange with soft constraints. The second search method 2A-H involves searching for all neighborhood solutions of the Or-opt and the CROSS-exchange with hard constraints. The difference between the search methods 2A-S and 2A-H is whether the instances are solved with hard or soft constraints. The reason why we introduced both soft and hard constraints is that we found that once a feasible solution of the mBSSRP is obtained, our conventional method [
20] searches for the feasible solutions successively. Then, 2A-S and 2A-H are the same as the conventional search method [
20]. The third search method 2B-S only searches for insertion in the normal order of the Or-opt with soft constraints. The fourth search method 2B-H only searches for insertion in the normal order of the Or-opt with hard constraints. The fifth search method 2C-S explores the insertion of the normal order of the Or-opt and the exchange of the normal order of the CROSS-exchange with soft constraints. Finally, the sixth search method 2C-H explores the insertion of the normal order of the Or-opt and the exchange of the normal order of the CROSS-exchange with hard constraints.
6. Conclusions
In this study, we proposed two search strategies with low computational costs while maintaining the performance of the mBSSRP: the first one is to reduce neighborhood solutions to obtain a feasible solution in a short time before finding a feasible solution of the mBSSRP, and the second one is to change the problem to be solved (mBSSRP or mBSSRP-S) after a feasible solution is obtained and to search good near-optimal solutions in a short time.
As the first search strategy before obtaining a feasible solution, we propose two search methods (the search methods 1B-S and 1C-S) and compare their performance with our previous method (the search method 1A-S) [
20]. The search method 1A-S searches for all neighborhood solutions of the Or-opt and the CROSS-exchange [
20]. The search method 1B-S searches only the normal order insertion of the Or-opt. The search method 1C-S searches for the normal order insertion of the Or-opt and the normal order exchange of the CROSS-exchange. The results of the numerical experiments showed that the search method 1C-S that only searches for the normal order insertion of the Or-opt and the normal order exchange of the CROSS-exchange is the best search method, because the search method 1C-S obtained many feasible solutions within a short time.
Next, as the second search strategy after a feasible solution of mBSSRP is obtained, we propose four search methods (2B-S, 2B-H, 2C-S and 2C-H) and compare their performance with our previous methods (2A-S and 2A-H) [
20]. The search method 2A-S searches for all neighborhood solutions of the Or-opt and the CROSS-exchange with soft constraints. The search method 2A-H searches for all neighborhood solutions of the Or-opt and the CROSS-exchange with hard constraints. The search method 2B-S searches for insertion in the normal order of the Or-opt with soft constraints. The search method 2B-H searches for insertion in the normal order of the Or-opt with hard constraints. The search method 2C-S searches for the insertion of the normal order of the Or-opt and the exchange of the normal order of the CROSS-exchange with soft constraints. Finally, the search method 2C-H searches for the insertion of the normal order of the Or-opt and the exchange of the normal order of the CROSS-exchange with hard constraints. The results of the numerical experiments showed that the search method 2C-H that searches for the normal order insertion of the Or-opt and the normal order exchange of the CROSS-exchange with hard constraints is the best search method, because the search method 2C-H obtained short tours within a short time. These results indicate that the search method 1C-S reduces the neighborhood solutions in the Or-opt and the CROSS-exchange and that the search space was shifted from soft constraints to hard constraints after a feasible solution is obtained is effective.
In this study, we used a dynamically changing method for weight parameters (
and
) to find good solutions of mBSSRP from the solution of the mBSSRP-S. To adjust these weight parameters
and
, we used an adjustment method (Equations (
3)–(
6)). Then, from the results of numerical experiments (
Figure 9a), it takes a long time to update the better solution for solving the mBSSRP-S (for example, 38.82 [s] for instance No. 1). The results of other instances are shown in
Table 9. From
Table 9, the search method 2C-S takes about 30∼80 [s] to find better feasible solutions. In this temporal duration of 30∼80 [s], the weight values
and
must be adjusted to appropriate values by Equations (3)–(6). If we can set optimal parameters for an initial state, the search method 2C-S can find a better solution more efficiently and quickly than the search method 2C-H. Thus, one of the future works is to determine the parameters by using a parameter optimization approach.
In the numerical experiment, we used randomly distributed ports. However, some bicycle sharing repositioning studies have used real BSS data. For instance, Rainer-Harbach et al. (2013) [
13], Raidl et al. (2013) [
14] and Papazek et al. (2013) [
15] use data from Vienna, and Chiariotti et al. (2018) [
17] use data from New York. Thus, it is an important work to use real BSS data in a future study.
Moreover, we need to address the dynamic cases because the number of bicycles in each port changes while restoring bicycles in the real BSS. The studies of static BSS (Rainer-Harbach et al. (2013) [
13], Raidl et al. (2013) [
14] and Papazek et al. (2013) [
15]) have been extended to dynamic BSS by Kloimüllner et al. (2014) [
16]. Thus, we also aim to extend this study to a dynamic case. The completion of the dynamic case would have a great impact on solving an issue inherent in BSS: users cannot rent or return a bicycle at an empty or a full port. A solution to this issue would lead to the development of BSS in the future, and it is expected that BSS will be used for a wide range of purposes such as commuting to work, daily life, and even in times of disaster. It will also have an effect to promote the spread of the next generation system called MaaS (Mobility as a Service), which has been attracting attention recently: all public transportation systems are connected by IT technology. To address the dynamic case, we have to consider forecasting future demands of bicycles from actual data and repositioning them based on the forecast data by investigating statistical characteristics of real data of BSSs [
32].
In our previous study [
20], we proposed a heuristic method by using the tabu search. Then, in this paper, we also used the tabu search to propose search strategies with low computational costs while maintaining the performance of the mBSSRP. However, it is important to use not only the tabu search but also the other approaches (for example, simulated annealing) and compare these approaches. In addition, we only focus on the operation of the rebalancing bicycles [
19,
20] and statistical analysis of BSS usage history [
32]. However, Kabak et al. (2018) [
33] focus on the location of bicycle ports. It is an important idea to evaluate the existing BSS ports and suggest alternative places for additional BSS ports.