Computational Complexity of the $α$-Ham-Sandwich Problem

Man-Kwun Chiu, Aruni Choudhary, Wolfgang Mulzer

Introduction

There are many alternative and more general variants of both the continuous and the discrete Ham-Sandwich Theorem. For example, Bárány and Matoušek [BM01] derived a version where measures in the plane can be divided into any (possibly different) ratios by fans instead of hyperplanes (lines). A discrete variant of this result was given by Bereg [Ber05]. Schnider [Sch19a] studied a generalization in higher dimensions. Recently Barba, Pilz, and Schnider [BPS19] showed that four measures in the plane can be bisected with two lines. Zivaljević and Vrećica [ZV90] proved a result that interpolates between the Ham-Sandwich Theorem and the Centerpoint Theorem [Rad46], of which there is also a no-dimensional version [CM20]. Schnider [Sch19b] presented a generalization based on this result among others.

Steiger and Zhao called their result the Generalized Ham-Sandwich Theorem, yet it is not a strict generalization of the classic Ham-Sandwich Theorem. Their result requires that the point sets obey well-separation and weak general position, while the classic theorem always holds without these assumptions. Therefore, we call this result the α\alpha-Ham-Sandwich theorem, for a clearer distinction. Set n=∑i∈[d]∣Si∣n=\sum_{i\in[d]}|S_{i}|. Steiger and Zhao gave an algorithm that computes the dividing hyperplane in O(n(log⁡n)d−3)O\left(n(\log n)^{d-3}\right) time, which is exponential in dd. Later, Bereg [Ber12] improved this algorithm to achieve a running time of n2O(d)n2^{O(d)}, which is linear in nn but still exponential in dd. We denote the associated computational search problem of finding the dividing hyperplane as Alpha-HS.

No polynomial algorithms are known for Ham-Sandwich and for Alpha-HS if the dimension is not fixed, and the notion of approximation is also not well-explored. Despite their superficial similarity, it is not immediately apparent whether the two problems are comparable in terms of their complexity. Due to the additional requirements on an input for Alpha-HS, an instance of Ham-Sandwich may not be reducible to Alpha-HS in general.

Since a dividing hyperplane for Alpha-HS is guaranteed to exist if the sets satisfy the conditions of well-separation and (weak) general position, Alpha-HS is a total search problem. In general, such problems are modelled by the complexity class TFNP\mathsf{TFNP} (Total Function Nondeterministic Polynomial) of NP\mathsf{NP}-search problems that always admit a solution. Two popular subclasses of TFNP\mathsf{TFNP}, originally defined by Papadimitriou [Pap94], are PPA\mathsf{PPA} (Polynomial Parity Argument) its sub-class PPAD\mathsf{PPAD}. These classes contain total search problems where the existence of a solution is based on a parity argument in an undirected or in a directed graph, respectively. Another sub-class of TFNP\mathsf{TFNP} is PLS\mathsf{PLS} (polynomial local search). It models total search problems where the solutions can be obtained as minima in a local search process, while the number of steps in the local search may be exponential in the input size. The class PLS\mathsf{PLS} was introduced by Johnson, Papadimitriou, and Yannakakis [JPY88]. A noteworthy sub-class of PPAD∩PLS\mathsf{PPAD}\cap\mathsf{PLS} is CLS\mathsf{CLS} (continuous local search) [DP11]. It models similar local search problems over a continuous domain using a continuous potential function.

Up to very recently, these complexity classes have mostly been studied in the context of algorithmic game theory. However, there have been increasing efforts towards mapping the complexity landscape of existence theorems in high-dimensional discrete geometry. Computing an approximate solution for the search problem associated with the Borsuk-Ulam Theorem is in PPA\mathsf{PPA}. In fact, this problem is complete for this class. The discrete analogue of the Borsuk-Ulam Theorem, Tucker’s Lemma [Tuc46], is also PPA\mathsf{PPA}-complete [ABB20]. Therefore, since the traditional proof of the Ham-Sandwich Theorem goes through the Borsuk-Ulam Theorem, it follows that Ham-Sandwich lies in PPA\mathsf{PPA}. In fact, Filos-Ratsikas and Goldberg [FRG19] recently showed that Ham-Sandwich is complete for PPA\mathsf{PPA}. The (presumably smaller) class PPAD\mathsf{PPAD} is associated with fixed-point type problems: computing an approximate Brouwer fixed point is a prototypical complete problem for PPAD\mathsf{PPAD}. The discrete analogue of Brouwer’s Fixed Point Theorem, Sperner’s Lemma, is also complete for PPAD\mathsf{PPAD}. In a celebrated result, the relevance of PPAD\mathsf{PPAD} for algorithmic game theory was made clear when it turned out that computing a Nash-equilibrium in a two player game is PPAD\mathsf{PPAD}-complete [CDT09]. In discrete geometry, finding a solution to the Colorful Carathéodory problem [Bár82] was shown to lie in the intersection PPAD∩PLS\mathsf{PPAD}\cap\mathsf{PLS} [MMSS17, MS18]. This further implies that finding a Tverberg partition (and computing a centerpoint) also lies in the intersection [Tve66, Sar92, LGMM19]. The problem of computing the (unique) fixed point of a contraction map is known to lie in CLS\mathsf{CLS} [DP11].

Recently, at ICALP 2019, Fearley, Gordon, Mehta, and Savani defined a sub-class of CLS\mathsf{CLS} that represents a family of total search problems with unique solutions [FGMS19]. They named the class Unique End of Potential Line (UEOPL\mathsf{UEOPL}) and defined it through the canonical complete problem UniqueEOPL. This problem is modelled as a directed graph. There are polynomially-sized Boolean circuits that compute the successor and predecessor of each node, and a potential value that always increases on a directed path. There is supposed to be only a single vertex with no predecessor (start of line). Under these conditions, there is a unique path in the graph that ends on a vertex (called end of line) with the highest potential along the path. This vertex is the solution to UniqueEOPL. Since the uniqueness of the solution is guaranteed only under certain assumptions, such a formulation is called a promise problem. Since there seems to be no efficient way to verify the assumptions, the authors allow two possible outcomes of the search algorithm: either report a correct solution, or provide any solution that was found to be in violation of the assumptions. This formulation turns UniqueEOPL into a non-promise problem and places it in TFNP\mathsf{TFNP}, since a correct solution is bound to exist when there are no violations, and otherwise a violation can be reported as a solution. Fearley et al. [FGMS19] also introduced the concept of a promise-preserving reduction between two problems AA and BB, such that if an instance of AA has no violations, then the reduced instance of BB is also free of violations. This notion is particularly meaningful for non-promise problems.

Contributions.

We provide the first non-trivial containment in a complexity class for the α\alpha-Ham-Sandwich problem by locating it in UEOPL\mathsf{UEOPL}. More precisely, we formulate Alpha-HS as a non-promise problem in which we allow for both valid solutions representing the correct dividing hyperplane, as well as violations accounting for the lack of well-separation and/or (weak) general position of the input point sets. A precise formulation of the problem is given in Definition 4 in Section 2. We then show a promise-preserving reduction from Alpha-HS to UniqueEOPL. This implies that Alpha-HS lies in UEOPL\mathsf{UEOPL}, and hence in CLS⊆PPAD∩PLS\mathsf{CLS}\subseteq\mathsf{PPAD}\cap\mathsf{PLS}. See Figure 2 for a pictorial description.

It is not surprising to discover that Alpha-HS lies in PPAD\mathsf{PPAD}, since the proof of the continuous version in [BHJ08] was based on Brouwer’s Fixed Point Theorem. The observation that it also lies in PLS\mathsf{PLS} is new and noteworthy, putting Alpha-HS into the reach of local search algorithms. In contrast, given our current understanding of total search problems, it is unlikely that the problem Ham-Sandwich would be in PLS\mathsf{PLS}.

Since Alpha-HS lies in PPAD⊆PPA\mathsf{PPAD}\subseteq\mathsf{PPA}, it is computationally easier than Ham-Sandwich, which is PPA\mathsf{PPA}-complete. This implies the existence of a polynomial-time reduction from Alpha-HS to Ham-Sandwich. A reduction in the other direction is unlikely. It thus turns out that well-separation brings down the complexity of the problem by a significant amount.

Often, problems in TFNP\mathsf{TFNP} come in the guise of a polynomial-size Boolean circuit with some property. In contrast, Alpha-HS is a purely geometric problem that has no circuit in its problem definition. This is the second problem in UEOPL\mathsf{UEOPL} apart from the PP-Matrix Linear complementarity problem and one of the few in CLS\mathsf{CLS} that does not have a description in terms of circuits.

Our local-search formulation is based on the intuition of rotating a hyperplane until we reach the desired solution. We essentially start with a hyperplane that is tangent to the convex hull of each input set, and we deterministically rotate the hyperplane until it hits a new point. This rotation can be continued whenever the hyperplane hits a new point, until we reach the correct dividing hyperplane. In other words, we can follow a local-search argument to find the solution. We show that this sequence of rotations can be modelled as a canonical path in a grid graph, and we give a potential function that guides the rotation and always increases along this path. Every violation of well-separation and (weak) general position can destroy this path. Furthermore, no efficient methods to verify these two assumptions are known. This poses a major challenge in handling the violations. One of our main technical contributions is to handle the violation solutions concisely.

An alternative approach would have been to look at the dual space of points where we get an arrangement of hyperplanes. The dividing hyperplane could then be found by looking at the correct level sets of the arrangement. However, this approach has the problem that the orientations of the hyperplanes in the original space and the dual space are not consistent. This complicates the arguments on the level sets, so we found it more convenient to use our notion of rotating hyperplanes. We show that we can maintain a consistent orientation throughout the rotation, and an inconsistent rotation is detected as a violation of the promise.

