Abstract 1 Introduction 2 Color Voronoi diagrams and arrangements 3 The colorful Clarkson–Shor framework 4 Color Voronoi diagrams under general distance functions 5 Iterative algorithms for color Voronoi diagrams References

Higher-Order Color Voronoi Diagrams and the Colorful Clarkson–Shor Framework

Sang Won Bae ORCID Division of AI Computer Science and Engineering, Kyonggi University, Suwon, South Korea Nicolau Oliver ORCID Faculty of Informatics, Università della Svizzera italiana, Lugano, Switzerland Evanthia Papadopoulou ORCID Faculty of Informatics, Università della Svizzera italiana, Lugano, Switzerland
Abstract

Given a set S of n colored sites, each sS associated with a distance-to-site function δs:2, we consider two distance-to-color functions for each color: one takes the minimum of δs for sites sS in that color and the other takes the maximum. These two sets of distance functions induce two families of higher-order Voronoi diagrams for colors in the plane, namely, the minimal and maximal order-k color Voronoi diagrams, which include various well-studied Voronoi diagrams as special cases. In this paper, we derive an exact upper bound 4k(nk)2n on the total number of vertices in both the minimal and maximal order-k color diagrams for a wide class of distance functions δs that satisfy certain conditions, including the case of point sites S under convex distance functions and the Lp metric for any 1p. For the L1 (or, L) metric, and other convex polygonal metrics, we show that the order-k minimal diagram of point sites has O(min{k(nk),(nk)2}) complexity, while its maximal counterpart has O(min{k(nk),k2}) complexity. To obtain these combinatorial results, we extend the Clarkson–Shor framework to colored objects, and demonstrate its application to several fundamental geometric structures, including higher-order color Voronoi diagrams, colored j-facets, and levels in the arrangements of piecewise linear/algebraic curves/surfaces. We also present iterative algorithms to compute higher-order color Voronoi diagrams.

Keywords and phrases:
higher-order Voronoi diagrams, color Voronoi diagrams, Hausdorff Voronoi diagrams, colored j-facets, arrangements, Clarkson–Shor technique
Funding:
Sang Won Bae: Supported by the National Research Foundation of Korea(NRF) grant funded by the Korea government(MSIT) (No. RS-2023-00251168), and in part by the University of Bayreuth Centre of International Excellence “Alexander von Humboldt” (Senior Fellowship 2023).
Evanthia Papadopoulou: Supported in part by the Swiss National Science Foundation (SNF), project 200021E-201356.
Copyright and License:
[Uncaptioned image] © Sang Won Bae, Nicolau Oliver, and Evanthia Papadopoulou; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Computational geometry
Related Version:
Full Version: https://arxiv.org/pdf/2504.06960 [10]
Acknowledgements:
The authors would like to thank Otfried Cheong, Christian Knauer, and Fabian Stehn for valuable discussions and comments, in particular, at the beginning of this work. The research by the first author has been done partly during his visits to Universität Bayreuth, Bayreuth, Germany and to Università della Svizzera italiana, Lugano, Switzerland.
Editors:
Oswin Aichholzer and Haitao Wang

1 Introduction

Let S be a set of n sites, each of which is assigned a color from a set K={1,,m} of m colors. Let SiS be the set of sites in color iK. We consider two distance-to-color functions from any point x2 to each color iK:

di(x):=minsSiδs(x) and d¯i(x):=maxsSiδs(x),

where δs(x) denotes the prescribed distance-to-site function from x2 to site sS.

Based on the two sets of point-to-color distances, we can define two different Voronoi diagrams of m colors in K. For each 1km1 and subset HK with |H|=k, define

Rk(H;S) :={x2di(x)<dj(x) for all iH and jKH}
R¯k(H;S) :={x2d¯i(x)>d¯j(x) for all iH and jKH},

called the minimal and maximal color Voronoi regions of H with respect to S. We then define the order-k minimal color Voronoi diagram of S, 𝖢𝖵𝖣k(S), to be the partition of 2 into the minimal regions Rk(H;S), for all HK with |H|=k; the order-k maximal color Voronoi diagram of S, 𝖢𝖵𝖣¯k(S), to be the partition of 2 into the maximal regions R¯k(H;S). In other words, 𝖢𝖵𝖣k(S) partitions 2 by k nearest colors under the minimal distances {di}iK, while 𝖢𝖵𝖣¯k(S) partitions 2 by k farthest colors under the maximal distances {d¯i}iK.

Figure 1: Special cases of 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) for colored points S in the Euclidean plane.

These two families of color Voronoi diagrams generalize various conventional counterparts.

  • For k=1, 𝖢𝖵𝖣1(S) and 𝖢𝖵𝖣¯1(S) correspond to the ordinary nearest-site and farthest-site Voronoi diagrams of S, 𝖵𝖣(S) and 𝖥𝖵𝖣(S), respectively, where adjacent faces of a common color are merged. See Figure 1(a) and (b).

  • If m=n, that is, each site in S has a distinct color, then 𝖢𝖵𝖣k(S)=𝖵𝖣k(S), the ordinary order-k Voronoi diagram without colors [29, 30, 13, 40]. In this case, it holds that 𝖢𝖵𝖣¯k(S)=𝖵𝖣nk(S), also known as the order-k farthest-site Voronoi diagram or the order-k maximal Voronoi diagram [5].

  • The farthest color Voronoi diagram 𝖥𝖢𝖵𝖣(S) [24, 1, 32, 9] partitions the plane 2 into regions of colors iK that consist of all points p2 whose farthest color is i with respect to the minimal distances {di}iK. The Hausdorff Voronoi diagram 𝖧𝖵𝖣(S), also known as the min-max or cluster Voronoi diagram [22, 37, 38], similarly partitions 2 into regions of nearest colors with respect to the maximal distances {d¯i}iK. We thus have 𝖥𝖢𝖵𝖣(S)=𝖢𝖵𝖣m1(S) and 𝖧𝖵𝖣(S)=𝖢𝖵𝖣¯m1(S). These two diagrams are often refined so that the region of each color iK is subdivided into subregions of sites in Si that determine the distance-to-color di or d¯i. We denote these refined versions by 𝖥𝖢𝖵𝖣(S) and 𝖧𝖵𝖣(S), respectively. See Figure 1(c) and (d).

In this paper, we initiate the study of higher-order color Voronoi diagrams, 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S), with general distance functions δs for sS. Our main results are as follows:

  1. (1)

    We prove an exact upper bound 4k(nk)2n on the total number of vertices of both order-k diagrams 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) for each 1km1, when the sites sS are points and δs(x) is given as the convex distance from x2 to s based on any convex and compact subset of 2, which includes the Lp distance for any 1p. This result is in fact a corollary of our more general result: we derive an exact equation under certain conditions on functions δs, which only requires relations on the numbers of vertices and unbounded edges in 𝖵𝖣(S) and 𝖥𝖵𝖣(S) for SS. The bound 4k(nk)2n can be realized, for example, when m=n, so it is tight and best possible.

  2. (2)

    Under the L1 or the L metric, we prove that 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) consist of at most min{4k(nk)2n,4(nk)2} and min{4k(nk)2n,2k2} vertices, respectively. Similar bounds are derived for any convex distance function based on a convex polygon.

  3. (3)

    We present an iterative algorithm that computes color Voronoi diagrams of order 1 to k in O(k2n+nlogn) expected or O(k2nlogn) worst-case time, provided that S and {δs}sS satisfy the requirements of abstract Voronoi diagrams [28] and an additional condition (see Section 5), which includes colored points S under any smooth convex distance function. For points S in the Euclidean plane, it can be reduced to O(k2n+nlogn) worst-case time.

