A deterministic sublinear-time nonadaptive algorithm for metric $1$ -median selection

Ching-Lueh Chang ¹¹1Department of Computer Science and Engineering, Yuan Ze University, Taoyuan, Taiwan. Email: clchang@saturn.yzu.edu.tw ²²2Innovation Center for Big Data and Digital Convergence, Yuan Ze University, Taoyuan, Taiwan. ³³3Supported in part by the Ministry of Science and Technology of Taiwan under grant 103-2221-E-155-026-MY2.

Abstract

We give a deterministic $O(hn^{1+1/h})$ -time $(2h)$ -approximation nonadaptive algorithm for $1$ -median selection in $n$ -point metric spaces, where $h\in\mathbb{Z}^{+}\setminus\{1\}$ is arbitrary. Our proof generalizes that of Chang [2].

1 Introduction

A metric space $(M,d)$ is a nonempty set $M$ endowed with a function $d\colon M\times M\to[\,0,\infty\,)$ such that for all $x$ , $y$ , $z\in M$ ,

•
$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).

The metric $1$ -median problem asks for a point in an $n$ -point metric space $(M,d)$ with the minimum average distance to other points. For $c\geq 1$ , a point $\hat{p}\in M$ is said to be $c$ -approximate for metric $1$ -median if

\sum_{x\in M}\,d\left(\hat{p},x\right)\leq c\cdot\min_{p\in M}\,\sum_{x\in M}\,d\left(p,x\right).

An algorithm for metric $1$ -median is nonadaptive if the sequence of distances that it inspects depends only on $M$ but not on $d$ . Because there are $n(n-1)/2$ nonzero distances, “sublinear-time” will mean “ $o(n^{2})$ -time.”

Indyk [5, 6] shows that metric $1$ -median has a Monte-Carlo $O(n/\epsilon^{2})$ -time $(1+\epsilon)$ -approximation algorithm for each $\epsilon>0$ . In $\mathbb{R}^{D}$ , where $D\geq 1$ , metric $1$ -median has a Monte-Carlo $O(2^{\text{poly}(1/\epsilon)}D)$ -time $(1+\epsilon)$ -approximation algorithm for all $\epsilon>0$ [7]. Many other algorithms are known for $k$ -median selection [4, 7, 1]. For example, Guha et al. [4] give a deterministic, $O(n^{1+\epsilon})$ -time, $O(n^{\epsilon})$ -space, $2^{O(1/\epsilon)}$ -approximation and one-pass algorithm as well as a Monte-Carlo algorithm for $k$ -median selection in metric spaces, where $\epsilon>0$ .

We show that metric $1$ -median has a deterministic $O(hn^{1+1/h})$ -time $(2h)$ -approximation nonadaptive algorithm for all $h\in\mathbb{Z}^{+}\setminus\{1\}$ , generalizing the following theorems:

Theorem 1 ([2]).

Metric $1$ -median has a deterministic $O(n^{1.5})$ -time $4$ -approximation nonadaptive algorithm.

Theorem 2 ([8]).

For each $h\in\mathbb{Z}^{+}\setminus\{1\}$ , metric $1$ -median has a deterministic $O(hn^{1+1/h})$ -time $(2h)$ -approximation (adaptive) algorithm.⁴⁴4The time complexity of $O(hn^{1+1/h})$ is originally presented as $O(n^{1+1/h})$ because $h$ is independent of $n$ . We include the $O(h)$ factor, which is implicit in the original proof, for ease of comparison.

When $n$ is a perfect square and $h=2$ , our proof is equivalent to that of Theorem 1 [2]. Chang [3] shows that metric $1$ -median has no deterministic $o(n^{2})$ -query $(4-\Omega(1))$ -approximation algorithms (where an algorithm’s query complexity is the number of distances that it inspects). So the approximation ratio of $4$ in Theorem 1 cannot be improved to a smaller constant.

2 Our algorithm

