Tropical bisectors and Voronoi diagrams

TL;DR

This paper introduces the concept of tropical Voronoi diagrams, exploring their structure and algorithms by studying bisectors under polyhedral norms, with a focus on tropical geometry.

Contribution

It initiates the study of tropical Voronoi diagrams, providing structural insights and algorithms based on bisectors of points with respect to polyhedral norms.

Findings

Bisectors of three points are homeomorphic to open subsets of Euclidean space under certain conditions.

Structural properties of tropical bisectors are established.

Algorithms for constructing tropical Voronoi diagrams are proposed.

Abstract

In this paper we initiate the study of tropical Voronoi diagrams. We start out with investigating bisectors of finitely many points with respect to arbitrary polyhedral norms. For this more general scenario we show that bisectors of three points are homeomorphic to a non-empty open subset of Euclidean space, provided that certain degenerate cases are excluded. Specializing our results to tropical bisectors then yields structural results and algorithms for tropical Voronoi diagrams.

Equations

$$\mathrm{dist}(a,b)=\underset{i\in [d+1]}{\max }\left(a_i-b_i\right)-\underset{j\in [d+1]}{\min }\left(a_j-b_j\right)=\underset{i,j\in [d+1]}{\max }(a_i-b_i-a_j+b_j).$$

$$\mathrm{dist}(a,b)=\underset{i\in [d+1]}{\max }\left(a_i-b_i\right).$$

$$\begin{array}{rl}{\mathbb{B}}^d& =\left\{x\in {\mathbb{R}}^{d+1}/\mathbb{R}1\right|\mathrm{dist}(x,0)\le 1\}=\bigcap _{i\ne j}\left\{x\in {\mathbb{R}}^{d+1}/\mathbb{R}1\right|x_i-x_j\le 1\}\\ & =\frac{1}{2}\mathrm{conv}\left(\{\pm 1{\}}^{d+1}\setminus \{\pm 1\}\right)+\mathbb{R}1.\end{array}$$

$$F=\sum _{i=1}^{d+1}{s}_i,$$

$$F_{-}=$$

$$F_{\ast }=$$

$$F_{+}=$$

$$\mathrm{bis}(S):=\left\{x\in {\mathbb{R}}^d\right|\mathrm{dist}(a,x)=\mathrm{dist}(b,x)\text{ for }a,b\in S\}.$$

$${\mathrm{bis}}_{(F_1,\dots ,F_k)}(\{a_1,\dots ,a_k\})=\mathrm{bis}(\{a_1,\dots ,a_k\})\cap a_1+F_1\cap \cdots \cap a_k+F_k,$$

$${\mathrm{bis}}_{(F_1,\dots ,F_k)}(a_1,\dots ,a_k)$$

$${\mathrm{bis}}_{(F_1,\dots ,F_k)}(a_1,\dots ,a_k),\text{ with }{F}_1,\dots ,F_k\in \mathcal{F}(K),$$

$${\mathrm{bis}}_{(F_1,\dots ,F_k)}(S)\cap {\mathrm{bis}}_{(F_1^{\prime },\dots ,F_k^{\prime })}(S)={\mathrm{bis}}_{(F_1\cap F_1^{\prime },\dots ,F_k\cap F_k^{\prime })}(S).$$

$$v_{1}=\left(-\sin\frac{2\pi}{3},\,\cos\frac{2\pi}{3}\right)\,,\ v_{2}=\left(\sin\frac{2\pi}{3},\,\cos\frac{2\pi}{3}\right)\,,\ v_{3}=\bigl{(}0,1\bigr{)}\enspace.$$

$${\mathrm{bis}}_{(F_1,\dots ,F_k)}(a_1,\dots ,a_k)=\varnothing ⟺{\mathrm{bis}}_{(F_1,\dots ,F_k)}(a_1^{\prime },\dots ,a_k^{\prime })=\varnothing .$$

$${\mathrm{bis}}_{(F,F)}(a,b)=(a+F)\cap (b+F),$$

$${\mathrm{bis}}_{(F_1,\dots ,F_k)}(a_1,\dots ,a_k)=Q\cap H$$

$${\mathrm{bis}}_{(F_1,\dots ,F_{d+1})}(a_1^{\prime },\dots ,a_{d+1}^{\prime })=Q^{\prime }\cap H.$$

$$({\lambda }_{F_i}-{\lambda }_{F_1})(x)={\lambda }_{F_i}(a_i)-{\lambda }_{F_1}(a_1)\text{for }i=2,\dots ,k.$$

$${\lambda }_{F_1}(x-a_1^{\prime })=\cdots ={\lambda }_{F_k}(x-a_k^{\prime }),$$

$${\mathrm{bis}}_{(F_1,\dots ,F_k)}(a_1^{\prime },\dots ,a_k^{\prime })=Q^{\prime }\cap H^{\prime }.$$

$$a_1-x=(0,2,2,4),$$

$$a_2-x=(0,4,1,3),$$

$$a_3-x=(2,4,0,2),$$

$$a_4-x=(4,0,1,3),$$

$$H_S(a)=\bigcap _{b\in S\setminus \{a\}}H(b-a).$$

$$H_S(a)=\mathrm{relint}\left(\bigcup \left\{F_i\right|{\lambda }_{F_i}(a)<{\lambda }_{F_i}(b)\text{ for all }b\in S\backslash a\}\right).$$

$$\bigcup _{a\in S}\overline{H_S(a)}=\partial K.$$

$${\varphi }^{-1}(C)\subset (a+\mathrm{cone}(K))\cap (b+\mathrm{cone}(\overline{H(a-b)})),$$

$$\mathrm{dist}(a+0v,a)=0,\mathrm{dist}(a+0v,b)>0,$$

$$\underset{\alpha \to \infty }{\lim }\mathrm{dist}(a+\alpha v,b)<\underset{\alpha \to \infty }{\lim }\mathrm{dist}(a+\alpha v,a).$$

$$\begin{array}{rl}\mathrm{bis}(a,b)=& \left\{x\right|x_3+1\le x_1,x_2,\le x_4-1\}\\ & \cup \left\{x\right|x_4+1\le x_1,x_2,\le x_3-1\}\cup \left\{x\right|x_1=x_2\}.\end{array}$$

$$\begin{array}{rl}\mathrm{bis}(c,d)=& \left\{x\right|x_1+2\le x_3,x_4,\le x_2-2\}\\ & \cup \left\{x\right|x_2+2\le x_3,x_4,\le x_1-2\}\cup \left\{x\right|x_3=x_4\}.\end{array}$$

$$\mathrm{bis}(a,b,c,d)\subseteq \left\{x\right|x_1=x_2,x_3=x_4\}.$$

Full text

F. Criado has been supported by Berlin Mathematical School and Einstein Foundation Berlin (EVF-2015-230). M. Joswig has been supported by Deutsche Forschungsgemeinschaft (EXC 2046: \csq@thequote@oinit\csq@thequote@oopenMATH+\csq@thequote@oclose, SFB-TRR 109: \csq@thequote@oinit\csq@thequote@oopenDiscretization in Geometry and Dynamics\csq@thequote@oclose, SFB-TRR 195: \csq@thequote@oinit\csq@thequote@oopenSymbolic Tools in Mathematics and their Application\csq@thequote@oclose, and GRK 2434: \csq@thequote@oinit\csq@thequote@oopenFacets of Complexity\csq@thequote@oclose). F. Santos has been supported by the Einstein Foundation Berlin (EVF-2015-230) and by Grants MTM2017-83750-P/AEI/10.13039/ 501100011033 and PID2019-106188GB-I00/AEI/10.13039/501100011033 of the Spanish State Research Agency

1. Introduction

One early route to the success of tropical geometry is based on the tropicalization of classical algebraic varieties defined over some valued field. Key examples include Mikhalkin’s correspondence principle, which relates tropical plane curves with classical complex algebraic curves [24], or the tropical Grassmannians of Speyer and Sturmfels [25]. In all of this the focus lies on the combinatorial properties of tropical varieties, which are ordinary polyhedral complexes.

More recently, however, tropical semi-algebraic sets and their intrinsic geometry came into the picture; see [1], [15]. For instance, their metric properties appear in [2] as a tool to show that standard versions of the interior point method of linear programming exhibit an exponential complexity in the unit cost model. The proof of that result is based on translating metric data on a family of tropical linear programs into curvature information about the central paths of their associated ordinary linear programs. Similarly, tropical analogs of isoperimetric (or isodiametric) inequalities have been studied in [8], where a tropical volume is defined that corresponds to an “energy gap” in mathematical physics [19]. Another example is the statistical analysis of phylogenetic trees by Lin, Monod and Yoshida [21].

We feel that all this calls for a more systematic investigation of metric properties of tropical varieties. Starting from first principles, this naturally leads to tropical Voronoi diagrams. The tropical distance between two points $a,b\in\mathbb{R}^{d+1}$ is

[TABLE]

It does not depend on choosing $\min$ or $\max$ as the tropical addition. The map $\operatorname{dist}:\mathbb{R}^{d+1}\times\mathbb{R}^{d+1}\to\mathbb{R}$ is non-negative, symmetric, and it satisfies the triangle inequality. Moreover, it is homogeneous, so it induces a norm on the tropical $d$-torus $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}\cong\mathbb{R}^{d}$, where $\mathbb{1}=(1,\dots,1)$ denotes the all ones vector. The tropical Voronoi region of a site $s\in S$ with respect to a set $S$ comprises those points in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ to which $s$ is the nearest among all sites in $S$, with respect to $\operatorname{dist}$. The tropical Voronoi diagram $\operatorname{Vor}(S)$ is the cell decomposition of $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ into Voronoi regions. Tropical Voronoi diagrams are a special case of Voronoi diagrams for polyhedral norms, a classical topic in convexity and computational geometry; cf. [4, Sect. 7.2] or [23, Sect. 4].

The intersection of two or more Voronoi regions is part of a bisector, i.e., the locus of points which are equidistant to a given set. For instance, in the Euclidean case the bisector of two points is a degenerate quadric which agrees with an affine hyperplane as a set. In the tropical setting, the bisector of two points can also be described as part of a tropical hypersurface, but this is now of degree $d+1$; see Proposition 4. Further, in the tropical case two points may already produce degenerate bisectors (which may contain, e.g., full-dimensional pieces), whereas the first degenerate case in the Euclidean metric arises for three points. So tropical Voronoi diagrams behave quite differently from Euclidean Voronoi diagrams.

Yet there are also similarities. A key structural result is that the tropical Voronoi regions are star convex and can be described as unions of finitely many ordinary polyhedra; see Proposition 1 and Theorem 6. We prove a second main result, Theorem 3, for the more general case of an arbitrary polyhedral norm in $\mathbb{R}^{d}$: the bisector of any three points in weak general position is homeomorphic to an open subset of $\mathbb{R}^{d-2}$. Our proof generalizes the arguments from [13], [14], where a similar result was proved for smooth norms in $d=2,3$. However, the global topology of tropical bisectors of three or more points can be radically different from the topology of the classical bisectors. For instance, tropical bisectors are sometimes disconnected and, more strongly, $d+1$ points can have more than one circumcenter. This may happen even in general position; see Examples 1 and 3. We do not know if bisectors may have nontrivial higher Betti numbers, but we suspect they can; see Theorem 4.

Another contribution in our paper is a randomized incremental algorithm for computing the tropical Voronoi diagram of $n$ points in general position in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ with an expected running time of order $O(n^{d}\log n)$, for fixed dimension $d$; see Theorem 10. Euclidean Voronoi diagrams of finite point sets can be explained fully in terms of ordinary convex polyhedra and convex hull algorithms; see [7], [4]. We do not know if there is a tropical analog.

Amini and Manjunath [3] study the Voronoi diagram of a lattice with respect to the following asymmetric version of the tropical distance:

[TABLE]

As they show in [3, Lemma 4.7], this is the polyhedral distance obtained taking as unit ball the standard simplex. Their motivation comes from work of Baker and Norine [5] on a Riemann–Roch theorem for graphs, which implies a Riemann–Roch theorem for tropical curves.

Our paper is organized as follows. The short Section 2, in which we verify that the tropical distance is induced by a polyhedral norm and discuss the combinatorics of the tropical unit ball, sets the stage. In Section 3 we collect our general structural results on bisectors and Voronoi diagrams. The results in this section are proved for general polyhedral norms, but all our examples address the tropical case. A subtle point is the right concept of “general position”. In fact, we distinguish between weak general position which prevents bisectors to contain full-dimensional parts (see Proposition 2), and a stronger general position which is defined via stability of facets in the bisector under perturbation of the sites. For instance, the bisector of any number $k$ of points in general position in $\mathbb{R}^{d}$ is either empty or a polyhedral complex of pure dimension $d+1-k$; see Corollary 2. As a special case, the bisector of $d+1$ points in $\mathbb{R}^{d}$ in general position is finite. Section 4 returns to the tropical case. We specialize our results on bisectors in general polyhedral norms, and we show that the combinatorial types of tropical bisectors of two points are classified in terms of a certain polyhedral fan related to the tropical unit ball and the braid arrangement; see Theorems 5 and 6. This is related to work of Develin [9] on the moduli of tropically collinear points. Finally, in Section 5 we discuss algorithms. This includes a tropical variant of Fortune’s beach line algorithm [10] for planar Voronoi diagrams as well as the aforementioned algorithm in arbitrary dimension.

Acknowledgment.

We thank Günter Rote and three anonymous referees for useful comments on a first version of this paper. Especially, one referee brought the reference [3] to our attention.

2. The tropical unit ball

The unit ball with respect to the tropical distance function defined in (1) is

[TABLE]

