Thomas Jefferson High School for Science and Technology, Virginia, USA and ~[email protected] California State University Channel Islands, California, USA and ~[email protected] Montgomery Blair High School, Silver Spring, Maryland, USA and ~[email protected] Department of Computer Science, University of Maryland, College Park, USA and ~[email protected] Department of Computer Science, University of Maryland, College Park, USA and ~[email protected]://orcid.org/0000-0002-5704-312X Mount Holyoke College, Massachusetts, USA and ~[email protected] Department of Computer Science, University of Maryland, College Park, USA and ~[email protected] Princeton University, New Jersey, USA and ~[email protected] Department of Computer Science, University of Maryland, College Park, Maryland, USA and https://www.cs.umd.edu/~mount/[email protected]://orcid.org/0000-0002-3290-8932 \CopyrightH. Banerjee, C. I. Day, M. Hunleth, S. Hwang, A. H. Gezalyan, M. Golovatskaia, N. Parepally, L. Wang, and D. Mount \ccsdesc[500]Theory of computation Computational geometry

Acknowledgements.
The REU grant CNS-2150382 for REU-CAAR funded part of this work. \hideLIPIcs

On The Heine-Borel Property and Minimum Enclosing Balls

Hridhaan Banerjee    Carmen Isabel Day    Megan Hunleth    Sarah Hwang    Auguste H. Gezalyan    Olya Golovatskaia    Nithin Parepally    Lucy Wang    David M. Mount
(March 3, 2025)
Abstract

In this paper, we contribute a proof that the problem of determining the minimum radius balls over metric spaces with the Heine-Borel property is always LP-type. Additionally, we prove that weak metric spaces, those without symmetry, also have this property if we fix the direction in which we take their distances from the centers of the balls. We use this to prove that the minimum radius ball problem is also LP-type in the Thompson metrics and Funk weak metric. We finally examine the LP algorithm and explicit primitives for computing the minimum radius ball in the Hilbert metric.

keywords:
Hilbert metric, Thompson metric, Heine-Borel property, convexity, LP-type problems

1 Introduction

The concept of an LP-type problem originates with the work of Micah Sharir and Emo Welzl in 1992 in their paper "A combinatorial bound for linear programming and related problems"[26]. Since then, a variety of problems have been shown to be LP-type, including minimum radius balls in the Euclidean metric [17] and spaces with Bregman divergences [18], finding the closest distance between two convex polygons [17], various game-theoretic games [11] such as simple stochastic games and parity games. The analysis of LP-type problems has been a continuous area of study in computational geometry since their inception. We contribute to this by proving that determining the minimum radius ball in any metric space satisfying the Heine-Borel property is LP-type. We provide an example with the Hilbert and Thompson metrics. We also show that this holds for weak metric spaces, that is, where the distance function fails to be symmetric.

The Hilbert metric originates from the work of David Hilbert in 1895 in relation to Hilbert’s fourth problem [12]. It presents a non-Euclidean metric in which the triangle inequality is not strict. The Hilbert metric generalizes the Cayley-Klein model of hyperbolic geometry to arbitrary convex bodies in an n𝑛nitalic_n-dimensional space. Definitions will be given in Section 2. Given a convex body, ΩnΩsuperscript𝑛\Omega\subset\mathbb{R}^{n}roman_Ω ⊂ blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT, the Hilbert metric has many desirable properties, such as the fact that straight lines are geodesics and are preserved under projective transformations. In the probability simplex, it provides a natural distance between discrete probability distributions, as shown by the works of Nielsen and Sun [20, 21]. For an excellent resource on Hilbert geometry, see "From Funk to Hilbert Geometry" or the "Handbook of Hilbert Geometry", both by Papadopoulos and Troyanov[23, 22].

The Hilbert metric has seen recent use in a diverse set of fields, especially in that of convexity approximation. This is in particular due to its relationship with Macbeath regions (which are equivalent to Hilbert balls up to a scaling factor) [1, 2] and the flag approximability of polytopes [29]. Other fields in which the Hilbert metric has been used are quantum information theory on convex cones defined by various operators[24], machine learning in the form of clustering and graph embeddings [20, 21], optimal mass transport [5], and a variety of situations in real analysis [13]. Due to its many uses, various algorithms from classical computational geometry have been modified for use in the Hilbert metric, including Voronoi diagrams [10, 4] and Delaunay triangulations [9]. We expand on these works by contributing an algorithm for minimum enclosing balls in the Hilbert polygonal geometry.

The Thompson metric was defined by A. C. Thompson in 1963 as an alternative to the Hilbert metric for its applications in analysis[27]. Like the Hilbert, it provides a metric space over convex bodies and has a similar geometry. Its primary uses are in analysis, in particular as a metric on cones [6, 14, 15, 16]. We contribute a proof that minimum radius balls in this metric are LP-type. Additionally, we contribute the fact that Thompson balls are, like Macbeath Regions, also equivalent to Hilbert balls up to a scaling factor, and therefore induce the same topology on convex bodies.

The Funk weak metric was defined by Paul Funk in 1929 [8]. The Funk metric is a weak metric space that can be used to define both the Hilbert and the Thompson metrics, and is often studied in their context [22, 28] and in the context of flags of polytopes [7]. Because the Funk metric is non-symmetric, it induces a reverse metric called the reverse Funk metric. We contribute a proof that the minimum radius ball problem in these weak metrics are LP-type.

2 Preliminaries

2.1 Metric Spaces

A metric space is the generalization of a distance d:𝐗×𝐗0:𝑑𝐗𝐗superscriptabsent0d:\mathbf{X}\times\mathbf{X}\rightarrow\mathbb{R}^{\geq 0}italic_d : bold_X × bold_X → blackboard_R start_POSTSUPERSCRIPT ≥ 0 end_POSTSUPERSCRIPT on a set 𝐗𝐗\mathbf{X}bold_X. It is a fundamental concept in geometry, real/complex analysis, and topology. We define it here:

Definition 2.1 (Metric Space).

The pair (𝐗,d)𝐗𝑑(\mathbf{X},d)( bold_X , italic_d ) is a metric space if, for all a,b,c𝐗𝑎𝑏𝑐𝐗a,b,c\in\mathbf{X}italic_a , italic_b , italic_c ∈ bold_X:

  1. 1.

    d(a,b)=0𝑑𝑎𝑏0d(a,b)=0italic_d ( italic_a , italic_b ) = 0 iff a=b𝑎𝑏a=bitalic_a = italic_b

  2. 2.

    d(a,b)=d(b,a)𝑑𝑎𝑏𝑑𝑏𝑎d(a,b)=d(b,a)italic_d ( italic_a , italic_b ) = italic_d ( italic_b , italic_a )

  3. 3.

    d(a,c)d(a,b)+d(b,c)𝑑𝑎𝑐𝑑𝑎𝑏𝑑𝑏𝑐d(a,c)\leq d(a,b)+d(b,c)italic_d ( italic_a , italic_c ) ≤ italic_d ( italic_a , italic_b ) + italic_d ( italic_b , italic_c )

