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 s∈S 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 s∈S 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 4⁢k⁢(n−k)−2⁢n 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 1⩽p⩽∞. 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⁢(n−k),(n−k)2}) complexity, while its maximal counterpart has O⁢(min⁡{k⁢(n−k),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 Si⊆S be the set of sites in color i∈K. We consider two distance-to-color functions from any point x∈ℝ2 to each color i∈K:

di⁢(x):=mins∈Si⁡δs⁢(x) and d¯i⁢(x):=maxs∈Si⁡δs⁢(x),

where δs⁢(x) denotes the prescribed distance-to-site function from x∈ℝ2 to site s∈S.

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

Rk⁢(H;S) :={x∈ℝ2∣di⁢(x)<dj⁢(x)⁢ for all i∈H and j∈K∖H}
R¯k⁢(H;S) :={x∈ℝ2∣d¯i⁢(x)>d¯j⁢(x)⁢ for all i∈H and j∈K∖H},

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 H⊂K 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}i∈K, while 𝖢𝖵𝖣¯k⁢(S) partitions ℝ2 by k farthest colors under the maximal distances {d¯i}i∈K.

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)=𝖵𝖣n−k⁢(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 i∈K that consist of all points p∈ℝ2 whose farthest color is i with respect to the minimal distances {di}i∈K. 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}i∈K. We thus have 𝖥𝖢𝖵𝖣⁢(S)=𝖢𝖵𝖣m−1⁢(S) and 𝖧𝖵𝖣⁢(S)=𝖢𝖵𝖣¯m−1⁢(S). These two diagrams are often refined so that the region of each color i∈K 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 s∈S. Our main results are as follows:

  1. (1)

    We prove an exact upper bound 4⁢k⁢(n−k)−2⁢n on the total number of vertices of both order-k diagrams 𝖢𝖵𝖣k⁢(S) and 𝖢𝖵𝖣¯k⁢(S) for each 1⩽k⩽m−1, when the sites s∈S are points and δs⁢(x) is given as the convex distance from x∈ℝ2 to s based on any convex and compact subset of ℝ2, which includes the Lp distance for any 1⩽p⩽∞. 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 S′⊆S. The bound 4⁢k⁢(n−k)−2⁢n 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⁡{4⁢k⁢(n−k)−2⁢n,4⁢(n−k)2} and min⁡{4⁢k⁢(n−k)−2⁢n,2⁢k2} 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⁢(k2⁢n+n⁢log⁡n) expected or O⁢(k2⁢n⁢log⁡n) worst-case time, provided that S and {δs}s∈S 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⁢(k2⁢n+n⁢log⁡n) 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⁢(n−k)) 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⁢(n−k),(n−k)2}) for a set S of points or line segments [30, 40].

