user: text: Write a concise one-sentence TLDR summary of the scientific paper content below.

Title: 

Paper content:
Sequential decision problems for real-world applications often need to be solved in real-time, requiring algorithms to perform well with a restricted computational budget.
Width-based lookaheads have shown state-of-the-art performance in classical planning problems as well as over the Atari games with tight budgets.
In this work we investigate width-based lookaheads over Stochastic Shortest paths (SSP).
We analyse why width-based algorithms perform poorly over SSP problems, and overcome these pitfalls proposing a method to estimate costs-to-go.
We formalize width-based lookaheads as an instance of the rollout algorithm, give a definition of width for SSP problems and explain its sample complexity.
Our experimental results over a variety of SSP benchmarks show the algorithm to outperform other state-of-the-art rollout algorithms such as UCT and RTDP.
Model-based lookahead algorithms provide the ability to autonomously solve a large variety of sequential decision making problems.
Lookaheads search for solutions by considering sequences of actions that can be made from the current state up to a certain time into the future.
For realworld applications decisions often need to be computed in real-time, requiring algorithms to perform with a restricted computational budget.
Limiting search in this way can result in considering states and trajectories which do not provide useful information.
To address this, lookaheads can be augmented with heuristics that estimate costs-to-go to prioritise states and trajectories, and have been shown to perform well where computation budgets are restricted BID8 .This paper is concerned with Stochastic Shortest Path (SSP) problems which are often used to compare and evaluate search algorithms.
We consider the width-based family of planning algorithms, first introduced by BID15 , which aim to prioritise the exploration of novel areas of the state space.
Two width-based planners, Lipovetzky and Geffner's breadth-first search, IW(1), and the depth-first search, Rollout-IW(1) BID1 , are investigated on SSP problems.
We first provide the necessary background for SSP problems and width-based algorithms, while also formalising width-based algorithms as instances of the rollout algorithm BID4 .
We then show the motive to augment width-based lookaheads with cost estimates on SSP problems, define the width of SSP problems and propose a novel width-based algorithm that estimates costs-to-go by simulating a general base policy.
Our experimental study shows that the algorithm compares favourably to the original Rollout-IW(1) algorithm and to other state-of-the-art instances of the rollout algorithm.
We concern ourselves with the problem of decision under stochastic uncertainty over a finite number of stages, which we characterise following closely the presentation of BID4 .
We are given a discrete-time dynamic system x k+1 = f k (x k , u k , w k ), k = 0, 1, . . .
, N − 1 (1) where the state x k is an element of a space S k ⊂ R d , the control u k is an element of space C k ⊂ N, and the random disturbance w k is an element of a space D k ⊂ R m 1 .
The control u k is constrained to take values in a given non-empty subset U (x k ) ⊂ C k , which depends on the current state x k , so that u k ∈ U k (x k ) for all x k ∈ S k and k. The random disturbance w k is characterised by a probability distribution P k (·|x k , u k ) that may depend explicitly on x k and u k but not on the values of previous disturbances w k−1 , . . .
, w 0 .
We consider the class of policies, or control laws, corresponding to the sequence of functions DISPLAYFORM0 where µ k maps states x k into controls u k = µ k (x k ) and is such that µ k (x k ) ∈ U (x k ) for all x k ∈ S k .
Such policies will be called admissible.
Given an initial state x 0 and admissible policy π, the states x k and disturbances w k are random variables with distributions defined through the system equation DISPLAYFORM1 Thus, for given functions g f (terminal cost) and g the expected cost of π starting at x 0 is DISPLAYFORM2 where the expectation is taken over the random variables w k and x k .
An optimal policy π * is one that minimises this cost DISPLAYFORM3 where Π is the set of all admissible policies.
The optimal cost J * (x 0 ) depends on x 0 and is equal to J π * (x 0 ).
We will refer to J * as the optimal cost or optimal value function that assigns to each initial state x 0 the cost J * (x 0 ).
We use BID4 definition, that formulates Stochastic Shortest Path (SSP) problems as the class of optimal control problems where we try to minimize DISPLAYFORM0 with α set to 1 and we assume there is a cost-free termination state t which ensures that J π (x 0 ) is finite.
Once the system reaches that state, it remains there at no further cost, that is, f (t, u, w) = t with probability 1 and g(t, u, w) = 0 for all u ∈ U (t).
We note that the optimal control problem defined at the beginning of this section is a special case where states are pairs (x k , k) and all pairs (x N , N ) are lumped into termination state t.
In order to guarantee termination with probability 1, we will assume that there exists an integer m such that there is a positive probability that t will be reached in m stages or less, regardless of what π is being used and the initial state x 0 .
That is, for all admissible policies and i = 1, ..., m it holds DISPLAYFORM1 A policy π will be proper if the condition above is satisfied for some m, and improper otherwise.
A particularly effective on-line approach to obtain suboptimal controls is rollout, where the optimal cost-to-go from current state x k is approximated by the cost of some suboptimal policy and a d-step lookahead strategy.
The seminal RTDP BID2 algorithm, is an instance of the rollout strategy where the lookahead is uniform, d = 1, and controlsμ(x k ) selected at stage k and for state x k are those that attain the minimum DISPLAYFORM0 whereJ k+1 is an approximation on the optimal cost-to-go J * k+1 .
If the approximation is from below, we will refer to it as a base heuristic, and can either be problem specific BID8 , domain independent BID5 BID22 or learnt from interacting with a simulator BID18 .
Alternatively, J k+1 can be defined as approximating the cost-to-go of a given suboptimal policy π, referred to as a base policy, where estimates are obtained via simulation BID19 .
We will denote the resulting estimate of cost-to-go as H k (x k ) 2 .
The result of combining the lookahead strategy and the base policy or heuristic is the rollout policy,π {μ 0 , µ 1 , . .
.,μ N −1 } with associated costJ(x k ).
Such policies have the property that for all x k and k DISPLAYFORM1 when H k is approximating from above the cost-to-go of a policy, as shown by BID4 from the DP algorithm that defines the costs of both the base and the rollout policy.
To compute at time k the rollout controlμ(x k ), we compute and minimize over the values of the Q-factors of state and control pairs (x l , u l ), DISPLAYFORM2 for admissible controls u l ∈ U (x l ), l = k + i, with i = 0, ..., d − 1, and DISPLAYFORM3 for l = k +d. In this paper we make a number of assumptions to ensure the viability of lookaheads with d > 1.
We will assume that we can simulate the system in Eq. 3 under the base policy, so we can generate sample system trajectories and corresponding costs consistent with probabilistic data of the problem.
We further assume that we can reset the simulator of the system to an arbitrary state.
Performing the simulation and calculating the rollout control still needs to be possible within the real-time constraints of the application, which is challenging as the number of Q-factors to estimate and minimizations to perform in Equations 9-10 is exponential on the average number of controls available per stage and d, the maximum depth of the lookahead.
We avoid the blowup of the size of the lookahead by cutting the recursion in Equation 9 and replacing the right hand side by that of Equation 10.
As detailed in the next section, we will do this when reaching states x l that are deemed not to be novel according to the notion of structural width by BID15 .
This results in a selective strategy alternative to the upper confidence bounds BID0 used in popular instances of Monte-Carlo Tree Search (MCTS) algorithms like Kocsis and Szepesvari's (2006) UCT, that also are instances of the rollout algorithm BID4 ).
We instantiate the rollout algorithm with an l-step, depth-selective lookahead policy using Width-based Search BID15 .
These algorithms both focus the lookahead and have good any-time behaviour.
When it comes to prioritisation of expanding states, width-based methods select first states with novel valuations of features defined over the states BID16 BID9 DISPLAYFORM0 is the domain of variable x i , has not been seen before in the current search.
Note that novel states are independent of the objective function used, as the estimated cost-to-go J is not used to define the novelty of the states.
IW(1) has recently been integrated as an instance of a rollout algorithm, and has been shown to perform well with respect to learning approaches with almost real-time computation budgets over the Atari games BID1 .
The breadth-first search strategy underlying IW(1) ensures a state variable x i is seen for the first time through the shortest sequence of control steps, i.e. the shortest path assuming uniform costs g(x, u, w).4 On the other hand, depth-first rollout algorithms cannot guarantee this property in general.
Rollout IW (RIW) changes the underlying search of IW into a depth-first rollout.
In order to ensure that RIW(1) considers a state to be novel iff it reaches at least one value of a state variable x i l through a shortest path, we need to adapt the definition of novelty.
Intuitively, we need to define a set of state features to emulate the property of the breadth-first search strategy.
Let d(x i , v) be the best upper bound known so far on the shortest path to reach each value v ∈ D(x i ) of a state variable from the root state x k .
Initially d(x i , v) = N for all state variables, where N is the horizon which is the maximum search depth allowed for the lookahead, thus denoting no knowledge initially.
When a state x l is generated, d(x i , v) is set to l for all state variables where l < d(x i , v).
Since RIW(1) always starts each new rollout from the current state x k , in order to prove a state x l to be novel we have to distinguish between x l being already in the lookahead tree and x l being new.
If x l is new in the tree, to conclude it is novel, it is sufficient to show that there exists a state variable x i whose known shortest path value d(x i , v) >
l. If x l is already in the tree, we have to prove the state contains at least one state variable value x i whose shortest path is l = d(x i , v), i.e. state x l is still novel and on the shortest path to x i .
Otherwise the rollout is terminated.
In order to ensure the termination of RIW (1) , non-novel states are marked with a solved label.
The label is backpropagated from a state x l+1 to x l if all the admissible control inputs u ∈ U (x l ) yield states x l+1 = f l (x l , u, w l ) already labeled as solved.
RIW(1) terminates once the root state is labeled as solved BID1 .
Nonnovel states x l are treated as terminals and their cost-to-go is set to 0.
This can induce a bias towards non-novel states rather than true terminal states.
In the next section we investigate how to overcome the ill-behaviour of RIW(1) when a state x l is non-novel.
We discuss the importance of estimating upper-bounds on the cost-to-go H l (x l ) instead of assigning termination costs.
This turns out to be essential for RIW (1) over SSPs.
Despite the successes of width-based algorithms on a variety of domains including the Atari-2600 games BID16 BID1 , the algorithms, as will be shown, have poor performance on SSP problems.
We illustrate this with two scenarios.
First, width-based lookaheads prefer trajectories leading to non-novel states over longer ones that reach a goal.
Second, and specific to depth-first width-based lookaheads, we show that useful information is ignored.
We can demonstrate these scenarios using a simple SSP problem with uniform and unitary action costs, shown in FIG0 .
The task is to navigate to a goal location using the least number of left, right, up or down actions.
Any action that would result in the agent moving outside of the grid produces no change in its position.
The features used by the width-based planners are the coordinates for the current agent position.
Both IW(1) and RIW(1) algorithms, given a sufficient budget, would result in the lookahead represented by yellow lines in FIG0 .
As expected, both lookaheads contain the shortest paths to make each feature of the problem true.
For both IW(1) and RIW(1), we back up the costs found in the lookahead starting from terminal and non-novel states.
In this instance a move down or left from the agent's initial state has no effect, thus immediately producing a non-novel state.
When backing up values, down and left have an expected cost of 1, which is less than the optimal cost of 2 for up, the action that leads to the top left goal state.
This prevents both IW(1) and RIW(1) from ever achieving the goal, as they keep selecting those useless actions.
Furthermore, if the goal is the top right location in FIG0 , RIW(1)'s random action selection can generate a trajectory that reaches the goal.
Yet, trajectories leading to the goal are pruned away, as non-novel states in later considered trajectories are treated as terminals, again resulting in the lookahead represented by the yellow lines in FIG0 .
BID1 introduced the algorithm RIW in the context of deterministic transition functions.
In this section we discuss its properties in the context of SSPs.
The set of features used to evaluate the novelty of a state is DISPLAYFORM0 is the domain of variable x i , and d is a possible shortest path distance.
Note that the horizon N is the upper-bound of d. The maximum number of novel states is O(|F |), as the maximum number of shortest paths for a feature (v, i, ·) ∈ F is N .
That is, in the worst case we can improve the shortest path for (v, i, ·) by one control input at a time.
The labeling of nodes ensures the number of rollouts from the initial state in RIW (1) is at most O(|F | × b), where b = max x l |U (x l )| is the maximum number of applicable control variables in a state, i.e. maximum branching factor.
When the labeling is applied to stochastic shortest path problems, the resulting lookahead tree is a relaxation of the original SSP, as it allows just one possible outcome of a control input.
Alternatively, one can back-propagate the label solved to a state x l iff 1) all admissible control inputs u ∈ U (x l ) have been applied resulting in states labeled as solved, and 2) the tree contains all the possible resulting states of each control input u ∈ U (x l ).
We refer to this new strategy to backpropagate labels as λ-labeling.
We denote as λ the maximum number of states that can result from applying u ∈ U (x l−1 ) in a state x l .
That is, λ = max x,u,w |f (x, u, w)|.
RIW(1) with λ-labeling will terminate after at most DISPLAYFORM1 Furthermore, we can reconcile the notion of width over classical planning problems BID15 with SSPs.
A terminal state t made of features f ∈ F has width 1 iff there is a trajectory x 0 , u 0 , . . . , u n−1 , x n for n ≤ N where x n = t, such that for each x j in the trajectory 1) the prefix x 0 , u 0 , . . .
, u j−1 , x j reaches at least one feature DISPLAYFORM2 , it is a shortest path possible to reach a value in x i , 2) any shortest path to f j can be extended with a single control input u into a shortest path for a feature f j+1 complying with property 1) in state x j+1 , and 3) the shortest path for f n is also a shortest path for termination state t. RIW(1) with the new labeling strategy is guaranteed to reach all width 1 terminal states t. Theorem 1.
Rollout IW(1) with λ-labeling is guaranteed to reach every width 1 terminal state t in polynomial time in the number of features F if λ = 1.
If λ = ∞, RIW(1) will not propagate any solved label, and terminate only when the computational budget is exhausted.
For simplicity, we assumed shortest paths are equivalent to the shortest sequence of control inputs.
To generalize to positive non-uniform costs, the distance d in the features should keep track of the cost of a path i g(x i , u i , w i ) instead of its length, and the horizon be applied to the cost of the path.
The most successful methods for obtaining cost-to-go approximations have revolved around the idea of running a number of Monte Carlo simulations of a suboptimal base policy π BID10 BID7 .
This amounts to generating a given number of samples for the expression minimized in Equation 7 starting from the states x l over the set of admissible controls u l ∈ U (x l ) in Equation 10, averaging the observed costs.
Three main drawbacks BID4 follow from this strategy.
First, the costs associated with the generated trajectories may be wildly overestimating J * (x l ) yet such trajectories can be very rare events for the given policy.
Second, some of the controls u l may be clearly dominated by the rest, not warranting the same level of sampling effort.
Third, initially promising controls may turn out to be quite bad later on.
MCTS algorithms aim at combining lookaheads with stochastic simulations of policies and aim at trading off computational economy with a small risk of degrading performance.
We add two new methods to the MCTS family, by combining the multi-step, width-based lookahead strategy discussed in the previous section with two simulation-based cost-to-go approximations subject to a computational budget that limits the number of states visited by both the lookahead and base policy simulation.
The first method, which we call RIW-RW, uses as the base policy a random walk, a stochastic policy that assigns the same probability to each of the controls u admissible for state x, and is generally regarded as the default choice when no domain knowledge is readily available.
A rolling horizon H is set for the rollout algorithm that follows from combining the RIW(1) lookahead with the simulation of random walks.
The maximal length of the latter is set to H − l, where l is the depth of the non-novel state.
Both simulations and the unrolling of the lookahead are interrupted if the computational budget is exhausted.
While this can result in trajectories sometimes falling short from a terminal, it keeps a lid on the possibility of obtaining extremely long trajectories that eat into the computational budget allowed and preclude from further extending the lookahead or sampling other potentially more useful leaf states x l .
One of the most striking properties of rollout algorithms is the cost improvement property in Equation 8, suggesting that upper bounds on costs-to-go can be used effectively to approximate the optimal costs J * .
Inspired by this, the second width-based MCTS method we discuss leverages the sampling techniques known as stochastic enumeration (SE) BID19 to obtain an unbiased estimator for upper bounds on costs-to-go, or in other words, estimates the maximal costs a stochastic rollout algorithm with a large depth lookahead can incur.
SE methods are inspired by a classic algorithm by D. E. Knuth to estimate the maximum search effort by backtracking search (1975).
Knuth's algorithm estimates the total cost of a tree T with root u keeping track of two quantities, C the estimate of the total cost, and D the expectation on the number of nodes in T at any given level of the tree, and the number of terminal nodes once the algorithm terminates.
Starting with the root vertex u and D ← 1, the algorithm proceeds by updating D to be D ← |S(u)|D and choosing randomly a vertex v from the set of successors S(u) of u. The estimate C is then updated C ← C + c(u, v)D using the cost of the edge between vertices u and v. These steps are then iterated until a vertex v is selected s.t.
S(v ) = 0.
We observe that Knuth's C quantity would correspond to the worst-case cost-to-goJ(x) k of a rollout algorithm using a lookahead strategy with d set to the rolling horizon H and the trivial base heuristic that assigns 0 to every leaf state.
Furthermore, we assume that the algorithm either does not find any terminals within the limits imposed by the computational budget assigned, or if it finds one such state, it is too the very last one being visited.
Lookaheads define trees over states connected by controls, edge costs c(u, v) correspond directly with realisations of the random variable g(x, u, w) and the set of successors S(v) of a vertex corresponds with the set of admissible controls U (x).
While Knuth's algorithm estimates are an unbiased estimator, the variance of this estimator can be exponential on the horizon, as the accuracy of the estimator lies on the assumption that costs are evenly distributed throughout the tree BID19 .
In the experiments discussed next, we use Knuth's algorithm directly to provide H k (x k ), adding the stopping conditions to enforce the computational budget and limiting the length of trajectories to H − l as above.
In comparison with simulating the random walk policy, the only overhead incurred is keeping up-to-date quantities C and D with two multiplications and an addition.
To evaluate the different methods we use a number of GridWorld (Sutton and Barto 2018) domains, an instance of a SSP problem.
The goal in GridWorld is to move from an initial position in a grid to a goal position.
In each state 4 actions are available: to move up, down, left or right.
Any action that causes a move outside of the grid results in no change to the agent's position.
Actions have a cost of 1, with the exception of actions that result in reaching the goal state, that have a cost of 0.
The complexity of GridWorld can be scaled through the size of the grid and the location and number of goals.
GridWorld also allows for extensions, which we use to have domains with a stationary goal, moving goals, obstacles and partial observability.
For each instance of the GridWorld domain we have a d 0 × d 1 grid, and the state is the current location of the agent, x = (a 0 , a 1 ) where a i is the agent's position in dimension i. The transition function is formalised as DISPLAYFORM0 where, ef specifies the change in the agent's position for each action, T k ⊂ S k is the set of goal states and x k+1 = x k where the condition in FIG0 is not met.
The cost of a transition is defined as DISPLAYFORM1 For the stationary goal setting we have a single goal state which is positioned in the middle of the grid by dividing 223.2 ± 10.4 199.9 ± 13.6 145.5 ± 12.9 DISPLAYFORM2 where δ t k gives the relative change of the goal state t k for the time step k + 1 and DISPLAYFORM3 Resulting in two goals starting at opposite corners of the grid moving back and forth on the same diagonal.
The obstacles setting, uses the stationary goal, but modifies S k such that, DISPLAYFORM4 where O ⊂ N 2 and is the set of obstacles, that is grid cells in which the agent is not allowed.
Having partially observable obstacles in GridWorld provides an instance of the stochastic Canadian Traveller Problem (CTP) (Papadimitriou and Yannakakis 1991).
The objective in CTP is to find the shortest path between one location in a road map to another, however, there is a known probability for each road in the map that due to weather conditions the road is blocked.
A road in CTP can only be observed as being blocked or unblocked by visiting a location connected to it, and once a road status is observed the status remains unchanged.
In terms of the GridWorld problem, each grid cell has a known probability of being a member of the obstacle set, O. The agent can only observe cells as being obstacles or not when it is in a neighbouring cell.
Once a grid cell is observed it is then known that it is either an obstacle or not for the remaining duration of the episode.
John Langford designed two MDP problems 5 described as Reinforcement Learning (RL) "Acid" intended to be RIW NA 212.9 ± 11.3 215.9 ± 10.8 223.5 ± 9.8 Rnd.236.2 ± 6.5 200.4 ± 11.4 110.6 ± 11.8 Table 2 : Same experimental setting as TAB2 over GridWorld with a moving goal.difficult to solve using common RL algorithms, such as Qlearning.
Langford's two problems allow two actions from every state.
The state space for the problems is S k = {i | 0 ≤ i < N } where the number of states value, N , allows the complexity of the problem to be controlled.
Langford originally specified the problems as reward-based, here we modify them to be SSP cost-based problems.
Reward shaping is commonly used to make Reinforcement Learning easier by encouraging actions, through higher rewards, towards a goal state or states.
The first of Langford's problems is named Antishaping and uses reward shaping to encourage actions away from the goal state.
Antishaping has the transition function DISPLAYFORM5 otherwise, the state remains unchanged, x k+1 = x k .
The set containing the goal state is T k = {N − 1}, which can be achieved by continuously selecting u k = 0.
The cost of each transition in Antishaping is 0.25 divided by N -x k+1 , except when x k+1 = N − 1 where the cost is 0.
The problem becomes a large plateau where longer paths become more costly at larger rates.
The motivation behind Langford's second problem, Combolock, is if many actions lead back to a start state, random exploration is inadequate.
The Combolock problem has the transition function DISPLAYFORM6 otherwise x k+1 is equal to the initial position of 0.
The goal state is T k = {N − 1} and sol x k is either 0 or 1 assigned to state x k which remains constant.
For each state x ∈ S, sol x has an equal chance of either being 0 or 1.
The cost of each transition in Combolock is 1 except for the transition that leads to the terminal state N − 1 where the cost is 0.
While common Reinforcement Learning algorithms such as Q-Learning methods will struggle to solve these domains, it is claimed by Langford that the E 3 (Kearns and Singh 2002) family of algorithms, whose exploration do not rely solely on random policies or reward feedback but on exploring the maximum number of states, will perform well.
43.9 ± 1.9 35.9 ± 2.5 25.3 ± 2.4 UCT Man.
37.0 ± 2.9 36.4 ± 2.9 36.4 ± 2.9 Rnd.43.8 ± 1.9 38.5 ± 2.5 25.9 ± 1.9 RIW Man.
36.4 ± 2.9 36.4 ± 2.9 36.4 ± 2.9 NA 49.8 ± 0.4 48.8 ± 1.0 49.3 ± 0.8 Rnd.44.9 ± 1.7 34.5 ± 2.8 19.3 ± 2.1 20 1Stp Man.
76.7 ± 5.5 77.1 ± 5.4 76.7 ± 5.5 Rnd.97.9 ± 1.6 88.0 ± 3.4 62.7 ± 4.8 UCT Man.
78.0 ± 5.4 78.4 ± 5.3 73.4 ± 5.7 Rnd.98.7 ± 1.2 96.4 ± 1.9 77.2 ± 4.2 RIW Man.
79.7 ± 4.9 76.7 ± 5.5 76.7 ± 5.5 NA 100.0 ± 0.0 100.0 ± 0.0 99.6 ± 0.8 Rnd.98.5 ± 1.3 88.0 ± 3.4 29.3 ± 3.1
Man.194.6 ± 13.4 191.5 ± 13.6 196.6 ± 13.2 Rnd.249.2 ± 1.1 244.4 ± 3.7 216.8 ± 9.3 UCT Man.194.6 ± 13.4 195.6 ± 13.3 184.4 ± 13.9 Rnd.249.0 ± 1.9 243.1 ± 4.3 231.6 ± 7.9 RIW Man.
208.6 ± 10.7 210.6 ± 11.1 193.5 ± 13.4 NA 250.0 ± 0.0 250.0 ± 0.0 250.0 ± 0.0 Rnd.247.9 ± 2.6 242.9 ± 4.3 196.1 ± 11.3 Table 3 : Same settings as TAB2 over GridWorld with fully observable obstacles and a stationary goal.
We evaluate the depth-first width-based rollout algorithm, RIW(1), with and without being augmented using base policies.
λ = 1 is used for the labels back-propagation.
We did not observe statistically significant changes with λ = ∞. For the GridWorld domain we define the features on which RIW(1) plans over as F = {(a, i, d) | a ∈ D(x i )} where d is the length of the control input path from the initial state, a is the agent's position in the grid in dimension i and D(x i ) is the domain of the agent's position, a, in dimension i. For Antishaping and Combolock the feature set will be F = {(i, d) | i ∈ N )} where i is the state number the agent is in and N is the number of states of the domain.
Two additional rollout algorithms are also considered, the one-step rollout algorithm, RTDP BID2 and the multi-step, selective, regret minimisation, rollout algorithm Upper Confidence bounds applied to Trees (UCT) BID14 .
The exploration parameter of UCT is set to 1.0 for all experiments.
For all the methods that use a base policies the maximum depth of a simulated trajectory is equal to H − l, where l is the depth at which the simulated trajectory began and H is the horizon value of the lookahead.
Also, a single, as opposed to multiple, simulated trajectory for the cost-to-go approximation is used, as initial results indicated it is favourable.
We also report the algorithms using a Manhattan distance heuristic for the GridWorld domains that use obstacles.
Using the Manhattan distance for the GridWorld problems with obstacles provides a lower bound on the cost-to-go.
Each method on the domains is evaluated at different levels of complexity by varying the number of states.
The methods are evaluated using different simulator budgets.
The simulator budgets are the maximum simulator calls allowed for the evaluation at each time step.
For each algorithm and domain Dim.
Heu.
Simulator Budget 100 1000 10000 10 1Stp Man.
36.6 ± 2.9 35.0 ± 3.0 35.9 ± 2.9 Rnd.40.4 ± 2.1 27.2 ± 2.5 15.4 ± 1.6 UCT Man.
28.4 ± 2.9 29.5 ± 3.1 18.5 ± 2.7 Rnd.41.5 ± 2.0 36.1 ± 2.3 22.2 ± 2.0 RIW Man.
28.1 ± 3.0 28.3 ± 3.0 26.5 ± 3.0 NA 49.5 ± 0.7 49.1 ± 0.9 49.8 ± 0.4 Rnd.43.5 ± 1.9 23.4 ± 2.6 11.1 ± 1.3
Man.
71.8 ± 5.8 74.0 ± 5.7 71.9 ± 5.8 Rnd.97.0 ± 1.7 82.3 ± 3.9 49.8 ± 4.6 UCT Man.
53.8 ± 5.7 60.9 ± 6.2 38.0 ± 5.5 Rnd.98.0 ± 1.3 87.4 ± 4.0 63.0 ± 4.8 RIW Man.
53.5 ± 5.6 44.1 ± 5.5 44.1 ± 5.6 NA 100.0 ± 0.0 99.5 ± 1.0 99.5 ± 1.0 Rnd.97.6 ± 1.5 79.7 ± 4.4 20.7 ± 1.9 Results are also presented for the Atari-2600 game Skiing, which is a SSP problem.
We use the OpenAI gym's BID6 interface of the Arcade Learning Environment (ALE) BID3 ) and use the slalom game mode of Skiing.
In the slalom mode the aim is to ski down the course as fast as possible while going through all the gates.
Once the finish line is reached, a 5 second time penalty is applied for each gate that is missed.
The reward values provided by ALE for Skiing is the negative value of the total time taken plus any time penalties in centiseconds, which we use as a cost.
We use the environment settings as described by BID17 with a frame skip of 5 and a probability of 0.25 of repeating the previous action sent to environment instead of the current one, which Machado et al. call sticky actions.
For evaluation we use a simulator budget of 100 and partial caching as described by BID1 , in that we cache simulator state-action transitions, thus assuming determinism, but clear the cached transitions when executing an action in the environment.
However, the lookahead tree itself is not cleared when executing an action in the environment as is done for the other domains trialed.
The maximum episode length is capped at 18,000 frames with a frame skip of 5 this equals 3,600 actions.
Using a simulation based cost-to-go approximation is infeasible with a simulator budget of 100 and the maximum episode length of 3,600 actions.
Therefore we report the algorithms using a heuristic cost-to-go estimate, which is the the number of gates that have either been missed or are still left times the time penalty of 500 centiseconds.
For the RIW(1) algorithms we use the pixel values from the current gray scaled screen at full resolution, that is 180 by 210 pixels, as features.
All experiments were run within the OpenAI gym frame- 1.7 ± 0.1 1.3 ± 0.1 1.1 ± 0.1 work BID6 ) and the code used for the algorithms and domains is available through GitHub 6 .
The different H functions reported here are H NA = 0, the random policy H Rnd , and the Manhattan distance H Man .
The algorithms were also evaluated using Knuth's algorithm with a different range of rollouts for the cost-to-go estimate, however, the results are not reported here as they are either statistically indifferent or dominated by the results using H Rnd with a single rollout.
BID4 suggests that MCTS algorithms should readily benefit from stronger algorithms to estimate costs-to-go by simulation of stochastic policies.
Our experiments showed that if synergies exist these do not manifest when using off-the-shelf stochastic estimation techniques like the ones discussed by BID19 .
TAB2 , 2 and 3 report the results of the different lookahead algorithms on the GridWorld domain variants with a stationary goal, moving goals and obstacles respectively.
For these three domains, results were also collected for a 100x100 grid, however, the results were omitted from the tables as the simulator budgets used were not sufficient to find anything meaningful.
The results on the stationary goal GridWorld domain shown in TAB2 provide a number of insights about the rollout algorithms reported.
First, we can see RIW(1) benefits from using H Rnd rather than H NA where the larger simulator budgets are used.
As the simulator budget increases, as could be expected, so does the performance of all the methods using H Rnd .
On the contrary, with H NA RIW(1)'s performance remains constant across the different budgets.
The explanation for this can be found in the motivating example we gave previously in this paper with the agent preferring the shorter trajectory of driving into the boundary of the grid.
TAB2 also shows that given the largest budget and H Rnd , RIW(1) statistically outperforms the other algorithms on the three domains of different size.
Table 2 for GridWorld with moving goals has similar patterns as the stationary goal domain in that RIW(1) with H Rnd dominates performance for the largest budget.
Also, the majority of performances for the smaller budgets, excluding RIW (1) 199.0 ± 1.9 193.1 ± 5.0 190.2 ± 6.0 Table 6 : Same settings as TAB6 over Combolock.
GridWorld with a stationary goal and obstacles results displayed in Table 3 continues the trend of results.
Using the largest budget RIW(1) with H Rnd outperforms all methods on the 10x10 and 20x20 domains.
For the 50x50 a number of results are statistically indifferent.
For this domain the algorithms using H Man as the base heuristic are also reported.
While using the largest budget on the 10x10 grid H Rnd dominates H Man , for the larger 50x50 we see H Man dominates H Rnd for UCT, and is competitive for the other methods.
For the smallest simulator budget on CTP reported in Table 4 using H Man with RIW(1) and UCT are the dominate methods.
For the largest simulator budget RIW(1) using H Rnd is dominant over all other methods for both sized domains.
We also see that in most cases for the two smaller budgets H Man dominates the H Rnd methods.
TAB6 and 6 show on the smaller 10 state domains RIW(1) with H Rnd is statistically dominant over all other methods on Antishaping and Combolock for the 500 and 1000 simulator budgets.
However, for the more complex 50 state domains, the results of all algorithms using H Rnd are statistically indifferent.
It can be observed that using H Rnd with RIW(1) does improve its performance compared with H NA across all the domain settings with simulator budgets of 500 and 1000, besides Antishaping with 50 states.
For the Skiing Atari-2600 game results in Table 7 H Heu is the heuristic value based on the number of gates missed and remaining as described in the previous section.
RIW(1) using H Heu dominates all other methods.
Comparing RIW(1) using H Heu results with those reported by BID17 , it has similar performance to the DQN algorithm BID18 after 100 million frames of training.
Since the simulation budget per action we use here is equivalent to 500 frames, and given that the maximum episode duration spans 3,600 actions, RIW(1) achieves the performance in Table 7 considering only 1.8 million frames.
Bertsekas (2017) considers AlphaGo Zero BID20 to be state-of-the-art in MCTS algorithms.
It combines the reasoning over confidence intervals first introduced with UCT BID14 and the classic simulation of base policies BID10 , adding to both supervised learning algorithms to obtain, offline, parametric representations of costs-to-go which are efficient to evaluate.
The resulting algorithm achieves super-human performance at the game of Go, long considered too hard for AI agents.
Rather Table 7 : Average and 95% confidence interval for the cost on the Atari-2600 Skiing game over 100 episodes.than using descriptions of game states directly as we do, AlphaZero uses a CNN to extract automatically features that describe spatial relations between game pieces.
Like us, AlphaZero's lookahead uses a stochastic policy to select what paths to expand, but rather than Q-factors, uses estimated win probabilities to prioritise controls, and simulates the opponent strategy via self-play to generate successor states.
Our simulators are always given and remain unchanged, rather than being dynamic as is the case for AlphaZero.
BID11 have recently presented a hybrid planning and learning approach that integrates BID1 rollout, with a deep neural network.
Similarly to AlphaZero, they use it to both guide the search, and also to extract automatically the set of state features F .
Interestingly, Junyent et al.'s work does not use deep neural networks to approximate costs-to-go as AlphaZero does.
A significant improvement in performance over Bandres et al. original rollout algorithm is reported with policies learnt after 40 million interactions with the simulator, using an overall computational budget much bigger than the one used by Bandres et al.
MCTS approaches typically combine lookaheads and cost-togo approximations, along with statistical tests to determine what are the most promising directions and focus their sampling effort.
The width-based methods described in this paper do so too, but in ways which are, at first sight, orthogonal to existing strategies.
It remains an area of active research to map out exactly how the width-based methods described in this paper, and those elsewhere by BID11 too, provide alternatives to the limitations of existing MCTS approaches.
Having said this, there is no general theory guiding the design of MCTS algorithms BID4 , and to avoid generating ad-hoc, problem dependent solutions involuntarily it is important to follow strict protocols that alert of potential lack of statistical significance in results, while relying on a diverse set of benchmarks that can be both easily understood, and highlight limitations of existing state-of-theart methods and overcome them.

assistant: text: We propose a new Monte Carlo Tree Search / rollout algorithm that relies on width-based search to construct a lookahead.