When all of the above properties are satisfied except for symmetry, the space is called a weak metric. Given a metric space, (𝐗,d)𝐗𝑑(\mathbf{X},d)( bold_X , italic_d ), the closed ball of radius r𝑟ritalic_r around a point p𝑝pitalic_p generalizes the Euclidean idea of a circle to an arbitrary metric space. We define it here:

Definition 2.2 (Closed Ball).

A closed ball around a point p𝑝pitalic_p in a metric space (𝐗,d)𝐗𝑑(\mathbf{X},d)( bold_X , italic_d ) of radius r𝑟ritalic_r is defined:

B(p,r)={q𝐗d(p,q)r}𝐵𝑝𝑟conditional-set𝑞𝐗𝑑𝑝𝑞𝑟B(p,r)=\{q\in\mathbf{X}\mid d(p,q)\leq r\}italic_B ( italic_p , italic_r ) = { italic_q ∈ bold_X ∣ italic_d ( italic_p , italic_q ) ≤ italic_r }

See Figure 1 for the closed balls in the four mentioned metrics.

A large set of metric spaces have a useful property known as the Heine-Borel property. A metric space is said to have the Heine-Borel property property if every closed and bounded set in the metric space is compact [30].

2.2 The Hilbert and Thompson Metrics, and the Funk Weak Metric

The Hilbert and Thompson metrics, and the Funk weak metrics, are defined on the points in the interior of a convex body ΩΩ\Omegaroman_Ω in dsuperscript𝑑\mathbb{R}^{d}blackboard_R start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT, where ΩΩ\Omegaroman_Ω is a closed and bounded convex set. Unless otherwise stated, we assume that ΩΩ\Omegaroman_Ω is a convex polygon with m𝑚mitalic_m edges. We let ΩΩ\partial\kern 1.0pt\Omega∂ roman_Ω refer to the boundary of ΩΩ\Omegaroman_Ω. Given any pair of points p𝑝pitalic_p and q𝑞qitalic_q within the interior of ΩΩ\Omegaroman_Ω, we define pq¯¯𝑝𝑞\overline{pq}over¯ start_ARG italic_p italic_q end_ARG to be the cord of ΩΩ\Omegaroman_Ω through those two points, and χ(p,q)𝜒𝑝𝑞\chi(p,q)italic_χ ( italic_p , italic_q ) to be the directed ray from p𝑝pitalic_p through q𝑞qitalic_q. Unless otherwise stated, when we take p,qΩ𝑝𝑞Ωp,q\in\Omegaitalic_p , italic_q ∈ roman_Ω, we mean the interior of ΩΩ\Omegaroman_Ω.

Refer to caption

Figure 1: (a) Funk ball, (b) Thompson ball (c) Hilbert ball, and (d) reverse Funk ball.
Definition 2.3 (Funk weak metric).

Given two points p,q𝑝𝑞p,qitalic_p , italic_q in a convex polygon ΩΩ\Omegaroman_Ω in nsuperscript𝑛\mathbb{R}^{n}blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT such that χ(p,q)𝜒𝑝𝑞\chi(p,q)italic_χ ( italic_p , italic_q ) intersects ΩΩ\partial\kern 1.0pt\Omega∂ roman_Ω at a point qsuperscript𝑞q^{\prime}italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT on ΩΩ\partial\kern 1.0pt\Omega∂ roman_Ω, we define the Funk weak metric to be:

FΩ(p,q)=lnpqqq,subscript𝐹Ω𝑝𝑞norm𝑝superscript𝑞norm𝑞superscript𝑞F_{\Omega}(p,q)\leavevmode\nobreak\ =\leavevmode\nobreak\ \ln\frac{\|p-q^{% \prime}\|}{\|q-q^{\prime}\|},italic_F start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_q ) = roman_ln divide start_ARG ∥ italic_p - italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG start_ARG ∥ italic_q - italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG ,

where FΩ(p,q)=0subscript𝐹Ω𝑝𝑞0F_{\Omega}(p,q)=0italic_F start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_q ) = 0.

The above definition is also sometimes called the forward Funk metric. Note however that the Funk weak metric is an asymmetric metric. Therefore, its reverse, the reverse Funk metric, is defined to be rFΩ(p,q)=FΩ(q,p)𝑟subscript𝐹Ω𝑝𝑞subscript𝐹Ω𝑞𝑝rF_{\Omega}(p,q)=F_{\Omega}(q,p)italic_r italic_F start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_q ) = italic_F start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_q , italic_p ). The Hilbert metric can be defined as the average of the forward and reverse Funk metrics.

Definition 2.4 (Hilbert metric).

Given two distinct points p,q𝑝𝑞p,qitalic_p , italic_q in a convex polygon ΩΩ\Omegaroman_Ω in nsuperscript𝑛\mathbb{R}^{n}blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT, let psuperscript𝑝p^{\prime}italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT and qsuperscript𝑞q^{\prime}italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT denote the endpoints of the chord pq¯¯𝑝𝑞\overline{pq}over¯ start_ARG italic_p italic_q end_ARG on (Ω)Ω\partial\kern 1.0pt(\Omega)∂ ( roman_Ω ), so that the points lie in order p,p,q,qsuperscript𝑝𝑝𝑞superscript𝑞\langle p^{\prime},p,q,q^{\prime}\rangle⟨ italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_p , italic_q , italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ⟩, the Hilbert distance between p𝑝pitalic_p and q𝑞qitalic_q, HΩ(p,q)subscript𝐻Ω𝑝𝑞H_{\Omega}(p,q)italic_H start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_q ), is:

HΩ(p,q)=12ln(qppppqqq),subscript𝐻Ω𝑝𝑞12norm𝑞superscript𝑝norm𝑝superscript𝑝norm𝑝superscript𝑞norm𝑞superscript𝑞H_{\Omega}(p,q)\leavevmode\nobreak\ =\leavevmode\nobreak\ \frac{1}{2}\ln\left(% \frac{\|q-p^{\prime}\|}{\|p-p^{\prime}\|}\frac{\|p-q^{\prime}\|}{\|q-q^{\prime% }\|}\right),italic_H start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_q ) = divide start_ARG 1 end_ARG start_ARG 2 end_ARG roman_ln ( divide start_ARG ∥ italic_q - italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG start_ARG ∥ italic_p - italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG divide start_ARG ∥ italic_p - italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG start_ARG ∥ italic_q - italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG ) ,

where HΩ(p,p)=0subscript𝐻Ω𝑝𝑝0H_{\Omega}(p,p)=0italic_H start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_p ) = 0.

