A Fissile Ripple Spreading Algorithm to Solve Time-Dependent Vehicle Routing Problem via Coevolutionary Path Optimization

,e time-dependent vehicle routing problems have lately received great attention for logistics companies due to their crucial roles in reducing the time and economic costs, as well as fuel consumption and carbon emissions. However, the dynamic routing environment and traffic congestions have made it challenging to make the actual travelling trajectory optimal during the delivery process. To overcome this challenge, this study proposed an unconventional path optimization approach, fissile ripple spreading algorithm (FRSA), which is based on the advanced structure of coevolutionary path optimization (CEPO). ,e objective of the proposed model is to minimize the travelling time and path length of the vehicle, which are the popular indicators in path optimization. Some significant factors usually ignored in other research are considered in this study, such as congestion evolution, routing environment dynamics, signal control, and the complicated correlation between delivery sequence and the shortest path. ,e effectiveness of the proposed approach was demonstrated well in two sets of simulated experiments.,e results prove that the proposed FRSA can scientifically find out the optimal delivery trajectory in a single run via global research, effectively avoid traffic congestion, and decrease the total delivery costs. ,is finding paves a new way to explore a promising methodology for addressing the delivery sequence and the shortest path problems at the same time.,is study can provide theoretical support for the practical application in logistics delivery.


Introduction
e vehicle routing problem (VRP) and its extended variants have recently received great attention in many fields, which are required to improve delivery efficiency and reduce the delivery costs of logistics companies [1][2][3][4][5][6]. In its basic form, some essential elements are often used to construct the foundation of the VRP resolution, for instance, several depots and a set of customers [7]. Each customer is located in a different geographic position and demands a specific amount of goods, which are carried to the designative destinations by a fleet of vehicles [8]. e vehicles will depart from the depots and finally return to the same depots, if the given delivery assignments are completed. e typical constraints of VRP are to pose a strict requirement on the vehicle travelling trajectory. at is, all matched sites can be visited only once along with the planned vehicle route. e common optimization objective of VRP is to seek the optimal delivery sequence and the shortest path performed by the specific vehicle. It requires that all customers can receive the associated goods, aiming at cutting down the related travelling time and path length of the overall routes.
Due to the boosting number of vehicles in urban networks, regular traffic congestion has been a common phenomenon that imposes some negative impacts on social development and daily life [9]. Especially for the logistics industry, several inevitable challenges are brought to impede the further research of the VRP and its multiple variants. How to effectively solve VRP plays a considerable role in minimizing delivery costs and ecological environment protection [10]. For example, the daily traffic congestions may cause huge uncertainties for vehicle route planning, when the vehicles intend to travel across a planned path with shorter length but will waste a long time waiting at the blocked areas [9]. e fuel consumption highly depends on the length and travelling time of the route. Hence, traffic congestion tends to result in high fuel consumption and air pollution with CO 2 emissions because the engine must keep running during waiting periods. According to reliable data reported by the International Energy Agency, the transportation sector has contributed to the second-largest CO 2 emissions in 2015 [11]. In recent decades, how to reduce carbon emissions and fuel consumption has become a significant issue that should be paid more attention to VRP [12][13][14][15][16]. Under the pressure of additional costs and environmental protection, the current concerns of logistics companies are to exploit an appropriate approach for addressing the VRP [17][18][19].
Most of the research studies investigated the VRP proposed some ideal hypotheses. For example, the vehicle travel speed between two nodes is assumed to be fixed, which means that the travelling time is simply regarded as a constant value. Under the background of the time-dependent vehicle routing problem (TDVRP), it is impractical to be supposed in the real world because the vehicle speeds may vary with the changes of traffic condition. In other words, the travel speed between two nodes is timedependent and usually relies on the starting time. When the vehicle reaches a node, the related link and node that will be visited soon may transform into blocked from the accessible state, influenced by sudden traffic congestion. As a result, an extra detour strategy or waiting behaviour both make the travelling time very long, even though the planned path possesses the shortest length. e longer path between two sites may not experience traffic congestion during the same period. A long path is not consistent with large travelling time [20].
erefore, the vehicle route should take into account the changes in traffic conditions. Malandraki was the first to investigate TDVRP via exploring the impacts of traffic congestion on vehicle speeds [21]. Eglese et al. designed a timetable of the shortest paths based on a practical network, to present the expected route distance and travel time between two nodes [22]. Donati et al. developed an advanced multi-ant-colony system to minimize the total travel time, during a fleet of vehicles travelling at time-dependent speeds with a fixed capacity [23]. Kok et al. employed a modified Dijkstra algorithm and dynamic programming heuristic to implement an appropriate strategy to avoid traffic congestion in TDVRP [24]. Huang et al. developed the TDVRP with path flexibility (TDVRPPF) under the stochastic traffic conditions and gave the path decisions that relied on the departure time [25]. Kuo adopted simulated annealing to minimize the fuel consumption used in the route for TDVRP [26]. Ma et al. proposed a combined order selection to address the TDVRP with time windows for perishable product delivery [27]. Taş et al. studied the TDVRP considering soft time window and stochastic travel time, and applied tabu search and adaptive large neighbourhood search to prove the effectiveness of the proposed model [28]. Andres Figliozzi proposed a promising algorithm that can well solve TDVRP with a soft or hard time window without any alteration [29].
To address TDVRP in the real world, this study aims to minimize the travelling time and path length of the route by solving the optimal delivery sequence and the best path at the same time. A novel fissile ripple spreading algorithm (FRSA) was proposed to find the optimal solution by avoiding daily traffic congestion in advance, or waiting to pass through a rapidly shrinking congestion. e coevolutionary path optimization (CEPO) technology was introduced to apply in the proposed method, according to the significant characteristic of traffic operational rules and delivery demand. e effectiveness of the proposed approach was illustrated in two sets of simulation experiments. e results proved that the proposed FRSA can effectively reduce the travelling time and path length of the vehicle route, which can decline fuel consumption and avoid regular traffic congestion. e rest of this study is organized as follows: Section 2 explains the problem description and the details of CEPO technology introduced in this work. Section 3 presents the basic idea of the proposed FRSA and shows how to integrate CEPO into it. Section 4 shows the experimental results by comparing the performance of FRSA with another compound method. Section 5 discusses some possible reasons for the outstanding results derived from FRSA and indicates its potential application. Section 6 outlines the conclusions of this work.

