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 -Ham-Sandwich theorem, for a clearer distinction. Set . Steiger and Zhao gave an algorithm that computes the dividing hyperplane in time, which is exponential in . Later, Bereg [Ber12] improved this algorithm to achieve a running time of , which is linear in but still exponential in . 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 (Total Function Nondeterministic Polynomial) of -search problems that always admit a solution. Two popular subclasses of , originally defined by Papadimitriou [Pap94], are (Polynomial Parity Argument) its sub-class . 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 is (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 was introduced by Johnson, Papadimitriou, and Yannakakis [JPY88]. A noteworthy sub-class of is (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 . In fact, this problem is complete for this class. The discrete analogue of the Borsuk-Ulam Theorem, Tucker’s Lemma [Tuc46], is also -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 . In fact, Filos-Ratsikas and Goldberg [FRG19] recently showed that Ham-Sandwich is complete for . The (presumably smaller) class is associated with fixed-point type problems: computing an approximate Brouwer fixed point is a prototypical complete problem for . The discrete analogue of Brouwer’s Fixed Point Theorem, Sperner’s Lemma, is also complete for . In a celebrated result, the relevance of for algorithmic game theory was made clear when it turned out that computing a Nash-equilibrium in a two player game is -complete [CDT09]. In discrete geometry, finding a solution to the Colorful Carathéodory problem [Bár82] was shown to lie in the intersection [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 [DP11].
Recently, at ICALP 2019, Fearley, Gordon, Mehta, and Savani defined a sub-class of that represents a family of total search problems with unique solutions [FGMS19]. They named the class Unique End of Potential Line () 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 , 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 and , such that if an instance of has no violations, then the reduced instance of 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 -Ham-Sandwich problem by locating it in . 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 , and hence in . See Figure 2 for a pictorial description.
It is not surprising to discover that Alpha-HS lies in , since the proof of the continuous version in [BHJ08] was based on Brouwer’s Fixed Point Theorem. The observation that it also lies in 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 .
Since Alpha-HS lies in , it is computationally easier than Ham-Sandwich, which is -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 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 apart from the -Matrix Linear complementarity problem and one of the few in 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 -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 -Ham-Sandwich Theorem [SZ10] here. The continuous version [BHJ08] follows a similar formulation.
We say that has very weak general position [SZ10], if for every choice of points , the affine hull of the set is a -flat and does not contain any other point of . 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 that contain multiple points from the same color. To simplify our proofs we need a slightly stronger form of general position. We say that has weak general position if the above restriction also applies to sets having exactly 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 points whose affine hull has dimension at most , with at least colors in the set. Testing whether a planar point set is in general position can be shown to be -Hard, using the result in [FKNN17]. It is easy to see that when , weak general position is equivalent to general position.
Well-separation.
A certificate for checking violations of well-separation is a colorful set whose affine hull has dimension at most . Another certificate is a partition 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 satisfying , , an -cut is an oriented hyperplane that contains one point from each color and satisfies for , where is the closed positive half-space defined by .
If an -cut exists, then it is unique.
If has weak general position, then a cut exists for each choice of , .
That means, every colorful -tuple of corresponds to exactly one -vector. Steiger and Zhao [SZ10] also presented an algorithm to compute the cut in time, where . The algorithm proceeds inductively in dimension and employs a prune-and-search technique. Bereg [Ber12] improved the pruning step to improve the runtime to .
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 be positive integers. The input consists of
The complexity class represents the class of problems that can be reduced in polynomial time to UniqueEOPL. This has been shown to lie in 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 and be two problems both having a formulation that allows for valid and violation solutions. A reduction from to is said to be promise-preserving, if whenever it is promised that has no violations, then the reduced instance of 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 -Ham-Sandwich in a non-promise setting:
A subset of of size and at least 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 to another colorful hyperplane through a sequence of local changes of the points on the hyperplanes such that the -vector of increases in some coordinate by one from that of . We next define the rotation operation in a little more detail. An anchor is a colorful -tuple of which spans a -flat. The following procedure takes as input an anchor and some point and determines the next hyperplane obtained by a rotation. The output is , where is an anchor and is some point.
Procedure
Let denote the hyperplane defined by and be the missing color in .
If the orientation of is not well-defined, report a violation of weak general position and well-separation.
Let be the subset of that lies in the closed halfspace and be the subset of that lies in the open halfspace . Let be the highest ranked point according to the order and be the highest ranked point according to .
If has color and , report out of range.
We rotate around the anchor in a direction such that the hyperplane is moving away from along the segment until it hits some point .
If the hyperplane hits multiple points at the same time, report a violation of weak general position.
If is not color , set and , where is a point in with the same color as . Otherwise, set and .
Figure 3 shows an application of this procedure, rotating to through .
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 is incident to a colorful set of points. This set of points defines possible anchors, and each anchor can be used to rotate 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 whose -vector is is given (we show later how this assumption can be removed). We start at and pick the anchor that excludes the first color, then apply a sequence of rotations until we hit another colorful hyperplane with -vector . Similarly, we move to a colorful hyperplane with -vector and so on until we reach . Then, we repeat this for the other colors in order to reach and so on until we reach the target -vector. This pattern of -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 to , with -vectors and respectively. During the rotation, we encounter a hyperplane for which its -vector is , which differs from our desired sequence of . This makes it difficult to define a potential function in the graph that strictly increases along the path where is the vertex representing hyperplane . One way to alleviate this problem is to not use as a vertex directly, but the double-wedge that is traced out by the rotation from to . If the -vector is now measured using the hyperplane that bisects the double-wedge, then we get the desired sequence of . 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 -vector is not sufficient to define the potential function, since the sequence of rotations between two colorful hyperplanes may have the same -vector. For instance, the bisectors of the rotations in in Figure 3 all have the same -vector. Hence, we need an additional measurement in order to determine the direction of rotation that increases the -vector.
Similar to how we define the orientation for a non-colorful hyperplane, let denote a hyperplane that passes through points of colors. Let denote the missing color in . Let be the highest ranked points under in and respectively. Let denote the intersection of and . We define a distance parameter called dist-value of to be the distance . In Figure 3, we can see that rotating from to sweeps the segment 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 -vector and dist-value for the tie-breaker.
Correctness.
We show that if there are no violations, we can always apply Procedure to increment the -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 , which may be swapped to form a new anchor in Step .
The well-separation condition guarantees that the potential function always increases along the rotation. Let denote a pair of hyperplanes that are the input and output of Procedure respectively. Let denote any intermediate hyperplane during the rotation from to through the common anchor. Let be the color missing from the anchor and be the highest ranked point under in . We say that the orientation of (resp. ) is consistent with that of if (resp. ). Lemma 5 shows that the orientations are always consistent when and are non-colorful hyperplanes even without the assumption of well-separation.
Assume that weak general position holds. Let be the input and output of Procedure respectively. Let denote any intermediate hyperplane within the rotation. The orientations of (resp. ) and are consistent when (resp. ) is a non-colorful hyperplane.
Since is a non-colorful hyperplane, let denote the color missing from . and give the same partition of into two sets because the continuous rotation from to does not hit any point in . Let and be the highest ranked points under in each set. Since we have weak general position, the segment cannot pass through the anchor of the rotation so that the orientations of and are well-defined by the colored points in the anchor and the intersections of the hyperplanes with the segment . Thus, the determinant defining the normal of the rotating hyperplane from to for the orientation is always non-zero. Since the intersection of the rotating hyperplane from to and the segment moves continuously along , by a continuity argument, the normal of the hyperplane does not flip during the rotation. Without loss of generality, assume that . This implies that is always in the positive half-space of and hence has a consistent orientation as . The same proof holds for . ∎
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 be a colorful hyperplane and be the first colorful hyperplane obtained by a sequence of rotations by Procedure . We denote be the non-colorful hyperplanes obtained from the above sequence of rotations. The dist-values of is strictly increasing.
Let denote the color missing from . Then, all miss the color , otherwise is not the first colorful hyperplane obtained by the rotations. Therefore, each gives the same partition of into two sets for because the continuous rotations from to does not hit any point in . Let and be the highest ranked points under in each set. Without loss of generality, assume that . Since are non-colorful hyperplanes, by Lemma 5, the consistent of the orientation can carry from to and so on. Then we have and . Let . According to Step of Procedure , each rotation is performed by moving away from along the segment . Hence we have . ∎
The last step for proving that the potential function always increases along the canonical path is to show that the -vector increases in some coordinate from one colorful hyperplane to another colorful hyperplane through Procedure . This requires the assumption of well-separation. Lemma 7 shows that if the orientations of and 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 to another colorful hyperplane through a sequence of non-colorful hyperplanes that miss color , we have and contains one additional point in that is hit by the last rotation. Therefore, is increased by and other 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 and are inconsistent, must be a colorful hyperplane by Lemma 5. Therefore, the point in that is not in the anchor is in , denoted by .
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 -vector 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 -vector 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 colors may contain additional points. This in turn implies that some -cut is missing, so that the correct solution for the target may not exist. In addition, the (highest-ranked) points from the missing color that we choose to define the orientation of a non-colorful hyperplane may form a segment that passes through the -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 -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 may be inconsistent, which may not guarantee that the -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 -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 is the number of bits used to represent each coordinate of points of .
Let denote the colorful points on and denote the colorful points on . Throughout this proof, we consider to be a closed halfspace while is an open halfspace. We prove the claim by induction on the dimension .
Without loss of generality, suppose that the open segment . Since and , there exists at least one point in order for to hold. If the open segment also lies in (resp. ), then there exists at least one point in (resp. ). We can see that the intersection point of and lies inside the triangles and (see Figure 11(b)).
Suppose that the open segment lies in (resp. ). In order to assign correct orientations to and , the order in which points of and appears on the hyperplanes along any direction must be the same for both. This is only feasible when lies between (resp. ) and the intersection point . Hence, (resp. ) lies inside the triangles (see Figure 11(c)).
The idea is to transform to a point set , 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 with the same -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 denote one of the non-colorful hyperplanes that passes through some and denote another non-colorful hyperplane that passes through some . Recall that they have the same -vector and dist-value. Let (resp. ) be the highest ranked points of under on either side of (resp. ) and let (resp. ) be the intersection of segment (resp. ) with (resp. ). By the definition and assumption of dist-value, the dist-value of and is .
In the following, we consider the case of with three sub-cases:
: the argument is symmetrical to the case above.
: we move towards along segment until it hits and move towards along segment until it hits .
Next, we consider and .
: we have . We move towards along segment until it hits and move towards along segment until it hits .
: we move towards along segment until it hits and move towards along segment until it hits .
: we have . We move towards along segment until it hits and move towards along segment until it hits .
: we have . We move towards along segment until it hits and move towards along segment until it hits .
The case for and are symmetrical to those above.
The last case is and . Basically, the sub-cases are the same as those above except one, which happens for and . In this case, we have and . We move and towards each other along segment and until they hit and . Note that the segment may intersect the -flat , but this case is also handled by Lemma 8. ∎
Alpha-HS .
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 -tuple of which spans a -flat. Let denote the missing color in the anchor. Then the tuple for the anchor is ordered as . An anchor along with a pair of points such that is called a double-wedge if all of the following hold:
the hyperplane through does not contain . This implies that the hyperplane through does not contain .
if are the highest ordered points of under on either sides of , then does not lie on a hyperplane.
the intersection of the open halfspaces is empty, that is, it does not contain any point of . Similarly must also be empty.
We visualize the anchor as a -ridge through which pass through. A rotation around the anchor changes one hyperplane to the other without passing through any other point of . Intuitively the double-wedge refers to the space and we use this interpretation several times. See Figure 5 for a simple example.
For a double-wedge , we define a representative hyperplane as the hyperplane that is the angular bisector of the double-wedge. Since a double-wedge is empty, does not contain any point of apart from . We define an orientation for based on and the color missing from . Let be points from the missing color as defined before. Without loss of generality, let , and . We call the first hyperplane among that intersects the directed segment as the upper hyperplane of and the other hyperplane as the lower hyperplane of . A simple example can be found in Figure 6.
The -vector of any oriented hyperplane is a -tuple of integers where is the number of points of in the closed halfspace for . The -vector of a double-wedge is defined as the -vector of its representative hyperplane. We say that a double-wedge is non-colorful, if both and are non-colorful, and colorful, if exactly one of and is colorful, and very colorful, if both and are colorful.
Under the assumption of weak general position, we additionally have that if is non-colorful, then are non-colorful, and if is colorful, exactly one of is colorful, and if is very colorful, both 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 is in weak general position. Later we show how to handle the cases when these assumptions fail.
We call two double-wedges neighboring if both share a common hyperplane, that is, , with an exception that we elaborate below. A double-wedge has different number of neighbors depending on how colorful its hyperplanes are. The anchor can be written in the form .
Let be non-colorful. Then both share their colors with those of . Suppose has the same color as . Then there are at most three neighboring double-wedges that share . One of them use the same anchor , and as an exception we do not count this as a neighboring double-wedge. For the two remaining neighbors, the anchor is where has replaced . The two rotational directions determine the two double-wedges. Only one of them has the same -vector as , since the representative hyperplanes contain on opposite sides. With a similar argument, there are at most two neighboring double-wedges that share and at most one of them has the same -vector as .
Let be colorful, where is colorful and is non-colorful, without loss of generality. By replacing some by we get an anchor that is contained in and which may define a double-wedge for each of the two rotational directions. Since there are possible anchors formed by replacement, there are at most double-wedges that share . Additionally, keeping the anchor fixed, there is at most one neighboring double-wedge. So there are at most neighboring double-wedges of that share . The case for is similar to case (1).
Let be very colorful. Similar to case (2), there are at most double-wedges sharing . The case for is similar.
We build a graph where each vertex corresponds to a double-wedge. Let be any double-wedge. For simplicity, we denote the vertex in corresponding to also by . If is colorful, we add an edge in between and the vertex of each neighboring double-wedge that shares . If is non-colorful, we add an edge only with the vertex of the double-wedge that shares its -vector with that of . Thus, non-colorful double-wedges have degree two in , while colorful and very colorful double-wedges have degrees at most and , respectively.
We transfer each attribute of a double-wedge to its vertex in . For instance, we call vertices of as non-colorful, colorful or very colorful corresponding to the color of the double-wedge representing the vertex. Similarly, each vertex has an -vector that corresponds to the -vector of its double-wedge, and so on. See Figure 7 for an elementary example.
Let be any vertex and let denote the -vector of . The largest -vector for any hyperplane is that occurs on a unique tangent hyperplane whose half-space contains . With our definition of the -vector of double-wedges using the representative hyperplanes, for any double-wedge , the -vector of is smaller in at least one coordinate from the maximum.
Let be a colorful hyperplane with -vector . For each , if (resp. ), then there is a path in such that the double-wedge is incident to , the double-wedge is incident to another colorful hyperplane , share the same -vector and the -vector for differs only in the -th component, where the value is (resp. ).
For the case that , we set an anchor in that excludes the -th colored point, say . Then, we apply Procedure starting from until we get another colorful hyperplane . Let be the double-wedge created by the first rotation. Note that is on the upper hyperplane of . During this sequence of rotations, we also get a sequence of double-wedges. Before the rotating hyperplane hits a point of repeated color , assume that is in the negative half-space of . Once is on , we swap with another point of the same color in the anchor and keep the rotation towards the opposite direction of the orientation so that is in the positive half-space of and from the positive half-space moves to the negative half-space. This is true because the orientation is consistent by Lemma 5. If is in the positive half-space of before is hit by , then remains in the positive half-space and as well (see Figure 3). Both cases maintain during the rotation. Thus, all non-colorful double-wedges in this sequence of rotations have the same -vector as . In the last rotation, the rotating hyperplane hits the first point of color . By Lemma 7, well-separation guarantees the consistency of the orientation of the rotating hyperplane so that moves from the negative half-space to the positive half-space and other points of color remain in the same sides of the hyperplane. Thus, is increased by one. The same argument also works for the case that by using the inverse of Procedure .
Since all the double-wedges created by the first rotation for each are incident to , they are also connected in by definition. ∎
By Lemma 12, we know that all (very) colorful double-wedges are connected. For non-colorful double-wedges, we apply Procedure 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 imitates a grid in a coarse sense. There is a ”vertex” for every colorful -tuple of , and there are paths connecting these grid vertices. We showed in Lemma 13 that is connected. Therefore, given a target -vector , the correct -tuple can be found by starting from some vertex and walking towards the solution. See Figure 8 for an illustration.
If 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 -cut, and no solutions for other cuts. The former case will manifest as multiple vertices with the same -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 that corresponds to the -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 in corresponds to in , where is a colorful point set of size from , and . We are only interested in the case when is a double-wedge, as per the definition in Section 4.1. Otherwise, we create a self loop on in . Furthermore, if there are no violations in , we can define a canonical path from the vertex with -vector to the unique vertex with -vector (shown in Section 3.1), which is the unique path in . For other vertices not on the path, we also create a self loop on . For instance, when fails weak general position assumption, then lie on the same hyperplane, where . In this case we create a cycle on with the same potential value on each , so that this violation may be reported as the violation (UV1) in . When 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 . 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.
contains bits representing the index of the missing color in ,
each contain bits for the indices of the points in ordered by ,
contain and bits for the index of the color and the index of respectively,
and we use the same idea to represent by . Altogether, we need at most bits in the encoding. Let denote the function that encodes to . If is a double-wedge, then also represents the same double-wedge. Since we do not want to create two valid vertices in corresponding to the same double-wedge, we only pick the one with on the upper hyperplane as a valid double-wedge. The following lemma details how we can verify whether a given encoding is a double-wedge:
The rotation process is handled by procedures and , which we describe next. Let be the current double-wedge. Some abnormal cases may happen in the output of or when the segment that defines the orientation of passes through (see Figure 10), or contains more than points, or the orientations of are not consistent. For these cases, the path will end or start at . When UniqueEOPL outputs , we can compute 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
Let and let be the missing color in .
Let be the subset of that lies in and be the subset of that lies in . Then, let be the highest ranked point according to the order and be the highest ranked point according to the order . As we describe previously, lies on the lower hyperplane.
Since is supposed to be on the canonical path and is not the end point, the -vector of is in the form of with .
If shares the same color of a point in , then set and .
If is in color and , then set and .
If is in color and , then let be the point in with color and set and .
We rotate around the anchor in a direction such that the hyperplane is moving away from along the segment until it hits a point .
Return .
Procedure
Let and let be the missing color in .
Let be the subset of that lies in and be the subset of that lies in . Then, let be the highest ranked point according to the order and be the highest ranked point according to the order . As we describe previously, lies on the upper hyperplane.
Since is supposed to be on the canonical path and is not the starting point, the -vector of is in the form of with . When , .
If shares the same color of a point in , then set and .
If is in color and , then set and .
If is in color and , then let be the point in with color and set and .
We rotate around the anchor in a direction such that the hyperplane is moving closer to along the segment until it hits a point .
Return .
The points and 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 , then Return , where is a colorful point set on a hyperplane that has -vector , and is the point that creates a double-wedge with the anchor and .
If is not the encoding of a double-wedge, then Return .
Let denote the double-wedge for which . If the orientations of , and are not consistent, then Return .
Recall that is the missing color in . If the -vector of is not in the form of with , then Return .
If the -vector of the lower hyperplane is , then Return .
If contains some points of other than , then let be those extra points ordered by , and let be the point just after in the order (if is after , then ). We Return .
Similarly, Return if contains some points other than , where is the point just after according to the order among those extra points.
Let be the output of . If the orientations of , and are well-defined and consistent, then Return .
If , then Return .
If is not a double-wedge, then Return .
Let be a double-wedge such that . If the orientations of , and are not consistent, then Return .
If and , then Return .
If is colorful, and the -vector of is , then Return .
Recall that is the missing color in . If the -vector of is not in the form of with , then Return .
If contains some points of other than , then let be those extra points ordered by , and let be the point just before in (if is before , then ), and Return .
Similarly, Return if contains some points other than , where is the point just before in among those extra points.
Let be the output of . If the orientations of , and are well-defined and consistent, then Return .
Given a double-wedge , let be any point of in and let be any point of in . Suppose that the segment does not pass through the affine hull of so that the intersections of with and are two distinct points. Let be the Euclidean distance between the two intersection points of with and . Let be the Euclidean distance between and the intersection point of with .
Let denote the number of bits needed to represent each coordinate of any point of . Then, is at less and is at most , where and .
Without loss of generality, we assume that the missing color in is color 1, i.e., is a set of points . Let (resp. ) be the intersection of and (resp. ). Since lies on and the affine hull of , we can represent as the convex combination of and and the linear combination of as follows:
where are column vectors. Since we assume that exists, we have . According to Cramer’s rule, we have , where is the matrix obtained by replacing the -th column of with . Using Leibniz formula for determinants, we can bound the denominator:
For , it is less than the longest possible line segment, so . ∎
If is not a double-wedge, then Return . Otherwise, let be a double-wedge such that .
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 -cut from the unique solution of the UniqueEOPL instance.
If there are no violations in , then the constructed UniqueEOPL instance only contains a type (U1) solution whose lower hyperplane is colorful and has -vector , which is a type (G1) solution of .
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 in the open halfspace ,
The proof for points in the open halfspace is similar. ∎
By Theorem 2, we already know that if there are multiple -cuts, then 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 and , for which no consistent orientation can be defined. This hyperplane gives the certificate of violation.
Let denote the colorful points on and denote the colorful points on .
The open segment is contained in either or for each .
We now extend Lemma 8 to the case where the two hyperplanes are non-colorful but have the same -vector and dist-value. On any path of non-colorful vertices of 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 -vector , where is the missing color for those non-colorful vertices and . If we find another non-colorful vertex with the same missing color and -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 and is color . Let be the colorful points on and be the colorful points on . Let (resp. ) be the highest ranked points of under on either side of (resp. ) 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 and . We can apply the same analysis from Lemma 9 for and in other cases.
When and , we have and . Since is the representative hyperplane, is the lower hyperplane and the dist-value of is larger than that of , the directed segment intersects them in the following order: , then . Since is a double-wedge, the orientations of and are consistent. That means, . There are two possibilities:
: using a similar argument as in the proof of Lemma 9, there must exist another point because and have the same -vector. Then, we move towards along segment until it hits and move towards along segment until it hits (Figure 12(a)). We apply Lemma 8 to the resulting point set and that completes the proof.
: in this case both and lie in , but since that would make the intersection order along as , which is a contradiction. Along the directed segment , the order of intersection is (Figure 12(b)) while on the directed segment it is . So we move towards along until it hits , and move towards along until it hits . Then we re-use the proof from Lemma 8.
In the following lemmas, we show how to reduce any violation solution of to a violation solution of 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 be the output of . The orientations of 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 . 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 . 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 -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 -cut exists for the point set. A given -cut may exist even when both assumptions are violated. On a related note, deciding whether the Linear Complementarity problem has a solution is -complete [Chu89]. The solution is unique if the problem involves a -matrix, but checking this condition is -complete [Cox94]. However, using witnesses to verify whether a matrix is P-matrix or not, a total search version is shown to be in . 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 . Faster algorithms for computing the -cut can also be explored.