Reverse-search algorithm

<h2 id="principles">Principles</h2>
A reverse-search algorithm generates the combinatorial objects in a <a href="/facts/State_space/H1vjb3LL">state space</a>, an <a href="/facts/Implicit_graph/ERbgYekg">implicit graph</a> whose vertices are the objects to be listed and whose edges represent certain "local moves" connecting pairs of objects, typically by making small changes to their structure. It finds each objects using a <a href="/facts/Depth-first_search/luxJKVph">depth-first search</a> in a rooted <a href="/facts/Spanning_tree/qQPN9taR">spanning tree</a> of this state space, described by the following information:<a class="footnote-ref" id="fnref:3" href="#fn:3">3</a>

<ul><li>The root of the spanning tree, one of the objects</li>
<li>A subroutine for generating the parent of each object in the tree, with the property that if repeated enough times it will eventually reach the root</li>
<li>A subroutine for listing all of the neighbors in the state space (not all of which may be neighbors in the tree)</li></ul>
From this information it is possible to find the children of any given node in the tree, reversing the links given by the parent subroutine: they are simply the neighbors whose parent is the given node. It is these reversed links to child nodes that the algorithm searches.<a class="footnote-ref" id="fnref:4" href="#fn:4">4</a>
A classical depth-first search of this spanning tree would traverse the tree recursively, starting from the root, at each node listing all of the children and making a recursive call for each one. Unlike a depth-first search of a graph with cycles, it is not necessary to maintain the set of already-visited nodes to avoid repeated visits; such repetition is not possible in a tree. However, this recursive algorithm may still require a large amount of memory for its <a href="/facts/Call_stack/1sMt713v">call stack</a>, in cases when the tree is very deep. Instead, reverse search traverses the spanning tree in the same order while only storing two objects: the current object of the traversal, and the previously traversed object. Initially, the current object is set to the root of the tree, and there is no previous object. From this information, it is possible to determine the next step of the traversal by the following case analysis:<a class="footnote-ref" id="fnref:5" href="#fn:5">5</a>

<ul><li>If there is no previous object, or the previous object is the parent of the current object, then this is the first time the traversal has reached the current object, so it is output from the search. The next object is its first child or, if it has no children, its parent.</li>
<li>In all other cases, the previous object must be a child of the current object. The algorithm lists the children (that is, state-space neighbors of the current object that have the current object as their parent) one at a time until reaching this previous child, and then takes one more step in this list of children. If another child is found in this way, it is the next object. If there is no next child and the current object is not the root, the next object is the parent of the current object. In the remaining case, when there is no next child and the current object is the root, the reverse search terminates.</li></ul>
This algorithm involves listing the neighbors of an object once for each step in the search. However, if there are 
 
 
 
 N
 
 
 {\displaystyle N}
 
 objects to be listed, then the search performs 
 
 
 
 2
 N
 −
 1
 
 
 {\displaystyle 2N-1}
 
 steps, so the number of times it generates neighbors of objects is within a factor of two of the number of times the recursive depth-first search would do the same thing.<a class="footnote-ref" id="fnref:6" href="#fn:6">6</a>

<h2 id="applications">Applications</h2>
Examples of the problems to which reverse search has been applied include the following combinatorial generation problems:

Vertices of simple convex polytopes
If a 
 
 
 
 d
 
 
 {\displaystyle d}
 
