Beating binary powering for polynomial matrices
Abstract
The th power of a polynomial matrix of fixed size and degree can be computed by binary powering as fast as multiplying two polynomials of linear degree in~. When Fast Fourier Transform (FFT) is available, the resulting complexity is \emph{softly linear} in~, i.e.~linear in~ with extra logarithmic factors. We show that it is possible to beat binary powering, by an algorithm whose complexity is \emph{purely linear} in~, even in absence of FFT. The key result making this improvement possible is that the entries of the th power of a polynomial matrix satisfy linear differential equations with polynomial coefficients whose orders and degrees are independent of~. Similar algorithms are proposed for two related problems: computing the th term of a C-finite sequence of polynomials, and modular exponentiation to the power for bivariate polynomials.
Cite
@article{arxiv.2302.04299,
title = {Beating binary powering for polynomial matrices},
author = {Alin Bostan and Vincent Neiger and Sergey Yurkevich},
journal= {arXiv preprint arXiv:2302.04299},
year = {2023}
}
Comments
10 pages, 3 figures, 2 tables, 5 algorithms