Fast evaluation of multiple zeta sums (Q4396455)

We show that the multiple zeta sum \[ \zeta(s_{1}, s_{2}, \dots, s_{d}) = \sum_{n_{1} > n_{2} > \dots > n_{d}} \frac{1 }{n_{1}^{s_{1}} n_{2}^{s_{2}} \dots n_{d}^{s_d}}, \] for positive integers \(s_{i}\), with \(s_{1}>1\), can always be written as a finite sum of products of rapidly convergent series. Perhaps surprisingly, one may develop fast summation algorithms of such efficiency that the overall complexity can be brought down essentially to that of one-dimensional summation. In particular, for any dimension \(d\) one may resolve \(D\) good digits of \(\zeta\) in \[ O(D \log D / \log \log D) \] arithmetic operations, with the implied big-\(O\) constant depending only on the set \(\{s_{1},\dots,s_{d}\}\).