Traversing Combinatorial 0/1-Polytopes via Optimization

SIAM Journal on Computing, Volume 53, Issue 5, Page 1257-1292, October 2024.
Abstract. In this paper, we present a new framework that exploits combinatorial optimization for efficiently generating a large variety of combinatorial objects based on graphs, matroids, posets, and polytopes. Our method is based on a simple and versatile algorithm for computing a Hamilton path on the skeleton of a 0/1-polytope [math], where [math]. The algorithm uses as a black box any algorithm that solves a variant of the classical linear optimization problem [math], and the resulting delay, i.e., the running time per visited vertex on the Hamilton path, is larger than the running time of the optimization algorithm only by a factor of [math]. When [math] encodes a particular class of combinatorial objects, then traversing the skeleton of the polytope [math] along a Hamilton path corresponds to listing the combinatorial objects by local change operations; i.e., we obtain Gray code listings. As concrete results of our general framework, we obtain efficient algorithms for generating all ([math]-optimal) bases and independent sets in a matroid; ([math]-optimal) spanning trees, forests, matchings, maximum matchings, and [math]-optimal matchings in a graph; vertex covers, minimum vertex covers, [math]-optimal vertex covers, stable sets, maximum stable sets, and [math]-optimal stable sets in a bipartite graph; as well as antichains, maximum antichains, [math]-optimal antichains, and [math]-optimal ideals of a poset. Specifically, the delay and space required by these algorithms are polynomial in the size of the matroid ground set, graph, or poset, respectively. Furthermore, all of these listings correspond to Hamilton paths on the corresponding combinatorial polytopes, namely the base polytope, matching polytope, vertex cover polytope, stable set polytope, chain polytope, and order polytope, respectively. As another corollary from our framework, we obtain an [math] delay algorithm for the vertex enumeration problem on 0/1-polytopes [math], where [math] and [math], and [math] is the time needed to solve the linear program [math]. This improves upon the 25-year-old [math] delay algorithm due to Bussieck and Lübbecke.