Shortest Paths in Graphs of Convex Sets I: The Perspective Mixed-Integer Convex Program (5.5) Has the Same Optimal Value as the Shortest-Path Problem in a Graph of Convex SetsResearch Paper
Motivation
Many planning problems combine a discrete choice with a continuous one. A robot moving through a cluttered room chooses which free regions to traverse and where exactly to cross between them; a vehicle chooses a sequence of operating modes and a trajectory within each mode. Marcucci, Umenberger, Parrilo and Tedrake (arXiv:2101.11565, SIAM J. Optim. 2024, doi:10.1137/22M1523790) abstract this as the shortest-path problem in a graph of convex sets (SPP in GCS): a shortest-path problem in which each vertex is not a point but a convex set, and the length of an edge depends on where in its two sets the path enters and leaves. The paper shows that the problem is NP-hard even in simple cases, and then gives a mixed-integer convex formulation that is exact and whose convex relaxation is strong in practice. That formulation underlies the GCS motion planners developed later by the same group.
Setting
A graph of convex sets consists of a finite directed graph , a source vertex , a target vertex , a nonempty compact convex set for each vertex , and for each edge an edge length , which is proper, closed and convex. No symmetry or triangle inequality is assumed. An - path is a sequence of distinct vertices with , and ; is the set of its edges. The SPP in GCS (2.1) is
The perspective of a set is the closed cone (Definition 4.1). The perspective of a function is , the function whose epigraph is the perspective of the epigraph of (Definition 4.4). For it equals .
With and the edges entering and leaving , the MICP (5.5) has a flow and two vectors per edge:
Its optimal value is . The paper derives (5.5) from a biconvex program (5.2), the shortest-path LP (5.1) together with vertex positions and the bilinear constraints , . Lemma 5.4 turns every valid linear inequality on the flows at a vertex into a valid perspective constraint.
Formalization targets
Goal: Theorem 5.7 (p. 10)
Under the standing assumptions of §2 and the assumption of p. 7 that no edge enters and no edge leaves ,
Both sides are infima in the extended reals and equal when there is no - path.
Milestones
- Remark 4.2 (p. 5): for bounded closed , .
- Remark 4.5 (p. 6): for , for , for proper closed convex .
- Lemma 5.4 (p. 9): a valid inequality for (5.2) lifts to the valid constraint .
- Remark 5.5 (p. 9): the equality version of Lemma 5.4.
- Path and cycles (proof of Theorem 5.7): binary flows satisfying (5.5b)–(5.5d) form the vertex-disjoint union of an - path and cycles.
- Simplification of (5.5e) (proof of Theorem 5.7): forces at cost ; forces , at cost .
- Flow along the path (proof of Theorem 5.7): for consecutive path edges , .
A companion item states the recovery part of Theorem 5.7 (second and third sentences): an optimal MICP solution yields an optimal path and optimal positions , .
Significance
Theorem 5.7 is what makes (5.5) usable: it says that the mixed-integer program, whose only nonconvexity is , solves the SPP in GCS exactly. The SPP in GCS is NP-hard by Theorem 3.1, so exactness is the most one can ask of a polynomial-size formulation, and dropping integrality gives a convex program whose value is a lower bound for branch and bound. The perspective construction of Lemma 5.4 is a general recipe, closely related to the reformulation-linearization technique, for convexifying products of flows and positions.
On the formal side, the paper's proof is half a page and leaves several points to the reader: why the integer flows decompose into a path and cycles, why the closure in the perspective cone disappears for compact sets, and why the perspective function vanishes at the origin. A machine-checked proof makes these steps explicit. The result has not been formalized before, and the mission also produces reusable statements about perspectives of sets and functions in extended-real arithmetic and about the decomposition of unit - flows with degree bounds.
Difficulty
The inequality is the routine half: a path and its positions give a feasible point of (5.5) of the same cost. The other half requires reading a path out of an arbitrary feasible point. Integer flows satisfying (5.5b)–(5.5d) need not be a path: they may also contain disjoint cycles, and the degree constraint (5.5c) is what prevents a cycle from sharing a vertex with the path (Example 5.10 shows the value drops without it). The cost of a cycle must be shown to be nonnegative, and the perspective constraints must be shown to collapse to on zero-flow edges. That step uses compactness of . For unbounded sets the zero slice of is the recession cone, and the argument fails. The cost at zero flow is , defined by a closure, and its value uses properness and convexity of , not just the formula .
Formalization scope
- is
Fin n → ℝ. Edge lengths and all optimal values areEReal; optimal values are infima (⨅/sInf), equal to on an empty feasible set. Attainment is not assumed. - The graph is a structure with a finite edge set
E : Finset (V × V)over a finite vertex type, a source, a target, setsX vand lengthsℓ e. The standing assumptions of §2 (; nonempty, compact, convex; proper closed convex) are the predicateIsGCS. The assumption of p. 7 is the separate predicateNoInSNoOutT, a hypothesis of the goal, of the path-and-cycles milestone and of the recovery item. - Paths are duplicate-free vertex lists from to along edges; is the list of consecutive pairs.
- The perspective of a function is Definition 4.4's infimum over the closed perspective cone of the epigraph. The formula is not used as a definition: in Lean it would evaluate to at for every argument and erase the content of (5.5e). The integrality constraint is part of the feasible set; without it the goal is false (§9.4 reports relaxation gaps).
- Disclosed additions: Lemma 5.4 and Remark 5.5 assume no self-loop at , which the page's partition of into incoming and outgoing edges presupposes. The recovery item takes the - path contained in , since the page's claim that this set is itself a path fails when an optimal solution carries a zero-length cycle.
- Needed infrastructure: closures of cones over compact sets, extended-real infima, and the combinatorics of unit flows on finite digraphs. Proofs of any milestone, and reusable lemmas on perspectives, are welcome.
Selected references
- T. Marcucci, J. Umenberger, P. A. Parrilo, R. Tedrake, Shortest Paths in Graphs of Convex Sets, SIAM Journal on Optimization 34(1), 2024; arXiv:2101.11565v5. https://arxiv.org/abs/2101.11565
- R. T. Rockafellar, Convex Analysis, Princeton University Press, 1970 (perspectives and recession cones, §8). https://doi.org/10.1515/9781400873173
- R. K. Ahuja, T. L. Magnanti, J. B. Orlin, Network Flows: Theory, Algorithms, and Applications, Prentice Hall, 1993 (flow decomposition, shortest-path LP).
- T. Marcucci, M. Petersen, D. von Wrangel, R. Tedrake, Motion Planning around Obstacles with Convex Optimization, Science Robotics 8(84), 2023. https://doi.org/10.1126/scirobotics.adf7843