Algebraic diagonals and walks: Algorithms, bounds, complexity. (November 2017)

• Record Type: Journal Article
• Authors: Bostan, Alin; Dumont, Louis; Salvy, Bruno
• Is Part Of: Journal of symbolic computation. Volume 83(2017)
• Page Start: 68
• Page End: 92
• Publication Date: 2017-11
• Subjects: 33F10 -- 68W30 -- 05A15; Diagonals -- Walks -- Algorithms; Mathematics -- Data processing -- Periodicals; Numerical analysis -- Data processing -- Periodicals; Automatic programming (Computer science) -- Periodicals; Mathématiques -- Informatique -- Périodiques; Analyse numérique -- Informatique -- Périodiques; Programmation automatique -- Périodiques; Automatic programming (Computer science); Mathematics -- Data processing; Numerical analysis -- Data processing; Periodicals; Electronic journals; 510.285
• DOI: 10.1016/j.jsc.2016.11.006
• Languages: English
• ISSNs: 0747-7171

Abstract

Abstract: The diagonal of a multivariate power series F is the univariate power series Diag F generated by the diagonal terms of F . Diagonals form an important class of power series; they occur frequently in number theory, theoretical physics and enumerative combinatorics. We study algorithmic questions related to diagonals in the case where F is the Taylor expansion of a bivariate rational function. It is classical that in this case Diag F is an algebraic function. We propose an algorithm that computes an annihilating polynomial for Diag F . We give a precise bound on the size of this polynomial and show that, generically, this polynomial is the minimal polynomial and that its size reaches the bound. The algorithm runs in time quasi-linear in this bound, which grows exponentially with the degree of the input rational function. We then address the related problem of enumerating directed lattice walks. The insight given by our study leads to a new method for expanding the generating power series of bridges, excursions and meanders. We show that their first N terms can be computed in quasi-linear complexity in N, without first computing a very large polynomial equation.