Pizza Sharing is PPA-hard

Argyrios Deligkas, John Fearnley, Themistoklis Melissourgos

Introduction

Mass partition problems ask to fairly divide measurable objects that are embedded into Euclidean space [RS20]. Perhaps the most popular mass partition problem is the ham sandwich problem, in which three masses are given in three-dimensional Euclidean space, and the goal is to find a single plane that cuts all three masses in half. Recently, there has been interest in pizza sharing problems, which are mass partition problems in the two-dimensional plane, and in this paper we study the computational complexity of such problems.

In the straight-cut pizza sharing problem, we are given 2n2n two-dimensional masses in the plane, and we are asked to find nn straight lines that simultaneously bisect all of the masses. See Figure 1(a) for an example. It has been shown that this problem always has a solution: the first result on the topic showed that solutions always exist when n=2n=2 [BPS19], and this was subsequently extended to show existence for all nn [HK20].

R^{+} and R−R^{-} (shaded and non-shaded areas respectively). Another related problem is the square-cut pizza sharing. In this problem, there are nn masses in the plane, and the task is to simultaneously bisect all masses using cuts, but the method of generating the cuts is different. Specifically, we seek a square-cut, which consists of a single path that is the union of horizontal and vertical line segments. See Figure 1(b) and Figure 1(c) for two examples of square-cuts. Intuitively, we can imagine that a pizza cutter is placed on the plane, and is then moved horizontally and vertically without being lifted in order to produce the cut. Note that the path is allowed to wrap around on the horizontal axis: if it exits the left or right boundary, then it re-appears on the opposite boundary. So the cut in Figure 1(c) is still considered to be a single square-cut.

It has been shown by [KRPS16] that, given nn masses, there always exists a square-cut-path (termed SC-path) which makes at most n−1n-1 turns and simultaneously bisects all of the masses. This holds even if the SC-path is required to be yy-monotone, meaning that the path never moves downwards.

Two-dimensional fair division is usually called land division in the literature. Land division is a prominent topic of interest in the Economics and AI communities that studies ways of fairly allocating two-dimensional objects among nn agents [Cha05, SHNHA17, SNHA20, ESS21, AD15, IH09, Hüs11]. The first popular appearance of such problems in a mathematical description was done by [Ste48], and since then, the existence of allocations under various fairness criteria have been extensively studied, together with algorithms that achieve them. These problems find applications from division of resources on land itself, to the Law of the Sea [SS03], to redistricting [LRY09, LS14].

Consensus halving is a problem that asks us to split a one-dimensional resource into two parts such that nn agents have equal value in both parts. Here, we study the same fairness criterion for nn agents, but for a two-dimensional resource. One can see that when we have the same fairness criterion at hand for any kk-dimensional resource, k≥2k\geq 2, we can always translate the problem into its one-dimensional version, by integrating each agent’s measure to a single dimension. Then a solution can be given by applying consensus halving. However, the solutions we get by doing so, are not taking into account the dimensionality of the problem, and as a result they might produce very unnatural solutions to a high-dimensional problem. For example, in land division, applying consensus halving would produce two parts, each of which can possibly be a union of ⌈n/2⌉\left\lceil n/2\right\rceil disjoint land strips. Can we get better solutions by exploiting all the dimensions of the problem?

In this work we investigate different cutting methods of the two-dimensional objects, and in particular, two pizza sharing methods for which a solution is guaranteed. While based on intuition one might assume that exploiting the two dimensions would allow the complexity of finding a solution to be lower, our results show that this is not the case. We present polynomial time reductions from the one-dimensional problem to the two-dimensional problems showing that the latter are at least as hard as the former, i.e., PPA-hard. Apart from the hardness results themselves, we believe that our reductions are interesting from another aspect too. They show ways to efficiently turn a problem into one of higher dimension, a task that has no standardised methods to be achieved (even for non-efficient reductions), and whose inverse is trivially achievable.

Computational complexity of fair division problems. There has been much interest recently in the computational complexity of fair division problems. In particular, the complexity class PPA has risen to prominence, because it appears to naturally capture the complexity of solving these problems. For example, it has recently been shown by [FRG18, FRG19] that the consensus halving problem, the ham sandwich problem, and the well-known necklace splitting problem are all PPA-complete.

The other class of relevance here is the class BU, which consists of all problems that can be polynomial-time reduced to finding an exact solution to a Borsuk-Ulam function. This class was defined by [DFMS21] and is believed to be substantially harder than the class PPA, because it is possible to construct a Borsuk-Ulam function that only has irrational solutions. Due to this, it is not currently expected that BU will be contained in FNP, whereas the containment of PPA in FNP is immediate.

Unfortunately, it is not currently known whether BU has complete problems. However, in [DFMS21] it was shown that exact consensus halving is in BU and also FIXP-hard, implying that FIXP⊆BU\textup{{FIXP}}\subseteq\textup{{BU}}. FIXP, defined by Etessami and Yannakakis [EY10], is the class of problems that can be reduced to finding an exact fixed point of a Brouwer function. It is known by the aforementioned work, that FIXP contains the problem Square Root Sum, which has as input positive integers a1,…,ana_{1},\dots,a_{n} and kk, and asks whether ∑i=1nai≤k\sum_{i=1}^{n}\sqrt{a_{i}}\leq k. The question of whether Square Root Sum is in NP has been open for more than 40 years ([GGJ76, Pap77, Tiw92]). Furthermore, since there exist Brouwer functions that only have irrational fixed points, it is likewise not expected that FIXP will be contained in FNP.

Our contribution. We study the computational complexity of the straight-cut and square-cut pizza sharing problems, and we specifically study the cases where (i) all mass distributions are unions of weighted polygons, and (ii) we are given unweighted point sets. We show that it is PPA-complete to find approximate solutions of the approximate versions of all the problems, while their decision variants are NP-complete. We also show that the exact square-cut pizza sharing problem with mass distributions is FIXP-hard and in BU, while its decision variant is ETR-complete. All of our hardness results are summarized in Table 1 and Table 2.

These results represent, to the best of our knowledge, the first PPA-hardness results for problems arising from computational geometry. We also note that pizza sharing problems do not need a circuit as part of the input, which makes them in some sense more “natural” than problems that are specified by circuits. Other known “natural” PPA-hard problems are one-dimensional, such as consensus halving [FHSZ20] and necklace splitting [FRG19]. Here we show the first known PPA-hardness result for a “natural” two-dimensional problem. Let us also mention here that shortly after the appearance of our result, Schnider in [Sch21] proved that the discrete version of straight-cut pizza sharing where each mass is represented by unweighted points is PPA-complete, while its continuous version for a more general input representation is FIXP-hard.

Recently, [BHH21] made great improvement towards showing BU-hardness of exact consensus halving. They introduced a class named BBU whose typical problem is an equivalent, alternative definition of the Borsuk-Ulam theorem, and they considered the strong approximation version of the classes BU and BBU, named BUa\textup{{BU}}_{a} and BBUa\textup{{BBU}}_{a}, respectively. Some of the most notable results of the aforementioned paper is that BUa=BBUa\textup{{BU}}_{a}=\textup{{BBU}}_{a} and that the strong approximation version of consensus halving is complete for BUa\textup{{BU}}_{a}. We believe that some of our reductions will be able to be translated into the framework of strong approximation and yield analogous BUa\textup{{BU}}_{a}-hardness results for the strong approximation version of pizza sharing problems. However, we remark that BU-completeness of either exact consensus halving or any of the pizza sharing problems is yet to be proven.

For both the straight-cut and the square-cut pizza sharing problems, we show that it is PPA-complete to find an ε\varepsilon-approximate solution for any constant ε∈(0,1/5)\varepsilon\in(0,1/5). This holds even when n+n1−δn+n^{1-\delta} lines are permitted in a straight-cut pizza sharing instance with 2n2n mass distributions, and when n−1+n1−δn-1+n^{1-\delta} turns of the square-cut path are permitted in a square-cut pizza sharing instance with nn mass distributions, for constant δ>0\delta>0. Furthermore, the PPA-hardness holds even when each mass distribution is uniform over polynomially many axis-aligned rectangles, and there is no overlap between any two mass distributions. The inapproximability for such high values of ε\varepsilon is possible due to a recent advancement in the inapproximability of consensus halving [DFHM22]. The PPA-inclusion holds even for inverse polynomial and inverse exponential ε\varepsilon, and for weighted polygons with holes (that is, the most general type of allowed input). Furthermore, we show that there exists an ε>0\varepsilon>0 such that deciding whether an ε\varepsilon-approximate solution of straight-cut pizza sharing with at most n−1n-1 lines (resp. an ε\varepsilon-approximate solution of square-cut pizza sharing with at most n−2n-2 turns) exists or, is NP-complete. All of these results hold also for the discrete version of the problems.

We then turn our attention to the computational complexity of finding an exact solution to the square-cut problem. We show that the problem of finding an SC-path with at most n−1n-1 turns that exactly bisects nn masses lies in BU, and is FIXP-hard. This hardness result applies even if all mass distributions are unions of weighted axis-aligned squares and right-angled triangles. In order to prove containment in BU, we provide a simpler existence proof for a solution to the square-cut pizza sharing problem that follows the lines of the original proof by [KRPS16], while to prove the FIXP-hardness, we reduce from the problem of finding an exact Consensus-Halving solution [DFMS21]. Regarding the decision version of the square-cut problem, we show that deciding whether there exists an exact solution with at most n−2n-2 turns is ETR-complete, where ETR consists of every decision problem that can be formulated in the existential theory of the reals (see Section 2 for its definition).

From a technical viewpoint, our PPA containment result for straight-cut pizza sharing is based on a reduction that transforms mass distributions to point sets in general position and then employs a recent result by [Sch21]. On the contrary, our containment results for square-cut pizza sharing are shown by directly reducing it to the Borsuk-Ulam problem. Our hardness results are obtained by reducing from the consensus halving problem, historically the first fair-division problem shown to be PPA-complete [FRG18]. It is worth mentioning that, if in the future, consensus halving with linear valuation densities (see definitions of kk-block-triangle valuations in Section 2) is shown to be BU-complete, then our work implies BU-completeness of exact square-cut pizza sharing.

Further related work. Since mass partitions lie in the intersection of topology, discrete geometry, and computer science there are several surveys on the topic; [BFHZ18, DLGMM19, Mat08, Živ17] focus on the topological point of view, while [AE+99, Ede12, KK03, Mat02] focus on computational aspects. Consensus halving [SS03] is the mass partition problem that received the majority of attention in Economics and Computation so far [DFH21, FRFGZ18, FRG19, FHSZ20, FRHSZ21]. Recently, Haviv [Hav22] showed PPA-completeness of finding fair independent sets on cycle graphs, having as a starting point the latter problem.

Preliminaries

Mass distributions. A mass distribution μ\mu on 2^{2} is a measure on the plane such that all open subsets of 2^{2} are measurable, 0<μ(2)<∞0<\mu\left(^{2}\right)<\infty, and μ(S)=0\mu(S)=0 for every subset of 2^{2} with dimension lower than 2. A mass distribution μ\mu is finite-separable, or simply separable, if it can be decomposed into a finite set of non-overlapping areas a1,a2,…,ada_{1},a_{2},\ldots,a_{d} such that μ(2)=∑j=1dμ(aj)\mu\left(^{2}\right)=\sum_{j=1}^{d}\mu(a_{j}). In addition, a separable mass distribution μ\mu is piece-wise uniform, if for every jj and every S⊆ajS\subseteq a_{j} it holds that μ(aj∩S)=wj⋅area(aj∩S)\mu(a_{j}\cap S)=w_{j}\cdot\text{area}(a_{j}\cap S) for some weight wj>0w_{j}>0 independent of SS. When additionally wj=wkw_{j}=w_{k} for all j,k∈[d]j,k\in[d] then the mass distribution is called uniform. Finally, a mass distribution is normalised if μ(2)=1\mu(^{2})=1. The support of mass distribution ii, denoted by supp(i)supp(i), is the area Ai⊆2A_{i}\subseteq^{2} which has the property that for every S⊆AiS\subseteq A_{i} with non-zero Lebesgue measure we have μi(S)>0\mu_{i}(S)>0. Let N:={I⊆[n]:⋂i∈Isupp(i)≠∅}N:=\{I\subseteq[n]:\bigcap_{i\in I}supp(i)\neq\emptyset\}. A set of mass distributions μ1,…,μn\mu_{1},\dots,\mu_{n}, or colours, has overlap kk if max⁡I∈N∣I∣=k\max_{I\in N}|I|=k. In order to easily present our results that concern additive approximations, throughout this work we consider all mass distributions to be normalised, which is without loss of generality.

Set of straight-cuts. A set of straight-cuts, or cut-lines, or simply lines defines subdivisions of the plane RR. Figure 1(a) shows an example of a set of straight-cuts. Each line creates two half-spaces, and arbitrarily assigns number “” to one and “11” to the other. A subdivision of RR is labeled “++” (and belongs to R+R^{+}) if its parity is odd (according to the labels given to the half-spaces) and “−-” (and belongs to R−R^{-}) otherwise. Observe that by flipping the numbers of two half-spaces defined by a line, we flip all the subdivisions’ labels. Thus, there are only two possible labelings of the subdivisions.

Square-cut-path. A square-cut-path, denoted for brevity SC-path, is a non-crossing directed path that is formed only by horizontal and vertical line segments and in addition it is allowed to “wrap around” in the horizontal dimension. Figure 1(b), Figure 1(c) show two examples of SC-paths. A turn of the path is where a horizontal segment meets with a vertical segment. An SC-path is yy-monotone if all of its horizontal segments are monotone with respect to the yy axis. Any SC-path partitions the plane into two regions, that we call R+R^{+} and R−R^{-}.

Pizza sharing. A set of lines (resp. a SC-path) ε\varepsilon-bisects a mass distribution μ\mu, if ∣μ(R+)−μ(R−)∣≤ε|\mu(R^{+})-\mu(R^{-})|\leq\varepsilon. It simultaneously ε\varepsilon-bisects a set of mass distributions MM if ∣μi(R+)−μi(R−)∣≤ε|\mu_{i}(R^{+})-\mu_{i}(R^{-})|\leq\varepsilon for every μi∈M\mu_{i}\in M.

R^{+} and R−R^{-} using at most nn lines such that for each i∈[2n]i\in[2n], it holds that ∣μi(R+)−μi(R−)∣≤ε|\mu_{i}(R^{+})-\mu_{i}(R^{-})|\leq\varepsilon. Definition 3. For any n≥1n\geq 1, the problem ε\varepsilon-SC-Pizza-Sharing is defined as follows: • Input: ε≥0\varepsilon\geq 0 and mass distributions μ1,μ2,…,μn\mu_{1},\mu_{2},\ldots,\mu_{n} on 2^{2}. • Output: A partition of 2^{2} to R+R^{+} and R−R^{-} using a SC-path with at most n−1n-1 turns such that for each i∈[n]i\in[n], it holds that ∣μi(R+)−μi(R−)∣≤ε|\mu_{i}(R^{+})-\mu_{i}(R^{-})|\leq\varepsilon. In [KRPS16], it was proved that ε\varepsilon-SC-Pizza-Sharing always admits a solution for arbitrary continuous measures with respect to the Lebesgue measure, and for any ε≥0\varepsilon\geq 0. In this work, we are interested in the computational aspect of the problem, hence we need to specify its input representation. For simplicity, we deal with measures determined by sets of polygons with holes, since even these simply-defined measures suffice to yield PPA-hardness. The precise input representation is described in Appendix A.

Although it is not hard to construct instances of ε\varepsilon-SC-Pizza-Sharing where n−1n-1 turns are necessary for any SC-path in order to constitute a solution, there might be cases where a solution can be achieved with an SC-path with fewer turns. Hence, we also study the decision version of the problem, in which we ask whether we can find a solution with kk turns, where k<n−1k<n-1.

Consensus halving. The hardness results that we will show in this paper will be shown by a reduction from the consensus halving problem.

In the ε\varepsilon-Consensus-Halving problem, there is a set of nn agents with valuation functions viv_{i} over the interval $,andthegoalistofindapartitionoftheintervalintosubintervalslabelledeither“, and the goal is to find a partition of the interval into subintervals labelled either “+”or“” or “-”,usingatmost”, using at mostncuts.Thispartitionshouldsatisfythatforeveryagentcuts. This partition should satisfy that for every agenti,thetotalvaluefortheunionofsubintervals, the total value for the union of subintervals\mathcal{I}^{+}labelled“labelled “+”andthetotalvaluefortheunionofsubintervals” and the total value for the union of subintervals\mathcal{I}^{-}labelled“labelled “-”isthesameupto” is the same up to\varepsilon,i.e.,, i.e.,|v_{i}(\mathcal{I}^{+})-v_{i}(\mathcal{I}^{-})|\leq\varepsilon.Wewillconsiderthefollowingtypesforavaluationfunction. We will consider the following types for a valuation functionv_{i}$; see Figure 2 for a visualization.

