English

Non asymptotic distributional bounds for the Dickman Approximation of the running time of the Quickselect algorithm

Probability 2018-10-22 v3 Numerical Analysis

Abstract

Given a non-negative random variable WW and θ>0\theta>0, let the generalized Dickman transformation map the distribution of WW to that of W=dU1/θ(W+1), W^*=_d U^{1/\theta}(W+1), where UU[0,1]U \sim {\cal U}[0,1], a uniformly distributed variable on the unit interval, independent of WW, and where =d=_d denotes equality in distribution. It is well known that WW^* and WW are equal in distribution if and only if WW has the generalized Dickman distribution Dθ{\cal D}_\theta. We demonstrate that the Wasserstein distance d1d_1 between WW, a non-negative random variable with finite mean, and DθD_\theta having distribution Dθ{\cal D}_\theta obeys the inequality d1(W,Dθ)(1+θ)d1(W,W). d_1(W,D_\theta) \le (1+\theta)d_1(W,W^*). The specialization of this bound to the case θ=1\theta=1 and coupling constructions yield d_1(W_{n,1},D) \le \frac{8\log (n/2)+10}{n} \quad \mbox{for all $n \ge 1$, where} \quad W_{n,1}=\frac{1}{n}C_{n,1}-1, and Cn,mC_{n,m} is the number of comparisons made by the Quickselect algorithm to find the mthm^{th} smallest element of a list of nn distinct numbers. A similar bound holds for m2m \ge 2, and together recover the results of [12] that show distributional convergence of WnW_n to the standard Dickman distribution in the asymptotic regime m=o(n)m=o(n). By developing an exact expression for the expected running time E[Cn,m]E[C_{n,m}], lower bounds are provided that show the rate is not improvable for all m2m \not = 2.

Keywords

Cite

@article{arxiv.1703.00505,
  title  = {Non asymptotic distributional bounds for the Dickman Approximation of the running time of the Quickselect algorithm},
  author = {Larry Goldstein},
  journal= {arXiv preprint arXiv:1703.00505},
  year   = {2018}
}

Comments

Proof of Lemma 2.5 simplified, minor correction to proof of Theorem 1.3