A solution approach to the orienteering problem with time windows and synchronisation constraints.

The orienteering problem with time windows and synchronisation constraints, known as the Cooperative Orienteering Problem with Time Windows (COPTW), is a class of problems with some important applications such as in home health care and emergency logistics management, and yet has received relatively little attention. In the COPTW, a certain number of team members are required to collect the associated reward from each node simultaneously and cooperatively. This requirement to have one or more team members simultaneously available at a vertex to collect the reward poses a challenging task. It means that while multiple paths need to be determined as in the team orienteering problem with time-windows (TOPTW), there is the additional requirement that certain paths must meet at some of the vertices. Exact methods are too slow for operational purposes and they are not able to handle large scale instances of the COPTW. In this paper, we address the problem of finding solutions to the COPTW in times that make the approach suitable for use in certain emergency response situations. This is achieved by developing new merit-based heuristics as elements of an Adaptive Large Neighbourhood Search (ALNS) algorithm. We validate the performance of this new approach through an extensive computational study. The computational results show that the proposed method is effective in obtaining high quality solutions in times that are suitable for operational purposes.

Introduction

The Orienteering Problem (OP) is a well-known problem in combinatorial optimisation, introduced by Golden et al. (1987). The OP emerges from a combination of the travelling salesman problem and the knapsack problem. It is a routing problem where travelling to all vertices is often not feasible due to a time constraint. The objective of the OP is to find the combination of nodes that maximise the total reward collected. The design of tourist itineraries is an example where the OP has been applied (Vansteenwegen and Van Oudheusden (2007), Vansteenwegen et al. (2011a), De Falco et al. (2015)). The tourist trip design problem seeks to select the most interesting combination of attractions to visit within some available time span (Vansteenwegen and Van Oudheusden (2007), Borràs et al. (2014)). An extension of this problem to cycle trip planning was proposed by Verbeeck et al. (2014). Further applications can be found in the review by Vansteenwegen et al. (2011b). There are several variants of the OP Gunawan et al. (2016). For example, the time-dependent orienteering problem considers the case where the travel time between two vertices depends on the leaving time of the first vertex Abbaspour and Samadzadegan (2011). This situation arises, amongst others, when travel times lengthen during peak traffic times. The time of arrival at a rail station, for example, can have a significant effect on waiting times when using public transport (Garcia et al. (2010)). The Team Orienteering Problem with Time Window (TOPTW) has received a great deal of attention (Labadie et al. (2012), Kim et al. (2013), Souffriau et al. (2013), Duque et al. (2015), Gunawan et al. (2015), Zhang et al. (2018)). The TOPTW arises when the sport of orienteering is played by teams of several people. The TOPTW extends the OP to identify multiple paths that maximise the total score while meeting the time window constraint at visited locations or finishing within a fixed time at the final destination.

The Orienteering problem is NP-hard (Golden et al. (1987)). An extensive number of exact and heuristic approaches have been proposed to solve this problem. A few papers have focused on exact methods (Keshtkaran et al. (2016), Poggi et al. (2010), Dang et al. (2013a), Bianchessi et al. (2018)) but generally they are not efficient for large-sized problems. As a result, the main body of the TOPTW literature is dominated by heuristic approaches (Bouly et al. (2010), Liang et al. (2013), Marinakis et al. (2015), Dang et al. (2013b), Gunawan et al. (2015), Labadie et al. (2012), Gunawan et al. (2017), Roozbeh et al. (2018b)). As the TOPTW has become a popular topic, several state-of-the-art papers on the solution methods can be found in the recent literature. Gambardella et al. (2012) enhanced the Ant Colony algorithm to overcome certain drawbacks of the algorithm that they developed earlier (Montemanni and Gambardella (2009)). Vansteenwegen et al. (2009) developed an iterated local search algorithm for the TOPTW, where high quality solutions are attained in much less computational time than exact methods. More complicated benchmark instances for the multi-period orienteering problem, which is a generalisation of the TOPTW, are solved by a variable neighbourhood search algorithm (Tricoire et al. (2010)). Lin and Vincent (2012) developed the slow and fast simulated annealing algorithms, where the latter performs efficiently in terms of computational effort and the former is more concerned about the solution quality. Gavalas et al. (2013) focused on tourist trip design problems. They introduced two cluster-based algorithms by considering the limitations of the iterative local search algorithm proposed by Vansteenwegen et al. (2009). Gunawan et al. (2015) extended the existing iterative local search algorithm by including more local search operations. Vidal et al. (2015) used variable neighbourhood search with an efficient select algorithm, equipped with iterated local search and a hybrid genetic algorithm, to solve the team orienteering problem. Also, Mei et al. (2016) investigated the multi-objective, time-dependent orienteering problem using a memetic algorithm. Finally, a recent article by Santini (2018) revealed the significance of using an ALNS algorithm for solving OPs.