Outline of the paper.

We discuss the background about the α\alpha-Ham-sandwich Theorem and UniqueEOPL in Section 2. In Section 3, we describe our instance of Alpha-HS and give an overview of the reduction and violation-handling. The technical details of the reduction are presented in Section 4 and Section 5. We conclude in Section 6.

Preliminaries

For conciseness, we describe the discrete version of α\alpha-Ham-Sandwich Theorem [SZ10] here. The continuous version [BHJ08] follows a similar formulation.

We say that PP has very weak general position [SZ10], if for every choice of points x1∈P1,…,xd∈Pdx_{1}\in P_{1},\dots,x_{d}\in P_{d}, the affine hull of the set {x1,…,xd}\{x_{1},\dots,x_{d}\} is a (d−1)(d-1)-flat and does not contain any other point of PP. This definition is sufficient for the result of Steiger and Zhao, where they simply call it as weak general position. Of course, this definition of weak general position has no restriction on sets {x1,…,xd}\{x_{1},\dots,x_{d}\} that contain multiple points from the same color. To simplify our proofs we need a slightly stronger form of general position. We say that PP has weak general position if the above restriction also applies to sets having exactly d−1d-1 colors. That means, each color may contribute at most one point to the set, except perhaps one color which is allowed to contribute two points. A certificate for checking violations of weak general position is a set of d+1d+1 points whose affine hull has dimension at most d−1d-1, with at least d−1d-1 colors in the set. Testing whether a planar point set is in general position can be shown to be NP\mathsf{NP}-Hard, using the result in [FKNN17]. It is easy to see that when d=2d=2, weak general position is equivalent to general position.

Well-separation.

A certificate for checking violations of well-separation is a colorful set {x1,…,xd}\{x_{1},\dots,x_{d}\} whose affine hull has dimension at most d−2d-2. Another certificate is a partition I,J⊂[d]I,J\subset[d] such that the convex hulls of the indexed sets are not separable. Due to Lemma 1, both certificates are equivalent and either can be converted to the other in polynomial time. To the best of our knowledge, the complexity of testing well-separation is unknown.

Given any set of positive integers {α1,…,αd}\{\alpha_{1},\dots,\alpha_{d}\} satisfying 1≤αi≤ni1\leq\alpha_{i}\leq n_{i}, i∈[d]i\in[d], an (α1,…,αd)(\alpha_{1},\dots,\alpha_{d})-cut is an oriented hyperplane HH that contains one point from each color and satisfies ∣H+∩Pi∣=αi|H^{+}\cap P_{i}|=\alpha_{i} for i∈[d]i\in[d], where H+H^{+} is the closed positive half-space defined by HH.

If an α\alpha-cut exists, then it is unique.

If PP has weak general position, then a cut exists for each choice of α\alpha, αi∈[ni]\alpha_{i}\in[n_{i}].

That means, every colorful dd-tuple of PP corresponds to exactly one α\alpha-vector. Steiger and Zhao [SZ10] also presented an algorithm to compute the cut in O(n(log⁡n)d−3)O(n(\log n)^{d-3}) time, where n=∑i=1dnin=\sum_{i=1}^{d}n_{i}. The algorithm proceeds inductively in dimension and employs a prune-and-search technique. Bereg [Ber12] improved the pruning step to improve the runtime to n2O(d)n2^{O(d)}.

2 Unique End of Potential Line

We briefly explain the Unique end of potential line problem that was introduced in [FGMS19]. More details about the problem and the associated class can be found in the above reference.

Let n,mn,m be positive integers. The input consists of

The complexity class UEOPL\mathsf{UEOPL} represents the class of problems that can be reduced in polynomial time to UniqueEOPL. This has been shown to lie in CLS\mathsf{CLS} in [FGMS19] and contains three classical problems: finding the fixed point of a contraction map, solving the P-Matrix Linear complementarity problem, and finding the unique sink of a directed graph (with arbitrary edge orientations) on the 1-skeleton of a hypercube.

A notion of promise-preserving reductions is also defined in [FGMS19]. Let XX and YY be two problems both having a formulation that allows for valid and violation solutions. A reduction from XX to YY is said to be promise-preserving, if whenever it is promised that XX has no violations, then the reduced instance of YY also has no violations. Thus a promise-preserving reduction to UniqueEOPL would mean that whenever the original problem is free of violations, then the reduced instance always has a single line that ends at a valid solution.

3 Formulating the search problem

We formalize the search problem for α\alpha-Ham-Sandwich in a non-promise setting:

A subset of PP of size d+1d+1 and at least d−1d-1 colors that lies on a hyperplane.

Here a solution of type (G1) corresponds to a solution representing a valid cut, while solutions of type (GV1) and (GV2) refer to violations of weak general position and well-separation, respectively. From Theorem 2 we see that a valid solution is guaranteed if no violations are presented, which shows that Alpha-HS is a total search problem.

Alpha-HS is in UEOPL

In this section we describe our instance of Alpha-HS in more detail and briefly outline a reduction to UniqueEOPL.

1 An overview of the reduction

We give a short overview of the ideas used in the reduction from Alpha-HS to UniqueEOPL. The details are technical and we defer them to Section 5. We encourage the interested reader to go through the details of our reduction.

Our intuition is based on rotating a colorful hyperplane HH to another colorful hyperplane H′H^{\prime} through a sequence of local changes of the points on the hyperplanes such that the α\alpha-vector of H′H^{\prime} increases in some coordinate by one from that of HH. We next define the rotation operation in a little more detail. An anchor is a colorful (d−1)(d-1)-tuple of PP which spans a (d−2)(d-2)-flat. The following procedure takes as input an anchor RR and some point p∈P∖Rp\in P\setminus R and determines the next hyperplane obtained by a rotation. The output is (R′,p′)(R^{\prime},p^{\prime}), where R′R^{\prime} is an anchor and p′∈P∖R′p^{\prime}\in P\setminus R^{\prime} is some point.

Procedure (R′,p′)=NextRotate(R,p)(R^{\prime},p^{\prime})=NextRotate(R,p)

Let HH denote the hyperplane defined by R∪{p}R\cup\{p\} and t1t_{1} be the missing color in RR.

If the orientation of HH is not well-defined, report a violation of weak general position and well-separation.

Let Pt1+P^{+}_{t_{1}} be the subset of Pt1P_{t_{1}} that lies in the closed halfspace H+H^{+} and Pt1−P^{-}_{t_{1}} be the subset of Pt1P_{t_{1}} that lies in the open halfspace H−H^{-}. Let x∈Pt1+x\in P^{+}_{t_{1}} be the highest ranked point according to the order ≺t1\prec_{t_{1}} and y∈Pt1−y\in P^{-}_{t_{1}} be the highest ranked point according to ≺t1\prec_{t_{1}}.

If pp has color t1t_{1} and ∣Pt1+∣=nt1|P^{+}_{t_{1}}|=n_{t_{1}}, report out of range.

We rotate HH around the anchor RR in a direction such that the hyperplane is moving away from xx along the segment xyxy until it hits some point q∈Pq\in P.

If the hyperplane hits multiple points at the same time, report a violation of weak general position.

If p′p^{\prime} is not color t1t_{1}, set R′:=R∪{q}∖{r}R^{\prime}:=R\cup\{q\}\setminus\{r\} and p′=rp^{\prime}=r, where rr is a point in RR with the same color as p′p^{\prime}. Otherwise, set R′=RR^{\prime}=R and p′=qp^{\prime}=q.

Figure 3 shows an application of this procedure, rotating H0H_{0} to H4H_{4} through H1,H2,H3H_{1},H_{2},H_{3}.

This rotation function can be interpreted as a function that assigns each hyperplane to the next hyperplane. The set of colorful hyperplanes can be interpreted as vertices in a graph with the rotation function determining the connectivity of the graph.

Each colorful hyperplane HH is incident to a colorful set of dd points. This set of points defines dd possible anchors, and each anchor can be used to rotate HH in a different fashion. To define a unique sequence of rotations, we pick a specific order as follows: first, we assume that the colorful hyperplane HH whose α\alpha-vector is (1,…,1)(1,\ldots,1) is given (we show later how this assumption can be removed). We start at HH and pick the anchor that excludes the first color, then apply a sequence of rotations until we hit another colorful hyperplane with α\alpha-vector (2,1,…,1)(2,1,\ldots,1). Similarly, we move to a colorful hyperplane with α\alpha-vector (3,1,…,1)(3,1,\ldots,1) and so on until we reach (α1,1,…,1)(\alpha_{1},1,\ldots,1). Then, we repeat this for the other colors in order to reach (α1,α2,1,…,1)(\alpha_{1},\alpha_{2},1,\ldots,1) and so on until we reach the target α\alpha-vector. This pattern of α\alpha-vectors helps in defining a potential function that strictly increases along the path. We can encode this sequence of rotations as a unique path in the UniqueEOPL instance, and we call it canonical path.