In this way, $\mathbb{B}^{d}$ is a polytope in the tropical torus $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$. We also write $\mathbb{B}^{d}(a,r)$ for the tropical ball with center $a$ and radius $r$. All tropical balls result from scaling and translating $\mathbb{B}^{d}$. In fact, the tropical norm agrees with polyhedral norm with respect to the tropical unit ball, in the sense of Section 3. Such distances are called convex distance functions in [4, Sect. 7.2]; see also [12, 14, 13].

Both the inequality and the vertex descriptions of $\mathbb{B}^{d}$ in Eq. (2) are non-redundant:

• •

$\mathbb{B}^{d}$ has $d(d+1)$ facets. Each facet corresponds to a choice of coordinates achieving the maximum and the minimum.

• •

$\mathbb{B}^{d}$ has $2^{d+1}-2$ vertices. Each vertex corresponds to a (nontrivial) partition of the coordinates into maxima and minima. For example, $\mathbb{B}^{2}$ is a hexagon and $\mathbb{B}^{3}$ is a rhombic dodecahedron.

The vertex description also shows that $\mathbb{B}^{d}$ equals the projection of the $(d{+}1)$-dimensional regular cube $[-1,1]^{d+1}$ in $\mathbb{R}^{d+1}$ along the direction $\mathbb{1}$. That is, $\mathbb{B}^{d}$ is a zonotope with $d+1$ generators in general position, and all its faces are parallelepipeds. These generators correspond to the $d+1$ coordinate directions in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$. This suggests a combinatorial way to specify the faces of $\mathbb{B}^{d}$: Each face $F$ can be written as a Minkowski sum

[TABLE]

where each $s_{i}$ is one of $\{-e_{i}\}$, $[-e_{i},e_{i}]$ or $\{+e_{i}\}$. We say that $F$ is of type $(F_{-},F_{*},F_{+})$ if

[TABLE]

Conversely, a partition of $[d+1]$ into three parts $F_{-}$, $F_{*}$, $F_{+}$ corresponds to a face of $\mathbb{B}^{d}$ if and only if neither $F_{-}$ nor $F_{+}$ is empty. Moreover, the dimension of $F$ equals the cardinality of $F_{*}$. In particular, the vertices of $\mathbb{B}^{d}$ correspond to the $2^{d+1}-2$ ways of partitioning $[d+1]$ into two non-empty subsets. The facets of $\mathbb{B}^{d}$ correspond to the $d(d+1)$ ways of choosing an ordered pair from $[d+1]$, without repetition.

Remark 1**.**

The zonotope $\mathbb{B}^{d}$ is dual to an arrangement of $d+1$ linear hyperplanes in general position in $\mathbb{R}^{d}$, oriented so that the intersection of all positive half-spaces is empty. In particular, its face lattice is the same as the lattice of non-zero covectors of the unique totally cyclic oriented matroid of rank $d$ with $d+1$ elements. Covectors of an oriented matroid are usually written as $(V_{-},V_{0},V_{+})$ but in our context we prefer to use $*$ instead of zero meaning that the corresponding coordinate is not fixed.

Remark 2**.**

Another general description of $\mathbb{B}^{d}$ is that it equals the (ordinary) Voronoi cell of the lattice of type $A_{d}$ (i.e., the triangular lattice for $d=2$ and the face centered cubic lattice (FCC) for $d=3$). Similarly, $\mathbb{B}^{d}$ is the polytope polar to the difference body $T-T$ of a regular $d$-simplex $T$. This description shows that $\mathbb{B}^{d}$ is the same as the polytope $U_{d}$ that appears in Makeev’s conjecture. See, e.g., [26, Conjecture 21.3.2].

3. Bisectors in polyhedral norms

Throughout this section we work in the general framework of Minkowski norms; see [4, Sect. 7.2], [12], [23]. Consider a polytope $K\subset\mathbb{R}^{d}$ with the origin in its interior. Let $\operatorname{dist}(a,b)$ be the unique scaling factor $\alpha>0$ such that $b-a\in\partial(\alpha K)$. Then, $\operatorname{dist}$ satisfies the triangle inequality, is invariant under translation, and homogeneous under scaling. If $K=-K$ then $\operatorname{dist}(a,b)=\operatorname{dist}(b,a)$ and $\operatorname{dist}(0,\cdot)$ is a norm in $\mathbb{R}^{d}$ in the usual sense. We allow $K\neq-K$, whence $\operatorname{dist}(a,b)\neq\operatorname{dist}(b,a)$ in general, but we still call it a norm. Bisectors and Voronoi diagrams for these norms have been studied in computational geometry [4, Sect. 7.2], [23, Sect. 4].

For any finite point set $S$ we define its bisector:

[TABLE]

Following the computational geometry tradition we will often call the elements of $S$ the sites. Although some of our results also hold for general convex bodies, for simplicity we assume $K$ to be a polytope. We denote by $\mathcal{F}(K)$ the face fan of $K$. The norm $\operatorname{dist}(0,\cdot)$ is linear in each of these cones, so we write

[TABLE]

for the intersection of the bisector with a choice of cones $F_{i}\in\mathcal{F}(K)$. Each cell of the form $\operatorname{bis}_{(F_{1},\dots,F_{k})}(a_{1},\dots,a_{k})$ is the intersection of the polyhedron $(a_{1}+F_{1})\cap\dots\cap(a_{k}+F_{k})$ with an affine subspace, which implies it is itself a polyhedron. As a consequence:

Proposition 1**.**

Let $K$ be a polytope with the origin in its interior, and let $\operatorname{dist}$ be the corresponding Minkowski norm. Let $S=\{a_{1},\dots,a_{k}\}\subset\mathbb{R}^{d}$ be a finite point set. Then the bisector $\operatorname{bis}(\{a_{1},\dots,a_{k}\})$ is a polyhedral complex whose cells are the polyhedra

[TABLE]

for all choices of $F_{1},\dots,F_{k}\in\mathcal{F}(K)$.

Proof.

The family of polyhedra

[TABLE]

forms a polyhedral complex since

[TABLE]

That polyhedral complex covers the entire bisector since for each point $p\in\operatorname{bis}(S)$ and for each $i$, the point $a_{i}$ must lie in some face $F_{i}$ of $p-\operatorname{dist}(a_{i},p)K$. ∎

Our primary example is the case where $K=\mathbb{B}^{d}$ is the tropical ball. In Figure 2 the point $p$, which is generic within the bisector of $a$ and $b$, lies in the facet $\operatorname{bis}_{(-*+),(+-*)}(a,b)$. For the purpose of drawing pictures, notice that any three vectors $v_{1},v_{2},v_{3}\in\mathbb{R}^{2}$ with $v_{1}+v_{2}+v_{3}=0$ define a map from $\mathbb{R}^{3}/\mathbb{R}\mathbb{1}$ to $\mathbb{R}^{2}$ via $e_{i}\mapsto v_{i}$. While $v_{1}=(1,0)$, $v_{2}=(0,1)$, $v_{3}=(-1,-1)$ is a common base for diagrams in tropical geometry, for our pictures we settle for the more symmetric isometric view where our base is:

[TABLE]

Remark 3**.**

If $\operatorname{bis}(a_{1},\dots,a_{k})=\emptyset$ then $\operatorname{bis}(a^{\prime}_{1},\dots,a^{\prime}_{k})=\emptyset$ for every choice of $a^{\prime}_{1},\dots,a^{\prime}_{k}$ sufficiently close to $a_{1},\dots,a_{k}$, by continuity of $\operatorname{dist}(\cdot,\cdot)$.

3.1. Weak general position and general position

Definition 1** (General position).**

A finite point set $S\subset\mathbb{R}^{d}$ is in general position with respect to $K$, if for every subset $\{a_{1},\dots,a_{k}\}\subset S$ there are neighborhoods $U_{i}$ of each $a_{i}$ such that for every choice of $\{a^{\prime}_{1},\dots,a^{\prime}_{k}\}$ with $a^{\prime}_{i}\in U_{i}$ and for every choice of maximal cones $F_{1},\dots,F_{k}\in\mathcal{F}(K)$ we have

[TABLE]

Moreover, the set $S$ is in weak general position if no pair of points $a,b\in S$ lies in a hyperplane parallel to a facet of $K$.

Remark 4**.**

As the name suggests, “weak general position” is implied by “general position” (see Corollary 1). A yet stronger notion would arise requiring stability not only of facets but also of lower-dimensional cells in bisectors; that is, allowing lower-dimensional cones from $\mathcal{F}(K)$ for the $F_{i}$ in the definition of general position. But this intermediate notion of “general position” is the most appropriate for our purposes and the algorithms in Section 5. As a first indication, Theorem 1 provides a local characterization.

By Proposition 1 bisectors are polyhedral complexes, and thus each cell has a dimension.

Proposition 2**.**

Two points $a,b$ are in weak general position if and only if $\operatorname{bis}(a,b)$ does not contain full-dimensional cells.

Proof.

Since $\operatorname{bis}_{(F,F^{\prime})}(a,b)\subset(a+F)\cap(b+F^{\prime})$, for it to be $d$-dimensional we need $F$ and $F^{\prime}$ to be cones of facets. We also need $F=F^{\prime}$, so that $\operatorname{dist}(a,\cdot)$ and $\operatorname{dist}(b,\cdot)$ have the same gradient on $(a+F)\cap(b+F)$, and we need $b-a$ to be parallel to the facet, so that $\operatorname{dist}(a,\cdot)=\operatorname{dist}(b,\cdot)$ on $(a+F)\cap(b+F)$.

Conversely, if $b-a$ is parallel to a facet of $K$ with cone $F$ then

[TABLE]

and this is $d$-dimensional. ∎

Corollary 1**.**

General position implies weak general position.

Proof.

Suppose $a-b$ is parallel to a facet of $K$, so that $\operatorname{bis}_{(F,F)}(a,b)$ is full-dimensional, where $F$ is the cone of that facet. Taking $b^{\prime}$ close to $b$ but away from the hyperplane parallel to the facet makes $\operatorname{bis}_{(F,F)}(a,b^{\prime})$ empty. Now the claim follows from Proposition 2. ∎

Theorem 1**.**

Let $S=\{a_{1},\dots,a_{k}\}\subseteq\mathbb{R}^{d}$ and for each $a_{i}$ choose a maximal cone $F_{k}\in\mathcal{F}(K)$. Let $Q:=(a_{1}+F_{1})\cap\dots\cap(a_{k}+F_{k})$, let $\lambda_{F_{i}}(x)$ be the linear function that restricts to $\operatorname{dist}(0,x)$ on $F_{i}$, and let $H$ be the affine subspace defined by $\lambda_{F_{1}}(x-a_{1})=\dots=\lambda_{F_{k}}(x-a_{k})$.

Then, the following conditions are equivalent:

1. (1)

There are neighborhoods $U_{i}$ of each $a_{i}$ such that for any choice of $a^{\prime}_{i}\in U_{i}$ the polyhedron $\operatorname{bis}_{(F_{1},\dots,F_{k})}(a^{\prime}_{1},\dots,a^{\prime}_{k})$ is not empty. 2. (2)

1. (a)

$Q$* is full-dimensional and $H$ intersects its interior; and* 2. (b)

the $k-1$ functions $\lambda_{F_{i}}-\lambda_{F_{1}}$ for $i=2,\dots,k$ are linearly independent.

Since

[TABLE]

condition (a) is equivalent to “$\operatorname{bis}_{(F_{1},\dots,F_{k})}(a_{1},\dots,a_{k})$ meets the interior of $Q$”.

Proof.

For the implication from “1” to “2”, let us first show that “1” forces $Q$ to be full-dimensional. Aiming for a contradiction, we assume that $Q$ is contained in the boundary of one of the cones $a_{i}+F_{i}$. That is, the polyhedron $Q_{i}:=\bigcap_{j\neq i}(a_{j}+F_{j})$ does not meet the interior of $a_{i}+F_{i}$, for some $i$. Then any $a^{\prime}_{i}$ in the interior of $a_{i}+F_{i}$ yields $(a^{\prime}_{i}+F_{i})\cap Q_{i}=\emptyset$. Hence $\operatorname{bis}_{(F_{1},\dots,F_{k})}(a^{\prime}_{1},\dots,a^{\prime}_{k})=\emptyset$, where $a^{\prime}_{j}=a_{j}$ for $j\neq i$. This contradicts “1” and shows that $Q$ is full-dimensional.

To see that $\operatorname{bis}_{(F_{1},\dots,F_{k})}(a_{1},\dots,a_{k})$ must intersect the interior of $Q$, suppose we are given neighborhoods $U_{i}$ as in “1”; recall (4). For each $i$, choose $v_{i}$ in the interior of $F_{i}$ and such that $\lambda_{i}(v_{i})=1$. Let $a^{\prime}_{i}:=a_{i}+\varepsilon v_{i}$, where $\varepsilon>0$ is taken small enough so that $a^{\prime}_{i}\in U_{i}$. Observe that our choice of $v_{i}$ makes the affine subspace defined by $\lambda_{F_{1}}(x-a^{\prime}_{1})=\dots=\lambda_{F_{k}}(x-a^{\prime}_{k})$ agree with $H$. Let $Q^{\prime}:=(a^{\prime}_{1}+F_{1})\cap\dots\cap(a^{\prime}_{k}+F_{k})$, which lies in the interior of $Q$. We have

[TABLE]

By “1” this is not empty, and thus $H$ intersects the interior of $Q$.

For condition (b) observe that $H$ can equivalently be defined by the $k-1$ affine equalities

[TABLE]

