English

Lower Bounds for Adaptive Sparse Recovery

Data Structures and Algorithms 2012-10-23 v2

Abstract

We give lower bounds for the problem of stable sparse recovery from /adaptive/ linear measurements. In this problem, one would like to estimate a vector xRnx \in \R^n from mm linear measurements A1x,...,AmxA_1x,..., A_mx. One may choose each vector AiA_i based on A1x,...,Ai1xA_1x,..., A_{i-1}x, and must output xx* satisfying |x* - x|_p \leq (1 + \epsilon) \min_{k\text{-sparse} x'} |x - x'|_p with probability at least 1δ>2/31-\delta>2/3, for some p{1,2}p \in \{1,2\}. For p=2p=2, it was recently shown that this is possible with m=O(1ϵkloglog(n/k))m = O(\frac{1}{\epsilon}k \log \log (n/k)), while nonadaptively it requires Θ(1ϵklog(n/k))\Theta(\frac{1}{\epsilon}k \log (n/k)). It is also known that even adaptively, it takes m=Ω(k/ϵ)m = \Omega(k/\epsilon) for p=2p = 2. For p=1p = 1, there is a non-adaptive upper bound of O~(1ϵklogn)\tilde{O}(\frac{1}{\sqrt{\epsilon}} k\log n). We show: * For p=2p=2, m=Ω(loglogn)m = \Omega(\log \log n). This is tight for k=O(1)k = O(1) and constant ϵ\epsilon, and shows that the loglogn\log \log n dependence is correct. * If the measurement vectors are chosen in RR "rounds", then m=Ω(Rlog1/Rn)m = \Omega(R \log^{1/R} n). For constant ϵ\epsilon, this matches the previously known upper bound up to an O(1) factor in RR. * For p=1p=1, m=Ω(k/(ϵlogk/ϵ))m = \Omega(k/(\sqrt{\epsilon} \cdot \log k/\epsilon)). This shows that adaptivity cannot improve more than logarithmic factors, providing the analog of the m=Ω(k/ϵ)m = \Omega(k/\epsilon) bound for p=2p = 2.

Keywords

Cite

@article{arxiv.1205.3518,
  title  = {Lower Bounds for Adaptive Sparse Recovery},
  author = {Eric Price and David P. Woodruff},
  journal= {arXiv preprint arXiv:1205.3518},
  year   = {2012}
}

Comments

19 pages; appearing at SODA 2013

R2 v1 2026-06-21T21:04:43.279Z