Since the product in the definition of the Hilbert metric is the cross ratio, which is preserved under projective transformations, it follows that the Hilbert metric is invariant under projective transformations[19]. It is also worth noting that straight lines are geodesics in the Hilbert metric, though not all geodesics are straight lines, and that the Hilbert metric satisfies all the properties of a metric, including symmetry and the triangle inequality. Nielsen and Shao showed how to algorithmically compute Hilbert balls with the aid of spokes.

Definition 2.5 (Spoke).

Given pΩ𝑝Ωp\in\Omegaitalic_p ∈ roman_Ω a spoke through p𝑝pitalic_p from a vertex v𝑣vitalic_v of ΩΩ\Omegaroman_Ω is pv¯¯𝑝𝑣\overline{pv}over¯ start_ARG italic_p italic_v end_ARG.

Succinctly, a Hilbert Ball of radius r𝑟ritalic_r around a point p𝑝pitalic_p can be constructed by computing all the points that are a distance of r𝑟ritalic_r away from p𝑝pitalic_p along the spokes of p𝑝pitalic_p from the vertices of ΩΩ\Omegaroman_Ω and forming their convex hull (see Figure 1(c)).

Lemma 2.6 (Nielsen and Shao [19]).

Hilbert balls have at most O(m)𝑂𝑚O(m)italic_O ( italic_m ) sides.

The Thompson metric is similar to the Hilbert in its construction. It is the maximum of both the Funk and reverse Funk metrics.

Definition 2.7 (Thompson metric).

Given a bounded closed convex body ΩΩ\Omegaroman_Ω in nsuperscript𝑛\mathbb{R}^{n}blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT and two points p,qint(Ω)𝑝𝑞intΩp,q\in\operatorname{int}(\Omega)italic_p , italic_q ∈ roman_int ( roman_Ω ), let psuperscript𝑝p^{\prime}italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT and qsuperscript𝑞q^{\prime}italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT denote the endpoints of the chord pq¯¯𝑝𝑞\overline{pq}over¯ start_ARG italic_p italic_q end_ARG on (Ω)Ω\partial\kern 1.0pt(\Omega)∂ ( roman_Ω ), so that the points lie in order p,p,q,qsuperscript𝑝𝑝𝑞superscript𝑞\langle p^{\prime},p,q,q^{\prime}\rangle⟨ italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_p , italic_q , italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ⟩, the Thompson distance between p𝑝pitalic_p and q𝑞qitalic_q TΩ(p,q)subscript𝑇Ω𝑝𝑞T_{\Omega}(p,q)italic_T start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_q ) is:

TΩ(p,q)=max(lnqppp,lnpqqq),subscript𝑇Ω𝑝𝑞𝑚𝑎𝑥norm𝑞superscript𝑝norm𝑝superscript𝑝norm𝑝superscript𝑞norm𝑞superscript𝑞T_{\Omega}(p,q)\leavevmode\nobreak\ =\leavevmode\nobreak\ max\left(\ln\frac{\|% q-p^{\prime}\|}{\|p-p^{\prime}\|},\ln\frac{\|p-q^{\prime}\|}{\|q-q^{\prime}\|}% \right),italic_T start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_q ) = italic_m italic_a italic_x ( roman_ln divide start_ARG ∥ italic_q - italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG start_ARG ∥ italic_p - italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG , roman_ln divide start_ARG ∥ italic_p - italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG start_ARG ∥ italic_q - italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∥ end_ARG ) ,

where TΩ(p,p)=0subscript𝑇Ω𝑝𝑝0T_{\Omega}(p,p)=0italic_T start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ( italic_p , italic_p ) = 0.

As it is the maximum of the two Funk weak metrics in a convex body, certain facts about this metric are clear.

Lemma 2.8.

Balls in the Thompson metric are convex polygons with O(m)𝑂𝑚O(m)italic_O ( italic_m ) sides.

Proof 2.9.

This follows directly from the fact that balls in the Funk metric are homotheties (the reverse is flipped) of ΩΩ\Omegaroman_Ω with O(m)𝑂𝑚O(m)italic_O ( italic_m ) sides [22] (see Figure 1 (a) and (d)) and Thompson balls are their intersections (see Figure 1(b)).

Theorem 2.10.

The topology induced by the Hilbert or Thompson metrics, as well as the Funk weak metrics in a bounded convex domain ΩnΩsuperscript𝑛\Omega\subset\mathbb{R}^{n}roman_Ω ⊂ blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT coincides with the Euclidean topology in that domain.

Proof 2.11.

This is a well known result in the field. For reference, this follows from Theorem 2.12.12.12.1 in [3] where the authors show that the Hilbert and Euclidean metrics induce the same topology on bounded convex domains in nsuperscript𝑛\mathbb{R}^{n}blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT. In addition, refer to Proposition 6.16.16.16.1 in "From Funk to Hilbert Geometry" for Funk weak metrics.

From this we immediately have the following corollary.

Corollary 2.12.

Any Hilbert or Thompson metric space, as well as any Funk weak metric space, has the Heine-Borel property.

3 Minimum Enclosing Radius Balls and the Heine-Borel Property

We begin by defining several concepts such as the minimum radius ball, and a set of criterion for being an LP-type problem[26].

Refer to caption

Figure 2: This is 2superscript2\mathbb{R}^{2}blackboard_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT with several closed balls removed. There is no Euclidean minimum radius ball through x,y,z𝑥𝑦𝑧x,y,zitalic_x , italic_y , italic_z.
Definition 3.1 (Minimum Radius enclosing Ball).

A closed ball B𝐵Bitalic_B is a minimum radius enclosing ball of S𝐗𝑆𝐗S\subset\mathbf{X}italic_S ⊂ bold_X when B𝐵Bitalic_B has radius inf{rSB(p,r)}infimumconditional-set𝑟𝑆𝐵𝑝𝑟\inf\{r\mid S\subset B(p,r)\}roman_inf { italic_r ∣ italic_S ⊂ italic_B ( italic_p , italic_r ) }. Note that in a general metric space, the minimum radius enclosing ball need not exist, nor be unique.

Definition 3.2 (Minimum Ball Property).

A metric space (𝐗,d)𝐗𝑑(\mathbf{X},d)( bold_X , italic_d ) satisfies the minimum ball property if for every finite H𝐗𝐻𝐗H\subset\mathbf{X}italic_H ⊂ bold_X there exists a minimum radius enclosing ball of H𝐻Hitalic_H.

For an example of a metric space without the minimum ball property see Figure 2. In this situation several closed disks were removed from 2superscript2\mathbb{R}^{2}blackboard_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT with the Euclidean distance.

Proposition 3.3.

Let (X,d)𝑋𝑑(X,d)( italic_X , italic_d ) be a metric space. If (X,d)𝑋𝑑(X,d)( italic_X , italic_d ) satisfies the Heine-Borel property, then (X,d)𝑋𝑑(X,d)( italic_X , italic_d ) satisfies the minimum ball property.

