Metric random matchings with applications¹¹1A preliminary version of this paper appears in https://arxiv.org/abs/1702.03106.

A metric space is a nonempty set $M$ endowed with a metric, i.e., a function $d\colon M\times M\to[\,0,\infty\,)$ such that

• •

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

• •

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

• •

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

for all $x$, $y$, $z\in M$ [14].

For all $n\in\mathbb{Z}^{+}$, define $[n]\equiv\{1,2,\ldots,n\}$. Given $n\in\mathbb{Z}^{+}$ and oracle access to a metric $d\colon[n]\times[n]\to[\,0,\infty\,)$, metric $1$-median asks for $\mathop{\mathrm{argmin}}_{y\in[n]}\,\sum_{x\in[n]}\,d(y,x)$, breaking ties arbitrarily. It generalizes the classical median selection on the real line and has a brute-force $\Theta(n^{2})$-time algorithm. More generally, metric $k$-median asks for $c_{1}$, $c_{2}$, $\ldots$, $c_{k}\in[n]$ minimizing $\sum_{x\in[n]}\,\min_{i=1}^{k}\,d(x,c_{i})$. Because $d(\cdot,\cdot)$ defines $\binom{n}{2}=\Theta(n^{2})$ nonzero distances, only $o(n^{2})$-time algorithms are said to run in sublinear time [9]. For all $\alpha\geq 1$, an $\alpha$-approximate $1$-median is a point $p\in[n]$ satisfying

For all $\epsilon>0$, metric $1$-median has a Monte Carlo $(1+\epsilon)$-approximation $O(n/\epsilon^{2})$-time algorithm [9, 10]. Guha et al. [8] show that metric $k$-median has a Monte Carlo, $O(\exp(O(1/\epsilon)))$-approximation, $O(nk\log n)$-time, $O(n^{\epsilon})$-space and one-pass algorithm for all small $k$ as well as a deterministic, $O(\exp(O(1/\epsilon)))$-approximation, $O(n^{1+\epsilon})$-time, $O(n^{\epsilon})$-space and one-pass algorithm. Given $n$ points in $\mathbb{R}^{D}$ with $D\geq 1$, the Monte Carlo algorithms of Kumar et al. [11] find a $(1+\epsilon)$-approximate $1$-median in $O(D\cdot\exp(1/\epsilon^{O(1)}))$ time and a $(1+\epsilon)$-approximate solution to metric $k$-median in $O(Dn\cdot\exp((k/\epsilon)^{O(1)}))$ time. All randomized $O(1)$-approximation algorithms for metric $k$-median take $\Omega(nk)$ time [12, 8]. Chang [3] shows that metric $1$-median has a deterministic, $(2h)$-approximation, $O(hn^{1+1/h})$-time and nonadaptive algorithm for all constants $h\in\mathbb{Z}^{+}\setminus\{1\}$, generalizing the results of Chang [2] and Wu [16]. On the other hand, he disproves the existence of deterministic $(2h-\epsilon)$-approximation $O(n^{1+1/(h-1)}/h)$-time algorithms for all constants $h\in\mathbb{Z}^{+}\setminus\{1\}$ and $\epsilon>0$ [4, 5].

In social network analysis, the closeness centrality of a point $v$ is the reciprocal of the average distance from $v$ to all points [15]. So metric $1$-median asks for a point with the maximum closeness centrality. Given oracle access to a graph metric, the Monte-Carlo algorithms of Goldreich and Ron [7] and Eppstein and Wang [6] estimate the closeness centrality of a given point and those of all points, respectively.

All known sublinear-time algorithms for metric $1$-median are either deterministic or Monte Carlo, the latter having a positive probability of failure. For example, Indyk’s Monte Carlo $(1+\epsilon)$-approximation algorithm outputs with a positive probability a solution without approximation guarantees. In contrast, we show that metric $1$-median has a randomized algorithm that always outputs a $(2+\epsilon)$-approximate solution in expected $O(n/\epsilon^{2})$ time for all $\epsilon\in(0,1)$. So, excluding the known deterministic algorithms (which are Las Vegas only in the degenerate sense), this paper gives the first Las Vegas approximation algorithm for metric $1$-median with an expected sublinear running time. Note that deterministic sublinear-time algorithms for metric $1$-median can be $4$-approximate but not $(4-\epsilon)$-approximate for any constant $\epsilon>0$ [2, 5]. So our approximation ratio of $2+\epsilon$ beats that of any deterministic sublinear-time algorithm. Inheriting Indyk’s algorithm, our algorithm outputs a $(1+\epsilon)$-approximate $1$-median in $O(n/\epsilon^{2})$ time with probability $\Omega(1)$ for all $\epsilon\in(0,1)$.