kk-block uniform. viv_{i} is kk-block and the density of every interval is cic_{i}.

Complexity classes. ε\varepsilon-SC-Pizza-Sharing is an example of a total problem, which is a problem that always has a solution. The complexity class TFNP (Total Function NP) defined in [MP91], contains all total problems whose solutions can be verified in polynomial time.

There are several well-known subclasses of TFNP that we will use in this paper. The class PPA, defined by Papadimitriou [Pap94], captures problems whose totality is guaranteed by the parity argument on undirected graphs. The complexity class PPAD⊆PPA\textup{{PPAD}}\subseteq\textup{{PPA}} is the subclass of PPA containing all problems whose totality is guaranteed by the parity argument on directed graphs.

The complexity class ETR consists of all decision problems that can be formulated in the existential theory of the reals [Mat14, Sch09]. It is known that NP⊆ETR⊆FNP\textup{{NP}}\subseteq\textup{{ETR}}\subseteq\textup{{FNP}} [Can88], and it is generally believed that ETR is distinct from the other two classes. The class FETR (Function ETR) consists of all search problems whose decision version is in ETR. The class TFETR is the subclass of FETR which contains only problems that admit a solution (i.e. all the instances of their decision version are “yes” instances). Both FETR and TFETR were introduced in [DFMS21] as the natural analogues of FNP and TFNP in the real RAM model of computation. For a definition of the real RAM model we refer the reader to the detailed work of Erickson, van der Hoog, and Miltzow [EvM20].

In this paper, our focus regarding computational complexity will be on the following two classes contained in TFETR. As mentioned earlier, the class BU⊆TFETR\textup{{BU}}\subseteq\textup{{TFETR}} was introduced in [DFMS21] and captures problems whose totality is guaranteed by the Borsuk-Ulam theorem. The class FIXP⊆BU\textup{{FIXP}}\subseteq\textup{{BU}} was defined in [EY10] and captures problems whose totality is guaranteed by Brouwer’s fixed point theorem.

Hardness results

Here we show all hardness results regarding the exact and approximate versions of our pizza sharing problems for mass distributions, as well as for point sets.

We start by proving that ε\varepsilon-Straight-Pizza-Sharing is PPA-hard for any ε<1/5\varepsilon<1/5, even for very simple mass distributions. We prove our result via a reduction from ε\varepsilon-Consensus-Halving with kk-block valuations, which for the special case of 33-block uniform valuations has been shown to be PPA-complete [DFHM22]. In addition, we explain how to combine the machinery of our reduction with that of [FRFGZ18] in order to get NP-hardness for ε\varepsilon-Straight-Pizza-Sharing, where ε>0\varepsilon>0 is a small constant.

We reduce from Consensus-Halving with 2n2n agents, and for each agent we create a corresponding mass in Straight-Pizza-Sharing. Firstly, we finely discretize the $intervalintoblocksandweplacetheblocksoninterval into blocks and we place the blocks ony=x^{2},where, wherex\geq 0.So,the. So, theintervalcorrespondstoapartofthequadraticequation.Thisguaranteesthateverylinecancutthis“bent”intervalatmosttwiceandinadditionthepartofeachmassthatisininterval corresponds to a part of the quadratic equation. This guarantees that every line can cut this “bent” interval at most twice and in addition the part of each mass that is inR^{+}isalmostthesameasvalueofthecorrespondingagentforthepieceofis almost the same as value of the corresponding agent for the piece oflabelledwith“labelled with “+”.Nextweshowhowtoconstructaninstance”. Next we show how to construct an instanceI_{P}ofof(\varepsilon-\varepsilon^{\prime})−Straight−Pizza−Sharingwith-Straight-Pizza-Sharing with2nmassdistributions,foranymass distributions, for any\varepsilon^{\prime}=\frac{1}{\textup{poly}(n)}<\varepsilon,givenaninstance, given an instanceI_{\textsc{CH}}ofof\varepsilon−Consensus−Halvingwith-Consensus-Halving with2nagentswithagents withk$-block valuations.

Let cmax⁡:=max⁡i∈[2n],m∈[k]cimc_{\max}:=\max_{i\in[2n],m\in[k]}c_{im}, where cimc_{im} is the value density of agent ii’s mm-th block in I\textscCHI_{\textsc{CH}}, and observe that cmax⁡≥1c_{\max}\geq 1 since the total valuation of any agent over $isis1.Inwhatfollows,itwillhelpustothinkoftheinterval. In what follows, it will help us to think of the intervalininI_{\textsc{CH}}asbeingdiscretizedinincrementsofas being discretized in increments ofd:=\frac{1}{\left\lceil 4\cdot n^{2}\cdot c_{\max}\right\rceil}.Werefertothesubinterval. We refer to the subinterval\left[(j-1)\cdot d,j\cdot d\right]astheas thej−th-thd−blockofinterval-block of intervalininI_{\textsc{CH}},for, forj\in[1/d]$.

We now describe the instance IPI_{P}. For ease of presentation, the space of the instance is inflated to [0,(6d2)2+1]2\left[0,\left(\frac{6}{d^{2}}\right)^{2}+1\right]^{2}, but by scaling the construction down, the results are attained. We consider two kinds of square tiles; 1/d1/d large square tiles of size 1×11\times 1, each of which contains 2n2n smaller square tiles of size 12n×12n\frac{1}{2n}\times\frac{1}{2n} on its diagonal. We will call the former type big-tile and denote it by tjt_{j} and the latter one small-tile and denote it by tijt_{ij} for some i∈[2n]i\in[2n], j∈[1/d]j\in[1/d].

The centers of consecutive big-tiles are positioned with 6d\frac{6}{d} distance apart in the xx-axis. For j∈[1/d]j\in[1/d] we define sj=6d⋅js_{j}=\frac{6}{d}\cdot j. For every agent i∈[2n]i\in[2n] we will create a uniform mass distribution μi\mu_{i} that consists of at most 1/d1/d many axis-aligned small-tiles. Each big-tile tjt_{j}, j∈[1/d]j\in[1/d], is centered at (sj,sj2)\left(s_{j},s_{j}^{2}\right) while, in it, each small-tile tijt_{ij}, i∈[2n]i\in[2n], belonging to mass distribution μi\mu_{i} has its bottom left corner at (sj+i−12n,sj2+i−12n)\left(s_{j}+\frac{i-1}{2n},s_{j}^{2}+\frac{i-1}{2n}\right). Each small-tile tijt_{ij} contains total mass (belonging to μi\mu_{i}) of exactly the same total value vijv_{ij} as that of agent ii’s jj-th dd-block in I\textscCHI_{\textsc{CH}}. Observe that, by definition, the value vijv_{ij} of the jj-th dd-block for agent ii is at most d⋅cmax⁡≤12n⋅2nd\cdot c_{\max}\leq\frac{1}{2n\cdot 2n}, and therefore it fits inside the small tile of size 12n×12n\frac{1}{2n}\times\frac{1}{2n}. In particular, the mass inside tijt_{ij} has width 12n\frac{1}{2n} and height vij⋅2nv_{ij}\cdot 2n. Figure 3 and Figure 4 depict our construction.

The proof of the following lemma can be found in LABEL:app:_{lem:lines-to-cuts}.

This, together with the fact that ε\varepsilon-Consensus-Halving is PPA -hard for any ε<1/5\varepsilon<1/5, due to [DFHM22], implies the main theorem of this section.

ε\varepsilon-Straight-Pizza-Sharing with 2n2n mass distributions is PPA-hard for any constant ε<1/5\varepsilon<1/5, even when n+n1−δn+n^{1-\delta} lines are allowed for any given constant δ>0\delta>0, every mass distribution is uniform over polynomially many rectangles, and there is no overlap between any two mass distributions.

We will now shift our attention to studying the decision version of the problem, where we are asking to find a solution that uses at most n−1n-1 straight lines, and notice that there is no guarantee for such a solution. We employ the NP-hard instances of ε\varepsilon-Consensus-Halving for their constant ε\varepsilon from [FRFGZ18], and reduce them according to the above reduction procedure to (ε−ε′)(\varepsilon-\varepsilon^{\prime})-Straight-Pizza-Sharing instances for some ε′\varepsilon^{\prime} that is inverse polynomial in the input size, e.g., ε′=1/n2\varepsilon^{\prime}=1/n^{2}. This gives the following.

There exists a constant ε>0\varepsilon>0 for which it is NP-hard to decide if an ε\varepsilon-Straight-Pizza-Sharing instance with 2n2n mass distributions admits a solution with at most n−1n-1 lines, even when every mass distribution is uniform over polynomially many rectangles, and there is no overlap between any two mass distributions.

2 Hardness of approximate SC-Pizza-Sharing

In this section we prove hardness results for ε\varepsilon-SC-Pizza-Sharing. We provide a reduction from ε\varepsilon-Consensus-Halving with kk-block valuations, which was shown to be PPA-complete in [DFHM22] for any constant ε<1/5\varepsilon<1/5 even for 33-block uniform valuations. Our construction shows that ε\varepsilon-SC-Pizza-Sharing remains PPA-hard even when there is no overlap between any two mass distributions. Notice that the case of non-overlapping mass distributions is the most simple type of an instance, since we can easily reduce it to one where an arbitrarily large number of masses overlap: make “dummy” exact copies of an existing mass distribution. Also, the machinery that we present, combined with the reduction by [FRFGZ18] implies NP-hardness for the decision version of ε\varepsilon-SC-Pizza-Sharing when ε\varepsilon is a small constant.

We reduce from a general ε\varepsilon-Consensus-Halving instance to an ε\varepsilon-SC-Pizza-Sharing instance, and the idea is to create a mass for each agent. For any ε′=1poly(n)<ε\varepsilon^{\prime}=\frac{1}{\textup{poly}(n)}<\varepsilon, given an instance I\textscCHI_{\textsc{CH}} of ε\varepsilon-Consensus-Halving with nn agents and kk-block valuations we will show a polynomial time construction to an (ε−ε′)(\varepsilon-\varepsilon^{\prime})-SC-Pizza-Sharing instance I\textscSCI_{\textsc{SC}}.

For our construction, we will use the same components as those in the proof of Lemma 4. In particular, let cmax⁡:=max⁡i∈[n],m∈[k]cimc_{\max}:=\max_{i\in[n],m\in[k]}c_{im}, where cimc_{im} is the value density of agent ii’s mm-th block in I\textscCHI_{\textsc{CH}} (and again note that cmax⁡≥1c_{\max}\geq 1 since the total valuation of the agent over $isis1).Similarlytotheaforementionedproof,wewilldiscretizethe). Similarly to the aforementioned proof, we will discretize theintervalofinterval ofI_{\textsc{CH}}inincrementsofin increments ofd:=\frac{1}{\left\lceil n^{2}\cdot c_{\max}\right\rceil}.Also,letusrestatethatforanygiven. Also, let us restate that for any givenj\in[1/d],thesubinterval, the subinterval\left[(j-1)\cdot d,j\cdot d\right]iscalledis calledj−th-thd−blockofinterval-block of intervalininI_{\textsc{CH}}$.

To simplify the presentation, all the analysis of instance I\textscSCI_{\textsc{SC}} will be for the scaled version of it, that is, in [0,1d]2\left[0,\frac{1}{d}\right]^{2}. It is easy to verify, though, that by proper re-scaling we can put the instance in 2^{2} as required. We will be using the same gadgets that were constructed for the proof of Lemma 4, namely the big-tiles, which contain small-tiles. In particular, we have 1/d1/d square big-tiles of size 1×11\times 1, each of which contains nn square small-tiles of size 1n×1n\frac{1}{n}\times\frac{1}{n} on its diagonal. For any i∈[n]i\in[n], j∈[1/d]j\in[1/d], we denote them by tjt_{j} and tijt_{ij}, respectively.

In this construction, however, the positioning of big-tiles will be different than that of the aforementioned proof. In particular, we will be placing them on the diagonal of [0,1d]2\left[0,\frac{1}{d}\right]^{2} as shown in Figure 4(b). For every agent ii we will create a uniform mass distribution μi\mu_{i} that consists of at most 1/d1/d many axis-aligned small-tiles. For each j∈[1/d]j\in[1/d], the bottom-left corner of big-tile tjt_{j} is at (j−1,j−1)\left(j-1,j-1\right). In it, each small-tile tijt_{ij}, i∈[n]i\in[n], belonging to mass distribution μi\mu_{i} has its bottom left corner at (j−1+i−1n,j−1+i−1n)\left(j-1+\frac{i-1}{n},j-1+\frac{i-1}{n}\right).

Inside a given big-tile tjt_{j}, the total mass of each small-tile tijt_{ij} contains total mass (belonging to μi\mu_{i}) of exactly the same total value vijv_{ij} as that of agent ii’s jj-th dd-block in I\textscCHI_{\textsc{CH}}. The shape of that mass is rectangular, and has width 1/n1/n and height vij⋅nv_{ij}\cdot n. Using identical arguments to those of Lemma 4 we conclude that the total value of the dd-block fits inside the small tile of size 1n×1n\frac{1}{n}\times\frac{1}{n}. Figure 4 depicts our construction.

We will now define how a solution to I\textscSCI_{\textsc{SC}}, i.e., an SC-path, is mapped back to a solution of I\textscCHI_{\textsc{CH}}, i.e., a set of cuts. This is identical to that of the aforementioned lemma’s translation. In particular, we consider again the big-tiles in sequential order and we add one cut at j⋅dj\cdot d whenever we find two big-tiles tjt_{j} and tj+2t_{j+2} that belong to different regions. Suppose that, following the aforementioned procedure, the next I\textscCHI_{\textsc{CH}} cut falls at j⋅d′j\cdot d^{\prime} for some d′>dd^{\prime}>d. If tjt_{j} belongs to region “++” (resp. “−-”) and tj+2t_{j+2} belongs to “−-” (resp. “++”), then the interval [j⋅d,j⋅d′][j\cdot d,j\cdot d^{\prime}] gets label “++” (resp. “−-”), and vice versa. This translation obviously takes polynomial time.

Fix a constant δ>0\delta>0, and some ε′∈[2nr,ε)\varepsilon^{\prime}\in[\frac{2}{n^{r}},\varepsilon), where r≥1r\geq 1 and 2nr<ε<1\frac{2}{n^{r}}<\varepsilon<1. Let an SC-path with at most n−1+n1−δn-1+n^{1-\delta} turns be a solution to (ε−ε′)(\varepsilon-\varepsilon^{\prime})-SC-Pizza-Sharing instance I\textscSCI_{\textsc{SC}}. Then we can find in polynomial time a solution to ε\varepsilon-Consensus-Halving instance I\textscCHI_{\textsc{CH}} with at most n+n1−δn+n^{1-\delta} cuts.

The above translation of the SC-path to I\textscCHI_{\textsc{CH}} cuts indicates that each cut introduces a discrepancy between the “++” and “−-” regions of I\textscCHI_{\textsc{CH}}, of value at most cmax⁡⋅dc_{\max}\cdot d for each valuation viv_{i}, i∈[n]i\in[n]. This results in total discrepancy of (n+n1−δ)⋅cmax⁡⋅d(n+n^{1-\delta})\cdot c_{\max}\cdot d for each agent ii. We now denote by I+\mathcal{I}^{+} and I−\mathcal{I}^{-} the regions in I\textscCHI_{\textsc{CH}} which correspond to the regions R+R^{+} and R−R^{-} of I\textscSCI_{\textsc{SC}} respectively, induced by the SC-path after disregarding the intersected big-tiles. Then for the discrepancy in the valuation of agent ii in I\textscCHI_{\textsc{CH}} we get

where the last inequality is due to the fact that in our reduction we can pick d≤ε′2n⋅cmax⁡d\leq\frac{\varepsilon^{\prime}}{2n\cdot c_{\max}}. This is always possible since, by assumption, ε′≥2/nr\varepsilon^{\prime}\geq 2/n^{r} for some r≥1r\geq 1, and additionally, according to our reduction, dd is required to be at most 1n2⋅cmax⁡\frac{1}{n^{2}\cdot c_{\max}}. ∎

The above, together with the fact that ε\varepsilon-Consensus-Halving is PPA -hard for any ε<1/5\varepsilon<1/5 ([DFHM22]) implies the main theorem of this section.

ε\varepsilon-SC-Pizza-Sharing with nn mass distributions is PPA-hard for any constant ε<1/5\varepsilon<1/5, even when n−1+n1−δn-1+n^{1-\delta} turns are allowed in the SC-path for any given constant δ>0\delta>0, every mass distribution is uniform over polynomially many rectangles, and there is no overlap between any two mass distributions.

Similarly to the case of the decision variant of ε\varepsilon-Straight-Pizza-Sharing (Theorem 6), we can get the following result by reducing from the NP-hard instances of ε\varepsilon-Consensus-Halving for constant ε\varepsilon.

There exists a constant ε>0\varepsilon>0 for which it is NP-hard to decide if an ε\varepsilon-SC-Pizza-Sharing instance with nn mass distributions admits a solution consisting of an SC-path with at most n−2n-2 turns, even when every mass distribution is uniform over polynomially many rectangles, and there is no overlap between any two mass distributions.