-dimensional <a href="/facts/Convex_polytope/2pRHcRgC">convex polytope</a> is defined as an intersection of <a href="/facts/Half-space_(geometry)/AGknBgnw">half-spaces</a>, then its <a href="/facts/Vertex_(geometry)/ID3cea9z">vertices</a> can be described as the points of intersection of 
 
 
 
 d
 
 
 {\displaystyle d}
 
 or more <a href="/facts/Hyperplane/2s2ycYYd">hyperplanes</a> bounding the halfspaces; it is a simple polytope if no vertex is the intersection of more than 
 
 
 
 d
 
 
 {\displaystyle d}
 
 of these hyperplanes. The <a href="/facts/Vertex_enumeration_problem/VT6MJ27t">vertex enumeration problem</a> is the problem of listing all of these vertices. The edges of the polytope connect pairs of vertices that have 
 
 
 
 d
 −
 1
 
 
 {\displaystyle d-1}
 
 hyperplanes in common, so the vertices and edges form a state space in which each vertex has 
 
 
 
 d
 
 
 {\displaystyle d}
 
 neighbors. The <a href="/facts/Simplex_algorithm/9fndF6cq">simplex algorithm</a> from the theory of <a href="/facts/Linear_programming/GduXFQxT">linear programming</a> finds a vertex maximizing a given linear function of the coordinates, by walking from vertex to vertex, choosing at each step a vertex with a greater value of the function; there are several standard choices of "pivot rule" that specify more precisely which vertex to choose. Any such pivot rule can be interpreted as defining the parent function of a spanning tree of the polytope, whose root is the optimal vertex. Applying reverse search to this data generates all vertices of the polytope. A similar algorithm can also enumerate all bases of a linear program, without requiring that it defines a polytope that is simple.<a class="footnote-ref" id="fnref:7" href="#fn:7">7</a><a class="footnote-ref" id="fnref:8" href="#fn:8">8</a>
Cells of hyperplane arrangements
A <a href="/facts/Hyperplane_arrangement/78c23rYy">hyperplane arrangement</a> decomposes <a href="/facts/Euclidean_space/R2UbzmzM">Euclidean space</a> into cells, each described by a "sign vector" that describes whether its points belong to one of the hyperplanes (sign 0), are on one side of the hyperplane (sign +1), or are on the other side (sign −1). The cells form a connected state space under local moves that change a single sign by one unit, and it is possible to check that this operation produces a valid cell by solving a linear programming feasibility problem. A spanning tree can be constructed for any choice of root cell by defining a parent operator that makes the first possible change that would bring the sign vector closer to that of the root. Using reverse search for this state space and parent operator produces an algorithm for listing all cells in polynomial time per cell.<a class="footnote-ref" id="fnref:9" href="#fn:9">9</a><a class="footnote-ref" id="fnref:10" href="#fn:10">10</a>
Point-set triangulations
The <a href="/facts/Point-set_triangulation/e8u1mrxf">triangulations of a planar point set</a> are connected by "flip" moves that remove one diagonal from a triangulation and replace it by another. If the <a href="/facts/Delaunay_triangulation/6VRk7cS9">Delaunay triangulation</a> is chosen as the root, then every triangulation can be flipped to the Delaunay triangulation by steps in which the triangulation of some subset of four points is replaced by its Delaunay triangulation.<a class="footnote-ref" id="fnref:11" href="#fn:11">11</a><a class="footnote-ref" id="fnref:12" href="#fn:12">12</a> Choosing the first Delaunay flip as the parent of each triangulation, and applying local search, produces an algorithm for listing all triangulations in polynomial time per triangulation.<a class="footnote-ref" id="fnref:13" href="#fn:13">13</a>
Connected subgraphs
The connected subgraphs, and connected <a href="/facts/Induced_subgraph/G3EUC0MQ">induced subgraphs</a>, of a given connected graph form a state space whose local moves are the addition or removal of a single edge or vertex of the graph, respectively. A spanning tree of this state space can be obtained by adding the first edge or vertex (in some ordering of the edges or vertices) whose addition produces another connected subgraph; its root is the whole graph. Applying local search to this state space and parent operator produces an algorithm for listing all connected subgraphs in polynomial time per subgraph.<a class="footnote-ref" id="fnref:14" href="#fn:14">14</a>
Other applications include algorithms for generating the following structures:

