[
Abstract
This paper describes Picat’s planner, its implementation, and planning models for several domains used in International Planning Competition (IPC) 2014. Picat’s planner is implemented by use of tabling. During search, every state encountered is tabled, and tabled states are used to effectively perform resourcebounded search. In Picat, structured data can be used to avoid enumerating all possible permutations of objects, and term sharing is used to avoid duplication of common state data. This paper presents several modeling techniques through the example models, ranging from designing state representations to facilitate data sharing and symmetry breaking, encoding actions with operations for efficient precondition checking and state updating, to incorporating domain knowledge and heuristics. Broadly, this paper demonstrates the effectiveness of tabled logic programming for planning, and argues the importance of modeling despite recent significant progress in domainindependent PDDL planners.
Planning as Tabled Logic Programming]
Planning as Tabled Logic Programming
N.F. Zhou]
NENGFA ZHOU
CUNY Brooklyn College and Graduate Center
and Roman Barták
Charles University
and Agostino Dovier
Univ. di Udine
\jdateFebruary 2014
2014
\pagerange[–LABEL:lastpage
\doiXXX
1 Introduction
Planning and logic programming are two close areas of research. PLANNER [Hewitt (1969)], which was designed as a language for proving theorems and manipulating models in a robot, is perceived as the first logic programming language. Planning has been an important problem domain for Prolog [Kowalski (1979), Warren (1974)]. Despite the amenability of Prolog to planning, Prolog is no longer a competitive tool for planning.
Tabling [Michie (1968), Tamaki and Sato (1986), Warren (1992)] is a technique used in logic and functional programming systems, which caches the results of certain calculations in memory and reuses them in subsequent calculations through a quick table lookup. Like state marking used in search algorithms, tabling can prevent the same state from being expanded more than once during search. Tabling has been found useful in many search problems, including theorem proving [Nielson et al. (2004), Pientka (2003)], program analysis [Dawson et al. (1996)] and model checking [Ramakrishna et al. (1997)]. Recently, tabled logic programming has been successfully employed to solve specific planning problems [Bartak and Zhou (2014), Zhou and Dovier (2013), Zhou (2014)], and has been shown to be significantly faster than the stateoftheart ASP (Answer Set Programming) planners on some problems [Zhou and Dovier (2013), Zhou (2014)].
This paper describes Picat’s planner and its implementation. This paper also presents planning models in Picat for several domains used in International Planning Competition 2014 (IPC’14) [Chrpa et al. (2014)], and demonstrates broader applicability of tabled logic programming to planning. Picat is a logicbased multiparadigm language that provides logic variables, pattern matching, nondeterminism through backtracking, loops, functions, constraints, and tabling as its core modeling and solving features. As a modeling language for planning, Picat differs from PDDL (Plan Domain Description Language) [McDermott (1998)] and ASP [Brewka et al. (2011), Gebser et al. (2012), Lifschitz (2002)] in several aspects: (1) Picat allows use of structures to represent states; (2) Picat supports explicit commitment and nondeterministic actions, which enables users to have better control over action applications; (3) Picat provides facilities for describing domain knowledge and heuristics for pruning search space.
As a solving system, Picat’s planner implements several techniques for better performance. First, it tables every state encountered during search and avoids repeating the exploration of the same state. Second, it adopts the hashconsing technique [Zhou and Have (2012)] to share common state data and to speed up the equality testing of states. Third, it utilizes tabled states to effectively perform resourcebounded search. For optimal planning, Picat offers builtins to perform iterative search, but unlike IDA* [Korf (1985)], Picat also reuses results tabled in early iterations [Zhou (2014)].
This paper shows that the abovementioned features of Picat make Picat a more appropriate language than PDDL for modeling and solving planning problems. To that end, this paper presents examples in Picat for several domains used in IPC’14. These examples illustrate several modeling techniques on how to design state representations to facilitate data sharing and symmetry breaking, on how to translate PDDL operators into Picat actions, and on how to incorporate domain knowledge and heuristics to reduce search spaces. This paper also gives the experimental results of the presented models and several other models encoded in the same way. The experimental results demonstrate the effectiveness of tabling and the importance of modeling.
2 A Brief Overview of Picat
Picat is a dynamicallytyped language. The basic types are taken from Prolog, except for arrays and maps. An array takes the form {,,}. The index of the first array element is 1, and the index notation X[I] can be used to access array elements. Picat also borrows the basic logical operators from Prolog, including conjunction (), negation (not ), disjunction (;), and ifthenelse (>;).
Picat allows function calls in arguments. For this reason, it requires structures to be preceded with a dollar symbol $ in order for them to be treated as data, unless the structure is special, or it occurs in a head pattern.
For each type, Picat provides a set of builtin functions and predicates. Many builtin predicates are taken from Prolog, including member/2, nth/3, and select/3. The function insert_ordered(,) inserts into the ordered list such that the resulting list remains ordered.
In Picat, predicates and functions are defined with patternmatching rules. Picat has two types of rules: the nonbacktrackable rule =>
, and the backtrackable rule ?=>
. In a predicate definition, the takes the form , where is a predicate name, and is the arity. The condition , which is an optional goal, specifies a condition under which the rule is applicable. For a call , if matches and succeeds, then the rule is said to be applicable to . When applying a rule to call , Picat rewrites into . If the used rule is nonbacktrackable, then the rewriting is a commitment, and the program can never backtrack to . However, if the used rule is backtrackable, then the program will backtrack to once fails, meaning that will be rewritten back to , and the next applicable rule will be tried on . In a function definition, the takes the form where is a function name and is a result to be returned. All of the rules in a function definition must be nonbacktrackable.
A pattern can contain aspatterns of the form @, where is a new variable in the head, and is a nonvariable term. The aspattern @ is the same as in pattern matching, but after pattern matching succeeds, is made to reference the term that matched .
Picat supports loops and list comprehensions. For example, the loop
foreach(E in L) Goal end
is true if Goal is true for each E in L. Picat adopts the following simple scoping rule: variables that occur only in a loop, but do not occur before the loop in the outer scope, are local to each iteration of the loop. Loops are compiled into tailrecursive predicates, and list comprehensions are compiled into tailrecursive predicates through foreach loops.
Picat supports tabling for dynamic programming solutions. Other features of Picat include assignments, list comprehensions, global maps for storing permanent data, higherorder functions, action rules for defining eventdriven actors, and modules for modeling and solving constraint satisfaction problems with CP, SAT, and MIP.
3 Tabling in Picat
Both predicates and functions can be tabled. In order to have all calls and answers of a predicate or function tabled, users just need to add the keyword table before the first rule. For a predicate definition, the keyword table can be followed by a tuple of table modes, including + (input),  (output), min, max, and nt (not tabled). For a predicate with a table mode declaration that contains min or max, Picat tables one optimal answer for each tuple of the input arguments. The last mode can be nt, which indicates that the corresponding argument will not be tabled.
Linear tabling [Zhou et al. (2008)] is used in Picat, and table modes are taken from [Guo and Gupta (2008)], except for the nt mode, which was initially proposed by [Zhou et al. (2010)]. Ground structured terms are hashconsed [Zhou and Have (2012)] so that common ground terms are tabled only once. For example, for the three lists [1,2,3], [2,3], and [3], the shared sublists [2,3] and [3] are reused from [1,2,3].
Modedirected tabling has been successfully used to solve specific planning problems such as Sokoban [Zhou and Dovier (2013)], and the Petrobras planning problem [Bartak and Zhou (2014)]. A planning problem is modeled as a pathfinding problem over an implicitly specified graph. The following gives the framework used in all these solutions.
table (+,,min) path(S,Path,Cost), final(S) => Path = [],Cost = 0. path(S,Path,Cost) => action(S,S1,Action,ActionCost), path(S1,Path1,Cost1), Path = [ActionPath1], Cost = Cost1+ActionCost.
The call path(S,Path,Cost) binds Path to an optimal path from S to a final state. The predicate final(S) succeeds if S is a final state, and the predicate action encodes the set of actions in the problem.
Consider using the path/3 predicate to find a shortest path from node a to node d in the DAG shown in Figure 1. The final state is d. The call path(a,Path,Cost) initiates the search from node a. In order to resolve the call, Picat applies the transition ab to node a, and generates a new call to path for finding a shortest path from b to d. After the answer bcd is found for the call, Picat tables the answer abcd, which has cost 8, for the inital call. After that, Picat backtracks, trying the alternative transition ac to node a. After the shortest path cd is found, Picat finds another path acd, which has cost 5, for the initial call. Since this path is shorter than the tabled path, Picat replaces the tabled path with this new path. Since the graph has no cycles, Picat returns the new path as the final answer. In general, calls need to be reevaluated until fixed points are reached if there are looping calls [Zhou et al. (2008)].
When applied to the singlesource shortest path problem, linear tabling is similar to Dijkstra’s algorithm, except that linear tabling tables shortest paths from the encountered states to the goal state rather than shortest paths to the encountered states from the initial state.
The above framework performs depthunbounded search. For many planning problems, branch & bound and IDA* [Korf (1985)] are useful for finding optimal solutions. The planner module of Picat provides builtins for planning with different types of search.
4 The planner Module of Picat
The planner module is based on tabling but it abstracts away tabling from users. For a planning problem, users only need to define the predicates final/1 and action/4, and call one of the search predicates in the module on an initial state in order to find a plan or an optimal plan.

final(): This predicate succeeds if is a final state.

action(,,,): This predicate encodes the state transition diagram of a planning problem. The state can be transformed to by performing . The cost of is , which must be nonnegative. If the plan’s length is the only interest, then = 1.
These two predicates are called by the planner. The action predicate specifies the precondition, effect, and cost of each of the actions. This predicate is normally defined with nondeterministic patternmatching rules. As in Prolog, the planner tries actions in the order they are specified. When a nonbacktrackable rule is applied to a call, the remaining rules will be discarded for the call.
The following predicates constitute the core of the planner module.

best_plan_unbounded(,,,): This predicate finds an optimal plan by performing depthunbounded search. This predicate is implemented based on the pathfinding framework shown above. The argument is not utilized to limit the depth of search. It is compared with after an optimal plan has been found. During depthunbounded search, once a state has failed, it will not be explored again.

plan(,,,): This predicate searches for a plan by performing resourcebounded search, in which a state is expanded only if it is new and its resource limit is nonnegative, or if the state has previously failed but the current occurrence has a higher resource limit than before. This predicate is defined as a tabled predicate that tables the argument but does not use it in variant checking [Zhou (2014)]. The implementation of this predicate is described in the next subsection.

best_plan(,,,): This predicate finds an optimal plan by performing resourcebounded iterativedeepening search. It calls the plan/4 predicate to find a plan, using 0 as the initial cost limit and gradually relaxing the cost limit until a plan is found. Unlike IDA*, which starts a new round from scratch, Picat also reuses the states that were tabled in the previous rounds.

best_plan_bb(,,,) This predicate finds an optimal plan using branch & bound. First, it calls plan/4 to find a plan. Then, it tries to find a better plan by imposing a stricter limit. This step is repeated until no better plan can be found. It returns the last plan that was found.

current_resource(): This function returns the resource limit argument of the latest call to plan/4. In order to retrieve the argument, the implementation has to traverse the callstack until it reaches a call to plan/4. This function can be used to check against a heuristic value. If the heuristic estimate of the cost to travel from the current state to a final state is greater than the resource limit, then the current state should fail.
5 The Implementation of ResourceBounded Search
Figure 2 sketches Picat’s implementation of the predicate plan/4. The following array is passed from a state to the next state as an nt argument:
{Limit,IPlan,IPlanCost}
In order to perform resourcebounded search, Picat treats the second argument of the predicate plan_bounded_aux differently from other nt arguments. Picat stores the Limit argument of each failed call to plan_bounded_aux, and uses this information to decide whether the same state should fail when it recurs.
The first rule of plan_bounded_aux throws the inherited plan and its cost as an exception if final(S) succeeds. The exception will be caught by the catch call in the rule body of plan/4.
The second rule of plan_bounded_aux calls action/4 to select an action, which produces a new state, NextS. Then, the rule computes the new resource limit, Limit1, by subtracting the cost of the selected action from Limit. If Limit1 >= 0 succeeds, then the rule continues with the tabled search by recursively calling plan_bounded_aux on the new state NextS. Otherwise, if Limit1 >= 0 fails, then Picat backtracks to select an alternative action.
The idea of resourcebounded search is to utilize tabled states and their resource limits to effectively decide when a state should be expanded and when a state should fail. Let denote a state with an associated resource limit, , as depicted in Figure 3. If is negative, then immediately fails. If is nonnegative and has never been encountered before, then is expanded by using a selected action. Otherwise, if the same state has failed before and was the resource limit when it failed, then is only expanded if , i.e., if the current resource limit is larger than the resource limit was at the time of failure.
6 Modeling Examples
This section describes the Picat models for four of the problem domains used in the sequential optimal track of IPC’14: Transport, Tetris, Floortile, and Parking. Each of the following models will show the state representation, the encoding of the actions, the search predicate that is used (best_plan or best_plan_unbounded), and the domain knowledge and heuristics that are employed. In these models, states are typically represented by lists, and preconditions and state updates are handled by standard list operations. Sometimes, arrays are used when no list suffixes can be shared. The Picat encodings of actions are mostly straightforward translations from the PDDL encodings. These models do not use sophisticated domain knowledge or heuristics that would hurt the readability or compromise the optimality of answers. These four domains, as well as five other domains, are described in Appendix A.
6.1 Example1: Transport
This problem is a variant of the popular logistics domain in planning. Given a weighted directed graph, a set of trucks each of which has a capacity for the number of packages it can carry, and a set of packages each of which has an initial location and a destination, the objective of the problem is to find an optimal plan to transport the packages from their initial locations to their destinations. This problem is more challenging than the Nomystery problem that was used in IPC’11, because of the existence of multiple trucks, and because an optimal plan normally requires trucks to cooperate. This problem degenerates into the shortest path problem if there is only one truck and only one package. The PDDL and Picat encodings of the problem are given in Appendix B.
Basic Encoding
A state is represented by an array of the form {Trucks,Packages}, where Trucks is an ordered list of trucks, and Packages is an ordered list of waiting packages. A package in Packages is a pair of the form (Loc,Dest) where Loc is the source location and Dest is the destination of the package. A truck in Trucks is a list of the form [Loc,Dests,Cap], where Loc is the current location of the truck, Dests is an ordered list of destinations of the loaded packages on the truck, and Cap is the capacity of the truck. At any time, the number of loaded packages must not exceed the capacity.
Note that keeping Cap as the last element of the list facilitates sharing, since the suffix [Cap], which is common to all the trucks that have the same capacity, is tabled only once. Also note that the names of the trucks and the names of packages are not included in the representation. Two packages in the waiting list that have the same source and the same destination are indistinguishable, and as are two packages loaded on the same truck that have the same destination. This representation breaks symmetries. Two configurations that only differ by a truck’s name or a package’s name are treated as the same state.
A state is final if all of the packages have been transported.
final({Trucks,[]}) => foreach([_Loc,Dests_] in Trucks) Dests == [] end.
The PDDL rules for the actions are straightforwardly translated into Picat as follows.
action({Trucks,Packages},NextState,Action,ACost) ?=> Action = $load(Loc), ACost = 1, select([Loc,Dests,Cap],Trucks,TrucksR), length(Dests) < Cap, select((Loc,Dest),Packages,PackagesR), NewDests = insert_ordered(Dests,Dest), NewTrucks = insert_ordered(TrucksR,[Loc,NewDests,Cap]), NextState = {NewTrucks,PackagesR}, action({Trucks,Packages},NextState,Action,ACost) ?=> Action = $unload(Loc), ACost = 1, select([Loc,Dests,Cap],Trucks,TrucksR), select(Dest,Dests,DestsR), NewTrucks = insert_ordered(TrucksR,[Loc,DestsR,Cap]), NewPackages = insert_ordered(Packages,(Loc,Dest)), NextState = {NewTrucks,NewPackages}. action({Trucks,Packages},NextState,Action,ACost) => Action = $move(Loc,NextLoc), select([LocTail],Trucks,TrucksR), road(Loc,NextLoc,ACost), NewTrucks = insert_ordered(TrucksR,[NextLocTail]), NextState = {NewTrucks,Packages}.
For the load action, the rule nondeterministically selects a truck that still has room for another package, and nondeterministically selects a package that has the same location as the truck. After loading the package to the truck, the rule inserts the package’s destination into the list of loaded packages of the truck. Note that the rule is nondeterministic. Even if a truck passes by a location that has a waiting package, the truck may not pick it. If this rule is made deterministic, then the optimality of plans is no longer guaranteed, unless there is only one truck and the truck’s capacity is infinite.
Domain Knowledge and Heuristics
Domain knowledge can be used to reduce the nondeterminism and avoid unnecessary applications of actions. In the complete Picat encoding given in Appendix B, the predicate action begins with a rule that deterministically unloads a package if the package’s destination is the same as the truck’s location.
Resourcebounded search is used to find an optimal plan. After each new state is generated, the following condition is checked to ensure that the current path is viable.
current_resource()  ACost >= estimated_cost(NewState).
Let , , be the remaining packages, and be the minimum cost of moving package to its destination by using any truck. The moving cost of a state can be safely estimated as max(). The estimated total cost is the estimated moving cost plus the loading and unloading costs of all of the remaining packages. This heuristic is admissible.
6.2 Example2: Tetris
The problem is a simplified version of the wellknown Tetris. There are three kinds of pieces: 1cell boxes, 2cell rectangles, and 3cell Lshaped pieces. Initially, all the pieces are distributed on a grid board. The pieces can freely move on the board to any cells as long as the cells are not occupied by other pieces. There is a region on the board that is designated as the target region. The goal of the game is to move all the pieces to the target region. In the IPC setting, rectangles can rotate but Lshaped pieces cannot.
Basic Encoding
Pieces can be represented as follows: A square is represented by the cell it occupies; a rectangle by a term of the form rect(C1,C2) where C1 and C2 are cell locations; and an Lshaped piece by a term of the form ell(C1,C2,C3). In the IPC setting, the ordering of the cells that are occupied by a piece is important. So rect(C1,C2) and rect(C2,C1) represent two different rectangle pieces. A state is represented by an array of the form {Squares,Rects,Ls}, where each argument gives a sorted list of a single kind of pieces.
The PDDL rules for actions can be translated into Picat in a straightforward manner. For example, the following rule selects a rectangle piece to move.
action(State@{Squares,Rects,Ls},NewState,Action,ACost) ?=> Action = $move(Rect,NewRect), ACost = 1, select(Rect,Rects,RectsR), Rect = $rect(C1,C2), connected(C2,C3), not_occupied(C3,State), %% C3 is free NewRect = $rect(C2,C3), NewState = {Squares,insert_ordered(Rects,NewRect),Ls}.
The predicate not_occupied(C3,State) is true if C3 is not occupied by any of the pieces in State.
Heuristics
Resourcebounded search is used to find an optimal plan for this problem. The goal of the problem is to move all the pieces to the target region. For each piece, assume that it’s the only piece on the board, disregarding all other pieces. The estimated cost of moving the piece to the target region is its shortest distance to the nearest cell in the target region times the cost of each move. The estimated cost of transforming a state into a final state is the sum of the estimated costs of all the pieces. This heuristic function is admissible and easy to compute, but very conservative because a cell in the target region can serve as a target for multiple pieces.
6.3 Example3: Floortile
A set of robots use two different colors (black and white) to paint patterns in floor tiles. The robots can move around the floor tiles in four directions (up, down, left and right). Robots paint with one color at a time, but can change their spray guns to the other color. However, robots can only paint the tile that is in front (up) and behind (down) them, and once a tile has been painted no robot can stand on it.
Basic Encoding
A state is represented by a list of the form [Robots,WTiles,BTiles], where Robots is an ordered list of robots and WTitles (BTiles) is an ordered list of locations of the painted white (black) tiles. A robot in Robots is a pair (Color,Loc) where Color is the color of the paint that the robot is holding and Loc is the robot’s location.
A state is final if all the tiles that are required to be painted are all painted.
final({_,WTiles,BTiles}) => painted_w_tiles_in_goal(WTiles), painted_b_tiles_in_goal(BTiles).
The actions are encoded according to the state representation. For example, once a tile is painted, the tile’s location is added into WTiles or BTiles depending on the color of the paint; once a robot (Color,Loc) moves from Loc to Loc1, the pair is changed to (Color,Loc1). The current graph is determined by the initial graph, the set of painted tiles, and the set of robots. A robot can move from its current location Loc to a new location NextLoc if NextLoc is connected to Loc in the graph, NextLoc is not painted, and NextLoc is not occupied by another robot.
Macro Actions and Domain Knowledge
The colorchanging action, if needed, can be forced to take place right before a painting action. This macro action helps reduce the search space.
For IPC’14, the used instances were generated in the fashion that robots should only paint tiles in front of them. This special condition can be exploited to reduce the nondeterminism of the painting rules: if a robot is at location Loc, trying to paint the up location ULoc, and the up location of ULoc has already been painted, then ULoc can be painted deterministically; similarly, if a robot is at location Loc, trying to paint the down location DLoc, and the down location of DLoc has already been painted, then DLoc can be painted deterministically.
Depthunbounded search is used for this problem. During depthunbounded search, dead ends are tabled and are not reexplored when they are encountered again. For this problem, a state becomes a dead end if there is an unpainted tile that is not reachable by any robot. No extra code is needed to detect dead ends. Because depthunbounded search is used, no heuristics are needed.
6.4 Example4: Parking
This domain involves parking cars on a street with N curb locations, and where cars can be doubleparked but not tripleparked. The goal is to find a plan to move from one configuration of parked cars to another configuration, by driving cars from one curb location to another. For each curb location, there are two parking spots: the curb side and the road side. For a car to move from a spot at a curb into a spot at a different curb, the destination spot must be clear; and if the spot is on the curb side, then the roadside spot must also be clear. In some ways, this problem is similar to the Blocks World and the Tower of Hanoi.
Basic Encoding
A configuration is represented by an array of curbs , where each is a list of up to two cars. When a curb has two cars [,], it is assumed that is the roadside car and is the curbside car. This representation allows the roadside car to be removed without touching the curbside car .
A state is represented by a pair of the form (CurrConfig,GoalConfig), where a CurrConfig is the current configuration, and GoalConfig is the goal configuration. As shown below, the inclusion of the goal configuration in the state representation facilitates the encoding of domain knowledge.
A state is final if the current configuration is the same as the goal configuration.
final((Config,Config)) => true.
All possible moves are specified with one rule as follows:
action((Config,GConfig),NextS,Action,ACost) => Action = $move(I,J), ACost = 1, nth(I,Config,ICars),ICars = [ClearCarNewICars], nth(J,Config,JCars), J !== I, JCars != [_,_], N = length(Config), NewConfig = new_array(N), foreach (K in 1..N) (K == I > NewConfig[K] = NewICars ;K == J > NewConfig[K] = [ClearCarJCars] ; NewConfig[K] = Config[K]) end, NextS = (NextConfig,GConfig).
The builtin predicate nth(I,Array,Arg) is true if the Ith argument of Array is Arg. When I is a variable, this predicate nondeterministically searches for an argument that unifies with Arg. The rule finds a nonempty curb I and a different curb J that has less than two cars, and moves the clear car, ClearCar, from curb I to curb J.
Domain Knowledge and Heuristics
Two knowledge rules can be incorporated into the basic encoding in order to speed up the search. First, when a car is moved into a spot that is its final spot in the goal configuration, such a move should be made deterministically. Note that if a car is moved to a roadside final spot then the curbside spot must be occupied by a correct car. Second, a car that has been placed in its final spot should not be moved away. This kind of domain knowledge has been used in solving the Tower of Hanoi [Alford et al. (2009)] and the Blocks World [Bacchus and Kabanza (2000)].
Resourcebounded search is used to find an optimal plan. The cost of transforming a state to a final state can be simply estimated as the number of incorrectlypositioned cars. This estimation is admissible, since at least one action is needed to move each incorrectlypositioned car. The following improved heuristic function, which takes special curb configurations into account, is used by the model: For each curb, if the curb has two cars [A,B] in the current state, but is required to have [B,A] in the final state, then its cost is 4; if the curb has two cars [A,B], but is required to have [C,A] or [B,C] (A C, B C), then its cost is 3; otherwise, the cost is the number of incorrectlypositioned cars.
7 Experimental Results
In addition to the four domains presented in this paper, we have encoded in Picat several other domains used in the deterministic sequential track of IPC’14. The domains and their Picat encodings are given in Appendix A. The Picat encodings are simple, compact, and comparable in size with the PDDL encodings used in IPC’14. Most of the encodings can be further improved by incorporating sophisticated domain knowledge and heuristics. In order to evaluate the effectiveness of the use of tabling and the use of heuristics, we have built two separate sets of encodings, namely, Picatnt and Picatnh, which use the same state representation as the original Picat encodings but have some component removed. Picatnt performs Prologstyle nontabled iterativedeepening search. Since Picatnt does not table any state, it may explore the same state multiple times during search. The Picatnh encodings do not use any heuristics. We have compared these Picat encodings with the IPC’14 PDDL encodings solved with Symba [Torralba et al. (2014)], a domainindependent bidirectional A* planner which won the optimal sequential track of IPC’14. This comparison offers a glimpse of how well Picat compares with the best domainindependent planner. A comparison of Picat’s planner and several domaindependent planners also shows the promise of tabled planning [Bartak et al. (2015)].
Domain  # insts  Picat  Picatnt  Picatnh  Symba 

Barman  14  14  0  14  6 
Cave  20  20  0  20  3 
Childsnack  20  20  20  20  3 
Citycar  20  20  17  18  17 
Floortile  20  20  0  20  20 
GED  20  20  19  13  19 
Parking  20  11  4  0  1 
Tetris  17  13  13  9  10 
Transport  20  9  0  4  8 
Table 1 shows the number of instances (#insts) in the domains used in IPC’14 and the number of (optimally) solved instances by each planner. The results were obtained on a Cygwin notebook computer with 2.4GHz Intel i5 and 4GB RAM. Both Picat and Symba were compiled using g++ version 4.8.3. For Symba, a setting suggested by one of Symba’s developers was used. A time limit of 30 minutes was used for each instance as in IPC. For every instance that was solved by both Symba and Picat, the plan quality is the same.
A comparison of Picat and Picatnt shows the effectiveness of the use of tabling. For every domain, except for Childsnack and Tetris, Picat solved more instances than Picatnt. For Barman, Cave, Floortile and Transport, Picatnt could not solve any of the instances. The Picat encodings for five of the domains (Citycar, GED, Parking, Tetris, and Transport) use heuristics. The use of heuristics is helpful for these domains, especially for Parking, for which Picatnh did not solve any of the instances.
Picat solved more instances than Symba for every domain except for Floortile, for which both systems solved all of the instances. The running times of the instances are not given, but the total runs for Picat were finished within 24 hours, while the total runs for Symba took more than 72 hours.
8 Conclusion and Future Work
This paper has presented Picat’s planner, its implementation, and example models for several domains from IPC’14. The example models illustrate several modeling techniques in Picat. One key task of modeling is finding an efficient state representation. While classical planning frameworks such as PDDL are based on a factored representation of states, Picat uses a structured representation. A structured state representation can leave out unnecessary information that is not needed for planning and can break symmetries by avoiding enumerating all possible permutations of objects. Another key task of modeling is utilizing domain knowledge to reduce search spaces. In the past, a lot of work has been done on the use of domain knowledge in planning [Bacchus and Kabanza (2000), Haslum and Scholz (2003), Kautz and Selman (1998)], but recently, this part of modeling has been put aside, because of the advancement of domainindependent PDDL planners. This paper has shown that, even with simple domain knowledge, the declarative encodings of Picat significantly outperform Symba, a stateoftheart domainindependent planner.
This paper has demonstrated for the first time that tabled logic programming is competitive with the cuttingedge PDDL planners. The key to the success is tabling. Tabling avoids repeating the exploration of the same state and facilitates performing resourcebounded search. The Picat planner does not do prior grounding that is typical for most current statespace search planners, and hence Picat has no problem with exploded memory consumption due to grounding. Nevertheless, memory consumption can be demanding during search since every encountered state is tabled. This is why careful modeling that removes symmetries and uses domain control knowledge to prune useless state transitions is important.
A good state representation should also exploit the underlying termsharing technique that is used in the tabling system. In the examples that are presented in this paper, ordered lists are used to represent collections of objects. It takes linear time to perform the basic operations. One direction for future work is to design data structures for state representations which are compact, efficient, and good for sharing.
A plethora of action languages have been designed for modeling and solving planning problems (e.g., A and its successors [Gelfond and Lifschitz (1998)], Golog [Levesque et al. (1997)], and K [Eiter et al. (2004)]). The focus of these languages has been on the modeling power rather than efficiency and scalability. These languages have been implemented by translation into SAT, ASP, CP, or PDDL [Baier et al. (2011), Dovier et al. (2011)], but no implementation has been shown to be competitive with the cuttingedge PDDL planners. Picat can be used as an implementation language for these action languages. Like in PDDL, a state is represented as a set of flat facts in all of these action languages. Another direction for future work is to devise an efficient translation from these action languages into Picat that automatically exploits structural representation, symmetries, domain control knowledge, and heuristics.
References
 Alford et al. (2009) Alford, R., Kuter, U., and Nau, D. S. 2009. Translating HTNs to PDDL: A small amount of domain knowledge can go a long way. In IJCAI. 1629–1634.
 Bacchus and Kabanza (2000) Bacchus, F. and Kabanza, F. 2000. Using temporal logics to express search control knowledge for planning. Artif. Intell. 116, 12, 123–191.
 Baier et al. (2011) Baier, J. A., Fritz, C., and McIlraith, S. A. 2011. Gologstyle search control for planning. In Knowing, Reasoning, and Acting: Essays in Honour of Hector J. Levesque, G. Lakemeyer and S. A. McIlraith, Eds. College Publications.
 Bartak et al. (2015) Bartak, R., Dovier, A., and Zhou, N.F. 2015. On modeling planning problems in logic programming. In PPDP. to appear.
 Bartak and Zhou (2014) Bartak, R. and Zhou, N. 2014. Using tabled logic programming to solve the Petrobras planning problem. TPLP 14, 45, 697–710.
 Brewka et al. (2011) Brewka, G., Eiter, T., and Truszczyński, M. 2011. Answer set programming at a glance. Commun. ACM 54, 12, 92–103.
 Chrpa et al. (2014) Chrpa, L., Vallati, M., and McCluskey, L. 2014. International planning competition.
 Dawson et al. (1996) Dawson, S., Ramakrishnan, C. R., and Warren, D. S. 1996. Practical program analysis using general purpose logic programming systems — A case study. ACM SIGPLAN Notices 31, 5, 117–126.
 Dovier et al. (2011) Dovier, A., Formisano, A., and Pontelli, E. 2011. Perspectives on logicbased approaches for reasoning about actions and change. In LNCS. Vol. 6565. 259–279.
 Eiter et al. (2004) Eiter, T., Faber, W., Leone, N., Pfeifer, G., and Polleres, A. 2004. A logic programming approach to knowledgestate planning: Semantics and complexity. ACM Trans. Comput. Log. 5, 2, 206–263.
 Gebser et al. (2012) Gebser, M., Kaminski, R., Kaufmann, B., and Schaub, T. 2012. Answer Set Solving in Practice. Morgan and Claypool Publishers.
 Gelfond and Lifschitz (1998) Gelfond, M. and Lifschitz, V. 1998. Action languages. Electron. Trans. Artif. Intell. 2, 193–210.
 Guo and Gupta (2008) Guo, H.F. and Gupta, G. 2008. Simplifying dynamic programming via modedirected tabling. Softw., Pract. Exper. 38, 1, 75–94.
 Haslum and Scholz (2003) Haslum, P. and Scholz, U. 2003. Domain knowledge in planning: Representation and use. In ICAPS Workshop on PDDL.
 Hewitt (1969) Hewitt, C. 1969. Planner: A language for proving theorems in robots. In IJCAI. 295–302.
 Kautz and Selman (1998) Kautz, H. and Selman, B. 1998. The role of domainspecific knowledge in the planning as satisfiability framework. In AIPS98. 181–189.
 Korf (1985) Korf, R. E. 1985. Depthfirst iterativedeepening: An optimal admissible tree search. Artif. Intell. 27, 1, 97–109.
 Kowalski (1979) Kowalski, R. 1979. Logic for Problem Solving. North Holland, Elsevier.
 Levesque et al. (1997) Levesque, H. J., Reiter, R., Lespérance, Y., Lin, F., and Scherl, R. B. 1997. GOLOG: A logic programming language for dynamic domains. J. Log. Program. 31, 13, 59–83.
 Lifschitz (2002) Lifschitz, V. 2002. Answer set programming and plan generation. Artif. Intell. 138, 12, 39–54.
 McDermott (1998) McDermott, D. 1998. The planning domain definition language manual. CVC Report 98003, Yale Computer Science Report 1165.
 Michie (1968) Michie, D. 1968. “memo” functions and machine learning. Nature, 19–22.
 Nielson et al. (2004) Nielson, F., Nielson, H. R., Sun, H., Buchholtz, M., Hansen, R. R., Pilegaard, H., and Seidl, H. 2004. The succinct solver suite. In Proc. Tools and Algorithms for the Construction and Analysis of Systems: 10th International Conference (TACAS), LNCS 2988. 251–265.
 Pientka (2003) Pientka, B. December 2003. Tabled higherorder logic programming. Ph.D. thesis, Technical Report CMUCS03185.
 Ramakrishna et al. (1997) Ramakrishna, Y. S., Ramakrishnan, C. R., Ramakrishnan, I. V., Smolka, S. A., Swift, T., and Warren, D. S. 1997. Efficient model checking using tabled resolution. In Computer Aided Verification. 143–154.
 Tamaki and Sato (1986) Tamaki, H. and Sato, T. 1986. OLD resolution with tabulation. In ICLP. 84–98.
 Torralba et al. (2014) Torralba, A., Alcazar, V., and Borrajo, D. 2014. Symbaâ: A symbolic bidirectional aâ planner. In The 2014 International Planning Competition. 105–109.
 Warren (1974) Warren, D. H. D. 1974. WARPLAN: A system for generating plans. Tech. Rep. DCL Memo 76, University of Edinburgh.
 Warren (1992) Warren, D. S. 1992. Memoing for logic programs. Comm. of the ACM, Special Section on Logic Programming 35, 93–111.
 Zhou (2014) Zhou, N.F. 2014. Combinatorial search with Picat. ICLP, invited talk, http://arxiv.org/abs/1405.2538.
 Zhou and Dovier (2013) Zhou, N.F. and Dovier, A. 2013. A tabled Prolog program for solving Sokoban. Fundam. Inform. 124, 4, 561–575.
 Zhou and Have (2012) Zhou, N.F. and Have, C. T. 2012. Efficient tabling of structured data with enhanced hashconsing. TPLP 12, 45, 547–563.
 Zhou et al. (2010) Zhou, N.F., Kameya, Y., and Sato, T. 2010. Modedirected tabling for dynamic programming, machine learning, and constraint solving. In ICTAI. 213–218.
 Zhou et al. (2008) Zhou, N.F., Sato, T., and Shen, Y.D. 2008. Linear tabling strategies and optimizations. TPLP 8, 1, 81–109.
Appendix A Benchmarks used in the paper
In this section we summarize, for reader’s convenience, the descriptions of all the domains used as benchmarks. Descriptions are drawn from https://helios.hud.ac.uk/scommv/IPC14/domains_sequential.html; Picat’s complete encodings for these benchmarks are available at http://picatlang.org/ipc14/.
a.1 Barman
There is a robot barman that manipulates drink dispensers, glasses, and a shaker. The goal is to find a plan of robot’s actions that serves a desired set of drinks. Robot hands can grasp at most one object at a time. Glasses need to be empty and clean to be filled. The benchmark was proposed by Sergio Jiménez Celorrio.


In Figure 4 we represent the initial configuration and the corresponding input specifications.
Actions available are:

grasp(OBJ) that executes the grasping either of a specific shot or shaker (OBJ)

leave(OBJ) that allows us to leave the shot or shaker (OBJ)

fill_shot(SHOT,ING) that allows us to fill the shot SHOT with the ingredient ING

empty_shot(SHOT) (resp., empty_shaker(SHAKER)) that allows us to empty the shot SHOT (resp., the skaker SHAKER)

clean_shot(SHOT) (resp., clean_shaker(SHAKER)) that allows us to clean the shot SHOT (resp., the skaker SHAKER)

pour_shot_to_shaker(SHOT,SHAKER) (resp., pour_shaker_to_shot(SHAKER,SHOT)) that allows us to pour the content of the shot SHOT in the shaker SHAKER (resp., vice versa).

shake(SHAKER) that executes that shaking of the shaker to mix the ingredients.

reduce (remove a shot from the state once it contains a required cocktail
All actions have cost 1 but reduce that has cost 0.
a.2 Cave Diving
There is a set of divers, each of who can carry four tanks of air. These divers must be hired to go into an underwater cave and either take photos or prepare the way for other divers by dropping full tanks of air. The cave is too narrow for more than one diver to enter at a time. Divers have a single point of entry. Certain rooms of the cave branches are objectives that the divers must photograph. Swimming and photographing both consume air tanks. Divers must exit the cave and decompress at the end. They can therefore only make a single trip into the cave. Certain divers have no confidence in other divers and will refuse to work if someone they have no confidence in has already worked. Divers have hiring costs inversely proportional to how hard they are to work with. This domain was proposed by Nathan Robinson, Christian Muise, and Charles Gretton.
The cave system is represented by an undirected acyclic graph. Divers can carry an amount of tanks according to their capacity. Rooms that need to be reached are among the leaves of the graph. In Figure 5 we represent an instance of the problem.


The actions available are:

hire_diver(Diver) that requires the availability of hiring cost and should satisfy the compatibility constraints among divers,

prepare_tank(T) that prepares the tank T for the current diver if his capacity allows it,

enter_water that requires the diver to be in the cave entrance,

photograph(Loc) that requires the diver to be in the target location Loc,

drop_tank(Loc) that allows the diver to leave a tank in the location Loc (the tank can be either full or empty),

swim(Loc1,Loc2) that allows the diver to swim between two locations that are adjacent in the graph,

pickup_tank(Loc) that allows the diver to collect a tank stored in the location Loc,

decompress should be made at the end of diving in the cave entrance.
Each action swim and photograph consumes (empties) one air tank. All actions but the first one have unitary cost.
a.3 Childsnack
This domain is to plan how to make and serve sandwiches for a group of children in which some are allergic to gluten. There are two actions for making sandwiches from their ingredients. The first one makes a sandwich and the second one makes a sandwich taking into account that all ingredients are glutenfree. There are also actions to put a sandwich on a tray and to serve sandwiches. Problems in this domain define the ingredients to make sandwiches at the initial state. Goals consist of having selected kids served with a sandwich to which they are not allergic. This domain was proposed by Raquel Fuentetaja, Tomàs de la Rosa Turbides.


Available actions are the following:

make_sandwich_no_gluten(Sw,B,Co) and make_sandwich(Sw,B,Co) where SW is a sandwich, B is a (nogluten) bread, and Co is a (nogluten) content allows us to make the sandwiches.

put_on_tray(Sw,T) puts the sandwich Sw on the tray T

serve_sandwich_no_gluten(Sw,Ch,T,Loc) and serve_sandwich(Sw,Ch,T,Loc) serves the (nogluten) sandwich Sw which is on tray T to the children Ch at the location Loc

move_tray(T,Loc1,Loc2), where T is a tray, Loc1 and Loc2 is a location (i.e., a table, the kitchen)
Each action has cost 1. make_sandwich (nogluten) consumes ingredients.
a.4 Citycar
This model aims to simulate the impact of road building/demolition on traffic flows. A city is represented as an acyclic graph, in which each node is a junction and edges are “potential” roads. Some cars start from different positions and have to reach their final destination as soon as possible. The agent has a finite number of roads available, which can be built for connecting two junctions and allowing a car to move between them. Roads can also be removed, and placed somewhere else, if needed. In order to place roads or to move cars, the destination junction must be clear, i.e., no cars should be in there. The domain was proposed by Mauro Vallati.


Allowed actions are the following:

car_arrived(Dest), which has cost 0. It allows to remove a car from the network and to remove the occurrence of the destination Dest (a junction) from the list of all final destinations.

car_start(Loc): A car is put in the road from the garage of location Loc: it has cost 1.

move_car_in_road(FromLoc) allows us to move a car in a road from the junction FromLoc (cost 1—the road is a straight line or a diagonal road starting in FromLoc).

move_car_out_road(ToLoc) allows us to move a out of a road as soon as the junction ToLoc is reached by the car (cost 1—the road is a straight line or a diagonal road ending in ToLoc).

These actions allow us to build diagonal, straight roads or of deleting one road:

build_diagonal_oneway(FromLoc,ToLoc) (cost 30),

build_straight_oneway(FromLoc,ToLoc) (cost 20),

destroy_road(FromLoc,ToLoc) (cost 10).

Let us observe that search symmetries are eliminated by considering the cars equivalent during the search. It is trivial to label them aposteriori given a correct plan.
a.5 Ged
The GED problem is to find a mincost sequence of operations that transforms one genome (signed permutation of genes) into another. The purpose of this is to use this cost as a measure of the distance between the two genomes, which is used to construct hypotheses about the evolutionary relationship between the organisms. The domains was proposed by Patrik Haslum.
This problem can be stated at several abstraction levels. A general version could include gene insertions and deletions. Let us focus on the abstraction level and on the three rules required by the competition benchmarks.
A gene is identified by a symbolic name. The connection between genes is stated by a binary predicate cw that encodes a linear graph. Each gene can occur in a regular direction (normal) or in reverse direction (inverted).
The three rules allowed are cut (of a substring) from the main genome, and then a splice of the cut substring directed or reversed in a selected point of the main genome. The reverse of a single gene is also allowed. Just to fix the ideas, let us consider the example in figure 8. Reversed genes are overlined.



Each complex action (cut and splice) is split in some subactions as done by Patrik Haslum in his PDDL encoding (http://picatlang.org/ipc14/ged.pddl).
a.6 Floortile, Parking, and Tetris
For the three domains discussed extensively in the core of paper we only show here an instance both in concrete form and as a picture (see Figures 9–11). The Transport domain is discussed in detail in the next section.






Appendix B The Transport Domain
b.1 PDDL Encoding of the Transport Domain
(define (domain transport) (:requirements :typing :actioncosts) (:types location target locatable  object vehicle package  locatable capacitynumber  object ) (:predicates (road ?l1 ?l2  location) (at ?x  locatable ?v  location) (in ?x  package ?v  vehicle) (capacity ?v  vehicle ?s1  capacitynumber) (capacitypredecessor ?s1 ?s2  capacitynumber) ) (:functions (roadlength ?l1 ?l2  location)  number (totalcost)  number ) (:action drive :parameters (?v  vehicle ?l1 ?l2  location) :precondition (and (at ?v ?l1) (road ?l1 ?l2) ) :effect (and (not (at ?v ?l1)) (at ?v ?l2) (increase (totalcost) (roadlength ?l1 ?l2)) ) ) (:action pickup :parameters (?v  vehicle ?l  location ?p  package ?s1 ?s2  capacitynumber) :precondition (and (at ?v ?l) (at ?p ?l) (capacitypredecessor ?s1 ?s2) (capacity ?v ?s2) ) :effect (and (not (at ?p ?l)) (in ?p ?v) (capacity ?v ?s1) (not (capacity ?v ?s2)) (increase (totalcost) 1) ) ) (:action drop :parameters (?v  vehicle ?l  location ?p  package ?s1 ?s2  capacitynumber) :precondition (and (at ?v ?l) (in ?p ?v) (capacitypredecessor ?s1 ?s2) (capacity ?v ?s1) ) :effect (and (not (in ?p ?v)) (at ?p ?l) (capacity ?v ?s2) (not (capacity ?v ?s1)) (increase (totalcost) 1) ) ) )
b.2 Picat Encoding of the Transport Domain
final({Trucks,[]}) => % no waiting packages and no loaded packages foreach([_Loc,Dests_] in Trucks) Dests == [] end. % unload a package action({Trucks,Packages},NextState,Action,ActionCost), select([Loc,Dests,Cap],Trucks,TrucksR), select(Loc,Dests,DestsR) % unload it deterministically => Action = $unload(Loc), ActionCost = 1, NewTrucks = insert_ordered(TrucksR,[Loc,DestsR,Cap]), NextState = {NewTrucks,Packages}. action({Trucks,Packages},NextState,Action,ActionCost) ?=> Action = $unload(Loc), ActionCost = 1, select([Loc,Dests,Cap],Trucks,TrucksR), select(Dest,Dests,DestsR), NewTrucks = insert_ordered(TrucksR,[Loc,DestsR,Cap]), NewPackages = insert_ordered(Packages,(Loc,Dest)), NextState = {NewTrucks,NewPackages}. % load a package onto a truck if the truck and the package are at the same location action({Trucks,Packages},NextState,Action,ActionCost) ?=> Action = $load(Loc), ActionCost = 1, select([Loc,Dests,Cap],Trucks,TrucksR), length(Dests) < Cap, select((Loc,Dest),Packages,PackagesR), % the package is at the same location as the truck NewTrucks = insert_ordered(TrucksR,[Loc,insert_ordered(Dests,Dest),Cap]), NextState = {NewTrucks,PackagesR}. % drive a truck from Loc to NextLoc action({Trucks,Packages},NextState,Action,ActionCost) => Action = $move(Loc,NextLoc), select([LocTail],Trucks,TrucksR), road(Loc,NextLoc,ActionCost), NewTrucks = insert_ordered(TrucksR,[NextLocTail]), NextState = {NewTrucks,Packages}, estimate_cost(NextState) =< current_resource()ActionCost. table estimate_cost({Trucks,Packages}) = Cost => LoadedPackages = [(Loc,Dest) : [Loc,Dests,_] in Trucks, Dest in Dests], NumLoadedPackages = length(LoadedPackages), TruckLocs = [Loc : [Loc_] in Trucks], travel_cost(TruckLocs,LoadedPackages,Packages,0,TCost), Cost = TCost+NumLoadedPackages+length(Packages)*2. % includes load and unload costs % the maximum of the minimum cost of transporting each single package travel_cost(_Trucks,[],[],Cost0,Cost) => Cost=Cost0. travel_cost(Trucks,[(PLoc,PDest)Packages],Packages2,Cost0,Cost) => Cost1 = min([D1+D2 : TLoc in Trucks, shortest_dist(TLoc,PLoc,D1), shortest_dist(PLoc,PDest,D2)]), travel_cost(Trucks,Packages,Packages2,max(Cost0,Cost1),Cost). travel_cost(Trucks,[],Packages2,Cost0,Cost) => travel_cost(Trucks,Packages2,[],Cost0,Cost).
b.3 An Instance of the Transport Domain
Initial state  Goal state 
b.3.1 Solving the instance with Picat
main => Facts = $[road(c3,c1,40),road(c1,c3,40),road(c3,c2,18), road(c2,c3,18),road(c4,c1,36),road(c1,c4,36), road(c4,c3,37),road(c3,c4,37),road(c5,c2,24), road(c2,c5,24),road(c5,c3,26),road(c3,c5,26)], cl_facts(Facts,[$road(+,,)]), Trucks = [[c2,[],3],[c1,[],2]], Packages = [(c1,c2),(c1,c2),(c3,c1),(c2,c5)], best_plan({sort(Trucks),sort(Packages)},Plan,PlanCost), foreach ({I,Action} in zip(1..len(Plan),Plan)) printf("%3d. %w\n",I,Action) end, println(plan_cost=PlanCost).
b.3.2 An Optimal Plan for the Instance
1. load(c1) 2. load(c1) 3. load(c2) 4. move(c1,c3) 5. move(c2,c5) 6. unload(c5) 7. move(c3,c2) 8. unload(c2) 9. unload(c2) 10. move(c2,c3) 11. load(c3) 12. move(c3,c1) 13. unload(c1) plan_cost = 148