3 Hardness of discrete Straight-Pizza-Sharing and SC-Pizza-Sharing

In this section, we study the discrete versions of Straight-Pizza-Sharing and SC-Pizza-Sharing.

R^{+} and R−R^{-} using at most nn lines such that for each i∈[2n]i\in[2n] it holds that ∣∣Pi∩R+∣−∣Pi∩R−∣∣≤ε⋅∣Pi∣\left||P_{i}\cap R^{+}|-|P_{i}\cap R^{-}|\right|\leq\varepsilon\cdot|P_{i}|. A point that is intersected by a line does not belong to any of R+R^{+}, R−R^{-}. Definition 11. For any n≥1n\geq 1, the problem ε\varepsilon-Discrete-SC-Pizza-Sharing is defined as follows: • Input: ε≥0\varepsilon\geq 0 and nn point sets P1,P2,…,PnP_{1},P_{2},\dots,P_{n} on 2^{2}. • Output: One of the following. (a) Two points with the same xx- or yy-coordinate. (b) A partition of 2^{2} to R+R^{+} and R−R^{-} using a yy-monotone SC-path with at most n−1n-1 turns such that for each i∈[n]i\in[n] it holds that ∣∣Pi∩R+∣−∣Pi∩R−∣∣≤ε⋅∣Pi∣\left||P_{i}\cap R^{+}|-|P_{i}\cap R^{-}|\right|\leq\varepsilon\cdot|P_{i}|. A point that is intersected by a line does not belong to any of R+R^{+}, R−R^{-}. Notice that the first kind of allowed output for both problems (Definition 10(a), Definition 11(a)) is a witness that the input points are not in general position or that their xx- or yy-coordinates are not unique, respectively, which can be checked in polynomial time. The second kind of output (Definition 10(b), Definition 11(b)) is the one that is interesting and can encode the hard instances studied here. In case the first kind of output does not exist, the other one is guaranteed to exist due to [Sch21] for ε\varepsilon-Discrete-Straight-Pizza-Sharing, while for ε\varepsilon-Discrete-SC-Pizza-Sharing its existence is guaranteed for every ε∈\varepsilon\in due to a reduction we present in Section 4.4 which shows containment in PPA.

The definition of ε\varepsilon-Discrete-Straight-Pizza-Sharing is a slightly modified form of the one that appears in [Sch21], where it is referred to as DiscretePizzaCutting. In our definition, we avoid having to “promise” an input of points that are in general position, by allowing as output a witness of an inappropriate input. Furthermore, the definition we present is more general, since it accommodates an approximation factor ε\varepsilon; in particular, for ε<min⁡i{1/∣Pi∣}\varepsilon<\min_{i}\{1/|P_{i}|\} we get the definition of the aforementioned paper. Similarly, we define ε\varepsilon-Discrete-SC-Pizza-Sharing, which to the best of our knowledge, has not been stated in previous work. Note that, if the input consists of points that are in general position, in an ε\varepsilon-Discrete-Straight-Pizza-Sharing solution only up to two points are allowed to be intersected by the same line, while in an ε\varepsilon-Discrete-SC-Pizza-Sharing solution only up to two points are allowed to be intersected by the same line segment.

Recall that in ε\varepsilon-SC-Pizza-Sharing we are given nn mass distributions, while in ε\varepsilon-Straight-Pizza-Sharing we are given 2n2n mass distributions as input. In Appendix C, we describe a general construction that takes as an input q∈{n,2n}q\in\{n,2n\} mass distributions μ1,…,μq\mu_{1},\dots,\mu_{q} normalized on 2^{2}, represented by weighted polygons with holes (see Appendix A for the detailed description), and turns it into qq sets of points P1,…,PqP_{1},\dots,P_{q} on 2^{2} whose union is in general position and furthermore, they have unique xx- and yy-coordinates. We prove that, given a set of m≤2nm\leq 2n lines or an SC-path of at most m′≤2n−1m^{\prime}\leq 2n-1 turns, that partition 2^{2} into R+R^{+} and R−R^{-} such that ∣∣Pi∩R+∣−∣Pi∩R−∣∣≤(ε−ε′)⋅∣Pi∣\left||P_{i}\cap R^{+}|-|P_{i}\cap R^{-}|\right|\leq(\varepsilon-\varepsilon^{\prime})\cdot|P_{i}|, the same lines and SC-path, respectively, separate the mass distributions such that ∣μi(R+)−μi(R−)∣≤ε\left|\mu_{i}(R^{+})-\mu_{i}(R^{-})\right|\leq\varepsilon.

Let N≥2nN\geq 2n be the input size of any of our two pizza sharing problems, and let the smallest area triangle in the mass distributions’ triangulation be α\alpha. For any ε′<ε\varepsilon^{\prime}<\varepsilon, where ε\varepsilon and α\alpha are at least inverse polynomial in NN, the construction results to a polynomial time reduction from ε\varepsilon-Straight-Pizza-Sharing to (ε−ε′)(\varepsilon-\varepsilon^{\prime})-Discrete-Straight-Pizza-Sharing and from ε\varepsilon-SC-Pizza-Sharing to (ε−ε′)(\varepsilon-\varepsilon^{\prime})-Discrete-SC-Pizza-Sharing. The reduction can be performed in time polynomial in the input size and in 1/α1/\alpha. It consists of two parts: (i) First we “pixelate” the mass distributions finely enough so that they are represented by a sufficiently large number of pixels. This will ensure a high enough “resolution” of the pixelated distributions. (ii) The pixels will then be turned into points, which we have to perturb in order to guarantee they are in general position and with unique coordinates, as required. In particular, in Appendix C, we prove the following.

Let NN be the input size of an approximate pizza sharing problem (either ε\varepsilon-Straight-Pizza-Sharing or ε\varepsilon-SC-Pizza-Sharing) whose triangulation has no triangle with area less than α>0\alpha>0. Also, let ε′∈[6Nc,ε)\varepsilon^{\prime}\in\left[\frac{6}{N^{c}},\varepsilon\right), where c>0c>0 is a fixed constant, and 6Nc<ε<1\frac{6}{N^{c}}<\varepsilon<1. Then, the instance can be reduced in time poly(N,1/α)\textup{poly}(N,1/\alpha) to its approximate discrete version, that is, (ε−ε′)(\varepsilon-\varepsilon^{\prime})-Discrete-Straight-Pizza-Sharing or (ε−ε′)(\varepsilon-\varepsilon^{\prime})-Discrete-SC-Pizza-Sharing.

Given Theorem 5 and Theorem 8, and since their instances are constructed such that α\alpha is an at least inverse polynomial function of the input size, the above lemma implies the following hardness results.

ε\varepsilon-Discrete-Straight-Pizza-Sharing with 2n2n point sets is PPA-hard for any constant ε<1/5\varepsilon<1/5, even when n+n1−δn+n^{1-\delta} lines are allowed for any given constant δ>0\delta>0.

ε\varepsilon-Discrete-SC-Pizza-Sharing with nn point sets is PPA-hard for any constant ε<1/5\varepsilon<1/5, even when n−1+n1−δn-1+n^{1-\delta} turns are allowed in the SC-path for any given constant δ>0\delta>0.

We note that PPA-hardness for ε\varepsilon-Discrete-Straight-Pizza-Sharing was so far known only for any ε∈[0,min⁡i{1/∣Pi∣})\varepsilon\in\left[0,\min_{i}\{1/|P_{i}|\}\right) (which is equivalent to ε=0\varepsilon=0), due to [Sch21]. Since, in the aforementioned paper’s constructions, min⁡i{1/∣Pi∣}∈O(1/poly(N))\min_{i}\{1/|P_{i}|\}\in O(1/\textup{poly}(N)), our result strengthens the hardness of the problem significantly.

Lemma 12 is general enough that allows us to derive NP-hardness results for the decision variants of the two discrete versions of the pizza sharing problems. If we ask for a solution with at most n−1n-1 straight lines or n−2n-2 turns in (ε−ε′)(\varepsilon-\varepsilon^{\prime})-Discrete-Straight-Pizza-Sharing and (ε−ε′)(\varepsilon-\varepsilon^{\prime})-Discrete-SC-Pizza-Sharing, respectively, then we can easily reduce to them from the instances of Theorem 6 and Theorem 9, picking ε′\varepsilon^{\prime} to be some inverse polynomial function of NN, e.g., ε′=1/N\varepsilon^{\prime}=1/N. In particular, we get the following.

There exists a constant ε>0\varepsilon>0 for which it is NP-hard to decide whether a solution of ε\varepsilon-Discrete-Straight-Pizza-Sharing with 2n2n point sets and at most n−1n-1 lines exists.

There exists a constant ε>0\varepsilon>0 for which it is NP-hard to decide whether a solution of ε\varepsilon-Discrete-SC-Pizza-Sharing with nn point sets and an SC-path with at most n−2n-2 turns exists.

4 Hardness of exact SC-Pizza-Sharing

In this section, we show hardness results for exact SC-Pizza-Sharing. We prove that solving SC-Pizza-Sharing is FIXP-hard and that deciding whether there exists a solution for SC-Pizza-Sharing with fewer than n−1n-1 turns is ETR-hard.

As mentioned earlier, computing an exact solution of a FIXP-hard problem may require computing an irrational number. To showcase this for SC-Pizza-Sharing, consider the following simple instance. Let us have a single mass distribution in the shape of a right-angled triangle, whose corners are on (0,1)(0,1), (1,1)(1,1), and (1,0)(1,0). It is normalised, i.e., its total mass is 11, therefore its weight is 22. An exact solution of SC-Pizza-Sharing is either a horizontal or a vertical straight line (with turns), which cuts the triangle such that each half-space has half of the mass, that is, 1/21/2. One can easily check that the solution is either the horizontal line y=2/2y=\sqrt{2}/2, or the vertical line x=2/2x=\sqrt{2}/2.

We provide a main reduction from (exact) Consensus-Halving to (exact) SC-Pizza-Sharing, whose gadgets are then employed to show the ETR-hardness of the decision version. This time, as a starting point, we will use the instances of Consensus-Halving produced in [DFMS21], which we will denote by I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}}. When clear from context, by I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}} we will denote a particular instance of the aforementioned family. We note that here the input consists of sets of points, i.e., we describe polygons by their vertices. For a detailed description of the input representation, see Appendix A. In [DFMS21], the FIXP-hard family of instances we reduce from has as input nn arithmetic circuits capturing the cumulative valuation of nn agents on $$, which were piece-wise polynomials of maximum degree 2. However, since their (density) valuation functions are piece-wise linear, the input of SC-Pizza-Sharing suffices to consist of only rectangles and triangles, and no other shape. Therefore, there is no need for extra translation of the input of Consensus-Halving to the input of SC-Pizza-Sharing.

The reduction. Here we show the main reduction, which will conclude with the proof that finding an exact solution to SC-Pizza-Sharing is FIXP-hard, and then will be used to show that its decision version is ETR-hard. We reduce from a Consensus-Halving instance I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}} with nn agents and kk-block-triangle valuations to an SC-Pizza-Sharing instance I\textscSCI_{\textsc{SC}} with nn mass distributions. When requesting a solution in I\textscSCI_{\textsc{SC}} with at most n−1n-1 turns in the SC-path, then we get FIXP-hardness, while when requesting to decide if I\textscSCI_{\textsc{SC}} is solvable with n−2n-2 turns, then we get ETR-hardness. Both of these results are due to [DFMS21], and hold even for 66-block-triangle valuations.

As in our previous proofs, for ease of presentation, the space of the instance is inflated to [0,2n(k+1)]2\left[0,2n(k+1)\right]^{2}, and by scaling the construction down to 2^{2} the correctness is attained. The key difference between this reduction and the previous reductions on the approximate versions is that the starting point of the reduction, i.e., instance I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}}, besides rectangular, also contains triangular-shaped valuations for agents. More specifically, all of the following hold:

the valuation function of every agent is 4-block-triangle, or 6-block;

every triangle has height 2 and belongs to exactly one interval of interest of the form [a,a+1][a,a+1];

for every agent i∈[n]i\in[n] there exists an interval [ai,bi][a_{i},b_{i}] that contains more than half of their total valuation, and in addition, for every i′≠ii^{\prime}\neq i we have (ai,bi)∩(ai′,bi′)=∅(a_{i},b_{i})\cap(a_{i^{\prime}},b_{i^{\prime}})=\emptyset.

Also, in this reduction, the resulting SC-Pizza-Sharing instances will contain weighted mass distributions (see definition in Section 2).

For each subinterval j∈[m]j\in[m] of I\textscCHI_{\textsc{CH}}, we will construct a tile tjt_{j} of size 1×11\times 1. Let the total value of agent ii in the jj-th interval of I\textscCHI_{\textsc{CH}} be vijv_{ij}. If for agent ii the valuation in the jj-th interval has constant density (i.e., the total valuation has a rectangular shape), then we create a square mass that belongs to μi\mu_{i}, of size 1×11\times 1 inside tjt_{j}, with weight wij=vijw_{ij}=v_{ij}. If her valuation has linear density (i.e., the total valuation has a triangular shape), then we create a mass of right-triangular shape that belongs to μi\mu_{i}, with two of its sides being of length 11 and touching the top and right sides of tjt_{j}, while its right angle touches the top-right corner of tjt_{j}, and its weight is wij=2w_{ij}=2. Observe that, for any given i∈[n]i\in[n], the total value of a subinterval in I\textscCHI_{\textsc{CH}} is equal to its total mass in the corresponding tile. Finally, we place the tiles sequentially in a diagonal manner, i.e., each tile tjt_{j} is axis-aligned and placed with its bottom-left corner at point (j,j)(j,j). See Figure 5 and Figure 6 for a depiction.

Now we need to show how an exact solution of I\textscSCI_{\textsc{SC}}, that is, an SC-path with n−1n-1 many turns is mapped back to a solution of I\textscCHI_{\textsc{CH}} with nn cuts. Consider the first tile tjt_{j} that is intersected by the SC-path, resulting in a part of it belonging to R+R^{+} and the rest of it to R−R^{-}. Let aj+:=μi(tj∩R+)a_{j}^{+}:=\mu_{i}(t_{j}\cap R^{+}), and by our definition above, it is implied that μi(tj∩R−)=vij−aj+\mu_{i}(t_{j}\cap R^{-})=v_{ij}-a_{j}^{+}. We then place a cut in the jj-th subinterval of I\textscCHI_{\textsc{CH}} at point xj′=xj+(xj+1−xj)⋅aj+vijx^{\prime}_{j}=x_{j}+(x_{j+1}-x_{j})\cdot\frac{a_{j}^{+}}{v_{ij}}, and label the interval [0,xj′][0,x^{\prime}_{j}] with “++”, and the interval [xj′,xj+1][x^{\prime}_{j},x_{j+1}] with “−-”. Next, consider the second tile tj′t_{j^{\prime}} that is intersected by the SC-path, resulting in a part of it belonging to R+R^{+} and the rest of it to R−R^{-}. For this tile, we will place a cut in the j′j^{\prime}-th subinterval of I\textscCHI_{\textsc{CH}} at point xj′′=xj′+(xj′+1−xj′)⋅vij′−aj′+vij′x^{\prime}_{j^{\prime}}=x_{j^{\prime}}+(x_{j^{\prime}+1}-x_{j^{\prime}})\cdot\frac{v_{ij^{\prime}}-a_{j^{\prime}}^{+}}{v_{ij^{\prime}}}, and label the interval [xj+1,xj′′][x_{j+1},x^{\prime}_{j^{\prime}}] with “−-”, and the interval [xj′′,xj′+1][x^{\prime}_{j^{\prime}},x_{j^{\prime}+1}] with “++”. We continue in this fashion with the translation of the SC-path into cuts and labels of I\textscCHI_{\textsc{CH}}, where, the rr-th in order intersected tile will be of the kind tjt_{j} if rr is odd, and of the kind tj′t_{j^{\prime}} if rr is even. This will result in a set of cuts that define regions with alternating signs.

It remains to prove that this is a solution to I\textscCHI_{\textsc{CH}}. To do this, we will use the following crucial observation.

In any solution of I\textscSCI_{\textsc{SC}} created by I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}}, there can be no turn inside a tile.

The truth of the statement can become apparent if one considers Item 3 of the above facts on I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}}, together with the fact that the endpoints of each [ai,bi][a_{i},b_{i}] are points of interest. It is implied then, that in any of the I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}} solutions, there needs to be at least one cut in each interval [ai,bi][a_{i},b_{i}] for each i∈[n]i\in[n]. And since we are allowed to draw at most nn cuts, there will be a single cut in each of those intervals. Also, due to the fact that ai,bia_{i},b_{i} are points of interest, each cut in [ai,bi][a_{i},b_{i}] belongs to a different subinterval, and therefore, there will be exactly nn cuts in nn distinct subintervals. Focusing now on our I\textscSCI_{\textsc{SC}} construction, the SC-path with n−1n-1 turns consists of a total of nn horizontal and vertical line segments. If any of those does not intersect any tile, then this SC-path will correspond to a set of at most n−1n-1 cuts in I\textscCHI_{\textsc{CH}}, which cannot be a solution. Therefore, every line segment intersects some tile.

