Algorithms for computing the sparsest shifts of polynomials via the Berlekamp/Massey algorithm

  • Authors:
  • Mark Giesbrecht;Erich Kaltofen;Wen-shin Lee

  • Affiliations:
  • University of Waterloo, Waterloo, Canada;North Carolina State University, Raleigh, North Carolina;University of Waterloo, Waterloo, Canada

  • Venue:
  • Proceedings of the 2002 international symposium on Symbolic and algebraic computation
  • Year:
  • 2002

Quantified Score

Hi-index 0.00

Visualization

Abstract

As a sub-procedure our algorithm executes the Berlekamp/Massey algorithm on a sequence of large integers or polynomials. We give a fraction-free version of the Berlekamp/Massey algorithm, which does not require rational numbers or functions and GCD operations on the arising numerators and denominators. The relationship between the solution of Toeplitz systems, Padé approximations, and the Euclidean algorithm is classical. Fraction-free versions [3] can be obtained from the subresultant PRS algorithm [2]. Dornstetter [6] gives an interpretation of the Berlekamp/Massey algorithm as a partial extended Euclidean algorithm. We map the subresultant PRS algorithm onto Dornstetter's formulation. We note that the Berlekamp/Massey algorithm is more efficient than the classical extended Euclidean algorithm.