Our combinatorial results generalize previously known bounds for the ordinary higher-order Voronoi diagrams 𝖵𝖣k(S) of uncolored sites S, which is a special case of m=n in our setting. The asymptotically tight bound O(k(nk)) on the complexity of 𝖵𝖣k(S) has been proved not only for points [29, 30] under the Lp metric but also for line segments [40] and even in the abstract setting [13]. Under the L1 (or L) metric, the complexity of 𝖵𝖣k(S) is known to be O(min{k(nk),(nk)2}) for a set S of points or line segments [30, 40].

In particular, the same exact number 4k(nk)2n can be derived for the ordinary order-k diagram 𝖵𝖣k(S) and order-(nk) diagram 𝖵𝖣nk(S) under the Euclidean metric from previous results [21, 29, 20]. In his book [21, Chapter 13], Edelsbrunner showed that the number of vertices of 𝖵𝖣k(S) is exactly (4k2)n2k2ek2i=1k1ei, if S is in general position, where ek denotes the number of unbounded edges in 𝖵𝖣k(S) (or, equivalently, k-sets in S). One can easily verify that the total number of vertices in 𝖵𝖣k(S) and 𝖵𝖣nk(S) is exactly 4k(nk)2n by using the identities: ek=enk and k=1n1ek=n(n1). This can also be verified from the inductive approach of Lee [29]. As a recent result related to ours, Biswas et al. [12] derived exact relations for the complexity of 3D Euclidean higher-order Voronoi diagrams with Morse theory, extending the inductive argument by Lee.

Another remarkable special case of our results is when k=m1, which yields the O(m(nm+1)) bound for the farthest color Voronoi diagram 𝖥𝖢𝖵𝖣(S)=𝖢𝖵𝖣m1(S) and the Hausdorff Voronoi diagram 𝖧𝖵𝖣(S)=𝖢𝖵𝖣¯m1(S). The worst-case complexity of 𝖥𝖢𝖵𝖣(S) and 𝖧𝖵𝖣(S) is known to be Θ(mn) or Θ(n2), if S is a set of points or line segments under the Euclidean or L1 metric [22, 24, 1, 38, 9]. While the upper bounds O(mn) and O(n2) are shown to be tight by matching lower bound constructions [24, 1, 22], they become significantly loose when n is close to m; if n=m, we have 𝖥𝖢𝖵𝖣(S)=𝖥𝖵𝖣(S) and 𝖧𝖵𝖣(S)=𝖵𝖣(S), where both have linear complexity. Hence, our new results not only prove tight upper bounds simultaneously on both diagrams, but also formally reveal a smooth extension upon the previous knowledge about the ordinary Voronoi diagrams for any mn. Prior to our results, it was known that the complexity of 𝖥𝖢𝖵𝖣(S) and 𝖧𝖵𝖣(S) can range from Θ(n) to Θ(m(nm+1)) expressed in terms of geometric parameters, called straddles for 𝖥𝖢𝖵𝖣(S) [33] and crossings for 𝖧𝖵𝖣(S) [38]. Conditions under which these diagrams have linear complexity have also been discussed [38, 33, 8].

Our combinatorial results are based on a color-augmented extension of the Clarkson–Shor framework [20], so-called the colorful Clarkson–Shor framework, which has its own interest with various applications. Our notions of colored objects and configurations are naturally inherited from any set system that fits in the original Clarkson–Shor framework, and yield a systematic scheme to deal with objects that are collections of primitive elements. Our new framework provides a unifying approach to derive the complexity of higher-order color Voronoi diagrams under general distances δs, including the well-studied diagrams 𝖥𝖢𝖵𝖣(S), 𝖧𝖵𝖣(S), and the ordinary higher-order diagram 𝖵𝖣k(S). In deriving our results, we make use of a close relation between the color diagrams 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) and the colored k-facets of S. An analogous relation for the Euclidean ordinary case (without colors) has been shown by Clarkson and Shor [20] and Edelsbrunner [21]. We also derive lower and upper bounds on the number of colored k-facets in 2. We also demonstrate more applications of the colorful Clarkson–Shor framework that result in several new bounds on levels of arrangements of piecewise linear or algebraic curves and surfaces; see [10].

Our algorithm follows the principle of iteratively computing all order-i Voronoi diagrams for 1ik, and compares to the O(k2nlogn)-time counterpart of Lee [29] for computing 𝖵𝖣1(S),,𝖵𝖣k(S) for points in the Euclidean metric. After a series of improvements [3, 7, 36, 2, 16, 41, 18], the first optimal O(nlogn+kn)-time algorithm that computes 𝖵𝖣k(S) for points S in the Euclidean plane was eventually presented in 2024 by Chan et al. [17]. Efficient algorithms of different approaches are also known for computing 𝖵𝖣k(S) of line segments S [40] or under the model of abstract Voronoi diagrams [14, 15]; and for computing 𝖥𝖢𝖵𝖣(S) [24, 1, 9, 33] and 𝖧𝖵𝖣(S) [22, 38, 4].

Due to space limit, proofs are omitted but can be found in the full version [10].

2 Color Voronoi diagrams and arrangements

Let S be a set of n sites, which can be any abstract objects, and δs:2 for sS be a given continuous function. We assume that the functions δs are in general position, similarly to [26]. More precisely, let γs be the graph of δs, that is, the xy-monotone surface {(p,δs(p))p2} in 3, and Γ:={γssS}. By the general position of functions δs, we mean the following: no more than three surfaces in Γ meet at a common point, no more than two surfaces in Γ meet at a one-dimensional curve, no two surfaces in Γ are tangent to each other, and none of the surfaces in Γ is tangent to the intersection curve of two others in Γ.

We denote the minimization diagram of Γ by 𝖵𝖣(S), the nearest-site Voronoi diagram of S, and the maximization diagram of Γ by 𝖥𝖵𝖣(S), the farthest-site Voronoi diagram of S. It is well known that the ordinary (uncolored) higher-order Voronoi diagrams 𝖵𝖣k(S) of S are determined by levels of the arrangement 𝒜(Γ) of Γ as established in earlier work [23, 5]. In the following, we discuss an analogous relation between order-k color Voronoi diagrams and levels of certain surfaces in 3.

Let us assume that the sites in S are colored with m colors in K={1,,m}. Let κ:SK denote a color assignment such that the color of sS is κ(s)K. Let Si:={sSκ(s)=i} and Γi:={γssSi} for iK. Define Ei and E¯i to be the lower and upper envelopes of surfaces in Γi, respectively, and consider the two arrangements 𝒜Γ=𝒜({E1,,Em}) and 𝒜¯Γ=𝒜({E¯1,,E¯m}). Note that Ei is the graph of the minimal distance function di of color iK, and E¯i is the graph of the maximal distance d¯i. We then consider the levels of the arrangements 𝒜Γ and 𝒜¯Γ. For 1km, let Lk be the k-level of 𝒜Γ from below and L¯k be the k-level of 𝒜¯Γ from above. So, L1 is the lower envelope of E1,,Em, and L¯1 is the upper envelope of E¯1,,E¯m. Thus, projecting L1 and L¯1 down onto 2 yields 𝖵𝖣(S) and 𝖥𝖵𝖣(S), respectively. On the other hand, Lm is the upper envelope of the lower envelopes E1,,Em, and L¯m is the lower envelope of the upper envelopes E¯1,,E¯m. Projecting Lm and L¯m down onto 2 yields the refined diagrams 𝖥𝖢𝖵𝖣(S) and 𝖧𝖵𝖣(S), respectively [24, 22].