Proof 3.4.

For a given finite H𝐗𝐻𝐗H\subset\mathbf{X}italic_H ⊂ bold_X define:

S:={s0 there exists a closed ball of radius s enclosing H}andr:=infSformulae-sequenceassign𝑆conditional-set𝑠superscriptabsent0 there exists a closed ball of radius 𝑠 enclosing 𝐻andassign𝑟infimum𝑆S:=\{s\in\mathbb{R}^{\geq 0}\mid\text{ there exists a closed ball of radius }s% \text{ enclosing }H\}\leavevmode\nobreak\ \leavevmode\nobreak\ \text{and}% \leavevmode\nobreak\ \leavevmode\nobreak\ r:=\inf Sitalic_S := { italic_s ∈ blackboard_R start_POSTSUPERSCRIPT ≥ 0 end_POSTSUPERSCRIPT ∣ there exists a closed ball of radius italic_s enclosing italic_H } and italic_r := roman_inf italic_S

Note that S𝑆Sitalic_S is bounded below by 0 and is nonempty because H𝐻Hitalic_H is bounded as a consequence of being finite. So the infimum is defined and finite. Thus, to show there exists a ball of minimal radius, we want to show that there exists a ball of radius r𝑟ritalic_r enclosing H𝐻Hitalic_H. Define the set of centers Cn:={p𝐗:HB(p,r+1/n)}assignsubscript𝐶𝑛conditional-set𝑝𝐗𝐻𝐵𝑝𝑟1𝑛C_{n}:=\{p\in\mathbf{X}:H\subset B(p,r+1/n)\}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT := { italic_p ∈ bold_X : italic_H ⊂ italic_B ( italic_p , italic_r + 1 / italic_n ) }. For all n𝑛n\in\mathbb{N}italic_n ∈ blackboard_N we have CnCn+1subscript𝐶𝑛1subscript𝐶𝑛C_{n}\supset C_{n+1}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ⊃ italic_C start_POSTSUBSCRIPT italic_n + 1 end_POSTSUBSCRIPT because for all pCn+1𝑝subscript𝐶𝑛1p\in C_{n+1}italic_p ∈ italic_C start_POSTSUBSCRIPT italic_n + 1 end_POSTSUBSCRIPT we have HB(p,r+1/(n+1))B(p,r+1/n)𝐻𝐵𝑝𝑟1𝑛1𝐵𝑝𝑟1𝑛H\subset B(p,r+1/(n+1))\subset B(p,r+1/n)italic_H ⊂ italic_B ( italic_p , italic_r + 1 / ( italic_n + 1 ) ) ⊂ italic_B ( italic_p , italic_r + 1 / italic_n ) and thus pCn𝑝subscript𝐶𝑛p\in C_{n}italic_p ∈ italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT as well.

Each Cnsubscript𝐶𝑛C_{n}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is bounded since for any hH𝐻h\in Hitalic_h ∈ italic_H and any pCn𝑝subscript𝐶𝑛p\in C_{n}italic_p ∈ italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT, d(p,h)r+1/n𝑑𝑝𝑟1𝑛d(p,h)\leq r+1/nitalic_d ( italic_p , italic_h ) ≤ italic_r + 1 / italic_n.

We will now show that each Cnsubscript𝐶𝑛C_{n}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is closed. Let hH𝐻h\in Hitalic_h ∈ italic_H and {pi}iCnsubscriptsubscript𝑝𝑖𝑖subscript𝐶𝑛\{p_{i}\}_{i\in\mathbb{N}}\subset C_{n}{ italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT } start_POSTSUBSCRIPT italic_i ∈ blackboard_N end_POSTSUBSCRIPT ⊂ italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT such that pipsubscript𝑝𝑖𝑝p_{i}\rightarrow pitalic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT → italic_p, p𝐗𝑝𝐗p\in\mathbf{X}italic_p ∈ bold_X. Then for any ε>0𝜀0\varepsilon>0italic_ε > 0, there exists i𝑖iitalic_i such that d(p,pi)<ε𝑑𝑝subscript𝑝𝑖𝜀d(p,p_{i})<\varepsilonitalic_d ( italic_p , italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) < italic_ε. So by the triangle inequality, d(p,h)d(p,pi)+d(pi,h)(r+1/n)+ε𝑑𝑝𝑑𝑝subscript𝑝𝑖𝑑subscript𝑝𝑖𝑟1𝑛𝜀d(p,h)\leq d(p,p_{i})+d(p_{i},h)\leq(r+1/n)+\varepsilonitalic_d ( italic_p , italic_h ) ≤ italic_d ( italic_p , italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) + italic_d ( italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_h ) ≤ ( italic_r + 1 / italic_n ) + italic_ε for any ε𝜀\varepsilonitalic_ε, so d(p,h)r+1/n𝑑𝑝𝑟1𝑛d(p,h)\leq r+1/nitalic_d ( italic_p , italic_h ) ≤ italic_r + 1 / italic_n. So hB(p,r+1/n)𝐵𝑝𝑟1𝑛h\in B(p,r+1/n)italic_h ∈ italic_B ( italic_p , italic_r + 1 / italic_n ) for any hH𝐻h\in Hitalic_h ∈ italic_H. Thus pCn𝑝subscript𝐶𝑛p\in C_{n}italic_p ∈ italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT. So Cnsubscript𝐶𝑛C_{n}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT contains all of its limit points and is thus closed.

Define D:=nCnassign𝐷subscript𝑛subscript𝐶𝑛D:=\bigcap_{n}C_{n}italic_D := ⋂ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT. Since we are assuming the Heine-Borel property and showed each Cnsubscript𝐶𝑛C_{n}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is closed and bounded, each Cnsubscript𝐶𝑛C_{n}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is compact. Because {Cn}subscript𝐶𝑛\{C_{n}\}{ italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT } is a decreasing nested sequence of compact sets of 𝐗𝐗\mathbf{X}bold_X with the metric topology, by Cantor’s intersection theorem D𝐷Ditalic_D is nonempty [25].

Fix pD𝑝𝐷p\in Ditalic_p ∈ italic_D. Then for any hH𝐻h\in Hitalic_h ∈ italic_H, d(p,h)r+1/n𝑑𝑝𝑟1𝑛d(p,h)\leq r+1/nitalic_d ( italic_p , italic_h ) ≤ italic_r + 1 / italic_n for all n𝑛n\in\mathbb{N}italic_n ∈ blackboard_N. Thus d(p,h)r𝑑𝑝𝑟d(p,h)\leq ritalic_d ( italic_p , italic_h ) ≤ italic_r. We conclude HB(p,r)𝐻𝐵𝑝𝑟H\subset B(p,r)italic_H ⊂ italic_B ( italic_p , italic_r ).