A natural way to define the UniqueEOPL graph would be to consider hyperplanes as the vertices in the graph. However, this leads to complications. Figure 3 shows a rotation from H0H_{0} to H4H_{4}, with α\alpha-vectors (2,3)(2,3) and (3,3)(3,3) respectively. During the rotation, we encounter a hyperplane H2H_{2} for which its α\alpha-vector is (2,4)(2,4), which differs from our desired sequence of (2,3),…,(2,3),(3,3)(2,3),\dots,(2,3),(3,3). This makes it difficult to define a potential function in the graph that strictly increases along the path vH0,…,vH4v_{H_{0}},\dots,v_{H_{4}} where vHiv_{H_{i}} is the vertex representing hyperplane HiH_{i}. One way to alleviate this problem is to not use HiH_{i} as a vertex directly, but the double-wedge that is traced out by the rotation from HiH_{i} to Hi+1H_{i+1}. If the α\alpha-vector is now measured using the hyperplane that bisects the double-wedge, then we get the desired sequence of (2,3),…,(2,3),(3,3)(2,3),\dots,(2,3),(3,3). See Figure 3 for an example.

With additional overhead, the rotation function can be extended to double-wedges. This in turn also leads to a neighborhood graph where the vertices are the double-wedges and the rotations can be used to define the edges. The graph is connected and has a grid-like structure that may be of independent interest. To simplify the exposition, we postpone the description of double-wedges and the associated graph to Section 4.

Distance parameter and potential function.

The α\alpha-vector is not sufficient to define the potential function, since the sequence of rotations between two colorful hyperplanes may have the same α\alpha-vector. For instance, the bisectors of the rotations in H0,…,H3H_{0},\dots,H_{3} in Figure 3 all have the same α\alpha-vector. Hence, we need an additional measurement in order to determine the direction of rotation that increases the α\alpha-vector.

Similar to how we define the orientation for a non-colorful hyperplane, let HH denote a hyperplane that passes through points of (d−1)(d-1) colors. Let PjP_{j} denote the missing color in HH. Let x,y∈Pjx,y\in P_{j} be the highest ranked points under ≺j\prec_{j} in H+H^{+} and H−H^{-} respectively. Let zz denote the intersection of xyxy and HH. We define a distance parameter called dist-value of HH to be the distance ∥x−z∥\|x-z\|. In Figure 3, we can see that rotating from H0H_{0} to H4H_{4} sweeps the segment xyxy in one direction, with the dist-value of the hyperplanes increasing strictly. This is sufficient to break ties and hence determine the correct direction of rotation. The precise statement is given in Lemma 6. We can extend this definition to the domain of double-wedges. We define a potential value for each vertex on the canonical path in UniqueEOPL using the sum of weighed components of α\alpha-vector and dist-value for the tie-breaker.

Correctness.

We show that if there are no violations, we can always apply Procedure NextRotateNextRotate to increment the α\alpha-vector until we find the desired solution, which implies that the canonical path exists. If the input satisfies weak general position, we can see that the rotating hyperplane always hits a unique point in Step 55, which may be swapped to form a new anchor in Step 77.

The well-separation condition guarantees that the potential function always increases along the rotation. Let H1,H2H_{1},H_{2} denote a pair of hyperplanes that are the input and output of Procedure NextRotateNextRotate respectively. Let HH denote any intermediate hyperplane during the rotation from H1H_{1} to H2H_{2} through the common anchor. Let PjP_{j} be the color missing from the anchor and xx be the highest ranked point under ≺j\prec_{j} in H1+H_{1}^{+}. We say that the orientation of H2H_{2} (resp. HH) is consistent with that of H1H_{1} if x∈H2+x\in H_{2}^{+} (resp. x∈H+x\in H^{+}). Lemma 5 shows that the orientations are always consistent when H1H_{1} and H2H_{2} are non-colorful hyperplanes even without the assumption of well-separation.

Assume that weak general position holds. Let H1,H2H_{1},H_{2} be the input and output of Procedure NextRotateNextRotate respectively. Let HH denote any intermediate hyperplane within the rotation. The orientations of H1H_{1} (resp. H2H_{2}) and HH are consistent when H1H_{1} (resp. H2H_{2}) is a non-colorful hyperplane.

Since H1H_{1} is a non-colorful hyperplane, let PjP_{j} denote the color missing from H1H_{1}. H1H_{1} and HH give the same partition of PjP_{j} into two sets because the continuous rotation from H1H_{1} to HH does not hit any point in PjP_{j}. Let xx and yy be the highest ranked points under ≺j\prec_{j} in each set. Since we have weak general position, the segment xyxy cannot pass through the anchor of the rotation so that the orientations of H1H_{1} and HH are well-defined by the (d−1)(d-1) colored points in the anchor and the intersections of the hyperplanes with the segment xyxy. Thus, the determinant defining the normal of the rotating hyperplane from H1H_{1} to HH for the orientation is always non-zero. Since the intersection of the rotating hyperplane from H1H_{1} to HH and the segment xyxy moves continuously along xyxy, by a continuity argument, the normal of the hyperplane does not flip during the rotation. Without loss of generality, assume that x∈H1+x\in H_{1}^{+}. This implies that xx is always in the positive half-space of HH and hence HH has a consistent orientation as H1H_{1}. The same proof holds for H2H_{2}. ∎

Next, we show that the dist-value is strictly increasing for all the intermediate hyperplanes in the sequence of rotations from one colorful hyperplane to another colorful hyperplane.

Assume that weak general position holds. Let H0H_{0} be a colorful hyperplane and HkH_{k} be the first colorful hyperplane obtained by a sequence of rotations by Procedure NextRotateNextRotate. We denote H1,…,Hk−1H_{1},\ldots,H_{k-1} be the non-colorful hyperplanes obtained from the above sequence of rotations. The dist-values of H1,…,Hk−1H_{1},\ldots,H_{k-1} is strictly increasing.

Let PjP_{j} denote the color missing from H1H_{1}. Then, H2,…,Hk−1H_{2},\ldots,H_{k-1} all miss the color PjP_{j}, otherwise HkH_{k} is not the first colorful hyperplane obtained by the rotations. Therefore, each HiH_{i} gives the same partition of PjP_{j} into two sets for i=1,…,k−1i=1,\ldots,k-1 because the continuous rotations from H1H_{1} to Hk−1H_{k-1} does not hit any point in PjP_{j}. Let xx and yy be the highest ranked points under ≺j\prec_{j} in each set. Without loss of generality, assume that x∈H1+x\in H_{1}^{+}. Since H1,…,Hk−1H_{1},\ldots,H_{k-1} are non-colorful hyperplanes, by Lemma 5, the consistent of the orientation can carry from H1H_{1} to H2H_{2} and so on. Then we have x∈H1+,…,x∈Hk−1+x\in H_{1}^{+},\dots,x\in H_{k-1}^{+} and y∈H1−,…,y∈Hk−1−y\in H_{1}^{-},\dots,y\in H_{k-1}^{-}. Let z1=xy∩H1,…,zk−1=xy∩Hk−1z_{1}=xy\cap H_{1},\dots,z_{k-1}=xy\cap H_{k-1}. According to Step 55 of Procedure NextRotateNextRotate, each rotation is performed by moving away from xx along the segment xyxy. Hence we have ∥x−z1∥<∥x−z2∥<⋯<∥x−zk−1∥\|x-z_{1}\|<\|x-z_{2}\|<\dots<\|x-z_{k-1}\|. ∎

The last step for proving that the potential function always increases along the canonical path is to show that the α\alpha-vector increases in some coordinate from one colorful hyperplane to another colorful hyperplane through Procedure NextRotateNextRotate. This requires the assumption of well-separation. Lemma 7 shows that if the orientations of H1,H2H_{1},H_{2} and HH are inconsistent, then well-separation is violated. By the contrapositive, if well-separation is satisfied, then all hyperplanes in the rotation always give consistent orientations. Then, it implies that rotating from a colorful hyperplane H0H_{0} to another colorful hyperplane HkH_{k} through a sequence of non-colorful hyperplanes that miss color PjP_{j}, we have H0+∩Pj⊂Hk+∩PjH_{0}^{+}\cap P_{j}\subset H_{k}^{+}\cap P_{j} and HkH_{k} contains one additional point in PjP_{j} that is hit by the last rotation. Therefore, αj\alpha_{j} is increased by 11 and other αi\alpha_{i}s keep the same value because of the way we swap the point of repeated color with the one in the anchor and the direction of rotation.

Since the orientations of H1H_{1} and HH are inconsistent, H1H_{1} must be a colorful hyperplane by Lemma 5. Therefore, the point in H1H_{1} that is not in the anchor is in PjP_{j}, denoted by pp.

In order to guarantee that there is no other path in UniqueEOPL apart from the canonical path, we introduce self-loops for vertices that are not on the canonical path. The detailed proof is given in Lemma 17 that if there are no violations, then the reduced instance of UniqueEOPL only gives a (U1) solution, which readily translates to a (G1) solution, so our reduction is promise-preserving, and this can be done in polynomial time.

Since we do not know the hyperplane with α\alpha-vector (1,…,1)(1,\dots,1) in advance, we split the problem into two sub-problems: in the first we start with any colorful hyperplane. We reverse the direction of the canonical path determined by the potential and construct an Alpha-HS instance for which the vertex with α\alpha-vector (1,…,1)(1,\dots,1) is the solution. In the second, we use this vertex as the input to the main Alpha-HS instance. If the input is free of violations, then both sub-problems give valid solutions and together they answer the original question.

Handling violations.

The reduction maps violations of Alpha-HS to those of the UniqueEOPL instance, and certificates for the violations can be recovered from additional processing. When a violation of weak general position is witnessed on a vertex that lies on the canonical path, a hyperplane incident to dd colors may contain additional points. This in turn implies that some α\alpha-cut is missing, so that the correct solution for the target may not exist. In addition, the (highest-ranked) points x,yx,y from the missing color that we choose to define the orientation of a non-colorful hyperplane may form a segment xyxy that passes through the (d−2)(d-2)-flat spanned by the anchor. In that case the orientation of the hyperplane is not well-defined. In the reduction, these problematic vertices are removed from the canonical path, thereby creating some additional starting points and end points in the reduced instance. These violations can be captured by (U1) with a wrong α\alpha-vector or (UV2). Furthermore, the hyperplanes that contains the degenerate point sets could be represented by different choices of anchors and a additional point on the plane. Each such pair represents a vertex in the reduced instance. We join these vertices in the form of a cycle in the UniqueEOPL instance with all vertices having the same potential value, so that the violations can also be captured by (UV1) and (UV3).