<ul><li><a href="/facts/Polyomino/XUlHFIRY">Polyominos</a>,<a class="footnote-ref" id="fnref:15" href="#fn:15">15</a> <a href="/facts/Polyiamond/sATfqa4Y">polyiamond</a> <a href="/facts/Prototile/zh5RBHy5">prototiles</a>,<a class="footnote-ref" id="fnref:16" href="#fn:16">16</a> and <a href="/facts/Polyhex_(mathematics)/QsdwYxde">polyhex (mathematics)</a> <a href="/facts/Hydrocarbon/Dcg95Ju5">hydrocarbon</a> molecules.<a class="footnote-ref" id="fnref:17" href="#fn:17">17</a></li>
<li><a href="/facts/Topological_ordering/VY7gNtYA">Topological orderings</a> of <a href="/facts/Directed_acyclic_graph/k6zq1os9">directed acyclic graphs</a>, using a state space whose local moves reverse the ordering of two elements.<a class="footnote-ref" id="fnref:18" href="#fn:18">18</a></li>
<li><a href="/facts/Spanning_tree/qQPN9taR">Spanning trees</a> of graphs, non-crossing spanning trees of planar point sets, and more generally bases of <a href="/facts/Matroid/5bLcGlRF">matroids</a>, using a state space that swaps one edge for another.<a class="footnote-ref" id="fnref:19" href="#fn:19">19</a></li>
<li><a href="/facts/Euler_tour/4l2G67c7">Euler tours</a> in graphs.<a class="footnote-ref" id="fnref:20" href="#fn:20">20</a></li>
<li>The <a href="/facts/Maximal_independent_set/eB4NQsny">maximal independent sets</a> of <a href="/facts/Sparse_graph/xFC7b7s9">sparse graphs</a>.<a class="footnote-ref" id="fnref:21" href="#fn:21">21</a></li>
<li><a href="/facts/Maximal_planar_graph/plREagr9">Maximal planar graphs</a><a class="footnote-ref" id="fnref:22" href="#fn:22">22</a> and <a href="/facts/Polyhedral_graph/IKPCDH6r">polyhedral graphs</a>.<a class="footnote-ref" id="fnref:23" href="#fn:23">23</a></li>
<li>Non-crossing <a href="/facts/Laman_graph/vdg4s2zR">minimally rigid graphs</a> on a given point set.<a class="footnote-ref" id="fnref:24" href="#fn:24">24</a></li>
<li><a href="/facts/Polygonalization/4YArjFyH">Surrounding polygons</a>, polygons that have some of a given set of points as vertices and surround the rest, using a state space that adds or removes one vertex of the polygon.<a class="footnote-ref" id="fnref:25" href="#fn:25">25</a></li>
<li>Vertices or facets of the <a href="/facts/Minkowski_sum/dEQfoJ1A">Minkowski sum</a> of convex polytopes.<a class="footnote-ref" id="fnref:26" href="#fn:26">26</a><a class="footnote-ref" id="fnref:27" href="#fn:27">27</a></li>
<li>The corners (multidegrees) of <a href="/facts/Monomial_ideal/UhogyF4z">monomial ideals</a>.<a class="footnote-ref" id="fnref:28" href="#fn:28">28</a></li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1">Avis, David; Fukuda, Komei (1992), "A pivoting algorithm for convex hulls and vertex enumeration of arrangements and polyhedra", Discrete & Computational Geometry, 8 (3): 295–313, doi:10.1007/BF02293050, MR 1174359; preliminary version in Seventh Annual Symposium on Computational Geometry, 1991, doi:10.1145/109648.109659 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></li>
<li id="fn:2">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></li>
<li id="fn:3">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></li>
<li id="fn:4">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></li>
<li id="fn:5">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></li>
<li id="fn:6">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></li>
<li id="fn:7">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></li>
<li id="fn:8">Avis, David (2000), "A revised implementation of the reverse search vertex enumeration algorithm", in Kalai, Gil; Ziegler, Günter M. (eds.), Polytopes—combinatorics and computation: Including papers from the DMV-Seminar "Polytopes and Optimization" held in Oberwolfach, November 1997, DMV Seminar, vol. 29, Basel: Birkhäuser, pp. 177–198, MR 1785299 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></li>
<li id="fn:9">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></li>
<li id="fn:10">Sleumer, Nora H. (1999), "Output-sensitive cell enumeration in hyperplane arrangements", Nordic Journal of Computing, 6 (2): 137–147, MR 1709978 <a href="/wiki/MR_(identifier)" target="_blank">/wiki/MR_(identifier)</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></li>
<li id="fn:11">Lawson, C. L. (1972), Generation of a triangular grid with applications to contour plotting, Memo 299, Jet Propulsion Laboratory <a href="#fnref:11" class="footnote-back-ref">↩</a></li>
<li id="fn:12">Sibson, R. (1973), "Locally equiangular triangulations", The Computer Journal, 21 (3): 243–245, doi:10.1093/comjnl/21.3.243, MR 0507358 <a href="/wiki/Robin_Sibson" target="_blank">/wiki/Robin_Sibson</a> <a href="#fnref:12" class="footnote-back-ref">↩</a></li>
<li id="fn:13">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:13" class="footnote-back-ref">↩</a></li>
<li id="fn:14">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:14" class="footnote-back-ref">↩</a></li>
<li id="fn:15">Liang, Xiaodong; Wang, Rui; Meng, Ji xiang (2017), "Code for polyomino and computer search of isospectral polyominoes", Journal of Combinatorial Optimization, 33 (1): 254–264, doi:10.1007/s10878-015-9953-z, MR 3595411, S2CID 254655722 <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:15" class="footnote-back-ref">↩</a></li>
<li id="fn:16">Horiyama, Takashi; Yamane, Shogo (2011), "Generation of polyiamonds for p6 tiling by the reverse search", in Akiyama, Jin; Jiang, Bo; Kano, Mikio; Tan, Xuehou (eds.), Computational Geometry, Graphs and Applications - 9th International Conference, CGGA 2010, Dalian, China, November 3-6, 2010, Revised Selected Papers, Lecture Notes in Computer Science, vol. 7033, Springer, pp. 96–107, doi:10.1007/978-3-642-24983-9_10, ISBN 978-3-642-24982-2, MR 2927314 <a href="978-3-642-24982-2" target="_blank">978-3-642-24982-2</a> <a href="#fnref:16" class="footnote-back-ref">↩</a></li>
<li id="fn:17">Caporossi, Gilles; Hansen, Pierre (May 1998), "Enumeration of polyhex hydrocarbons to 
 
 
 
 h
 =
 21
 
 
 {\displaystyle h=21}
 
