Square-root sum problem

From HandWiki
Question, Web Fundamentals.svg Unsolved problem in computer science:
What is the Turing run-time complexity of the square-root sum problem?
(more unsolved problems in computer science)

The square-root sum problem (SRS) is a computational decision problem from the field of numerical analysis, with applications to computational geometry.

Definitions

SRS is defined as follows:[1]

Given positive integers [math]\displaystyle{ a_1,\ldots,a_k }[/math] and an integer t, decide whether [math]\displaystyle{ \sum_{i=1}^k \sqrt{a_i} \leq t }[/math].

An alternative definition is:

Given positive integers [math]\displaystyle{ a_1,\ldots,a_k }[/math] and [math]\displaystyle{ b_1,\ldots,b_k }[/math], decide whether [math]\displaystyle{ \sum_{i=1}^k \sqrt{a_i} \leq \sum_{i=1}^k \sqrt{b_i} }[/math].

Run-time complexity

SRS can be solved in polynomial time in the Real RAM model.[2] However, its run-time complexity in the Turing machine model is open, as of 1997.[1] The main difficulty is that, in order to solve the problem, the square-roots should be computed to a high accuracy, which may require a large number of bits. The problem is mentioned in the Open Problems Garden.[3]

Blomer[4] presents a polynomial-time Monte Carlo algorithm for deciding whether a sum of square roots equals 0.

Allender, Burgisser, Pedersen and Miltersen[5] prove that SRS lies in the counting hierarchy (which is contained in PSPACE).

Separation bounds

One way to solve SRS is to prove a lower bound on the absolute difference [math]\displaystyle{ \left|t - \sum_{i=1}^k \sqrt{a_i} \right| }[/math] or [math]\displaystyle{ \left|\sum_{i=1}^k \sqrt{a_i} - \sum_{i=1}^k \sqrt{b_i}\right| }[/math]. Such lower bound is called a "separation bound" since it separates between the difference and 0. For example, if the absolute difference is at least 2-d, it means that we can round all numbers to d bits of accuracy, and solve SRS in time polynomial in d.

This leads to the mathematical problem of proving bounds on this difference. Define r(n,k) as the smallest positive value of the difference [math]\displaystyle{ \sum_{i=1}^k \sqrt{a_i} - \sum_{i=1}^k \sqrt{b_i} }[/math], where ai and bi are integers between 1 and n; define R(n,k) is defined as -log r(n,k), which is the number of accuracy digits required to solve SRS. Computing r(n,k) is open problem 33 in the open problem project.[6]

