Notes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.

Similar documents
5. Analytic Combinatorics

DISCRETE AND ALGEBRAIC STRUCTURES

ADDITIONAL NOTES. 3 Surjections 6

Unlabelled Structures: Restricted Constructions

DISCRETE AND ALGEBRAIC STRUCTURES

On (Effective) Analytic Combinatorics

Generating Functions

7 Asymptotics for Meromorphic Functions

4. Complex Analysis, Rational and Meromorphic Asymptotics

Considering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.

LECTURE 7. 2 Simple trees and Lagrange Inversion Restricting the out-degree Lagrange Inversion... 4

3. Generating Functions

Algorithms for Combinatorial Systems: Between Analytic Combinatorics and Species Theory.

De part en retraite de Pierre

Combinatorial Enumeration. Jason Z. Gao Carleton University, Ottawa, Canada

ASYMPTOTIC ENUMERATION OF PERMUTATIONS AVOIDING GENERALIZED PATTERNS

1 Solving Linear Recurrences

Fourth Week: Lectures 10-12

INTRODUCTION TO REAL ANALYTIC GEOMETRY

Complex Analysis Math 147 Winter 2008

A BIJECTION BETWEEN WELL-LABELLED POSITIVE PATHS AND MATCHINGS

Generating Functions

Qualifying Exam Complex Analysis (Math 530) January 2019

Part IB. Further Analysis. Year

PROBABILITY VITTORIA SILVESTRI

are the q-versions of n, n! and . The falling factorial is (x) k = x(x 1)(x 2)... (x k + 1).

Combinatorial/probabilistic analysis of a class of search-tree functionals

COMPLEX ANALYSIS Spring 2014

Generating Function Notes , Fall 2005, Prof. Peter Shor

Conformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.

Contents. Counting Methods and Induction

LECTURE Questions 1

1 Basic Combinatorics

F (z) =f(z). f(z) = a n (z z 0 ) n. F (z) = a n (z z 0 ) n

NATIONAL UNIVERSITY OF SINGAPORE Department of Mathematics MA4247 Complex Analysis II Lecture Notes Part II

MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5

ANALYTIC COMBINATORICS COMPLEX ASYMPTOTICS

Growing and Destroying Catalan Stanley Trees

k-protected VERTICES IN BINARY SEARCH TREES

Complex Analysis Slide 9: Power Series

Riemann sphere and rational maps

PROBABILITY. Contents Preface 1 1. Introduction 2 2. Combinatorial analysis 5 3. Stirling s formula 8. Preface

Definitions. Notations. Injective, Surjective and Bijective. Divides. Cartesian Product. Relations. Equivalence Relations

or E ( U(X) ) e zx = e ux e ivx = e ux( cos(vx) + i sin(vx) ), B X := { u R : M X (u) < } (4)

Foundations of Mathematics MATH 220 FALL 2017 Lecture Notes

Solutions to Complex Analysis Prelims Ben Strasser

2 Write down the range of values of α (real) or β (complex) for which the following integrals converge. (i) e z2 dz where {γ : z = se iα, < s < }

Complex Analysis Topic: Singularities

Analytic Combinatorics: A Primer

INTRODUCTION TO COMPLEX ANALYSIS W W L CHEN

Symbolic enumeration methods for unlabelled structures

Lecture 9. = 1+z + 2! + z3. 1 = 0, it follows that the radius of convergence of (1) is.

Mapping the Discrete Logarithm: Talk 2

4 Uniform convergence

Complex Analysis MATH 6300 Fall 2013 Homework 4

arxiv:math.cv/ v1 23 Dec 2003

MASTERS EXAMINATION IN MATHEMATICS

Random Boolean expressions

A = A U. U [n] P(A U ). n 1. 2 k(n k). k. k=1

Binomial Coefficient Identities/Complements

Reading 11 : Relations and Functions

GENERALIZED STIRLING PERMUTATIONS, FAMILIES OF INCREASING TREES AND URN MODELS

A Generalization of Wigner s Law

Topic 4 Notes Jeremy Orloff

Complex Analysis Homework 9: Solutions

Complex Analysis Important Concepts

Solutions to practice problems for the final

Derangements with an additional restriction (and zigzags)

Combinatorial Hopf algebras in particle physics I

An Introduction to Combinatorial Species

Last Update: March 1 2, 201 0

PHYS 3900 Homework Set #03

THE RESIDUE THEOREM. f(z) dz = 2πi res z=z0 f(z). C

