English

Beating binary powering for polynomial matrices

Symbolic Computation 2023-05-29 v2

Abstract

The NNth 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~NN. When Fast Fourier Transform (FFT) is available, the resulting complexity is \emph{softly linear} in~NN, i.e.~linear in~NN with extra logarithmic factors. We show that it is possible to beat binary powering, by an algorithm whose complexity is \emph{purely linear} in~NN, even in absence of FFT. The key result making this improvement possible is that the entries of the NNth power of a polynomial matrix satisfy linear differential equations with polynomial coefficients whose orders and degrees are independent of~NN. Similar algorithms are proposed for two related problems: computing the NNth term of a C-finite sequence of polynomials, and modular exponentiation to the power NN for bivariate polynomials.

Keywords

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