Coevolutionary Path Optimization
2.1. Problem Description. As a critical solution to reduce travel time and fuel consumption, TDVRP has been studied in extensive fields, which provides a great potential alternative to boost delivery efficiency and decrease economic cost. However, most of the existing studies ignored the practical network topology and separate the essential relationship between the shortest path of two sites and delivery sequence. Generally speaking, the delivery sequence decision is greatly affected by the weight values between two sites, which are obtained from the shortest path optimization. On the other hand, if a vehicle set out from a site at different periods, it tends to encounter different traffic conditions, and the shortest paths are likely to be changed. One can see that the path optimization results are also influenced by the delivery sequence. Furthermore, few studies focus on how to design a general algorithm to consider traffic congestion. Traffic congestion is the most common phenomenon in urban roads, which tends to cause a lot of time delay for vehicles as well as fuel consumption, along with some potential damage for emergency rescue. erefore, although some encouraging progress has been made, developing a better intelligent algorithm with excellent optimization performance and good applicability is less successful and still requires further research. To overcome these research gaps, this paper firstly proposed a universal approach called FRSA to address TDVRP, considering dynamic traffic congestion and the coupling relationship between the shortest path and delivery sequence.
anks to the foundation offered by the previous research, the proposed TDVRP in this study aims to make the actual delivery trajectory optimal, and it can be separated into two important issues. e one is to propose an appropriate path optimization method to find out the shortest path between two sites. e other is to exploit a feasible algorithm to determine the best delivery sequence with small travelling time and short path length on account of the path optimization method. In other words, the goal is to minimize the total travel time and path length by optimizing the vehicle delivery order and the delivery paths between any two sites.
According to the formal description of TDVRP in previous studies, several crucial rules should be established. A series of predetermined sites must be reached by a specific vehicle that is located in a single depot. e locations and the demands of each site are known, while the location and capacity of the related depot are given as well. Due to the complicated effects of traffic flow and road geometry, vehicle speeds may vary over time in a road section [30][31][32][33]. Additionally, when the vehicle travels on different roads or different driving directions, it tends to encounter various speeds and traffic conditions. To clarify the application scope of this study, two assumptions must be made. (1) e vehicle is allowed to wait at signalized intersections or congested links, and the delay caused by this inevitable procedure will be contained in total travel time. (2) e vehicle can get access to future traffic data from the traffic information prediction platform or big data prediction algorithm. e future traffic data include congestion scope, the average delay of signalized intersections, and average traffic speeds of all roads. (3) Any change in the routing environment will eventually be reflected by the certain variation of traffic speeds in related roads.

