Deterministic metric $1$-median selection with very few queries ¹¹1Part of this paper appears in Proceedings of the 27th International Computing and Combinatorics Conference (COCOON 2021).

An $n$-point metric space $(M,d)$ is a size-$n$ set $M$ endowed with a distance function $d\colon M\times M\to[0,\infty)$ such that

• •

  $d(x,y)=0$ if and only if $x=y$,

• •

  $d(x,y)=d(y,x)$, and

• •

  $d(x,y)+d(y,z)\geq d(x,z)$ (triangle inequality)

for all $x$, $y$, $z\in M$ [16]. Metric $1$-median asks for a point $p\in M$ minimizing $\sum_{x\in M}\,d(p,x)$. Clearly, it has a brute-force $O(n^{2})$-time algorithm. Furthermore, it generalizes the classical median selection [6] and can be generalized further to metric $k$-median clustering. In social network analysis, metric $1$-median asks for an actor with the maximum closeness centrality [17]. For all $\beta\geq 1$, a $\beta$-approximate $1$-median of $(M,d)$ is a point $p\in M$ satisfying $\sum_{y\in M}\,d(p,y)\leq\beta\cdot\min_{q\in M}\sum_{y\in M}\,d(q,y)$. By convention, a $\beta$-approximation algorithm for metric $1$-median must output a $\beta$-approximate $1$-median of $(M,d)$. A query inspects $d(x,y)$ for some $x$, $y\in M$. An algorithm is nonadaptive if its $i$th query $(x_{i},y_{i})\in M^{2}$ is independent of the answers to the first $i-1$ queries, for all $i>1$. Write $d_{G}$ for the distance function induced by an undirected graph $G$.

Indyk [11, 12] gives a Monte Carlo $O(n/\epsilon^{2})$-time $(1+\epsilon)$-approximation algorithm for metric $1$-median, where $\epsilon>0$. His time complexity is optimal w.r.t. $n$. When restricted to $\mathbb{R}^{D}$, metric $1$-median has a Monte Carlo $O(D\cdot\exp(\text{poly}(1/\epsilon)))$-time $(1+\epsilon)$-approximation algorithm [14]. The more general $k$-median clustering in metric spaces has streaming approximation algorithms [10], requires $\Omega(nk)$ time for $O(1)$-approximations [15] and is inapproximable to within $(1+2/e-\Omega(1))$ unless $\text{NP}\subseteq\text{DTIME}(n^{O(\log\log n)})$ [13]. For $\mathbb{R}^{D}$ and graph metrics, a well-studied problem is to find the average distance from a query point to a finite set of points [1, 8, 9].

Deterministic $\omega(n)$-query computation is almost completely understood for metric $1$-median: For all constants $\epsilon\in(0,1)$, the best approximation ratio achievable by deterministic $o(n^{2})$-query and $O(n^{1+\epsilon})$-query algorithms is $4$ and $2\lceil 1/\epsilon\rceil$, respectively [2, 4, 18]. The same holds with “query” replaced by “time” and regardless of whether the algorithms can be adaptive [2, 4]. In contrast, we study the largely unknown deterministic $O(n)$- or $o(n)$-query computation. An $o(n)$-query algorithm enjoys the strength of ignoring a $1-o(1)$ fraction of points.

It is folklore that every point is an $(n-1)$-approximate $1$-median. Surprisingly, this is the current best upper bound for deterministic $o(n)$-query algorithms. In particular, no deterministic $o(n)$-query $o(n)$-approximation algorithms are known for metric $1$-median. Instead, we give a deterministic, $o(n)$-query, $o(f(n)\cdot\log n)$-approximation and nonadaptive algorithm for each computable function $f\colon\mathbb{Z}^{+}\to\mathbb{Z}^{+}$ satisfying $f(n)=\omega(1)$. So, e.g., metric $1$-median has a deterministic $o(n)$-query $o(\alpha(n)\cdot\log n)$-approximation algorithm for the very slowly growing inverse Ackermann function $\alpha(\cdot)$. Our main technical discovery is that a $\beta$-approximate $1$-median of $(S,d|_{S\times S})$ (where $d|_{S\times S}$ denotes $d$ restricted to $S\times S$) is an $O(\beta n/|S|)$-approximate $1$-median of $(M,d)$, for all $\emptyset\subsetneq S\subseteq M$ and $\beta\geq 1$. When $S\subseteq M$ is a uniformly random set of a sufficiently large size, an approximate solution to metric $k$-median clustering for $(S,d|_{S\times S})$ is a good one for $(M,d)$ with high probability [7]. But our discovery is for any $S$ and is new.