In particular, the same exact number 4⁢k⁢(n−k)−2⁢n can be derived for the ordinary order-k diagram 𝖵𝖣k⁢(S) and order-(n−k) diagram 𝖵𝖣n−k⁢(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 (4⁢k−2)⁢n−2⁢k2−ek−2⁢∑i=1k−1ei, 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 𝖵𝖣n−k⁢(S) is exactly 4⁢k⁢(n−k)−2⁢n by using the identities: ek=en−k and ∑k=1n−1ek=n⁢(n−1). 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=m−1, which yields the O⁢(m⁢(n−m+1)) bound for the farthest color Voronoi diagram 𝖥𝖢𝖵𝖣⁢(S)=𝖢𝖵𝖣m−1⁢(S) and the Hausdorff Voronoi diagram 𝖧𝖵𝖣⁢(S)=𝖢𝖵𝖣¯m−1⁢(S). The worst-case complexity of 𝖥𝖢𝖵𝖣⁢(S) and 𝖧𝖵𝖣⁢(S) is known to be Θ⁢(m⁢n) 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⁢(m⁢n) 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 m⩽n. Prior to our results, it was known that the complexity of 𝖥𝖢𝖵𝖣⁢(S) and 𝖧𝖵𝖣⁢(S) can range from Θ⁢(n) to Θ⁢(m⁢(n−m+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 1⩽i⩽k, and compares to the O⁢(k2⁢n⁢log⁡n)-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⁢(n⁢log⁡n+k⁢n)-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 s∈S 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 x⁢y-monotone surface {(p,δs⁢(p))∣p∈ℝ2} in ℝ3, and Γ:={γs∣s∈S}. 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 κ:S→K denote a color assignment such that the color of s∈S is κ⁢(s)∈K. Let Si:={s∈S∣κ⁢(s)=i} and Γi:={γs∣s∈Si} for i∈K. 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 i∈K, and E¯i is the graph of the maximal distance d¯i. We then consider the levels of the arrangements 𝒜Γ and 𝒜¯Γ. For 1⩽k⩽m, 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 1⩽k⩽m−1 are the projections of Lk∩Lk+1 and of L¯k∩L¯k+1, respectively, onto ℝ2.

  • ■

    For each 1⩽k⩽m, 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 s∈Si such that, for any point p∈f, di⁢(p)=δs⁢(p) and i∈K is the k-th nearest color from p with respect to {di}i∈K (resp. d¯i⁢(p)=δs⁢(p) and i∈K is the k-th farthest color from p with respect to {d¯i}i∈K). 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)=𝖢𝖵𝖣m−1⁢(S) and 𝖧𝖵𝖣⁢(S)=𝖢𝖵𝖣¯m−1⁢(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 s6∈S3 and those with s9∈S2 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}i∈K or from {E¯i}i∈K, 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 𝖢𝖵𝖣k−1∗⁢(S) or in 𝖢𝖵𝖣¯k−1∗⁢(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 1⩽c⩽3 and 0⩽j⩽m−1 to be the number of c-chromatic vertices in 𝒜Γ below which there are exactly j surfaces from {Ei}i∈K; 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}i∈K.

Lemma 1.

For any 1⩽c⩽3 and 1⩽k⩽m, the following hold:

  1. (i)

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

  2. (ii)

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

  3. (iii)

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

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 d⩾1, 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 s∈S 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 κ:S→K to the objects in S, where K={1,2,…,m} denotes the set of m colors with m⩽n. For each color i∈K, let Si:={s∈S∣κ⁢(s)=i}. For f∈ℱ⁢(S) and set Df⊆S of d objects defining f, κ⁢(Df)={κ⁢(s)∣s∈Df} is called a set of colors defining f. We build a color-to-configuration conflict relation χκ⊆K×ℱ⁢(S) such that a color i∈K is in conflict with a configuration f∈ℱ⁢(S) if an object s∈Si is in conflict with f, that is, (i,f)∈χκ if and only if (s,f)∈χ for some s∈Si. 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 1⩽c⩽d, 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,⋃jℱc,j⁢(S,κ),χκ) for each 1⩽c⩽d 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 1⩽c⩽d, r⩾0, and a⩾0 be integers, and R⊆K be a random subset of r colors. Then,

(mr)⁢𝐄[|ℱc,a⁢(SR,κR)|]⩾∑j=0m−c|ℱc,j⁢(S,κ)|⁢(ja)⁢(m−c−jr−c−a),

where SR=⋃i∈RSi and κR:SR→R 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 S′⊆S, let ℱ0⁢(S′)⊆ℱ⁢(S′) be the set of (uncolored) configurations f of weight 0, that is, (s,f)∉χ for all s∈S′. 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 S′⊆S 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 1⩽c⩽d, the total number of (⩽c)-chromatic colored configurations is bounded by

∑b=1c∑j=0m−b(m−b−jc−b)⁢|ℱb,j⁢(S,κ)|⩽(mc)⋅1m⁢∑i∈KT0⁢(c⁢|Si|)=O⁢(mc−1⋅T0⁢(c⁢n)).

Also, for each 2⩽c⩽d and 0⩽k⩽⌊mc⌋−1, it holds that

∑j=0k|ℱc,j⁢(S,κ)|=O⁢((k+1)cm⋅∑i∈KT0⁢(m⁢|Si|k+1))=O⁢((k+1)cm⋅T0⁢(m⁢nk+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)c⋅T0⁢(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 i∈K for some ρ⩾1. Then, for 1⩽c⩽d,

∑b=1c∑j=0m−b(m−b−jc−b)⁢|ℱb,j⁢(S,κ)|⩽(mc)⁢T0⁢(c⁢ρ⋅nm),

and for 2⩽c⩽d and 0⩽k⩽⌊mc⌋−1,

∑j=0k|ℱc,j⁢(S,κ)|=O⁢((k+1)c⋅T0⁢(ρ⋅nk+1)).

3.1 Colored 𝒋-facets

Let S be a set of n points in ℝd. A j-facet in S is an oriented (d−1)-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⁢(k⌈d/2⌉⁢n⌊d/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 A⊂ℝd, we shall say that A intersects a color i∈K, if A contains some s∈Si. 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 i∈K; see the full version [10] for a more detailed discussion.

For 1⩽c⩽d and j⩾0, 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/3⁢n) in ℝ2 and ek⁢(S)=O⁢(k2/3⁢n2) 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 0⩽k⩽⌊md⌋−1, the number of (⩽k)-facets in S is ∑j=0kej⁢(S)=O⁢(m⌊d/2⌋−1⁢k⌈d/2⌉⁢n⌊d/2⌋). If there is a constant ρ⩾1 such that |Si|⩽ρ⋅nm for every i∈K, then the bound is improved to ∑j=0kej⁢(S)=O⁢(k⌈d/2⌉⁢n⌊d/2⌋).

Note that for large k with k⩾⌊md⌋, 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 2⩽r⩽m and R⊆K be a random set of r colors. It holds that

(mr)⁢𝐄[e0⁢(SR)]⩽2⁢(m−1r−1)⁢n−4⁢(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 2⩽r⩽m, resulting in:

Theorem 8.

Let S⊂ℝ3 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 0⩽j⩽m−2

e3,j⁢(S)+∑i=0je2,i⁢(S)+∑i=0j(j−i+1)⁢e1,i⁢(S)=2⁢(j+1)⁢(n−j−2).

Note that Theorem 8 reveals an exact equation on ej⁢(S)=e3,j⁢(S)+e2,j⁢(S)+e1,j⁢(S) for each 0⩽j⩽m−2. If m=n, that is, |Si|=1 for all i∈K, we have e1,j⁢(S)=e2,j⁢(S)=0 for every j, so the equality ej⁢(S)=2⁢(j+1)⁢(n−j−2) 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 κ:S→K={1,…,m}, and δs⁢(x)=‖x−s‖2 is the Euclidean distance for each s∈S and any x∈ℝ2. Consider 𝖢𝖵𝖣k⁢(S) and 𝖢𝖵𝖣¯k⁢(S) for 1⩽k⩽m−1 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 1⩽c⩽3 and 0⩽j⩽m−c, 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∪={s∪∣s∈S} 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 1⩽c⩽3 and 0⩽j⩽m−c, 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 1⩽k⩽m−1. The total number of vertices in 𝖢𝖵𝖣k⁢(S) and 𝖢𝖵𝖣¯k⁢(S) is exactly

4⁢k⁢(n−k)−2⁢n−2⁢∑i=0k−2e2,i⁢(S∪)−∑i=0k−1(2⁢k−2⁢i−1)⁢e1,i⁢(S∪).

Theorem 10 implies the O⁢(k⁢(n−k)) 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⁢(n−m+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 s∈S 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 S′⊆S, 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′|−2−u0⁢(S′) for any S′⊆S.

V2

v¯0⁢(S′)=u¯0⁢(S′)−2 for any S′⊆S.

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 j⩾0, let uc,j=uc,j⁢(S) be the number of c-chromatic unbounded edges in 𝒜Γ that lie above exactly j surfaces in {Ei}i∈K, 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}i∈K. 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 R⊆K,

∑c=13vc,0⁢(SR)=2⁢|SR|−2−∑c=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 0⩽j⩽m−1, define

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

With the above notations, let 0⩽j⩽m−2. Condition V1 implies Vj+Uj=(j+1)⁢(2⁢n−j−2); condition V2 implies V¯j−U¯j=−(j+1)⁢(j+2).

Theorem 13.

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

2⁢k⁢(2⁢n−k)−2⁢n−2⁢∑i=0k−2v2,i⁢(S)−∑i=0k−1(2⁢k−2⁢i−1)⁢v1,i⁢(S)−Uk−1−Uk−2;

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

U¯k−1+U¯k−2−2⁢k2−2⁢∑i=0k−2v¯2,i⁢(S)−∑i=0k−1(2⁢k−2⁢i−1)⁢v¯1,i⁢(S).

Remark that condition V1 already implies the asymptotic complexity O⁢(k⁢(n−k)) of 𝖢𝖵𝖣k⁢(S) for k⩽n2 as 2⁢n−k⩽3⁢(n−k), while we have O⁢(k⁢n) for k>n2. Further, to show the O⁢(k⁢(n−k)) 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)⁢(2⁢n−j−2)+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)⁢(2⁢n−j−2)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}s∈S 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⁢(n−k)) for 1⩽k⩽⌊n2⌋ and O⁢(k⁢n) for ⌊n2⌋+1⩽k⩽m.

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 S′⊆S.