Basic Idea of CEPO.
In this paper, an advanced framework of CEPO is introduced to improve the path optimization algorithm, which concentrates on how to find out the best path between different sites. As a newly emerging approach, CEPO is firstly proposed by Hu et al. who have certified its outstanding performance and feasibility in discrete urban networks [34]. Moreover, another advanced variant timing coevolutionary path optimization (TCEPO) was developed to prove the effectiveness of CEPO by several comparative simulation experiments [35]. Due to its coevolutionary optimization mechanism, it can achieve a global optimal result as long as the given prediction data are precise enough. e prediction data in this study are obtained from those intelligent prediction algorithms, which are used to predict the short-term traffic flow or traffic speed based on historical traffic data. e prediction algorithms tend to adopt machine learning theory to predict traffic dynamics and usually achieve more than 90% of mean prediction accuracy [36][37][38]. Some of them can achieve higher prediction accuracy, even more than 95% [38]. ese encouraging results can provide the predictive traffic data, which are assumption basis for the proposed method. To explain how it can implement a superior optimization process, the basic idea of CEPO is illustrated in Figure 1. e definition of CEPO is, under the predicted traffic environment dynamics, each time-unit-oriented optimization step coevolves with the traffic speed variation of associated links once the traffic environment changes in routing networks, to make actual travelling trajectory optimal in a single run. All optimization steps will be carried out before the vehicle sets out from the depot. Unlike other mainstream research studies based on online reoptimization (OLRO), CEPO has more advantages in utilizing future traffic information and optimization accuracy. is is because the future routing environment dynamics are considered into the optimization procedure, which brings an extra superiority to know the evolutionary trend of routing environment in advance. As a result, it can help a vehicle to avoid wrong detour decisions or escape from congestion areas in advance. No matter which kind of roles it plays, CEPO poses a positive effect on the best route decisions that can help vehicle travel across the temporary congestion areas as fast as possible.

Mathematic Expression of CEPO.
It is assumed that the predictable routing environment dynamics and vehicle delivery sequence are both given, CEPO aims to address the minimization problem in terms of time cost. Travelling time is usually a major concern in path optimization, particularly in the dynamic urban routing environment [39,40]. Moreover, as explained in Section 2.3, CEPO technology is a time-unitoriented optimization methodology. us, it is an effective way to use travelling time as the optimization objective of CEPO. e objective function can be established as follows: which is subject to the following equations: where F c is the function calculating the time cost of all optimization paths. n denotes the number of delivery sites. P k denotes the kth optimization path between two sites in the presupposed delivery sequence. L(P k ) denotes the number of all nodes in path P k . f k (P k , L(P k )) denotes the time to travel through all nodes in path P k . f k (P k , P k (i + 1)) denotes the time to pass through the node P k (i + 1) when the vehicle just reaches the node P k (i). D t (P k (i)) denotes the time delay to cross over the node P k (i) in path P k at time t. C t (P k (i), P k (i + 1)) denotes the travelling time associated with the link between nodes P k (i) and P k (i + 1) at time t.
In this section, the objective function in (1) is defined as minimizing the total travelling time of the complete delivery loop. Equations (2) and (3) clearly show the calculation derivation of the travelling time, consisting of the time costs that pass through the nodes and links, respectively. It is believed that the time cost of a subpath in P k tends to be dependent on when the subpath is to be travelled through, which is often overlooked in other studies.
us, if the routing environment dynamics are predictable, one-off optimization will be prospectively conducted to seek the best path P * in the time-varying traffic environment before the vehicle starts to perform the delivery assignments. is means that the reoptimization of the path P * required in OLRO-based methods is unnecessary to use at each time unit.

Algorithm Origin and Development.
To realize the CEPO method as mentioned in Section 2, it is essential to exploit an appropriate algorithm to adapt to the synchronous state changes of nodes and links. It has achieved a great success for deterministic and heuristic methods in solving static path optimization and time-dependent path optimization. For example, the Dijkstra algorithm, A * algorithm, genetic algorithm, particle swarm optimization algorithm, simulated annealing algorithm, and ant colony algorithm. It is unlikely to apply them to an effective realization of CEPO. One important reason is that their basic computational step is the link-oriented analysis that associates with a long period, which may cause serious scalability problems [30]. What makes things worse is that it is impossible to pre-know which link will be accessed or when the given links will be accessed. Additionally, it must be considered for the improvement of delivery strategy and the best-first search of the shortest paths between sites. According to the latest knowledge, it seems to be quite hard to synchronously perform this procedure for most methods. is is why in a time-dependent routing environment, an advanced optimization algorithm needs to be designed to achieve the CEPO, to cover all necessary changes in the route network. FRSA is an evolutionary branch of the ripple spreading algorithm (RSA), which is first reported by Hu and Liao in 2016 [41]. Like other evolutionary computation algorithms, the basic idea of FRSA takes inspiration from nature and intends to imitate the ripple spreading phenomenon. e core technology of FRSA presents excellent optimality and surprisingly robustness in resolving complex problems, which is conducted by a multilevel and fissionable ripple relay race. FRSA is an agent-based, deterministic, and discrete simulation model, which provides a rare opportunity to evolve with the changes in the predicted environment at each time unit. During a single run of FRSA, it delivers a great number of ripples and the spreading speeds of ripples reflected in the changes of links to program optimization steps. Furthermore, when a ripple activates a site from the last site, the fissile mechanism is permitted to carry out more independent ripple relay races that explore optimal resolution in global research. us, once a single run of FRSA has been done, optimal delivery sequence and the shortest paths of given OD pairs (site to site or depot to site) can be obtained.