", Journal of Chemical Information and Computer Sciences, 38 (4): 610–619, doi:10.1021/ci970116n <a href="/wiki/Journal_of_Chemical_Information_and_Computer_Sciences" target="_blank">/wiki/Journal_of_Chemical_Information_and_Computer_Sciences</a> <a href="#fnref:17" class="footnote-back-ref">↩</a></li>
<li id="fn:18">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:18" class="footnote-back-ref">↩</a></li>
<li id="fn:19">Avis, David; Fukuda, Komei (1996), "Reverse search for enumeration", Discrete Applied Mathematics, 65 (1–3): 21–46, doi:10.1016/0166-218X(95)00026-N, MR 1380066 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:19" class="footnote-back-ref">↩</a></li>
<li id="fn:20">Kurita, Kazuhiro; Wasa, Kunihiro (2022), "Constant amortized time enumeration of Eulerian trails", Theoretical Computer Science, 923: 1–12, arXiv:2101.10473, doi:10.1016/j.tcs.2022.04.048, MR 4436557 <a href="/wiki/Theoretical_Computer_Science_(journal)" target="_blank">/wiki/Theoretical_Computer_Science_(journal)</a> <a href="#fnref:20" class="footnote-back-ref">↩</a></li>
<li id="fn:21">Eppstein, David (2009), "All maximal independent sets and dynamic dominance for sparse graphs", ACM Transactions on Algorithms, 5 (4): A38:1–A38:14, arXiv:cs/0407036, doi:10.1145/1597036.1597042, MR 2571901, S2CID 2769046 <a href="/wiki/David_Eppstein" target="_blank">/wiki/David_Eppstein</a> <a href="#fnref:21" class="footnote-back-ref">↩</a></li>
<li id="fn:22">Avis, David (1996), "Generating rooted triangulations without repetitions", Algorithmica, 16 (6): 618–632, doi:10.1007/s004539900067, MR 1412663 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:22" class="footnote-back-ref">↩</a></li>
<li id="fn:23">Deza, Antoine; Fukuda, Komei; Rosta, Vera (1994), "Wagner's theorem and combinatorial enumeration of 3-polytopes", Proceedings of a symposium held at the Research Institute for Mathematical Sciences, Kyoto University, Kyoto, May 17–19, 1993, RIMS Kôkyûroku Bessatsu, vol. 872, pp. 30–34, MR 1330480 <a href="/wiki/Komei_Fukuda" target="_blank">/wiki/Komei_Fukuda</a> <a href="#fnref:23" class="footnote-back-ref">↩</a></li>
<li id="fn:24">Avis, David; Katoh, Naoki; Ohsaki, Makoto; Streinu, Ileana; Tanigawa, Shin-ichi (June 2007), "Enumerating non-crossing minimally rigid frameworks" (PDF), Graphs and Combinatorics, 23 (S1): 117–134, doi:10.1007/s00373-007-0709-0, S2CID 10874512 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:24" class="footnote-back-ref">↩</a></li>
<li id="fn:25">Yamanaka, Katsuhisa; Avis, David; Horiyama, Takashi; Okamoto, Yoshio; Uehara, Ryuhei; Yamauchi, Tanami (2021), "Algorithmic enumeration of surrounding polygons" (PDF), Discrete Applied Mathematics, 303: 305–313, doi:10.1016/j.dam.2020.03.034, MR 4310502 <a href="/wiki/David_Avis" target="_blank">/wiki/David_Avis</a> <a href="#fnref:25" class="footnote-back-ref">↩</a></li>
<li id="fn:26">Fukuda, Komei (2004), "From the zonotope construction to the Minkowski addition of convex polytopes", Journal of Symbolic Computation, 38 (4): 1261–1272, doi:10.1016/j.jsc.2003.08.007, MR 2094220 <a href="/wiki/Komei_Fukuda" target="_blank">/wiki/Komei_Fukuda</a> <a href="#fnref:26" class="footnote-back-ref">↩</a></li>
<li id="fn:27">Weibel, Christophe (2010), "Implementation and parallelization of a reverse-search algorithm for Minkowski sums", in Blelloch, Guy E.; Halperin, Dan (eds.), Proceedings of the Twelfth Workshop on Algorithm Engineering and Experiments, ALENEX 2010, Austin, Texas, USA, January 16, 2010, Society for Industrial and Applied Mathematics, pp. 34–42, doi:10.1137/1.9781611972900.4, ISBN 978-0-89871-931-4 <a href="978-0-89871-931-4" target="_blank">978-0-89871-931-4</a> <a href="#fnref:27" class="footnote-back-ref">↩</a></li>
<li id="fn:28">Bayer, Dave; Taylor, Amelia (2009), "Reverse search for monomial ideals", Journal of Symbolic Computation, 44 (10): 1477–1486, doi:10.1016/j.jsc.2009.05.002, MR 2543431 <a href="/wiki/Dave_Bayer" target="_blank">/wiki/Dave_Bayer</a> <a href="#fnref:28" class="footnote-back-ref">↩</a></li>
</ol>

Reverse-search algorithm open-in-new

Reverse-search algorithm