Note that the converse does not hold. As a counterexample, any infinite set equipped with the discrete metric does not satisfy the Heine-Borel property, yet it satisfies the minimum ball property. Additionally, note that in the above proof of Theorem 3.3 we did not use the symmetric property of the metric space, only the triangle inequality. As such, given the same definition of a ball Definition 2.2 and Definition 3.1 e.g. fixing the direction of the ball, we contribute the following corollary.

Corollary 3.5.

Let (X,d)𝑋𝑑(X,d)( italic_X , italic_d ) be a weak metric space. If (X,d)𝑋𝑑(X,d)( italic_X , italic_d ) satisfies the Heine-Borel property, then (X,d)𝑋𝑑(X,d)( italic_X , italic_d ) satisfies the minimum ball property.

We now introduce the characterization of an LP-type problem by Sharir and Welzl [26].

Definition 3.6 (LP-type).

A pair (H,f)𝐻𝑓(H,f)( italic_H , italic_f ) is called an LP-type problem if H𝐻Hitalic_H is a finite set and f𝑓fitalic_f is a function from subsets of H𝐻Hitalic_H to a totally ordered set such that f𝑓fitalic_f satisfies the following two properties:

Monotonicity:FGHf(F)f(G)f(H).:𝑀𝑜𝑛𝑜𝑡𝑜𝑛𝑖𝑐𝑖𝑡𝑦𝐹𝐺𝐻𝑓𝐹𝑓𝐺𝑓𝐻Monotonicity:F\subset G\subset H\implies f(F)\leq f(G)\leq f(H).italic_M italic_o italic_n italic_o italic_t italic_o italic_n italic_i italic_c italic_i italic_t italic_y : italic_F ⊂ italic_G ⊂ italic_H ⟹ italic_f ( italic_F ) ≤ italic_f ( italic_G ) ≤ italic_f ( italic_H ) .
Locality:FGH,xH,f(F)=f(G)=f(F{x})f(F)=f(G{x}).Locality:F\subset G\subset H,x\in H,f(F)=f(G)=f(F\cup\{x\})\implies f(F)=f(G% \cup\{x\}).italic_L italic_o italic_c italic_a italic_l italic_i italic_t italic_y : italic_F ⊂ italic_G ⊂ italic_H , italic_x ∈ italic_H , italic_f ( italic_F ) = italic_f ( italic_G ) = italic_f ( italic_F ∪ { italic_x } ) ⟹ italic_f ( italic_F ) = italic_f ( italic_G ∪ { italic_x } ) .

Let (𝐗,d)𝐗𝑑(\mathbf{X},d)( bold_X , italic_d ) be a metric space satisfying the minimum ball property. Let 𝐗𝐗\mathbf{X}bold_X have some well-ordering precedes-or-equals\preceq. Take a finite set H𝐗𝐻𝐗H\subset\mathbf{X}italic_H ⊂ bold_X. We define:

f:2H×𝐗,f(G)=inf{(r,p)GB(p,r),ordered lexicographically}:𝑓formulae-sequencesuperscript2𝐻𝐗𝑓𝐺infimumconditional-set𝑟𝑝𝐺𝐵𝑝𝑟ordered lexicographicallyf:2^{H}\rightarrow\mathbb{R}\times\mathbf{X},f(G)=\inf\{(r,p)\mid G\subset B(p% ,r),\text{ordered lexicographically}\}italic_f : 2 start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT → blackboard_R × bold_X , italic_f ( italic_G ) = roman_inf { ( italic_r , italic_p ) ∣ italic_G ⊂ italic_B ( italic_p , italic_r ) , ordered lexicographically }

Here, f(G)=(r^,p^)𝑓𝐺^𝑟^𝑝f(G)=(\hat{r},\hat{p})italic_f ( italic_G ) = ( over^ start_ARG italic_r end_ARG , over^ start_ARG italic_p end_ARG ) gives us a unique minimum radius enclosing ball of G𝐺Gitalic_G, by taking the minimum center with respect to the well-ordering. We call this ball BG=B(p^,r^)subscript𝐵𝐺𝐵^𝑝^𝑟B_{G}=B(\hat{p},\hat{r})italic_B start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT = italic_B ( over^ start_ARG italic_p end_ARG , over^ start_ARG italic_r end_ARG ). Note that if f(G)=(r,p)𝑓𝐺𝑟𝑝f(G)=(r,p)italic_f ( italic_G ) = ( italic_r , italic_p ), there may exist other points psuperscript𝑝p^{\prime}italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT such that GB(r,p)𝐺𝐵𝑟superscript𝑝G\subset B(r,p^{\prime})italic_G ⊂ italic_B ( italic_r , italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) as well; however, f𝑓fitalic_f defines only one minimum radius enclosing ball. Note also that if 𝐗=n𝐗superscript𝑛\mathbf{X}=\mathbb{R}^{n}bold_X = blackboard_R start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT then by Proposition 3.10 we can take the lexicographic usual total order and f𝑓fitalic_f will still be well-defined.

Theorem 3.7.

The problem of finding the minimum radius ball of a set of points in a metric space that satisfies the minimum ball property is LP-type.

Proof 3.8.

To see that f𝑓fitalic_f satisfies monotonicity, consider a finite set H𝐗𝐻𝐗H\subset\mathbf{X}italic_H ⊂ bold_X. Let FGH𝐹𝐺𝐻F\subset G\subset Hitalic_F ⊂ italic_G ⊂ italic_H and suppose f(H)=(r^,p^)𝑓𝐻^𝑟^𝑝f(H)=(\hat{r},\hat{p})italic_f ( italic_H ) = ( over^ start_ARG italic_r end_ARG , over^ start_ARG italic_p end_ARG ). Note that GHBH𝐺𝐻subscript𝐵𝐻G\subset H\subset B_{H}italic_G ⊂ italic_H ⊂ italic_B start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT and thus (r^,p^){(r,p)GB(p,r)}^𝑟^𝑝conditional-set𝑟𝑝𝐺𝐵𝑝𝑟(\hat{r},\hat{p})\in\{(r,p)\mid G\subset B(p,r)\}( over^ start_ARG italic_r end_ARG , over^ start_ARG italic_p end_ARG ) ∈ { ( italic_r , italic_p ) ∣ italic_G ⊂ italic_B ( italic_p , italic_r ) }. It follows from the definition of f𝑓fitalic_f that f(G)f(H)𝑓𝐺𝑓𝐻f(G)\leq f(H)italic_f ( italic_G ) ≤ italic_f ( italic_H ). By the same logic, f(F)f(G)𝑓𝐹𝑓𝐺f(F)\leq f(G)italic_f ( italic_F ) ≤ italic_f ( italic_G ).