Fissile Ripple Relay Race Based on Global Research.
To realize the CEPO methodology, the optimization principle of the FRSA is extended to be divided into two parts, that is, independent ripple relay race and fissile program. An original ripple relay race is applied in a static route network [42]. But in the actual application, especially for the urban network, traffic environment dynamics should be integrated into an independent ripple relay. e illustration of an independent ripple relay race is given in Figure 2.
In the ripple relay race, the initial ripple is activated to split into several ripples from the origin node, and new ripples travel toward the adjacent nodes along with the connected links. When a ripple reaches an unvisited node, it will trigger new ripples at that node, and the unvisited state will be transformed into an active state. ose new ripples will continue to spread out until arriving in the adjacent nodes, while their spreading speeds are consistent with the traffic speeds of the link. As the ripple relay race goes on, numerous ripples will be generated to compete with each other until a ripple reaches the given destination node. Eventually, the shortest path from the origin to the destination will be determined by tracking back the visiting node order, when the ripple reaches the destination firstly. erefore, the travelling time of the shortest path between sites can be calculated as the end of an independent ripple relay race.
Although the independent relay race just offers a supportive foundation for addressing TDVRP in a dynamic routing environment, it is still less successful. It seems to be an ongoing challenge that how to decide an optimal delivery sequence with the smallest travel time cost. Hence, how to apply ripple relay race to resolve the best visiting order for all preset sites may be a question, which is worth thoroughly pondering. However, thanks to the development of computer science technologies, stronger computational ability can be used to execute big data processing. It opens a door to the fissile program that can evolve more independent relay races, to search more loop paths associated with different visiting orders.
To explore the optimal delivery sequence, it is a critical way to search for all delivery schemes in parallel by performing extensive ripple relay races, which are independent of each other. e diagram of the fissile program in FRSA is given in Figure 3. When the nearest site will be triggered by a ripple from the depot, it will generate a unique list recording activation order of the site. At the same time, a new ripple relay race starts to run at the site. All ripples from the new ripple relay race will not be affected by any ripple from the previous one. In other words, the ripples of different relay races are independent of each other, although they may arrive in the same nodes or links simultaneously. All ripple relay races run in the original size network. What is more important, the previous ripple relay race will continue to program, even if it has given rise to a new ripple relay race at the unvisited site. e other ripples from the previous relay race still own similar ability to trigger other sites, as long as these sites are not contained in the list of activity orders  Journal of Advanced Transportation 5 recorded by the associated ripples. is strict rule of FRSA indicates that a site will not be triggered twice in the same ripple relay race but remains accessible to others. As time goes on, more and more independent ripple relay races are conducted to evolve with the routing environment dynamics. e optimization of FRSA will be terminated, when all sites are recorded in the list of activation order and the related ripple finally reaches the depot in the first place. e optimal loop path can be gained by a single run of FRSA after a large number of computations for time-unit-oriented analysis.