An extension to TOPTW is the Cooperative Orienteering Problem with Time Windows (COPTW), proposed by Van Der Merwe et al. (2014). This problem arises when some tasks that need to be undertaken at certain locations can only be accomplished by two or more individuals acting cooperatively and simultaneously. Thus, in the COPTW, each vertex has a unique resource requirement. This requirement specifies the number of team members that simultaneously must visit each vertex to collect the reward at that vertex. The reward is also conditional upon the requirement being met within the time windows specified for each vertex. This means that while multiple paths need to be determined as in the team orienteering problem with time-windows (TOPTW), there is the additional requirement that certain paths must meet at some of the nodes.

The COPTW arises in several applications which motivate this study. The first application arose from a problem involving asset protection during escaped wildfires (Roozbeh et al. (2018a)). Valuable assets such as bridges and hospitals are distributed over a landscape with a network of roads. A wildfire sweeps over this landscape in certain directions. For a wildfire that is beyond control, fire-fighters are deployed to visit each asset to undertake tasks that will mitigate the risk of each asset's destruction. The asset must be serviced before the fire reaches it, but must also not be serviced too early as hosed down structures will dry out again or areas cleared of debris can become littered again. Some assets might require a simultaneous visit by an aerial truck for accessing tall structures and a pumper. Others might need a tanker (own water) and a personnel vehicle.

Tour planning for UAV's is another application of the OP. In this case, rewards at a node correspond to priorities of target locations. Some of these problems require coordination and cooperation typical of the COPTW (Bertuccelli et al., 2004). Various OP applications of UAV's include operations management in the event of disaster, and for other military and civilian purposes. Issues arising from uncertainties in these problems are also addressed (Bertuccelli et al., 2004, Evers et al., 2015). To handle this complexity, the OP has been extended to include, for example, the OP with Stochastic Weights (Evers et al., 2014c), Robust OP (Evers et al., 2014b) and Online Stochastic OP with time windows and time-sensitive targets (Evers et al., 2014a).

A further application is in the home care services. In such problems, some operations may require more than one staff member. For example, where heavy lifting and specialist medical expertise are required (Bredström and Rönnqvist (2008)). A closely related problem is in the home cleaning services where some homes may require one or more services such as basic cleaning, window cleaning, and washing. Access to each house is within a given time-window. Further potential applications of the COPTW could be considered in post-disaster emergency services, where some services require visits from a supply vehicle to be synchronised with visits from personnel to distribute supplies.

For some applications of the COPTW, short solution times are essential. Due to the nature of the COPTW with its synchronisation constraints, the use of exact algorithms, as adopted by Van Der Merwe et al. (2014) is only suitable for relatively small instances. For large problems, exact methods become computationally prohibitive, especially where time-critical operational decisions have to be made. Thus, the aim of this paper is to develop an ALNS algorithm for the COPTW capable of handling large-scale problems and delivering good solutions within times suitable for operational purposes in emergency situations. To achieve this we build on the ideas given in Roozbeh et al. (2018b) and which were used to solve an asset protection problem (APP), a special case of the COPTW. Their method involved adapting a basic merit-based algorithm followed by various ruin and recreate heuristics. In this paper extensions to the proposed algorithm are introduced along with new insertion heuristics, and different problem-specific ruin algorithms are developed to efficiently search the problem space. These heuristics work together in an ALNS framework by detecting ineffective or redundant elements of the ALNS and deleting them. The performance of the approach was evaluated and tested for both small and large-size instances. For small instances, results are compared with optimal solutions obtained using exact methods. However, when it came to large scale instances CPLEX was unable to solve the problem to optimality. We also generate a large set of COPTW instances and report benchmark results for future studies.

In the following section the mathematical formulation of the COPTW problem is described. In section 3 we introduce the general scheme of the proposed approach. Our findings from an extensive experimental study are discussed in section 4. Finally, the conclusion and potential future research directions are stated in section 5.

Cooperative orienteering problem with time windows

The COPTW generalises the TOPTW formulation (Labadie et al., 2012), where a certain number of team members are required at vertices (Van Der Merwe et al., 2014). The location is considered served if a required team members arrive at the vertex within the time windows and start the service simultaneously at time for a duration of time units to collect a reward, . In the COPTW, a homogeneous fleet of P team members start their route from and must return to the depot by time , where both and represent the same location. Unlike the classical OP, the COPTW enforces time windows for service activities and therefore is represented as the closing time for the depot, . For any two vertices, indicates the required time for each team member to travel from i to j. The binary decision variable takes the value 1 if team members visit within the appropriate time windows for the required duration, 0 otherwise. Another decision variable is which takes 1 if is traversed, otherwise 0. Furthermore, represents the number of team members travelling from i to j. Lastly, ξ defines a set of arcs that can be traversed. For , if a team member departs from vertex i at and arrives at vertex j before the can be traversed. Following the definition of ξ, a set of feasible arcs that can be traversed to and from node i are shown by and , respectively.

Van Der Merwe et al. (2014) formulated a mixed integer programming model of the COPTW. For ease of reference, the model is repeated below with a few minor changes in notation.