Chang [3] shows that metric $1$-median has a deterministic, $O(\exp(O(1/\epsilon))\cdot n\log n)$-time, $O(\exp(O(1/\epsilon))\cdot n)$-query, $(\epsilon\log n)$-approximation and nonadaptive algorithm, for all $\epsilon>0$. So deterministic $O(n)$-query algorithms can be $(\epsilon\log n)$-approximate for each $\epsilon>0$. Currently, the best lower bound against deterministic $O(n)$-query algorithms is that they cannot be $O(1)$-approximate [4]. So there is a huge gap between Chang’s [3] approximation ratio of $\epsilon\log n$ and the current best lower bound. We close the gap by showing each deterministic $O(n)$-query algorithm for metric $1$-median to be not $(\delta\log n)$-approximate for a sufficiently small constant $\delta>0$ (depending on the algorithm). Our approach, sketched below, adversarially answers the queries of a deterministic $O(n)$-query algorithm Alg:

1. (I)

   Start with the complete graph on $M$.

2. (II)

   Mark all edges in an $O(1)$-regular expander graph as permanent.

3. (III)

   Repeat the following:

   1. (1)

      Upon receiving a query $(a,b)\in M^{2}$, find a shortest $a$-$b$ path $P$ and answer by the length of $P$.

   2. (2)

      Mark all edges of $P$ as permanent.

   3. (3)

      For each vertex $v$ incident to too many permanent edges, remove all non-permanent edges incident to $v$.

Intuitively, item (III3) keeps degrees small, thus forcing the output of Alg to have a large average distance to other points. Because item (III1) answers a query by the length of $P$, items (III2)–(III3) must preserve all edge of $P$ (by marking them as permanent and not removing them) for the consistency in answering future queries. Items (I) and (III1)–(III3) follow Chang’s [4] paradigm. To prove a lower bound against Alg, we shall make the output of Alg a lot worse than a $1$-median, presumably by identifying or planting a vertex with a sufficiently small average distance to other points. However, Chang fails in this respect. We overcome his problem by item (II), which allows a vertex to have an $O(1)$ average distance to other vertices.

An extension of our lower bound forbids each deterministic $o(n)$-query algorithm for metric $1$-median to be $o(f(n)\cdot\log n)$-approximate for some computable function $f\colon\mathbb{Z}^{+}\to\mathbb{Z}^{+}$ satisfying $f(n)=\omega(1)$. In particular, deterministic $o(n)$-query $O(\log n)$-approximation algorithms do not exist. Previously, the best lower bound against deterministic $o(n)$-query algorithms $A$ is folklore and forbids $A$ to be $h_{A}(n)$-approximate for some $h_{A}(n)=\omega(1)$.³³3For a sketch of proof, answer all queries of $A$ by $1$ and put all points not involved in the queries to be extremely close to one another but extremely far away from $A$’s output and from the points involved in the queries. So previous works do not yet refute the existence of deterministic $o(n)$-query $O(\alpha(n))$-approximation algorithms, where $\alpha(\cdot)$ is the very slowly growing inverse Ackermann function.

Chang [5]’s adversarial method shows that metric $1$-median has no deterministic $O(n)$-query $o(\log n)$-approximation algorithms that make each point involve in $O(1)$ queries to $d$. But his adversary is rather naïve and does not seem to yield any unconditional lower bound such as ours.

Take an $n$-point metric space $(M,d)$ and $\emptyset\subsetneq S\subseteq M$. Define

to be a $1$-median of $(M,d)$ and $(S,d|_{S\times S})$, respectively, breaking ties arbitrarily. Furthermore, pick $\boldsymbol{u}$ and $\boldsymbol{v}$ independently and uniformly at random from $S$. So

is the average distance in $(S,d|_{S\times S})$.

For all $x^{\prime}_{S}\in S$,

The next two lemmas constitute our main discovery.

Lemmas 3–4 imply the following.

The following theorem is due to Chang [3].

Below is our main theorem.

Taking a very slowly growing $f(\cdot)$ (e.g., the iterated logarithm or the inverse Ackermann function), Theorem 7 allows deterministic $o(n)$-query algorithms to be very close to being $O(\log n)$-approximate.

Fix any deterministic $q$-query algorithm Alg, where $q=q(n)=O(n)$. Then take a constant $C>2d+4q/n$, where $d=O(1)$ is such that $d$-regular expander graphs exist. By padding, assume the number of Alg’s queries to be exactly $q$. Adversary Adv in Fig. 1 answers the queries of Alg. All graphs are assumed to be undirected.

As a remark, whenever an edge of a graph is marked as permanent, that edge is considered to be permanent in all graphs. For example, an edge of $G^{\text{exp}}$ marked as permanent in line 3 of Adv is considered to be permanent in lines 11–13, even though the latter processes $G^{(i)}$ rather than $G^{\text{exp}}$. Similarly, although an edge marked as permanent by line 8 comes from $G^{(i-1)}$ by line 6, it is considered to be permanent in lines 11–13 as well.