Let $(\{0,1,\ldots,n-1\},d)$ be a metric space, $h\in\mathbb{Z}^{+}\setminus\{1\}$ and $t=\lceil n^{1/h}\rceil$ . For all $j\in\{0,1,\ldots,n-1\}$ , denote the (unique) $t$ -ary representation of $j$ by

\displaystyle\left(s_{h-1}(j),s_{h-2}(j),\ldots,s_{0}(j)\right)\in\left\{0,1,\ldots,t-1\right\}^{h},

i.e.,

\displaystyle\sum_{\ell=0}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}=j.

(1)

For all $i$ , $j\in\{0,1,\ldots,n-1\}$ ,

\displaystyle\tilde{d}\left(i,i+j\bmod{n}\right)\stackrel{{\scriptstyle\text{def.}}}{{=}}\sum_{k=0}^{h-1}\,d\left(i+\sum_{\ell=h-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n},i+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right).

(2)

By convention, empty sums vanish, e.g., $\sum_{\ell=h}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}=0$ .

Lemma 3.

For all $i$ , $j\in\{0,1,\ldots,n-1\}$ ,

d\left(i,i+j\bmod{n}\right)\leq\tilde{d}\left(i,i+j\bmod{n}\right).

Proof.

By equation (2) and the triangle inequality for $d$ ,

\tilde{d}\left(i,i+j\bmod{n}\right)\geq d\left(i,i+\sum_{\ell=0}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right).

This and equation (1) complete the proof. ∎

Lemma 4.

For all $\alpha\in\{0,1,\ldots,n-1\}$ with

\displaystyle\sum_{j=0}^{n-1}\,\tilde{d}\left(\alpha,\alpha+j\bmod{n}\right)=\min_{i=0}^{n-1}\,\sum_{j=0}^{n-1}\,\tilde{d}\left(i,i+j\bmod{n}\right),

(3)

we have

\displaystyle\sum_{j=0}^{n-1}\,d\left(\alpha,j\right)\leq 2h\cdot\min_{i=0}^{n-1}\,\sum_{j=0}^{n-1}\,d\left(i,j\right).

(4)

Proof.

Let $\boldsymbol{u}$ be a uniformly random element of $\{0,1,\ldots,n-1\}$ and

\displaystyle i^{\prime}=\mathop{\mathrm{argmin}}_{i=0}^{n-1}\,\sum_{j=0}^{n-1}\,d\left(i,j\right),

(5)

breaking ties arbitrarily. It is easy to see that

\displaystyle\sum_{j=0}^{n-1}\,d\left(\alpha,j\right)=\sum_{j=0}^{n-1}\,d\left(\alpha,\alpha+j\bmod{n}\right).

Furthermore,

	$\displaystyle\sum_{j=0}^{n-1}\,d\left(\alpha,\alpha+j\bmod{n}\right)$	$\displaystyle\stackrel{{\scriptstyle\text{Lemma~{}\ref{pseudodistanceislarger}}}}{{\leq}}$	$\displaystyle\sum_{j=0}^{n-1}\,\tilde{d}\left(\alpha,\alpha+j\bmod{n}\right)$
		$\displaystyle\stackrel{{\scriptstyle\text{equation~{}(\ref{theoptimalpointwithrespecttopseudodistance})}}}{{\leq}}$	$\displaystyle\mathop{\mathrm{E}}\left[\,\sum_{j=0}^{n-1}\,\tilde{d}\left(\boldsymbol{u},\boldsymbol{u}+j\bmod{n}\right)\,\right].$

By equation (2),

			$\displaystyle\mathop{\mathrm{E}}\left[\,\sum_{j=0}^{n-1}\,\tilde{d}\left(\boldsymbol{u},\boldsymbol{u}+j\bmod{n}\right)\,\right]$
		$\displaystyle=$	$\displaystyle\mathop{\mathrm{E}}\left[\,\sum_{j=0}^{n-1}\,\sum_{k=0}^{h-1}\,d\left(\boldsymbol{u}+\sum_{\ell=h-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n},\boldsymbol{u}+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right)\,\right].$