When a violation of well-separation is witnessed on a vertex on the canonical path, the orientations of the two hyperplanes paired by Procedure NextRotateNextRotate may be inconsistent, which may not guarantee that the α\alpha-vector is incremented in one component by one (See Figure 4). Hence, the canonical path is split into two paths that can be captured by (UV2). Furthermore, a violation of well-separation also creates multiple colorful hyperplanes with the same α\alpha-vector (See Figure 4, left). Two vertices in the UniqueEOPL graph with the same potential value, which could correspond to some colorful or non-colorful hyperplanes, can be reported by (UV3). We show that this gives a certificate of violation of well-separation in the following lemmata, where m0m_{0} is the number of bits used to represent each coordinate of points of PP.

Let p1∈P1,…,pd∈Pdp_{1}\in P_{1},\ldots,p_{d}\in P_{d} denote the colorful points on HpH_{p} and q1∈P1,…,qd∈Pdq_{1}\in P_{1},\ldots,q_{d}\in P_{d} denote the colorful points on HqH_{q}. Throughout this proof, we consider H+H^{+} to be a closed halfspace while H−H^{-} is an open halfspace. We prove the claim by induction on the dimension dd.

Without loss of generality, suppose that the open segment p1q1∈Hp−∩Hq+p_{1}q_{1}\in H_{p}^{-}\cap H_{q}^{+}. Since p1∈Hp+∩Hq+p_{1}\in H_{p}^{+}\cap H_{q}^{+} and q1∈Hp−∩Hq+q_{1}\in H_{p}^{-}\cap H_{q}^{+}, there exists at least one point r1∈P1∩(Hp+∩Hq−)r_{1}\in P_{1}\cap(H_{p}^{+}\cap H_{q}^{-}) in order for ∣P1∩Hp+∣=∣P1∩Hq+∣|P_{1}\cap H_{p}^{+}|=|P_{1}\cap H_{q}^{+}| to hold. If the open segment p2q2p_{2}q_{2} also lies in Hp−∩Hq+H_{p}^{-}\cap H_{q}^{+} (resp. Hp+∩Hq−H_{p}^{+}\cap H_{q}^{-}), then there exists at least one point r2r_{2} in P2∩(Hp+∩Hq−)P_{2}\cap(H_{p}^{+}\cap H_{q}^{-}) (resp. P2∩(Hp−∩Hq+)P_{2}\cap(H_{p}^{-}\cap H_{q}^{+})). We can see that the intersection point xx of HpH_{p} and HqH_{q} lies inside the triangles △p1q1r1\triangle p_{1}q_{1}r_{1} and △p2q2r2\triangle p_{2}q_{2}r_{2} (see Figure 11(b)).

Suppose that the open segment p2q2p_{2}q_{2} lies in Hp+∩Hq+H_{p}^{+}\cap H_{q}^{+} (resp. Hp−∩Hq−H_{p}^{-}\cap H_{q}^{-}). In order to assign correct orientations to HpH_{p} and HqH_{q}, the order in which points of P1P_{1} and P2P_{2} appears on the hyperplanes along any direction must be the same for both. This is only feasible when p2p_{2} lies between p1p_{1} (resp. q1q_{1}) and the intersection point x=Hp∩Hqx=H_{p}\cap H_{q}. Hence, p2p_{2} (resp. q2q_{2}) lies inside the triangles △p1q1r1\triangle p_{1}q_{1}r_{1} (see Figure 11(c)).

The idea is to transform PP to a point set P′P^{\prime}, in which we can find two points from the missing color that can each be moved onto one of the non-colorful hyperplanes. Then, the two non-colorful hyperplanes become colorful hyperplanes in P′P^{\prime} with the same α\alpha-vector so that we are in the setup of Lemma 8 and the claim follows.

Without loss of generality, we assume that the missing color of the two non-colorful hyperplanes is color 1. Let HpH_{p} denote one of the non-colorful hyperplanes that passes through some p2∈P2,…,pd∈Pdp_{2}\in P_{2},\ldots,p_{d}\in P_{d} and HqH_{q} denote another non-colorful hyperplane that passes through some q2∈P2,…,qd∈Pdq_{2}\in P_{2},\ldots,q_{d}\in P_{d}. Recall that they have the same α\alpha-vector and dist-value. Let xp,ypx_{p},y_{p} (resp. xq,yqx_{q},y_{q}) be the highest ranked points of P1P_{1} under ≺1\prec_{1} on either side of HpH_{p} (resp. HqH_{q}) and let zpz_{p} (resp. zqz_{q}) be the intersection of segment xpypx_{p}y_{p} (resp. xqyqx_{q}y_{q}) with HpH_{p} (resp. HqH_{q}). By the definition and assumption of dist-value, the dist-value of HpH_{p} and HqH_{q} is ∣∣xp−zp∣∣=∣∣xq−zq∣∣||x_{p}-z_{p}||=||x_{q}-z_{q}||.

In the following, we consider the case of xp,xq∈Hp+∩Hq+x_{p},x_{q}\in H_{p}^{+}\cap H_{q}^{+} with three sub-cases:

[yp∈Hp−∩Hq−\mboxandyq∈Hp+∩Hq−][y_{p}\in H_{p}^{-}\cap H_{q}^{-}\mbox{ and }y_{q}\in H_{p}^{+}\cap H_{q}^{-}]: the argument is symmetrical to the case above.

[yp∈Hp−∩Hq+\mboxandyq∈Hp+∩Hq−][y_{p}\in H_{p}^{-}\cap H_{q}^{+}\mbox{ and }y_{q}\in H_{p}^{+}\cap H_{q}^{-}]: we move ypy_{p} towards xpx_{p} along segment xpypx_{p}y_{p} until it hits HpH_{p} and move yqy_{q} towards xqx_{q} along segment xqyqx_{q}y_{q} until it hits HqH_{q}.

Next, we consider xp∈Hp+∩Hq+x_{p}\in H_{p}^{+}\cap H_{q}^{+} and xq∈Hp−∩Hq+x_{q}\in H_{p}^{-}\cap H_{q}^{+}.

[yp∈Hp−∩Hq+\mboxandyq∈Hp−∩Hq−][y_{p}\in H_{p}^{-}\cap H_{q}^{+}\mbox{ and }y_{q}\in H_{p}^{-}\cap H_{q}^{-}]: we have xq=ypx_{q}=y_{p}. We move xpx_{p} towards xqx_{q} along segment xpxqx_{p}x_{q} until it hits HpH_{p} and move xqx_{q} towards yqy_{q} along segment xqyqx_{q}y_{q} until it hits HqH_{q}.

[yp∈Hp−∩Hp−\mboxandyq∈Hp+∩Hq−][y_{p}\in H_{p}^{-}\cap H_{p}^{-}\mbox{ and }y_{q}\in H_{p}^{+}\cap H_{q}^{-}]: we move xpx_{p} towards xqx_{q} along segment xpxqx_{p}x_{q} until it hits HpH_{p} and move xqx_{q} towards ypy_{p} along segment xqypx_{q}y_{p} until it hits HqH_{q}.

[yp∈Hp−∩Hq+\mboxandyq∈Hp+∩Hq−][y_{p}\in H_{p}^{-}\cap H_{q}^{+}\mbox{ and }y_{q}\in H_{p}^{+}\cap H_{q}^{-}]: we have xq=ypx_{q}=y_{p}. We move ypy_{p} towards xpx_{p} along segment xpypx_{p}y_{p} until it hits HpH_{p} and move yqy_{q} towards xpx_{p} along segment xpyqx_{p}y_{q} until it hits HqH_{q}.

[yp∈Hp−∩Hq−\mboxandyq∈Hp−∩Hq−][y_{p}\in H_{p}^{-}\cap H_{q}^{-}\mbox{ and }y_{q}\in H_{p}^{-}\cap H_{q}^{-}]: we have yp=yqy_{p}=y_{q}. We move xpx_{p} towards xqx_{q} along segment xpxqx_{p}x_{q} until it hits HpH_{p} and move xqx_{q} towards yqy_{q} along segment xqyqx_{q}y_{q} until it hits HqH_{q}.

The case for xp∈Hp+∩Hq−x_{p}\in H_{p}^{+}\cap H_{q}^{-} and xq∈Hp+∩Hq+x_{q}\in H_{p}^{+}\cap H_{q}^{+} are symmetrical to those above.

The last case is xp∈Hp+∩Hq−x_{p}\in H_{p}^{+}\cap H_{q}^{-} and xq∈Hp−∩Hq+x_{q}\in H_{p}^{-}\cap H_{q}^{+}. Basically, the sub-cases are the same as those above except one, which happens for yp∈Hp−∩Hp+y_{p}\in H_{p}^{-}\cap H_{p}^{+} and yq∈Hp+∩Hq−y_{q}\in H_{p}^{+}\cap H_{q}^{-}. In this case, we have xp=yqx_{p}=y_{q} and xq=ypx_{q}=y_{p}. We move xpx_{p} and xqx_{q} towards each other along segment xpxqx_{p}x_{q} and until they hit HpH_{p} and HqH_{q}. Note that the segment xpxqx_{p}x_{q} may intersect the (d−2)(d-2)-flat Hp∩HqH_{p}\cap H_{q}, but this case is also handled by Lemma 8. ∎