Lemma 15.

Condition V3 implies Uj=U¯j for any 0⩽j⩽m−1.

Assuming conditions V1–V3, we obtain the same exact number as in Theorem 10.

Theorem 16.

Given S and {δs}s∈S in general position as above, if conditions V1–V3 hold, then the total number of vertices in 𝖢𝖵𝖣k⁢(S) and 𝖢𝖵𝖣¯k⁢(S) for 1⩽k⩽m−1 is exactly

4⁢k⁢(n−k)−2⁢n−2⁢∑i=0k−2(v2,i⁢(S)+v¯2,i⁢(S))−∑i=0k−1(2⁢k−2⁢i−1)⁢(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 B⊂ℝ2 be any convex and compact body whose interior contains the origin. Define δs⁢(x)=‖x−s‖B for point s∈S to be the convex distance from x∈ℝ2 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 s∈S as above, suppose B is smooth. For c∈{1,2} and 0⩽j⩽m−1, 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 4⁢k⁢(n−k)−2⁢n 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)=‖x−s‖B^ for s∈S are in general position (see Section 2), and (b) for any scaled and translated copy Bp⁢q⁢r of B having three points p,q,r∈S on its boundary, there is also a scaled and translated copy B^p⁢q⁢r of B^, whose boundary goes through p,q,r and the separation of S by B^p⁢q⁢r is the same as that by Bp⁢q⁢r.

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 1⩽k⩽m−1, the total number of vertices in 𝖢𝖵𝖣k⁢(S) and 𝖢𝖵𝖣¯k⁢(S) under any Lp metric for 1⩽p⩽∞ or any convex distance function is at most 4⁢k⁢(n−k)−2⁢n.

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 0⩽k⩽m−2,