Finally,

	$\displaystyle\mathop{\mathrm{E}}\left[\,\sum_{j=0}^{n-1}\,\sum_{k=0}^{h-1}\,d\left(\boldsymbol{u}+\sum_{\ell=h-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n},\boldsymbol{u}+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right)\,\right]$	(6)
$\displaystyle\leq$	$\displaystyle\mathop{\mathrm{E}}\left[\,\sum_{j=0}^{n-1}\,\sum_{k=0}^{h-1}\,d\left(i^{\prime},\boldsymbol{u}+\sum_{\ell=h-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right)+d\left(i^{\prime},\boldsymbol{u}+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right)\,\right]$
$\displaystyle=$	$\displaystyle\sum_{j=0}^{n-1}\,\sum_{k=0}^{h-1}\,\mathop{\mathrm{E}}\left[\,d\left(i^{\prime},\boldsymbol{u}+\sum_{\ell=h-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right)\,\right]+\mathop{\mathrm{E}}\left[\,d\left(i^{\prime},\boldsymbol{u}+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right)\,\right]$
$\displaystyle=$	$\displaystyle\sum_{j=0}^{n-1}\,\sum_{k=0}^{h-1}\,\left(\mathop{\mathrm{E}}\left[\,d\left(i^{\prime},\boldsymbol{u}\right)\,\right]+\mathop{\mathrm{E}}\left[\,d\left(i^{\prime},\boldsymbol{u}\right)\,\right]\right)$
$\displaystyle=$	$\displaystyle 2nh\cdot\mathop{\mathrm{E}}\left[\,d\left(i^{\prime},\boldsymbol{u}\right)\,\right],$

where the inequality follows from the triangle inequality for $d$ , and the second-to-last equality is true because $\boldsymbol{u}+\sum_{\ell=h-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}$ distributes uniformly at random over $\{0,1,\ldots,n-1\}$ for any $j\in\{0,1,\ldots,n-1\}$ and $k\in\{0,1,\ldots,h\}$ . Inequalities (5)–(6) imply inequality (4). ∎

For a predicate $P$ , let $\chi[\,P\,]=1$ if $P$ is true and $\chi[\,P\,]=0$ otherwise. Define $(s^{\prime}_{h-1},s^{\prime}_{h-2},\ldots,s^{\prime}_{0})\in\{0,1,\ldots,t-1\}^{h}$ to be the $t$ -ary representation of $n-1$ . So $\sum_{r=0}^{h-1}\,s^{\prime}_{r}\cdot t^{r}=n-1$ . For $i\in\{0,1,\ldots,n-1\}$ and $m\in\{0,1,\ldots,h-1\}$ ,

$\displaystyle f\left(i,m\right)$	$\displaystyle\stackrel{{\scriptstyle\text{def.}}}{{=}}$	$\displaystyle\sum_{s_{m},s_{m-1},\ldots,s_{0}=0}^{t-1}\,\chi\left[\sum_{r=0}^{m}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}\right]$	(7)
	$\displaystyle\cdot$	$\displaystyle\sum_{k=0}^{m}\,d\left(i+\sum_{\ell=m+1-k}^{m}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+\sum_{\ell=m-k}^{m}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right),$	(7)
$\displaystyle g\left(i,m\right)$	$\displaystyle\stackrel{{\scriptstyle\text{def.}}}{{=}}$	$\displaystyle\sum_{s_{m},s_{m-1},\ldots,s_{0}=0}^{t-1}\,\sum_{k=0}^{m}\,d\left(i+\sum_{\ell=m+1-k}^{m}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+\sum_{\ell=m-k}^{m}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right).\,\,\,\,\,\,\,\,\,\,$	(8)

Clearly,

	$\displaystyle f\left(i,0\right)$	$\displaystyle=$	$\displaystyle\sum_{s_{0}=0}^{s^{\prime}_{0}}\,d\left(i,i+s_{0}\bmod{n}\right),$		(9)
	$\displaystyle g\left(i,0\right)$	$\displaystyle=$	$\displaystyle\sum_{s_{0}=0}^{t-1}\,d\left(i,i+s_{0}\bmod{n}\right).$		(10)

Lemma 5.

For all $i\in\{0,1,\ldots,n-1\}$ ,

\displaystyle f\left(i,h-1\right)=\sum_{j=0}^{n-1}\,\tilde{d}\left(i,i+j\bmod{n}\right).

Proof.

As $\sum_{r=0}^{h-1}\,s^{\prime}_{r}\cdot t^{r}=n-1$ ,

	$\displaystyle f\left(i,h-1\right)$	$\displaystyle=$	$\displaystyle\sum_{s_{h-1},s_{h-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\,\sum_{r=0}^{h-1}\,s_{r}\cdot t^{r}\leq n-1\,\right]$		(11)
		$\displaystyle\cdot$	$\displaystyle\sum_{k=0}^{h-1}\,d\left(i+\sum_{\ell=h-k}^{h-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right).\,\,\,\,\,$		(11)

By the existence and uniqueness of a $t$ -ary representation of each $j\in\{0,1,\ldots,n-1\}$ ,

	$\displaystyle\sum_{j=0}^{n-1}\,\sum_{k=0}^{h-1}\,d\left(i+\sum_{\ell=h-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n},i+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}(j)\cdot t^{\ell}\bmod{n}\right)$	(12)
$\displaystyle=$	$\displaystyle\sum_{s_{h-1},s_{h-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\,\sum_{r=0}^{h-1}\,s_{r}\cdot t^{r}\leq n-1\,\right]$
$\displaystyle\cdot$	$\displaystyle\sum_{k=0}^{h-1}\,d\left(i+\sum_{\ell=h-k}^{h-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+\sum_{\ell=h-1-k}^{h-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right).$

Equations (2) and (11)–(12) complete the proof. ∎

Lemma 6.

For all $i\in\{0,1,\ldots,n-1\}$ and $m\in\{1,2,\ldots,h-1\}$ ,

	$\displaystyle g\left(i,m\right)$	$\displaystyle=$	$\displaystyle t^{m}\sum_{s_{m}=0}^{t-1}\,d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)$
		$\displaystyle+$	$\displaystyle\sum_{s_{m}=0}^{t-1}\,g\left(i+s_{m}\cdot t^{m}\bmod{n},m-1\right).$

Proof.

By equation (8),

			$\displaystyle g\left(i,m\right)$
		$\displaystyle=$	$\displaystyle\sum_{s_{m}=0}^{t-1}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\left(d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)\vphantom{\sum_{\ell=m-1-k}^{m-2}\,s_{\ell}\cdot t^{\ell}}\right.$
			$\displaystyle\left.+\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right)\right),$
			$\displaystyle g\left(i+s_{m}\cdot t^{m}\bmod{n},m-1\right)$
		$\displaystyle=$	$\displaystyle\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right)$

for $s_{m}\in\{0,1,\ldots,t-1\}$ . Furthermore,

\displaystyle\sum_{s_{m}=0}^{t-1}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)=t^{m}\sum_{s_{m}=0}^{t-1}\,d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right).