The objective function (1) maximises the sum of the rewards collected at each vertex i. Constraint (2) ensures all members depart from and return to a designated depot. Constraint (3) guarantees flow conservation by enforcing the equality of incoming and outgoing arcs to each node. Constraint (4) ensures that the collection of a reward (i.e. a score) at each location is dependent upon the condition of fulfilling its resource requirement. Constraint (5) makes sure that the number of travelling members of a fleet through an arc never exceeds P and it also ensures that is zero when the path ij is not traversed. Constraint (6) ensures that at each vertex the service can only be started when the previously visited location has been served completely and there is sufficient time to travel to the vertex, where M represents a large constant. Setting is sufficiently large for this purpose. Constraints (7) and (8) ensure each vertex is visited within its time window. Integer and binary conditions are defined in constraints (9) and (10). Note that the , in the time budget constraint in orienteering problems is assumed to be the closing time of the ending point () meaning that there is a time window for all nodes and routes should be finished by returning back to the final depot by .

A graphical representation of a simple solution for COPTW is sketched in Fig. 1 as shown below. Two of the main attributes of each node, the resource requirement () and corresponding reward (), are specified while all nodes have a time window constraint. A given number of team members leave the depot and members must be at vertex i before starting the service to collect the associated reward (). The numbers over each arc represent the number of members travelling through an arc. It is infeasible to visit all vertices within their time windows. Considering the set of all feasible solutions, a schedule that maximises the score (i.e. the sum of rewards collected) is identified. In Fig. 1, three members leave the depot by travelling through different arcs and attend nodes 4, 1 and 3 individually. As their associated resource requirements are all equal to one, they complete the service and collect the rewards. Thereafter, all three members travel to node number 6 to fulfil its resource requirement and collect the associated reward (10). After collection of maximum amount of possible scores, all members return to the depot. Note that both depots represent the same location.

Figure 1

A sample solution of the COPTW.

A solution approach for the COPTW

We propose an Adaptive Large Neighbourhood Search (ALNS) that brings together a powerful algorithmic framework to solve the COPTW. We extend the work by Roozbeh et al. (2018b) developed for the asset protection problem. Each of the ruin and recreate heuristics is re-examined for the COPTW and for opportunities for further efficiency enhancements and then modified accordingly. In addition some new insertion and removal algorithms are introduced. Significant changes from the APP algorithm are indicated with an asterisk (*) in Table 2. In this section, we define the ALNS algorithm and describe the heuristics implemented.

Table 2

Statistics for performance of the removal and insertion algorithms. The problem specific heuristics are indicated with an asterisk (*). Numbers in the Table 2 represent the average results for five instances after ten runs (5 instances×10 runs). The second column denotes the name of algorithms, the third column illustrates the percentage of total iterations that a corresponding algorithm is chosen. The fourth, fifth and sixth columns show the change in the objective values after exclusion of the corresponding algorithm. A positive value indicates an improvement in the objective value. The last column indicates heuristics chosen to omit (Y) or not (N) in subsequent use of the ALNS.

ALNS procedure

The proposed ALNS algorithm comprises the construction of an initial solution followed by ruin and recreate strategies to improve the solution. The latter process continues until some termination conditions are satisfied. Stopping conditions for the ALNS heuristic are based on the total number of iterations and the run time limit. Through every iteration () a ruin heuristic () is selected and partially destroys the current feasible solution (). Thereafter, a recreate heuristic () inserts unvisited nodes into enough routes to achieve a better solution (). In the ALNS, a wide ranges of removal and insertion heuristics are implemented and being used dynamically and adaptively according to their past performance. This leads to the exploration of large neighbourhoods and escaping from local optimal solutions. An adaptive weight and a score is defined for each ruin and recreate heuristic (). Initially, all weights are equal, and scores are set to be zero. The algorithm runs for N iterations, divided into k segments and therefore iterates times over each segment. The corresponding score of ruin and recreate heuristics are increased by , and , defined below, after each iteration. If the new solution improves the current best solution (), the value of is added to the score of the associated ruin and recreate heuristics. The value of the is added to scores if the new solution improves the best current solution through a segment. Finally, the value of is added if the new solution is accepted despite degrading the objective value. After each segment, the weight of algorithm are updated by using the formula . The parameter ρ is called the reaction factor. It effectively regulates the number of iterations that ineffective algorithms are utilised when they are not contributing substantially to improving the solution. Also, associated weights, scores and number of times a heuristic is selected are defined by , and , respectively. The probability of using an algorithm in segment is defined based on the formula where the denominator calculates the cumulative weight of ruin and recreate heuristics and the scores are set back to zero when a segment has been completed (Roozbeh et al., 2018b).

It has been shown by Santini et al. (2018) that acceptance criterion can significantly impact the solution quality. For this reason, we take advantage of Simulated Annealing (SA) (Van Laarhoven and Aarts (1987)) to escape from local optimal solutions by accepting solutions with worse objective values with probability , while solutions with improved objective values are always accepted. In the SA, T, and represent temperature, the objective value of the current solution and the best solution, respectively. The initial temperature needs to be defined in such a way as to accept solutions with probability where the objective is δ% worse than the initial solution . Therefore, can be set as . The temperature is updated every N iterations according to the cooling rate, () where .

Insertion feasibility