Realization of FRSA considering the CEPO Technology.
To make the actual loop path optimal, it is crucial to permit ripples to wait at a node that represents a signalized interaction in practical application. e waiting behaviour of ripples reflects that the vehicles have to obey legal traffic rules, such as the restriction of red light and forward vehicles. Another condition is that the ripples can wait at a link or a node when they are temporarily blocked due to the evolution of the congestion area. However, this waiting behaviour provides a potential likelihood of waiting at a link or node until they become accessible again, which may take less time than travelling around the block area. It can impose a positive impact on avoiding the wrong detour, which is helpful to find out the shortest path with the smallest travelling time. erefore, to extend the FRSA to realize the CEPO technology, it needs to introduce a waiting behaviour of ripple into the relay race, to imitate the real performing actions of vehicles in the urban roads.
To understand the outline of FRSA, some variables need to be calibrated to explain how it works. Let N denote the total number of nodes in the network; n denotes the number of sites. For the sake of simplicity but without losing generality, this paper supposes that there are one hundred nodes, three sites, and one depot in the network. Let node 1 represent depot. m is the number of adjacent nodes for any node and sets m � 4, which means each node is connected to four nodes. In a ripple relay race, node i can be triggered with unlimited times and its ripple is called ripple i, which can split into several arcs spreading to adjacent nodes. S r (i) is the state of the ripple i, S r (i) � 1, 2, 3 means ripple i remains inactive, active, and waiting, respectively. r(i, j) is the radius of ripple arcs from node i to node j (i ≠ j). V t (i, j) denotes the spreading speed of ripple arc from node i to node j at time t, which is related to expand the associated arcs radius r(i, j). F r (i) denotes a crucial variable that records which ripple triggered the node i. For instance, F r (i) � j means node i is triggered by the ripple from node j, which offers detailed information to track back the travelling trajectory. When the ripple relay races are terminated, the optimal loop path P * will be determined by iterative deduction of F r (i), and the expression is illustrated in the following equation: To avoid the repeated activation of the sites in the same ripple relay race, M r (i) is used to record the existing travelling trajectory of the ripple i via tracking back F r (i). One can see that any site will not be repeated on each M r (i) except for the depot, which is the origin as well as the destination. In the real world, the spreading speed of ripple i is subject to the speed limit of the urban roads. us, V t (i, j) must satisfy the following constraints, which are given in the following equations: where V lim (i, j) is the maximal speed limit of the link (i, j). t unit is the time unit that is small enough to capture the major dynamics of the routing environment. L(i, j) is the length of the link (i, j). As time goes on, the arc radius of the active ripple will be updated according to the associated spreading. e calculation formula of the arc radius is shown in the following equation: If the node i is reached by the ripple j before the ending of a time unit, the new ripple i would continue to spread out in the remaining time t re . e spreading radius of arcs from the ripple i to its adjacent node k in the remaining time t re can be obtained according to the following equations: where t re is the remaining time in a time unit when a ripple arrives at a node. r(i, k) denotes the spreading radius of arcs from the ripple i to its adjacent node k in the remaining time t re . e flowchart of the FRSA based on the CEPO methodology is demonstrated in Figure 4. One can see that Step 1 to Step 3 well describe how to construct a simulation network and initialize some parameters as the prior preparation.
Step 4 defines the termination condition of the ripple relay race; that is, all sites are visited by a specific ripple that finally returns to the depot.
Step 5 indicates the function that converts the state of the node from waiting into active if the node is accessible.
Step 6 gives the derivation formula of the arc radius for an active ripple. Step 7 means that the remaining time will be used to develop the arc radius of the new active ripple, which is triggered by other ripples before the ending of a time unit. Step 8 demonstrates the requirement of a waiting behaviour at a triggered node, while it is blocked.
Step 9 updates the state of the nodes to inactive again, when their all ripples have reached adjacent nodes.
Step 10 judges whether FRSA needs to start a new ripple relay race from a node. Step 11 explains the mandatory operation rules for a new ripple relay race.
Step 12 outputs the optimal travelling trajectory based on complete M r (1).

Static Delivery Sequence Determination Tests Compared with GA.
In this experiment, to achieve a fair understanding of the performance of FRSA, four tests with different site numbers are conducted to resolve the traditional VRP. A widely used GA is selected to output the results for the comparative analysis. A series of sites (red points) and depot (green point) are generated as follows: the location of the depot is set as the left side as an example and the locations of the sites are randomly distributed in the right district. All nodes have the associated links and connected to their adjacent nodes, and the weights of the links depend on their length. To show the generality of the proposed FRSA in the delivery application scenario, the depot is given only one, while there are four sites scales with n � 10, 12, 14, and 16, respectively. e diagrams of the node connection for all tests are shown in Figure 5.
Since the ripple relay race of FRSA is the time-unitoriented process, it will be terminated once a ripple firstly finds the shortest path in a certain time unit. After that, the FRSA will arrange all the complete routes in that time unit, according to the path length from small to large to select the optimal path. For FRSA, an appropriate spreading speed can Step 1: create N nodes network; give the location and number of depots and all sites. Set node 1 as depot.
Step 4: if all sites are recorded in M r (1)?
Step 5: convert the state of the nodes. For i = 1: N, whether the node i is accessible, S r (i) = 3 and F r (i) > 0. If so, then change the state of node i from waiting to active. Update S r (i) = 2.
Step 6: calculate the arc radii of the ripples. For i = 1: N, whether S r (i) = 2. If so, then update all arc radii of the ripple i over time according to formula (7). Let t = t + 1.
Step 7: manage the remaining time. For i = 1: N, whether S r (i) = 1 and node i is triggered by adjacent node j. If so, then update all arc radii of the ripple i in the remaining time according to formulas (8) and (9). Let F r (i) = j and S r (i) = 2.
Step 8: realize the waiting behavior. For i = 1: N, whether the node i is triggered by its adjacent ripple j but blocked, and S r (i) = 0. If so, then let F r (i) = j, r(i) = 0, and S r (i) = 3, where r(i) denotes all arc radii of ripple i.
Step 9: update the state of the ripple. For i = 1: N, whether the ripple i has triggered each adjacent node and S r (i) = 2. If so, then let S r (i) = 1.
Step 10: if any site is reached by the ripple?
Yes Yes No

