Computer Algebra for Lattice Path Combinatorics Alin Bostan

Total Page:16

File Type:pdf, Size:1020Kb

Computer Algebra for Lattice Path Combinatorics Alin Bostan Computer Algebra for Lattice path Combinatorics Alin Bostan To cite this version: Alin Bostan. Computer Algebra for Lattice path Combinatorics. Symbolic Computation [cs.SC]. Université Paris 13, 2017. tel-01660300 HAL Id: tel-01660300 https://hal.archives-ouvertes.fr/tel-01660300 Submitted on 12 Dec 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Université Paris 13 Laboratoire d’Informatique de Paris Nord Habilitation à Diriger des Recherches Spécialité : Sciences Calcul Formel pour la Combinatoire des Marches Soutenue le 15 décembre 2017 par Alin Bostan (Inria) devant le jury composé de : Mme. Frédérique Bassino Université Paris 13, Villetaneuse M. Olivier Bodini Université Paris 13, Villetaneuse Mme. Mireille Bousquet-Mélou CNRS, Université de Bordeaux Mme. Lucia Di Vizio CNRS, Université de Versailles M. Mark Giesbrecht Université de Waterloo, Canada (rapporteur) M. Florent Hivert Université Paris 11, Orsay M. Christian Krattenthaler Université de Vienne, Autriche (rapporteur) M. Gilles Villard CNRS, ENS de Lyon (rapporteur) COMPUTER ALGEBRA FOR LATTICE PATH COMBINATORICS ALIN BOSTAN∗ Abstract. Classifying lattice walks in restricted lattices is an important problem in enumerative combinatorics. Recently, computer algebra has been used to explore and to solve a number of diffi- cult questions related to lattice walks. We give an overview of recent results on structural properties and explicit formulas for generating functions of walks in the quarter plane, with an emphasis on the algorithmic methodology. Key words. Enumerative combinatorics, random walks in cones, lattice paths in the quarter plane, Gessel walks, generating functions, computer algebra, automated guessing, creative telescoping, diago- nals, binomial sums, algebraic functions, D-finite functions, hypergeometric functions, elliptic integrals. AMS subject classifications. Primary 05A10, 05A15, 05A16, 97N70, 33F10, 68W30, 14Q20; Secondary 33C05, 97N80, 13P15, 33C75, 12Y05, 13P05, 14Q20. This document is structured as follows. Section1 gives an overview of recent re- sults obtained in lattice path combinatorics with the help of computer algebra, with a focus on the exact enumeration of walks confined to the quarter plane. Sections2 and3 then go into more details of two classes of fruitful algorithmic approaches: guess-and-prove and creative telescoping. 1. General presentation. 1.1. Prelude. Consider the following innocent-looking problem. A tandem-walk is a path in Z2 taking steps from f", , &g only. Show that, for any integer n ≥ 0, the following quantities are equal: (i) the number an of tandem-walks of length n (i.e., using n steps), confined to the upper half-plane Z × N, that start and end at (0, 0); (ii) the number bn of tandem-walks of length n confined to the quar- ter plane N2, that start at (0, 0) and finish on the diagonal x = y. For instance, for n = 3, this common value is a3 = b3 = 3, as shown below. (i) (ii) The problem establishes a rather surprising connection between tandem-walks in the lattice plane, submitted to two different kinds of constraints: the evolution domain of the walk, and its ending point. The domain constraint is weaker for the first family of walks, while the ending constraint is relaxed for the second family. It appears that this problem is far from being trivial. Several solutions exist, but none of them is elementary. One of the main aims of the present text is to ∗Inria, Université Paris-Saclay, 91120 Palaiseau, France ([email protected]). 2 COMPUTER ALGEBRA FOR LATTICE PATH COMBINATORICS 3 convince the reader that this problem (and many others with a similar flavor) can be solved with the help of a computer. More precisely, Computer Algebra tools, extensively described in the following sections, can be used to discover and to prove the following equalities (3n)! (1) a = b = , and am = bm = 0 if 3 does not divide m. 3n 3n n!2 · (n + 1)! It goes without saying that such a simple and beautiful expression cannot be an element of chance. As it will turn out, closed forms are quite rare for this kind of enumeration problems. Nevertheless, even in absence of nice formulas, the struc- tural properties of the corresponding enumeration sequences reflect the symmetries of the step set and of the evolution domain. Equation (1) shows that the sequences (an) and (bn) are P-recursive, that is, they satisfy a linear recurrence with polyno- mial coefficients (in the index n). One of the messages that will emerge from the text is that this important property of the enumeration sequences is intimately related to the finiteness of a certain group, naturally attached to the step set f", , &g. 1.2. General context: lattice paths confined to cones. Let us put the previous problem into a more general framework. Let d ≥ 1 be an integer (dimension), let S d d be a finite subset (called step set, or model) of vectors in Z , and p0 2 Z (starting point). A S-path (or S-walk) of length n starting at p0 is a sequence (p0, p1,..., pn) d of elements in the lattice Z such that pi+1 − pi 2 S for all 0 ≤ i < n. Let C be a cone of Rd, that is a subset of Rd such that r · v 2 C for any v 2 C and r > 0, assumed to contain p0. We will be interested in the (exact and asymptotic) enumeration of S-walks confined to the cone C, and potentially subject to additional constraints. Example 1. Consider the model S = f(1, 0), (−1, 0), (1, −1), (−1, 1)g (called the Gouyou-Beauchamps model) in dimension d = 2, with starting point p0 = (0, 0) and 2 with cone C = R+ (the quarter plane). The picture below displays the step set of the model (on the left), and a S-walk of length n = 17 confined to C (on the right). (i, j) = (5, 1) The main typical questions in this context are then the following: • What is the number an of n-step S-walks contained in C and starting at p0? • For fixed i 2 C, what is the number an;i of such walks that end at i? • What is the nature of their generating functions n n i A(t) = ∑ ant and A(t; x) = ∑ an;it x ? n n,i As expected from the introductory example of tandem-walks, the answers to these questions are not simple, and heavily depend on the various parameters. The aim of this text is to provide a survey of recent results —notably classification results and closed form expressions— obtained using Computer Algebra. 4 ALIN BOSTAN 1.3. Why count walks in cones?. Lattice paths are fundamental objects in com- binatorics. They have been studied at least since the second half of the 19th century, in connection with the ballot problem (see §1.4). Even earlier, embryonic occurrences (around 1650) are in Pascal’s and Huygens’ solutions of the so-called problem of di- vision of the stakes (or, problem of points), and of the gambler’s ruin problem, which motivated the beginnings of modern probability theory [170, 226, 157]. Despite these historically important examples, the enumeration of lattice walks has long re- mained part of what may be called recreational mathematics. It is only in the late 1960s that their study really became an independent field of research, at the cross- roads of pure and applied mathematics. Since then, various approaches have been progressively involved, separately or in interaction, in the study of lattice walks. These methods arise from various fields of classical mathematics (algebra, combi- natorics, complex analysis, probability theory), and more recently from computer science. There are several reasons for the ubiquity of lattice walks, but the most solid one is that they encode several important classes of mathematical objects, in discrete mathematics (permutations, trees, words, urns, . ), in statistical physics (magnetism, polymers, . ), in probability theory (branching processes, games of chance, . ), in operations research (birth-death processes, queueing theory, . ). Therefore, many questions from all these various fields can be reduced to solving lattice path problems. For more motivations, the reader is referred to the introduc- tion of [26]. Nowadays, several books are entirely devoted to lattice paths and their applications [355, 312, 315, 160, 146, 384, 180, 388, 47, 284, 44], and an international conference titled Lattice path combinatorics and applications is entirely devoted to this field. We recommend Humphreys’ article [237] for a brief review of the history of lattice path enumeration and for a survey of the recent evolution of the field. Also, Krattenthaler’s recent survey [269] is an excellent overview of various results and methods in lattice path enumeration. 1.4. The ballot problem and the reflection principle. As mentioned before, the enumeration of lattice walks is an old topic. We want to illustrate this using Bertrand’s ballot problem [36, 10]. The aim is not only to provide the flavor of a nice piece of combinatorial reasoning, but especially to introduce the so-called reflection principle, seemingly invented by Aebly and Mirimanoff [5, 306], which contains the roots of a systematic method for lattice walks, to be presented later, and based on the notion of group of a walk, see §1.18. Bertrand’s problem is the following: Suppose that two candidates A and B are running in an election.
Recommended publications
  • Mathematisches Forschungsinstitut Oberwolfach Enumerative
    Mathematisches Forschungsinstitut Oberwolfach Report No. 12/2014 DOI: 10.4171/OWR/2014/12 Enumerative Combinatorics Organised by Mireille Bousquet-M´elou, Bordeaux Michael Drmota, Wien Christian Krattenthaler, Wien Marc Noy, Barcelona 2 March – 8 March 2014 Abstract. Enumerative Combinatorics focusses on the exact and asymp- totic counting of combinatorial objects. It is strongly connected to the proba- bilistic analysis of large combinatorial structures and has fruitful connections to several disciplines, including statistical physics, algebraic combinatorics, graph theory and computer science. This workshop brought together experts from all these various fields, including also computer algebra, with the goal of promoting cooperation and interaction among researchers with largely vary- ing backgrounds. Mathematics Subject Classification (2010): Primary: 05A; Secondary: 05E, 05C80, 05E, 60C05, 60J, 68R, 82B. Introduction by the Organisers The workshop Enumerative Combinatorics organized by Mireille Bousquet-M´elou (Bordeaux), Michael Drmota (Vienna), Christian Krattenthaler (Vienna), and Marc Noy (Barcelona) took place on March 2-8, 2014. There were over 50 par- ticipants from the US, Canada, Australia, Japan, Korea, and various European countries. The program consisted of 13 one hour lectures, accompanied by 17 shorter contributions and the special session of 5 presentations by Oberwolfach Leibniz graduate students. There was also a lively problem session led by Svante Linusson. The lectures were intended to provide overviews of the state of the art in various areas and to present relevant new results. The lectures and short talks ranged over a wide variety of topics including classical enumerative prob- lems, algebraic combinatorics, asymptotic and probabilistic methods, statistical 636 Oberwolfach Report 12/2014 physics, methods from computer algebra, among others.
    [Show full text]
  • Ilse Fischer Platz 1 Universität Wien 1090 Wien Ifischer/ T +43-1-4277-50454 Ilse.fischer @Univie.Ac.At U +43-1-4277-50450
    Fakultät für Mathematik Oskar-Morgenstern- Ilse Fischer Platz 1 Universität Wien 1090 Wien http://www.mat.univie.ac.at/˜ ifischer/ T +43-1-4277-50454 ilse.fischer @univie.ac.at u +43-1-4277-50450 Born on June 29, 1975, in Klagenfurt, Austria. Education 2006 Habilitation in Mathematics, Universität Wien. Thesis A polynomial method for the enumeration of plane partitions and alternating sign matrices 2000 Dr. rer. nat in Mathematics, Universität Wien, summa cum laude. Thesis Enumeration of perfect matchings: Rhombus tilings and Pfaffian graphs Supervisors Christian Krattenthaler and Franz Rendl 1998 Mag. rer. nat in Mathematics, Universität Wien, summa cum laude. Thesis Der Satz von Skolem, Mahler, Lech – ein Ergebnis über die Nullstellen linearer Rekursionsfolgen Supervisor Johann Cigler 1993 – 1998 Studies of Mathematics, Universität Wien. Academic Employment since 2017 Full Professor, Universität Wien. 2011 – 2017 Associate Professor, Universität Wien. 2007 and 2009 Visiting Professor, Alpen-Adria-Universität Klagenfurt. 2004 – 2011 Universitätsassistentin, Universität Wien, permanent position since 2008. 2002 Visiting Researcher, Universität Wien. 1999 – 2004 Universitätsassistentin, Alpen-Adria-Universität Klagenfurt. 1995 –1999 Teaching Assistant, Universität Wien. Awards David P. Robbins Prize, American Mathematical Society (AMS), jointly with Roger 2019 E. Behrend (Cardiff) and Matjaž Konvalinka (Ljubljana). 1/7 2009 START Award, Austrian Federal Ministry of Science and Research / Austrian Sci- ence Fund FWF. 2006 Dr.-Maria-Schaumayer-Preis. Research Grants 2017 – 2021 Principal Investigator in the SFB (Special Research Area) “Algorithmic and Enumerative Combinatorics”, Project part: “Alternating sign arrays of triangular shapes”, Austrian Science Fund (FWF), e 364.875. 2010 – 2018 Principal Investigator of the START project, “Compact enumeration formulas for generalized partitions”, Austrian Science Fund (FWF), e 1.1 Mill.
    [Show full text]
  • Curriculum Vitae
    Curriculum Vitae Personal Data Full name Priv.-Doz. Dipl.-Inf. Dr. Christoph Bernhard Koutschan Date of birth 12.12.1978 Place of birth Dillingen an der Donau, Germany Nationality German Marital status Married, two children Contact Address Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences Altenberger Straße 69, A-4040 Linz Phone +43 732 2468 5254 E-mail [email protected] Homepage http://www.koutschan.de ORCID 0000-0003-1135-3082 Education and Work Experience 05/2017 Habilitation in Mathematics, Johannes Kepler University Linz, Austria. Thesis: Quod Erat Demonstrandum: Proofs by Computer (evaluation committee: James H. Davenport, Christian Krattenthaler, Marko Petkovˇsek,Bruno Salvy, Carsten Schneider, Nobuki Takayama) 09/2012 { Research scientist at the Johann Radon Institute for Computational and Applied Mathematics (RICAM), Linz, Austria. 09/2011 { 08/2012 Postdoctoral researcher at Institut National de Recherche en Informatique et en Automatique (INRIA), MSR-INRIA Joint Centre, Orsay, France. 09/2010 { 08/2011 Postdoctoral researcher at the Research Institute for Symbolic Computa- tion (RISC), Linz, Austria. 10/2009 { 06/2010 Postdoctoral researcher at Tulane University, New Orleans, USA. 10/2005 { 09/2009 Ph.D. studies in symbolic computation at RISC, Johannes Kepler Univer- sity Linz, Austria. Ph.D. thesis: Advanced Applications of the Holonomic Systems Approach (advisor: Univ.-Prof. Dr. Peter Paule) Christoph Koutschan Curriculum Vitae 2 10/1999 { 07/2005 Undergraduate studies in computer science (minor subject: mathematics), Friedrich-Alexander University of Erlangen-N¨urnberg, Germany. Diploma thesis: Regular Languages and Their Generating Function: The Inverse Problem (advisor: Prof. Dr. Volker Strehl) 07/1998 { 07/1999 Civil Service 06/1998 High school graduation (Abitur) in Coburg, Germany.
    [Show full text]
  • Meetings and Conferences of The
    Meetings & Conferences of the AMS IMPORTANT INFORMATION REGARDING MEETINGS PROGRAMS: AMS Sectional Meeting programs do not appear in the print version of the Notices. However, comprehensive and continually updated meeting and program information with links to the abstract for each talk can be found on the AMS website.See http://www.ams.org/meetings/. Programs and abstracts will continue to be displayed on the AMS website in the Meetings and Conferences section until about three weeks after the meeting is over. Final programs for Sectional Meetings will be archived on the AMS website in an electronic issue of the Notices as noted below for each meeting. Christian Krattenthaler, University of Lyon-I, Exact and Mainz, Germany asymptotic enumeration of vicious walkers with a wall interaction. June 16–19, 2005 Frank Natterer, University of Muenster, Imaging and Thursday – Sunday inverse problems for partial differential equations. Meeting #1008 Horng-Tzer Yau, New York University and Stanford University, Dynamics of Bose-Einstein condensates. Joint International Meeting with the Deutsche Mathematiker- Vereinigung (DMV) and the Oesterreichische Mathematis- Special Sessions che Gesellschaft (OMG) Associate secretary: Susan J. Friedlander Affine Algebraic Geometry, Shreeram Abhyankar, Purdue Announcement issue of Notices: February 2005 University, Hubert Flenner, Ruhr University Bochum, and Makar Limanov, Wayne State University. Program first available on AMS website: Not applicable Program issue of electronic Notices: Not applicable Algebraic Combinatorics, Patricia Hersh, Indiana Univer- Issue of Abstracts: Not applicable sity-Bloomington, Christian Krattenthaler, University of Lyon-I, and Volkmar Welker, Philipps University Marburg. Deadlines Algebraic Cryptography, Dorian Goldfeld, Columbia For organizers: Expired University, Martin Kreuzer and Gerhard Rosenberger, For consideration of contributed papers in Special Sessions: Universität Dortmund, and Vladimir Shpilrain, The City Expired College of New York.
    [Show full text]
  • Notices of the American Mathematical Society
    ISSN 0002-9920 (print) ISSN 1088-9477 (online) of the American Mathematical Society September 2018 Volume 65, Number 8 Hispanic Heritage Month September 15–October 15, 2018 The Geometry of Matroids page 902 Movie Animation: A Continuum Approach for Frictional Contact page 909 September Sectional Sampler: Braids, Surfaces, and Homological Invariants page 914 2018 AMS Election page 917 AMS/MAA Grad School Fair . Meeting and Recruiting Students for Your Graduate School Programs! The 2019 Grad School Fair: • Held at the Baltimore Convention Center, Baltimore, MD, during the Joint Mathematics Meetings in January • Research-oriented special undergraduate programs • Connect with over 300 talented student attendees. • More than 60 graduate programs in the mathematical sciences represented! The event is free for registered students to attend. Schools will pay a small table fee to represent their programs. Learn more at: For further information: www.ams.org/gradfair phone: 800-321-4AMS, ext. 4060 email: [email protected] Please check the Joint Mathematics Meetings registration site for updated dates and times. Notices of the American Mathematical Society September 2018 FEATURED 902684 909 26914 The Geometry of Matroids Movie Animation: A AMS Fall Sectional Sampler Federico Ardila Continuum Approach for Braids, Surfaces, and Homological Frictional Contact Invariants J. Elisenda Grigsby Joseph Teran Ricardo Cortez presents our issue for Hispanic Heritage Month (page 901), which begins September 15. Federico Ardila shows how matroid theory has led to the solution of long-standing questions. Joseph Teran explains how movie animations model skin, clothing, and snow-covered ground. In a sampler of her address at the Eastern Sectional at Delaware this month, Elisenda Grigsby describes some invariants for knots and links, with possible application to the slice ribbon conjecture.
    [Show full text]
  • Cesar Ceballos –
    Institut für Geometrie TU Graz B [email protected] http://www.mat.univie.ac.at/~cesar/ Cesar Ceballos https://orcid.org/0000-0002-7348-3248 Personal information Born in Colombia on October 9, 1987. Married, two children (aged 6 & 8). Main areas of research Discrete geometry, polytope theory, algebraic, enumerative and geometric combinatorics, Hopf algebras, diagonal harmonics, representation theory, and tropical geometry. Academic career and positions held Current Privatdozent, Technische Universität Graz, Austria. May 2020 Privatdozent (Habilitation Degree), Universität Wien, Austria. 2015 – 2020 Postdoc, Universität Wien, Austria. 2013 – 2015 Visitor Member, Fields Institute, Canada. 2013 – 2015 Banting Postdoctoral Fellow, York University, Canada. 2012 – 2013 Postdoc, Freie Universität Berlin, Germany. 2009 – 2012 Dr.rer.nat., Freie Universität Berlin, Germany, Summa Cum Laude. Title – On associahedra and related topics. Supervisor – Günter Ziegler. 2008 – 2009 Master: Aspects of Calabi-Yau Geometries, Universiteit Utrecht, The Netherlands. 2007 – 2008 Master degree in mathematics, Universidad de Los Andes, Colombia. 2002 – 2006 Bachelor degree in mathematics, Universidad del Valle, Colombia. Selected awards 2021-2025 FWF Stand Alone Project Grant. 2018 Elsevier Mathematical Sciences Sponsorship Fund. 2013-2015 Banting Postdoctoral Fellowship, Government of Canada’s most prestigious award for postdoctoral researches. 2009-2013 Doctoral and Postdoctoral Fellowships, RTG Methods for Discrete Structures. 2008-2009 MRI Master Class Fellowship. 2008 Student Award for Art and Science, Mazda Foundation. 2008 Student Research Project Grant, Proyecto Semilla, Universidad de los Andes. 2007-2008 Master and Teaching Assistant Fellowship, Universidad de los Andes. CV Cesar Ceballos 2007 Second best graduating undergraduate student, Universidad del Valle. 2006 Bronze Medal, IX Ibero-American Mathematical Olympiads for University Students.
    [Show full text]
  • FWF Annual Report 2014
    Annual Report 2014 2014 Annual Report Annual Report 2014 We strengthen science and the humanities in Austria. Annual report submitted to the Austrian Federal Ministry of Science, Research and Economy in accordance with Art. 4 para. 1 of the Austrian Research and Technology Funding Act (FTFG). Vienna, March 2015. CONTENTS Table of contents 4 The FWF’s corporate policy 49 Bodies of the FWF Supervisory Board, FWF Management, FWF Board, 6 REPORT OF THE MANAGEMENT Assembly of Delegates, International START/ On the state of scholarly research in Austria Wittgenstein Jury, SciComm Jury, PEEK Board 20 GENERAL ACTIVITY REPORT 59 FWF Office Development of funding activities 63 Balance sheet, income statement, 28 FWF funding categories notes on the financial statements 31 APPENDIX / TABLES 68 Index of tables and illustrations ANNUAL REPORT 2014 3 INTRODUCTION The FWF’s corporate policy »We strengthen science and the humanities in Austria « The Austrian Science Fund (FWF) is Our objectives To emphasise and enhance the inter- Austria’s central funding organisa- To strengthen Austria’s international active effects of science and research tion for basic research. performance and capabilities in science with all other areas of culture, the and research as well as the country’s at- economy and society, and in particu- Our mission tractiveness as a location for high-level lar to increase the acceptance of sci- The purpose of the FWF is to support scientific activities, primarily by funding ence and research through concerted the ongoing development of Austrian top-quality research projects for individ- public relations activities. science and basic research in line with uals and teams and by enhancing the high international standards.
    [Show full text]