Alpha-HS ∈\in UEOPL\mathsf{UEOPL} ⊆\subseteq CLS\mathsf{CLS}.

Double-wedges and the neighborhood graph

In this section we formally define the notion of double-wedges and the underlying graph that is defined using rotations.

An anchor is a colorful (d−1)(d-1)-tuple of PP which spans a (d−2)(d-2)-flat. Let PiP_{i} denote the missing color in the anchor. Then the tuple for the anchor R=(p1,…,pd−1)R=(p_{1},\dots,p_{d-1}) is ordered as R={p1∈P1,…,pi−1∈Pi−1,pi∈Pi+1,…,pd−1∈Pd}R=\{p_{1}\in P_{1},\dots,p_{i-1}\in P_{i-1},p_{i}\in P_{i+1},\dots,p_{d-1}\in P_{d}\}. An anchor RR along with a pair of points p,q∈Pp,q\in P such that p,q∉Rp,q\not\in R is called a double-wedge if all of the following hold:

the hyperplane HpH_{p} through R∪{p}R\cup\{p\} does not contain qq. This implies that the hyperplane HqH_{q} through R∪{q}R\cup\{q\} does not contain pp.

if x,yx,y are the highest ordered points of PiP_{i} under ≺i\prec_{i} on either sides of Hp,HqH_{p},H_{q}, then R∪{x,y}R\cup\{x,y\} does not lie on a hyperplane.

the intersection of the open halfspaces Hq+∩Hp−H_{q}^{+}\cap H_{p}^{-} is empty, that is, it does not contain any point of PP. Similarly Hp+∩Hq−H_{p}^{+}\cap H_{q}^{-} must also be empty.

We visualize the anchor as a (d−2)(d-2)-ridge through which Hp,HqH_{p},H_{q} pass through. A rotation around the anchor changes one hyperplane to the other without passing through any other point of PP. Intuitively the double-wedge refers to the space (Hq+∩Hp−)∪(Hp+∩Hq−)(H_{q}^{+}\cap H_{p}^{-})\cup(H_{p}^{+}\cap H_{q}^{-}) and we use this interpretation several times. See Figure 5 for a simple example.

For a double-wedge w:=(R,p,q)w:=(R,p,q), we define a representative hyperplane HwH_{w} as the hyperplane that is the angular bisector of the double-wedge. Since a double-wedge is empty, HwH_{w} does not contain any point of PP apart from RR. We define an orientation for HwH_{w} based on RR and the color missing from RR. Let x,yx,y be points from the missing color as defined before. Without loss of generality, let x∈Hw+x\in H_{w}^{+}, and y∈Hw−y\in H_{w}^{-}. We call the first hyperplane among Hp,HqH_{p},H_{q} that intersects the directed segment xyxy as the upper hyperplane of ww and the other hyperplane as the lower hyperplane of ww. A simple example can be found in Figure 6.

The α\alpha-vector of any oriented hyperplane HH is a dd-tuple (a1,…,ad)(a_{1},\dots,a_{d}) of integers where aia_{i} is the number of points of PiP_{i} in the closed halfspace H+H^{+} for i∈[d]i\in[d]. The α\alpha-vector of a double-wedge w=(R,p,q)w=(R,p,q) is defined as the α\alpha-vector of its representative hyperplane. We say that a double-wedge w=(R,p,q)w=(R,p,q) is non-colorful, if both R∪{p}R\cup\{p\} and R∪{q}R\cup\{q\} are non-colorful, and colorful, if exactly one of R∪{p}R\cup\{p\} and R∪{q}R\cup\{q\} is colorful, and very colorful, if both R∪{p}R\cup\{p\} and R∪{q}R\cup\{q\} are colorful.

Under the assumption of weak general position, we additionally have that if ww is non-colorful, then Hp,HqH_{p},H_{q} are non-colorful, and if ww is colorful, exactly one of Hp,HqH_{p},H_{q} is colorful, and if ww is very colorful, both Hp,HqH_{p},H_{q} are colorful.

The definition of dist-value for hyperplanes can be extended to double-wedges by setting the dist-value of a double-wedge as that of its representative hyperplane. Consequently, the results of Lemmas 5, 6 and 7 extend to double-wedges with simple modifications.

2 Defining a neighborhood graph

We define a concept of neighborhood between double-wedges, and then we use this to define a graph whose vertices correspond to the double-wedges. We first describe the graph under the assumptions that the colors are well-separated and PP is in weak general position. Later we show how to handle the cases when these assumptions fail.

We call two double-wedges (R,p,q),(R′,p′,q′)(R,p,q),(R^{\prime},p^{\prime},q^{\prime}) neighboring if both share a common hyperplane, that is, {Hp,Hq}∩{Hp′,Hq′}≠∅\{H_{p},H_{q}\}\cap\{H_{p^{\prime}},H_{q^{\prime}}\}\neq\varnothing, with an exception that we elaborate below. A double-wedge w=(R,p,q)w=(R,p,q) has different number of neighbors depending on how colorful its hyperplanes are. The anchor can be written in the form R={x1,…,xd−1}R=\{x_{1},\dots,x_{d-1}\}.

Let ww be non-colorful. Then p,qp,q both share their colors with those of RR. Suppose pp has the same color as xix_{i}. Then there are at most three neighboring double-wedges that share HpH_{p}. One of them use the same anchor RR, and as an exception we do not count this as a neighboring double-wedge. For the two remaining neighbors, the anchor is R′=(x1,…,p,…,xd−1)R^{\prime}=(x_{1},\dots,p,\dots,x_{d-1}) where pp has replaced xix_{i}. The two rotational directions determine the two double-wedges. Only one of them has the same α\alpha-vector as ww, since the representative hyperplanes contain xix_{i} on opposite sides. With a similar argument, there are at most two neighboring double-wedges that share HqH_{q} and at most one of them has the same α\alpha-vector as ww.

Let ww be colorful, where HpH_{p} is colorful and HqH_{q} is non-colorful, without loss of generality. By replacing some xix_{i} by pp we get an anchor that is contained in HpH_{p} and which may define a double-wedge for each of the two rotational directions. Since there are (d−1)(d-1) possible anchors formed by replacement, there are at most 2(d−1)2(d-1) double-wedges that share HpH_{p}. Additionally, keeping the anchor RR fixed, there is at most one neighboring double-wedge. So there are at most 2d−12d-1 neighboring double-wedges of ww that share HpH_{p}. The case for HqH_{q} is similar to case (1).

Let ww be very colorful. Similar to case (2), there are at most 2d−12d-1 double-wedges sharing HpH_{p}. The case for HqH_{q} is similar.

We build a graph GG where each vertex corresponds to a double-wedge. Let w=(R,p,q)w=(R,p,q) be any double-wedge. For simplicity, we denote the vertex in GG corresponding to ww also by ww. If HpH_{p} is colorful, we add an edge in GG between ww and the vertex of each neighboring double-wedge that shares HpH_{p}. If HpH_{p} is non-colorful, we add an edge only with the vertex of the double-wedge that shares its α\alpha-vector with that of ww. Thus, non-colorful double-wedges have degree two in GG, while colorful and very colorful double-wedges have degrees at most 2d2d and 4d−24d-2, respectively.

We transfer each attribute of a double-wedge to its vertex in GG. For instance, we call vertices of GG as non-colorful, colorful or very colorful corresponding to the color of the double-wedge representing the vertex. Similarly, each vertex has an α\alpha-vector that corresponds to the α\alpha-vector of its double-wedge, and so on. See Figure 7 for an elementary example.

Let v∈Gv\in G be any vertex and let (α1,…,αd)(\alpha_{1},\dots,\alpha_{d}) denote the α\alpha-vector of vv. The largest α\alpha-vector for any hyperplane is (n1,…,nd)(n_{1},\dots,n_{d}) that occurs on a unique tangent hyperplane whose half-space contains PP. With our definition of the α\alpha-vector of double-wedges using the representative hyperplanes, for any double-wedge ww, the α\alpha-vector of ww is smaller in at least one coordinate from the maximum.

Let HH be a colorful hyperplane with α\alpha-vector (α1,…,αd)(\alpha_{1},\dots,\alpha_{d}). For each j∈[d]j\in[d], if αj≤nj−1\alpha_{j}\leq n_{j}-1 (resp. αj≥2\alpha_{j}\geq 2), then there is a path w,w1,w2,…,wk,w′w,w_{1},w_{2},\dots,w_{k},w^{\prime} in GG such that the double-wedge ww is incident to HH, the double-wedge w′w^{\prime} is incident to another colorful hyperplane H′H^{\prime}, w,w1,w2,…,wkw,w_{1},w_{2},\dots,w_{k} share the same α\alpha-vector and the α\alpha-vector for H′H^{\prime} differs only in the jj-th component, where the value is αj+1\alpha_{j}+1 (resp. αj−1\alpha_{j}-1).

