Tractability of Approximation for Weighted Korobov Spaces on Classical and Quantum Computers

Abstract We study the approximation problem (or problem of optimal recovery in the $L_2$-norm) for weighted Korobov spaces with smoothness parameter $\a$. The weights $\gamma_j$ of the Korobov spaces moderate the behavior of periodic functions with respect to successive variables. The nonnegative smoothness parameter $\a$ measures the decay of Fourier coefficients. For $\a=0$, the Korobov space is the $L_2$ space, whereas for positive $\a$, the Korobov space is a space of periodic functions with some smoothness and the approximation problem corresponds to a compact operator. The periodic functions are defined on $[0,1]^d$ and our main interest is when the dimension $d$ varies and may be large. We consider algorithms using two different classes of information. The first class $\lall$ consists of arbitrary linear functionals. The second class $\lstd$ consists of only function values and this class is more realistic in practical computations. We want to know when the approximation problem is tractable. Tractability means that there exists an algorithm whose error is at most $\e$ and whose information cost is bounded by a polynomial in the dimension $d$ and in $\e^{-1}$. Strong tractability means that the bound does not depend on $d$ and is polynomial in $\e^{-1}$. In this paper we consider the worst case, randomized, and quantum settings. In each setting, the concepts of error and cost are defined differently and, therefore, tractability and strong tractability depend on the setting and on the class of information. In the worst case setting, we apply known results to prove that strong tractability and tractability in the class $\lall$ are equivalent. This holds if and only if $\a>0$ and the sum-exponent $s_{\g}$ of weights is finite, where $s_{\g}= \inf\{s>0 : \xxsum_{j=1}^\infty\g_j^s\,<\,\infty\}$. In the worst case setting for the class $\lstd$ we must assume that $\a>1$ to guarantee that functionals from $\lstd$ are continuous. The notions of strong tractability and tractability are not equivalent. In particular, strong tractability holds if and only if $\a>1$ and $\xxsum_{j=1}^\infty\g_j<\infty$. In the randomized setting, it is known that randomization does not help over the worst case setting in the class $\lall$. For the class $\lstd$, we prove that strong tractability and tractability are equivalent and this holds under the same assumption as for the class $\lall$ in the worst case setting, that is, if and only if $\a>0$ and $s_{\g} < \infty$. In the quantum setting, we consider only upper bounds for the class $\lstd$ with $\a>1$. We prove that $s_{\g}<\infty$ implies strong tractability. Hence for $s_{\g}>1$, the randomized and quantum settings both break worst case intractability of approximation for the class $\lstd$. We indicate cost bounds on algorithms with error at most $\e$. Let $\cc(d)$ denote the cost of computing $L(f)$ for $L\in \lall$ or $L\in \lstd$, and let the cost of one arithmetic operation be taken as unity. The information cost bound in the worst case setting for the class $\lall$ is of order $\cc (d) \cdot \e^{-p}$ with $p$ being roughly equal to $2\max(s_\g,\a^{-1})$. Then for the class $\lstd$ in the randomized setting, we present an algorithm with error at most $\e$ and whose total cost is of order $\cc(d)\e^{-p-2} + d\e^{-2p-2}$, which for small $\e$ is roughly $$ d\e^{-2p-2}. $$ In the quantum setting, we present a quantum algorithm with error at most $\e$ that uses about only $d + \log \e^{-1}$ qubits and whose total cost is of order $$ (\cc(d) +d) \e^{-1-3p/2}. $$ The ratio of the costs of the algorithms in the quantum setting and the randomized setting is of order $$ \frac{d}{\cc(d)+d}\,\left(\frac1{\e}\right)^{1+p/2}. $$ Hence, we have a polynomial speedup of order $\e^{-(1+p/2)}$. We stress that $p$ can be arbitrarily large, and in this case the speedup is huge.

[1]  P. Shor INTRODUCTION TO QUANTUM ALGORITHMS , 2000, quant-ph/0005003.

[2]  G. Wahba Spline models for observational data , 1990 .

[3]  Arthur O. Pittenger,et al.  An Introduction to Quantum Computing Algorithms , 2000 .

[4]  Erich Novak,et al.  The Real Number Model in Numerical Analysis , 1995, J. Complex..

[5]  Henryk Wozniakowski,et al.  Path Integration on a Quantum Computer , 2002, Quantum Inf. Process..

[6]  H. Wozniakowski Efficiency of Quasi-Monte Carlo Algorithms for High Dimensional Integrals , 2000 .

[7]  Artur Ekert,et al.  Basic Concepts in Quantum Computation (量子情報理論とその応用論文小特集) , 1998 .

[8]  Henryk Wozniakowski,et al.  When Are Quasi-Monte Carlo Algorithms Efficient for High Dimensional Integrals? , 1998, J. Complex..

[9]  Henryk Wozniakowski,et al.  Information-based complexity , 1987, Nature.

[10]  N. Aronszajn Theory of Reproducing Kernels. , 1950 .

[11]  Jozef Gruska,et al.  Quantum Computing , 2008, Wiley Encyclopedia of Computer Science and Engineering.

[12]  Erich Novak,et al.  Complexity of Linear Problems with a Fixed Output Basis , 2000, J. Complex..

[13]  Lov K. Grover A fast quantum mechanical algorithm for database search , 1996, STOC '96.

[14]  H. Woxniakowski Information-Based Complexity , 1988 .

[15]  Henryk Wozniakowski,et al.  On the Power of Standard Information for Weighted Approximation , 2001, Found. Comput. Math..

[16]  F. J. Hickernell,et al.  Tractability of Multivariate Integration for Periodic Functions , 2001, J. Complex..

[17]  Stefan Heinrich Quantum Summation with an Application to Integration , 2002, J. Complex..

[18]  Henryk Wozniakowski,et al.  Weighted Tensor Product Algorithms for Linear Multivariate Problems , 1999, J. Complex..

[19]  Henryk Wozniakowski,et al.  Tractability of Multivariate Integration for Weighted Korobov Classes , 2001, J. Complex..

[20]  Fred J. Hickernell,et al.  Integration and approximation in arbitrary dimensions , 2000, Adv. Comput. Math..

[21]  E. Novak,et al.  Foundations of Computational Mathematics: When are integration and discrepancy tractable? , 2001 .

[22]  Henryk Wozniakowski,et al.  Intractability Results for Integration and Discrepancy , 2001, J. Complex..

[23]  Erich Novak Quantum Complexity of Integration , 2001, J. Complex..

[24]  Henryk Wozniakowski,et al.  Tractability and Strong Tractability of Linear Multivariate Problems , 1994, J. Complex..

[25]  Stefan Heinrich Quantum integration in Sobolev classes , 2003, J. Complex..

[26]  Erich Novak,et al.  Optimal linear randomized methods for linear operators in Hilbert spaces , 1992, J. Complex..

[27]  Ronald de Wolf,et al.  Quantum lower bounds by polynomials , 2001, JACM.

[28]  Lov K. Grover A framework for fast quantum mechanical algorithms , 1997, STOC '98.

[29]  Erich Novak,et al.  On a problem in quantum summation , 2003, J. Complex..

[30]  S. Heinrich,et al.  Optimal Summation and Integration by Deterministic, Randomized, and Quantum Algorithms , 2002 .

[31]  G. Brassard,et al.  Quantum Amplitude Amplification and Estimation , 2000, quant-ph/0005055.

[32]  Grace Wahba,et al.  Spline Models for Observational Data , 1990 .

[33]  Felix Wu,et al.  The quantum query complexity of approximating the median and related statistics , 1998, STOC '99.