∎

Lemma 7.

For all $i\in\{0,1,\ldots,n-1\}$ and $m\in\{1,2,\ldots,h-1\}$ ,

$\displaystyle f\left(i,m\right)$	$\displaystyle=$	$\displaystyle\left(1+\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\right)d\left(i,i+s^{\prime}_{m}\cdot t^{m}\bmod{n}\right)$
	$\displaystyle+$	$\displaystyle t^{m}\sum_{s_{m}=0}^{s^{\prime}_{m}-1}\,d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)$
	$\displaystyle+$	$\displaystyle f\left(i+s^{\prime}_{m}\cdot t^{m}\bmod{n},m-1\right)$
	$\displaystyle+$	$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}-1}\,g\left(i+s_{m}\cdot t^{m}\bmod{n},m-1\right).$

Proof.

Observe the following for all $s_{m}$ , $s_{m-1}$ , $\ldots$ , $s_{0}\in\{0,1,\ldots,t-1\}$ :

(i)
If $s_{m}=s^{\prime}_{m}$ , then $\sum_{r=0}^{m}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}$ if and only if $\sum_{r=0}^{m-1}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}$ ;
(ii)
If $s_{m}<s^{\prime}_{m}$ , then $\sum_{r=0}^{m}\,s_{r}\cdot t^{r}<\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}$ ;
(iii)
If $s_{m}>s^{\prime}_{m}$ , then $\sum_{r=0}^{m}\,s_{r}\cdot t^{r}>\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}$ .