In particular, it is interesting whether r(n,k) is in O(poly(k,log(n)). A positive answer would imply that SRS can be solved in polynomial time in the Turing Machine model. Some currently known bounds are:

  • Qian and Wang[7] prove by an explicit construction that, for any k and n, [math]\displaystyle{ r(n,k)\in O(n^{-2k+3/2}) }[/math], so [math]\displaystyle{ R(n,k)\geq (2k-3/2)\cdot \log{n} }[/math]. This number is optimal for k=2, and also for a wide range of integers.
  • Burnikel, Fleischer, Mehlhorn and Schirra[8] proved an upper bound on the number of digits: [math]\displaystyle{ R(n,k)\in O(2^{2 k}\cdot \log{n}) }[/math].
  • Cheng, Meng, Sun and Chen[9] showed that [math]\displaystyle{ R(n,k)\in 2^{O(n/\log{n})}\cdot \log{n} }[/math].
  • Cheng and Li[10] showed that [math]\displaystyle{ R(n,k)\in 2^{O(n/\log{n})} }[/math]. This implies an that SRS can be solved in time [math]\displaystyle{ 2^{o(k)} \cdot (\log{n})^{O(1)} }[/math], as long as n is in o(k log k). They also present an algorithm to compute r(n,k) in time [math]\displaystyle{ n^{k+o(k)} }[/math].
  • Eisenbrand, Haeberle and Singer[11] prove that [math]\displaystyle{ r(n,k)\geq \gamma\cdot n^{-2n} }[/math], where gamma is a constant that depends on the inputs a1,...,an, and steps from the Subspace theorem. This improves the previous bound [math]\displaystyle{ r(n,k)\geq \left( n\cdot \max_i (\sqrt{a_i})\right)^{-2^n} }[/math].

Applications

SRS is important in computational geometry, as Euclidean distances are given by square-roots, and many geometric problems (e.g. Minimum spanning tree in the plane and Euclidean traveling salesman problem) require to compute sums of distances.

Etessami and Yannakakis[12] show a reduction from SRS to the problem of termination of recursive concurrent stochastic games.

Relation to semidefinite programming

SRS also has a theoretic importance, as it is a simple special case of a semidefinite programming feasibility problem. Consider the matrix [math]\displaystyle{ \left( \begin{matrix} 1 & x \\ x & a \end{matrix} \right) }[/math]. This matrix is positive semidefinite iff [math]\displaystyle{ a - x^2 \geq 0 }[/math], iff [math]\displaystyle{ |x|\leq \sqrt{a} }[/math]. Therefore, to solve SRS, we can construct a feasibility problem with n constraints of the form [math]\displaystyle{ \left( \begin{matrix} 1 & x_i \\ x_i & a_i \end{matrix} \right) \succeq 0 }[/math], and additional linear constraints [math]\displaystyle{ x_i\geq 0, \sum_{i=1}^n x_i \geq k }[/math]. The resulting SDP is feasible if and only if SRS is feasible. As the runtime complexity of SRS in the Turing machine model is open, the same is true for SDP feasibility (as of 1997).

Extensions

Kayal and Saha[13] extend the problem from integers to polynomials. Their results imply a solution to SRS for a special class of integers.

References

  1. 1.0 1.1 Goemans, Michel X. (1997-10-01). "Semidefinite programming in combinatorial optimization" (in en). Mathematical Programming 79 (1): 143–161. doi:10.1007/BF02614315. ISSN 1436-4646. https://doi.org/10.1007/BF02614315. 
  2. Tiwari, Prasoon (1992-12-01). "A problem that is easier to solve on the unit-cost algebraic RAM". Journal of Complexity 8 (4): 393–397. doi:10.1016/0885-064X(92)90003-T. ISSN 0885-064X. https://dx.doi.org/10.1016/0885-064X%2892%2990003-T. 
  3. "Complexity of square-root sum | Open Problem Garden". http://garden.irmacs.sfu.ca/?q=op/complexity_of_square_root_sum. 
  4. "CSDL | IEEE Computer Society". https://www.computer.org/csdl/proceedings-article/focs/1991/0185434/12OmNAKuoSa. 
  5. Allender, Eric; Bürgisser, Peter; Kjeldgaard-Pedersen, Johan; Miltersen, Peter Bro (January 2009). "On the Complexity of Numerical Analysis" (in en). SIAM Journal on Computing 38 (5): 1987–2006. doi:10.1137/070697926. ISSN 0097-5397. http://epubs.siam.org/doi/10.1137/070697926. 
  6. Demaine, Erik D.. "TOPP: Problem 33: Sum of Square Roots". https://topp.openproblem.net/p33. 
  7. Qian, Jianbo; Wang, Cao An (2006-12-16). "How much precision is needed to compare two sums of square roots of integers?". Information Processing Letters 100 (5): 194–198. doi:10.1016/j.ipl.2006.05.002. ISSN 0020-0190. https://www.sciencedirect.com/science/article/pii/S0020019006001396. 
  8. Burnikel, C.; Fleischer, R.; Mehlhorn, K.; Schirra, S. (2000-05-01). "A Strong and Easily Computable Separation Bound for Arithmetic Expressions Involving Radicals" (in en). Algorithmica 27 (1): 87–99. doi:10.1007/s004530010005. ISSN 1432-0541. https://doi.org/10.1007/s004530010005. 
  9. Cheng, Qi; Meng, Xianmeng; Sun, Celi; Chen, Jiazhe (April 2010). "Bounding the sum of square roots via lattice reduction" (in en). Mathematics of Computation 79 (270): 1109–1122. doi:10.1090/S0025-5718-09-02304-7. ISSN 0025-5718. Bibcode2010MaCom..79.1109C. https://www.ams.org/mcom/2010-79-270/S0025-5718-09-02304-7/. 
  10. Cheng, Qi; Li, Yu-Hsin (2011-09-09). "On the minimum gap between sums of square roots of small integers". Theoretical Computer Science 412 (39): 5458–5465. doi:10.1016/j.tcs.2011.06.014. ISSN 0304-3975. 
  11. Eisenbrand, Friedrich; Haeberle, Matthieu; Singer, Neta (2023). "An improved bound on sums of square roots via the subspace theorem". arXiv:2312.02057 [cs.CG].
  12. Etessami, Kousha; Yannakakis, Mihalis (2008-11-11). "Recursive Concurrent Stochastic Games". Logical Methods in Computer Science 4 (4). doi:10.2168/LMCS-4(4:7)2008. ISSN 1860-5974. https://lmcs.episciences.org/1196. 
  13. Kayal, Neeraj; Saha, Chandan (2012-11-01). "On the Sum of Square Roots of Polynomials and Related Problems". ACM Transactions on Computation Theory 4 (4): 9:1–9:15. doi:10.1145/2382559.2382560. ISSN 1942-3454. https://doi.org/10.1145/2382559.2382560.