Notice that, due to the diagonal placement of tiles, any SC-path that is a solution has to be xx-monotone and yy-monotone, i.e., to have a “staircase” form. Also, this diagonal placement of tiles dictates that, if there was a turn of the SC inside a tile, then two line segments are used to intersect it instead of one. This means that at most n−1n-1 tiles will be intersected, and therefore this translates to a set of at most n−1n-1 cuts in I\textscCHI_{\textsc{CH}}, which cannot be a solution. ∎

Put differently, having a turn inside a tile would be a “waste”, and in our instances, all turns are needed in order for a solution to exist.

Let us now focus on a tile trt_{r}. We consider two cases that can appear in I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}}:

For all i∈[n]i\in[n], viv_{i} is nonnegative constant in the rr-th subinterval. Let trt_{r} be a tile as the aforementioned tjt_{j}. Now notice that the “++” part of the jj-th subinterval equals vij⋅xj′−xjxj+1−xj=aj+=μi(tj∩R+)v_{ij}\cdot\frac{x^{\prime}_{j}-x_{j}}{x_{j+1}-x_{j}}=a_{j}^{+}=\mu_{i}(t_{j}\cap R^{+}), and the “−-” part equals vij⋅xj+1−xj′xj+1−xj=vij−aj+=μi(tj∩R−)v_{ij}\cdot\frac{x_{j+1}-x^{\prime}_{j}}{x_{j+1}-x_{j}}=v_{ij}-a_{j}^{+}=\mu_{i}(t_{j}\cap R^{-}). If on the other hand, trt_{r} is a tile as the aforementioned tj′t_{j^{\prime}}, it is again easy to see that, for every i∈[n]i\in[n], the total “++” and “−-” parts of tj′t_{j^{\prime}} have mass equal to the corresponding “++” and “−-” parts of the j′j^{\prime}-th subinterval.

There is an i∈[n]i\in[n] such that viv_{i} is linear in the rr-th subinterval. Due to the previous case, for all i′≠ii^{\prime}\neq i, the total “++” and “−-” parts of the tile have mass equal to the corresponding “++” and “−-” parts of the corresponding subinterval in I\textscCHI_{\textsc{CH}}. For ii, due to the structure of I\textscCH\textscDFMSI_{\textsc{CH}}^{\textsc{DFMS}}, the slope of viv_{i} is 22, and the length of the subinterval is 11, therefore we have vij=1v_{ij}=1. Let trt_{r} be a tile as the aforementioned tjt_{j}, and without loss of generality, let the SC-path intersect it horizontally at j+sj+s (y-coordinate), and its bottom part is labelled “++”, while the top part is labelled “−-”. According to the previous case, a cut should be placed at point xj′=xj+sx^{\prime}_{j}=x_{j}+s in the jj-th subinterval, and its part to the left should be labelled “++” while the one to its right should be labelled “−-”. Now notice that for agent ii the “++” part of the jj-th subinterval equals s⋅2s2=s2\frac{s\cdot 2s}{2}=s^{2}, while her “++” mass in tjt_{j} is s⋅s2⋅wij=s2\frac{s\cdot s}{2}\cdot w_{ij}=s^{2}, since wij=2w_{ij}=2. Similarly, both the “−-” parts of I\textscCHI_{\textsc{CH}} and I\textscSCI_{\textsc{SC}} are equal in the subinterval and the tile, respectively. Finally, if, trt_{r} is a tile as the aforementioned tj′t_{j^{\prime}}, it is again easy to see that the total “++” and “−-” parts of tj′t_{j^{\prime}} have mass equal to the corresponding “++” and “−-” parts of the j′j^{\prime}-th subinterval.

The above analysis shows that, for any i∈[n]i\in[n], in each individual tile the total “++” mass is equal to the total “++” value of the corresponding subinterval. Therefore, given a solution to I\textscSCI_{\textsc{SC}} where the total mass of R+R^{+} will equal that of R−R^{-} for every i∈[n]i\in[n], the induced cuts on I\textscCHI_{\textsc{CH}} constitute a solution. Finally, due to the FIXP-hardness of exact Consensus-Halving shown in [DFMS21], we get the following.

SC-Pizza-Sharing is FIXP-hard even when every mass distribution is uniform, consists of at most six pieces that can be unit-squares or right-angled triangles, and have overlap at most 3.

We can also show that deciding whether there exists an exact SC-Pizza-Sharing solution with n−2n-2 turns is ETR-hard. To show this, we will use a result of [DFMS21], where it was shown that deciding whether there exists an exact Consensus-Halving solution with nn agents and n−1n-1 cuts is ETR-hard. We give a reduction from this version of Consensus-Halving to the decision problem for SC-Pizza-Sharing. The reduction uses the same ideas that we presented for the FIXP-hardness reduction. The full details are in Appendix D, where the following theorem is shown.

It is ETR-hard to decide if an exact SC-Pizza-Sharing instance admits a solution with a SC-path with at most n−2n-2 turns, even when every mass distribution consists of at most six pieces that can be unit-squares or right-angled triangles, and have overlap at most 3.

Containment results

In this section, we present containment results for the exact and approximate versions of Straight-Pizza-Sharing and SC-Pizza-Sharing that we study in this paper. Our containment result for the former problem is possible by reducing it to its discrete version, which was recently shown by [Sch21] to be in PPA. For SC-Pizza-Sharing, our results revolve around a proof that solutions exist, which utilizes the Borsuk-Ulam theorem. A proof of this kind was already presented in [KRPS16], but we present a new one which is conceptually simpler, as it does not use any involved topological techniques. An advantage of our proof is that it can be made algorithmic, and so it can be used to show that SC-Pizza-Sharing is contained in BU. Then, by making simple modifications to the BU containment proof, we show that the other containment results hold. Then, we study the corresponding decision variants of the problems, and acquire containment in NP for approximate and discrete versions of the problems and ETR containment of exact SC-Pizza-Sharing.

Here we prove that ε\varepsilon-Straight-Pizza-Sharing with 2n2n mass distributions is in PPA for any ε∈Ω(1/poly(N))\varepsilon\in\Omega(1/\textup{poly}(N)) and α∈Ω(1/poly(N))\alpha\in\Omega(1/\textup{poly}(N)), where NN is the input size and α\alpha is the smallest area among the triangles of the triangulated mass distributions of ε\varepsilon-Straight-Pizza-Sharing. This answers a big open question left from [DFM22], where PPA containment was elusive. We also show that deciding whether a solution with at most n−1n-1 straight lines exists is in NP. Both those results are derived by reducing the problem to its discrete version, where instead of mass distributions, the input consists of points, and the goal is to bisect (up to one point) each of the 2n2n point sets using at most nn straight lines. The latter was recently shown by [Sch21] to be in PPA.

PPA containment. We will employ Lemma 12 in a straightforward way. In particular, given an ε\varepsilon-Straight-Pizza-Sharing instance for some ε∈Ω(1/poly(N))\varepsilon\in\Omega(1/\textup{poly}(N)), we will pick an ε′<ε\varepsilon^{\prime}<\varepsilon as prescribed in the aforementioned lemma, and reduce our problem to Discrete-Straight-Pizza-Sharing in time poly(N,1/α)\textup{poly}(N,1/\alpha).

ε\varepsilon-Straight-Pizza-Sharing with 2n2n weighted mass distributions with holes is in PPA for any ε∈Ω(1/poly(N))\varepsilon\in\Omega(1/\textup{poly}(N)) and α∈Ω(1/poly(N))\alpha\in\Omega(1/\textup{poly}(N)), where NN is the input size and α\alpha is the smallest area among the triangles of the triangulated mass distributions.

NP containment. Observe that Lemma 12 shows how to turn a set of 2n2n mass distributions (weighted polygons with holes) into a set of (unweighted) points, which, if cut by at most 2n2n straight lines, will result to an approximate cut of the mass distributions relaxed by an extra additive ε′∈Ω(1/Nc)\varepsilon^{\prime}\in\Omega(1/N^{c}) for any c>0c>0. When the number of straight lines is at most n−1n-1, this gives straightforwardly a reduction to -Discrete-Straight-Pizza-Sharing, since we can check how many of the polynomially many points of the latter instance are in R+R^{+} and R−R^{-}.

Deciding whether ε\varepsilon-Straight-Pizza-Sharing with 2n2n weighted mass distributions with holes has a solution with at most n−1n-1 straight lines is in NP for any ε∈Ω(1/poly(N))\varepsilon\in\Omega(1/\textup{poly}(N)) and α∈Ω(1/poly(N))\alpha\in\Omega(1/\textup{poly}(N)), where NN is the input size and α\alpha is the smallest area among the triangles of the triangulated mass distributions.

2 Containment of exact SC-Pizza-Sharing

Existence of a SC-Pizza-Sharing solution. We begin by proving that a solution to exact SC-Pizza-Sharing always exists. This proof holds for arbitrary mass distributions, but for our algorithmic results, we will only consider the case where the mass distributions are unions of polygons with holes. Our proof is based on the proof of Karasev, Roldán-Pensado, and Soberón [KRPS16], but they use more involved techniques from topology, which we would like to avoid since our goal is to implement the result algorithmically.

Figure 7 (see Appendix E) gives an overview for our SC-path embedding. The embedding considers only yy-monotone SC-paths. The path itself is determined by the points x1x_{1}, x2x_{2}, …, and y1y_{1}, y2y_{2}, …, which define the points at which the path turns. The path begins on the boundary on the line y1y_{1}, and then moves to x1x_{1}, at which points it turns and moves upwards to y2y_{2}, and it then turns to move to x2x_{2}, and so on.

But these points alone do not fully specify the path, since the decision of whether to move to the right or to the left when moving from xix_{i} to xi+1x_{i+1} has not been specified. To do this, we also include signs that are affixed to each horizontal strip. We will call the two sides of the cut side A (non-shaded) and side B (shaded). The (i+1)(i+1)-st strip [yi,yi+1][y_{i},y_{i+1}] is split into two by the vertical line segment on x=xix=x_{i} that starts from point (xi,yi)(x_{i},y_{i}) and ends at (xi,yi+1)(x_{i},y_{i+1}), with one part being on side A, and the other being on side B. If the strip is assigned the sign ++ then the area on the left of the strip is on side A of the cut, while if the strip is assigned the sign −- then the areas on the right of the strip is on side A of the cut. Once the sides of the cut have been decided, there is then a unique way to move through the points in order to separate sides A and B.

So, an SC-path with n−1n-1 turns, can be represented by nn variables and ⌈n/2⌉\left\lceil n/2\right\rceil signs. We then embed these into the sphere SnS^{n} in the following way. We give an informal description here, and the full definition will be given in the proof of Theorem 22. We encode the yy values as variables zi∈z_{i}\in where yi+1=yi+∣zi∣y_{i+1}=y_{i}+|z_{i}|, while the xx values are encoded as values in $,where, where|x_{i}|definesthedefines theithvalueontheth value on thexaxis.Weusethesignsoftheaxis. We use the signs of thez_{i}variablestodefinethesignsforthestrips,andnotethatatleastoneofthevariables to define the signs for the strips, and note that at least one of thez_{i}’sisnon−zerosincethesumofthestrips’lengthsshouldequal1.Finally,wethenshrinkallvariablessothattheirsumliesin’s is non-zero since the sum of the strips’ lengths should equal 1. Finally, we then shrink all variables so that their sum lies in,andweembedthemasapointin, and we embed them as a point inS^{n}$ using one extra dimension as a slack variable in case the sum of the absolute values of the variables does not equal 1.

The key property is that, if a point \vvP∈Sn\vv{P}\in S^{n} represents an SC-path, then the point −\vvP-\vv{P} represents exactly the same SC-path but with all signs flipped, and thus sides A and B will be exchanged. Therefore, we can write down a function ff, where f(\vvP)if(\vv{P})_{i} outputs the amount of mass of the ii-th mass distribution that lies on the A side of the cut, and therefore any point \vvP\vv{P} satisfying f(\vvP)=f(−\vvP)f(\vv{P})=f(-\vv{P}) must cut all mass distributions exactly in half.

So, we get the following theorem that is proved formally in LABEL:app:sc-pizza_exists.

BU-containment. The next step is to turn this existence result into a proof that SC-Pizza-Sharing is contained in BU. We begin by recapping the definition of BU given in [DFMS21]. An arithmetic circuit is a circuit that operates on real numbers, and uses gates from the set {c,+,−,×c,×,max⁡,min⁡}\{c,+,-,\times c,\times,\max,\min\}, where a cc-gate outputs the constant cc, a ×c\times c gate multiplies the input by a constant cc, and all other gates behave according to their standard definitions. The class BU contains every problem that can be reduced in polynomial time to Borsuk-Ulam.

The key issue is how to determine how much of a polygon lies on a particular side of the cut. To do this, we first triangulate all polygons, as shown in Figure 8(a) (see Section F.1) to obtain mass distributions that are unions of triangles. We then break each individual triangle down into a sum of axis-aligned right-angled triangles. This is shown in Figure 8(b), where the triangle \trnglABC\trngl{ABC} is represented as the sum of \trnglXYB+\trnglXBZ−\trnglAYB−\trnglXAC−\trnglCBZ\trngl{XYB}+\trngl{XBZ}-\trngl{AYB}-\trngl{XAC}-\trngl{CBZ}.

We then explicitly build an arithmetic circuit that, given the point x∈Sn+1x\in S^{n+1} used in Theorem 3, and a specific axis-aligned right-angled triangle, can output the amount of mass of that triangle that lies on side A of the cut. Then the function ff used in Theorem 3 can be built simply by summing over the right-angled triangles that arise from decomposing each of the polygons. So we get the following theorem, which is proved formally in LABEL:app:BU-inclusion.

Exact SC-Pizza-Sharing for weighted polygons with holes is in BU.

Deciding whether there exists an SC-path with kk turns that is an exact solution for SC-Pizza-Sharing with nn mass distributions is in ETR.

3 Containment of approximate SC-Pizza-Sharing

PPA containment. The following theorem shows PPA containment of ε\varepsilon-SC-Pizza-Sharing via a reduction to the ε\textsc−Borsuk−Ulam\varepsilon\textsc{-Borsuk-Ulam} problem which is in PPA [DFMS21].

ε\varepsilon-SC-Pizza-Sharing for weighted polygons with holes is in PPA.

Deciding whether there exists an SC-path with kk turns that is a solution of ε\varepsilon-SC-Pizza-Sharing with nn mass distributions is in NP.

4 Containment of discrete SC-Pizza-Sharing

It has already been shown in [Sch21] that ε\varepsilon-Discrete-Straight-Pizza-Sharing is in PPA even for ε=0\varepsilon=0. We complete the picture regarding inclusion of discrete pizza sharing problems, by showing that ε\varepsilon-Discrete-SC-Pizza-Sharing is also in PPA for ε=0\varepsilon=0, and therefore, for every ε∈\varepsilon\in. We will reduce ε\varepsilon-Discrete-SC-Pizza-Sharing to ε′\varepsilon^{\prime}-SC-Pizza-Sharing for ε=0\varepsilon=0 and ε′=1/2N\varepsilon^{\prime}=1/2N, where NN is the input size. Finally, as one would expect from the discrete version, its decision variant is in NP since its candidate solutions are verifiable in polynomial time.

PPA containment. Consider an instance of ε\varepsilon-Discrete-SC-Pizza-Sharing with point sets P1,…,PnP_{1},\dots,P_{n}, and denote P:=P1∪⋯∪PnP:=P_{1}\cup\dots\cup P_{n}. We intend to turn each point into a mass of non-zero area. To do so, we need to first scale down the landscape of the points, to create some excess free space as a “frame” around them. It suffices to scale down by an order of 3 and center it in the middle of 2^{2}, that is, to map each point (x,y)(x,y) to ( 13+x3,13+y3)\left(\ \frac{1}{3}+\frac{x}{3},\frac{1}{3}+\frac{y}{3}\right). Now all our points are in [1/3,2/3]2[1/3,2/3]^{2}. For convenience, for each i∈[n]i\in[n], we will be still denoting by PiP_{i} the new set of points after scaling and centering.

Now, we check whether for any pair i≠ji\neq j we have Pi,j:=Pi∩Pj≠∅P_{i,j}:=P_{i}\cap P_{j}\neq\emptyset, which means that two points of two point sets have identical positions. Consider all Pi,j≠∅P_{i,j}\neq\emptyset and let their union be P′P^{\prime}. Now consider all points that do not belong in P′P^{\prime}, that is R:=P∖P′R:=P\setminus P^{\prime}. We want to find the minimum positive difference in the xx- and yy-coordinates between any pair of points in RR. Let a point of RR be denoted pt=(xt,yt)p_{t}=(x_{t},y_{t}), and let us denote