We have

			$\displaystyle f\left(i,m\right)$
		$\displaystyle\stackrel{{\scriptstyle\text{equation~{}(\ref{subsumlessthanorequalto})}}}{{=}}$	$\displaystyle\sum_{s_{m}=0}^{t-1}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\sum_{r=0}^{m}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}\right]\cdot\left(d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)\vphantom{\sum_{\ell=m-2-k}^{m-2}\,s_{\ell}\cdot t^{\ell}\bmod{n}}\right.$
			$\displaystyle\left.+\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right)\right)$
		$\displaystyle\stackrel{{\scriptstyle\text{item~{}(\ref{item3})}}}{{=}}$	$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\left(\chi\left[\left(s_{m}=s^{\prime}_{m}\right)\land\left(\sum_{r=0}^{m}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}\right)\right]\right.$
			$\displaystyle\left.+\chi\left[\left(s_{m}<s^{\prime}_{m}\right)\land\left(\sum_{r=0}^{m}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}\right)\right]\right)\cdot\left(d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)\vphantom{\sum_{\ell=m-2-k}^{m-2}\,s_{\ell}\cdot t^{\ell}\bmod{n}}\right.$
			$\displaystyle\left.+\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right)\right)$
		$\displaystyle\stackrel{{\scriptstyle\text{item~{}(\ref{item1})}}}{{=}}$	$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\left(\chi\left[\left(s_{m}=s^{\prime}_{m}\right)\land\left(\sum_{r=0}^{m-1}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\right)\right]\right.$
			$\displaystyle\left.+\chi\left[\left(s_{m}<s^{\prime}_{m}\right)\land\left(\sum_{r=0}^{m}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m}\,s^{\prime}_{r}\cdot t^{r}\right)\right]\right)\cdot\left(d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)\vphantom{\sum_{\ell=m-2-k}^{m-2}\,s_{\ell}\cdot t^{\ell}\bmod{n}}\right.$
			$\displaystyle\left.+\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right)\right)$
		$\displaystyle\stackrel{{\scriptstyle\text{item~{}(\ref{item2})}}}{{=}}$	$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\left(\chi\left[\left(s_{m}=s^{\prime}_{m}\right)\land\left(\sum_{r=0}^{m-1}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\right)\right]\right.$
			$\displaystyle\left.+\chi\left[s_{m}<s^{\prime}_{m}\right]\vphantom{\sum_{r=0}^{m-2}\,s^{\prime}_{r}\cdot t^{r}}\right)\cdot\left(d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)\vphantom{\sum_{\ell=m-2-k}^{m-2}\,s_{\ell}\cdot t^{\ell}\bmod{n}}\right.$
			$\displaystyle\left.+\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right)\right).$

By equation (7),

			$\displaystyle f\left(i+s^{\prime}_{m}\cdot t^{m}\bmod{n},m-1\right)$
		$\displaystyle=$	$\displaystyle\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\,\sum_{r=0}^{m-1}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\,\right]$
		$\displaystyle\cdot$	$\displaystyle\sum_{k=0}^{m-1}\,d\left(i+s^{\prime}_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s^{\prime}_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right)$
		$\displaystyle=$	$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\,\left(s_{m}=s^{\prime}_{m}\right)\land\left(\sum_{r=0}^{m-1}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\right)\,\right]$
		$\displaystyle\cdot$	$\displaystyle\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right).$