No
Step 11: generate new independent ripple relay race. For any ripple, i will not trigger the sites that are contained in M r (i) and the depot except for all sites contained in M r (i) but still can trigger other normal nodes. Assuming that rhe site k is reached by the ripple j, then let F r (k) = j, S r (k) = 2, M r (k) = {l, ... , j, k}.
Step 12: output the optimal loop path P * , which is recorded in M r (1). Journal of Advanced Transportation improve the computational efficiency and reduce computation or storage pressures. To ensure the computation efficiency of FRSA optimal, the spreading speed of ripples cannot be too large or too small in the network. If the spreading speed is relatively high, the path length travelled by ripples will be pretty long in a time unit. When the optimal route is solved in a time unit, a great number of travelling trajectories that are closed to the optimal path are likely to be found in the remaining time and still need to be arranged. e unexpected results pose extra work on the calculation module, which means that the computational time required for sorting path length of each complete route increases accordingly. If the spreading speed is relatively low, resulting in more time units required to extend the arc radius of the ripples. e workload of updating the ripple state significantly increases with the growth of time units used in a single run. Hence, the ripple spreading speed can neither be too large nor too small. After comparing multiple spreading speeds through simulated tests, this paper proves that selecting 60 km/h has better computing efficiency. Additionally, this speed also conforms to the speed limit management regulations of urban roads in China. As a result, this experiment sets ripple spreading speed v(i, j) � 60 km/h. e related results of FRSA and GA are given in Table 1 and Figure 6. One can see that, in terms of PL, it is apparent to observe that FRSA still outperforms GA in all kinds of nodes sizes. e distance of the optimized paths given by FRSA is approximately 33.5% to 57.7% smaller than that of GA. Compared with GA, the reduction rate of FRSA reaches the minimal value at n � 12 and peaks at n � 14 in terms of PL. One possible reason is that the core technology of FRSA is to globally explore all paths, and eventually ended when the ripple has visited all sites without repetitive activation and returns to the origin. Compared with GA, it possesses better search capacity and will not fall into a local optimum. Unfortunately, as a sacrifice of the huge advantage in solving the best path, the CT of FRSA remarkably goes up with the expansion of node sizes.  Reduction rate (%) The reduction rate of FRSA compared with GA in terms of PL  However, unlike other global searching algorithms, FRSA needs no requirements to traversal all results. is is because the optimization steps are divided into a lot of time units, and the ripple relay race will be terminated once the first group of paths is found at a certain time unit. It means that extensive data calculation used for searching other longer loop paths has been omitted. In this respect, FRSA has achieved great progress in terms of CT without losing global searching capacity, although it is relatively large compared with those heuristic algorithms.