and finally, d:=min⁡{xmin⁡,ymin⁡}d:=\min\{x_{\min},y_{\min}\}.

We now turn each point of PiP_{i} into an axis-aligned square of size d3×d3\frac{d}{3}\times\frac{d}{3}, with its bottom-left corner having the point’s coordinates. Notice that the total area of the squares is ∣Pi∣⋅d29|P_{i}|\cdot\frac{d^{2}}{9}, therefore, by setting the weight of each square to 9∣Pi∣d2\frac{9}{|P_{i}|d^{2}} we have the full description of a mass distribution μi\mu_{i}. Notice that, since d≤1/3d\leq 1/3, all of the mass distributions are in [1/3−1/3⋅3,2/3+1/3⋅3]=[2/9,7/9]2[1/3-1/3\cdot 3,2/3+1/3\cdot 3]=[2/9,7/9]^{2}.

Any SC-path that is a solution to the resulting ε\varepsilon-SC-Pizza-Sharing instance for ε=1/2N\varepsilon=1/2N, can be turned into a solution of ε′\varepsilon^{\prime}-Discrete-SC-Pizza-Sharing for ε′=0\varepsilon^{\prime}=0 in polynomial time.

If a horizontal (resp. vertical) line segment of the SC-path solution intersects two squares (of any two mass distributions), this means that their corresponding points in Discrete-SC-Pizza-Sharing had the same yy- (resp. xx-) coordinate. To see this, without loss of generality, suppose that the two squares are intersected by the same horizontal line segment, and that their corresponding points do not have the same yy-coordinate. Then, the distance between their bottom-left corners is positive but no greater than d/3d/3, which implies that d≤d/3d\leq d/3 (by definition of dd), a contradiction. A symmetric argument holds for two squares that are intersected by the same vertical line segment. In both the above cases, we return as a solution to -Discrete-SC-Pizza-Sharing the two corresponding points of the squares, which are of the kind of “Output (a)” in Definition 11.

Let Pmax⁡:=max⁡i∈[n]∣Pi∣P_{\max}:=\max_{i\in[n]}|P_{i}|. Suppose ∣Pi∣|P_{i}| is odd for every i∈[n]i\in[n]. Then, since we are asking for an 1/2N1/2N-SC-Pizza-Sharing solution, its SC-path cannot be non-intersecting with any of the squares, otherwise ∣μi(R+)−μi(R+)∣≥1/∣Pi∣≥1/Pmax⁡>1/2Pmax⁡≥1/2N|\mu_{i}(R^{+})-\mu_{i}(R^{+})|\geq 1/|P_{i}|\geq 1/P_{\max}>1/2P_{\max}\geq 1/2N, a contradiction. Therefore, at least one square of PiP_{i} is intersected, and this holds for every i∈[n]i\in[n]. If no line segment of the SC-path intersects two squares, we conclude that the SC-path with n−1n-1 turns and nn line segments will intersect at most nn squares. But since, as discussed above, at least one square of each mass distribution has to be intersected by SC, we get that SC intersects exactly one square of each mass distribution. Each side, R+R^{+} and R−R^{-} of the SC-path includes at least ⌊∣Pi∣/2⌋\left\lfloor|P_{i}|/2\right\rfloor entire squares for every i∈[n]i\in[n], and therefore, their bottom-left corners, i.e., the corresponding points of Discrete-SC-Pizza-Sharing. This is a solution to the -Discrete-SC-Pizza-Sharing instance.

Now suppose ∣Pi∣|P_{i}| is even for some i∈[n]i\in[n]. We can remove an arbitrary square from all mass distributions that come from point sets with even cardinality, and perform the aforementioned reduction to 1/2N1/2N-SC-Pizza-Sharing. Then, let us call “ii-th segment” the one that intersects one square of PiP_{i}, called ii-th square, and let it be a vertical segment, without loss of generality. By placing back the removed square, its bottom-left corner will be: (i) either on opposite sides with that of the ii-th square, (ii) or on the same side (notice that due to the allowed discrepancy 1/2N<1/2Pmax⁡1/2N<1/2P_{\max}, the ii-th segment cannot fall on the bottom-left corner of the ii-th square). Then, in case (i) each side contains exactly ∣Pi∣/2|P_{i}|/2 bottom-left corners of squares (i.e., points). By scaling up the positions of SC-path’s segments (recall that we have scaled down), this is a solution to the -Discrete-SC-Pizza-Sharing instance. In case (ii), suppose without loss of generality that both the inserted square and the bottom-left corner of the ii-th square are on the left side of the ii-th segment. We modify the SC-path by shifting the ii-th segment to the left such that it is now located d/3d/3 to the left of ii-th square’s bottom-left corner. Notice that this position is to the right of the inserted square since there is at least 2d/32d/3 distance between the two squares. We do the same for every i∈[n]i\in[n] has even number of points/squares. Then, each side of the SC-path for every i∈[n]i\in[n] has exactly ∣Pi∣/2|P_{i}|/2 bottom-left corners of squares which represent the initial points. After scaling up the positions of the SC-path’s segments, this is a solution to -Discrete-SC-Pizza-Sharing.

Finally, it is clear that the aforementioned operations can be performed in poly(N)\textup{poly}(N) time. ∎

By the PPA containment of 1/2N1/2N-SC-Pizza-Sharing (see Theorem 27), we get the following.

ε\varepsilon-Discrete-SC-Pizza-Sharing is in PPA for any ε∈\varepsilon\in.

NP containment. It is also easy to see that, by checking whether each of the points of each PiP_{i} is in R+R^{+} or R−R^{-} as defined by a candidate SC-path solution, we can decide in polynomial time if indeed it is a solution or not to -Discrete-SC-Pizza-Sharing (since the points are polynomially many in NN, by definition).

Deciding whether there exists an SC-path solution to ε\varepsilon-Discrete-SC-Pizza-Sharing is in NP for any ε∈\varepsilon\in.

Conclusions

For ε\varepsilon-Straight-Pizza-Sharing we have shown that finding a solution is PPA-complete for any ε∈[1/Nc,1/5)\varepsilon\in\left[1/N^{c},1/5\right), where NN is the input size and c>0c>0 is any constant. This result holds for both its continuous (even when the input contains only axis-aligned squares) and its discrete version. We have also shown that the same result holds for ε\varepsilon-SC-Pizza-Sharing, where the PPA containment holds even for inverse explonential ε\varepsilon. One open question that remains is “Can we prove containment in PPA of ε\varepsilon-Straight-Pizza-Sharing for inverse exponential ε\varepsilon?”. For the decision variant of both these problem, we show that there exists a small constant ε\varepsilon such that they are NP-complete. For both these problems and their search/decision variants, a big open question is “What is the largest constant ε\varepsilon for which the problem remains PPA-hard and NP-hard, respectively?”. The most interesting open question is “are there any good algorithms that guarantee a solution in polynomial time for some constant ε∈[1/5,1)\varepsilon\in[1/5,1), even when slightly more lines (resp. turns in an SC-path) are allowed?”.

We have also shown that exact SC-Pizza-Sharing is FIXP-hard and in BU. The interesting question that needs to be settled is “What is the complexity of the problem; is it complete for one of those classes or is it complete for some other class?”. Schnider in [Sch21] showed that exact Straight-Pizza-Sharing is FIXP-hard for a more general type of input than ours. So, a natural question is “When the input consists of weighted polygons with holes, is the problem FIXP-hard and inside BU, and if yes, is it complete for any of the two classes?”. For a strong approximation version of Consensus-Halving, [BHH21] showed that the problem is BUa\textup{{BU}}_{a}-complete. We conjecture that the same holds for the two pizza sharing problems studied here.

Another problem that remains open is the complexity of ε\varepsilon-Straight-Pizza-Sharing and ε\varepsilon-SC-Pizza-Sharing when every mass distribution consists of a constant number of non-overlapping rectangles. And finally, what is the complexity of the modified problem when instead of two parts we ask to fairly split the plane into d≥3d\geq 3 equal parts?

Appendix A Input representation

While the mathematical proof of Theorem 22 holds for general measures which are absolutely continuous with respect to the Lebesgue measure, for the computational problems Straight-Pizza-Sharing and SC-Pizza-Sharing we need a standardized way to describe the input, and therefore restrict to particular classes of measures. We consider the class of mass distributions that are defined by weighted polygons with holes. This class consists of mass distributions with the property that are succinctly represented in the input of a Turing machine.

Appendix B Proof of Lemma 4

Observe that the number of big-tiles intersected by a line in the set L\mathcal{L} is bounded by the number of times this line can cut y=x2y=x^{2}. Any straight line can cut y=x2y=x^{2} at most twice, and the same holds when instead of points we have big-tiles on y=x2y=x^{2} if we take care of the sparsity of the big-tiles sj=6d⋅js_{j}=\frac{6}{d}\cdot j.

Recall that we present the instance as if it was [0,(6d2)2+1]2\left[0,\left(\frac{6}{d^{2}}\right)^{2}+1\right]^{2} instead of 2^{2}. But furthermore, let us do the analysis of the proof in the artificial square [0,1d2+1]2\left[0,\frac{1}{d^{2}}+1\right]^{2} and show how from there we go to [0,(6d2)2+1]2\left[0,\left(\frac{6}{d^{2}}\right)^{2}+1\right]^{2}. In particular, in the artificial square the barycenters of consecutive big-tiles have distance 1 on the xx-axis, and their size is d6×d6\frac{d}{6}\times\frac{d}{6}. We will prove that for any big-tile of size 2r×2r2r\times 2r with r≤d/12r\leq d/12 there is no straight line that intersects more than two big-tiles.

For (x,y)=(j,j2)(x,y)=(j,j^{2}) the shortest distance is

We will now show that D(j,r)>r2D(j,r)>r\sqrt{2}, and therefore it does not intersect with the big-tile, since the greatest distance from a point of a big-tile to its barycenter is r2r\sqrt{2}. By the choice of r≤d/12r\leq d/12 we can see that 6jr+7r+r2≤16jr+7r+r^{2}\leq 1 and 2j+3−r>1+2r2j+3-r>1+2r, for every j∈[1d−2]j\in[\frac{1}{d}-2]. Therefore, we have

Let us denote by D(j,j′,r)D(j,j^{\prime},r) the distance of the statement. We have

From the latter two claims, we conclude the following.

For j=2j=2 and j′=1j^{\prime}=1, from Equation 1 for (x,y)=(j′,(j′)2)(x,y)=(j^{\prime},(j^{\prime})^{2}) we get

By the choice of r=d/12r=d/12, the numerator is strictly greater than 22, and 1+2r<7−r1+2r<7-r. Therefore, D(2,1,r)>2(7−r)2D(2,1,r)>\frac{2}{(7-r)\sqrt{2}}. Similarly to the proof of 32, we require D(2,1,r)>r2D(2,1,r)>r\sqrt{2}, which is true since 2(7−r)2≥r2\frac{2}{(7-r)\sqrt{2}}\geq r\sqrt{2}.

For j=3j=3 and j′=2j^{\prime}=2, from Equation 1 for (x,y)=(j′,(j′)2)(x,y)=(j^{\prime},(j^{\prime})^{2}) we get

Similarly to the previous case, the numerator is strictly greater than 22, and 1+2r<9−r1+2r<9-r. Therefore, D(3,2,r)>2(9−r)2D(3,2,r)>\frac{2}{(9-r)\sqrt{2}}. We require D(3,2,r)>r2D(3,2,r)>r\sqrt{2}, which is true since 2(9−r)2≥r2\frac{2}{(9-r)\sqrt{2}}\geq r\sqrt{2}.

For j=3j=3 and j′=1j^{\prime}=1, from Equation 1 for (x,y)=(j′,(j′)2)(x,y)=(j^{\prime},(j^{\prime})^{2}) we get

The numerator is strictly greater than 55, and 9−r>1+2r9-r>1+2r. Therefore, D(3,2,r)>5(9−r)2D(3,2,r)>\frac{5}{(9-r)\sqrt{2}}. We require D(3,1,r)>r2D(3,1,r)>r\sqrt{2}, which is true since 5(9−r)2≥r2\frac{5}{(9-r)\sqrt{2}}\geq r\sqrt{2}. ∎

Any line that intersects two big-tiles cannot intersect any big-tile located below the lowest of the two.

Now we are ready to prove the main claim.

Any line can intersect at most two big-tiles.

From Corollary 34, 35, 36 and the above paragraph, we have deduced Corollary 37. To prove the current claim it remains to prove that any line that intersects two big-tiles cannot intersect any big-tile between the two or above the highest of the two. However, this is easy to show by just shift in perspective. In particular, for the sake of contradiction, consider three big-tiles namely the jj-th, (j+k)(j+k)-th and (j+k′)(j+k^{\prime})-th with 1≤k<k′1\leq k<k^{\prime} that are intersected by a line. But then, by considering the latter two big-tiles and applying to them Corollary 37 we see that our assumption about the line intersecting the jj-th big-tile is contradicting the corollary. This completes the proof. ∎

The instance we have created is in square [0,1d2+1]2\left[0,\frac{1}{d^{2}}+1\right]^{2}. It is now easy to re-normalize it in [0,(6d2)2+1]2\left[0,\left(\frac{6}{d^{2}}\right)^{2}+1\right]^{2} by inflating it such that the big-tiles have size 1×11\times 1 instead of 2r×2r2r\times 2r. This results in the barycenters of consecutive big-tiles now having distance 1/(2r)=6/d1/(2r)=6/d instead of 11 in the xx-axis, since we picked r=d/12r=d/12. Therefore, the barycenters of the big-tiles for j∈[1/d]j\in[1/d] have coordinates (sj,sj2)(s_{j},{s_{j}}^{2}), where sj=6d⋅js_{j}=\frac{6}{d}\cdot j. Similarly, by further scaling the instance we can put it in 2^{2}.

We conclude that each cut can intersect at most two big-tiles, thus our available cuts will intersect at most 2(n+n1−δ)2(n+n^{1-\delta}) big-tiles. Let us disregard the big-tiles that are intersected by the set of lines L\mathcal{L}. By doing so, each cut introduces at most 2⋅cmax⁡⋅d2\cdot c_{\max}\cdot d discrepancy for each valuation viv_{i}, i∈[2n]i\in[2n], resulting to 2(n+n1−δ)⋅cmax⁡⋅d2(n+n^{1-\delta})\cdot c_{\max}\cdot d discrepancy in a total.

Now, we are ready to define the cuts for the Consensus-Halving instance I\textscCHI_{\textsc{CH}}. We consider the big-tiles in sequential order and we add one cut at j⋅dj\cdot d whenever we find two big-tiles tjt_{j} and tj+2t_{j+2} that belong to different regions, i.e. “++”, “−-”, and vice versa. This change of region can happen at most 2(n+n1−δ)2(n+n^{1-\delta}) times as argued earlier. Hence, we have at most 2(n+n1−δ)2(n+n^{1-\delta}) cuts in the instance I\textscCHI_{\textsc{CH}}. The labels of the pieces for I\textscCHI_{\textsc{CH}} follow the labels of the big-tiles of the instance IPI_{P} (excluding those who contain intersected small-tiles, which are arbitrarily given the same label as that of the next big-tile).

Let us denote by I+\mathcal{I}^{+} and I−\mathcal{I}^{-} the regions in I\textscCHI_{\textsc{CH}} corresponding (according to the above mapping) to the regions R+R^{+} and R−R^{-} of IPI_{P} respectively, induced by L\mathcal{L} after disregarding the intersected big-tiles. Then for the discrepancy in the valuation of agent ii in I\textscCHI_{\textsc{CH}} we get

where the last inequality comes from the fact that in our reduction we can pick d≤ε′4n⋅cmax⁡d\leq\frac{\varepsilon^{\prime}}{4n\cdot c_{\max}}. Recall that this is always possible since ε′≥1/nr\varepsilon^{\prime}\geq 1/n^{r} for some r≥1r\geq 1, and additionally, according to our reduction, dd is required to be at most 14n2⋅cmax⁡\frac{1}{4n^{2}\cdot c_{\max}}. This completes the proof.

Appendix C Proof of Lemma 12

Pixelation. We will start with the task of pixelating the mass distributions. Consider the input of an ε\varepsilon-SC-Pizza-Sharing or an ε\varepsilon-Straight-Pizza-Sharing instance, that is, q∈{n,2n}q\in\{n,2n\} mass distributions, respectively, on 2^{2} consisting of weighted polygons with holes (see Appendix A for details on the input representation). Let the instance’s input size be N≥2nN\geq 2n (by definition). As a first step, we perform a pixelation procedure: each polygon will be turned into a union of smaller squares that will have approximately the same total area as the polygon.

As shown in Section F.1, it is possible to decompose a polygon into a union of disjoint non-obtuse triangles. Each of those triangles’ area is rational since it can be computed by adding and subtracting the areas of five right-angled triangles with rational coordinates, which additionally, are axis-aligned. Furthermore, the cardinality of the non-obtuse triangles is a polynomial function in the input size of the polygon’s description, that is, the coordinates of the points that define its corners and the value that defines its weight. Therefore, the exact area of any mass distribution can be computed in polynomial time.