A survey of vehicle routing problems with different types of synchronisation constraints can be found in Drexl (2012). Synchronisation constraints exist in various areas such as home health care routing and scheduling (Mankowska et al. (2014)), pick-up and delivery (Clarke et al. (2017)), transportation and logistics, medical delivery and airport ground handling services (Fink et al. (2019)). Our study focuses on synchronised operations in which multiple team members are required to visit a set of nodes to accomplish various tasks. Synchronisation constraints significantly increase complexities in maintaining feasibility during insertion operations. To efficiently handle these complexities in a time efficient manner, the following approaches are implemented.

We first filter out infeasible arcs through the search process; therefore, when for all nodes, then will be removed in the initial solution. Also, due to the requirement of nodes that need to be visited by multiple routes, we filter visiting of node i after j if j has already been inserted before i in one of the existing routes. We finally define a “feasibility matrix” comprising binary elements to facilitate the insertion process. Most importantly, to evaluate the possibility of insertions in constant computational time within a route, we need to define . This is to determine the existing amount of spare time that can be invested to meet a node i before arrival. Arrival time at a node can simply be defined as below. Since a service must be accomplished simultaneously by a required number of team members (T), a synchronised start of service is calculated as below. Subsequently the time that a team member waits at node i to start the service is as follows.

For a given route τ, a visit at node i is defined by . The equation (14) defines the time that arrival can be delayed () while the solution feasibility is satisfied. To calculate the we start from the last visit in each route. As synchronisation visits are required, the minimum value of for each node in existing routes must be taken. Therefore, for the if there exists such that we have:

To determine the insertion feasibility of a node k between i and in route τ, calculation of the caused shift () is required. Thus: If the caused shift () is less than or equal to the , the insertion will be considered as valid.

Insertion algorithms

We introduce new variants to the classical merit-based (MB) heuristics (Roozbeh et al., 2016), where their speed and performance assist us in handling complexities of the problem. The MB heuristic is utilised with different approaches including the classical MB heuristic, the noised imposed MB heuristic, the cumulative MB and the noised imposed cumulative MB heuristic. In addition to the variants of the MB heuristics, our algorithm, which can be regarded as a greedy time heuristic is formulated with consideration for the relative travel time. Overall, the proposed algorithm includes five insertion heuristics, which are described below.

Classical merit-based heuristic

The MB algorithm constructs a list of tuples . The refers to the merit value which is calculated by using the following function: In equation (17), is the angle between the vectors from the depot to the nodes i and j, respectively. The score for each node and the average score of all vertices are represented by and , respectively. Also, in the above formulation, is the distance between nodes i and j and is used to scale the distance term. To ensure the assignment of nodes with the highest merit values are prioritised, we implemented parallel route construction in the algorithm. Moreover, best values for the three parameters, ϑ, λ and μ, are defined in the parameter tuning part of the algorithm. Please note that travel time and travel distance are assumed to be correlated and they can be used interchangeably while having no impact on the order of merit values.

Algorithm 1 represents the MB heuristic. The algorithm starts by ranking unvisited nodes according to the ratio where are the resources (team members) required to collect the reward at node m. By experimentation the square root in the denominator was found to yield better results than a straight ratio of reward to required resources. Next, the algorithm creates a Euclidean distance matrix and a feasibility matrix as explained in section 3.2. The Merit Pair List (MPL) is constructed and sorted in descending order of (lines 5 and 6). The algorithm iterates over each node in the set U (line 7) for each pair (line 9) as it attempts to assign unvisited nodes into an existing subroute. Otherwise if resources are sufficient a new route is created (line 17). Relevant information is then updated (lines 22-25).

Algorithm 1

Pseudocode for the merit-based heuristic.

Noise-imposed merit-based heuristic

The noise-imposed merit-based heuristic (NMB) heuristic is used to bring additional diversification to the search. We define a noise parameter , then is the amount of noise allowed. The NMB heuristic updates merit values by where .

Cumulative merit-based heuristic

Through the insertion process, using the classical MB heuristic we look at the merit value of inserting the unvisited node i where node j is the succeeding node. However, we may need to consider also the node preceding i to accomplish an efficient insertion. Thus, to incorporate both the succeeding and preceding nodes into the merit values, we make a tuple of four components () where is the merit value that can be achieved by insertion of the unvisited node k between nodes i and j.

The Cumulative merit-based (CMB) heuristic uses the calculated merit values for pairs in equation (17) and determines values of by equation (18). Fig. 2 illustrates insertion operation before (Fig. 2(a)) and after (Fig. 2(b)) using equation (18). Insertion by equation (17) may result in a route as defined in Fig. 2(a). In the generated merit pair list using the equation (18), as both succeeding and preceding nodes are considered, the gets a larger value than . Therefore, insertion will take place as shown in Fig. 2(b) which is the optimal assignment.

Figure 2

An illustrative improvement by the cumulative merit-based heuristic. Fig. 2(a) is an insertion made by the MB heuristic compared to the CMB heuristic in Fig. 3(b).

Figure 3

An illustrative insertion of the greedy time heuristic. Fig. 3(a) shows two alternatives for insertion, and Fig. 3(b) is the insertion made by the GT heuristic.

Noise-imposed cumulative merit-based heuristic

The same approach as in section 3.3.2 is used for the noise-imposed cumulative merit-based (NCMB) heuristic to broaden the search space.