If the left-hand sides are linearly dependent, then one of the $k-1$ equations, say the $i$th one, is redundant. But then choosing a point $a^{\prime}_{i}\in U_{i}$ with $\lambda_{F_{i}}(a^{\prime}_{i})\neq\lambda_{F_{i}}(a_{i})$ (and letting $a^{\prime}_{j}=a_{j}$ for the rest) renders the system of equations infeasible. Hence $\operatorname{bis}_{(F_{1},\dots,F_{d+1})}(a^{\prime}_{1},\dots,a^{\prime}_{d+1})=\emptyset$, contradicting “1”.

We now show that “2” implies “1”. Consider arbitrary points $a^{\prime}_{i}$, and let $Q^{\prime}_{i}=\bigcap_{i}(a^{\prime}_{i}+F_{i})$. Further let $H^{\prime}$ be the affine subspace defined by

[TABLE]

so that

[TABLE]

We want to show that if each $a^{\prime}_{i}$ is sufficiently close to the corresponding $a_{i}$ for all $i$ then $Q^{\prime}\cap H^{\prime}$ is not empty.

Condition (b) says that $H^{\prime}$ is $(d{+}1{-}k)$-dimensional (and parallel to $H$) for every choice of $a^{\prime}_{i}$s and that it varies continuously with the choice. Condition (a) says that $Q^{\prime}$ stays full-dimensional if the $a^{\prime}_{i}$s are sufficiently close to the original $a_{i}$s, and that it also varies continuously with the choice, in the following strong sense: consider a description of each $a_{i}+F_{i}$ by a finite system of linear inequalities. Then $a^{\prime}_{i}+F_{i}$ is defined by a system with the same linear functions and with right-hand sides varying continuously with the $a^{\prime}_{i}$s.

Thus, if each $a^{\prime}_{i}$ is close to $a_{i}$, then $Q^{\prime}$ and $H^{\prime}$ are a full-dimensional polyhedron and a $(d{+}1{-}k)$-dimensional affine subspace close to $Q$ and $H$ respectively. Since, by (a), $H$ intersects the interior of $Q$, continuity implies that $H^{\prime}$ still intersects the interior of $Q^{\prime}$ when $a^{\prime}_{i}$ is close enough to $a_{i}$ for each $i$. In particular, $\operatorname{bis}_{(F_{1},\dots,F_{k})}(a^{\prime}_{1},\dots,a^{\prime}_{k})$ is not empty. ∎

Corollary 2**.**

The bisector of $k$ points in general position is either empty or pure of dimension $d+1-k$. In particular, the bisector of $d+1$ points in general position is finite, and this is empty for more than $d+1$ points.

Proof.

Every maximal non-empty cell $\operatorname{bis}_{(F_{1},\dots,F_{k})}(a_{1},\dots,a_{k})$ is the intersection of the polyhedron $Q$ with the affine subspace $H$ of Theorem 1. That result implies that $H$ has dimension $d+1-k$ and meets the interior of the full-dimensional polyhedron $Q$. Thus, the cell has dimension $d+1-k$. ∎

Corollary 3**.**

If every subset of at most $d+2$ points in $S$ is in general position then so is $S$.

Corollary 4**.**

For any $n\geq 1$, the sequences of $n$ points in $\mathbb{R}^{d}$ in general position form an open dense subset of $(\mathbb{R}^{d})^{n}$.

Proof.

Let $S\subset\mathbb{R}^{d}$ be a set of cardinality $n$. Then $S$ is in general position if and only if, for each subset $\{a_{1},\dots,a_{k}\}\subset S$ and maximal cones $F_{1},\dots,F_{k}\in\mathcal{F}(K)$, the polyhedron $\operatorname{bis}_{(F_{1},\dots,F_{k})}(a_{1},\dots,a_{k})$ either is empty or satisfies condition (1) of Theorem 1. Since there are finitely many choices of $\{a_{1},\dots,a_{k}\}$ and $\{F_{1},\dots,F_{k}\}$, it suffices to prove the statement for one such choice.

For openness: Emptyness is an open condition because of the continuity of $\operatorname{dist}(\cdot,\cdot)$. On the other hand, Theorem 1 says that condition (1) is equivalent to (2), which is open by the arguments in the proof (namely, the fact that both $Q$ and $H$ depend continuously on the sites).

For density, we are going to show that if $Q\cap H$ is not empty but fails to satisfy condition (a) or (b) then the set $\{a_{1},\dots,a_{k}\}$ lies in one of finitely many linear hyperplanes in $(\mathbb{R}^{d})^{k}$. To emphasize that $Q$ and $H$ depend on the choice of sites we denote them $Q(a_{1},\dots,a_{k})$ and $H(a_{1},\dots,a_{k})$.

If (b) fails for $\{a_{1},\dots,a_{k}\}$ then the linear system (5) defining $H(a_{1},\dots,a_{k})$ is feasible but overdetermined, which implies a linear relation, depending solely on $F_{1},\dots,F_{k}$, among the right-hand sides $\lambda_{F_{i}}(a_{i})-\lambda_{F_{1}}(a_{1})$. The relation is not tautological on the $a_{i}$s since, as shown in the proof of Theorem 1, it is easy to construct a point set $\{a^{\prime}_{1},\dots,a^{\prime}_{k}\}$ with $H(a^{\prime}_{1},\dots,a^{\prime}_{k})=\emptyset$.

For (a), consider the inequality descriptions of the cones $F_{i}$ and translate them to obtain an inequality description of $Q(a_{1},\dots,a_{k})$ as the feasibility region of a system $\mathcal{S}(a_{1},\dots,a_{k})$ of affine inequalities with fixed gradients and with right-hand sides parameterized linearly by the $a_{i}$s. If (a) fails for $\{a_{1},\dots,a_{k}\}$ then one of two things happen:

• •

$Q(a_{1},\dots,a_{k})$ is non-empty but not full-dimensional. Consider a minimal subsystem of $\mathcal{S}(a_{1},\dots,a_{k})$ that already defines a non-full-dimensional feasibility region. Minimality implies that this feasibility region is an affine subspace of $\mathbb{R}^{d}$ and that turning the inequalities to equalities produces an over-determined subsystem. This implies, as in the previous case, a linear relation among the $a_{i}$s.

• •

$Q(a_{1},\dots,a_{k})$ is full-dimensional but $Q(a_{1},\dots,a_{k})\cap H(a_{1},\dots,a_{k})$ is contained in its boundary. This implies $Q(a_{1},\dots,a_{k})\cap H(a_{1},\dots,a_{k})$ to be contained in a facet of $Q(a_{1},\dots,a_{k})$. Let $H_{0}(a_{i})$ be the hyperplane containing that facet. Note that the hyperplane $H_{0}$ depends only on one $a_{i}$ since it comes from the description of one of the cones $a_{i}+F_{i}$. Then, adding to the the linear system (5) the equation defining $H_{0}(a_{i})$ produces again an over-determined system, hence there is a linear relation among the $a_{i}$s.

The relations for condition (a) are not tautological on the $a_{i}$s since we can easily make $Q(a^{\prime}_{1},\dots,a^{\prime}_{k})$ full-dimensional and $H(a^{\prime}_{1},\dots,a^{\prime}_{k})$ intersect its interior as follows: choose each $a^{\prime}_{i}$ in the interior of $-F_{i}$ and such that $\lambda_{F_{i}}(a_{i})=-1$, so that both $H$ and the interior of $Q$ contain the origin. The relations are finitely many since we have at most one for each subsystem (respectively equation) of the system $\mathcal{S}(a_{1},\dots,a_{k})$, which depends only on the choice of $F_{i}$s. ∎

For the case of the tropical norm, condition (b) admits a nice combinatorial characterization. Observe that a choice of facets $F_{1},\dots,F_{k}\in\mathcal{F}(\mathbb{B}^{d})$ can be encoded as a directed graph on the vertex set $[d+1]$ and with an arc $a_{i}$ going from the coordinate that is minimized at $F_{i}$ to the coordinate that is maximized at $F_{i}$, for $i=1,\dots,k$. We denote this graph $G(F_{1},\dots,F_{k})$.

Proposition 3**.**

For the case of the tropical norm, condition 2.(b) of Theorem 1 holds if and only if the graph $G(F_{1},\dots,F_{k})$ either has no (undirected) cycle or it has a unique cycle and it is unbalanced; that is, the number of arcs in one direction is different from the other direction.

Proof.

A cycle in $G(F_{1},\dots,F_{k})$ is equivalent to a linear dependence among the corresponding linear functions $\lambda_{F_{i}}$s, by simply adding them with signs corresponding to the direction of the arcs along the cycle. If the cycle is balanced then $\lambda_{F_{1}}$ can be subtracted from each $\lambda_{F_{i}}$ so that the corresponding functions $(\lambda_{F_{i}}-\lambda_{F_{1}})$ are also dependent. The same thing can be done if $G(F_{1},\dots,F_{k})$ has two different (unbalanced) cycles, since a linear combination of the two corresponding dependences can be made balanced.

Conversely, any linear dependence among the functions $(\lambda_{F_{i}}-\lambda_{F_{1}})$ corresponds to a balanced dependence among the corresponding $\lambda_{F_{i}}$. The latter either corresponds to a balanced circuit in the graph or decomposes into two (or more) linear dependences with distinct supports. ∎

As a consequence of Corollary 2 the bisector of a set $S$ of $d+1$ points in general position is a finite set of points, which we call circumcenters of $S$. In dimension two, three points in (weak) general position have a most one circumcenter, as we show in Corollary 6 below. In higher dimension the same is known to be false for other polytopal norms [13], and here is an example for the tropical norm:

Example 1** (Non-uniqueness of circumcenters).**

Let us consider the four points $a_{1}=(0,2,3,3)$, $a_{2}=(0,4,2,2)$, $a_{3}=(2,4,1,1)$ and $a_{4}=(4,0,2,2)$. Their bisector contains the points $x=(0,0,1,-1)$ and $y=(0,0,-1,1)$. Indeed, both $x$ and $y$ are at distance $4$ from all the $a_{i}$’s since we have

[TABLE]

The points $a_{1},\dots,a_{4}$ are not in weak general position, as they lie in the plane $x_{3}-x_{4}=0$. However, they satisfy conditions (a) and (b) of Theorem 1 for the polytopes $Q_{x}$ and $Q_{y}$ containing the circumcenters $x$ and $y$: For condition (b) observe that the digraphs corresponding to $x$ and $y$ are, respectively, $\{14,12,32,21\}$ and $\{13,12,42,21\}$. They both have a single cycle, $\{12,21\}$, which is unbalanced. Condition (a) follows from the fact that $x$ (resp. $y$) is in the interior of all the cones whose intersection defines $Q_{x}$ (resp. $Q_{y})$. This is equivalent to the fact that all the vectors $a_{i}-x$ and $a_{i}-y$ have a unique maximum and a unique minimum entries.

Since condition (b) does not depend on the sites and condition (a) is, by Theorem 1, open, any perturbation of the sites will still produce at least two circumcenters. In particular, by Corollary 4, there are sites in general position for which this happens.

3.2. Halfspheres, sectors, and the bisector of two points

The topology of a bisector is closely related to the following partition of $\partial K$. Let $S\subseteq\mathbb{R}^{d}$ be a finite set of sites, and pick $a,b\in S$ distinct. The open halfsphere in the direction of $b-a$, denoted as $H(b-a)$, is the set of points in $\partial K$ whose exterior normal cone is contained in $(b-a)^{\vee}:=\{\lambda\,|\,\lambda(b-a)>0\}$. Informally, $H(b-a)$ are the faces of $K$ that “can be seen” from the direction $(b-a)$. For a fixed site $a\in S$, the sector of $a$ is the set

[TABLE]

We denote $\mathcal{H}_{S}:=\{H_{S}(a)\,|\,a\in S\}$. Observe that $H(b-a)$ and, hence, $H_{S}(a)$, are open in $\partial K$.

Lemma 1**.**

Let $F_{1},\dots,F_{m}$ be the facets of $K$ and let $\lambda_{F_{i}}(x)\leq 1$ be the valid linear inequality defining $F_{i}$. Then, for each $a\in S$,

[TABLE]

In particular, $H_{S}(a)\cap H_{S}(b)=\emptyset$ for every $a,b\in S$ and, if $S$ is in weak general position,

[TABLE]

where $\overline{H_{S}(a)}$ denotes the topological closure of $H_{S}(a)$.

Proof.

It is clear from the definition that $H_{S}(a)$ contains the relative interior of every $F_{i}$ with $\lambda_{F_{i}}(a)<\lambda_{F_{i}}(b)$ for $b\in S\setminus a$. By convexity of the cones $\{\lambda\,|\,\lambda(b-a)>0\}$ for each $a,b$, $H(b-a)$ (hence $H_{S}(a)$) also contains the relative interior of every lower dimensional face contained only in such facets. This proves the first formula. The second part follows from the first and the fact that in weak general position the minimum of each $\lambda_{F_{i}}$ is attained at a single point of $S$. ∎

Remark 5**.**

Assuming weak general position, Lemma 1 allows us to think of $\mathcal{H}_{S}$ as a labeling of the facets of $K$ by the elements of $S$ or, equivalently, as a map $\mathcal{F}(K)\to S$. If $K$ is centrally symmetric, then each pair of opposite facets $F$ and $-F$ belong one to $H(b-a)$ and the other to $H(a-b)$. If $K$ is not, we can still guarantee that $H(a-b)$ is never empty, and always disjoint from $H(b-a)$. As a consequence, $H(b-a)$ (and hence $H_{S}(a)$) cannot contain all the (relative interiors of) facets of $K$.

For the case $K=\mathbb{B}^{d}$ of the tropical ball this partition of the facets translates into something more meaningful. Recall (see Proposition 3 and the paragraph before it) that facets of $\mathbb{B}^{d}$ can be represented as the arcs in the complete digraph on $d+1$ nodes. In particular, $\mathcal{H}_{S}$ colors the arcs by the points of $S$. Then:

1. (1)

Each coloring is a partial ordering of vertices, i.e. there is no monochromatic cycle. 2. (2)

For the case of two points in general position, the two colors are opposite acyclic tournaments. In particular, there is a bijection between the possible halfspheres $H(b-a)$ and the total orderings of $d+1$ elements.

Theorem 2** ([13] for $d=3$).**

Let $a,b\in\mathbb{R}^{d}$ be in weak general position. Then the central projection from $a$ induces a homeomorphism between $\operatorname{bis}(a,b)$ and $a+H(b-a)$. Hence, $\operatorname{bis}(a,b)$ is homeomorphic to $\mathbb{R}^{d-1}$.

Proof.

Let us first show that $\operatorname{bis}(a,b)$ is contained in $a+\operatorname{cone}(H(b-a))$. To seek a contradiction, let $c\in\operatorname{bis}(a,b)$ such that $c-a\notin\operatorname{cone}(H(b-a))$. This implies that the smallest ball centered at $a$ that contains $c$ touches it at a facet $F$ with functional $\lambda_{F}$ such that $\lambda_{F}(b-a)\leq 0$. Now, $c$ is equidistant to $a$ and $b$, and $a$ and $b$ cannot be in the same facet of the ball centered at $c$ (because they are in weak general position). Therefore, $\operatorname{dist}(c,\varepsilon a+(1-\varepsilon)b)<\operatorname{dist}(c,a)$, by convexity of the ball. This contradicts the fact that $\lambda_{F_{i}}(b-a)\leq 0$.

Hence, we have a well-defined map $\phi:\operatorname{bis}(a,b)\to a+H(b-a)$ given by central projection. The map $\phi$ is continuous since it is the restriction of central projection. It is also proper (that is, the inverse image of a compact set is compact) by a following argument: Let $C$ be a compact subset of $H(b-a)$. By continuity, $\phi^{-1}(C)$ is closed in $\operatorname{bis}(a,b)$, hence in $\mathbb{R}^{d}$ since $\operatorname{bis}(a,b)$ itself is closed (it is the zero set of the continuous function $d(x,a)-d(x,b)$). Thus, we only need to prove that $\phi^{-1}(C)$ is bounded. This follows from the fact that

[TABLE]

and that $\operatorname{cone}(C)$ and $\operatorname{cone}(\overline{H(a-b)})$ are two closed linear cones meeting only at the origin, since ${H(a-b)}$ and ${H(b-a)}$ are open and disjoint in $\partial K$.

Once we know $\phi$ is proper and continuous, we only need to check that this map is bijective in order for it to be a homeomorphism. To show this, we construct its inverse. For each $v\in H(b-a)$ we consider the ray $r_{v}=\{a+\alpha v:\alpha\geq 0\}$. Along $r_{v}$, the distance to $a$ is linear in $\alpha$, the distance to $b$ is convex in $\alpha$ and both functions are continuous. Observe also that

[TABLE]

The last inequality comes from the fact that as the we move farther away from $a$ along $r_{v}$, eventually $(a+\alpha v)-b$ will be in the same cone of $\mathcal{F}(K)$ as $(a+\alpha v)-a=\alpha v$ (by weak general position), and $\langle b-a,v\rangle>0$ since $v\in H(b-a)$.

Hence, the function $\alpha\mapsto\operatorname{dist}(a+\alpha v,b)-\operatorname{dist}(a+\alpha v,a)$ is negative at zero, positive at infinity, continuous, and convex. Therefore, it has exactly one root, which means $r_{v}$ intersects the bisector exactly once. We define $\psi(a+v)$ as this unique intersection point.

The maps $\phi$ and $\psi$ are clearly inverses of one another. ∎

Looking at the proof, the reader can check that central projection gives a proper and continuous map from $\operatorname{bis}(a,b)$ to $a+H(b-a)$ even without assuming weak general position. We only need weak general position to construct the inverse.

Corollary 5**.**

If $S$ is in weak general position and there is an empty sector $H_{S}(a)\in\mathcal{H}_{S}$ then the bisector of $S$ is empty.

Proof.

Assume that there is a point $c\in\operatorname{bis}(S)$. For a site $a\in S$ let us show that $H_{S}(a)\neq\emptyset$. By definition, $c\in\operatorname{bis}(a,b)$ for $b\in S\setminus\{a\}$. By Proposition 2, each $\operatorname{bis}(a,b)$ can be mapped to $H(b-a)$ by central projection. Since $c$ is in all of these bisectors, the central projection of $c$ into the ball $a+K$ lies in $H(b-a)$ for all $b$, and hence in $H_{S}(a)$. ∎

The converse of Proposition 5 is true for three points in arbitrary dimension (Theorem 3) but not for more, even in general position, as the following example shows:

Example 2** (Empty bisector, with non-empty sectors).**

Consider $a=(1,-1,$ $0,0)$, $b=(-1,1,0,0)$, $c=(0,0,2,-2)$ and $d=(0,0,-2,2)$. Then we have

[TABLE]

By symmetry, we also have

[TABLE]

Since $\operatorname{bis}(a,b,c,d)$ lies in the intersection of the two, we have

[TABLE]

So for $x\in\operatorname{bis}(a,b,c,d)$, we may assume $x_{3}=0$, which entails:

[TABLE]

with equality only when $x_{1}=0$ and

[TABLE]

with equality only if $|x_{1}|\geq 2$. This shows that $\operatorname{bis}(a,b,c,d)=\emptyset$.

However, $\mathcal{H}_{S}$ has no empty sector since the sectors of $a,b,c,d$ contain the facets with outer normals $(1,-1,0,0)$, $(-1,1,0,0)$, $(0,0,1,-1)$ and $(0,0,$ $-1,1)$, respectively. Both this property and the emptiness of $\operatorname{bis}(a,b,c,d)$ are preserved under perturbation, so this behavior happens also in general position.

3.3. Bisectors of three points

The goal of this section is to prove our first main result.

Theorem 3**.**

Let $S=\{a_{1},a_{2},a_{3}\}$ be a set of three distinct points in $\mathbb{R}^{d}$ which lie in weak general position with respect to a polytope $K$. If $H_{S}(a_{i})\neq\emptyset$ for $i=1,2,3$ then $\operatorname{bis}(a_{1},a_{2},a_{3})$ is homeomorphic to a non-empty open subset of $\mathbb{R}^{d-2}$.

Corollary 6**.**

For any three points in weak general position $\operatorname{bis}(a_{1},a_{2},a_{3})$ is either empty or pure of dimension $d-2$. If $d=2$ then $\operatorname{bis}(a_{1},a_{2},a_{3})$ is either empty or a single point.

We begin with the case $d=2$ of Theorem 3; this occurs in [13]:

Lemma 2**.**

Let $S=\{a_{1},a_{2},a_{3}\}\subset\mathbb{R}^{2}$ be in weak general position with respect to a polytope $K$. If $H_{S}(a_{i})\neq\emptyset$ for the three of them, then the bisector $\operatorname{bis}(S)$ is a single point.

Proof.

We first show that $\operatorname{bis}(S)$ cannot contain two distinct points. To seek a contradiction, suppose that there are two points $x\neq y$ in the bisector. This means that there exist $\alpha,\beta>0$ such that $B_{x}:=x-\alpha K$ and $B_{y}:=y-\beta K$ satisfy $S\subset\partial B_{x}\cap\partial B_{y}$. For brevity we call $B_{x}$ and $B_{y}$ negative balls centered at $x$ and $y$, respectively. The two balls are related to one another by a homothety $\rho:B_{x}\to B_{y}$. Let $a^{\prime}_{i}:=\rho(a_{i})$ and $S^{\prime}=\{a^{\prime}_{1},a^{\prime}_{2},a^{\prime}_{3}\}=\rho(S)$. Then, $S\cup S^{\prime}\subset\partial B_{y}$, but this is impossible: the vertex sets of two homothetic triangles can only lie in the boundary of a polytope in the plane if three of the points are collinear. Three collinear points of $\partial B_{y}$ necessarily lie on a single edge of $B_{y}$, and (at least) two of them would come from the same triangle $S$ or $S^{\prime}$, violating the weak general position of $S$.

It remains to exclude the case where $\operatorname{bis}(S)=\emptyset$. Suppose that this would hold. That is, the two-point bisectors are pairwise disjoint. Now the two-point bisectors are homeomorphic to lines by Theorem 2. So the fact that they do not meet implies that each of them appears either in full or not at all in the Voronoi diagram of $S$. But the three of them cannot appear, since then the diagram would have four regions, not three. Thus, one of them, say $\operatorname{bis}(a_{1},a_{3})$, does not appear at all in $\operatorname{Vor}_{S}$. Consequently, $\operatorname{bis}(a_{1},a_{3})$ is contained in the Voronoi region of the third point $a_{2}$. We will show that $H_{S}(a_{2})=\emptyset$, and this yields the desired contradiction.

To simplify the exposition, we call the line $a_{1}a_{3}$ “horizontal”. Let $u$ and $v$ be the points where the ball $K$ has a horizontal tangent, i.e., parallel to $a_{1}a_{3}$. Without loss of generality, $u$ is at the top of $K$, and $v$ is at the bottom; see Figure 3. The points $u$ and $v$ are unique and thus vertices of $K$, by weak general position. Consider the negative tangent cones $u+\operatorname{cone}(u-K)$ and $v+\operatorname{cone}(v-K)$ of $K$. Call $K_{u}$ and $K_{v}$ their respective translations that have $a_{1}$ and $a_{3}$ on the boundary, and let $u^{\prime}$ and $v^{\prime}$ be the apices of these translated cones. Such translations exist and are unique since the two boundary rays of $\operatorname{cone}(u-K)$ point upward, and the boundary rays of $\operatorname{cone}(v-K)$ point downward.

By construction the boundaries of $K_{u}$ and $K_{v}$ intersect precisely in $a_{1}$ and $a_{3}$. Now, every negative ball of the form $B_{u^{\prime},r}:=u^{\prime}-r(K-u)$ has $K_{u}$ as its tangent cone at $u^{\prime}$, and balls like $B_{v^{\prime},r}:=v^{\prime}-r(K-v)$ are similar. Taking $r$ sufficiently large, we can force

[TABLE]

Since $B_{u^{\prime},r}$ and $B_{v^{\prime},r}$ are negative balls with $a_{1}$ and $a_{3}$ in their boundaries, their centers lie on the bisector $\operatorname{bis}(a_{1},a_{3})$. Our initial assumption that the bisector is contained in the Voronoi region of $a_{2}$ implies that $a_{2}$ is in the interior of both balls, i.e., in the interior of both $K_{u}$ and $K_{v}$. This implies that $H_{S}(a_{2})$ is empty: every facet-defining functional of $K$ takes its minimum on $K_{u}\cap K_{v}$ at either $a_{1}$ or $a_{3}$.

We conclude that, if $H_{S}(a_{i})\neq\emptyset$ for $i=1,2,3$ then $\operatorname{bis}(S)\neq\emptyset$, and this finishes the proof. ∎

For the rest of the proof of Theorem 3 let $S=\{a_{1},a_{2},a_{3}\}\subset\mathbb{R}^{d}$ be three points in general position with respect to a convex polytope $K$. The idea is to reduce the general problem to two dimensions via the following construction. Let $\operatorname{pr}:\mathbb{R}^{d}\rightarrow\mathbb{R}^{d-2}$ be the affine projection that quotients out the $2$-plane $\Pi$ containing $S$. Next we exhibit relevant properties of that map.

Lemma 3**.**

With the above notation, let $x\in\operatorname{int}(\operatorname{pr}(K))\subset\mathbb{R}^{d-2}$. Further let $\Pi_{x}:=\operatorname{pr}^{-1}(x)$, which is a $2$-plane parallel to $\Pi$, and $K_{x}:=K\cap\Pi_{x}$. Then $K_{x}$ is a convex polygon.

Proof.

This is a general property of polytope projections: if $Q$ is the image of a polytope $P$ under an affine map then the fiber of every point $x\in\operatorname{int}(Q)$ is a polytope of dimension $\dim(P)-\dim(Q)$. ∎

Lemma 4**.**

With the above notation, again let $x\in\operatorname{int}(\operatorname{pr}(K))$. Further, let $H_{S}^{(x)}(a_{i})$ denote the sector of $a_{i}\in\Pi$ with respect to the polygon $K_{x}$. Then

[TABLE]

Proof.

Note that no facet in $K_{x}$ is parallel to any $a_{i}-a_{j}$ because if it were, the facet of $K$ containing it would be parallel, too.

Let $F^{\prime}$ be a facet of $K_{x}$, and let $F$ be the corresponding facet in $K$. Then, $n(F^{\prime})\in\Pi_{x}$, the normal vector to $F^{\prime}$, is the projection into $\Pi$ of $n(F)$, and,

[TABLE]

∎

Lemma 5**.**

Let $S=\{a_{1},a_{2},a_{3}\}$ as before. If $H_{S}(a_{i})\neq\emptyset$ for all $i$, then the intersection $\bigcap_{a_{i}\in S}\operatorname{pr}(H_{S}(a_{i}))\subset\mathbb{R}^{d-2}$ is open and not empty.

Proof.

First, observe that an $x\in\partial K$ with $\operatorname{pr}(x)\in\partial(\operatorname{pr}(K))$ cannot be in any of the $H_{S}(a_{i})$: indeed, $x\in\partial K$ implies that there is a normal vector of $K$ at $x$ orthogonal to $\Pi$, hence orthogonal to $a_{i}-a_{j}$ for every $a_{i},a_{j}$. As a consequence,