As we have discussed earlier, in order for our approximation parameter ε∈\varepsilon\in to make sense, we consider normalized mass distributions, meaning that μi(2)=1\mu_{i}\left(^{2}\right)=1 for all i∈[q]i\in[q]. Notice that it can be the case that some mass distributions have total area constant (in which case their weight is constant), while some others might have total area exponentially small (and therefore exponentially large weight). Therefore, in our analysis, we make sure that the “resolution” we provide to any polygon FF is relative to its actual area area(F)\text{area}(F) rather than its measure μi(F)\mu_{i}(F).

Consider some polygon FF of mass distribution μi\mu_{i} for i∈[q]i\in[q] and let it be triangulated into non-obtuse triangles, all with non-zero Lebesgue measure (strictly positive area). We will focus on one of FF’s non-obtuse triangles, \trnglABC\trngl{ABC} (see Figure 8(b)) with area S:=area(\trnglABC)>0S:=\text{area}\left(\trngl{ABC}\right)>0 and perimeter T>0T>0. Notice that since \trnglABC\trngl{ABC} is in 2^{2}, we have S≤1S\leq 1 and T≤3⋅2<5T\leq 3\cdot\sqrt{2}<5. Suppose that among all triangles of all μi\mu_{i}’s, the minimum area triangle has area α\alpha. The first step is to pixelate \trnglABC\trngl{ABC}. Let our pixels have size t×tt\times t for t=α/15N1+ct=\alpha/15N^{1+c}, where c>0c>0 is any fixed constant, and consider an axis-aligned square grid of pixels in our space, 2^{2}. We create a pixel for \trnglABC\trngl{ABC} if and only if the pixel’s intersection with \trnglABC\trngl{ABC} has non-zero Lebesgue measure. Then, the pixelated version of the triangle, denoted \trnglABCp\trngl{ABC}_{p}, has area S+S′S+S^{\prime}, where S′S^{\prime} is the excess area induced by the pixels intersected by the three sides of the triangle. By definition, S′S^{\prime} is at least , and at most the area of pixels that intersect the three sides of \trnglABC\trngl{ABC}. Therefore, by denoting the number of such pixels for each side by nAB,nBC,nCAn_{AB},n_{BC},n_{CA} and referring to Figure 8(b), we have nAB≤⌈AYt⌉+1+⌈YBt⌉+1≤AYt+YBt+4n_{AB}\leq\left\lceil\frac{AY}{t}\right\rceil+1+\left\lceil\frac{YB}{t}\right\rceil+1\leq\frac{AY}{t}+\frac{YB}{t}+4, and similarly for nBC,nCAn_{BC},n_{CA}. This gives

where the last strict inequality comes from the fact that t<T12t<\frac{T}{12}. To see this, we have to first notice that t=α15N1+c≤S15<S3Tt=\frac{\alpha}{15N^{1+c}}\leq\frac{S}{15}<\frac{S}{3T}. Then we also have to use the known formula that connects SS and TT, namely,

which implies ST<T4\frac{S}{T}<\frac{T}{4}, and therefore t<T12t<\frac{T}{12}.

We want to bound the proportion of excess area due to the pixelation compared to the triangle’s actual area, that is, S′/SS^{\prime}/S. We have

The pixelation of \trnglABC\trngl{ABC} results in \trnglABCp\trngl{ABC}_{p}, where area(\trnglABCp)<(1+αS⋅N1+c)⋅area(\trnglABC)\text{area}\left(\trngl{ABC}_{p}\right)<\left(1+\frac{\alpha}{S\cdot N^{1+c}}\right)\cdot\text{area}\left(\trngl{ABC}\right).

The total area of MM is at most α/5N1+c\alpha/5N^{1+c}.

For any disjoint ML,MRM_{L},M_{R} with ML∪MR=MM_{L}\cup M_{R}=M, we have ∣(area(L)−area(R))−(area(Lp∪ML)−area(Rp∪MR))∣≤2α/N1+c|(\text{area}(L)-\text{area}(R))-(\text{area}(L_{p}\cup M_{L})-\text{area}(R_{p}\cup M_{R}))|\leq 2\alpha/N^{1+c}.

It suffices to show that 0≤area(Lp∪ML)−area(L)≤2α/N1+c0\leq\text{area}(L_{p}\cup M_{L})-\text{area}(L)\leq 2\alpha/N^{1+c}. The first inequality is easy to see, since L⊆Lp∪MLL\subseteq L_{p}\cup M_{L}, which implies area(L)≤area(Lp∪ML)\text{area}(L)\leq\text{area}(L_{p}\cup M_{L}). For the second inequality, we have

Similarly, 0≤area(Rp∪MR)−area(R)≤2α/N1+c0\leq\text{area}(R_{p}\cup M_{R})-\text{area}(R)\leq 2\alpha/N^{1+c}, or equivalently, −2α/N1+c≤−(area(Rp∪MR)−area(R))≤0-2\alpha/N^{1+c}\leq-(\text{area}(R_{p}\cup M_{R})-\text{area}(R))\leq 0. Therefore,

Recall that, in an ε\varepsilon-Straight-Pizza-Sharing solution, at most 2n2n lines can intersect \trnglABC\trngl{ABC} and \trnglABCp\trngl{ABC}_{p} (even though the standardized version of the problem requires at most nn straight lines, as we showed in Theorem 5, PPA-hardness holds even for at most n+n1−δn+n^{1-\delta} lines for any constant δ>0\delta>0). Similarly for an ε\varepsilon-SC-Pizza-Sharing solution, since its SC-path comprises of at most n−1n-1 turns, i.e., nn straight line segments (and again, by Theorem 8 PPA-hardness holds even for at most n+n1−δn+n^{1-\delta} line segments for any constant δ>0\delta>0) By inductively applying 40 and 41 NN times (recall that N≥2nN\geq 2n), we get the following.

Let at most 2n2n straight lines intersect \trnglABCp\trngl{ABC}_{p}, and the side of each pixel be α/15N1+c\alpha/15N^{1+c} for any c>0c>0. Also, let MM be the set of its pixels that are intersected by the lines, and ML,MRM_{L},M_{R} be an arbitrary partition of MM. Then, area(M)≤α/5Nc\text{area}(M)\leq\alpha/5N^{c}, and furthermore, ∣(area(L)−area(R))−(area(Lp∪ML)−area(Rp∪MR))∣≤2α/Nc|(\text{area}(L)-\text{area}(R))-(\text{area}(L_{p}\cup M_{L})-\text{area}(R_{p}\cup M_{R}))|\leq 2\alpha/N^{c}, for any c>0c>0.

Turning pixels into points in general position. So far, for ε\varepsilon-SC-Pizza-Sharing and ε\varepsilon-Straight-Pizza-Sharing, with q∈{n,2n}q\in\{n,2n\} mass distributions, respectively, we have described how to turn each distribution μi\mu_{i}, i∈[q]i\in[q] into a pixelated version of it. We will now turn each of those pixels into a set of points. Let μi\mu_{i} consist of bb many polygons, and recall that each polygon has its own weight wi,j>0w_{i,j}>0, j∈[b]j\in[b]. Let Ti,j∈FiT^{i,j}\in\mathcal{F}_{i} be a non-obtuse triangle belonging to the jj-th polygon of μi\mu_{i}, and Fi\mathcal{F}_{i} be the set of these triangles composing μi\mu_{i}. Suppose a pixel contains non-zero Lebesgue measure of triangles {Ti,k}k∈DT^{i,k}\}_{k\in D} for some D⊆[b]D\subseteq[b], and let us denote wmax⁡i:=max⁡k∈Dwi,kw^{i}_{\max}:=\max_{k\in D}w_{i,k}. Observe that, due to our assumption that the mass distributions are normalised, we have 1=∑Ti,j∈Fiwi,j⋅area(Ti,j)≥∑Ti,j∈Fiwi,j⋅α1=\sum_{T^{i,j}\in\mathcal{F}_{i}}w_{i,j}\cdot\text{area}(T^{i,j})\geq\sum_{T^{i,j}\in\mathcal{F}_{i}}w_{i,j}\cdot\alpha, therefore,

We will place ⌈wmax⁡i⋅Nc⌉\left\lceil w^{i}_{\max}\cdot N^{c}\right\rceil points at the pixel’s bottom-left corner, that is, all having the same position. Each of the pixels of mass distribution ii has no weight, as desired, and they form a set PiP_{i}. Recall that each of the PiP_{i}’s we created contains at most 2Nc/t2α=1800N2+3c/a32N^{c}/t^{2}\alpha=1800N^{2+3c}/a^{3} points, i.e., polynomially many in the instance’s description size and 1/α1/\alpha. That is because its pixels can be at most ⌈1/t⌉⋅⌈1/t⌉≤4/t2=900N2+2c/a2\left\lceil 1/t\right\rceil\cdot\left\lceil 1/t\right\rceil\leq 4/t^{2}=900N^{2+2c}/a^{2}, with each pixel containing at most ⌈wmax⁡i⋅Nc⌉≤Ncα+1≤2Ncα\left\lceil w^{i}_{\max}\cdot N^{c}\right\rceil\leq\frac{N^{c}}{\alpha}+1\leq\frac{2N^{c}}{\alpha} points.

Observe, however, that in the discrete version of the pizza sharing instance we created, the points of the q∈{n,2n}q\in\{n,2n\} point sets lie on vertices of a square grid with edge length t=α/15N1+ct=\alpha/15N^{1+c}. These points are not in general position, and therefore, not all solutions of that instance (a set of three points intersected by the same straight line) can be translated back to a solution of the original corresponding continuous version of the instance. That is why we need to turn this instance into one where the points in P1∪⋯∪PqP_{1}\cup\dots\cup P_{q} are in general position.

First, we have to slightly shift the points inside each PiP_{i} which have identical positions. Notice that, by Equation 3, at most ⌈wmax⁡i⋅Nc⌉≤2Nc/α\left\lceil w^{i}_{\max}\cdot N^{c}\right\rceil\leq 2N^{c}/\alpha such points have identical positions. Let these points be p1,…,pmp_{1},\dots,p_{m} for m≤2Nc/αm\leq 2N^{c}/\alpha. We shift vertically each point pj=(xj,yj)p_{j}=(x_{j},y_{j}), j∈[m]j\in[m] to (xj,yj+(j−1)⋅α230N2+2c2N)\left(x_{j},y_{j}+(j-1)\cdot\frac{\alpha^{2}}{30N^{2+2c}2^{N}}\right). It is easy to see that none of these points have the same position. Observe also, that none of those points acquired the same position as another one from PiP_{i}, since the closest point of PiP_{i} to any of the points p1,…,pmp_{1},\dots,p_{m} for m≤2Nc/αm\leq 2N^{c}/\alpha has distance (in the yy-axis) at least (yj+α15N1+c)−(yj+m⋅α230N2+2c2N)≥α15N1+c(1−1N2N)>0\left(y_{j}+\frac{\alpha}{15N^{1+c}}\right)-\left(y_{j}+m\cdot\frac{\alpha^{2}}{30N^{2+2c}2^{N}}\right)\geq\frac{\alpha}{15N^{1+c}}\left(1-\frac{1}{N2^{N}}\right)>0.

Then, we have to take care of some points of Pi,PjP_{i},P_{j} for i≠ji\neq j, which might have identical positions. To prevent such a case, we create a small gap between the position of any two point sets’ points, by shifting their grid in the xx-direction. Formally, for each i∈[2n]i\in[2n], each point of PiP_{i} with coordinates (xi,yi)(x_{i},y_{i}) now acquires coordinates (xi+(i−1)⋅α230N2+2c2N,yi)\left(x_{i}+(i-1)\cdot\frac{\alpha^{2}}{30N^{2+2c}2^{N}},y_{i}\right). It is immediate that two points that used to have identical positions now have different positions; it is also straightforward that we have preserved the uniqueness of position among the points inside each PiP_{i}. What remains is to check whether a shifted point of PiP_{i} has the same position as a point of PjP_{j}. There are two cases: if they used to have the same yy-coordinate, then their minimum distance in the xx-coordinate is at least (xi+α15N1+c)−(xi+(2n−1)⋅α230N2+2c2N)>α15N1+c(1−αNc2N)>0\left(x_{i}+\frac{\alpha}{15N^{1+c}}\right)-\left(x_{i}+(2n-1)\cdot\frac{\alpha^{2}}{30N^{2+2c}2^{N}}\right)>\frac{\alpha}{15N^{1+c}}\left(1-\frac{\alpha}{N^{c}2^{N}}\right)>0; if they used to have the same xx-coordinate, then they must have had different yy-coordinates, which remained the case after the shifting. From the shifting we performed, we get the following.

The points of P1∪⋯∪PqP_{1}\cup\dots\cup P_{q} can lie only on vertices of a square grid in 2^{2} with edge length t′=α230N2+2c2Nt^{\prime}=\frac{\alpha^{2}}{30N^{2+2c}2^{N}}.

Now that each of the points of our instance has a unique position, we need to make sure that they are in general position, meaning that no triplet of them can be intersected by the same straight line. We will do this by further distorting the position of the points in a way, such that, even though the points will have a tiny distance from their original place (and therefore seem to remain on a square grid), they will be in general position.

To make the presentation of the analysis easier, instead of our fine grid (of 43), we will show the proof of correctness of our construction for the inflated square grid {0,…,k}×{0,…,k}\{0,\dots,k\}\times\{0,\dots,k\} with edge length 11. Observe, that in such a grid, the only lines that intersect at least two points have a slope of the form S=Y2−Y1X2−X1S=\frac{Y_{2}-Y_{1}}{X_{2}-X_{1}}, where (X1,Y1)(X_{1},Y_{1}), (X2,Y2)(X_{2},Y_{2}) is any pair of points. We will call those lines of interest. Since X1,Y1,X2,Y2∈{0,…,k}X_{1},Y_{1},X_{2},Y_{2}\in\{0,\dots,k\}, we deduce that the slopes of interest are actually of the form S=ΔYΔXS=\frac{\Delta Y}{\Delta X}, where ΔY∈{−k,…,0,…,k}\Delta Y\in\{-k,\dots,0,\dots,k\}, and ΔX∈{−k,…,−1,1,…,k}\Delta X\in\{-k,\dots,-1,1,\dots,k\}, or it is ∞\infty.

Consider now replacing each point in the grid, with a disk of radius d>0d>0 centered at its corresponding point’s position. Let us identify each disk by its center, that is, its initial point on the grid. For convenience, we will be calling the grid of points point-grid and the grid of disks disk-grid. Consider now a set RdR_{d} of at least two disks that can be intersected by the same straight line, called distorted line of interest, and let R′R^{\prime} be the collection of all such sets (and again, R′R^{\prime} is a finite collection since kk is finite).

S=S′S=S^{\prime}: Then we have ∣tan⁡ψ∣≥∣Y+1−d−(Y+d)k+d−(0−d)∣=1−2dk+2d≥7d>∣tan⁡θ∣|\tan{\psi}|\geq\left|\frac{Y+1-d-(Y+d)}{k+d-(0-d)}\right|=\frac{1-2d}{k+2d}\geq 7d>|\tan{\theta}|, where the last weak inequality holds for any d≤142k4d\leq\frac{1}{42k^{4}}.