Greedy time heuristic

Time spent travelling may curtail the number of nodes that can be visited. Thus a Greedy Time (GT) heuristic is developed. The GT heuristic focuses on travel time minimisation. A list of tuples is generated and sorted in ascending order of time ratio values, , defined by equation (19). Note that, although this equation focuses on travel time, it includes other problem specific attributes in the denominator. In equation (19), the numerator calculates the increase in the travel time as a result of inserting node k between i and j. The is the associated score for node k which is divided by the square root of the ratio of the requirements at node k () to total number of team members, (P).

In Fig. 3(a) two potential nodes are shown for an insertion between nodes i and j. Although visiting node L involves a longer travel time than node k, the GT heuristic inserts L into the route. This is because impacts on the total value of and ends up with a higher ranking in the list.

Ruin methods

The initial solution generated by the classical MB heuristic () is partially destroyed by one of the ruin (removal) algorithms. At each iteration, a ruin heuristic removes D nodes from the solution, where D is called the degree of destruction. The value for D is a random number in the range of when K nodes are visited. In the following, we explain all ruin methods that are implemented in our ALNS in details.

A Random Removal (RR) algorithm randomly removes D customers from the solution. Worst-Distance Removal (WDR) determines the distance of nodes from succeeding and preceding nodes and sorts them in descending order in a list named as O. Then, it removes nodes in position from the solution and the list. Parameters and bring randomness to select new nodes for removal. Worst-Time Removal (WTR) performs the same as the WDR, however the list O is defined based on values of . Shaw Removal (SR) measures relatedness between nodes using where and are the Shaw parameters and varies in the range of based on the number of routes that nodes i and j have in common. The SR heuristic picks a random node and measures relatedness of other nodes with the one selected. Thereafter, position will define the node that has to be removed from the relatedness list . In the SR algorithm parameters is the determinism factor and is a random number. The SR increase the chance of substitution of nodes with the ones that have a high level of similarity in terms of distance, routes in common and time window, so that the search can move to a new neighbourhood. Proximity Removal (PR), Time Removal (TR) and Requirement-Based Removal (RBR) algorithms are special cases of the SR heuristic where , and are set to a value of 1, while others are set at 0.

Waiting-Time Oriented Removal (WTOR) is similar to the WDR, however the cost is the waiting time calculated as . Worst-Requirement Removal (WRR) algorithm removes nodes with highest demand from the solution. Relative-Requirement Removal (RRR1) heuristic is designed for heterogeneous fleet of members (i.e. vehicles). The RRR1 algorithm follows the same logic as the WRR, however, considers available number of members from each type. Therefore, nodes with lowest value of ω defined as have higher chance of removal, where and are the number of member type q needs to visit node i and the assigned importance rate of member type q, respectively. The parameter gets a larger value when more of member type q is available. The RR1 is introduced for sensitivity analysis in Table 2 on the APP's benchmark instances and is not practical for the COPTW as it deals with a homogeneous fleet of members. Resource Requirement Removal (RRR2) algorithm removes nodes with similar demands to increase the chance of being reinserted in alternative positions.

The Historical-Node Removal (HNR) algorithm uses the cost function for removing nodes, where the first two terms are normalised values from WTOR and WDR algorithms. The HNR removes nodes with worst where is the best position cost of node j before iteration i. The Time Windows-Oriented Removal (TWR) algorithm removes nodes with similar time windows using the same approach as the HNR algorithm. The TWR algorithm divides the total time window length into four zones and removes nodes from each zone, where denote respectively the total nodes and the unvisited nodes in the zone i. Every time one of the WDR, WTR, WTOR, WRR and RRR1 algorithms is called, we use a random binary variable to decide whether the calculated cost values should be divided by their associated score or not. The Last removal algorithm is the Cluster Removal (CR) heuristic. We use a density-based algorithm for clustering nodes (Ester et al. (1996)) which have proved to be efficient for classical OPs (Santini (2018)). We select a random cluster and remove D nodes from the chosen cluster in the current solution. The process of random selection of clusters and removal of nodes from the intersection of the cluster and the solution continues until D nodes are being removed.

The overall description of the ALNS is given in Algorithm 2. The ALNS algorithm determines an initial solution by the classical MB heuristic. Then, probabilities for insertion and removal heuristics are initialised (line 5). The preliminary temperature needs to be initialised for the SA algorithm and is updated (line 23) through the search process. Using the roulette-wheel mechanism as an adaptive layer of the ALNS a recreate and ruin heuristic is selected in every iteration (lines 10 and 11). Note that all variants of the MB heuristic in addition to greedy time heuristic fall within the set of recreate (repair) algorithms (line 9). Depending on the improvements resulting from the selected algorithms, scores are updated for the heuristics (line 18). Finally, the ALNS returns the best solution () found after N iterations (line 24).

Algorithm 2

The overall structure for the ALNS algorithm with simulated annealing.

Computational study

Extensive numerical studies were conducted to evaluate the efficacy of the proposed solution approach. A set of benchmark instances was generated by adding the problem-specific attribute to the well-known existing benchmark sets (see Vansteenwegen et al. (2009)). The resource requirement attribute was added to each vertex by picking 1, 2 or 3 randomly, which indicates how many members of the team are required to collect the associated reward at each node.1