Figure 2: The minimal color Voronoi diagrams 𝖢𝖵𝖣k(S) and the refined diagrams 𝖢𝖵𝖣k(S).
Figure 3: The maximal color diagrams 𝖢𝖵𝖣¯k(S) and the refined diagrams 𝖢𝖵𝖣¯k(S).

From the viewpoint of [23, 5], we observe the following.

  • The order-k color Voronoi diagrams, 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S), for 1km1 are the projections of LkLk+1 and of L¯kL¯k+1, respectively, onto 2.

  • For each 1km, let 𝖢𝖵𝖣k(S) denote the planar map obtained by projecting Lk down onto 2; analogously, let 𝖢𝖵𝖣¯k(S) be the map obtained by projecting L¯k onto 2. By definition, 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) refine 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S), respectively. Each face f of 𝖢𝖵𝖣k(S) (or of 𝖢𝖵𝖣¯k(S)) is associated with a site sSi such that, for any point pf, di(p)=δs(p) and iK is the k-th nearest color from p with respect to {di}iK (resp. d¯i(p)=δs(p) and iK is the k-th farthest color from p with respect to {d¯i}iK). This way, the refined diagrams 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) partition the plane 2 by the k-th nearest and k-th farthest colors, respectively.

From the construction, it is clear that 𝖢𝖵𝖣1(S)=𝖵𝖣(S) and 𝖢𝖵𝖣¯1(S)=𝖥𝖵𝖣(S). Note also that 𝖢𝖵𝖣m(S)=𝖥𝖢𝖵𝖣(S) and 𝖢𝖵𝖣¯m(S)=𝖧𝖵𝖣(S), whereas 𝖥𝖢𝖵𝖣(S)=𝖢𝖵𝖣m1(S) and 𝖧𝖵𝖣(S)=𝖢𝖵𝖣¯m1(S).

Figures 2 and 3 illustrate an example under the Euclidean metric, where S consists of n=9 points and m=4 colors: S1={s1,s2,s5}, S2={s3,s9}, S3={s4,s6,s8}, and S4={s7}, in red, blue, purple, and black, respectively; selected regions of 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) are labeled in (a)–(c); faces associated with s6S3 and those with s9S2 in 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) are shaded in purple and blue, respectively, in (d)–(g).

Now, consider the vertices and edges of the arrangements 𝒜Γ and 𝒜¯Γ. By the general position assumption, each vertex v of 𝒜Γ or of 𝒜¯Γ is determined by exactly three sites s,s,s′′S in such a way that v is a common intersection of three surfaces γs,γs,γs′′Γ. Such a vertex v is called c-chromatic for c{1,2,3} if |{κ(s),κ(s),κ(s′′)}|=c. (Note that any 1-chromatic vertex is a vertex of some single envelope Ei or E¯i.) Similarly, each edge of 𝒜Γ and of 𝒜¯Γ is determined by exactly two sites in S, and is either 1- or 2-chromatic according to the number of involved colors. We identify each vertex or edge of 𝖢𝖵𝖣k(S) or of 𝖢𝖵𝖣¯k(S) by its original lifted copy in 𝒜Γ or in 𝒜¯Γ. Observe that any c-chromatic vertex or edge of 𝒜Γ or of 𝒜¯Γ appears in c consecutive levels, as it lies in the intersection of c surfaces from {Ei}iK or from {E¯i}iK, respectively. Thus, c-chromatic vertices appear in c consecutive orders of the refined diagrams. We call a vertex or an edge of 𝖢𝖵𝖣k(S) or of 𝖢𝖵𝖣¯k(S) new if it does not appear in 𝖢𝖵𝖣k1(S) or in 𝖢𝖵𝖣¯k1(S), respectively; and old, otherwise. By definition, the vertices and edges of 𝖢𝖵𝖣1(S) and 𝖢𝖵𝖣¯1(S) are all new. Note that every edge of 𝖢𝖵𝖣k(S) (or, of 𝖢𝖵𝖣¯k(S)) is 2-chromatic and new, and appears both in 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣k+1(S) (in 𝖢𝖵𝖣¯k(S) and 𝖢𝖵𝖣¯k+1(S), resp.), being first new and then old. See Figures 2 and 3, where new 2-chromatic edges are in black, old 2-chromatic edges in gray, 1-chromatic edges in their own color, and new vertices marked by small squares.

We define vc,j=vc,j(S) for 1c3 and 0jm1 to be the number of c-chromatic vertices in 𝒜Γ below which there are exactly j surfaces from {Ei}iK; v¯c,j=v¯c,j(S) to be the number of c-chromatic vertices in 𝒜¯Γ above which there are exactly j surfaces from {E¯i}iK.

Lemma 1.

For any 1c3 and 1km, the following hold:

  1. (i)

    The number of new c-chromatic vertices of 𝖢𝖵𝖣k(S) is exactly vc,k1, and the number of new c-chromatic vertices of 𝖢𝖵𝖣¯k(S) is exactly v¯c,k1.

  2. (ii)

    The number of vertices of 𝖢𝖵𝖣k(S) is exactly v3,k1+v3,k2+v3,k3+v2,k1+v2,k2+v1,k1, and the number of vertices of 𝖢𝖵𝖣¯k(S) is exactly v¯3,k1+v¯3,k2+v¯3,k3+v¯2,k1+v¯2,k2+v¯1,k1, where vc,j=v¯c,j=0 for j<0.

  3. (iii)

    The number of vertices of 𝖢𝖵𝖣k(S) is exactly v3,k1+v3,k2+v2,k1, and the number of vertices of 𝖢𝖵𝖣¯k(S) is exactly v¯3,k1+v¯3,k2+v¯2,k1.

3 The colorful Clarkson–Shor framework

The Clarkson–Shor technique [20] is based on a general framework dealing with so-called configurations or ranges defined by a set of objects. Specifically, the following three ingredients are supposed to be given with a constant integer parameter d1, see also Sharir [42]:

  • A set S of n objects.

  • A set (S) of configurations, each of which is defined by exactly d objects in S.

  • A conflict relation χS×(S) between objects sS and configurations f(S) with the requirement that none of the d objects defining f are in conflict with f.

In the original framework, the number of objects that define a configuration does not have to be exactly d, but at most d. This restriction can be achieved by adding dummy objects to S.

Let us call such a triplet (S,(S),χ) a CS-structure. Given any CS-structure (S,(S),χ) with parameter d, we now impose a color assignment κ:SK to the objects in S, where K={1,2,,m} denotes the set of m colors with mn. For each color iK, let Si:={sSκ(s)=i}. For f(S) and set DfS of d objects defining f, κ(Df)={κ(s)sDf} is called a set of colors defining f. We build a color-to-configuration conflict relation χκK×(S) such that a color iK is in conflict with a configuration f(S) if an object sSi is in conflict with f, that is, (i,f)χκ if and only if (s,f)χ for some sSi. We are then interested in those configurations f(S) such that none of its defining colors in κ(Df) are in conflict with f. Let (S,κ)(S) be the set of these configurations, called colored configurations with respect to κ. We call f(S,κ) c-chromatic if |κ(Df)|=c for 1cd, and let the weight of f be the number of colors in K that are in conflict with f. Let c,j(S,κ)(S,κ) be the set of c-chromatic weight-j colored configurations in (S,κ).