Delivery Trajectory Optimization
Tests considering Dynamic Routing Environment. As mentioned in Section 2.1, this study is especially interested in addressing a significant issue via designing an efficient methodology based on CEPO technology, which aims to make the actual delivery trajectory optimal. What makes things difficult is that the simultaneous resolution of the shortest paths between sites and the optimal delivery sequence poses a great challenge on the current methods and tends to be ignored. erefore, this study investigates the interaction relationship of them and tests the performance of FRSA based on CEPO, comparing to the mainstream heuristic algorithm based on OLRO. Different from Section 4.1, traffic dynamics of networks are introduced to imitate the actual traffic environment. However, travelling time is significantly affected by traffic flow in urban roads. When the traffic flow of the specific road section is large, the associated vehicle density also is high, which means the space headway between vehicles will be small. However, the smaller space headway between vehicles, along with the more serious disturbance effect between vehicles. is is because the drivers are subject to the front vehicles or the lane-changing behaviour of the adjacent vehicles. e frequency and severity of traffic flow turbulence will be significantly enlarged as the growth of traffic flow, such as dangerous deceleration behaviour, risky lane changes, and sudden accidents. us, the negative effects of traffic flow turbulence lead to the reduction of average traffic speed, which poses a direct impact on the increase of travelling time. Conversely, if the traffic flow on the road section is small, the speed of the vehicle tends to be less affected by that, and the travelling time is also reduced accordingly.
us, travelling time also needs to be an additional assessment indicator, which is used to analyze the temporal efficiency of the optimized routes. In this study, a set of experiments is conducted to verify the promising application of the CEPO-based FRSA in terms of travelling time (TT), path length (PL), and computational time (CT).
To obtain the shortest path between sites in GA, a popular Dijkstra algorithm is introduced to achieve the realization of the OLRO [43][44][45][46]. In OLRO, the basic application of the Dijkstra algorithm is to repeat calculating the current shortest path based on the dynamic routing network in each time unit until the vehicle reaches the destination. Based on the current shortest paths between sites, GA is used to find out the optimal delivery sequence via iterative computation and finally presents the complete delivery path. In CEPO, the proposed FRSA will give the optimal delivery trajectory in a single run according to the future routing environment.
In this part, different sizes of urban networks are randomly generated as following: the locations of N nodes are randomly evenly distributed in a rectangular area. ere are 4 network scales with N � 64, 144, 256, and 400, respectively. All nodes are connected to their neighbouring nodes, while each node has 3 or 4 connections. e shortest link is less than 150 m. e depot (green point) is always set as the central location of the network, while 8 sites (red points) are distributed around the depot randomly. Furthermore, to achieve the best computational efficiency, 1 min is selected to represent a time unit. However, the most important thing is to construct the routing environment dynamics, which can be described as traffic speed dynamics: where σ t (i, j) denotes the traffic speed variation of the link (i, j) at time t. ere are three hypotheses for processing nodes and links, as follows: (i) To restore the actual traffic operation in the urban network, it is needed to determine the maximum speed and speed variation of the links at each time unit. Suppose 0 < V t (i, j) ≤ 80 km/h and the absolute value of speed variation is less than 4 km/h for all links (|σ t (i, j)| ≤ 4) at each time unit. (ii) Due to the compulsive traffic laws in the signalized intersections, the vehicles must wait at the nodes until they get access to the right to pass through. us, the stopping delay caused by waiting at nodes should be contained in travelling time during path optimization. Let d(i) denote the average delay of node i, and its value randomly varies from 20 to 60 seconds. (iii) As a common phenomenon in the urban network, the evolution of the congestion area also should be considered in TDVRP. As shown in Figure 7, the central nodes of congestion areas are evenly distributed in the network and have three kinds of states: spreading, static, and shrinking. It is supposed that the spreading or shrinking congestion areas will expand or narrow in every two minutes, and this action can evolve three times at most. Once the congestion areas spread to the maximum scope, it will keep unchanged in two minutes and then begin to shrink in a regular period.
Some experimental results are clearly given in Table 2 and Figure 8. Table 2 shows the results of the actual delivery trajectories given by FRSA and GA, in terms of TT, PL, and CT. One can see that no matter which scale of the route networks, the proposed FRSA (CEPO-based) outperforms GA (OLRO-based) in terms of TTand PL. From Figure 8, the TT of FRSA is approximately 11.86% to 20.25% smaller than that of GA. Moreover, the PL of FRSA is approximately 0.86% to 27.16% shorter than that of GA. It is obvious to prove that the large PL is often inconsistent with the larger TT [34]. In small-size networks (N � 64, 144, and 256), the CT of FRSA is about 35.52% to 74.37% smaller than that of GA. Nevertheless, this trend tends to reverses when the network size becomes relatively large, and the CT of FRSA is 13.22% longer than that of GA when N � 400.
In terms of TT, PL, and CT, the reduction rate of FRSA gradually boosts from n � 64, then reaches the peak at n � 144, and finally declines smoothly over the boost of nodes. Compared with GA, surprisingly, the reduction percentage of FRSA is no-linear growth as the expansion of network size in terms of TT, PL, and CT.