In the first study, we investigate the performance of the proposed removal and insertion algorithms. Thereafter, we evaluate the efficiency of our ALNS approach. In the next study, truncated benchmark sets are designed to solve sufficiently small-sized instances by means of both the CPLEX commercial solver and our algorithm. The number of vertices in the small-sized instances are carefully chosen to obtain the optimal solutions using CPLEX. Furthermore, we explore the trade-off between an increased number of available members for service on the one hand and the computational time and objective value on the other hand. In addition, we demonstrate the performance of the proposed heuristic in terms of time and accuracy on the large-sized benchmark instances and present our results as benchmarks for future studies. All the above computations were performed on a single CPU with 16 GB of RAM on the Australian National Computational Infrastructure using a single thread. Each node is equipped with dual 8-core Intel Xeon (Sandy Bridge 2.6 GHz) processors. The algorithm was programmed in C++, using a GCC 6.2.0 compiler. Where applicable, MILP models were solved by the CPLEX 12.6 commercial solver in deterministic parallel optimisation mode. All tables show the execution times as elapsed time in seconds.

Parameter tuning and sensitivity analysis of algorithms

Many aspects of the COPTW are similar to those of the APP. Thus, it seemed logical to use the same parameter values identified in Roozbeh et al. (2018b). This also meant that we could perform a fair comparison on the performance of our new heuristics against the APP developed in that paper. The parameters had been determined using a similar tuning methodology to Emeç et al. (2016), Ropke and Pisinger (2006). Parameters and their tuned values are summarised below. Please note that similar parameter values to Roozbeh et al. (2018b) are used for validation purposes.

We investigate the performance of the insertion and removal heuristics in removing inefficient heuristics. This is done by conducting a sensitivity analysis on the heuristics using the same instances that were used for the parameter tuning: R104, R206, RC104, RC108 and RC206. Among the removal algorithms, WTOR, WDR and HNR are the most frequently used algorithms while CR, TWR and RR rarely contribute through the ALNS algorithm. The results show that all insertion heuristics contribute to the solution.

An algorithm with a minor contribution through the search process might still play a significant role in helping the ALNS to escape a local optimum Emeç et al. (2016). To identify whether omission of an algorithm might improve the final solution, we conduct further experiments omitting just one algorithm. The results in Table 2 indicate that exclusion of the MB heuristic has a negative impact on the worst, average and best values of the objective function in 10 runs. However, exclusion of the NCMB improves the ALNS performance. Thus, this heuristic, and others (indicated in the table with a “Y”), are deleted from the ALNS. The other heuristics (indicated with a “N”) make some positive contribution to the objective and remain in the ALNS.

Experiments on APP instances

We validate the performance of our algorithm against the results published in Roozbeh et al. (2018b) due to its similarity with the COPTW. In this paper an algorithm is proposed to solve an asset protection problem during escaped wildfires where community assets are endangered by the moving fire front. The objective function aims to maximise the total scores collected by protecting at-risk assets. An asset can be protected if a sufficient number of vehicles (i.e. members) from different types accomplish tasks within the time window. Our results are achieved by using the parameters in Table 1 after exclusion of inefficient heuristics. Table 3 shows a comparison of the results of our algorithm and the APP. We define two sets of vehicles as in the original paper named as Set1 and Set2 for the sake of a fair comparison. In order to compare results of each algorithm, the average and best solutions are noted. Results under the label “APP” referring to asset protection problem's results and ALNS refers to our new algorithm. The proposed ALNS for the COPTW has built on the work by Roozbeh et al. (2018b), and our newly developed heuristics along with the sensitivity analysis improved the ALNS performance significantly as can be seen in Table 3. Our algorithm is implemented to solve the asset protection instances and results compared to the adaptive large neighbourhood search algorithm developed for the APP by Roozbeh et al. (2018b).

Table 1

Parameter tuning results.

Description	Parameter	Value
Parameters for MB heuristic	(λ,μ,ϑ)	(2,1,3)
Improving solution score	σ2	12
Number of iterations	N	3000
Number of iterations over each segment	n	100
Roulette wheel reaction factor	ρ	0.1
Global solution score	σ1	35
Worse solution score	σ3	5
Shaw parameters	θ1, θ2, θ3	(3,13,7)
SA parameter	δ	0.05
Cooling rate	ϵ	0.9999
Noise parameter	αnoise	0.6
WDR determinism factor	κ	8
Shaw determinism factor	η	12

Table 3

A performance comparison of results for 100-nodes. APP represents results for asset protection problem and ALNS shows our results. Vehicle numbers are defined in two categories: Set1 = (V1=6, V2=5, V3=4) and Set2 = (V1=7, V2=6, V3=5). Results in the table are percentage values of scores collected and time is in seconds.

Instances	# Vehicles	APP			ALNS			Δavg	Δbest
		time	Avg	Best	time	Avg	Best
C100	Set1	138.47	60.17	61.84	135.83	62.34	63.18	3.48	2.12
	Set2	150.39	66.66	68.35	154.27	68.46	69.70	2.63	1.95
