Tighter Estimates for epsilon-nets for Disks

TL;DR

This paper improves the upper bounds on the size of epsilon-nets for disks in the plane, leading to better approximation algorithms for the geometric hitting set problem with practical implementations.

Contribution

It presents a new, constructive algorithm using Delaunay triangulations that reduces the epsilon-net size bound from 24/epsilon to 13.4/epsilon for disks.

Findings

Improved epsilon-net size bound to 13.4/epsilon for disks.

Implemented and publicly released the new algorithm.

Experimental results indicate actual epsilon-net sizes are around 9/epsilon.

Abstract

The geometric hitting set problem is one of the basic geometric combinatorial optimization problems: given a set $P$ of points, and a set $\mathcal{D}$ of geometric objects in the plane, the goal is to compute a small-sized subset of $P$ that hits all objects in $\mathcal{D}$. In 1994, Bronniman and Goodrich made an important connection of this problem to the size of fundamental combinatorial structures called $\epsilon$-nets, showing that small-sized $\epsilon$-nets imply approximation algorithms with correspondingly small approximation ratios. Very recently, Agarwal and Pan showed that their scheme can be implemented in near-linear time for disks in the plane. Altogether this gives $O(1)$-factor approximation algorithms in $\tilde{O}(n)$ time for hitting sets for disks in the plane. This constant factor depends on the sizes of $\epsilon$-nets for disks; unfortunately, the current state-of-the-art bounds are large -- at least $24/\epsilon$ and most likely larger than $40/\epsilon$. Thus the approximation factor of the Agarwal and Pan algorithm ends up being more than $40$. The best lower-bound is $2/\epsilon$, which follows from the Pach-Woeginger construction for halfspaces in two dimensions. Thus there is a large gap between the best-known upper and lower bounds. Besides being of independent interest, finding precise bounds is important since this immediately implies an improved linear-time algorithm for the hitting-set problem. The main goal of this paper is to improve the upper-bound to $13.4/\epsilon$ for disks in the plane. The proof is constructive, giving a simple algorithm that uses only Delaunay triangulations. We have implemented the algorithm, which is available as a public open-source module. Experimental results show that the sizes of $\epsilon$-nets for a variety of data-sets is lower, around $9/\epsilon$.