Core Concepts
We present a simple and efficient algorithm that computes the composition of two power series in near-linear time complexity, improving upon the previous best algorithms.
Abstract
The paper presents an efficient algorithm for computing the composition of two power series, f(g(x)) mod xn, where f(x) and g(x) are polynomials in A[x] of degrees less than m and n, respectively.
Key highlights:
The previous best algorithms for general power series composition had a time complexity of O(n^(1+o(1))) or O(n^1.43), whereas the new algorithm achieves a time complexity of O(M(n) log m + M(m)), where M(n) is the time complexity of polynomial multiplication.
The algorithm builds upon the Graeffe iteration approach introduced by Bostan and Mori to manipulate rational power series, reducing the problem to a series of polynomial multiplications with exponentially decreasing degrees.
The algorithm also has consequences for other computation models, such as boolean circuits and multitape Turing machines, depending on the fast arithmetic available in the underlying ring A.
The paper also presents an efficient algorithm for the transposed problem, power projection, and shows how the transposition principle can be used to derive the power series composition algorithm.
Stats
The previous best algorithms for general power series composition had a time complexity of O(n^(1+o(1))) or O(n^1.43).
The new algorithm achieves a time complexity of O(M(n) log m + M(m)), where M(n) is the time complexity of polynomial multiplication.
Quotes
"We present a simple algorithm that reduces the complexity of power series composition to near-linear time."
"Our algorithm also has consequences for other computation models measured by bit complexity, such as boolean circuits and multitape Turing machines."