[TABLE]

which implies it is open.

For any point $x\in\operatorname{int}(\operatorname{pr}(K))$, the preimage $\operatorname{pr}^{-1}(x)$ is a polygon, a slice of $K$. This slice has to intersect at least two of the classes of $\mathcal{H}$, because $\mathcal{H}$ is a partition (so at the slice intersects at least one class), but no class can contain a set of facets whose vectors are positively dependent (because each class is an intersection of half-spheres). Then, any point $x\in\operatorname{int}(K)$ lies in at least two sets $\operatorname{pr}(H_{S}(a_{i}))$.

Thus, the three open sets $\operatorname{pr}(H_{S}(a_{i}))$ cover each point of $\operatorname{int}(\operatorname{pr}(K))$ at least twice. At least two of these sets must intersect, say $\operatorname{pr}(H_{S}(a_{1}))$ and $\operatorname{pr}(H_{S}(a_{2}))$. Suppose that $(\operatorname{pr}(H_{S}(a_{1}))\cap\operatorname{pr}(H_{S}(a_{2})))\cap\operatorname{pr}(H_{S}(a_{3}))=\emptyset$. Then, $\operatorname{int}(\operatorname{pr}(K))$ would be disconnected, because it is covered by two disjoint open sets. Since this is not possible, there must be a point in the common intersection of the three $H_{S}(a_{i})$. ∎

Our next lemma finishes the proof of Theorem 3, and it actually gives more information.

Lemma 6**.**

The set $\operatorname{bis}(a_{1},a_{2},a_{3})$ is homeomorphic to $\bigcap_{i=1,2,3}\operatorname{pr}(H_{S}(a_{i}))$.

Proof.

Consider the map

[TABLE]

defined as follows: Let $p\in\operatorname{bis}(S)$, and let $v_{i}\in H_{S}(a_{i})$ be the central projection from $p$ to $a_{i}+\partial K$, for each $i=1,2,3$. Note that each $v_{i}$ lies in the corresponding $H_{S}(a_{i})$, by Theorem 2. Further, the three points $v_{1}$, $v_{2}$ and $v_{3}$ lie in a plane parallel to $\Pi$. In particular, $\operatorname{pr}(v_{1})=\operatorname{pr}(v_{2})=\operatorname{pr}(v_{3})$ lies in $\bigcap_{i=1,2,3}\operatorname{pr}(H_{S}(a_{i}))$ and we define

[TABLE]

To show that $\phi$ is a homeomorphism, let us construct its inverse $\psi$. Let $\gamma:\operatorname{pr}(\operatorname{int}(K))\to\operatorname{int}(K)$ be a continuous section of $\operatorname{pr}$ in $K$. For example, but not necessarily, for each $2$-plane $\Pi^{\prime}$ parallel to $\Pi$ and intersecting $K$ let $\gamma(\operatorname{pr}(\Pi^{\prime}))$ be the centroid of $\Pi^{\prime}\cap K$.

Now, let $x\in\bigcap_{i=1,2,3}\operatorname{pr}(H_{S}(a_{i}))$. Let $\Pi_{x}=\operatorname{pr}^{-1}(x)$ and let $w_{i}=\gamma(x)+a_{i}$, for each $i=1,2,3$. In the $2$-plane $\Pi_{x}$ we have a set $S_{x}=\{w_{1},w_{2},w_{3}\}$ and a unit ball $K\cap\Pi_{x}$. Lemma 3 gives that $H_{S_{x}}(w_{i})=H_{S}(a_{i})\cap\Pi_{x}$. By choice of $x$ we have $\bigcap_{i}H_{S_{x}}(w_{i})\neq\emptyset$ and Lemma 2 guarantees that the bisector of $S_{x}$ is a unique point $r\in\Pi_{x}$; see Figure 4.

Let $v_{i}$ be the central projection of $r$ to $w_{i}+\partial K_{x}$. Observe that $|w_{i}r|/|w_{i}v_{i}|$ is independent of $i$ since it equals $\operatorname{dist}_{K_{x}}(w_{i},r)$, where $\operatorname{dist}_{K_{x}}$ denotes the distance induced by $K_{x}$ in $\mathbb{R}^{2}$, and $r$ is in $\operatorname{bis}_{K_{X}}(w_{1},w_{2},w_{3})$. Since $w_{i}-a_{i}=\gamma(x)$ is also independent of $i$, the three rays $a_{i}v_{i}$ meet at the point

[TABLE]

and $\operatorname{dist}_{K}(a_{i},p)=\operatorname{dist}_{K_{X}}(w_{i},r)$. Thus, $p\in\operatorname{bis}(a_{1},a_{2},a_{3})$ and we define $p$ to be $\psi(x)$.

This gives us a well-defined map

[TABLE]

and by construction $\psi$ is the inverse of $\phi$ both ways. The map $\gamma$ is continuous, as the bisector of three points in the plane depends continuously on a continuous deformation of the unit ball. So $\psi$ is continuous. Thus, $\phi$ and $\psi$ are homeomorphisms between $\operatorname{bis}(a_{1},a_{2},a_{3})$ and $\bigcap_{i=1,2,3}\operatorname{pr}(H_{S}(a_{i}))$. Since the latter is not empty and open by Lemma 5, the former is homeomorphic to a non-empty open subset of $\mathbb{R}^{d-2}$. ∎

Our proof of Theorem 3 closely follows the proof of the $3$-dimensional case in [14]. There, it is additionally shown that the number of connected components of $\operatorname{bis}(a_{1},a_{2},a_{3})$ equals the total number of connected components of the three sectors $H_{S}(a_{1})$ minus two. We can extend this formula to higher dimension and to higher Betti numbers (over any field):

Theorem 4**.**

Let $a_{1},a_{2},a_{3}\in\mathbb{R}^{d}$ be three points in weak general position with respect to a polytope $K$ and assume that $H_{S}(a_{i})\neq\emptyset$ for all three. Then, for $j\in\{0,\dots,d-3\}$, we have

[TABLE]

Proof.

Consider the same projection $\operatorname{pr}:\mathbb{R}^{d}\to\mathbb{R}^{d-1}$ as before. Observe that, since $\Pi_{x}\cap H_{S}(a_{i})$ is empty or contractible for every plane $\Pi_{x}$ parallel to $\operatorname{pr}$, we have

[TABLE]

We now apply Alexander duality in the one-point compactification $\mathcal{S}$ of $\operatorname{int}(\operatorname{pr}(K))$, which is a sphere of dimension $d-2$. Alexander duality says that if $U$ is an open and locally contractible subset of $\mathcal{S}$ then

[TABLE]

In particular, if we let $C_{i}=\mathcal{S}\setminus\operatorname{pr}(H_{S}(a_{i}))$ we have

[TABLE]

and

[TABLE]

Yet $C_{1}$, $C_{2}$ and $C_{3}$ are pairwise disjoint except for the “point at infinity” of $\mathcal{S}$, because each point of $\operatorname{int}(\operatorname{pr}(K))$ lies in at least two of the sets $\operatorname{pr}(H_{S}(a_{i}))$. Thus, $\bigcup_{i=1}^{3}C_{i}$ is the topological wedge (or one-point sum) of $C_{1}$, $C_{2}$ and $C_{3}$, which makes the right-hand sides of the two last equations coincide. ∎

Remark 6**.**

One may ask how complicated the Betti numbers $\tilde{\beta}_{j}(H_{S}(a_{i}))$ in Theorem 4 can be. Equivalently, how complicated the topology of three point bisectors can be. Such a bisector is $(d-2)$-dimensional, so the relevant Betti numbers are $\tilde{\beta}_{0},\dots,\tilde{\beta}_{d-2}$. The last one, $\tilde{\beta}_{d-2}$, must vanish as $\tilde{\beta}_{j}(H_{S}(a_{i}))=\tilde{\beta}_{j}(\operatorname{pr}(H_{S}(a_{i})))$, and the latter is an open subset of $\mathbb{R}^{d-2}$. But $\tilde{\beta}_{d-3}$ can be non-zero, as the following example shows.

Example 3**.**

Let $K=\mathbb{B}^{3}$ be the tropical unit ball, i.e., $d=3$. The three sites $a=(0,0,4,4)$, $b=(-3,0,2,0)$ and $c=(0,-3,0,2)$ lie in weak general position in $\mathbb{R}^{4}/\mathbb{R}\mathbb{1}$. We describe each sector in terms of the facets whose relative interiors lie in that sector. The facet $F_{ij}$ is the one at which coordinate $i$ is minimized and $j$ is maximized. We obtain

[TABLE]

The sector $H(a)$ is disconnected; i.e., $\tilde{\beta}_{0}(H(a))$ is non-zero. By Theorem 4 $\operatorname{bis}(a,b,c)$ is disconnected, too.

Remark 7**.**

The results in this section and in Section 3.2 involve only weak general position. While we keep assuming that $K$ is a polytope for simplicity, these results generalize to an arbitrary convex body $K$. In that case $S$ is in weak general position if the boundary $\partial K$ contains no segment parallel to a difference $a-b$ with $a,b\in S$.

4. Classification of tropical bisectors of two points

4.1. Tropical bisectors and tropical hypersurfaces

In the classical case the bisector of two points is a degenerate quadric, namely the affine hyperplane perpendicular to the connecting line segment and which runs through the midpoint. The tropical analog is more interesting.

Proposition 4**.**

Let $a,b\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ be in weak general position. Then the homogeneous $\max$-tropical Laurent polynomial

[TABLE]

vanishes on the bisector $\operatorname{bis}(a,b)$. That is, the set $\operatorname{bis}(a,b)$ is contained in a $\max$-tropical hypersurface of degree $d+1$.

Proof.

Recall that a $\max$-tropical (Laurent) polynomial vanishes if the maximum is attained at least twice; see [22, §3.1]. First, we check that there are no duplicates among the terms in the representation (6) of $\phi(a,b)$. Assume the contrary, i.e., $x_{i}-a_{i}-x_{j}+a_{j}=x_{i}-b_{i}-x_{j}+b_{j}$ for some $i,j\in[d+1]$. Then $a_{j}-a_{i}=b_{j}-b_{i}$, which forces $\langle e_{j}-e_{i},b-a\rangle=(b_{j}-a_{j})-(b_{i}-a_{i})=0$. Thus $b-a$ is parallel to the facet of $\mathbb{B}^{d}$ with normal vector $e_{j}-e_{i}$. This was explicitly excluded in our assumption, and we arrive at the desired contradiction. We infer that the $2d(d+1)$ terms are pairwise distinct.

Let $x\in\operatorname{bis}(a,b)$. This means that $\operatorname{dist}(a,x)=\operatorname{dist}(b,x)$, and thus

[TABLE]

It follows that $\phi(a,b)$ vanishes at $x$.

The degree of the bisector tropical hypersurface can be read off any Laurent monomial like $x_{i}-x_{j}$ by adding $x_{1}+x_{2}+\dots+x_{d+1}$, which yields the true monomial $2x_{i}+\sum_{k\in[d+1]-\{i,j\}}x_{k}$. The latter has degree $2+(d+1-2)=d+1$. ∎

Proposition 4 yields a trivial algorithm to compute tropical bisectors in weak general position: enumerate the maximal cells of the tropical hypersurface defined by (6) and select those maximal cells that attain maxima in one monomial of type $x_{i}-a_{i}-x_{j}+a_{j}$ and one monomial of type $x_{k}-b_{k}-x_{\ell}+b_{\ell}$. This algorithm needs to go through the $\Theta(d^{4})$ choices of one monomial from the left and one from the right. This is worst case optimal, as we will prove in Corollary 8 that tropical bisectors can have $\Omega(d^{4})$ maximal cells.

Example 4**.**

The labeling of the faces of a tropical bisector does not need to be unique if $a$ and $b$ are not in weak general position. For instance, if $b-a=(-1,1,0)$ then

[TABLE]

is the only face; see Figure 5 (left).

4.2. The bisection fan

Normal equivalence of tropical bisectors is preserved by translation and scaling. In particular the equivalence class of $\operatorname{bis}(a,b)$ is uniquely determined by the direction of the vector $b-a$. The bisection fan $\mathcal{F}_{\text{bis}}^{d}$ is the complete polyhedral fan in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ whose relatively open cones are defined by “$x$ and $y$ lie in the same cone if and only if $\operatorname{bis}(0,x)$ and $\operatorname{bis}(0,y)$ are normally equivalent”. Put differently, two points $a,b\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ are in general position if and only if the difference $b-a$ lies in a maximal cone of $\mathcal{F}_{\text{bis}}^{d}$. In the rest of this section we show that $\mathcal{F}_{\text{bis}}^{d}$ is indeed a polyhedral fan and give an explicit description of it.

Recall that an ordered partition or total preorder on a finite set $S$ is a partition of $S$ into non-empty parts together with a total order on the parts. If the parts are denoted $S_{1},\dots,S_{k}$ (in this order), we can write $x\leq y$ meaning “$x\in S_{i}$ and $y\in S_{j}$ for some $i\leq j$”. In particular, for all $x,y\in S$ we have $x\leq y\leq x$ if and only if $x$ and $y$ lie in the same part.

Any real vector $v=(v_{1},\dots,v_{d+1})\in\mathbb{R}^{d+1}$ induces an ordered partition $S(v)$ of $[d+1]$ by putting together the coordinates that have the same value and ordering the groups according to their values. For example, the vector $v=(3,1,6,4,6,3,1)$ of length seven induces the partition $(\{2,7\},\{1,6\},\{4\}$, $\{3,5\})$ of the set $\{1,2,\dots,7\}$ into four parts. Note that the ordered partition $S(v)$ is constant on the class $v+\mathbb{R}\mathbb{1}$. Hence these ordered partitions are defined for points in the projective tropical torus $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$.