C200	Set1	133.48	59.04	60.72	130.64	61.06	61.92	3.31	1.94
	Set2	143.92	64.87	66.58	136.67	66.73	67.73	2.79	1.70
R100	Set1	134.47	61.19	62.50	125.48	63.39	64.49	3.47	3.09
	Set2	138.97	68.45	69.86	131.66	70.28	71.31	2.60	2.03
R200	Set1	135.75	63.64	65.30	129.17	65.80	66.75	3.28	2.17
	Set2	144.65	69.76	71.38	134.30	71.66	72.74	2.65	1.87
RC100	Set1	143.41	66.77	68.59	146.12	68.84	69.98	3.01	1.99
	Set2	149.53	73.21	75.17	152.94	74.94	75.85	2.31	0.90
RC200	Set1	142.97	67.17	69.13	130.35	69.35	70.45	3.14	1.87
	Set2	146.91	73.12	74.71	139.50	74.95	75.79	2.44	1.42

In Table 3, and are the percentage deviation from the reported results in the literature. Our ALNS improves results for all instances and enhances the average results for the and by around 2% and 3%, respectively. Table 3 shows that our approach performs more efficiently and validates our algorithm on the existing benchmarks. The computation times are quite close for both methods, while the ALNS invests more time on heuristics that bring improvement. It would appear that keeping the best performing heuristics, while discarding others, identified in the sensitivity analysis has yielded a positive improvement in performance.

Experiments on small-size instances of the COPTW

For further evaluation, a collection of small-sized benchmark instances was generated and each instance was solved by means of both CPLEX and the ALNS algorithm. Note that in this problem set, Euclidean distance is used and rounded down to the first decimal, while results are rounded to two decimal places.

A summary of the tests for 10 and 12 nodes with 3 and 4 team members on instance sets C100, R100 and RC100 is provided in the Table 4. The sizes of truncated instances are chosen in such a way as to investigate the correlation between the increase in problem size and the growth in computational effort. It is worthwhile to mention that infeasible edges are excluded in the MILP formulations to simplify models for the CPLEX implementation. Also, larger problems with more than 12 nodes could not be solved by CPLEX within the time limit of 5 hours.

Table 4

A summary of the ALNS performance for small-size instances for 10 and 12 nodes on C100, R100 and RC100 datasets. All computational times are in seconds.

Set	# Members	10				12
		CPLEX	ALNS	Δaverage	Δbest	CPLEX	ALNS	Δaverage	Δbest
C100	p=3	62.15	5.48	0.00	0.00	142.37	6.86	0.00	0.00
	p=4	80.16	61.24	0.00	0.00	61.75	113.21	0.00	0.00
R100	p=3	244.29	4.55	0.00	0.00	5809.95	6.91	0.00	0.00
	p=4	355.91	28.78	-0.02	0.00	6018.78	43.48	-0.20	0.00
RC100	p=3	356.87	26.33	0.00	0.00	4502.62	21.19	-0.28	0.00
	p=4	277.14	60.99	-0.13	0.00	3297.12	111.34	-0.10	0.00

Computational times for both CPLEX and ALNS are reported in Table 4. The optimality gaps are shown by “” and “” that represent the percentage gap between the CPLEX optimal solution and ALNS algorithm average and best results. The average gap for and are 0.06% and 0.00%, a promising performance by the proposed algorithm. Moreover, it can be seen that the computational time increases significantly with minor changes in the problem size for CPLEX compared with negligible changes for the ALNS.

Table 5, gives a summary of results for 24 and 26 vertices with the same number of available team members. One can see that CPLEX solves larger problems from the sets C200, R200 and RC200 compared to those in Table 4. This is due to the nature of the studied class of problems as the time window intervals are different in length and a larger portion of nodes can be covered by the same number of team members. In Table 5, the average computational time remains around two minutes for all instances, while it takes hours to solve some sets by CPLEX. In Table 5, the average deviation of the ALNS algorithm from optimal solutions is 0.00% for and just 0.03% for which are reasonable. Due to the CPLEX computational time we are unable to solve larger instances in Table 5 while for some instances large amounts of nodes are covered.

Table 5

A summary of the ALNS performance for small-size instances for 24 and 26 nodes on C200, R200 and RC200 dataset. All computational times are in seconds.

Set	# Members	24				26
		CPLEX	ALNS	Δaverage	Δbest	CPLEX	ALNS	Δaverage	Δbest
C200	p=3	140.95	51.79	0.00	0.00	217.15	55.87	0.00	0.00
	p=4	145.81	92.63	0.00	0.00	144.80	104.69	0.00	0.00
R200	p=3	351.03	117.21	0.00	0.00	473.92	128.09	0.00	0.00
	p=4	224.51	132.4	0.00	0.00	286.45	134.96	0.00	0.00
RC200	p=3	509.47	130.13	0.00	0.00	4516.54	147.05	-0.38	0.00
	p=4	313.45	143.61	0.00	0.00	336.61	168.53	-0.04	0.00

