The $N$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 $N$. When Fast Fourier Transform (FFT) is available, the resulting arithmetic complexity is \emph{softly linear} in $N$, i.e. linear in $N$ with extra logarithmic factors. We show that it is possible to beat binary powering, by an algorithm whose complexity is \emph{purely linear} in $N$, even in absence of FFT. The key result making this improvement possible is that the entries of the $N$th power of a polynomial matrix satisfy linear differential equations with polynomial coefficients whose orders and degrees are independent of $N$. Similar algorithms are proposed for two related problems: computing the $N$th term of a C-recursive sequence of polynomials, and modular exponentiation to the power $N$ for bivariate polynomials.
翻译:固定大小和度的多元矩阵的“美元”的功率可以通过以美元乘以两个线性多边矩阵的二元功率来计算。当快速傅里叶变形(FFT)可用时,由此产生的算术复杂性是以美元计算,即以美元线性计算,即以美元线性计算,加上额外的对数因素。我们表明,可以通过复杂程度为 empph{纯线性}的算法,以美元乘以美元来击败二元动力,即使没有FFFT。 使这一改进成为可能的关键性结果是,多元矩阵的“美元”的分数符合多数值的线性差异方程式,其顺序和水平独立于$N。为两个相关的问题提出了类似的算法:计算一个C-对数序列的“美元”术语,以及将“美元”的组合式推算出用于双变量的“美元”。