For $v,w\in\mathbb{R}^{d+1}$ with $S(v)=S(w)$ we have $S(v+w)=S(v)$ and, moreover, $S(\alpha v)=S(v)$ for any positive real $\alpha$. That is to say, the stratification of $\mathbb{R}^{d+1}$ by ordered partitions forms a complete polyhedral fan. In what follows we seek to refine that fan by recording which part, or which gap between parts, contains the midvalue

[TABLE]

The bisected ordered partition of $[d+1]$ induced by $v\in\mathbb{R}^{d+1}$ is the ordered partition $S(v)$ as defined above, together with the information of which (values of) parts are smaller, equal or greater than the midvalue $\mu(v)$; see also [9]. Equivalently, this is the ordered partition associated with the extended vector $(v,\mu(v))\in\mathbb{R}^{d+2}$. We denote by $\mathcal{F}_{\text{bop}}^{d}$ the fan of bisected ordered partitions of dimension $d$.

Remark 8**.**

The “finest” or “most generic” ordered partitions are the permutations, in which each part is a singleton. Hence, the fan of ordered partitions equals the normal fan of the permutahedron. This, in turn, coincides with the fan of regions in the braid arrangement or Coxeter arrangement of type $A_{d}$, which consists of the hyperplanes $\left\{x\,\vphantom{x_{i}=x_{j}}\right|\left.\vphantom{x}\,x_{i}=x_{j}\right\}$ for $1\leq i<j\leq d+1$. We denote this fan $\mathcal{F}(A_{d})$. It is intermediate between the central fan of the tropical ball (which is coarser) and the fan of bisected ordered partitions (which is finer):

[TABLE]

see Figure 6 for a visualization of the case $d=3$. Note that $\mathcal{F}(A_{d})$ is also the fan of weak general position: $a$ and $b$ are in weak general position if and only if $b-a$ lies in a full-dimensional cell.

Example 5**.**

In $d=2$, the fans $\mathcal{F}(A_{d})$ and $\mathcal{F}(\mathbb{B}^{d})$ coincide, they form the face fan of the regular hexagon $\mathbb{B}^{2}$. The bisection fan $\mathcal{F}_{\text{bop}}^{d}$ is the barycentric subdivision of it; see Figure 7. Excluding permutations of the coordinates, and sign inversion, we infer that there are three types of tropical bisectors in the plane, and these are shown in Figure 5. The type to the left is in weak general position but not in general position, the type to the left is in general position, and the type to the right is not even in weak general position.

Recall that the $\max$-tropical line segment spanned by two points, $a$ and $b$, is the set

[TABLE]

It is worth noting that the combinatorial types of tropical segments are classified by the braid fan $\mathcal{F}(A_{d})$; see [22, Prop. 5.11]. In this sense the fan of bisected ordered partitions $\mathcal{F}_{\text{bop}}^{d}$ classifies the combinatorial types of “bisected tropical segments”:

Proposition 5**.**

The bisected ordered partition of $b-a$ contains the same information as the combinatorial type of the tropical line segment $[a,b]$ together with the information of which part contains the midpoint.

Proof.

Suppose for simplicity that $a,b\in\mathbb{R}^{d+1}$ satisfy

[TABLE]

Then $[a,b]$ is the union of at most $d$ ordinary line segments, one for each subset of coordinates between a strict inequality in (7). That is, the combinatorics of the tropical segment is the same as the ordered partition of $b-a$. The midvalue $\mu(b-a)$ selects one of the ordinary segments. ∎

The goal of the rest of this section is to prove the following:

Theorem 5**.**

We have $\mathcal{F}_{\text{bis}}^{d}=\mathcal{F}_{\text{bop}}^{d}$. That is, given $a,b,a^{\prime},b^{\prime}\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ the polyhedra $\operatorname{bis}(a,b)$ and $\operatorname{bis}(a^{\prime},b^{\prime})$ are normally equivalent if and only if $b-a$ and $b^{\prime}-a^{\prime}$ induce the same bisected ordered partition of $[d+1]$.

4.3. Proof of Theorem 5

In (3) we defined the type $(F_{-},F_{*},F_{+})$ of a face $F$ of $\mathbb{B}^{d}$ or the corresponding face cone in $\mathcal{F}(\mathbb{B}^{d})$. For a pair of faces, $F$ and $G$, this gives rise to the following labeling partition of $[d+1]$:

[TABLE]

As a first step in the proof of Theorem 5, the following lemma characterizes when $\operatorname{bis}_{(F,G)}(a,b)$ is non-empty. Recall that this is the case if and only if there is a tropical ball touching $a$ and $b$ at faces $F$ and $G$, respectively.

Lemma 7**.**

Let $F$ and $G$ be a fixed pair of faces of $\mathbb{B}^{d}$ with the labeling partition defined as in (8). Further let $a,b\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$. Then the set $\operatorname{bis}_{(F,G)}(a,b)$ is not empty if and only if there exist $\gamma\in\mathbb{R}$ and $\delta\in[0,\infty)$ such that the following conditions are satisfied:

[TABLE]

Proof.

Let us assume that the face of the tropical bisector $\operatorname{bis}(a,b)$ defined by $(F,G)$ is non-empty. Then there is a point $x$ such that $\operatorname{dist}(a,x)=\operatorname{dist}(x,b)=\delta$ and $a-x\in F$ as well as $b-x\in G$. We set $\gamma_{a}=\min_{i\in[d+1]}(a_{i}-x_{i})$ and $\gamma_{b}=\min_{i\in[d+1]}(b_{i}-x_{i})$. The possible values for the coordinates of $a-x$ and $b-x$ are

[TABLE]

for some $\delta\geq 0$. Setting $\gamma=\gamma_{b}-\gamma_{a}$ the above translates into (9).

For the converse, note that going from (13) and (17) to (9) is the Fourier-Motzkin elimination of the variables $x_{i}$. Therefore, any $\gamma$ and $\delta\geq 0$ which are feasible for (9) can be lifted to a solution of (13) and (17). That is to say, we can set $\gamma_{a}=0$ and $\gamma_{b}=\gamma$, and the conditions in (13) and (17) yield a point $x\in\operatorname{bis}_{(F,G)}(a,b)$. ∎

Proposition 6**.**

Let $F$ and $G$ be a fixed pair of faces of $\mathbb{B}^{d}$. Then the set

[TABLE]

is both a polyhedral cone and a tropical cone, although perhaps not a tropical polyhedral cone.

Proof.

Let $a,a^{\prime},b,b^{\prime}\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ with $\operatorname{bis}_{(F,G)}(a,b)$ and $\operatorname{bis}_{(F,G)}(a^{\prime},b^{\prime})$ both non-empty. Since $\operatorname{bis}_{(F,G)}(a,b)\neq\emptyset$, by Lemma 7, there are scalars $\gamma$ and $\delta$ satisfying the conditions (9). Likewise there are certificates $\gamma^{\prime}$ and $\delta^{\prime}$ for $\operatorname{bis}_{(F,G)}(a^{\prime},b^{\prime})\neq\emptyset$. By linearity of the conditions (9) it follows that $\gamma+\gamma^{\prime}$ and $\delta+\delta^{\prime}$ certify that $\operatorname{bis}_{(F,G)}(a+a^{\prime},b+b^{\prime})\neq\emptyset$: for instance, we have $(b+b^{\prime}-a-a^{\prime})_{i}=\gamma+\gamma^{\prime}$ for $i\in L_{0}=(F_{-}\cap G_{-})\cup(F_{+}\cap G_{+})$. Since clearly $\alpha c\in C$ for all $c\in C$ and $\alpha\geq 0$ we conclude that $C$ is an ordinary cone. This cone is polyhedral as it is defined in terms of the finitely many linear conditions (9).

A similar argument shows that $C$ is also closed with respect to taking arbitrary $(\max,+)$-linear combinations: for instance, with the above notation we have $\max((b-a)_{i},(b^{\prime}-a^{\prime})_{i})=\max(\gamma,\gamma^{\prime})$ for $i\in L_{0}$. This shows that $C$ is a tropical cone. ∎

Corollary 7**.**

The bisection fan of tropical bisectors is a classical polyhedral fan, and a tropical (perhaps not tropical polyhedral) fan.

Proof.

We know that the feasibility region of a face $(F,G)$ is a tropical and classical polyhedral cone. Finite intersections of these cones are again tropical cones, classical cones, and polyhedral cones. Therefore, the feasibility region of a normal equivalence class, which is the intersection of the cones of its non-empty faces, is again a tropical and classical polyhedral cone. Hence, the whole fan has this structure. ∎

The following shows one direction of Theorem 5, namely, $\mathcal{F}_{\text{bis}}^{d}$ is coarser than $\mathcal{F}_{\text{bop}}^{d}$.

Lemma 8**.**

Let $F,G\in\mathcal{F}(\mathbb{B}^{d})$. Then, whether $\operatorname{bis}_{(F,G)}(a,b)$ is empty or not, for each $a,b\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ depends only on the bisected ordered partition of $b-a$.

Proof.

Consider the partition of $[d+1]$ into six sets $L_{0},L_{+},L_{-},L_{+1},L_{-1},L_{*}$ defined in (8). We want to show that feasibility of the system (9) for a given $a$ and $b$ depends only on the bisected ordered partition of $b-a$. Without loss of generality we assume $a=0$ and $b_{1}\leq b_{2}\leq\dots\leq b_{d+1}$ as in (7).

Let $\mu(b):=\tfrac{1}{2}(b_{d+1}+b_{1})=\tfrac{1}{2}\operatorname{dist}(0,b)+b_{1}$ be the midvalue of $b-0$ and $m=\max(0+\mu(b)\mathbb{1},b)$ the midpoint of the segment $[0,b]$. In particular, $\operatorname{dist}(0,m)=\operatorname{dist}(m,b)=\tfrac{1}{2}\operatorname{dist}(0,b)$.

We distinguish three cases, depending on whether both, none, or exactly one of $L_{+1}$ and $L_{-1}$ are empty.

Claim I: Suppose $L_{+1}\cup L_{-1}=\emptyset$. Then, (9) is feasible if and only if there are $k\leq\ell$ with

[TABLE]

Indeed, in this case feasibility of (9) is equivalent to feasibility of

[TABLE]

which implies the ordered partition to satisfy (I.1). Conversely, if the ordered partition satisfies (I.1), then let $\gamma$ be chosen to satisfy (19) and let $\delta\geq\min(\gamma-b_{1},b_{d+1}-\gamma)$. This yields a feasible solution to (9). In other words, in this case we can tell if $C$ is empty or not by just looking at the the ordered partition of $b$; the relative position of the midpoint is irrelevant.

Claim II: Suppose $L_{+1}\neq\emptyset=L_{-1}$. Then, (9) is feasible if and only if, in addition to (I.1), we have

[TABLE]

If (I.1), (II.1), (II.2), and (II.3) hold, take $\gamma=\max\left\{b_{i}\,\vphantom{i\in L_{0}\cup L_{-}\cup\{1\}}\right|\left.\vphantom{b_{i}}\,i\in L_{0}\cup L_{-}\cup\{1\}\right\}$ and $\delta=b_{d+1}-\gamma$.

Conversely, if $(\gamma,\delta)$ is feasible for (9) then $\gamma+\delta=b_{i}=b_{d+1}$ for all $i\in L_{+1}$. In particular, $\operatorname{bis}_{(F,G)}(a,b)$ is empty unless $\{b_{i}\,|\,i\in L_{+1}\}=\{b_{d+1}\}$ is a singleton. Since the coefficients of $b$ are in ascending order, it follows that

[TABLE]

and $b_{i}=b_{d+1}$ for all $i\in L_{+1}$. This shows that (II.1) holds. Now the constraints of (9) translate into (19) as in the previous case, which implies (I.1). Additionally

[TABLE]

Adding (20) and (21) now yields

[TABLE]

which, by (19), gives (II.2). Finally, (II.3) follows from the fact that the only possible value in $\left\{b_{i}\,\vphantom{i\in L_{0}}\right|\left.\vphantom{b_{i}}\,i\in L_{0}\right\}$ is $\gamma$.

The case where $L_{-1}\neq\emptyset$ and $L_{+1}=\emptyset$ is analogous.

Claim III: Suppose $L_{+1}\neq\emptyset\neq L_{-1}$. Then, (9) is feasible if and only if

[TABLE]

Indeed, in this case the only candidate solution for (9) is $\gamma=\mu(b)=(b_{d+1}+b_{1})/2$ and $\delta=(b_{d+1}-b_{1})/2$. This is a solution or not depending only on whether Equations (III.1) and (III.2) are satisfied. ∎

The following gives the other direction of Theorem 5: $\mathcal{F}_{\text{bis}}^{d}$ refines $\mathcal{F}_{\text{bop}}^{d}$.

Lemma 9**.**

Let $a,b,a^{\prime},b^{\prime}\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ be two pairs of points. If the bisected ordered partitions of $b-a$ and $b^{\prime}-a^{\prime}$ are not the same, then there is a pair of faces $F,G\in\mathcal{F}(\mathbb{B}^{d})$ such that exactly one of $\operatorname{bis}_{F,G}(a,b)$ or $\operatorname{bis}_{F,G}(a^{\prime},b^{\prime})$ is empty.

Proof.

As before, we assume without loss of generality that $a=a^{\prime}=0$. We know that the bisected ordered partitions of $b$ and $b^{\prime}$ are different. Our goal is to find a pair of faces $(F,G)$ that lies in one and only one of the bisectors. We do this in three cases, depending on the difference between the bisected ordered partitions

Case I: Suppose that

[TABLE]