Indyk [9, 10] gives a Monte-Carlo $O(n/\epsilon^{3.5})$-time algorithm that approximates the average distance in any metric space $([n],d)$ to within a multiplicative factor in $[\,1-\epsilon,1+\epsilon\,]$, for all $\epsilon>0$. Barhum, Goldreich and Shraibman [1] improve Indyk’s time complexity of $O(n/\epsilon^{3.5})$ to $O(n/\epsilon^{2})$. This paper gives a Monte-Carlo $O(n/\epsilon)$-time algorithm that approximates the average distance in $([n],d)$ to within a multiplicative factor in $[\,1/2-\epsilon,1\,]$, for all $\epsilon>0$. For all $\epsilon=\omega(1/n^{1/4})$, we present a Monte-Carlo $O(n)$-time algorithm approximating the average distance of any graph metric to within a multiplicative factor in $[\,1-\epsilon,1+\epsilon\,]$. But for general metrics, we do not know whether the $O(n/\epsilon^{2})$ running time of Barhum, Goldreich and Shraibman can be improved to $O(n/\epsilon^{2-\Omega(1)})$.

For a metric space $([n],d)$,

breaking ties arbitrarily in equation (2). So $\bar{r}$ is the average distance in $([n],d)$, and $p^{*}$ is a $1$-median.

An algorithm $A$ with oracle access to $d\colon[n]\times[n]\to[\,0,\infty\,)$ is denoted by $A^{d}$ and may query $d$ on any $(x,y)\in[n]\times[n]$ for $d(x,y)$. In this paper, all Landau symbols (such as $O(\cdot)$, $o(\cdot)$, $\Omega(\cdot)$ and $\omega(\cdot)$) are w.r.t. $n$. The following result is due to Indyk.

Henceforth, denote Indyk’s algorithm in Fact 1 by Indyk median. It is given $n\in\mathbb{Z}^{+}$, $\epsilon>0$ and oracle access to a metric $d\colon[n]\times[n]\to[\,0,\infty\,)$. The following fact on estimating the average distance is due to Barhum, Goldreich and Shraibman.

This section presents a randomized algorithm that always outputs a $(2+\epsilon)$-approximate $1$-median, where $\epsilon\in(0,1)$. Clearly,

For each permutation $\pi\colon[n]\to[n]$,

where the first and the second inequalities follow from the triangle inequality and the injectivity of $\pi$.

Inequalities (3)–(4) yield the following.

By Fact 1, Indyk median satisfies condition (2) in Theorem 8. But it does not satisfy condition (1).

We now justify the optimality of the ratio of $2+\epsilon$ in Theorem 8. Let $A$ be a randomized algorithm that always outputs a $(2-\epsilon)$-approximate $1$-median. Furthermore, denote by $p\in[n]$ (resp., $Q\subseteq[n]\times[n]$) the output (resp., the set of queries as unordered pairs) of $A^{d_{1}}(n)$, where $d_{1}$ is the discrete metric (i.e., $d_{1}(x,y)=1$ and $d_{1}(x,x)=0$ for all distinct $x$, $y\in[n]$). Without loss of generality, assume $(p,y)\in Q$ for all $y\in[n]\setminus\{p\}$ by adding dummy queries. So the queries in $Q$ witness that

Assume without loss of generality that $A$ never queries for the distance from a point to itself.

In the sequel, consider the case that $|Q|<\epsilon\cdot(n-1)^{2}/8$. By the averaging argument, there exists a point $\hat{p}\in[n]\setminus\{p\}$ involved in at most $2\cdot|Q|/(n-1)$ queries in $Q$ (note that each query involves two points). Because every function $f\colon[n]\times[n]\to[\,0,\infty\,)$ with