Considering the colors in K as new objects, each of which is a collection of objects in S, observe that this color-augmented structure (K,jc,j(S,κ),χκ) for each 1cd is again a CS-structure with parameter c. Therefore, the main lemma of Clarkson and Shor [20, Lemma 2.1] automatically implies the following.

Lemma 2.

With the above notations, let 1cd, r0, and a0 be integers, and RK be a random subset of r colors. Then,

(mr)𝐄[|c,a(SR,κR)|]j=0mc|c,j(S,κ)|(ja)(mcjrca),

where SR=iRSi and κR:SRR denotes the restriction of κ to SR. The equality holds if each configuration in (S,κ) is defined by a unique set of d objects in S.

In probabilistic arguments dealing with CS-structures, it is usually necessary to have an upper bound on the number of weight-0 configurations. For any subset SS, let 0(S)(S) be the set of (uncolored) configurations f of weight 0, that is, (s,f)χ for all sS. Let T0(n) be any nondecreasing function with T0(0)=0 that upper bounds the number |0(S)| of these configurations for any set S of n uncolored objects. The following is obvious by definition.

Lemma 3.

Let SS be any subset and κ be any color assignment for S. Then,

c=1d|c,0(S,κ)|=|0(S)|T0(|S|).

In many applications, such an upper bound function T0 is either known from previous work or relatively easy to obtain. Once we have T0, we can derive general upper bounds on the number of corresponding colored configurations of weight at most k in such a procedural way as done for uncolored cases [43, 20, 34].

Theorem 4.

With the above notations, suppose T0 is a convex function. For each 1cd, the total number of (c)-chromatic colored configurations is bounded by

b=1cj=0mb(mbjcb)|b,j(S,κ)|(mc)1miKT0(c|Si|)=O(mc1T0(cn)).

Also, for each 2cd and 0kmc1, it holds that

j=0k|c,j(S,κ)|=O((k+1)cmiKT0(m|Si|k+1))=O((k+1)cmT0(mnk+1)).

Remark that the left-hand side of the first bound can be seen as a “weighted” count of (c)-chromatic colored configurations, and that the second bound in Theorem 4 for n=m implies Clarkson and Shor’s original bound, O((k+1)cT0(n/(k+1))), for uncolored cases [20, Theorem 3.1]. Note that Theorem 4 still implies the same Clarkson–Shor bound if T0 is linear. Moreover, if the colors are assigned in a favorably uniform way, we can derive similar bounds as well without assuming the convexity of T0.

Theorem 5.

With the above notations, suppose |Si|ρnm for every iK for some ρ1. Then, for 1cd,

b=1cj=0mb(mbjcb)|b,j(S,κ)|(mc)T0(cρnm),

and for 2cd and 0kmc1,

j=0k|c,j(S,κ)|=O((k+1)cT0(ρnk+1)).

3.1 Colored 𝒋-facets

Let S be a set of n points in d. A j-facet in S is an oriented (d1)-simplex σ with its vertices chosen from S such that the open half-space on its positive side σ+ contains exactly j points of S. Among the first applications of the original Clarkon–Shor framework was the tight upper bound O(kd/2nd/2) on the number of (k)-facets [20], implying the same asymptotic bound on the (k)-level in the arrangement of hyperplanes via the point-to-hyperplane duality [21]. Many variants of j-facets have been discussed in the literature; we refer to a survey article by Wagner [44].

Figure 4: Colored j-facets in colored points in 2: a 1-chromatic 2-facet σ1 and a 2-chromatic 2-facet σ2 are shown, which choose half-planes σ1+ and σ2+ on their right side.

Now, we assume that the points in S are colored by a color assignment κ with m colors K. For any subset Ad, we shall say that A intersects a color iK, if A contains some sSi. According to our notion of colored configurations, a colored j-facet σ in S with respect to κ is an oriented simplex defined by d points DσS such that exactly j colors, but none of the defining colors in κ(Dσ), are intersected by σ+. See Figure 4. Notice that colored j-facets correspond to vertices of the arrangement of m lower/upper envelopes of hyperplanes dual to Si for iK; see the full version [10] for a more detailed discussion.

For 1cd and j0, let ec,j(S) be the number of c-chromatic j-facets in S and ej(S):=cec,j(S) be the number of all j-facets. Katoh and Tokuyama [27, Proposition 15] proved that ek(S)=O(k1/3n) in 2 and ek(S)=O(k2/3n2) in 3 based on a generalized Lovász’s Lemma. Theorems 4 and 5 directly imply the following bounds.

Corollary 6.

For a set S of n colored points in d with m colors and any 0kmd1, the number of (k)-facets in S is j=0kej(S)=O(md/21kd/2nd/2). If there is a constant ρ1 such that |Si|ρnm for every iK, then the bound is improved to j=0kej(S)=O(kd/2nd/2).

Note that for large k with kmd, the above bound on the number of (k)-facets is asymptotically the same as the total number of configurations (see Theorems 4 and 5).

We continue our discussion for the case of d=3. Lemma 3 implies that e0(S) counts the number of facets of the convex hull of S in 3. Using this fact, we observe the following.

Lemma 7.

Let 2rm and RK be a random set of r colors. It holds that

(mr)𝐄[e0(SR)]2(m1r1)n4(mr)

with equality when the points in S are in convex and general position.

Now, suppose that S is in convex and general position. Then, Lemmas 2 and 7 provide two different ways of exactly counting (mr)𝐄[e0(SR)] for 2rm, resulting in:

Theorem 8.

Let S3 be a set of n points in convex and general position, each of which is colored by one of m colors. Then, it holds that for each 0jm2

e3,j(S)+i=0je2,i(S)+i=0j(ji+1)e1,i(S)=2(j+1)(nj2).

Note that Theorem 8 reveals an exact equation on ej(S)=e3,j(S)+e2,j(S)+e1,j(S) for each 0jm2. If m=n, that is, |Si|=1 for all iK, we have e1,j(S)=e2,j(S)=0 for every j, so the equality ej(S)=2(j+1)(nj2) holds. This exact number for the case of m=n was proved earlier by Clarkson and Shor [20, Theorem 3.5].

3.2 Euclidean color Voronoi diagrams

Suppose that S consists of n points in general position in 2, with a given color assignment κ:SK={1,,m}, and δs(x)=xs2 is the Euclidean distance for each sS and any x2. Consider 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) for 1km1 in this setting.

Figure 5: Selected new vertices (small squares) in 𝖢𝖵𝖣2(S) and 𝖢𝖵𝖣¯2(S) and their corresponding circles. Green vertices are 2-chromatic, while orange ones are 3-chromatic. Observe that exactly one color is intersected by the interior C^ of each circle C in (a) and the exterior C¯ of each circle C in (b).