For the case that αj≤nj−1\alpha_{j}\leq n_{j}-1, we set an anchor RR in HH that excludes the jj-th colored point, say pjp_{j}. Then, we apply Procedure NextRotateNextRotate starting from (R,pj)(R,p_{j}) until we get another colorful hyperplane H′H^{\prime}. Let ww be the double-wedge created by the first rotation. Note that pjp_{j} is on the upper hyperplane of ww. During this sequence of rotations, we also get a sequence of double-wedges. Before the rotating hyperplane H′′H^{\prime\prime} hits a point pip_{i} of repeated color ii, assume that pip_{i} is in the negative half-space of H′′H^{\prime\prime}. Once pip_{i} is on H′′H^{\prime\prime}, we swap pip_{i} with another point pi′p^{\prime}_{i} of the same color in the anchor and keep the rotation towards the opposite direction of the orientation so that pip_{i} is in the positive half-space of H′′H^{\prime\prime} and pi′p^{\prime}_{i} from the positive half-space moves to the negative half-space. This is true because the orientation is consistent by Lemma 5. If pip_{i} is in the positive half-space of H′′H^{\prime\prime} before pip_{i} is hit by H′′H^{\prime\prime}, then pip_{i} remains in the positive half-space and pi′p^{\prime}_{i} as well (see Figure 3). Both cases maintain αi\alpha_{i} during the rotation. Thus, all non-colorful double-wedges in this sequence of rotations have the same α\alpha-vector as ww. In the last rotation, the rotating hyperplane hits the first point pj′p^{\prime}_{j} of color jj. By Lemma 7, well-separation guarantees the consistency of the orientation of the rotating hyperplane so that pj′p^{\prime}_{j} moves from the negative half-space to the positive half-space and other points of color jj remain in the same sides of the hyperplane. Thus, αj\alpha_{j} is increased by one. The same argument also works for the case that αj≥2\alpha_{j}\geq 2 by using the inverse of Procedure NextRotateNextRotate.

Since all the double-wedges created by the first rotation for each jj are incident to HH, they are also connected in GG by definition. ∎

By Lemma 12, we know that all (very) colorful double-wedges are connected. For non-colorful double-wedges, we apply Procedure NextRotateNextRotate on its lower hyperplane until the rotating hyperplane hits some point of the missing color, which implies that non-colorful double-wedges also connect to some (very) colorful double-wedge. ∎

The neighborhood graph GG imitates a grid in a coarse sense. There is a ”vertex” for every colorful dd-tuple of PP, and there are paths connecting these grid vertices. We showed in Lemma 13 that GG is connected. Therefore, given a target α\alpha-vector (α1,…,αd)(\alpha_{1},\dots,\alpha_{d}), the correct dd-tuple can be found by starting from some vertex and walking towards the solution. See Figure 8 for an illustration.

If PP violates well-separation or weak general position, then many nice properties of the neighborhood graph are destroyed. Double-wedges may fail to have consistent orientations by Lemma 7. There may be multiple solutions for the same α\alpha-cut, and no solutions for other cuts. The former case will manifest as multiple vertices with the same α\alpha-vector, but they may lie in different connected components, so Lemma 13 will fail, making the graph disconnected. For the latter case, there will be no vertex in GG that corresponds to the α\alpha-cut. The grid-like structure exhibited in Lemma 12 is also not applicable anymore, meaning that the canonical path may not exist. See Figure 9 for a graph that contains violations.

The formal reduction

As shown in Section 3, a vertex vv in I′\cal I^{\prime} corresponds to (R,p,q)(R,p,q) in I\cal I, where RR is a colorful point set of size (d−1)(d-1) from PP, and p,q∈Pp,q\in P. We are only interested in the case when (R,p,q)(R,p,q) is a double-wedge, as per the definition in Section 4.1. Otherwise, we create a self loop on vv in I′\cal{I^{\prime}}. Furthermore, if there are no violations in I\cal I, we can define a canonical path from the vertex v0v_{0} with α\alpha-vector =(1,…,1)=(1,\ldots,1) to the unique vertex vαv_{\alpha} with α\alpha-vector =(α1,…,αd)=(\alpha_{1},\dots,\alpha_{d}) (shown in Section 3.1), which is the unique path in I′\cal I^{\prime}. For other vertices vv not on the path, we also create a self loop on vv. For instance, when I\cal I fails weak general position assumption, then R∪{p1,p2,…,pm}R\cup\{p_{1},p_{2},\ldots,p_{m}\} lie on the same hyperplane, where p1≺⋯≺pm∈Pp_{1}\prec\dots\prec p_{m}\in P. In this case we create a cycle on v1=(R,p1,q),v2=(R,p2,q),…,vm=(R,pm,q)v_{1}=(R,p_{1},q),v_{2}=(R,p_{2},q),\ldots,v_{m}=(R,p_{m},q) with the same potential value on each viv_{i}, so that this violation may be reported as the violation (UV1) in I′\cal I^{\prime}. When I\cal I fails the well-separated assumption, the graph we constructed may contain more than one path, which may be reported as the violations (UV2) or (UV3) in I′\cal I^{\prime}. In particular, if a hyperplane witnesses both the violations of weak general position and well-separation, then the cycle may become a path with the same potential value, so any violation could be possible.

t1t_{1} contains ⌈log⁡d⌉\lceil\log d\rceil bits representing the index of the missing color in RR,

t2,…,tdt_{2},\ldots,t_{d} each contain ⌈log⁡n0⌉\lceil\log n_{0}\rceil bits for the indices of the points in RR ordered by ≺\prec,

(td+1,td+2)(t_{d+1},t_{d+2}) contain ⌈log⁡d⌉\lceil\log d\rceil and ⌈log⁡n0⌉\lceil\log n_{0}\rceil bits for the index of the color and the index of p∈Ptd+1p\in P_{t_{d+1}} respectively,

and we use the same idea to represent qq by (td+3,td+4)(t_{d+3},t_{d+4}). Altogether, we need at most κ=3⋅⌈log⁡d⌉+(d+1)⋅⌈log⁡n0⌉\kappa=3\cdot\lceil\log d\rceil+(d+1)\cdot\lceil\log n_{0}\rceil bits in the encoding. Let fv:{P}d+1→{0,1}κf_{v}:\{P\}^{d+1}\rightarrow\{0,1\}^{\kappa} denote the function that encodes (R,p,q)(R,p,q) to (t1,…,td+4)(t_{1},\ldots,t_{d+4}). If (R,p,q)(R,p,q) is a double-wedge, then (R,q,p)(R,q,p) also represents the same double-wedge. Since we do not want to create two valid vertices in GG corresponding to the same double-wedge, we only pick the one with (td+1,td+2)(t_{d+1},t_{d+2}) on the upper hyperplane as a valid double-wedge. The following lemma details how we can verify whether a given encoding (t1,…,td+4)(t_{1},\ldots,t_{d+4}) is a double-wedge:

The rotation process is handled by procedures NextNeighborNextNeighbor and PrevNeighborPrevNeighbor, which we describe next. Let w=(R,p,q)w=(R,p,q) be the current double-wedge. Some abnormal cases may happen in the output w′=(R′,p′,q′)w^{\prime}=(R^{\prime},p^{\prime},q^{\prime}) of NextNeighborNextNeighbor or PrevNeighborPrevNeighbor when the segment xyxy that defines the orientation of Hw′H_{w^{\prime}} passes through R′R^{\prime} (see Figure 10), Hp′H_{p^{\prime}} or Hq′H_{q^{\prime}} contains more than dd points, or the orientations of Hp′,Hq′,Hw′H_{p^{\prime}},H_{q^{\prime}},H_{w^{\prime}} are not consistent. For these cases, the path will end or start at ww. When UniqueEOPL outputs ww, we can compute w′w^{\prime} and find the certificate of a violation as follows: for the first or second case, it is easy to see that it violates well-separation and/or weak general position. For the third case, we show that it violates well-separation and that we can obtain a certificate for the violation in Lemma 7.

Procedure (R′,p′,q′)=NextNeighbor(R,p,q)(R^{\prime},p^{\prime},q^{\prime})=NextNeighbor(R,p,q)

Let w=(R,p,q)w=(R,p,q) and let t1t_{1} be the missing color in RR.

Let Pt1+P^{+}_{t_{1}} be the subset of Pt1P_{t_{1}} that lies in Hw+H^{+}_{w} and Pt1−P^{-}_{t_{1}} be the subset of Pt1P_{t_{1}} that lies in Hw−H^{-}_{w}. Then, let x∈Pt1+x\in P^{+}_{t_{1}} be the highest ranked point according to the order ≺t1\prec_{t_{1}} and y∈Pt1−y\in P^{-}_{t_{1}} be the highest ranked point according to the order ≺t1\prec_{t_{1}}. As we describe previously, qq lies on the lower hyperplane.

Since ww is supposed to be on the canonical path and is not the end point, the α\alpha-vector of HqH_{q} is in the form of (α1,…,αt1−1,bt1,1…,1)(\alpha_{1},\ldots,\alpha_{t_{1}-1},b_{t_{1}},1\ldots,1) with bt1≤αt1b_{t_{1}}\leq\alpha_{t_{1}}.

If qq shares the same color of a point rr in RR, then set R′:=R∪{q}∖{r}R^{\prime}:=R\cup\{q\}\setminus\{r\} and p′:=rp^{\prime}:=r.

If qq is in color t1t_{1} and bt1<αt1b_{t_{1}}<\alpha_{t_{1}}, then set R′:=RR^{\prime}:=R and p′:=qp^{\prime}:=q.

If qq is in color t1t_{1} and bt1=αt1b_{t_{1}}=\alpha_{t_{1}}, then let rr be the point in RR with color t1+1t_{1}+1 and set R′:=R∪{q}∖{r}R^{\prime}:=R\cup\{q\}\setminus\{r\} and p′=rp^{\prime}=r.

We rotate HqH_{q} around the anchor R′R^{\prime} in a direction such that the hyperplane is moving away from xx along the segment xyxy until it hits a point q′∈Pq^{\prime}\in P.

Return (R′,p′,q′)(R^{\prime},p^{\prime},q^{\prime}).