Lecture 17: Trees and Merge Sort 10:00 AM, Oct 15, 2018

Asymptotic Enumeration of Compacted Binary Trees

Part II. Riemann Surfaces. Year

What you learned in Math 28. Rosa C. Orellana

13 Maximum Modulus Principle

Random Boolean expressions

Enumerative Combinatorics 7: Group actions

On the parity of the Wiener index

and the compositional inverse when it exists is A.

Numerical Sequences and Series

Metric spaces and metrizability

INTRODUCTION TO REAL ANALYSIS II MATH 4332 BLECHER NOTES

Exercises for Part 1

III.2. Analytic Functions

Chapter 30 MSMYP1 Further Complex Variable Theory

MORE CONSEQUENCES OF CAUCHY S THEOREM

INTEGRATION WORKSHOP 2003 COMPLEX ANALYSIS EXERCISES

Equivalence Relations and Partitions, Normal Subgroups, Quotient Groups, and Homomorphisms

LECTURE-13 : GENERALIZED CAUCHY S THEOREM

Statistical Mechanics and Combinatorics : Lecture IV

Notes on Complex Analysis

MS 3011 Exercises. December 11, 2013

Lesson 8: Why Stay with Whole Numbers?

means is a subset of. So we say A B for sets A and B if x A we have x B holds. BY CONTRAST, a S means that a is a member of S.

Problems for MATH-6300 Complex Analysis

Review problems solutions

Problem 1. Arfken