We consider all circles through any three points in S with no regards of colors and let (S) and ¯(S) be the sets of the interiors and exteriors, respectively, of these circles. Also, define two conflict relations χS×(S) and χ¯Sׯ(S) to be the inclusion relation. We then consider colored configurations (S,κ) and ¯(S,κ) with respect to the given color assignment κ. Observe that each colored configuration of weight j in (S,κ) or in ¯(S,κ) corresponds to a new vertex of 𝖢𝖵𝖣j+1(S) or of 𝖢𝖵𝖣¯j+1(S), respectively, by Lemma 1 and the discussions in Section 2, see Figure 5 illustrating the case of j=1. Therefore, for each 1c3 and 0jmc, we have vc,j(S)=|c,j(S,κ)| and v¯c,j(S)=|¯c,j(S,κ)|.

Now, consider the well-known lifting that maps points p=(x,y) in 2 onto the unit paraboloid U={z=x2+y2} in 3: p=(x,y)p=(x,y,x2+y2). Let S={ssS} be the set of lifted colored points in 3. (The horizontal plane {z=0} is identified as the original plane 2.) Consider colored j-facets in S as in the previous section, and recall that ec,j(S) denotes the number of c-chromatic j-facets in S. We then observe the following.

Lemma 9.

For 1c3 and 0jmc, we have vc,j(S)+v¯c,j(S)=ec,j(S).

Since S is in convex and general position, Theorem 8, together with Lemmas 1 and 9, implies an exact equation on the number of vertices in 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S).

Theorem 10.

Let S be a set of n points with m colors in general position in the Euclidean plane, and 1km1. The total number of vertices in 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) is exactly

4k(nk)2n2i=0k2e2,i(S)i=0k1(2k2i1)e1,i(S).

Theorem 10 implies the O(k(nk)) bound on the complexity of 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) by Lemma 1. An interesting special case of the above result is the following.

Corollary 11.

Given a set S of n colored points with m colors in the Euclidean plane, the complexity of both 𝖥𝖢𝖵𝖣(S) and 𝖧𝖵𝖣(S) is bounded by O(m(nm+1)).

4 Color Voronoi diagrams under general distance functions

We extend our results for the Euclidean case to general distance functions. We continue the discussions from Section 2, so S is a set of n sites, colored with m colors from K by a color assignment κ, and the functions δs:2 for sS are in general position.

Notice that the vertices of the arrangements 𝒜Γ and 𝒜¯Γ form two set systems that fit in our colorful framework. More precisely, let (S) be the set of vertices of the arrangement of n surfaces in Γ, and consider two conflict relations χ,χ¯S×(S) such that (s,v)χ if v(S) lies above surface γsΓ and (s,v)χ¯ if v lies below γs. Based on two CS-structures (S,(S),χ) and (S,(S),χ¯), we consider their colored configurations with respect to κ, denoted by (S,κ) and ¯(S,κ), respectively. By this construction, we have vc,j(S)=|c,j(S,κ)| and v¯c,j(S)=|¯c,j(S,κ)|, counting c-chromatic weight-j colored configurations in (S,κ) and in ¯(S,κ), respectively, and, simultaneously, counting new c-chromatic vertices in 𝖢𝖵𝖣j+1(S) and in 𝖢𝖵𝖣¯j+1(S) by Lemma 1.

For each SS, let v0(S) and u0(S) denote the numbers of vertices and unbounded edges111Hereafter, by counting unbounded edges, we mean counting vertices at infinity. So, if an unbounded edge separates the plane 2, then it is counted twice. , respectively, in 𝖵𝖣(S); let v¯0(S) and u¯0(S) denote the numbers of vertices and unbounded edges, respectively, in 𝖥𝖵𝖣(S). We consider the following conditions.

V1

v0(S)=2|S|2u0(S) for any SS.

V2

v¯0(S)=u¯0(S)2 for any SS.

Note that if every region in 𝖵𝖣(S) is nonempty and simply connected, then Euler’s formula and the general position assumption imply condition V1. If 𝖥𝖵𝖣(S) forms a tree, then every face of 𝖥𝖵𝖣(S) is unbounded and thus condition V2 holds by Euler’s formula.

In addition, for c{1,2} and j0, let uc,j=uc,j(S) be the number of c-chromatic unbounded edges in 𝒜Γ that lie above exactly j surfaces in {Ei}iK, and u¯c,j=u¯c,j(S) be the number of c-chromatic unbounded edges in 𝒜¯Γ that lie below exactly j surfaces in {E¯i}iK. From the discussion in Section 2, observe that uc,j and u¯c,j are equal to the numbers of new c-chromatic unbounded edges in 𝖢𝖵𝖣j+1(S) and in 𝖢𝖵𝖣¯j+1(S), respectively. Further, as for vc,j and v¯c,j, note that uc,j and u¯c,j indeed count c-chromatic weight-j colored configurations based on two CS-structures for unbounded edges in the arrangement of n surfaces in Γ. Hence, assuming V1 and V2, Lemma 3 implies: for any subset RK,

c=13vc,0(SR)=2|SR|2c=12uc,0(SR) and c=13v¯c,0(SR)=c=12u¯c,0(SR)2,

since 𝖢𝖵𝖣1(SR)=𝖵𝖣(SR) and 𝖢𝖵𝖣¯1(SR)=𝖥𝖵𝖣(SR). Combining these equations and the others obtained by Lemma 2, we have two systems of linear equations that can be solved in a similar way as done in Theorem 8. For 0jm1, define

Vj:=v3,j+i=0j(v2,i+(ji+1)v1,i), and  Uj:=i=0j(u2,i+(ji+1)u1,i),
V¯j:=v¯3,j+i=0j(v¯2,i+(ji+1)v¯1,i), and  U¯j:=i=0j(u¯2,i+(ji+1)u¯1,i).
Lemma 12.

With the above notations, let 0jm2. Condition V1 implies Vj+Uj=(j+1)(2nj2); condition V2 implies V¯jU¯j=(j+1)(j+2).

Theorem 13.

Given S and {δs}sS in general position as above, let 1km1 be an integer. If condition V1 is true, then the number of vertices in 𝖢𝖵𝖣k(S) is exactly

2k(2nk)2n2i=0k2v2,i(S)i=0k1(2k2i1)v1,i(S)Uk1Uk2;

if condition V2 is true, the number of vertices in 𝖢𝖵𝖣¯k(S) is exactly

U¯k1+U¯k22k22i=0k2v¯2,i(S)i=0k1(2k2i1)v¯1,i(S).

Remark that condition V1 already implies the asymptotic complexity O(k(nk)) of 𝖢𝖵𝖣k(S) for kn2 as 2nk3(nk), while we have O(kn) for k>n2. Further, to show the O(k(nk)) bound for any value of k for both 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S), it suffices to show that Uj(j+1)(j+2)o(j2) and U¯j(j+1)(2nj2)+o(j2). In this way, Theorem 13 reduces the problem of bounding the complexity of higher-order color Voronoi diagrams to that of bounding the number of their unbounded edges. Also, note that Lemma 12 implies Uj(j+1)(2nj2)andU¯j(j+1)(j+2), if conditions V1 and V2 hold.

Remark also that if 𝖵𝖣(S) and 𝖥𝖵𝖣(S) fall under the umbrella of abstract Voronoi diagrams, then conditions V1 and V2 hold [28, 35], so Theorem 13 implies:

Corollary 14.

Suppose that S and {δs}sS imply a bisector system that satisfies the conditions of abstract Voronoi diagrams [28]. Then, the complexity of 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣k(S) is O(k(nk)) for 1kn2 and O(kn) for n2+1km.