Procedure (R′,p′,q′)=PrevNeighbor(R,p,q)(R^{\prime},p^{\prime},q^{\prime})=PrevNeighbor(R,p,q)

Let w=(R,p,q)w=(R,p,q) and let t1t_{1} be the missing color in RR.

Let Pt1+P^{+}_{t_{1}} be the subset of Pt1P_{t_{1}} that lies in Hw+H^{+}_{w} and Pt1−P^{-}_{t_{1}} be the subset of Pt1P_{t_{1}} that lies in Hw−H^{-}_{w}. Then, let x∈Pt1+x\in P^{+}_{t_{1}} be the highest ranked point according to the order ≺t1\prec_{t_{1}} and y∈Pt1−y\in P^{-}_{t_{1}} be the highest ranked point according to the order ≺t1\prec_{t_{1}}. As we describe previously, pp lies on the upper hyperplane.

Since ww is supposed to be on the canonical path and is not the starting point, the α\alpha-vector of HpH_{p} is in the form of (α1,…,αt1−1,bt1,1…,1)(\alpha_{1},\ldots,\alpha_{t_{1}-1},b_{t_{1}},1\ldots,1) with bt1≤αt1b_{t_{1}}\leq\alpha_{t_{1}}. When t1=1t_{1}=1, b1>1b_{1}>1.

If pp shares the same color of a point rr in RR, then set R′:=R∪{p}∖{r}R^{\prime}:=R\cup\{p\}\setminus\{r\} and q′:=rq^{\prime}:=r.

If pp is in color t1t_{1} and bt1>1b_{t_{1}}>1, then set R′:=RR^{\prime}:=R and q′:=pq^{\prime}:=p.

If pp is in color t1t_{1} and bt1=1b_{t_{1}}=1, then let rr be the point in RR with color t1−1t_{1}-1 and set R′:=R∪{p}∖{r}R^{\prime}:=R\cup\{p\}\setminus\{r\} and q′=rq^{\prime}=r.

We rotate HpH_{p} around the anchor R′R^{\prime} in a direction such that the hyperplane is moving closer to xx along the segment xyxy until it hits a point p′∈Pp^{\prime}\in P.

Return (R′,p′,q′)(R^{\prime},p^{\prime},q^{\prime}).

The points xx and yy from the missing color can be found in linear time. We can also check which point will hit the hyperplane first during the rotation by a prune-and-search technique in polynomial time. ∎

If v=0κv=0^{\kappa}, then Return fv({p2′,…,pd′},p1′,q′)f_{v}(\{p^{\prime}_{2},\ldots,p^{\prime}_{d}\},p^{\prime}_{1},q^{\prime}), where {p1′,…,pd′}\{p^{\prime}_{1},\ldots,p^{\prime}_{d}\} is a colorful point set on a hyperplane H0H_{0} that has α\alpha-vector =(1,…,1)=(1,\ldots,1), and q′q^{\prime} is the point that creates a double-wedge with the anchor {p2′,…,pd′}\{p^{\prime}_{2},\ldots,p^{\prime}_{d}\} and p1′p^{\prime}_{1}.

If (t1,…,td+4)(t_{1},\ldots,t_{d+4}) is not the encoding of a double-wedge, then Return (t1,…,td+4)(t_{1},\ldots,t_{d+4}).

Let (R,p,q)(R,p,q) denote the double-wedge for which fv(R,p,q)=(t1,…,td+4)f_{v}(R,p,q)=(t_{1},\ldots,t_{d+4}). If the orientations of HwH_{w}, HpH_{p} and HqH_{q} are not consistent, then Return fv(R,p,q)f_{v}(R,p,q).

Recall that t1t_{1} is the missing color in RR. If the α\alpha-vector of (R,p,q)(R,p,q) is not in the form of (α1,…,αt1−1,bt1,1,…,1)(\alpha_{1},\ldots,\alpha_{t_{1}-1},b_{t_{1}},1,\ldots,1) with bt1<αt1b_{t_{1}}<\alpha_{t_{1}}, then Return fv(R,p,q)f_{v}(R,p,q).

If the α\alpha-vector of the lower hyperplane HqH_{q} is (α1,…,αd)(\alpha_{1},\ldots,\alpha_{d}), then Return fv(R,p,q)f_{v}(R,p,q).

If HpH_{p} contains some points of PP other than R∪{p}R\cup\{p\}, then let p1,p2,…,pmp_{1},p_{2},\ldots,p_{m} be those extra points ordered by ≺\prec, and let pip_{i} be the point just after pp in the order ≺\prec (if pp is after pmp_{m}, then pi=p1p_{i}=p_{1}). We Return fv(R,pi,q)f_{v}(R,p_{i},q).

Similarly, Return fv(R,p,qi)f_{v}(R,p,q_{i}) if HqH_{q} contains some points other than R∪{q}R\cup\{q\}, where qiq_{i} is the point just after qq according to the order ≺\prec among those extra points.

Let w′=(R′,p′,q′)w^{\prime}=(R^{\prime},p^{\prime},q^{\prime}) be the output of NextNeighbor(R,p,q)NextNeighbor(R,p,q). If the orientations of Hw′H_{w^{\prime}}, Hp′H_{p^{\prime}} and Hq′H_{q^{\prime}} are well-defined and consistent, then Return fv(R′,p′,q′)f_{v}(R^{\prime},p^{\prime},q^{\prime}).

If v=0κv=0^{\kappa}, then Return 0κ0^{\kappa}.

If (t1,…,td+4)(t_{1},\ldots,t_{d+4}) is not a double-wedge, then Return (t1,…,td+4)(t_{1},\ldots,t_{d+4}).

Let (R,p,q)(R,p,q) be a double-wedge such that fv(R,p,q)=(t1,…,td+4)f_{v}(R,p,q)=(t_{1},\ldots,t_{d+4}). If the orientations of HwH_{w}, HpH_{p} and HqH_{q} are not consistent, then Return fv(R,p,q)f_{v}(R,p,q).

If R={p2′,…,pd′}R=\{p^{\prime}_{2},\ldots,p^{\prime}_{d}\} and p=p1′p=p^{\prime}_{1}, then Return 0κ0^{\kappa}.

If HpH_{p} is colorful, t1=1t_{1}=1 and the α\alpha-vector of HpH_{p} is (1,…,1)(1,\ldots,1), then Return fv(R,p,q)f_{v}(R,p,q).

Recall that t1t_{1} is the missing color in RR. If the α\alpha-vector of (R,p,q)(R,p,q) is not in the form of (α1,…,αt1−1,bt1,1,…,1)(\alpha_{1},\ldots,\alpha_{t_{1}-1},b_{t_{1}},1,\ldots,1) with bt1<αt1b_{t_{1}}<\alpha_{t_{1}}, then Return fv(R,p,q)f_{v}(R,p,q).

If HpH_{p} contains some points of PP other than R∪{p}R\cup\{p\}, then let p1,p2,…,pmp_{1},p_{2},\ldots,p_{m} be those extra points ordered by ≺\prec, and let pip_{i} be the point just before pp in ≺\prec (if pp is before p1p_{1}, then pi=pmp_{i}=p_{m}), and Return fv(R,pi,q)f_{v}(R,p_{i},q).

Similarly, Return fv(R,p,qi)f_{v}(R,p,q_{i}) if HqH_{q} contains some points other than R∪{q}R\cup\{q\}, where qiq_{i} is the point just before qq in ≺\prec among those extra points.

Let w′=(R′,p′,q′)w^{\prime}=(R^{\prime},p^{\prime},q^{\prime}) be the output of PrevNeighbor(R,p,q)PrevNeighbor(R,p,q). If the orientations of Hw′H_{w^{\prime}}, Hp′H_{p^{\prime}} and Hq′H_{q^{\prime}} are well-defined and consistent, then Return fv(R′,p′,q′)f_{v}(R^{\prime},p^{\prime},q^{\prime}).

Given a double-wedge (R,p,q)(R,p,q), let x′x^{\prime} be any point of PP in Hp+∩Hq+H^{+}_{p}\cap H^{+}_{q} and let y′y^{\prime} be any point of PP in Hp−∩Hq−H^{-}_{p}\cap H^{-}_{q}. Suppose that the segment x′y′x^{\prime}y^{\prime} does not pass through the affine hull of RR so that the intersections of x′y′x^{\prime}y^{\prime} with HpH_{p} and HqH_{q} are two distinct points. Let dmin⁡d_{\min} be the Euclidean distance between the two intersection points of x′y′x^{\prime}y^{\prime} with HpH_{p} and HqH_{q}. Let dmax⁡d_{\max} be the Euclidean distance between xx and the intersection point of x′y′x^{\prime}y^{\prime} with HqH_{q}.

Let m0m_{0} denote the number of bits needed to represent each coordinate of any point of PP. Then, dmin⁡d_{\min} is at less 1/N21/N^{2} and dmax⁡d_{\max} is at most MM, where N=d!2dm0N=d!2^{dm_{0}} and M=d2m0M=\sqrt{d}2^{m_{0}}.

Without loss of generality, we assume that the missing color in RR is color 1, i.e., RR is a set of points {p2∈P2,…,pd∈Pd}\{p_{2}\in P_{2},\ldots,p_{d}\in P_{d}\}. Let zpz_{p} (resp. zqz_{q}) be the intersection of x′y′x^{\prime}y^{\prime} and HpH_{p} (resp. HqH_{q}). Since zpz_{p} lies on x′y′x^{\prime}y^{\prime} and the affine hull of HpH_{p}, we can represent zpz_{p} as the convex combination of x′x^{\prime} and y′y^{\prime} and the linear combination of R∪{p}R\cup\{p\} as follows:

where y′−x′,p2−p1,…y^{\prime}-x^{\prime},p_{2}-p_{1},\dots are column vectors. Since we assume that zpz_{p} exists, we have det⁡A≠0\det A\neq 0. According to Cramer’s rule, we have λi=det⁡Aidet⁡A\lambda_{i}=\frac{\det A_{i}}{\det A}, where AiA_{i} is the matrix obtained by replacing the ii-th column of AA with bb. Using Leibniz formula for determinants, we can bound the denominator:

For dmax⁡d_{\max}, it is less than the longest possible line segment, so dmax⁡≤d2m0d_{\max}\leq\sqrt{d}2^{m_{0}}. ∎

If (t1,…,td+4)(t_{1},\ldots,t_{d+4}) is not a double-wedge, then Return . Otherwise, let (R,p,q)(R,p,q) be a double-wedge such that fv(R,p,q)=(t1,…,td+4)f_{v}(R,p,q)=(t_{1},\ldots,t_{d+4}).

The following lemma shows that if there are no violations in Alpha-HS, there are also no violations in the constructed UniqueEOPL instance. This makes the reduction promise-preserving. In particular, we can find the (α1,…,αd)(\alpha_{1},\ldots,\alpha_{d})-cut from the unique solution of the UniqueEOPL instance.

If there are no violations in I\cal I, then the constructed UniqueEOPL instance I′\cal I^{\prime} only contains a type (U1) solution whose lower hyperplane is colorful and has α\alpha-vector =(α1,…,αd)=(\alpha_{1},\ldots,\alpha_{d}), which is a type (G1) solution of I\cal I.

There are different types of violations that may return from a UniqueEOPL instance. To obtain the certificate of violation (GV1) in Alpha-HS, we need the following lemmas to convert the violation solutions of UniqueEOPL to the (GV1) certificate.

Similarly, for any point yy in the open halfspace H−H^{-},

The proof for points in the open halfspace H−H^{-} is similar. ∎

By Theorem 2, we already know that if there are multiple α\alpha-cuts, then PP is not well-separated. The following two lemmas give the same result, but we provide a constructive proof so that we can find the certificate of the violation.

The proof is based on the idea of finding a hyperplane that ”interpolates” between HpH_{p} and HqH_{q}, for which no consistent orientation can be defined. This hyperplane gives the certificate of violation.

Let p1∈P1,…,pd∈Pdp_{1}\in P_{1},\ldots,p_{d}\in P_{d} denote the colorful points on HpH_{p} and q1∈P1,…,qd∈Pdq_{1}\in P_{1},\ldots,q_{d}\in P_{d} denote the colorful points on HqH_{q}.

The open segment piqip_{i}q_{i} is contained in either Hp+∩Hq+H_{p}^{+}\cap H_{q}^{+} or Hp−∩Hq−H_{p}^{-}\cap H_{q}^{-} for each i∈[d]i\in[d].

We now extend Lemma 8 to the case where the two hyperplanes are non-colorful but have the same α\alpha-vector and dist-value. On any path of non-colorful vertices of GG that starts and ends at colorful vertices, the dist-values of the vertices along the path is bounded by the dist-values of the two end points. Hence, we have an interval of the dist-values along the path with respect to an α\alpha-vector (α1,…,αi−1,bi,1,…,1)(\alpha_{1},\ldots,\alpha_{i-1},b_{i},1,\ldots,1), where ii is the missing color for those non-colorful vertices and bi<αib_{i}<\alpha_{i}. If we find another non-colorful vertex with the same missing color and α\alpha-vector but its dist-value is outside the interval, we show that it implies a violation of well-separation in the following lemma.

Without loss of generality, we assume that the missing color of HwH_{w} and HH is color 11. Let p2∈P2,…,pd∈Pdp_{2}\in P_{2},\ldots,p_{d}\in P_{d} be the colorful points on HH and q2∈P2,…,qd∈Pdq_{2}\in P_{2},\ldots,q_{d}\in P_{d} be the colorful points on HwH_{w}. Let xp,ypx_{p},y_{p} (resp. xq,yqx_{q},y_{q}) be the highest ranked points of P1P_{1} under ≺1\prec_{1} on either side of HH (resp. HwH_{w}) As we saw in the proof of Lemma 9, the condition that the hyperplanes share the dist-value is required only in the first case when xp,xq∈H+∩Hw+x_{p},x_{q}\in H^{+}\cap H_{w}^{+} and yp,yq∈H−∩Hw−y_{p},y_{q}\in H^{-}\cap H_{w}^{-}. We can apply the same analysis from Lemma 9 for HH and HwH_{w} in other cases.

When xp,xq∈H+∩Hw+x_{p},x_{q}\in H^{+}\cap H_{w}^{+} and yp,yq∈H−∩Hw−y_{p},y_{q}\in H^{-}\cap H_{w}^{-}, we have xp=xqx_{p}=x_{q} and yp=yqy_{p}=y_{q}. Since HwH_{w} is the representative hyperplane, HqH_{q} is the lower hyperplane and the dist-value of HH is larger than that of HqH_{q}, the directed segment xqyqx_{q}y_{q} intersects them in the following order: HwH_{w}, HqH_{q} then HH. Since ww is a double-wedge, the orientations of HwH_{w} and HqH_{q} are consistent. That means, q∈Hq∩Hw−q\in H_{q}\cap H_{w}^{-}. There are two possibilities:

q∈H+∩Hw−q\in H^{+}\cap H_{w}^{-}: using a similar argument as in the proof of Lemma 9, there must exist another point r′∈P1∩H−∩Hw+r^{\prime}\in P_{1}\cap H^{-}\cap H_{w}^{+} because q∈P1q\in P_{1} and H,HwH,H_{w} have the same α\alpha-vector. Then, we move qq towards xpx_{p} along segment xpqx_{p}q until it hits HwH_{w} and move r′r^{\prime} towards xpx_{p} along segment xpr′x_{p}r^{\prime} until it hits HH (Figure 12(a)). We apply Lemma 8 to the resulting point set and that completes the proof.

q∈H−∩Hw−q\in H^{-}\cap H_{w}^{-}: in this case both yqy_{q} and qq lie in q∈H−∩Hw−q\in H^{-}\cap H_{w}^{-}, but yq≠qy_{q}\neq q since that would make the intersection order along xqyqx_{q}y_{q} as Hw,H,HqH_{w},H,H_{q}, which is a contradiction. Along the directed segment xqqx_{q}q, the order of intersection is H,HwH,H_{w} (Figure 12(b)) while on the directed segment xqyqx_{q}y_{q} it is Hw,HH_{w},H. So we move yqy_{q} towards xqx_{q} along xqyqx_{q}y_{q} until it hits HH, and move qq towards xqx_{q} along xqqx_{q}q until it hits HwH_{w}. Then we re-use the proof from Lemma 8.

In the following lemmas, we show how to reduce any violation solution of I′\cal I^{\prime} to a violation solution of I\cal I to complete the reduction.

The next lemma is to handle type (UV1) solutions that capture a vertex at which the potential value is not increasing. From the way we construct the graph, it only happens when the weak general position fails.

A solution of type (UV2) means that there is another starting point of some other path. The proof of the following lemma is basically the same as Lemma 22.

For Case 1, we apply the same argument in the proof of Lemma 22 to show that the case is impossible. For Case 2, let w′=(R′,p′,q′)w^{\prime}=(R^{\prime},p^{\prime},q^{\prime}) be the output of PrevNeighbor(v)PrevNeighbor(v). The orientations of w′w^{\prime} are not well-defined or inconsistent, then we can find a type (GV2) solution. ∎

In (UV3), either we have two representative hyperplanes with the same potential value, or we find a hyperplane whose potential value is between the potential values of two consecutive representative hyperplanes in GG. We use Lemmas 9 and 21 to find a violation solution of Alpha-HS.

Conclusion and future work

We gave a complexity-theoretic upper bound for Alpha-HS. No hardness results are known for this search problem, and the next question is determining if this is hard for UEOPL\mathsf{UEOPL}. One challenge is that UniqueEOPL is formulated as Boolean circuits, whereas Alpha-HS is purely geometric. Emulating circuits using purely geometric arguments is highly non-trivial. Filos-Ratsikas and Goldberg showed a reduction of this form in [FRG19]. They reduced the PPA\mathsf{PPA}-complete Tucker circuit to Ham-Sandwich, going via the Consensus-Halving [SS03], and the Necklace-splitting problems [AW86]. It could be a worthwhile exercise to investigate if their techniques can provide insights for hardness of Alpha-HS.

Some related problems are determining the complexity of answering whether a point set is well-separated, whether it is in weak general position, or whether a given α\alpha-cut exists for the point set. A given α\alpha-cut may exist even when both assumptions are violated. On a related note, deciding whether the Linear Complementarity problem has a solution is NP\mathsf{NP}-complete [Chu89]. The solution is unique if the problem involves a PP-matrix, but checking this condition is coNP\mathsf{coNP}-complete [Cox94]. However, using witnesses to verify whether a matrix is P-matrix or not, a total search version is shown to be in UEOPL\mathsf{UEOPL}. Our result for Alpha-HS would go in a similar vein, if the complexities of the above problems were better determined.

Another line to work could be to determine the computational complexities of other extensions of the Ham-Sandwich theorem. For other geometric problems that are total and admit unique solutions, it could be worthwhile to explore their place in the class UEOPL\mathsf{UEOPL}. Faster algorithms for computing the α\alpha-cut can also be explored.

References