S≠S′S\neq S^{\prime}: Without loss of generality, let φ,φ′∈[0,2π)\varphi,\varphi^{\prime}\in[0,2\pi), and let ρ:=min⁡{∣φ−φ′∣,∣φ−φ′−π∣}\rho:=\min\{|\varphi-\varphi^{\prime}|,|\varphi-\varphi^{\prime}-\pi|\}. Then ψ≥ρ−2⋅θ\psi\geq\rho-2\cdot\theta, therefore, ∣tan⁡ψ∣≥∣tan⁡(ρ−2θ)∣=∣tan⁡ρ−tan⁡2θ∣∣1+tan⁡ρ⋅tan⁡2θ∣|\tan{\psi}|\geq|\tan{(\rho-2\theta)}|=\frac{|\tan{\rho}-\tan{2\theta}|}{|1+\tan{\rho}\cdot\tan{2\theta}|}. Furthermore, we have that tan⁡ρ≤1k−0=1k\tan{\rho}\leq\frac{1}{k}-0=\frac{1}{k}, and that tan⁡2θ=2⋅tan⁡θ1−tan⁡2θ≤14d1−49d2≤41d≤4142k2<1\tan{2\theta}=\frac{2\cdot\tan{\theta}}{1-\tan^{2}{\theta}}\leq\frac{14d}{1-49d^{2}}\leq 41d\leq\frac{41}{42k^{2}}<1. This implies that ∣tan⁡ψ∣≥∣tan⁡ρ−tan⁡2θ∣1+1/k≥∣tan⁡ρ−tan⁡2θ∣2≥∣tan⁡ρ∣2|\tan{\psi}|\geq\frac{|\tan{\rho}-\tan{2\theta}|}{1+1/k}\geq\frac{|\tan{\rho}-\tan{2\theta}|}{2}\geq\frac{|\tan{\rho}|}{2}. Also, for finite S:=tan⁡φ,S′:=tan⁡φ′S:=\tan{\varphi},S^{\prime}:=\tan{\varphi^{\prime}}, we know that ∣tan⁡φ∣≤k1|\tan{\varphi}|\leq\frac{k}{1}, ∣tan⁡φ′∣≤k1|\tan{\varphi^{\prime}}|\leq\frac{k}{1}, therefore ∣1+tan⁡φ⋅tan⁡φ′∣≤1+k2|1+\tan{\varphi}\cdot\tan{\varphi^{\prime}}|\leq 1+k^{2}. We also have that ∣S−S′∣=∣ΔYΔX−ΔY′ΔX′∣=∣ΔY⋅ΔX′−ΔY′⋅ΔXΔX⋅ΔX′∣≥1k2|S-S^{\prime}|=\left|\frac{\Delta Y}{\Delta X}-\frac{\Delta Y^{\prime}}{\Delta X^{\prime}}\right|=\left|\frac{\Delta Y\cdot\Delta X^{\prime}-\Delta Y^{\prime}\cdot\Delta X}{\Delta X\cdot\Delta X^{\prime}}\right|\geq\frac{1}{k^{2}}, since ∣ΔY⋅ΔX′−ΔY′⋅ΔX∣|\Delta Y\cdot\Delta X^{\prime}-\Delta Y^{\prime}\cdot\Delta X| is a positive integer, and ∣ΔX⋅ΔX∣≤k⋅k|\Delta X\cdot\Delta X|\leq k\cdot k. By definition of ρ\rho, we have

and therefore, ∣tan⁡ψ∣≥14k4≥7d>∣tan⁡θ∣|\tan{\psi}|\geq\frac{1}{4k^{4}}\geq 7d>|\tan{\theta}|, where the last weak inequality holds for any d≤142k4d\leq\frac{1}{42k^{4}}.

For both of the above cases, we have shown that the maximum difference θ\theta in the angle of distorted lines of interest, as compared to their non-distorted version, is strictly smaller than the minimum angle ψ\psi between any two distorted lines of interest that intersect two pairs of disks, where each pair defines a different line of interest. ∎

The above claim shows that, for small enough disk radius dd, the only disks that can be intersected by the same line, are only those whose centers can be intersected by a common line. Recall now that our ultimate goal is to move each point of the {0,…,k}×{0,…,k}\{0,\dots,k\}\times\{0,\dots,k\} grid slightly, so that they are in general position. To do that, we will properly move each point inside the first quadrant of its respective disk with radius dd, so that this extra distortion results in a general position.

Since each “skewed” point will be in its respective disk, by 44, if there are three points that can be intersected by a common line, the centers of their respective disks should lie on the same line of interest. Therefore, we only need to prove that no three skewed points whose original position was on the same line of interest can be intersected by a common line. We will denote by Xi∈{0,…,k}X_{i}\in\{0,\dots,k\} and Yi∈{0,…,k}Y_{i}\in\{0,\dots,k\} the coordinates of an original point on the grid, and by xi,yix_{i},y_{i} the skewed points after moving the original ones inside their respective disk. Let r:=12kr:=\frac{1}{2k}, and δi:=(Xi2+Yi2)⋅d4k2≤d2\delta_{i}:=(X_{i}^{2}+Y_{i}^{2})\cdot\frac{d}{4k^{2}}\leq\frac{d}{2}. The “skewing” will be as follows: xi:=Xi+δix_{i}:=X_{i}+\delta_{i}, and yi:=Yi+r⋅δiy_{i}:=Y_{i}+r\cdot\delta_{i}.

The skewed points are in general position.

We define the auxiliary quantities W1:=Y1−Y0W_{1}:=Y_{1}-Y_{0}, W2:=Y2−Y0W_{2}:=Y_{2}-Y_{0}, Z1:=X1−X0Z_{1}:=X_{1}-X_{0}, Z2:=X2−X0Z_{2}:=X_{2}-X_{0}, D1:=Z12+2X0Z1+W12+2Y0W1D_{1}:=Z_{1}^{2}+2X_{0}Z_{1}+W_{1}^{2}+2Y_{0}W_{1}, D2:=Z22+2X0Z2+W22+2Y0W2D_{2}:=Z_{2}^{2}+2X_{0}Z_{2}+W_{2}^{2}+2Y_{0}W_{2}. Let us analyse the case where SS is finite. Then S=W1Z1=W2Z2=Y2−Y1X2−X1S=\frac{W_{1}}{Z_{1}}=\frac{W_{2}}{Z_{2}}=\frac{Y_{2}-Y_{1}}{X_{2}-X_{1}}. We have

where the last equalities hold, since S<∞S<\infty, and therefore, Z1≠0Z_{1}\neq 0, Z2≠0Z_{2}\neq 0. We have assumed that the three skewed points are intersected by a common line, which implies that y1−y0x1−x0=y2−y0x2−x0\frac{y_{1}-y_{0}}{x_{1}-x_{0}}=\frac{y_{2}-y_{0}}{x_{2}-x_{0}}. From Equation 4, Equation 5, we get

which, after simplification, gives (S−r)⋅D1Z1=(S−r)⋅D2Z2(S-r)\cdot\frac{D_{1}}{Z_{1}}=(S-r)\cdot\frac{D_{2}}{Z_{2}}, or equivalently, D1Z1=D2Z2\frac{D_{1}}{Z_{1}}=\frac{D_{2}}{Z_{2}}, since r=1/2kr=1/2k, and the smallest positive slope is S=1/kS=1/k. The latter equality implies

We have assumed that the three skewed points are intersected by a common line, therefore, by equating the left-hand sides of Equation 6 and Equation 7, after simplification, we get W1=W2W_{1}=W_{2}, or equivalently, Y1=Y2Y_{1}=Y_{2}. Since also X1=X2X_{1}=X_{2}, the points (X1,Y1)(X_{1},Y_{1}), (X2,Y2)(X_{2},Y_{2}) are identical, a contradiction. ∎

Recall now that 44 and Lemma 45 were proven for a general (original) point-grid {0,…,k}×{0,…,k}\{0,\dots,k\}\times\{0,\dots,k\}, for any k≥1k\geq 1. To use them in the context of our reduction, we need to set the proper values for kk and then re-scale the entire grid so that it is inside 2^{2}. From 43, we need k=⌈1t′⌉=⌈30N2+2c2Nα2⌉k=\left\lceil\frac{1}{t^{\prime}}\right\rceil=\left\lceil\frac{30N^{2+2c}2^{N}}{\alpha^{2}}\right\rceil. Now we have a very large grid with edge length 1, so in order to be in 2^{2} as required, consider the grid’s edge length to be t′=α230N2+2c2Nt^{\prime}=\frac{\alpha^{2}}{30N^{2+2c}2^{N}}. Finally, to ensure that all our auxiliary results go through, we have to also scale down rr by further multiplying it with t′t^{\prime}.

Consider the input of an ε\varepsilon-SC-Pizza-Sharing or an ε\varepsilon-Straight-Pizza-Sharing instance, meaning, the description of q∈{n,2n}q\in\{n,2n\} sets, respectively, of weighted polygons with holes on 2^{2}. By definition, its size is N∈Ω(n)N\in\Omega(n). In time poly(N)\textup{poly}(N), we perform a triangulation of each polygon into non-obtuse triangles, therefore, in total we require poly(N)\textup{poly}(N) time for this task. Then, we perform the “pixelation” procedure, which requires, for each of qq mass distributions, checking whether each of ⌈225N2+2c/α2⌉\left\lceil 225N^{2+2c}/\alpha^{2}\right\rceil pixels has a non-empty intersection with a triangle. This task can be performed in poly(N,1/α)\textup{poly}(N,1/\alpha) time, since c>0c>0 is a fixed constant. Next, the points of each point set PiP_{i} created from the respective pixels are shifted positively by an exponentially small value in their yy- and xx-coordinate, so that there is no overlap between any two points. This takes again poly(N,1/α)\textup{poly}(N,1/\alpha) time. Finally, these points are skewed so that they are in general position, by adding very small values in their xx- and yy-coordinates, which takes once again poly(N,1/α)\textup{poly}(N,1/\alpha) time. Also, we remark that the skewing procedure, results in points such that no two of them have the same xx- or yy-coordinate, thus Output (a) of Definition 11 cannot be produced. Finally, notice that, even though we have ensured that each of the created points’ position before skewing is on a vertex of an exponentially large grid, the number of points we have is poly(N,1/α)\textup{poly}(N,1/\alpha).

Recall that Ti,j∈FiT^{i,j}\in\mathcal{F}_{i} is a non-obtuse triangle which belongs to the jj-th polygon of μi\mu_{i}, and Fi\mathcal{F}_{i} is the set of such triangles that compose μi\mu_{i}. By Tpi,jT^{i,j}_{p} we denote the pixelated version of Ti,jT^{i,j}, while M+i,j,M−i,jM^{i,j}_{+},M^{i,j}_{-} are Tpi,jT^{i,j}_{p}’s respective parts of the pixels intersected by lines, that join the R+R^{+} and the R−R^{-} sides, respectively. To bound ∣Pi∣|P_{i}| we will use the fact that ∑Ti,j∈Fiwi,j⋅area(Ti,j)=1\sum_{T^{i,j}\in\mathcal{F}_{i}}w_{i,j}\cdot\text{area}(T^{i,j})=1, or equivalently, ∑Ti,j∈Fiwi,j⋅Nc⋅area(Ti,j)=Nc\sum_{T^{i,j}\in\mathcal{F}_{i}}w_{i,j}\cdot N^{c}\cdot\text{area}(T^{i,j})=N^{c}. We have

and since ∑Ti,j∈Fiarea(Ti,j)≤1\sum_{T^{i,j}\in\mathcal{F}_{i}}\text{area}(T^{i,j})\leq 1, we get

Now observe that, after pixelation, only the pixels at the boundary of each triangle Ti,jT^{i,j} can correspond to ⌈wmax⁡i⋅Nc⌉\left\lceil w^{i}_{\max}\cdot N^{c}\right\rceil points instead of ⌈wi,j⋅Nc⌉\left\lceil w_{i,j}\cdot N^{c}\right\rceil. Therefore, using the notation of Appendix C, where S:=area(Ti,j)S:=\text{area}(T^{i,j}), only at most a fraction S′/SS^{\prime}/S can correspond to ⌈wmax⁡i⋅Nc⌉\left\lceil w^{i}_{\max}\cdot N^{c}\right\rceil points. So,

where the second to last inequality comes from the fact that 1≤Nc/α1\leq N^{c}/\alpha.

where the first inequality is acquired by the reverse triangle inequality in combination with Lemma 42, and the second inequality is due to the fact that the number of points in a pixel of Tpi,jT^{i,j}_{p} is ⌈wmax⁡i⋅Nc⌉\left\lceil w^{i}_{\max}\cdot N^{c}\right\rceil, where wmax⁡i≥wi,jw^{i}_{\max}\geq w_{i,j}, by definition. ∎

Appendix D Proof of Theorem 19

Before we prove the theorem, let us give a brief sketch of the ETR-hardness proof of the aforementioned Consensus-Halving decision version of [DFMS21]. We are given an instance of the following problem which was shown to be ETR-complete (Lemma 15 of the aforementioned paper). Definition 46 (\textscFeasible\textsc{Feasible}_{}). Let p(x1,…,xm)p(x_{1},\dots,x_{m}) be a polynomial. We ask whether there exists a point (x1,…,xm)∈m(x_{1},\dots,x_{m})\in^{m} that satisfies p(x1,…,xm)=0p(x_{1},\dots,x_{m})=0.

We will use exactly the same technique up to the point where we have a Consensus-Halving instance that checks whether q1=q2q_{1}=q_{2}. Then, we use the gadgets described in the FIXP-hardness reduction of Section 3.4 that reduce the valuation functions of nn agents in Consensus-Halving into mass distributions of a SC-Pizza-Sharing instance with nn colours. According to 17 in any solution of the resulting SC-Pizza-Sharing instance, the SC-path does not have turns inside any unit-square. This means that each horizontal/vertical segment of the SC-path that cuts a unit square in SC-Pizza-Sharing has a 1-1 correspondence to a cut of a Consensus-Halving solution, thus a Consensus-Halving solution that uses n−1n-1 cuts would correspond to a SC-path with n−1n-1 line segments, i.e., n−2n-2 turns. Therefore, if and only if there is a SC-path that solves SC-Pizza-Sharing with nn colours and n−2n-2 turns, there is a (n−1)(n-1)-cut that solves Consensus-Halving with nn agents. Equivalently, there is a \vvx∈m\vv{x}\in^{m} such that p(\vvx)=0p(\vv{x})=0, making the \textscFeasible\textsc{Feasible}_{} instance satisfiable. ∎

Appendix E Proof of Theorem 22

We also define the variables z1,z2,…,zm+1z_{1},z_{2},\dots,z_{m+1} where ∣zi∣=yi−yi−1|z_{i}|=y_{i}-y_{i-1}, i∈[m+1]i\in[m+1]. The sign of ziz_{i}, i∈[m+1]i\in[m+1] indicates the sign of the leftmost part of slice [yi−1,yi][y_{i-1},y_{i}] that xi−1x_{i-1} defines, and we set the whole slice [0,y1][0,y_{1}] to have the sign of z1z_{1}. Clearly, (z1,z2,…,zm+1)∈Sm(z_{1},z_{2},\dots,z_{m+1})\in S^{m}, since ∑i=1m+1∣zi∣=1\sum_{i=1}^{m+1}|z_{i}|=1, and (x1,…,xm)∈m(x_{1},\dots,x_{m})\in^{m}, by definition. A feasible solution of SC-Pizza-Sharing is then the vector (z1,…,zm+1,x1,…,xm)(z_{1},\dots,z_{m+1},x_{1},\dots,x_{m}); this defines a directed path that consists of horizontal and vertical line segments with at most 2m−12m-1 turns. We can recover this path using Algorithm 1. Note that in the algorithm we implicitly suggest to navigate using the following trick: if we are at a horizontal part of the path, say y=y′y=y^{\prime}, and we hit the boundary x=0x=0 or x=1x=1 before we reach a vertical cut, then we wrap around in the horizontal dimension and continue from point (1,y′)(1,y^{\prime}) or (0,y′)(0,y^{\prime}) respectively.

Finally, we highlight that in the map-back process, RR is only used to indicate the sign of variable zm+1z_{m+1}, while its absolute value has no other purpose than to serve as a remainder: it ensures that (Z1,…,Zm,R,X1,…,Xm)∈S2m(Z_{1},\dots,Z_{m},R,X_{1},\dots,X_{m})\in S^{2m}, i.e., it is ∣R∣=2m−(∑i=1m∣Zi∣+∑i=1m∣Xi∣)|R|=2m-(\sum_{i=1}^{m}|Z_{i}|+\sum_{i=1}^{m}|X_{i}|).

For any given point \vvP=(Z1,…,Zm,R,X1,…,Xm)∈S2m\vv{P}=(Z_{1},\dots,Z_{m},R,X_{1},\dots,X_{m})\in S^{2m}, the Borsuk-Ulam function is defined to be the total “++” (positive) measure on 2^{2} induced by \vvP\vv{P}, and we denote it by μ(R+;\vvP)\mu(R^{+};\vv{P}), that is, f(\vvP)=μ(R+;\vvP)f(\vv{P})=\mu(R^{+};\vv{P}). The total positive measure is a continuous function of the variables: ff is a continuous function of (z1,…,zm+1,x1,…,xm)(z_{1},\dots,z_{m+1},x_{1},\dots,x_{m}), by the interpretation of the variables on 2^{2}; and the mapping we defined from \vvP\vv{P} to (z1,…,zm+1,x1,…,xm)(z_{1},\dots,z_{m+1},x_{1},\dots,x_{m}) is also continuous. To see the latter, note that for fixed signs of the variables of \vvP\vv{P}, the variables (z1,…,zm+1,x1,…,xm)(z_{1},\dots,z_{m+1},x_{1},\dots,x_{m}) are continuous by definition, therefore we only need to check what happens when a variable from the latter changes sign. By the mapping we defined above, one can verify that whenever a variable changes its sign it has to necessarily pass from value 0 first, and hence the mapping is continuous.

The proof is similar when n=2m−1n=2m-1. The horizontal cuts are again 0≤y1≤y2≤⋯≤ym≤10\leq y_{1}\leq y_{2}\leq\dots\leq y_{m}\leq 1, but the vertical cuts are x1,x2,…,xm−1∈x_{1},x_{2},\dots,x_{m-1}\in, meaning that the top slice is not vertically cut. Also, it is easy to see that one could consider the path to be again yy-monotone but in the opposite direction, meaning that there is no line segment pointing upwards.

