1 COMPUTER ALGEBRA FOR LATTICE PATH COMBINATORICS
2 ALIN BOSTAN∗
3 Abstract. Classifying lattice walks in restricted lattices is an important problem in enumerative 4 combinatorics. Recently, computer algebra has been used to explore and to solve a number of diffi- 5 cult questions related to lattice walks. We give an overview of recent results on structural properties 6 and explicit formulas for generating functions of walks in the quarter plane, with an emphasis on the algorithmic methodology.7
8 Key words. Enumerative combinatorics, random walks in cones, lattice paths, generating functions, 9 computer algebra, automated guessing, creative telescoping, Gessel walks, algebraic functions, D-finite functions, elliptic functions.10
11 AMS subject classifications. Primary 05A10, 05A15, 05A16, 97N70, 33F10, 68W30, 14Q20; Secondary 33C05, 97N80, 13P15, 33C75, 12Y05, 13P05, 14Q20.12
13 This document is structured as follows. Section1 gives an overview of recent re- 14 sults obtained in lattice path combinatorics with the help of computer algebra, with 15 a focus on the exact enumeration of walks confined to the quarter plane. Sections2 16 and3 then go into more details of two classes of fruitful algorithmic approaches: 17 guess-and-prove and creative telescoping.
18 1. General presentation.
19 1.1. Prelude. Consider the following innocent-looking problem.
20 A tandem-walk is a path in Z2 taking steps from {↑, ←, &} only. 21 Show that, for any integer n ≥ 0, the following quantities are equal:
22 (i) the number an of tandem-walks of length n (i.e., using n steps), 23 confined to the upper half-plane Z × N, that start and end at (0, 0);
24 (ii) the number bn of tandem-walks of length n confined to the quar- 25 ter plane N2, that start at (0, 0) and finish on the diagonal x = y.
26 For instance, for n = 3, this common value is a3 = b3 = 3, as shown below.
(i)
27
(ii)
28 The problem establishes a rather surprising connection between tandem-walks 29 in the lattice plane, submitted to two different kind of constraints: the evolution 30 domain of the walk, and its ending point. The domain constraint is weaker for the 31 first family of walks, while the ending constraint is relaxed for the second family. 32 It appears that this problem is far from being trivial. Several solutions exist, 33 but none of them is elementary. One of the main aims of the present text is to
∗Inria and Université Paris-Saclay, 91120 Palaiseau, France ([email protected]). 1
This manuscript is for review purposes only. 2 ALIN BOSTAN
34 convince the reader that this problem (and many others with a similar flavor) can 35 be solved with the help of a computer. More precisely, Computer Algebra tools, 36 extensively described in the following sections, can be used to discover and to prove 37 the following equalities (3n)! 38 (1) a = b = , and am = bm = 0 if 3 does not divide m. 3n 3n n!2 · (n + 1)!
39 It goes without saying that such a simple and beautiful expression cannot be an 40 element of chance. As it will turn out, closed forms are quite rare for this kind of 41 enumeration problems. Nevertheless, even in absence of nice formulas, the struc- 42 tural properties of the corresponding enumeration sequences reflect the symmetries 43 of the step set and of the evolution domain. Equation (1) shows that the sequences 44 (an) and (bn) are P-recursive, that is, they satisfy a linear recurrence with polyno- 45 mial coefficients (in the index n). One of the messages that will emerge from the text 46 is that this important property of the enumeration sequences is intimately related 47 to the finiteness of a certain group, naturally attached to the step set {↑, ←, &}.
48 1.2. General context: lattice paths confined to cones. Let us put the previous 49 problem into a more general framework. Let d ≥ 1 be an integer (dimension), let S d d 50 be a finite subset (called step set, or model) of vectors in Z , and p0 ∈ Z (starting 51 point). A S-path (or S-walk) of length n starting at p0 is a sequence (p0, p1,..., pn) d 52 of elements in the lattice Z such that pi+1 − pi ∈ S for all 0 ≤ i < n. Let C be a 53 cone of Rd, that is a subset of Rd such that r · v ∈ C for any v ∈ C and r > 0, assumed 54 to contain p0. We will be interested in the (exact and asymptotic) enumeration of 55 S-walks confined to the cone C, and potentially subject to additional constraints. 56 Example 1. Consider the model S = {(1, 0), (−1, 0), (1, −1), (−1, 1)} (called the 57 Gouyou-Beauchamps model) in dimension d = 2, with starting point p0 = (0, 0) and 2 58 with cone C = R+ (the quarter plane). The picture below displays the step set of 59 the model (on the left), and a S-walk of length n = 17 confined to C (on the right).