Discussion
e main goal of this study was to validate the extensive applicability of FRSA based on the framework of CEPO. With this aim, this study performed traditional tests for FRSA and GA in the static networks, as presented in Section 4.1. For the sake of fairness, the reduction rate was used to measure the difference between FRSA and GA in terms of the TT, PL, and CT. e proposed FRSA exhibits superior searching ability and is far better than GA in terms of the cost of the related solutions. However, the encouraging findings also bring corresponding disadvantages. As regards the most serious defect, a massive amount of data is handled in FRSA, along with a sharp increase of CT as the number of sites goes on. An inspiring discovery thing is that, compared to other accurate approaches that achieved a similar effect, FRSA has cut down a tremendous number of computational works without seeking other possible solutions.
As discussed in the previous section, this study aims to propose a commonly used algorithm to address various popular branches of VRP, especially for TDVRP. erefore, the time-varying traffic environment needs to be considered in testing FRSA as an inevitable development. As a crucial indicator, the TT of the solutions becomes a current focus. In a dynamic routing environment, FRSA is always superior to GA in terms of TT, PL, and CT, except for a special case (N � 400) that the CT of FRSA is larger than that of GA.
ere are three main reasons giving rise to the results that FRSA outperforms the GA in terms of TT and PL. One reason for this is that, in FRSA, the optimal delivery sequence and the shortest path between given OD pairs are solved effectively at the same time. GA just focuses on how to ensure the delivery sequence better, while Dijkstra only concentrates on how to find the best path between OD pairs. e second reason is that FRSA has more advantages to derive the optimal solution via global search, while many heuristic algorithms can only find an approximate solution through parallel computing. e last reason contributes to most to this result, the time-unit-oriented optimization steps of FRSA can coevolve with the changes of routing   environment and achieve the realization of CEPO. It owns unique features that avoid the coming congestion in advance, or allow waiting behaviour to prevent the wrong detour measure. erefore, under the coupling effect of three critical factors, the actual travelling trajectory with a shorter TT can be obtained in the nature of things. Except for the big size network like n � 400, the CT of FRSA is smaller than that of GA in addressing TDVRP. As a result, a smaller CT is spent in FRSA to solve the path optimization in dynamic routing network. But, in n � 400 network, there are no more advantages for FRSA compared with GA. e main reason contributing to this is that each node provides an opportunity to trigger the ripples, and one node is usually triggered many times during optimization. More nodes also imply the boosting number of ripples will be generated, which asks for a larger CT to update the state transition and ripple spreading situation in a time unit.
However, FRSA possesses some distinct advantages in dealing with VRP, compared with the previous methods. e most outstanding one is that it shows a reliable guarantee of optimality in a time-varying routing environment. e other promising features, like a one-off run before a complicated trip begins and the positive avoidance of congestion areas, also proved that it is superior to the other. ese results suggest that, in some urban networks (less than 400 nodes), it would be advisable to execute the FRSA to search the best solution as a practical implementation. As for large network scenarios, it needs to pay attention to whether there is a strict requirement for computational time. If not, FRSA will be a good choice to make the actual travelling trajectory better for addressing the existing VRP.

Conclusions
is work paves a new way to address various VRP variants, particularly for TDVRP whose main features have been studied for several decades. A mathematical model of FRSA is formulated based on the CEPO technology, which makes use of the future traffic environment to implement path optimization. Due to the effects of traffic flow, network structure, and driver characteristics, some critical factors like timedependent traffic speeds, network topology, temporal congestion area evolution, different vehicle departure times, vehicle waiting behaviour, and network size are considered in the proposed model [47][48][49]. e objective of the proposed method is to minimize the sum of the delivery costs, such as travelling time and path length. Unlike the existing research studies, this study pays extra attention to explain the internal correlation between determining the delivery sequence and searching the shortest paths between sites, which tends to be ignored. However, FRSA owns several outstanding features, including a congestion avoidance approach, a single run optimization, waiting measure, parallel computing, coevolutionary planning, and global research. e effectiveness of the proposed approach was demonstrated in two sets of simulation experiments. e experimental results show the following: (1) FRSA presents an excellent improvement in reducing actual travelling time and path length, effectively avoiding the wrong detour and regular traffic congestions. It can help to decrease the total distribution cost, fuel consumption, and carbon emissions, which can effectively protect the environment. (2) e optimal delivery sequence and the shortest paths can be derived simultaneously in the dynamic routing networks. is encouraging characteristic is significantly important and beneficial to make the actual travelling trajectory better, which is more likely to be adopted in the practical application compared with other methods. (3) e application of CEPO methodology considering the future routing environment dynamics will decline the travelling time. By contrast, other OLRO-based methods tend to take wrong coping strategies to avoid traffic congestions, which will increase the time and economic costs for logistics companies.
Despite some inspiring results that have been achieved, there still exist some limitations that need to be overcome. Further research may consider the extension of the proposed algorithmic framework to address other VRP variants in special scenarios, by taking into account time window, service time, vehicle capacity, site demand, and pickup. Additionally, the proposed FRSA can be further developed to resolve scalability problems and improve computational efficiency. As a new and effective deterministic algorithm, FRSA opens the door to the development of precise exploration.

Data Availability
e data used to support the findings of this study are currently under embargo while the research findings are commercialized. Requests for data, 12 months after publication of this article, will be considered by the corresponding author.