Note that a symmetric proof exists, where the slices are vertical instead of horizontal, and the cuts within the slices are horizontal instead of vertical. The analysis is similar to the one we give here, and it guarantees the existence of an SC-path which is allowed to wrap around in the vertical dimension, it bisects all nn measures and is xx-monotone with no line segment heading left (or right).

Appendix F Proof of Theorem 24

Consider an arbitrary instance of exact SC-Pizza-Sharing, i.e., the one of Definition 3, where we additionally restrict the mass distributions to be weighted polygons (with holes). For ease of presentation, we will give distinct colours to the mass distributions μi\mu_{i}, i∈[n]i\in[n] and call them colours 1,2,…,n1,2,\dots,n. We are given a set of polygons with holes in the input form described in Appendix A. We will first do a preprocessing of the input: (a) normalization in 2^{2}, and (b) triangulation. The former is needed in order to simulate the space of the mathematical existence proof of Theorem 22, while the latter allows us to construct the Borsuk-Ulam function of the aforementioned proof via circuits. For the computation of the Borsuk-Ulam function described in Appendix E we need a way of computing the “positive” measure of a polygon, as dictated by a given feasible solution \vvP\vv{P}. To simplify this computation, we will triangulate further the polygon to end up with only non-obtuse triangles, which can be decomposed into axis-aligned right-angled triangles (Section F.1).

For task (a) we first check among all polygons of all colours, what the smallest and greatest coordinates of their vertices are, which defines a rectangle that includes all our polygons. Then, we find the square with smallest sides within which the aforementioned rectangle can be inscribed. Finally, by shifting the rectangle so that its bottom left vertex is at (0,0)(0,0), and then scaling it down so that each side is of length 11, we get a normalized input where each polygon is in 2^{2} and their relative positions and areas are preserved. Task (b) is easily taken care of via already known algorithms (e.g., [GJPT78, AAP86, Meh84, GM91]) that triangulate polygons with holes in O(nlog⁡n)O(n\log n) time (which is optimal) without inserting additional vertices.

Here we first show how an arbitrary triangle can be decomposed into right-angled triangles whose right angle is additionally axis-aligned. Then, by using only the allowed gates of BU, we present a way to construct a Borsuk-Ulam function. We show that any solution of Borsuk-Ulam whose input is the aforementioned function can be mapped back in polynomial time to an exact SC-Pizza-Sharing solution.

By the definition of the Borsuk-Ulam function of Section 4, it is apparent that we need to be able to compute parts of the area of a polygon, depending on where square-cuts fall. The function we provide can compute parts of the area of axis-aligned right-angled triangles. For that reason, after a standard triangulation (e.g., using the technique of [AAP86]) we further preprocess it and decompose each obtuse triangle into two right-angled triangles, by adding an extra line segment.

In particular, we check the obtuseness of a triangle \trnglABC\trngl{ABC} by computing the squared lengths of its sides AB2,BC2,AC2AB^{2},BC^{2},AC^{2} (each is rational; a sum of squares of rationals), taking the largest one, w.l.o.g. AC2AC^{2} and then checking whether AB2+BC2<AC2AB^{2}+BC^{2}<AC^{2}. If the inequality is not true then \trnglABC\trngl{ABC} is non-obtuse and we proceed. Otherwise, we add the line segment BDBD that starts from BB and ends at DD on side ACAC, where BDC^=ADB^=90∘\widehat{BDC}=\widehat{ADB}=90^{\circ}. We will denote by \trnglABC\trngl{ABC} a triangle with vertices A,B,CA,B,C and, when clear from context, we will also use the same notation to indicate the area of the triangle. Two intersecting line segments ABAB, BCBC define two angles, denoted ABC^\widehat{ABC} and CBA^\widehat{CBA}. The order of the vertices implies a direction of the segments, i.e., in the former angle we have ABAB, BCBC and in the latter we have CBCB, BABA. We consider the direction of the segments and define the angle to be the intersection of the left halfspaces of the segments. Therefore ABC^=360∘−CBA^\widehat{ABC}=360^{\circ}-\widehat{CBA}. This order will not matter if clear from context (e.g., in triangles). The coordinates (xD,yD)(x_{D},y_{D}) of DD are rationals since they are the solution of the following two equations: (a) one that dictates that DD is on ACAC: yD−yAxD−xA=yA−yCxA−xC\frac{y_{D}-y_{A}}{x_{D}-x_{A}}=\frac{y_{A}-y_{C}}{x_{A}-x_{C}}, and (b) one that captures the fact that BDBD and ACAC are perpendicular: yB−yDxB−xD⋅yA−yCxA−xC=−1\frac{y_{B}-y_{D}}{x_{B}-x_{D}}\cdot\frac{y_{A}-y_{C}}{x_{A}-x_{C}}=-1. In fact,

At this point the triangulation of each polygon consists of non-obtuse triangles. The following proposition makes the computation of areas of these triangles’ parts possible via a decomposition into axis-aligned right-angled triangles.

The area of any non-obtuse triangle can be computed by the areas of five axis-aligned right-angled triangles.

A proof by picture is presented in Figure 8(b), where we draw a segment from the top-left corner to the bottom-right one, and \trnglXYB+\trnglXBZ−\trnglAYB−\trnglXAC−\trnglCBZ\trngl{XYB}+\trngl{XBZ}-\trngl{AYB}-\trngl{XAC}-\trngl{CBZ}.

The proof is immediate if we show that every non-obtuse triangle \trnglABC\trngl{ABC} can be tightly inscribed inside a rectangle, meaning that all of its vertices touch the rectangle’s perimeter. In particular, the rectangle has coordinates

First, we argue that one of \trnglABC\trngl{ABC}’s vertices is on a corner of the rectangle. W.l.o.g. suppose arg⁡min⁡{xA,xB,xC}=A\arg\min\{x_{A},x_{B},x_{C}\}=A. If arg⁡min⁡{yA,yB,yC})=A\arg\min\{y_{A},y_{B},y_{C}\})=A or arg⁡max⁡{yA,yB,yC})=A\arg\max\{y_{A},y_{B},y_{C}\})=A then AA is the bottom-left or top-left corner of the rectangle, respectively. Otherwise, suppose w.l.o.g. that arg⁡min⁡{yA,yB,yC})=B\arg\min\{y_{A},y_{B},y_{C}\})=B. If arg⁡max⁡{xA,xB,xC}=B\arg\max\{x_{A},x_{B},x_{C}\}=B then BB is the bottom-right corner of the rectangle. Otherwise, arg⁡max⁡{xA,xB,xC}=C\arg\max\{x_{A},x_{B},x_{C}\}=C. Then, if arg⁡max⁡{yA,yB,yC})=C\arg\max\{y_{A},y_{B},y_{C}\})=C, CC is the top-right corner of the rectangle. Otherwise, arg⁡max⁡{yA,yB,yC})=A\arg\max\{y_{A},y_{B},y_{C}\})=A and AA is the top-left corner of the rectangle. We conclude that a vertex of the triangle, w.l.o.g. AA, is on a corner of the rectangle.

Now we need to show that BB and CC lie on the perimeter of the rectangle. By the definition of the aforementioned rectangle’s vertices, it cannot be that both CC and DD are not on the perimeter; if, for example, AA is the top-right corner then arg⁡min⁡{yA,yB,yC}∈{B,C}\arg\min\{y_{A},y_{B},y_{C}\}\in\{B,C\}, therefore vertex arg⁡min⁡{yA,yB,yC}\arg\min\{y_{A},y_{B},y_{C}\} touches the lower side of the rectangle, and similarly if AA is any of the other corners. Now, for the sake of contradiction, suppose that the remaining vertex, w.l.o.g. CC does not touch the perimeter of the rectangle. Then, by definition of the rectangle’s vertices, BB is on another corner. If AA and BB are on the same side of the rectangle, then it is clear by the definition of the rectangle’s coordinates that CC must be on the perimeter of it. Otherwise AA and BB are diagonal corners of the rectangle. Then, ABAB is the largest side of the triangle and if CC is not on the boundary, it holds that AC2+BC2<AB2AC^{2}+BC^{2}<AB^{2}, meaning that \trnglABC\trngl{ABC} is obtuse, a contradiction. Therefore, any non-obtuse triangle can be tightly inscribed in a rectangle. ∎

F.2 Constructing the Borsuk-Ulam function

Consider the τ≥1\tau\geq 1 weighted polygons of the ii-th colour, and let us focus on a particular polygon t∈[τ]t\in[\tau]. We have triangulated the polygon into mtm_{t} non-obtuse triangles. Consider one such triangle j∈[mt]j\in[m_{t}] and the virtual triangles Tj1,Tj2,Tj3,Tj4,Tj5T_{j}^{1},T_{j}^{2},T_{j}^{3},T_{j}^{4},T_{j}^{5}, which are the corresponding five axis-aligned right-angled triangles described in the previous subsection (also see Figure 8(b)). W.l.o.g. we consider Tj1T_{j}^{1} and Tj2T_{j}^{2} to be the positively contributing triangles and the rest to be the negatively contributing triangles. For each of them we will be computing the positive measure that the cut \vvp\vv{p} of SC-Pizza-Sharing defines (see Appendix E). By the axis-aligned right-angled triangle decomposition described in the proof of Proposition 48, it suffices to show how to compute parts of areas of such a triangle, for all of its four possible orientations: QI,QII,QIII,QIVQ_{I},Q_{II},Q_{III},Q_{IV}, where QoQ_{o} is the orientation when, by shifting the triangle so that the vertex of the right angle is on (0,0)(0,0), the whole triangle is in the oo-th quadrant.

First, we test the orientation of our triangle. Observe that for Tj1T_{j}^{1} and Tj2T_{j}^{2} we know that the orientation is QIQ_{I} and QIIIQ_{III}, respectively (see \trnglXYB\trngl{XYB} and \trnglXZB\trngl{XZB} in Figure 8(b)). For Tj3,Tj4,Tj5T_{j}^{3},T_{j}^{4},T_{j}^{5}, we check the coordinates and identify what kind of triangles they are.

For a fixed colour i∈[n]i\in[n], for each possible category QoQ_{o}, o∈{I,II,III,IV}o\in\{I,II,III,IV\} we show how to compute the term that an axis-aligned right-angled triangle Tjr=\trnglABCT_{j}^{r}=\trngl{ABC} (as in Figure 8(b)), j∈[m]j\in[m], r∈r\in contributes to the Borsuk-Ulam function f(\vvp)if(\vv{p})_{i}. Let us denote by gj,r(x)g_{j,r}(x) the linear function of the line segment of \trnglABC\trngl{ABC}’s hypotenuse (with respect to the horizontal axis xx), and note that it is easy to derive from the segment’s end-points. We also denote by (xj,r,yj,r)(x_{j,r},y_{j,r}) and (xj,r′,yj,r′)(x^{\prime}_{j,r},y^{\prime}_{j,r}) the bottom and top points that define the hypotenuse of \trnglABC\trngl{ABC}. The vertex of the right angle is defined by these four coordinates depending on QoQ_{o}.

For any given point (feasible solution) \vvp\vv{p} that defines a SC-path of SC-Pizza-Sharing (see Appendix E), recall that there are ⌈n/2⌉\left\lceil n/2\right\rceil slices 0=y0≤y1≤⋯≤y⌈n/2⌉≤y⌈n/2⌉+1=10=y_{0}\leq y_{1}\leq\dots\leq y_{\left\lceil n/2\right\rceil}\leq y_{\left\lceil n/2\right\rceil+1}=1 that determine ⌈n/2⌉+1\left\lceil n/2\right\rceil+1 slices ∣z1∣,…,∣z⌈n/2⌉+1∣|z_{1}|,\dots,|z_{\left\lceil n/2\right\rceil+1}| that partition 2^{2}. One can see that, given the zsz_{s}’s, s∈[⌈n/2+1⌉]s\in[\left\lceil{n/2}+1\right\rceil] from \vvp\vv{p}, the ysy_{s}’s can be computed by the recursive expression ys=∣zs∣−ys−1y_{s}=|z_{s}|-y_{s-1}.

We first define the following auxiliary functions Ajr(zs)A_{j}^{r}(z_{s}), Bjr(−zs)B_{j}^{r}(-z_{s}) for each QoQ_{o} and r∈r\in:

For an example of case QIIQ_{II} see Figure 9.

Then, we get the total area that all slices up to zsz_{s} (i.e., [0,ys][0,y_{s}]) define together with xsx_{s} from

which has the property that if zs>0z_{s}>0 (resp. zs<0z_{s}<0), then the computed area considers the thickness ∣zs∣|z_{s}| only in the part of the slice that is left (resp. right) to xsx_{s}. For zs=0z_{s}=0, ∣zs∣|z_{s}| has no thickness, its measure is 0, and consistently vanishes from the above functions.

Then, to isolate the part that only slice zsz_{s} contributes to the positive measure, for r∈r\in (see definition of virtual triangles TjrT_{j}^{r} above), we define

Consequently, the positive measure that an (unweighted) non-obtuse triangle jj contributes to the Borsuk-Ulam function according to the SC-path \vvp\vv{p} is

Finally, given that colour i∈[n]i\in[n] has τ\tau many weighted polygons, each of weight wtw_{t}, t∈[τ]t\in[\tau], which has been decomposed into mtm_{t} many non-obtuse triangles, ii’s positive measure (i.e., the ii-th coordinate of the Borsuk-Ulam function) is

Appendix G Proof of Theorem 25

The proof is by showing that the aforementioned problem can be formulated as an ETR problem. We will use the machinery of the BU containment proof (Theorem 24) to build the piece-wise polynomial function ff. In the proof of the aforementioned theorem we construct the function f:Sn↦Rnf:S^{n}\mapsto R^{n} so that it computes the “positive” part of measure i∈[n]i\in[n] in its ii-th coordinate, and we require (at most) n−1n-1 turns in the SC-path. It is easy to modify this construction for any given number kk of turns: by the proof of Theorem 22, what we need to do is to cut ⌈(k+1)/2⌉+1\left\lceil(k+1)/2\right\rceil+1 slices (i.e., place ⌈(k+1)/2⌉\left\lceil(k+1)/2\right\rceil horizontal cuts) in 2^{2} and another k+1−⌈(k+1)/2⌉k+1-\left\lceil(k+1)/2\right\rceil vertical cuts on them starting with the second slice from the bottom. Then, by the same mapping as in the aforementioned proof, we have a function f:Sk+1↦Rnf:S^{k+1}\mapsto R^{n} for which, if f(\vvP∗)=f(−\vvP∗)f(\vv{P}^{*})=f(-\vv{P}^{*}) for some \vvP∗\vv{P}^{*} we have a solution to exact SC-Pizza-Sharing with kk turns in its SC-path. One can see that the SC-Pizza-Sharing problems of Theorem 22 and Theorem 24 were special cases for k≥n−1k\geq n-1, for which the aforementioned equality is guaranteed by the Borsuk-Ulam theorem.

Appendix H Proof of Theorem 27

Appendix I Proof of Theorem 28

As shown in the proof of Theorem 27, the function ff we use (proof of Theorem 24) is continuous, piece-wise polynomial, and therefore λ\lambda-Lipschitz continuous, where λ=max⁡j=1n+1{sup⁡x∣∣∂f(x)∂xj∣∣∞}\lambda=\max_{j=1}^{n+1}\left\{\sup_{x}\left|\left|\frac{\partial f(x)}{\partial x_{j}}\right|\right|_{\infty}\right\} (again, note that λ\lambda is constant, and points where fi(x)f_{i}(x), i∈[n]i\in[n] is non-differentiable do not matter for Lipschitzness). Suppose that a point x∈Sk+1x\in S^{k+1} satisfies ∣∣f(x)−f(−x)∣∣∞≤ε||f(x)-f(-x)||_{\infty}\leq\varepsilon for a given ε>0\varepsilon>0, and therefore it is a solution to the decision version of ε\textsc−Borsuk−Ulam\varepsilon\textsc{-Borsuk-Ulam} when parameterized by the number of turns kk. Then, the following inequalities hold for any other point y∈Sk+1y\in S^{k+1}:

Consider now a number MM that is upper-bounded by an inverse-polynomial of the input size. Then, for any yy such that ∣∣x−y∣∣∞≤M2λ||x-y||_{\infty}\leq\frac{M}{2\lambda} we have from Equation 9 that

meaning that yy is also a solution to ε′\varepsilon^{\prime}-Borsuk-Ulam parameterized by kk, for ε′=M+ε\varepsilon^{\prime}=M+\varepsilon. Therefore, yy is a solution to an instance of the same problem when additively relaxed by at most an inverse polynomial (in the input size) quantity.

By the above, we conclude that there are always polynomial-size solutions to the problem (and thus to the initial ε\varepsilon-SC-Pizza-Sharing problem) given that there are exact solutions to it. Therefore, a candidate solution can be verified in polynomial time.

Acknowledgements

The second author was supported by EPSRC grant EP/W014750/1 “New Techniques for Resolving Boundary Problems in Total Search”. The third author was supported by the Alexander von Humboldt Foundation with funds from the German Federal Ministry of Education and Research (BMBF).

References