The quantities Uj and U¯j, related to the number of unbounded edges, often turn out to be equal; the very typical example is the Euclidean case for point sites S, where the equality uc,j(S)=u¯c,j(S)=ec,j(S) holds. This inspires us to consider the following third condition:

V3

u0(S)=u¯0(S) for any SS.

Lemma 15.

Condition V3 implies Uj=U¯j for any 0jm1.

Assuming conditions V1V3, we obtain the same exact number as in Theorem 10.

Theorem 16.

Given S and {δs}sS in general position as above, if conditions V1V3 hold, then the total number of vertices in 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) for 1km1 is exactly

4k(nk)2n2i=0k2(v2,i(S)+v¯2,i(S))i=0k1(2k2i1)(v1,i(S)+v¯1,i(S)).

Below, we discuss some specific cases of functions δs for a set S of points in the plane 2, in which new results are derived by applying Theorems 13 and 16.

Convex distance functions.

From now on, suppose S consists of n colored points in 2. Let B2 be any convex and compact body whose interior contains the origin. Define δs(x)=xsB for point sS to be the convex distance from x2 to s based on B [19, 31]. Since Voronoi diagrams of point sites under a convex distance function fall under the model of abstract Voronoi diagrams [28, 35], conditions V1 and V2 hold.

Condition V3, however, is not guaranteed in general; a popular example is the L1 or L metric, under which 𝖵𝖣(S) may have Θ(|S|) parallel unbounded edges while 𝖥𝖵𝖣(S) has at most four. In the following, we first assume that B is smooth that is, there is a unique line tangent to B at each point on its boundary [31]. We then make the following observation for a smooth B, stronger than condition V3, which has been known for the Euclidean metric even in higher dimensions [11].

Lemma 17.

Given S and δs for sS as above, suppose B is smooth. For c{1,2} and 0jm1, we have uc,j(S)=u¯c,j(S)=ec,j(S), the number of c-chromatic j-facets in S.

By Lemma 17, Theorem 16 implies the same upper bound 4k(nk)2n on the total number of vertices in both 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) under any smooth convex distance function.

We then relax the smoothness of B by a limit argument, so let B be any convex and compact body. Consider a sequence of smooth and convex bodies B0,B1, that converges to B. Obviously, there exists B^=Bi sufficiently close to B so that: (a) the functions δ^s(x)=xsB^ for sS are in general position (see Section 2), and (b) for any scaled and translated copy Bpqr of B having three points p,q,rS on its boundary, there is also a scaled and translated copy B^pqr of B^, whose boundary goes through p,q,r and the separation of S by B^pqr is the same as that by Bpqr.

This implies that the vertices in 𝖢𝖵𝖣k(S) and in 𝖢𝖵𝖣¯k(S) under the convex distance function based on B are preserved in their counterpart diagrams under the convex distance function based on B^. Furthermore, the general position assumption can be relaxed, as it does not decrease the number of vertices in the diagrams. Hence, we conclude:

Corollary 18.

Let S be a set of n colored points in 2 with m colors. For 1km1, the total number of vertices in 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) under any Lp metric for 1p or any convex distance function is at most 4k(nk)2n.

It is worth noting that Lemmas 12 and 17 imply bounds for colored j-facets in 2.

Corollary 19.

Let S be a set of n colored points in 2 with m colors. For 0km2,

(k+1)(k+2)j=0ke2,j(S)+j=0k(kj+1)e1,j(S)(k+1)(2nk2).

More on polygonal convex distance functions.

First, we consider the L metric, so B is the unit square centered at the origin. Liu, Papadopoulou, and Lee [30] proved an upper bound of O((nk)2) for ordinary order-k Voronoi diagrams of n points under the L metric using the Hanan grid. An analogous argument can also be applied to our color diagrams.

Lemma 20.

Under the L metric, the number of vertices in 𝖢𝖵𝖣k(S) is at most 4(nk)2.

Hence, the complexity of 𝖢𝖵𝖣k(S) under the L metric is O(min{k(nk),(nk)2}). For the maximal counterpart 𝖢𝖵𝖣¯k(S), we prove the following.

Lemma 21.

Under the L metric, for any 0jm2, we have U¯j2(j+1)(j+2). Therefore, the number of vertices of 𝖢𝖵𝖣¯k(S) for 1km1 is at most 2k2.

Summarizing, we obtain:

Theorem 22.

Let S be a set of n colored points with m colors in the L or L1 plane. For 1km1, the number of vertices in 𝖢𝖵𝖣k(S) is at most min{4k(nk)2n,4(nk)2} and the number of vertices in 𝖢𝖵𝖣¯k(S) is at most min{4k(nk)2n,2k2}.

The above approach also works for polygonal convex distances, concluding the following.

Corollary 23.

Let B be a convex 2b-gon with 2b4, centrally symmetric around the origin, and S a set of n colored points with m colors. For 1km1, 𝖢𝖵𝖣k(S) and 𝖢𝖵𝖣¯k(S) under the convex distance function based on B consist of at most min{4k(nk)2n,2(b2b)(nk)2} and min{4k(nk)2n,2(b2b1)k2} vertices, respectively.

Remark that a more careful analysis could reduce the factor depending on b, and relax the central symmetry of B.

5 Iterative algorithms for color Voronoi diagrams

In this section, we present an iterative approach to compute the order-k color Voronoi diagrams and their refined counterparts for an increasing order of k. Recall that S is a set of n sites associated with distance functions δs for sS. We assume the general position assumption on {δs}sS given in Section 2. We first establish some key structural properties, which add the concept of color to well-known properties of order-k Voronoi diagrams.

Consider a face f of an order-k Voronoi region Rk(H;S) of 𝖢𝖵𝖣k(S), or a region R¯k(H;S) of 𝖢𝖵𝖣¯k(S), where HK with |H|=k. Recall that SHS is the set of sites whose colors are included in H. Let SfS be the set of sites that, together with sites in SH, define the edges along the boundary of f. The following property is directly derived from the definitions.

Lemma 24.

Let fRk(H;S) be a face of 𝖢𝖵𝖣k(S) for 1km1. It holds that:

  1. (i)

    𝖢𝖵𝖣1(Sf)f=𝖢𝖵𝖣k+1(S)f and 𝖵𝖣(Sf)f=𝖢𝖵𝖣k+1(S)f.

  2. (ii)

    𝖥𝖢𝖵𝖣(SH)f=𝖢𝖵𝖣k1(S)f and 𝖥𝖢𝖵𝖣(SH)f=𝖢𝖵𝖣k(S)f.

Analogous claims hold for the maximal diagrams, however, for an unbounded face f of 𝖢𝖵𝖣¯k(S), the set Sf is no longer adequate to derive the portion of 𝖢𝖵𝖣¯k+1(S) that lies within f. We need the set Sf+SSH that defines the unbounded faces of 𝖢𝖵𝖣¯k+1(S)f.

Lemma 25.