Transcription:

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Notes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.. Tuesday, June 2, 202 Combinatorics is the study of finite structures that combine via a finite set of rules. Algebraic combinatorics uses algebraic methods to help you solve counting problems. Often algebraic problems are aided by combinatorial tools; combinatorics thus becomes quite interdisciplinary. Analytic combinatorics uses methods from mathematical analysis, like complex and asymptotic analysis, to analyze a class of objects. Instead of counting, we can get a sense for their rate of growth. Source: Analytic Combinatorics, free online at http://algo.inria.fr/flajolet/publications/ book.pdf. Goal.. Use methods from complex and asymptotic analysis to study qualitative properties of finite structures. Example.2. Let S be a set of objects labeled by integers {,..., n}; how many ways can they form a line? Answer: n!. How fast does this quantity grow? How can we formalize this notion of quickness of growth? We take a set of functions that we will use as yardsticks, e.g. line, exponential, etc.; we will see what function approximates our sequence. Stirling s formula tells us: ( n ) n n! 2πn. e a where a n b n implies that lim n b. This formula could only arrive via analysis; indeed Stirling s formula appeared in the decades after Newton & Leibniz. Almost instantly, we know that (00)! has 58 digits. Similarly, (000)! has 0 2568 digits. Example.3. Let an alternating permutation σ S n be a permutation such that σ() < σ(2) > σ(3) < σ(4) > > σ(n); by definition, n must be odd. Let T n be the number of alternating permutations in S n. Listing the first few terms of this sequence we have: n T n n! 3 2 2 5 6 20 7 272 5040 9 7936 362, 880 It appears that n! grows much faster than T n. How do we quantify this notion of much faster? () Find a formula to approximate T n. (2) Compare it to n! or the Stirling formula, i.e. we look at the ratio T n /n! as n. If the limit is, then the growth is as fast, if the limit is 0, then we say it is much slower. We want to find a formula that approximates T n ; in 88, Desire Andre found that the coefficient of zn n! in the function tan(z) is precisely T n. Therefore, tan(z) is the exponential generating function (e.g.f) for T n. We now describe a bijection between these permutations and the set of decreasing binary trees, i.e. a rooted binary tree where the value of a child is always less than the value of a parent. For example, taking the permutation (23456789 ), 57283649

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá we translate it into a decreasing binary tree by starting with the largest number, and letting the left child and the right child be the largest number to its left and its largest child to the right, then iterate. So, an alternating permutation looks like (σ L, max, σ R ). Let τ be the class of alternating permutations. Then, τ = () (τ, max, τ). T (z) = z + z 0 T (ω) 2 dω. dt (z) = + T (z) 2, T (0) = 0, dz a differential equation that we solve and find T (z) = tan(z). Remark.4. What we have learned: () A simple algorithm to find T n. (2) g.f admits an explicit expression in terms of a well known mathematical object, i.e. tan(z). (3) the probability of randomly obtaining an alternating permutation: T (n) n! is the coefficient of z n in Taylor expansion of tan(z). There is an amazing relationship between the growth of the coefficients of a counting sequence and the singularities of the generating function when viewed as an analytic function C C. Example.5. The sequence,,,,,... are the coefficients of the generating function: z n =, z <, z with singularity at z =. Example.6. The sequence, 2, 4, 8, 6,... are the coefficients of the generating function: 2 n z n = 2z, z < 2, with singularity at z = /2. Based on this, the singularities give us information about the growth of the sequence. Therefore, in Example.3, we want to examine the singularities of tan(z) near zero. The closest singularities are ± π 2. Let us approximate tan(z) near ± π 2, we have: tan(z) 8z π 2 4z 2. Extracting the coefficient of z n, we obtain: ( ) T n 2 n+ n! 2 0. π Therefore, we can say with more rigor that T n grows much slower than n!. Remark.7. Here, we used only real singularities. However, in general, we must study complex values of z. It is singularities in the complex plane that matter and complex analysis that is needed to draw conclusion regarding the asymptotic behavior of the counting sequence. 2

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Example.8. Let c n = the number of all unlabeled ordered binary trees with n internal nodes. n c n 0 2 2 3 5 4 4 5 42 You may recognize these as the Catalan numbers. We can obtain this via the symbolic method as well: C = (C,, C). C(z) = + zc(z) 2. C(z) = 4z. 2z = ( ) 2n z n. n + n As this is the ordinary generating function, we are concerned with the coefficient of z n, which gives us the formula: c n = ( ) 2n = (2n)! n + n n + n!n!. We can approximate this last expression using Stirling s formula: c n c n = 4n πn 3. Here the growth is related to a non-pole singularity, i.e. z = 4. 2. Wednesday, June 3, 202 / Symbolic Method (Note that all topics discussed today deal with unlabeled structures.) There are many set-theoretic operations that have a nice translation into generating functions. Today, we will discuss: Admissible constructions on combinatorial classes Symbolic Method Constructible combinatorial classes O.G.F. Definition 2.. A Combinatorial Class is a finite or enumerable set with a size function, such that: () The size of any element 0. (2) The number of elements of a given size is finite. A n counts number of elements of size n. Definition 2.2. An ordinary generating function µ A(z) = A n z n, where A is a combinatorial class. Equivalently, A(z) = α A z α. 3

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Admissible Constructions: m-ary construction: Associates to any collection of m combinatorial classes B (),..., B (m), a new class A = Φ(B (),..., B (m) ). It is admissible iff A n depends only on (B n () ) n,..., (B n (m) ) n. If we have an admissible construction, then there exists ψ such that A(z) = ψ[b () (z),..., B (m) (z)]. Example 2.3 (Cartesian Product). Example 2.4 (Disjoint Union of Sets). A = B C A = {α = (β, γ) β B, γ C}. n α A = β B + γ C, A n = B k C n k A(z) = B(z)C(z). k=0 A = B C, B C =. The size of ω A is just its size as an element of B or C. Remark 2.5. Even if A n = B n + C n A(z) = B(z) + C(z). Definition 2.6. A constructible class is a class that can be defined from primal elements by means of admissible constructions. Example 2.7. The neutral class E has one object of size 0. A = E A = A E A(z) =. The atomic class Z has one object of size. Its generating function is Z(z) = z. Basic admissibility: The combinatorial union, cartesian product, SEQ, PSET, MSET, CYC are all admissible with the operators on ordinary generating functions: () A = B + C. If the sets are not disjoint, we simply color them by taking the Cartesian product each with a particular color. Then, perform the usual disjoint union. The ordinary generating functions satisfy: A(z) = B(z) + C(z). (2) A = B C. As we saw before, A(z) = B(z)C(z). (3) SEQ(B) = {ɛ} + B + B B +. A = SEQ(B) = {α = (β,..., β l ) l 0}. The size function has α A = l i= β i B is admissible iff B has no objects of size 0. SEQ(Z) \ ɛ =,,,, simply the unary representation of integers. As for the ordinary generating functions, SEQ(B) : A(z) = + B(z) + B(z) 2 + = (4) A = P SET (B) := set of all finite subsets of B. The size is given by α = {{β,..., β l } β i B}. α A = 4 l β i. i= B(z).

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá A = P SET (B) := β B({ɛ} + {β}). A(z) = ( + z β ) = + z β B n ( n ) Bn. exp log A(z) = exp( B n log( + z n )) = exp( n k (5) MSET (B) := all finite multi sets of B. B is finite. A = MSET (B) = β B SEQ({β}). ( ) k B(z k )). k A(z) = z β. β B (6) Cycle construction CY C(B) := sequences taken up to circular shift. CY C(B) := SEQ(B)/S circular shift. Theorem 2.8 (Symbolic Method, Unlabelled). The generating function of a constructible class is a completion of a system of functional equations whose terms are built from:, z, +,, Q, Exp, Exp, log, where the non-obvious operations are defined by: Q[f] = f Exp[f] = exp( k Exp[f] = exp( k log[f] = k Iterative Combinatorial Classes: (SEQ). f(z k ) ) (MSET ). k ( ) k f(z k ) ) (P SET ). k ϕ(k) k log f(z k ). Example 2.9 (Binary Words). Let W be the set of binary words, i.e. words on the alphabet {a, b}. W = SEQ(z + z) W (z) = 2z. Example 2.0 (General Plane Trees). A plane tree has a root, and its children could be considered an ordered set of roots for new trees. Therefore, we solve recursively. G = ZSEQ(G) G(z) = z G(z) 2 ( 4z). Example 2. (Polygon Triangulations). An (m + 2)-gon can be tiled by m triangles. T = {ɛ} + (T,, T ). T (z) = + zt 2 (z) T (z) = 4z. 2z T n = C n. 5

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá 3. Thursday, June 4, 202 Symbolic Method for Labeled Structures. Recall that the number of trees on n + internal vertices satisfies: C n where the 4 is related to the singularity at /4. 4n πn 3, Definition 3.. A labelled class has labelled objects. () Weakly-labelled objects: labeled graph whose vertices I. (2) Well-labelled objects of size [n] if the set of labels is [... n]. A labelled class is a combinatorial class of well-labelled objects. Example 3.2 (Labelled graphs). The graph on four vertices 2 3 4 is equivalent to 4 3 2. Therefore, we have 2 labelled graphs on four vertices. Definition 3.3. The exponential generating function (E.G.F.) of a class A is given by: A(z) := z n A n n! = z α α!. α A The neutral class E and atomic class Z are the same as in the unlabeled universe. Three Basic Classes of Labels Example 3.4. The class of Permutations P, is given by: ( ) 2 P = {ɛ, (),,...}. 2 Clearly, P n = n!, so the exponential generating function is P (z) = Example 3.5. The class of Urns U, is given by: z n P n n! = z n = z. U = {ɛ, (), (2), (23),...}. U n =, so the exponential generating function is U(z) = z n U n n! = z n n! = ez. Example 3.6. The class of Circular graphs C, is given by: ( ) 23 C = {ɛ, (), (2),,...}. 32 C n = (n )!, so the exponential generating function is C(z) = z n C n n! = z n n = log z. Admissible Labeled Constructions 6

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Remark 3.7. Product of A(z) and B(z): A(z)B(z) = C(Z) C(z) = z n a n n! z n b n n!. C(z) = Remark 3.8. Product of multiple functions: A(z) = B () (z) B (r) (z) A(z) = ( n k=0 ( ) n z )a n k b n k k n!. k +...+k r=n ( n k, k 2,..., k r Definition 3.9. β γ := {(β, γ ) the well-labeled objects ρ(β ) = β, ρ(γ ) = γ}. Example 3.0. Product of a labeled triangle graph and a labeled edge graph. ) b () k b (r) zn k r Definition 3.. The labelled product B C is obtained by forming ordered pairs from B C and performing all possible order-preserving relabeling. A = B C = (β γ). A n = β + γ =n The list of admissible labeled constructions: β B,γ C ( ) β + γ A(z) = B(z)C(z). β () Labeled Product. (2) The set of k-sequences (k-seq), and the set of all sequences (SEQ). (3) The set of k-sets (k-set), and the set of all sets (SET). (4) The set of k-cycles (k-cyc), and the set of all cycles (CYC). (5) Combinatorial Sum. The corresponding generating functions are: () (2) (3) (4) (5) A(z) B(z) log B(z) e B(z) B(z) A(z) + B(z) 7 n!.

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Example 3.2 (Surjective functions). Consider R as the set of surjections f : A B, where B = r. Let R n (r) = all surjections f : [n] [r]. We can write this as the sequence (f(), f(2),..., f(n)) all of which are elements of [r]. We can also think of it as the sequence of subsets of [n] that map to k. For example, given ( ) 23456789 ({2}, {, 3}, {5, 7}, {9}, {4, 6, 8}). 22535354 Therefore, R = SEQ r (SET (Z)). This returns a generating function: Example 3.3 (Set Partitions). Let S (r) n R(z) = (e z ) r. = the number of ways to partition [n] into r blocks. S r = n S (r) n = SET r (SET (Z)). S (r) (z) = r! (ez ) r. Example 3.4 (Words). Let X = {a,..., a r } be an alphabet. Let W be the class of all words, and W n be the class of all words of size n. W = U r = SEQ r (U) W (z) = e rz. Example 3.5 (Non-Plane Labeled Rooted Trees). Let T be the class of all non-plane (i.e. subtrees are not ordered) labeled rooted trees. T = Z SEQ(T ). 4. Friday, June 5, 202 Example 4.. Recall that Alternating Permutations led us to the generating function tan(z) : C C. We also found that lim T (z) z ±π/2 8z π 2 4z 2. These coefficients grow similarly to 2 ( ) 2 n+ ( π = 4 2 ) n. π π Why does this tell us anything about the growth of the coefficients of tan(z)? The Symbolic Method has asked us to think as follows: Example 4.2. Sequence Formal Power Series Functions {A n } A n z n A(z). {,,,..., } z n z. {, 2, 4, 8,...} 2 n z n 2z. Today we will concern ourselves more with the functions on the right. 8

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá In general the sequence gives us a function that is analytic on a neighborhood of 0 in C. The analytic structure of A(z) tells us about the asymptotics of A n. Let f be an analytic function on a neighborhood of 0. Then f has a power series representation f(z) = f n z n. Singularities of f: Let f : Ω C. f is analytic on Ω if at every z 0 Ω, it has a power series representations f(z) = c n (z z 0 ) n on z z 0 < δ. Equivalently, f is differentiable at each z 0 Ω. A singularity of f is a point where f cannot be extended analytically. Example 4.3. {,,...} z n z. This power series converges only on Ω = z <, and the function is defined and analytic on C\{}. This leads to a singularity at z 0 =. Example 4.4. The exponential generating function for derangements: f n = [z n ]f(z), This has a singularity at z 0 =. f(z) = e z z D n = n!f n. Example 4.5. The exponential generating function for surjections: g(z) = 2 e z. This has a singularity when e z = 2, i.e. for all k Z, is a singular point. z = log 2 + 2πik Example 4.6. The ordinary generating function for Catalan numbers C n : h(z) = 4z. 2 ω ω has a singularity at ω 0 = 0. Therefore, this O.G.F. has a singularity at z 0 = 4. is a singular point. First Principle: The exponential growth of f n tells us the location of singularities of f. By location, we mean how close is the closest singularity to zero? Second Principle: The sub-exponential growth of f n tells us the nature of singularities of f. More precisely, suppose f(z) = f n z n has radius of convergence R such that 0 < R <. Then: ( ) n f n = Θ(n), where lim sup Θ(n) /n =. R n More explicitly, Θ(n) ( + ɛ) n for finitely many n, Θ(n) ( ɛ) n for infinitely many n. Fact 4.7. The closest singularity to 0 of f is on the circle z = R. One can justify this by assuming the negative and applying Cauchy s Integral Formula. Exercise 4.8. Prove: If Θ(n) := R n f n, then lim sup Θ(n) /n =. Hint: Root test. 9

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Example 4.9. f d = R Exp. z n 2z /2 2 n 4z ( log 2 2 /4 4 n 2 e z log 2 Transfer The second principle is called Transfer by Flajolet-Sedgwick. We want to know ) n f(z) σ(z) as z z 0? fn σ n as n. Strategy: () Find singularities of f(z) closest to zero. This gives you the dominant distance R. Then, ( ) n f n = Θ(n). R (2) Normalize to R = and z 0 =. (3) Find σ in catalog such that f σ as z z 0. (4) Conclude f n σ n. Example 4.0. Suppose f is analytic on C \ (, ); the closest singularity is at. Suppose also that f ( z) α = σ(z), for α =, 2, 3.... Claim: f n nα (α )!. We can evaluate a contour integral f n = 2πi γ f(z) dz. zn+ where γ is a Pacman shape with a little circle around and a big circle going out to infinity, with rays at ±π/4, as in Figure. π 4 Figure. Contour for Integral. Example 4.. Let σ = z, then σ n asymptotically approaches 2 πn 3. Therefore, for f = 2 ( 4z), we have f n 4 n 4 πn 3. 0

Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Example 4.2. Let σ = z, then σ n asymptotically approaches (and always equals). Let σ = ( z) 2, then σ n asymptotically approaches (and always equals) n.