Without loss of generality, there is an index $i$ such that $b_{i}$ is maximum and $b^{\prime}_{i}$ is not, or $b_{i}$ is minimum and $b^{\prime}_{i}$ is not. Let $F$ be the face with

[TABLE]

and let $G=-F$. This choice makes

[TABLE]

The cell $\operatorname{bis}_{F,G}(a,b)$ is not empty by lemma 7, since the following is a solution for (9):

[TABLE]

However, $\operatorname{bis}_{(F,G)}(a^{\prime},b^{\prime})$ is empty: in order for it not to be empty we would need

[TABLE]

Case II: Suppose that $b$ and $b^{\prime}$ have exactly the same maxima and minima but the ordered partitions of $b$ and $b^{\prime}$ do not coincide. That is, there is a pair of indices, $i,j\in[d+1]\setminus\{1,d+1\}$ such that $b_{i}\geq b_{j}$ but $b^{\prime}_{i}<b^{\prime}_{j}$.

We may assume that $1$ and $d+1$ are a minimum and a maximum, respectively, of both $b$ and $b^{\prime}$. Let $F$ be the face with $F_{+}=\{d+1\}$ and $F_{-}=\{1\}$. Let $G$ be the face with $G_{+}=\{i\}$, and $G_{-}=\{j\}$. Then, (8) gives

[TABLE]

Then, $\operatorname{bis}_{F,G}(a,b)$ is not empty since $\gamma=(b_{i}+b_{j})/2$, and $\delta=\operatorname{dist}(a,b)$ is a solution of (9). However, the system for $b^{\prime}$ is infeasible, since $b^{\prime}_{i}<b^{\prime}_{j}$. Thus, $\operatorname{bis}_{F,G}(a^{\prime},b^{\prime})$ is empty.

Case III: Suppose that $b$ and $b^{\prime}$ have exactly the same maxima and minima and the same ordered partitions but the midvalue does not coincide.

As before, we assume without loss of generality that $1$ and $d+1$ are a minimum and a maximum, respectively, of both $b$ and $b^{\prime}$. Then, there is an index $i\in[d+1]\setminus\{1,d+1\}$ such that $\mu(b)\leq b_{i}$ but $\mu(b^{\prime})>b^{\prime}_{i}$ (or vice-versa, but that would give an equivalent case).

In this case, we let $F$ and $G$ be the faces with $F_{+}=\{d+1\}$, $F_{-}=\{1\}$, $G_{+}=\{1\}$ and $G_{-}=\{i\}$. These faces produce

[TABLE]

Then, $\operatorname{bis}_{F,G}(a,b)$ is not empty since $\gamma=\mu(b)$, $\delta=\operatorname{dist}(a,b)/2$ is a solution of (9) for $b$

However, the system for $b^{\prime}$ is infeasible. This is because (9) specifies that $\gamma+\delta\geq b^{\prime}_{d+1}$, and $\gamma-\delta=b^{\prime}_{1}$. Adding them together and dividing by two we get $\gamma\geq(b^{\prime}_{1}+b^{\prime}_{d+1})/2=\mu(b^{\prime})$. We also need by (9) and $i\in L_{+}$ that $\gamma\leq b^{\prime}_{i}$. Then, $\mu(b)\leq b^{\prime}_{i}$, which contradicts our assumption. ∎

Corollary 8**.**

The tropical bisector of two points in general position has $\Theta(d^{4})$ maximal cells.

Proof.

The upper bound is trivial, each maximal cell corresponds to a choice of a pair of facets $(F,G)$ from the tropical ball. For the lower bound, assume without loss of generality that $a=0$ and $b_{1}<b_{2}<\dots<b_{d+1}$. Then, for each choice of $i,j,k,\ell\in\{1,\dots,d+1\}$, all different and with $\max\{j,\ell\}<\min\{i,k\}$, let $F_{+}=\{i\}$, $F_{-}=\{j\}$, $G_{+}=\{k\}$, $G_{-}=\{\ell\}$. By Claim 1 in the proof of Lemma 8 the set $\operatorname{bis}_{(F,G)}(a,b)$ is not empty. Since $a$ and $b$ are in general position and $F$ and $G$ are facets of $\mathbb{B}^{d}$, we have $\dim\operatorname{bis}_{(F,G)}(a,b)=d-1$. There are $4\binom{d+1}{4}$ ways of choosing such $\{i,j,k,\ell\}$. ∎

4.4. The structure of tropical Voronoi regions

A polytrope is an ordinary polytope which is also convex in the tropical sense (with respect to $\min$ and $\max$ simultaneously); see [17]. These are precisely the ordinary polytopes whose facets normals are roots of type $A_{d}$, i.e., $e_{i}-e_{j}$ for $i\neq j$; they generalize the “alcoved polytopes” of Lam and Postnikov [20]. Here we relax this notion by also calling a not necessarily bounded ordinary polyhedron a polytrope if its facets normals are roots of type $A_{d}$; this was called a “weighted digraph polyhedron” in [18].

The tropical unit ball $\mathbb{B}^{d}$ is a polytrope. But a more important example for us are the polytropes $Q=\bigcap_{a\in S}(a_{+}F_{a})$, where $F_{a}\in\mathbb{B}^{d}$ for each $a\in S$. Recall from Section 3 that in such a $Q$ bisectors of subsets of $S$ agree with affine subspaces. Thus:

Lemma 10**.**

For each polytrope $Q$ as above and $a\in S$, the set $Q\cap\operatorname{Vor}_{S}(a)$ is the intersection of $Q$ with ordinary affine halfspaces with facet normal $e_{i}-e_{j}-e_{k}+e_{\ell}$, where $i$ and $j$ are fixed.

Proof.

Let $i$ and $j$ be the coordinates maximized and minimized in $F_{a}$, respectively. For each $b\in S\setminus a$, the condition for $x$ to be closer to $a$ than to $b$ is that $x_{i}-a_{i}-x_{j}+a_{j}\leq x_{k}-b_{k}-x_{\ell}+b_{\ell}$, where $k$ and $\ell$ are the coordinates corresponding to $F_{b}$; see Proposition 4. ∎

We call the intersection of a (possibly unbounded) polytrope with ordinary affine halfspaces with facet normal $e_{i}-e_{j}-e_{k}+e_{\ell}$, where $i$ and $j$ are fixed, a semi-polytrope of type $(i,j)$. A semi-polytrope in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}\cong\mathbb{R}^{d}$ has at most $2\tbinom{d+1}{2}$ facets, since there are at most $(d+1)d$ vectors $e_{i}-e_{j}-e_{k}+e_{\ell}$ for $k\neq\ell$ and fixed $(i,j)$, plus the (at most) $(d+1)d$ facets of a polytrope.

A set $X\subset\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}\cong\mathbb{R}^{d}$ is star convex with center $c$ if for any point $x\in X$ the ordinary line segment $[c,x]$ is contained in $X$. Clearly any convex set is star convex, but the converse does not hold. Star convex sets are contractible. Despite the many differences to Euclidean Voronoi diagrams, the following result expresses a key similarity.

Theorem 6**.**

Let $S\subset\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ be a finite set in weak general position. Then each tropical Voronoi region of $S$ is the star convex union of finitely many (possibly unbounded) semi-polytropes.

Proof.

That Voronoi regions for polyhedral norms are star-convex is a well-known fact (see [4, p. 133] or [23, p. 127]), which follows for example from Theorem 2. By Lemma 10, $\operatorname{Vor}_{S}(a)$ decomposes as finitely many semi-polytropes, by intersecting it with the individual polyhedra $Q=\bigcap_{a\in S}(a_{+}F_{a})$, for all choices of $\{F_{a}\}_{a\in S}$. ∎

Semi-polytropes are not necessarily tropically convex, and this entails that the regions of a tropical Voronoi diagram are not necessarily tropically convex either; see Figure 8 for an example. The tropical torus $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ is compactified by the tropical projective space $\mathbb{TP}^{d}$; the latter is the $\max$-tropical convex hull of the $d+1$ $\max$-tropical unit vectors

[TABLE]

In this way, $\mathbb{TP}^{d}$ may be seen as an infinitely scaled tropical unit ball, which is a polytrope; see [18, §3.5]. Similarly for arbitrary (semi-)polytropes the line between bounded and unbounded is blurred in the compactification.

5. Computing tropical Voronoi diagrams

We will discuss several algorithms. Some of these methods are similar to their classical Euclidean counterparts, others rely on tailored data structures, which are based on Theorem 6. For the complexity analysis of our algorithms we will consider the dimension as constant.

5.1. The planar case

There are several methods for computing Euclidean Voronoi diagrams in $\mathbb{R}^{2}$ with the optimal time complexity $O(n\log n)$ and linear space; see [7, §7.2]. This agrees with the situation for planar tropical convex hull computations; see [16, §5]. Chew and Drysdale [6] gave a divide-and-conquer algorithm with the same complexity for planar Voronoi diagrams with respect to arbitrary norms. Here we sketch a tropical analog of Fortune’s beach line algorithm [10]; see also [27].

Suppose that we are given a set $S$ of $n$ sites in $\mathbb{R}^{3}/\mathbb{R}\mathbb{1}$. In view of Theorem 6 the tropical Voronoi diagram of $S$ gives rise to a planar graph where vertices are circumcenters of triples of points in $S$, edges are two point bisectors, and faces are Voronoi regions. We can make this planar embedding piecewise linear by subdividing each bisector into at most five segments; see Figure 5. The relevant data structure, as in the classical setting, is a doubly-connected edge list which requires $O(n)$ space; see [7, §2.2].

The beach line algorithm is based on a line sweep. The tropical sweep line at time $t$ in $\mathbb{R}^{3}/\mathbb{R}\mathbb{1}$ is the set $L(t)=(0,t,0)+\mathbb{R}(0,0,1)+\mathbb{R}\mathbb{1}$. Note that $L(t)$ is an ordinary line which is also tropically convex (with respect to $\min$ and $\max$). For an arbitrary point $x$ we call the set

[TABLE]

the parabola spanned by $x$ and $L(t)$; here $\operatorname{dist}(x,L(t))=\min\{\operatorname{dist}(x,y)\,|\,y-(0,t,0)\in\mathbb{R}(0,0,1)+\mathbb{R}\mathbb{1}\}$. This is a $1$-dimensional polyhedral complex, which is homeomorphic with $L(t)$ via orthogonal projection, consisting of five segments.

We will assume that our set $S$ of sites is in general position and hence, in particular, each sweep line contains at most one site. A point $a=(a_{1},a_{2},a_{3})$ is said to have been visited by the sweep line $L(t)$ if $a_{2}-a_{1}\leq t$.

The beach line $B(t)$ of $S$ at time $t$ is formed by the points $(b_{1},b_{2},b_{3})$ which lie on a parabola $P(s,t)$ for a visited point $s$ such that $b_{2}-b_{1}$ is maximal among all such points for a fixed value $b_{3}-b_{1}$. That is, the beach line is formed by the right-most points on the parabolas spanned by the visited points and the sweep line; see Figure 9. So $B(t)$ is a union of parabolic arcs; it is easy to see that each parabola contributes at most two arcs to the beach line at any time. Like a single parabola also the beach line $B(t)$ is homeomorphic to $L(t)$ via orthogonal projection. In the portion of $\mathbb{R}^{3}/\mathbb{R}\mathbb{1}$ left to $B(t)$ the tropical Voronoi diagram of $S$ is known at time $t$.

Observation 7**.**

The beach line is a polygonal line with $O(n)$ segments.

The actual algorithm works as in the classical case. We maintain a priority queue of site events (when the sweep line visits a site) and circle events (when there is a candidate for a new vertex of the tropical Voronoi diagram). The total number of events is linear in $n$. As in the classical case, it is possible to relax the condition on general position by means of symbolic perturbation.

Theorem 8**.**

The beach line algorithm computes a tropical Voronoi diagram of $n$ sites in $\mathbb{R}^{3}/\mathbb{R}\mathbb{1}$ in $O(n\log n)$ time and $O(n)$ space.

For the output we can choose, with the same complexity, between an abstract planar graph (encoding $\operatorname{Vor}(S)$ topologically) and its piecewise linear embedding resulting from Theorem 6.

5.2. Polytrope partitions

Let $S\subseteq\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ be a finite set of sites. From Theorem 6 we know that the tropical Voronoi diagram can be described in terms of (semi-)polytropes. For the definition and basic facts on polytropes, see Section 4.4 and [17]. The following takes inspiration from the trapezoid map data structure; see [7, §6].

Definition 2**.**

A polytrope partition for $S$ is a finite collection $\mathcal{C}$ of (perhaps unbounded) non-degenerate polytropes with disjoint interiors, covering $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$, such that:

1. (1)

each facet-defining hyperplane of any cell in $\mathcal{C}$ lies in the hyperplane arrangement $S+A_{d}$. 2. (2)

for each cell $P$ in $\mathcal{C}$ and site $a\in S$ the restricted Voronoi region $\operatorname{Vor}_{S}(a)\cap P$ is contained in a maximal cone $a+F$ of $a+\mathcal{F}(\mathbb{B}^{d})$.

A valid labeling for $\mathcal{C}$ assigns to each cell $P\in\mathcal{C}$ a (partial) matching $\mathcal{L}_{\mathcal{C}}(P)$ of $S\times\mathcal{F}(\mathbb{B}^{d})$ containing

[TABLE]

The distance function $x\mapsto\operatorname{dist}(x,a)$ to a fixed site $a$ is piecewise linear, as it is linear in each translated cone $a+F$ for $F\in\mathcal{F}(\mathbb{B}^{d})$. Polytrope partitions are designed to exploit this fact, so that the distances to the relevant sites in each cell are linear. In particular:

Observation 9**.**

Let $P$ be a cell in a polytrope partition $\mathcal{C}$ for $S$. Then for all $x\in P$ we have

[TABLE]

where $\lambda_{F_{a}}$ is the linear function defined by restricting the distance to $a$ on some maximal cone $a+F_{a}$ of $a+\mathcal{F}(\mathbb{B}^{d})$ which contains $\operatorname{Vor}_{S}(a)\cap P$. Thus computing the restriction of $\operatorname{Vor}(S)$ to the polytrope $P$ amounts to finding the regions of linearity of the tropical polynomial $\min_{a\in S}\lambda_{F_{a}}(x)$. The latter can be obtained via an ordinary dual convex hull computation.

Note that the maximal cone $F_{a}$ in the above is irrelevant if $\operatorname{Vor}_{S}(a)\cap P=\emptyset$, and it is unique, by axiom (2), if $\operatorname{Vor}_{S}(a)\cap P\neq\emptyset$. If the polytrope partition is equipped with a valid labeling, then this tells us the choice of the right cone $F_{a}$ for each site, if it exists.

The following shows that polytrope partitions exist.

Example 6**.**

The braid arrangement $A_{d}$ consists of the $\tbinom{d+1}{2}$ ordinary hyperplanes $\{x\,|\,x_{i}=x_{j}\}$, where $i\neq j$. This gives rise to the standard polytrope partition $S+A_{d}$, which is finer than any other polytrope partition for $S$; Figure 8 shows an example for $d=2$. This construction occurs in planar tropical convex hull algorithms; see [16, Figure 3].

And finally, the following lemma shows that valid labelings exist for every polytrope partition, if the sites $S$ are in weak general position:

Lemma 11**.**

Let $\mathcal{C}$ be a polytrope partition for $S$. If $S$ is in weak general position then there is a valid labeling of $\mathcal{C}$. Moreover, for $d$ considered constant, a labeling of each polytropal cell has constant size, and it can be computed in $O(n)$ time.

Proof.

Suppose that a valid labeling does not exist for some cell $P$. Either the set of pairs in (2) matches two cones with the same site, or matches two sites with the same cone of $\mathcal{F}(\mathbb{B}^{d})$. The former cannot happen since $P$ is full-dimensional and the cones in $\mathcal{F}(\mathbb{B}^{d})$ only intersect in lower-dimensional polyhedral cones.

Then there are sites $a,b\in S$ and a maximal cone $F\in\mathcal{F}(\mathbb{B}^{d})$ such that the sets $(a+F)\cap P\cap\operatorname{Vor}_{S}(s)$ and $(b+F)\cap P\cap\operatorname{Vor}_{S}(t)$ both are non-empty.

With the notation of (22) we have $\lambda_{F_{a}}=\lambda_{F_{b}}$; and we shortly write $\lambda$. Since the sites are in weak general position, we may assume that $\lambda(b)>\lambda(a)$. Picking $y\in(b+F)\cap P\cap\operatorname{Vor}_{S}(b)$ yields

[TABLE]

where the last inequality follows from (22). The resulting inequality $\operatorname{dist}(y,b)>\operatorname{dist}(y,S)$ implies that $y\not\in\operatorname{Vor}_{S}(b)$, which is a contradiction. Hence a valid labeling does exist.

To compute such labeling, we iterate through all the sites. For each site $a$, the candidate facet of $F_{a}\in\mathcal{F}(\mathbb{B}^{d})$ is known by definition of the polytrope partition. To check if $(a,F_{a})$ is a labeling candidate, we need to determine if $(a+F_{a})\cap P\cap\operatorname{Vor}_{S}(a)$ is empty or not. This amounts to solving a linear program that has constant size (as $d$ is a constant). It follows that the entire labeling can be computed in $O(n)$ time. ∎

We aim at a first algorithm for computing a tropical Voronoi diagram in arbitrary dimension. This will employ the standard polytrope partition from Example 6.

Lemma 12**.**

If $S$ is in weak general position and has size $n$ then the trivial polytrope partition has

[TABLE]

maximal cells, if we consider $d$ a fixed constant.

Proof.

Pick a generic direction $v\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$. The cells of the polytrope partition that are bounded in the direction of $v$ are in a one-to-one correspondence with the vertices of the arrangement, by associating each polytrope with the optimum of the linear program maximizing $v^{T}x$. Since the number of vertices equals

[TABLE]

by Cayley’s formula, it suffices to show that the number of unbounded cells is in $O(n^{d-1})$, as the bases of $A_{d}$ correspond to spanning trees in the complete graph with $d$ vertices. The unbounded cells intersect a hyperplane, $H$, normal to $v$ that is far enough in the $v$ direction. The cells intersecting $H$ are the same as the cells in the restricted hyperplane arrangement, which is a $(d-1)$-dimensional arrangement with $N={d+1\choose 2}n$ hyperplanes. The number of such cells is known to be in $O\left(N^{d-1}\right)$, which agrees with $O\left(n^{d-1}\right)$ as $N$ depends linearly on $n$. ∎

Remark 9** (Standard polytrope partition algorithm).**

This directly yields a first algorithm for computing a tropical Voronoi diagram of $n$ sites in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ in $O(n^{d+1})$ time, as follows: First, we sort $S$ along each of the ${(d+1)\choose 2}$ directions $e_{i}+e_{j}$, in $O(n\log n)$ time. As in the proof of Lemma 12 we pick a generic direction $v\in\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$. We can compute the vertices of the hyperplane arrangement $A_{d}+S$ in time $O(n^{d})$ by enumerating all $d$-sets of independent directions, which can be derived from the the oriented spanning tree of $K_{d+1}$, in constant time. For each of the $d$ directions we choose an index $i\in[n]$.

Next we perturb each such vertex $p$ by a small multiple of $-v$, and we collect the intersection of bands of contiguous parallel hyperplanes of $A_{d}+S$ that contain the perturbed point. This can be done in time $O(\log n)$ for each direction. In this way, we find those cells which are bounded in the direction of this particular $v$ in linear output-dependent time. We repeat the same procedure for a set of directions $v_{1},\dots,v_{d+1}$ which positively span the entire space $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$. Each polytropal cell will be bounded in at least one of these directions, and thus their enumeration is still in $O(n^{d})$ for $d$ fixed.

Then, for each polytrope $P$, we compute a corresponding labeling in time $O(n)$, by Lemma 11. Therefore, we can compute the standard polytrope partition, including labels, in time $O(n^{d+1})$ for fixed $d$. The tropical Voronoi diagram in each cell is an ordinary dual convex hull problem of constant size. This computation splits each polytrope in the partition into semi-polytropes. The convex hull problem can be solved in constant time, and hence this algorithm takes $O(n^{d+1})$ time, if $d$ is considered a fixed constant.

We implemented this procedure in polymake [11], version 3.6.

Question 1**.**

In the plane $\mathbb{R}^{3}/\mathbb{R}\mathbb{1}$, we believe that ideas similar to the “trapezoidal maps” used in point location, see [7, §6.1], should yield polytrope partitions of linear size but we did not work out the details. More generally: Is there a polytrope partition of complexity better than $\Theta(n^{d})$ in arbitrary dimension? One could hope for something in $O(n^{d/2})$, which is the worst-case complexity of Euclidean Voronoi diagrams.

5.3. An $O(n^{d}\log n)$ randomized incremental algorithm in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$

We can improve the algorithm from Remark 9 by constructing a polytrope partition incrementally. The idea is to update an existing polytrope partition by including a new point and to employ randomization to improve the efficiency. Moreover, we will also produce a coarser polytrope partition than the standard one, but only by a constant in $d$.

A key ingredient is a new data structure that we call a polytrope tree. Throughout we assume that the set $S$ of $n$ sites forms a subset of $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ in general position. We fix the polytrope partition $\mathcal{C}:=S+\mathcal{F}(\mathbb{B}^{d})$, which is coarser than the standard polytrope partition but only by a factor which is constant in $d$; see Example 6.

Definition 3**.**

A polytrope tree for $S$ is a (rooted) tree $T$ such that

1. (1)

for each leaf $\ell$ there is a polytropal cell $P(\ell)$ of $\mathcal{C}$; 2. (2)

for each interior node $i$ there is a site $a(i)\in S$ and a polytrope $P(i)$.

These satisfy the following consistency conditions:

• •

for the root node $r$ of $T$ we have $P(r)=\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$, which may be seen as an unbounded polytrope;

• •

the map $\ell\mapsto P(\ell)$ is a bijection between the set of leaves of $T$ and the set of polytropes in $\mathcal{C}$;

• •

the map $i\mapsto a(i)$ is a surjection from the interior nodes onto the set $S$;

• •

if $i$ is an interior node with children $c_{1},\dots,c_{k}$, then $P(c_{1}),\dots,P(c_{k})$ form the maximal cells of $(a(i)+\mathcal{F}(\mathbb{B}^{d}))\cap P(i)$.

It is easy to construct a polytrope tree for $S$, and its purpose is to speed up the computation of a valid labeling. This will reduce the algorithmic complexity from $O(n^{d+1})$ to $O(n^{d}\log n)$.

For the incremental update to insert a new site $b\not\in S$ we maintain a stack $\Sigma$ of unvisited nodes in a given polytrope tree for $\Sigma$ and process it as follows:

• •

the stack $\Sigma$ is initialized with the root node $r$;

• •

we remove the top node $q$ from the stack $\Sigma$ unless it is empty;

• •

if $q$ is an interior node such that $P(q)$ intersects more than one maximal cone of $s+\mathcal{F}(\mathbb{B}^{d})$, then we push the children of $q$ onto the stack $\Sigma$;

• •

if $q$ is a leaf such that $P(\ell)$ intersects more than one maximal cone of $p+\mathcal{F}(\mathbb{B}^{d})$, then we create the intersections of $P(\ell)$ with $s+\mathcal{F}(\mathbb{B}^{d})$ as new leaves, which now become children of $q$, and we set $a(q)\leftarrow b$.

Note that an interior node $q$ with $P(q)$ contained in a unique maximal cone of $a+\mathcal{F}(\mathbb{B}^{d})$ is kept unchanged, and its children will not be visited. The following is the essential part of the complexity analysis.

Proposition 7**.**

Let $T$ be a polytrope tree created in the way explained above, where the $n$ sites in $S$ are processed in uniformly random order. Then the expected height of $T$ is of order $O(\log n)$, if $d$ is considered a fixed constant.

Proof.

Let $P$ be a polytrope in the polytrope partition $\mathcal{C}(S)$. For each ordering $\pi:[n]\to S$ of $S$ we have a polytrope tree $T(S,\pi)$ with $P$ as a leaf. By induction on $n$ we will show:

[TABLE]

where the expectation $E[\cdot]$ is taken uniformly over all $n!$ orderings of $S$, and $h_{T(S,\pi)}(P)$ is the depth of the leaf of $P$ in $T(S,\pi)$.

We proceed by backwards analysis. Let $S^{\prime}\subset S$ be the subset of sites that lie in some facet-defining hyperplane of $P$. Since $P$ has at most $d(d+1)$ facets and (by general position) their corresponding hyperplanes contain each exactly one point of $S$, we have $|S^{\prime}|\leq d(d+1)$. Thus, the probability that the height $h_{T}(P)$ increases in the last insertion is at most $d(d+1)/n$. Since the increase is by exactly one, we have

[TABLE]

where $P^{\prime}$ is the polytrope containing $P$ in the polytrope partition before the last insertion. By induction hypothesis

[TABLE]

The last two formulas give Eq. (23). ∎

Corollary 9**.**

The above method constructs a polytrope tree for the polytrope partition $S+\mathcal{F}(\mathbb{B}^{d})$ in expected time $O(n^{d}\log n)$ and space $O(n^{d})$, for $d$ constant.

Proof.

The algorithm that inserts a new site $a$ into the tree only visits nodes that are above some leaf requiring an update. For each such leaf $\ell$ the polytrope $P(\ell)$ intersects one of the $d(d+1)$ hyperplanes in $a+A_{d}$. This implies that there are $O(n^{d-1})$ of them. Since the expected depth of every leaf is $O(\log n)$ it requires expected time $O(n^{d-1}\log n)$ for inserting $a$. Hence the total complexity for $n$ sites amounts to $O(n^{d}\log n)$. ∎

In order to compute the tropical Voronoi diagram, we also need to compute the labeling of this polytrope partition. The naive way is to compute the labeling for each leaf as we did in Remark 9.

A slight improvement is to compute the labeling during the depth-first-search (DFS) exploration of the tree at each insertion of a new site. But in this way, even if an interior node is completely contained in only one cone of the fan $a+\mathcal{F}(\mathbb{B}^{d})$, we need to descend to its subtree in order to update the labels. This would slow the algorithm down to $\Theta(n^{d+1})$ because each insertion will have to iterate through all the leaves.

A better way is to compute the labeling lazily. To this end we equip each interior node $i$ with a partial labeling $\mathcal{L}_{\mathcal{C}}(i)$. With each new insertion, we proceed as we just explained, but we do not cascade down the label updates. Only once all sites in $S$ have been inserted we cascade the partial labels, updating them in DFS order. This takes $O(n^{d}\log n)$ time to compute the polytropes and the lazy labelings, plus $O(n^{d})$ time to cascade the lazy labelings down in the tree, for a total time complexity of $O(n^{d}\log n)$ time. This gives our final result.

Theorem 10**.**

There is a randomized incremental algorithm for computing tropical Voronoi diagrams of $n$ sites in $\mathbb{R}^{d+1}/\mathbb{R}\mathbb{1}$ in general position with expected time complexity $O(n^{d}\log n)$ and space complexity $O(n^{d})$, for $d$ constant.