Approximating Iterated Multiplication of Stochastic Matrices in Small Space

Matrix powering, and more generally iterated matrix multiplication, is a fundamental linear algebraic primitive with myriad applications in computer science. Of particular interest is the problem's space complexity as it constitutes the main route towards resolving the BPL vs. L problem. The seminal work by Saks and Zhou [32] gives a deterministic algorithm for approximating the product of 𝑛 stochastic matrices of dimension 𝑤 ×𝑤 in space 𝑂 (log 3/2 𝑛 + √︁ log 𝑛 •log 𝑤). The first improvement upon [32] was achieved by Hoza [15] who gave a logarithmic improvement in the 𝑛 = poly(𝑤) regime, attaining 𝑂 (

We give the first polynomial improvement over [32]. Our algorithm achieves space complexity of

In particular, in the regime log 𝑛 > log 2 𝑤, our algorithm runs in nearly-optimal 𝑂 (log 𝑛) space, improving upon the previous best 𝑂 (log 3/2 𝑛).

To obtain our result for the special case of matrix powering, we harness recent machinery from time-and space-bounded Laplacian solvers to the framework of [32] and devise an intricate precisionalternating recursive scheme. This enables us to bypass the bottleneck of paying log 𝑛-space per recursion level. The general case of iterated matrix multiplication poses several additional challenges, the substantial of which is handled by devising an improved shift and truncate mechanism. The new mechanism is made possible by a novel use of the Richardson iteration.