To see that f𝑓fitalic_f satisfies locality, consider a finite set H𝐗𝐻𝐗H\subset\mathbf{X}italic_H ⊂ bold_X. Take F,GH𝐹𝐺𝐻F,G\subset Hitalic_F , italic_G ⊂ italic_H and xH𝑥𝐻x\in Hitalic_x ∈ italic_H such that FGH𝐹𝐺𝐻F\subset G\subset Hitalic_F ⊂ italic_G ⊂ italic_H and f(F)=f(G)=f(F{x})𝑓𝐹𝑓𝐺𝑓𝐹𝑥f(F)=f(G)=f(F\cup\{x\})italic_f ( italic_F ) = italic_f ( italic_G ) = italic_f ( italic_F ∪ { italic_x } ). It follows that BF=BG=BF{x}subscript𝐵𝐹subscript𝐵𝐺subscript𝐵𝐹𝑥B_{F}=B_{G}=B_{F\cup\{x\}}italic_B start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT italic_F ∪ { italic_x } end_POSTSUBSCRIPT. Hence, xBG𝑥subscript𝐵𝐺x\in B_{G}italic_x ∈ italic_B start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT and FG{x}BG=BF𝐹𝐺𝑥subscript𝐵𝐺subscript𝐵𝐹F\subset G\cup\{x\}\subset B_{G}=B_{F}italic_F ⊂ italic_G ∪ { italic_x } ⊂ italic_B start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT. Thus, by monotonicity, f(F)f(G{x})f(G)=f(F)𝑓𝐹𝑓𝐺𝑥𝑓𝐺𝑓𝐹f(F)\leq f(G\cup\{x\})\leq f(G)=f(F)italic_f ( italic_F ) ≤ italic_f ( italic_G ∪ { italic_x } ) ≤ italic_f ( italic_G ) = italic_f ( italic_F ), and we conclude that f(F)=f(G{x})𝑓𝐹𝑓𝐺𝑥f(F)=f(G\cup\{x\})italic_f ( italic_F ) = italic_f ( italic_G ∪ { italic_x } ).

Note that the only properties of the metric space used for Theorem 3.7 is the minimum ball property, which we’ve shown can also hold for weak metric spaces in Corollary 3.5, and that if GH𝐺𝐻G\subset Hitalic_G ⊂ italic_H then GBH𝐺subscript𝐵𝐻G\subset B_{H}italic_G ⊂ italic_B start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT, which still holds for weak metrics so long as the direction of the ball is fixed. As such, we contribute the following corollary.

Corollary 3.9.

Finding the minimum radius ball of a set of points in a weak metric space that satisfies the minimum ball property is LP-type.

Note that a well-ordering is not always necessary and can be replaced with a total ordering, such as a lexicographic ordering.

Proposition 3.10.

If X=k𝑋superscript𝑘X=\mathbb{R}^{k}italic_X = blackboard_R start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT for n<𝑛n<\inftyitalic_n < ∞ then the well-order precedes-or-equals\preceq on 𝐗𝐗\mathbf{X}bold_X can be replaced with the lexicographical usual total order \leq from ksuperscript𝑘\mathbb{R}^{k}blackboard_R start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT.

Proof 3.11.

Recall, the following notation and results from Theorem 3.3.

Cn:={p𝐗:HB(p,r+1/n)}andD:=nCnformulae-sequenceassignsubscript𝐶𝑛conditional-set𝑝𝐗𝐻𝐵𝑝𝑟1𝑛andassign𝐷subscript𝑛subscript𝐶𝑛C_{n}:=\{p\in\mathbf{X}:H\subset B(p,r+1/n)\}\qquad\text{and}\qquad D:=\bigcap% _{n}C_{n}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT := { italic_p ∈ bold_X : italic_H ⊂ italic_B ( italic_p , italic_r + 1 / italic_n ) } and italic_D := ⋂ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT

We saw that each Cnsubscript𝐶𝑛C_{n}italic_C start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is closed and bounded with respect to the metric topology. Thus it is compact in the standard (metric) topology in ksuperscript𝑘\mathbb{R}^{k}blackboard_R start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT, and thus D𝐷Ditalic_D, the set of centers of minimum radius enclosing balls is also nonempty and compact in this topology.

Let pi:kd:subscript𝑝𝑖superscript𝑘superscript𝑑p_{i}:\mathbb{R}^{k}\rightarrow\mathbb{R}^{d}italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT : blackboard_R start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT → blackboard_R start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT be the projection onto the ithsuperscript𝑖𝑡i^{th}italic_i start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT component. Since D𝐷Ditalic_D is compact in the standard topology, its projection onto each component is compact and thus compact in the standard topology on \mathbb{R}blackboard_R. We want to show that there is a minimal point with respect to the lexicographic order on ksuperscript𝑘\mathbb{R}^{k}blackboard_R start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT.

The body p1(D)subscript𝑝1𝐷p_{1}(D)italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_D ) is compact, and thus it has a minimal element, x1subscript𝑥1x_{1}italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT.

Then D1:=p11(x1)Dassignsubscript𝐷1superscriptsubscript𝑝11subscript𝑥1𝐷D_{1}:=p_{1}^{-1}(x_{1})\bigcap Ditalic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT := italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ⋂ italic_D is non-empty and compact (in standard topology) since D𝐷Ditalic_D is compact and p11(x1)superscriptsubscript𝑝11subscript𝑥1p_{1}^{-1}(x_{1})italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) is closed.