(k+1)⁢(k+2)⩽∑j=0ke2,j⁢(S)+∑j=0k(k−j+1)⁢e1,j⁢(S)⩽(k+1)⁢(2⁢n−k−2).

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⁢((n−k)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⁢(n−k)2.

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

Lemma 21.

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

Summarizing, we obtain:

Theorem 22.

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

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

Corollary 23.

Let B be a convex 2⁢b-gon with 2⁢b⩾4, centrally symmetric around the origin, and S a set of n colored points with m colors. For 1⩽k⩽m−1, 𝖢𝖵𝖣k⁢(S) and 𝖢𝖵𝖣¯k⁢(S) under the convex distance function based on B consist of at most min⁡{4⁢k⁢(n−k)−2⁢n,2⁢(b2−b)⁢(n−k)2} and min⁡{4⁢k⁢(n−k)−2⁢n,2⁢(b2−b−1)⁢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 s∈S. We assume the general position assumption on {δs}s∈S 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 H⊆K with |H|=k. Recall that SH⊆S is the set of sites whose colors are included in H. Let Sf⊆S 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 f⊆Rk⁢(H;S) be a face of 𝖢𝖵𝖣k⁢(S) for 1⩽k⩽m−1. It holds that:

  1. (i)

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

  2. (ii)

    𝖥𝖢𝖵𝖣⁢(SH)∩f=𝖢𝖵𝖣k−1⁢(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+⊆S∖SH that defines the unbounded faces of 𝖢𝖵𝖣¯k+1⁢(S)∩f.

Lemma 25.

Let f⊆R¯k⁢(H;S) be a face of 𝖢𝖵𝖣¯k⁢(S) for 1⩽k⩽m−1. Let Sf+⊆S∖SH 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⁢(Sf∪Sf+)∩f=𝖢𝖵𝖣¯k+1⁢(S)∩f and 𝖥𝖵𝖣⁢(Sf∪Sf+)∩f=𝖢𝖵𝖣¯k+1∗⁢(S)∩f.

  2. (ii)

    𝖧𝖵𝖣⁢(SH)∩f=𝖢𝖵𝖣¯k−1⁢(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 s∈S 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 S′⊆S.

Figure 6: Illustration of condition V3′ and Lemma 26 with the same set S={s1,…,s9} as in Figures 2–3 under the Euclidean metric. In (a)(b), 𝖵𝖣⁢(S′) and 𝖥𝖵𝖣⁢(S′) have the same sequence of sites that define unbounded faces for any S′⊆S, 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 1⩽k⩽m. (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}s∈S be a set of n colored sites with m colors and distance functions that satisfy the conditions of abstract Voronoi diagrams. Then, for 1⩽k⩽m, in O⁢(k2⁢n+n⁢log⁡n) expected time or in O⁢(k2⁢n⁢log⁡n) 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 s∈S, the time bound is reduced to O⁢(k2⁢n+n⁢log⁡n) 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 1⩽k⩽m, we can compute 𝖢𝖵𝖣1∗⁢(S),…,𝖢𝖵𝖣k∗⁢(S) in O⁢(k2⁢n+n⁢log⁡n) expected time or in O⁢(k2⁢n⁢log⁡n) 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.