Algorithms for Linearly Recurrent Sequences of Truncated Polynomials

Abstract

Linear recurrent sequences are those whose elements are defined as linear combinations of preceding elements, and finding recurrence relations is a fundamental problem in computer algebra. In this paper, we focus on sequences whose elements are vectors over the ring A = K[ ]/⟨ ⟩ of truncated polynomials. Finding the ideal of their recurrence relations has applications such as the computation of minimal polynomials and determinants of sparse matrices over A. We present three methods for finding this ideal: a Berlekamp-Massey-like approach due to Kurakin, one which computes the kernel of some block-Hankel matrix over A via a minimal approximant basis, and one based on bivariate Padé approximation. We propose complexity improvements for the first two methods, respectively by avoiding the computation of redundant relations and by exploiting the Hankel structure to compress the approximation problem. Then we confirm these improvements empirically through a C++ implementation, and we discuss the above-mentioned applications. CCS CONCEPTS• Computing methodologies → Algebraic algorithms; • Theory of computation → Design and analysis of algorithms.

Rank-Sensitive Computation of the Rank Profile of a Polynomial Matrix

“…Furthermore, there are specific situations where rank deficiency is actually expected by design, and one would like to be able to take advantage of this in algorithms. Recently, in a context of computing generators of linearly recurrent sequences, minimal approximant bases of rank-deficient, structured matrices F have been encountered [13,Sec. 5].…”

Rank-Sensitive Computation of the Rank Profile of a Polynomial Matrix

“…Furthermore, there are specific situations where rank deficiency is actually expected by design, and one would like to take advantage of this in algorithms. Recently, in a context of computing generators of linearly recurrent sequences, minimal approximant bases of rankdeficient, structured matrices F have been encountered [13,Sec. 5].…”