By equation (8),

			$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}-1}\,g\left(i+s_{m}\cdot t^{m}\bmod{n},m-1\right)$
		$\displaystyle=$	$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\,s_{m}<s^{\prime}_{m}\,\right]$
		$\displaystyle\cdot$	$\displaystyle\sum_{k=0}^{m-1}\,d\left(i+s_{m}\cdot t^{m}+\sum_{\ell=m-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n},i+s_{m}\cdot t^{m}+\sum_{\ell=m-1-k}^{m-1}\,s_{\ell}\cdot t^{\ell}\bmod{n}\right).$

Because each number in $\{0,1,\ldots,\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\}$ can be written uniquely as $\sum_{r=0}^{m-1}\,s_{r}\cdot t^{r}$ , where $s_{m-1}$ , $s_{m-2}$ , $\ldots$ , $s_{0}\in\{0,1,\ldots,t-1\}$ ,

			$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\,\left(s_{m}=s^{\prime}_{m}\right)\land\left(\sum_{r=0}^{m-1}\,s_{r}\cdot t^{r}\leq\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\right)\,\right]$
		$\displaystyle\cdot$	$\displaystyle d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)$
		$\displaystyle=$	$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\left(1+\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\right)\cdot\chi\left[\,s_{m}=s^{\prime}_{m}\,\right]\cdot d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)$
		$\displaystyle=$	$\displaystyle\left(1+\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r}\right)d\left(i,i+s^{\prime}_{m}\cdot t^{m}\bmod{n}\right).$

Finally,

			$\displaystyle\sum_{s_{m}=0}^{s^{\prime}_{m}}\,\sum_{s_{m-1},s_{m-2},\ldots,s_{0}=0}^{t-1}\,\chi\left[\,s_{m}<s^{\prime}_{m}\,\right]\cdot d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right)$		(14)
		$\displaystyle=$	$\displaystyle t^{m}\sum_{s_{m}=0}^{s^{\prime}_{m}-1}\,d\left(i,i+s_{m}\cdot t^{m}\bmod{n}\right).$		(14)

Equations (2)–(14) complete the proof. ∎

t\leftarrow\lceil n^{1/h}\rceil

;

2: Find the

t

-ary representation of

n-1

, denoted

(s^{\prime}_{h-1},s^{\prime}_{h-2},\ldots,s^{\prime}_{0})\in\{0,1,\ldots,t-1\}^{h}

;

3: for

i\in\{0,1,\ldots,n-1\}

f[i][0]\leftarrow\sum_{s_{0}=0}^{s^{\prime}_{0}}\,d(i,i+s_{0}\bmod{n})

;

g[i][0]\leftarrow\sum_{s_{0}=0}^{t-1}\,d(i,i+s_{0}\bmod{n})

;

6: end for

7: for

m=1

up to

h-1

8: for

i\in\{0,1,\ldots,n-1\}

f[i][m]\leftarrow(1+\sum_{r=0}^{m-1}\,s^{\prime}_{r}\cdot t^{r})\,d(i,i+s^{\prime}_{m}\cdot t^{m}\bmod{n})

;

10:

f[i][m]\leftarrow f[i][m]+t^{m}\sum_{s_{m}=0}^{s^{\prime}_{m}-1}\,d(i,i+s_{m}\cdot t^{m}\bmod{n})

;

11:

f[i][m]\leftarrow f[i][m]+f[i+s^{\prime}_{m}\cdot t^{m}\bmod{n}][m-1]

;

12:

f[i][m]\leftarrow f[i][m]+\sum_{s_{m}=0}^{s^{\prime}_{m}-1}\,g[i+s_{m}\cdot t^{m}\bmod{n}][m-1]

;

13:

g[i][m]\leftarrow t^{m}\sum_{s_{m}=0}^{t-1}\,d(i,i+s_{m}\cdot t^{m}\bmod{n})

;

14:

g[i][m]\leftarrow g[i][m]+\sum_{s_{m}=0}^{t-1}\,g[i+s_{m}\cdot t^{m}\bmod{n}][m-1]

;

15: end for

16: end for

17: Output

\mathop{\mathrm{argmin}}_{i=0}^{n-1}\,f[i][h-1]

, breaking ties arbitrarily;

Figure 1: Algorithm find-median with input a metric space

(\{0,1,\ldots,n-1\},d)

and

h\in\mathbb{Z}^{+}\setminus\{1\}

Lemma 8.

Algorithm find-median in Fig. 1 is $(2h)$ -approximate for metric $1$ -median.

Proof.

By equations (9)–(10), lines 4–5 of find-median compute $f(i,0)$ and $g(i,0)$ . Then, by Lemmas 6–7, $f(i,m)$ and $g(i,m)$ can be found by dynamic programming as in lines 9–14. So line 17 outputs $\mathop{\mathrm{argmin}}_{i=0}^{n-1}\,f(i,h-1)$ , which equals

\mathop{\mathrm{argmin}}_{i=0}^{n-1}\,\sum_{j=0}^{n-1}\,\tilde{d}\left(i,i+j\bmod{n}\right)

by Lemma 5. Now Lemma 4 gives the approximation ratio of $2h$ . ∎

We now state our main theorem.

Theorem 9.

Metric $1$ -median has a deterministic $O(hn^{1+1/h})$ -time $(2h)$ -approximation nonadaptive algorithm for each $h\in\mathbb{Z}^{+}\setminus\{1\}$ .

Proof.

Clearly, find-median is deterministic and nonadaptive. Furthermore, it is $(2h)$ -approximate for metric $1$ -median by Lemma 8. As $s^{\prime}_{i}\leq t-1$ for all $i\in\{0,1,\ldots,h-1\}$ , the loop in lines 3–6 of find-median takes $O(nt)$ time. By precomputing $t^{i}$ and $\sum_{r=0}^{i}\,s^{\prime}_{r}\cdot t^{r}$ for all $i\in\{0,1,\ldots,h-1\}$ , each iteration of the loop in lines 8–15 takes $O(t)$ time. ∎

References

[1] M. R. Ackermann, J. Blömer, and C. Sohler. Clustering for metric and nonmetric distance measures. ACM Transactions on Algorithms, 6(4), 2010. Article 59.
[2] C.-L. Chang. Deterministic sublinear-time approximations for metric $1$ -median selection. Information Processing Letters, 113(8):288–292, 2013.
[3] C.-L. Chang. A lower bound for metric $1$ -median selection. Technical Report arXiv: 1401.2195, 2014.
[4] S. Guha, A. Meyerson, N. Mishra, R. Motwani, and L. O’Callaghan. Clustering data streams: Theory and practice. IEEE Transactions on Knowledge and Data Engineering, 15(3):515–528, 2003.
[5] P. Indyk. Sublinear time algorithms for metric space problems. In Proceedings of the 31st Annual ACM Symposium on Theory of Computing, pages 428–434, 1999.
[6] P. Indyk. High-dimensional computational geometry. PhD thesis, Stanford University, 2000.
[7] A. Kumar, Y. Sabharwal, and S. Sen. Linear-time approximation schemes for clustering problems in any dimensions. Journal of the ACM, 57(2):5, 2010.
[8] B.-Y. Wu. On approximating metric $1$ -median in sublinear time. Information Processing Letters, 114(4):163–166, 2014.

A deterministic sublinear-time nonadaptive algorithm for metric 111-median selection

Abstract

1 Introduction

Theorem 1 ([2]).

Theorem 2 ([8]).

2 Our algorithm

Lemma 3.

Proof.

Lemma 4.

Proof.

Lemma 5.

Proof.

Lemma 6.

Proof.

Lemma 7.

Proof.

Lemma 8.

Proof.

Theorem 9.

Proof.

References

A deterministic sublinear-time nonadaptive algorithm for metric $1$ -median selection