An instance where the ALNS achieved an optimal solution on a small set is demonstrated in Fig. 4. For the sake of better presentation, nodes are assigned into the cells of arrays. The times for starting the service and the corresponding time windows for each succeeding vertex are provided. Consider, for example Fig. 4, three members leave the depot and arrive at node number 5 within its time windows [15,67], where they start the service simultaneously to collect the associated score at time 15.1. After that, second and third members leave node 5 toward vertex number 3. Finally, all the team members finish their tour by returning to the depot before , where is the closing time at the final depot, time budget constraint. It can be seen that the illustrated routes are highly dependent such that any minor changes in sequence of nodes at any tour, requires reconstruction of other routes.

Figure 4

Scheduled tour by ALNS.

Experiments on large-size instances of the COPTW

We performed further experiments on larger instances. Table 6 shows the computational results with 4 and 6 team members for 100 nodes on instance sets C100, C200, R100, R200, RC100 and RC200. We utilise the ALNS benchmark since no benchmark exists for the COPTW. Proportional to the problem size more team members are considered in order to cover a substantial percentage of available scores. Our results demonstrate the worst, average and best percentage value of scores collected for each set. In the last two columns of Table 6, we report on the most effective insertion and removal algorithms for each set of instances. Merit-based heuristics and noise-imposed cumulative merit-based are recognised as the best insertion algorithms for all members and sets. It shows that using the merit-based heuristic directly, improves the value of objective function.

Table 6

Computational results for the larger instances with 100 vertices. Results are shown as the percentage of rewards collected.

Set	# Members	100
		worst (%)	average (%)	best (%)	best removal algorithms	best insertion algorithms
C100	p=4	35.36	35.98	36.71	WDR/HNR/WTR/WTOR	MB/NMB/GT/NCMB
	p=6	46.41	46.95	47.82	WDR/WTR/TR/SR	MB/NCMB/NMB/GT
C200	p=4	72.17	73.25	74.65	WDR/WTR/HNR/RR1	MB/NMB/NCMB/GT
	p=6	85.84	87.00	88.26	WDR/HNR/WTR/WTOR	MB/NMB/NCMB/GT
R100	p=4	32.48	33.48	34.44	WDR/WTR/WRR/RRR1	MB/NCMB/NMB/GT
	p=6	42.26	43.40	44.71	WDR/RRR1/WTR/WTOR	MB/NMB/GT/NCMB
R200	p=4	74.26	75.99	78.32	WDR/WTR/SR/TR	MB/NCMB/NMB/GT
	p=6	87.00	88.57	90.61	WDR/WTR/RRR1/SR	MB/NMB/GT/NCMB
RC100	p=4	30.12	30.78	31.51	HNR/WTR/WTOR/WDR	MB/NMB/GT/NCMB
	p=6	39.96	40.74	41.51	WDR/WTOR/HNR/RRR1	MB/NMB/GT/NCMB
RC200	p=4	66.63	68.37	70.44	WDR/WTR/WTOR/TR	MB/NMB/GT/NCMB
	p=6	81.70	83.61	86.09	WTR/WDR/RRR1/SR	MB/NCMB/NMB/GT

Conclusion

The COPTW is an important class of the orienteering problem that arises naturally in many important applications, such as home health care problems, the asset protection problem, UAV problems, and in other emergency logistics management. In practical problems many of these applications require solutions as a matter of urgency. Current algorithms are not designed to handle the complexities resulting from the requirements of synchronised visits in the COPTW. Thus, we initially explained handling the synchronous visits within an algorithmic framework in Section 3.2. Then, we developed an ALNS algorithm equipped with variants of the MB heuristic as insertion operators to deal with the various issues that arise from this problem. Together with the MB heuristic 19 heuristics were implemented using the problem specific attributes, where five of them are newly introduced. Thereafter, we evaluated the effectiveness of the removal and insertion algorithms to improve the efficiency of the ALNS by investing more time on effective heuristics. The tailored ALNS algorithm can solve large-sized problems in computational times suitable for operational practice. The average time spent to solve each instance with 100 nodes in Table 6 was 193.6 seconds, which is a time suitable for operational purposes.

In order to validate our solution approach, we modified the algorithm to implement tests on the existing benchmarks. We conducted tests on the APP's instances, a special case of the COPTW. In the APP any minor improvements in the objective function is significant and important as it means that more assets can be protected. Our algorithm improved the existing benchmarks using modest computational resources. For further evaluation, a new benchmark set was generated for the COPTW. The performance of the algorithm was validated successfully for small-sized instances by comparing solutions with the optimal results obtained by the CPLEX solver. Further experiments with large scale problems demonstrated the efficacy as well as accuracy of the ALNS in terms of achieved objective values and computation times.

As the applications grow in the cooperative orienteering class of problems, other developments may be needed to handle complexities such as soft time windows and real-time changes to routes. The ALNS algorithm is a powerful tool which further constraints might be added to deal with such problems.

Declarations

Author contribution statement

I. Roozbeh: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data; Wrote the paper.

J.W. Hearne, D. Pahlevani: Conceived and designed the experiments; Analyzed and interpreted the data; Wrote the paper.

Funding statement

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Competing interest statement

The authors declare no conflict of interest.

Additional information

Data associated with this study has been deposited at https://sites.google.com/site/imanrzbh/datasets.