Let fR¯k(H;S) be a face of 𝖢𝖵𝖣¯k(S) for 1km1. Let Sf+SSH be the set of sites that define unbounded faces in 𝖢𝖵𝖣¯k+1(S)f; Sf+= if f is bounded. The following hold:

  1. (i)

    𝖢𝖵𝖣¯1(SfSf+)f=𝖢𝖵𝖣¯k+1(S)f and 𝖥𝖵𝖣(SfSf+)f=𝖢𝖵𝖣¯k+1(S)f.

  2. (ii)

    𝖧𝖵𝖣(SH)f=𝖢𝖵𝖣¯k1(S)f and 𝖧𝖵𝖣(SH)f=𝖢𝖵𝖣¯k(S)f.

Using Lemma 24(i), we can iteratively compute 𝖢𝖵𝖣k+1(S), given 𝖢𝖵𝖣k(S). However, to iteratively compute 𝖢𝖵𝖣¯k+1(S), we first need to identify the sites that define the new unbounded edges of 𝖢𝖵𝖣¯k+1(S). This information, however, is not encoded in 𝖢𝖵𝖣¯k(S). We give a condition, related to condition V3, under which we can use 𝖢𝖵𝖣k+1(S) to derive the information missing from 𝖢𝖵𝖣¯k+1(S). This condition is satisfied if, for example, S is a set of points and δs for sS is a convex distance based on a smooth body.

V3

The unbounded faces of 𝖵𝖣(S) and 𝖥𝖵𝖣(S) are defined by the same sequence of sites, for any SS.

Figure 6: Illustration of condition V3 and Lemma 26 with the same set S={s1,,s9} as in Figures 23 under the Euclidean metric. In (a)(b), 𝖵𝖣(S) and 𝖥𝖵𝖣(S) have the same sequence of sites that define unbounded faces for any SS, so condition V3 holds. The shaded region in (a) is a face f of 𝖢𝖵𝖣1(S) corresponding to face f of 𝖢𝖵𝖣¯1(S) shaded in (b). In (c)(d), shaded regions show how the portions of f and f at infinity are subdivided in 𝖢𝖵𝖣2(S) and 𝖢𝖵𝖣¯2(S).
Lemma 26.

Condition V3 implies that the unbounded faces of 𝖢𝖵𝖣k(S) and of 𝖢𝖵𝖣¯k(S) are defined by the same sequence of sites for any 1km. (See Figure 6.)

Now, we assume that the sites S and their distance functions δs fall under the model of abstract Voronoi diagrams [28]. Furthermore, we assume that any distance or bisector can be computed in O(1) time. We then conclude the following using procedures in [6, 39, 3, 25].

Theorem 27.

Let S and {δs}sS be a set of n colored sites with m colors and distance functions that satisfy the conditions of abstract Voronoi diagrams. Then, for 1km, in O(k2n+nlogn) expected time or in O(k2nlogn) worst-case time, we can compute 𝖢𝖵𝖣1(S),,𝖢𝖵𝖣k(S). If in addition condition V3 holds, then we can also compute 𝖢𝖵𝖣¯1(S),,𝖢𝖵𝖣¯k(S) in the same time bound. If S consists of points and δs(x) is the Euclidean distance to sS, the time bound is reduced to O(k2n+nlogn) in the worst case.

The convex distance functions satisfy the conditions of Theorem 27, hence we have:

Corollary 28.

Let B be a convex and compact body in 2 of a constant complexity that contains the origin in its interior. Given a set S of n colored points in 2 with m colors and an integer 1km, we can compute 𝖢𝖵𝖣1(S),,𝖢𝖵𝖣k(S) in O(k2n+nlogn) expected time or in O(k2nlogn) worst-case time. If B is smooth, then we can also compute 𝖢𝖵𝖣¯1(S),,𝖢𝖵𝖣¯k(S) in the same time bound.