We similarly recursively define xjsubscript𝑥𝑗x_{j}italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT to be the minimal element of pj(Dj1)subscript𝑝𝑗subscript𝐷𝑗1p_{j}(D_{j-1})italic_p start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_D start_POSTSUBSCRIPT italic_j - 1 end_POSTSUBSCRIPT ) and find Dj:=pj1(xj)Dj1assignsubscript𝐷𝑗superscriptsubscript𝑝𝑗1subscript𝑥𝑗subscript𝐷𝑗1D_{j}:=p_{j}^{-1}(x_{j})\bigcap D_{j-1}italic_D start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT := italic_p start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) ⋂ italic_D start_POSTSUBSCRIPT italic_j - 1 end_POSTSUBSCRIPT to be compact for j{2,..,k}j\in\{2,..,k\}italic_j ∈ { 2 , . . , italic_k }. It follows that (x1,x2,xk)subscript𝑥1subscript𝑥2subscript𝑥𝑘(x_{1},x_{2},...x_{k})( italic_x start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , … italic_x start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is the minimal element of D𝐷Ditalic_D with respect to the lexicographic order.

Corollary 3.12.

The minimum radius ball problem in the Hilbert and Thompson metrics, as well as the Funk weak metrics, is LP-type.

Proof 3.13.

Note that Corollary 2.12, Proposition 3.3, and Theorem 3.7 together imply this result for the Hilbert and Thompson metrics. Note secondly that Corollary 10, Corollary 3.5, and Corollary 3.9 together imply the result for the Funk weak metrics.

4 Hilbert Radius Minimum Enclosing Ball

We have shown in Corollary 3.12 that the minimum radius enclosing ball problem is LP-type in a Hilbert metric, so we will now focus on how to implement the LP algorithm. By Proposition 3.10 the well-order precedes-or-equals\preceq on 𝐗𝐗\mathbf{X}bold_X can be replaced with the lexicographical usual total order \leq from ksuperscript𝑘\mathbb{R}^{k}blackboard_R start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT. This is useful as we can comprehend and compute using the lexicographic total order.

Running an LP algorithm on an LP-type problem requires solving two primitive operations [17]. These are the violation test (given a basis BS𝐵𝑆B\subset Sitalic_B ⊂ italic_S and element xS𝑥𝑆x\in Sitalic_x ∈ italic_S, whether f(B)=f(B{x})𝑓𝐵𝑓𝐵𝑥f(B)=f(B\cup\{x\})italic_f ( italic_B ) = italic_f ( italic_B ∪ { italic_x } )) and basis computations (how to find a basis of B{x}𝐵𝑥B\cup\{x\}italic_B ∪ { italic_x }). Instead of a well-order, we will let δ𝛿\deltaitalic_δ be the lexicographic order 2superscript2\mathbb{R}^{2}blackboard_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for the Hilbert convex polygons. We need the two following supporting lemmas:

Lemma 4.1.

The combinatorial dimension for minimum Hilbert radius balls is 3.

Proof 4.2.

This follows immediately from the fact that, in general position at most three points defines a Hilbert ball [9](see Lemma 14). This is because Hilbert balls intersect along line segment edges in general position. Three points cannot intersect along the same edge at the same distance without two of them lying along the same spoke.

Lemma 4.3.

Given an m𝑚mitalic_m sided convex polygon ΩΩ\Omegaroman_Ω in 2superscript2\mathbb{R}^{2}blackboard_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, the center of the minimum Hilbert radius ball around two points, p,qintΩ𝑝𝑞intΩp,q\in\operatorname{int}\Omegaitalic_p , italic_q ∈ roman_int roman_Ω, can be computed in time O(logm)𝑂𝑚O(\log m)italic_O ( roman_log italic_m ).

Proof 4.4.

The balls of minimum radius for a set of two points are always centered on a section of the bisector in the sector directly between the two points. This is because Hilbert balls are polygons and the balls around two points meet at a segment across this sector[19] (see Lemma 10101010 page 4444). This piece of the bisector is a line [4]. We binary search the boundary of ΩΩ\Omegaroman_Ω to determine which edges define this sector, and then use the bisector equation there [9, 4] (see Section 3 in [9]). We choose the lexicographically smallest point along the bisector in this sector to serve as the center of the minimum Hilbert radius ball around points p,q𝑝𝑞p,qitalic_p , italic_q.

Lemma 4.5.

The violation test in the Hilbert metric for radius balls can be computed in time O(log3m)𝑂superscript3𝑚O(\log^{3}m)italic_O ( roman_log start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_m ).

Proof 4.6.

Given a basis B𝐵Bitalic_B, we are interested in if f(B)=f(B{x})𝑓𝐵𝑓𝐵𝑥f(B)=f(B\cup\{x\})italic_f ( italic_B ) = italic_f ( italic_B ∪ { italic_x } ). We assume that x𝑥xitalic_x is not already in B𝐵Bitalic_B. We have two nontrivial cases based on the size of the basis.

Case 1: B𝐵Bitalic_B has two elements y,z𝑦𝑧y,zitalic_y , italic_z. To determine if f(B)=f(B{x})𝑓𝐵𝑓𝐵𝑥f(B)=f(B\cup\{x\})italic_f ( italic_B ) = italic_f ( italic_B ∪ { italic_x } ), it suffices to check whether x𝑥xitalic_x is contained in the minimum radius ball of y,z𝑦𝑧y,zitalic_y , italic_z. If x𝑥xitalic_x is contained, then f(B)=f(B{x})𝑓𝐵𝑓𝐵𝑥f(B)=f(B\cup\{x\})italic_f ( italic_B ) = italic_f ( italic_B ∪ { italic_x } ). If x𝑥xitalic_x is not contained, either the center must move or the radius must increase, so f(B)<f(B{x})𝑓𝐵𝑓𝐵𝑥f(B)<f(B\cup\{x\})italic_f ( italic_B ) < italic_f ( italic_B ∪ { italic_x } ). Since we can compute the center of a 2222 point Hilbert ball in O(logm)𝑂𝑚O(\log m)italic_O ( roman_log italic_m ) time by Lemma 4.3, we can check this case in O(logm)𝑂𝑚O(\log m)italic_O ( roman_log italic_m ) time.

Case 2: B𝐵Bitalic_B has three elements. In this case, we calculate the center of the Hilbert ball around the three points in time O(log3m)𝑂superscript3𝑚O(\log^{3}m)italic_O ( roman_log start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_m ) using the algorithm from "Delaunay Triangulations in the Hilbert Metric" [9]. If a center exists, we compute the distance between x𝑥xitalic_x and the center. If a center does not exist, one of the three elements is contained in a ball defined by the other two, so we check all pairs of elements using Lemma 4.3 in O(logm)𝑂𝑚O(\log m)italic_O ( roman_log italic_m ) time to find the minimum enclosing ball. We take the resulting ball and check the center’s distance to x𝑥xitalic_x. Either x𝑥xitalic_x is in the ball or it is not, in which case f(B)𝑓𝐵f(B)italic_f ( italic_B ) must increase because either the radius or center must move.

Lemma 4.7.

The basis computation in the Hilbert metric for radius balls can be done in time O(log3m)𝑂superscript3𝑚O(\log^{3}m)italic_O ( roman_log start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_m ).

Proof 4.8.

Suppose we have a previous basis B𝐵Bitalic_B and a new element x𝑥xitalic_x that is not contained inside the ball formed by the basis. We would like to compute the basis of B{x}𝐵𝑥B\cup\{x\}italic_B ∪ { italic_x }. This gives us a few cases. We will consider the non-trivial cases where B𝐵Bitalic_B contains two or three elements. In either situation, we compute the minimum enclosing balls of all pairs and triples of points in B{x}𝐵𝑥B\cup\{x\}italic_B ∪ { italic_x }, and check for containment of the remaining points in those balls. This can be done in O(log3m)𝑂superscript3𝑚O(\log^{3}m)italic_O ( roman_log start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_m ) as described in Case 3 of Lemma 4.5.

Theorem 4.9.

Hilbert minimum radius balls, of n𝑛nitalic_n points, can be computed in time O(nlog3m)𝑂𝑛superscript3𝑚O(n\log^{3}m)italic_O ( italic_n roman_log start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_m ).

Proof 4.10.

This follows directly from the running time of our primitive operations and Theorem 7 and Corollary 8 from Sharir and Welzl’s "A combinatorial bound for linear programming and related problems"[26].

5 Conclusion

In this paper, we contributed a criterion for showing that minimum radius balls are LP-type. This criterion being that finite sets in the metric space always have at least one minimum radius ball. We showed that if a metric space, or weak metric space, has the Heine-Borel property, it has this property. We used this to contribute a minimum radius ball algorithm for the Hilbert metric and proved that the minimum radius ball problem is LP-type for the Thompson metric and Funk weak metrics.

References

  • [1] Ahmed Abdelkader and David M. Mount. Economical Delone sets for approximating convex bodies. In Proc. 16th Scand. Workshop Algorithm Theory, pages 4:1–4:12, 2018. doi:10.4230/LIPIcs.SWAT.2018.4.
  • [2] Ahmed Abdelkader and David M Mount. Convex approximation and the Hilbert geometry. In 2024 Symposium on Simplicity in Algorithms (SOSA), pages 286–298. SIAM, 2024.
  • [3] AF Beardon. The Klein, Hilbert and Poincaré metrics of a domain. Journal of computational and applied mathematics, 105(1-2):155–162, 1999.
  • [4] Madeline Bumpus, Caesar Dai, Auguste H. Gezalyan, Sam Munoz, Renita Santhoshkumar, Songyu Ye, and David M. Mount. Software and analysis for dynamic Voronoi diagrams in the Hilbert metric, 2023. arXiv:2304.02745.
  • [5] Yongxin Chen, Tryphon Georgiou, and Michele Pavon. Entropic and displacement interpolation: A computational approach using the Hilbert metric. SIAM J. Appl. Math., 76:2375–2396, 2016. doi:10.1137/16M1061382.
  • [6] Ştefan Cobzaş and Mircea-Dan Rus. Normal cones and Thompson metric. Topics in Mathematical Analysis and Applications, pages 209–258, 2014.
  • [7] Dmitry Faifman, Constantin Vernicos, and Cormac Walsh. Volume growth of Funk geometry and the flags of polytopes. arXiv preprint arXiv:2306.09268, 2023.
  • [8] Paul Funk. Über geometrien, bei denen die geraden die kürzesten sind. Mathematische Annalen, 101(1):226–237, 1929.
  • [9] Auguste H. Gezalyan, Soo H. Kim, Carlos Lopez, Daniel Skora, Zofia Stefankovic, and David M. Mount. Delaunay triangulations in the hilbert metric. In 19th Scandinavian Symposium and Workshops on Algorithm Theory (SWAT 2024). Schloss Dagstuhl–Leibniz-Zentrum für Informatik, 2024.
  • [10] Auguste H. Gezalyan and David M. Mount. Voronoi diagrams in the Hilbert metric. In 39th International Symposium on Computational Geometry (SoCG 2023). Schloss Dagstuhl-Leibniz-Zentrum für Informatik, 2023.
  • [11] Nir Halman. Simple stochastic games, parity games, mean payoff games and discounted payoff games are all lp-type problems. Algorithmica, 49:37–50, 2007.
  • [12] D. Hilbert. Ueber die gerade Linie als kürzeste Verbindung zweier Punkte. Math. Annalen, 46:91–96, 1895. doi:10.1007/BF02096204.
  • [13] Bas Lemmens and Roger Nussbaum. Birkhoff’s version of Hilbert’s metric and its applications in analysis, 2013. arXiv:1304.7921.
  • [14] Bas Lemmens and Mark Roelands. Unique geodesics for Thompson’s metric. Annales de l’Institut Fourier, 65(1):315–348, 2015.
  • [15] Bas Lemmens and Mark Roelands. Midpoints for Thompson’s metric on symmetric. Osaka Journal of Mathematics, 54(1):197–208, 2017.
  • [16] Yongdo Lim. Geometry of midpoint sets for Thompson’s metric. Linear Algebra and its Applications, 439(1):211–227, 2013.
  • [17] Jiří Matoušek, Micha Sharir, and Emo Welzl. A subexponential bound for linear programming. In Proceedings of the eighth annual symposium on Computational geometry, pages 1–8, 1992.
  • [18] Frank Nielsen and Richard Nock. On the smallest enclosing information disk. Information Processing Letters, 105(3):93–97, 2008.
  • [19] Frank Nielsen and Laetitia Shao. On balls in a Hilbert polygonal geometry. In Proc. 33rd Internat. Sympos. Comput. Geom., pages 67:1–67:4, 2017. (Multimedia contribution). doi:10.4230/LIPIcs.SoCG.2017.67.
  • [20] Frank Nielsen and Ke Sun. Clustering in Hilbert’s projective geometry: The case studies of the probability simplex and the elliptope of correlation matrices. In Frank Nielsen, editor, Geometric Structures of Information, pages 297–331. Springer Internat. Pub., 2019. doi:10.1007/978-3-030-02520-5_11.
  • [21] Frank Nielsen and Ke Sun. Non-linear embeddings in Hilbert simplex geometry. In Topological, Algebraic and Geometric Learning Workshops 2023, pages 254–266. PMLR, 2023.
  • [22] Athanase Papadopoulos and Marc Troyanov. From Funk to Hilbert geometry. In Handbook of Hilbert geometry, volume 22 of IRMA Lectures in Mathematics and Theoretical Physics, pages 33–68. European Mathematical Society Publishing House, 2014. doi:10.4171/147-1/2.
  • [23] Athanase Papadopoulos and Marc Troyanov. Handbook of Hilbert geometry, volume 22 of IRMA Lectures in Mathematics and Theoretical Physics. European Mathematical Society Publishing House, 2014. doi:10.4171/147.
  • [24] David Reeb, Michael J. Kastoryano, and Michael M. Wolf. Hilbert’s projective metric in quantum information theory. J. Math. Physics, 52(8), 2011. doi:10.1063/1.3615729.
  • [25] Walter Rudin et al. Principles of mathematical analysis, volume 3. McGraw-hill New York, 1964.
  • [26] Micha Sharir and Emo Welzl. A combinatorial bound for linear programming and related problems. In STACS 92: 9th Annual Symposium on Theoretical Aspects of Computer Science Cachan, France, February 13–15, 1992 Proceedings 9, pages 567–579. Springer, 1992.
  • [27] Anthony C Thompson. On certain contraction mappings in a partially ordered vector space. Proceedings of the American Mathematical Society, 14(3):438–443, 1963.
  • [28] Marc Troyanov. Funk and Hilbert geometries from the Finslerian viewpoint. arXiv preprint arXiv:1311.2508, 2013.
  • [29] Constantin Vernicos and Cormac Walsh. Flag-approximability of convex bodies and volume growth of Hilbert geometries, 2018. arXiv:1809.09471.
  • [30] Robert Williamson and Ludvik Janos. Constructing metrics with the Heine-Borel property. Proceedings of the American Mathematical Society, 100(3):567–573, 1987.