References

  • [1] Manuel Abellanas, Ferran Hurtado, Christian Icking, Rolf Klein, Elmar Langetepe, Lihong Ma, Belén Palop, and Vera Sacristán. The farthest color Voronoi diagram and related problems. Technical Report 002, Rheinische Friedrich–Wilhelms–Universität Bonn, 2006.
  • [2] Pankaj K. Agarwal, Mark de Berg, Jiří Matoušek, and Otfried Schwarzkopf. Constructing levels in arrangements and higher order Voronoi diagrams. SIAM J. Comput., 27(3):654–667, 1998. doi:10.1137/S0097539795281840.
  • [3] Alok Aggarwal, Leonidas J. Guibas, James B. Saxe, and Peter W. Shor. A linear-time algorithm for computing the Voronoi diagram of a convex polygon. Discrete Comput. Geom., 4(1):591–604, 1989. doi:10.1007/BF02187749.
  • [4] Elena Arseneva and Evanthia Papadopoulou. Randomized incremental construction for the Hausdorff Voronoi diagram revisited and extended. Journal of Combinatorial Optimization, 37(2):579–600, February 2019. doi:10.1007/s10878-018-0347-x.
  • [5] Franz Aurenhammer. Power diagrams: properties, algorithms and applications. SIAM J. Comput., 16(1):78–96, 1987. doi:10.1137/0216006.
  • [6] Franz Aurenhammer, Rolf Klein, and Der-Tsai Lee. Voronoi Diagrams and Delaunay Triangulations. World Scientific, Singapore, 2013. doi:10.1142/8685.
  • [7] Franz Aurenhammer and Otfried Schwarzkopf. A simple on-line randomized incremental algorithm for computing higher order Voronoi diagram. Int. J. Comput. Geom. Appl., 2(4):363–381, 1992. doi:10.1142/S0218195992000214.
  • [8] Sang Won Bae. On linear-sized farthest-color Voronoi diagrams. IEICE Trans. Inf. Syst., 95-D(3):731–736, 2012. doi:10.1587/transinf.E95.D.731.
  • [9] Sang Won Bae. Tight bound and improved algorithm for farthest-color Voronoi diagrams of line segments. Comput. Geom.: Theory Appl., 47(8):779–788, 2014. doi:10.1016/j.comgeo.2014.04.005.
  • [10] Sang Won Bae, Nicolau Oliver, and Evanthia Papadopoulou. Higher-order color Voronoi diagrams and the colorful Clarkson–Shor framework, 2025. arXiv:2504.06960.
  • [11] Gill Barequet, Evanthia Papadopoulou, and Martin Suderland. Unbounded regions of high-order Voronoi diagrams of lines and line segments in higher dimensions. Discrete & Computational Geometry, 72(3):1304–1332, May 2023. doi:10.1007/s00454-023-00492-2.
  • [12] Ranita Biswas, Sebastiano Cultrera di Montesano, Herbert Edelsbrunner, and Morteza Saghafian. Counting cells of order-k Voronoi tessellations in 3 with Morse theory. In Kevin Buchin and Éric Colin de Verdière, editors, Proceedings of the 37th International Symposium on Computational Geometry (SoCG 2021), volume 189 of LIPIcs, pages 16:1–16:15. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2021. doi:10.4230/LIPIcs.SoCG.2021.16.
  • [13] Cecilia Bohler, Panagiotis Cheilaris, Rolf Klein, Chih-Hung Liu, Evanthia Papadopoulou, and Maksym Zavershynskyi. On the complexity of higher order abstract Voronoi diagrams. Comput. Geom.: Theory Appl., 48(8):539–551, 2015. doi:10.1016/j.comgeo.2015.04.008.
  • [14] Cecilia Bohler, Rolf Klein, and Chih-Hung Liu. An efficient randomized algorithm for higher-order abstract Voronoi diagrams. Algorithmica, 81(6):2317–2345, 2019. doi:10.1007/s00453-018-00536-7.
  • [15] Cecilia Bohler, Chih-Hung Liu, Evanthia Papadopoulou, and Maksym Zavershynskyi. A randomized divide and conquer algorithm for higher-order abstract Voronoi diagrams. Comput. Geom.: Theory Appl., 59:26–38, 2016. doi:10.1016/j.comgeo.2016.08.004.
  • [16] Timothy M. Chan. Random sampling, halfspace range reporting, and construction of (k)-levels in three dimensions. SIAM J. Comput., 30(2):561–575, 2000. doi:10.1137/S0097539798349188.
  • [17] Timothy M. Chan, Pingan Cheng, and Da Wei Zheng. An optimal algorithm for higher-order Voronoi diagrams in the plane: The usefulness of nondeterminism. In David P. Woodruff, editor, Proceedings of the 2024 ACM-SIAM Symposium on Discrete Algorithms, SODA 2024, Alexandria, VA, USA, January 7-10, 2024, pages 4451–4463. SIAM, 2024. doi:10.1137/1.9781611977912.156.
  • [18] Timothy M. Chan and Konstantinos Tsakalidis. Optimal deterministic algorithms for 2-d and 3-d shallow cuttings. Discrete Comput. Geom., 56(4):866–881, 2016. doi:10.1007/s00454-016-9784-4.
  • [19] L. Paul Chew and Robert L.S. Drysdale III. Voronoi diagrams based on convex distance functions. In Proceedings of the First Annual Symposium on Computational Geometry, Baltimore, Maryland, USA, June 5-7, 1985, pages 235–244. ACM, 1985. doi:10.1145/323233.323264.
  • [20] Kenneth L. Clarkson and Peter W. Shor. Application of random sampling in computational geometry, II. Discrete Comput. Geom., 4:387–421, 1989. doi:10.1007/BF02187740.
  • [21] Herbert Edelsbrunner. Algorithms in Combinatorial Geometry. Springer-Verlag, 1987.
  • [22] Herbert Edelsbrunner, Leonidas J. Guibas, and Micha Sharir. The upper envelope of piecewise linear functions: Algorithms and applications. Discrete Comput. Geom., 4(4):311–336, 1989. doi:10.1007/BF02187733.
  • [23] Herbert Edelsbrunner and Raimund Seidel. Voronoi diagrams and arrangements. Discrete Comput. Geom., 1:25–44, 1986. doi:10.1007/BF02187681.
  • [24] Daniel P. Huttenlocher, Klara Kedem, and Micha Sharir. The upper envelope of Voronoi surfaces and its applications. Discrete Comput. Geom., 9:267–291, 1993. doi:10.1007/BF02189323.
  • [25] Kolja Junginger and Evanthia Papadopoulou. Deletion in abstract Voronoi diagrams in expected linear time and related problems. Discrete & Computational Geometry, 69(4):1040–1078, March 2023. doi:10.1007/s00454-022-00463-z.
  • [26] Haim Kaplan, Wolfgang Mulzer, Liam Roditty, Paul Seiferth, and Micha Sharir. Dynamic planar Voronoi diagrams for general distance functions and their algorithmic applications. Discrete Comput. Geom., 64(3):838–904, 2020. doi:10.1007/s00454-020-00243-7.
  • [27] Naoki Katoh and Takeshi Tokuyama. k-Levels of concave surfaces. Discrete Comput. Geom., 27:567–584, 2002. doi:10.1007/s00454-001-0086-z.
  • [28] Rolf Klein. Concrete and Abstract Voronoi Diagrams, volume 400 of Lecture Notes in Computer Science. Springer-Verlag, Berlin, Germany, 1989. doi:10.1007/3-540-52055-4.
  • [29] Der-Tsai Lee. On k-nearest neighbor Voronoi diagrams in the plane. IEEE Trans. Computers, 31(6):478–487, 1982. doi:10.1109/TC.1982.1676031.
  • [30] Chih-Hung Liu, Evanthia Papadopoulou, and Der-Tsai Lee. The k-nearest-neighbor Voronoi diagram revisited. Algorithmica, 71(2):429–449, 2015. doi:10.1007/s00453-013-9809-9.
  • [31] Lihong Ma. Bisectors and Voronoi Diagams for Convex Distance Functions. PhD thesis, Fern Unversität Hagen, 2000.
  • [32] Ioannis Mantas, Evanthia Papadopoulou, Vera Sacristán, and Rodrigo I. Silveira. Farthest color Voronoi diagrams: Complexity and algorithms. In Proc. LATIN 2020, volume 12118 of Lecture Notes in Computer Science, pages 283–295, 2020. doi:10.1007/978-3-030-61792-9_23.
  • [33] Ioannis Mantas, Evanthia Papadopoulou, Rodrigo I. Silveira, and Zeyu Wang. The farthest color Voronoi diagram in the plane. Algorithmica to appear, preprint available at Research Square. doi:10.21203/rs.3.rs-4644060/v1.
  • [34] Jiří Matoušek. Lectures on Discrete Geometry. Springer, 2002.
  • [35] Kurt Mehlhorn, Stefan Meiser, and Ronald Rasch. Furthest site abstract Voronoi diagrams. Int. J. Comput. Geom. Appl., 11(6):583–616, 2001. doi:10.1142/S0218195901000663.
  • [36] Ketan Mulmuley. On levels in arrangements and Voronoi diagrams. Discrete Comput. Geom., 6(1):307–338, 1991. doi:10.1007/BF02574692.
  • [37] Evanthia Papadopoulou. Critical area computation for missing material defects in VLSI circuits. IEEE Trans. Comput. Aided Des. Integr. Circuits Syst., 20(5):583–597, 2001. doi:10.1109/43.920683.
  • [38] Evanthia Papadopoulou. The Hausdorff Voronoi diagram of point clusters in the plane. Algorithmica, 40(2):63–82, 2004. doi:10.1007/s00453-004-1095-0.
  • [39] Evanthia Papadopoulou. Abstract Voronoi-like Graphs: Extending Delaunay’s Theorem and Applications. In Proc. 39th International Symposium on Computational Geometry (SoCG), volume 258 of LIPIcs, pages 52:1–52:16, 2023. doi:10.4230/LIPIcs.SoCG.2023.52.
  • [40] Evanthia Papadopoulou and Maksym Zavershynskyi. The higher-order Voronoi diagram of line segments. Algorithmica, 74(1):415–439, 2016. doi:10.1007/s00453-014-9950-0.
  • [41] Edgar A. Ramos. On range reporting, ray shooting, and k-level construction. In Proceedings of the Fifteenth Annual Symposium on Computational Geometry, Miami Beach, Florida, USA, June 13-16, 1999, pages 390–399. ACM, 1999. doi:10.1145/304893.304993.
  • [42] Micha Sharir. The Clarkson–Shor technique revisited and extended. Comb. Probab. Comput., 12(2):191–201, 2003. doi:10.1017/S0963548302005527.
  • [43] Micha Sharir and Pankaj K. Agarwal. Davenport-Schinzel Sequences and Their Geometric Applications. Cambridge University Press, 1995.
  • [44] Uli Wagner. k-sets and k-facets. In Jacob E. Goodman, János Pach, and Richard Pollack, editors, Surveys on Discrete and Computational Geometry, volume 453 of Contemporary Mathematics, pages 443–513. Amer. Math. Soc., 2008.