Scripture Combinatorics
Scripture Combinatorics
Markus Fulmek
Summer Term 2018
The present (draft of) lecture notes is based on lectures “Combinatorics” held
since summer term 2012. It contains selected material from the following text
books:
‚ Species: Combinatorial Species and tree–like Structures by F. Bergeron,
G. Labelle and P. Leroux [1],
‚ Partially ordered sets: Enumerative Combinatorics I by R. Stanley [7],
‚ Asymptotics: Analytic Combinatorics by P. Flajolet and R. Sedgewick
[3].
The reader is assumed to be familiar with basic concepts from Discrete Math-
ematics (see, for instance, lecture notes [4]). Moreover, some knowledge from
Analysis (power series), Linear Algebra, Algebra (basic group theory) and Com-
plex Analysis (contour integrals) is a necessary precondition for the material
presented here.
This draft was translated to English by Ilse Fischer, Markus Fulmek, Christian
Krattenthaler and Michael Schlosser. Many thanks in advance for any hints
regarding typos and other flaws: most likely, this draft contains many of them!
Vienna, summer term 2019
Contents
Note that the condition b p0q “ 0 is necessary, since otherwise we could en-
counter infinite sums for the coefficients of the composition: we do not deal
with infinite sums in the calculus of formal power series.
T HEOREM 1.1.3. Let a pzq “ a1 ¨ z ` a2 ¨ z2 ` ¨ ¨ ¨ be a formal power series with a p0q “
0. Then there is a unique compositional inverse b pzq “ b1 ¨ z ` b2 ¨ z2 ` ¨ ¨ ¨ , i.e.
pa ˝ bq pzq “ pb ˝ aq pzq “ z,
if and only if a1 ‰ 0.
But how can we find this inverse? In fact, there is a quite general and very useful
formula which we shall derive in the following: For that, we consider a slightly
more general problem, which involves an extension of the calculus of formal
power series.
D EFINITION 1.2.1. A formal Laurent series (in some variable z) is a (formal) sum
ÿ
an ¨ zn
něN
with coefficients an P C, for some N P Z: N may be smaller than zero, i.e., there may
be negative powers of z.
We denote the set of all formal Laurent series by Cppzqq
Addition, multiplication and composition for formal Laurent series are defined
exactly as for formal power series.
T HEOREM 1.2.2. A formal Laurent series a pzq possesses a (unique) multiplicative
inverse if and only if a pzq ‰ 0.
This implies that the algebra of formal Laurent series is, in fact, a field.
In the following we shall always assume that the formal power series f pzq (with
f pzq “ 0) under consideration “starts with z1 , i.e., is of the form
f pzq “ f 1 ¨ z1 ` f 2 ¨ z2 ` ¨ ¨ ¨
f pzq
where f 1 ‰ 0 (we might express this as z ‰ 0): The general case
(Here, we made use of Iverson’s notation: r“some assertion”s equals 1 if “some asser-
tion” is true, otherwise it equals 01.)
` ˘
P ROOF. Observe that for n ‰ 0 we have n ¨ f n´1 pzq ¨ f 1 pzq “ p f n pzqq1 : For
n ą 0, this is a formal power series; and since the coefficient of z´1 equals 0 for
every formal power series, the assertion is true for all n ą 0.
If n “ ´m ă 0, then f n pzq “ c´m ¨ z´m ` ¨ ¨ ¨ ` c´1 ¨ z´1 ` c0 ` c1 ¨ z1 ` . . . and
p f n pzqq1 “ p´mq ¨ c´m ¨ z´m´1 ` ¨ ¨ ¨ ` p´1q ¨ c´1 ¨ z´2 ` c1 ` . . . ,
so the assertion is also true for all n ă 0.
1Iverson’s notation is a generalization of Kronecker’s delta: δ “ ri “ js.
i,j
1.2. LAGRANGE’S INVERSION FORMULA 3
P ROOF. For the first claim (1.2) we take the derivative of (1.1),which gives
ÿ
g1 pzq “ k ¨ ck ¨ f 1 pzq ¨ f k´1 pzq .
k
and the claim follows, again, by taking the coefficients of z´1 on both sides.
Now Lagrange’s inversion formula is the following special case of Theorem 1.2.4:
4 1. GENERATING FUNCTIONS AND LAGRANGE’S INVERSION FORMULA
k r ´k z ´n
Jzn K Fk pzq “ z f pzq for n ‰ 0. (1.4)
n
An alternative expression for these coefficients is
r z
´k´1
n k
Jz K F pzq “ z f 1 pzq f ´n´1 pzq . (1.5)
pkq 1 r ´1 z k r ´k z ´n
Fn “ z k ¨ zk´1 ¨ f ´n pzq “ z f pzq ,
n n
which is equivalent to (1.4). In the same way, applying (1.3) gives (1.5).
E XAMPLE 1.2.6 (Catalan numbers). Let Fn be the number of triangulations of an
pn ` 1q–gon (with vertices labelled 1, 2, . . . , n successively): Setting F1 “ 1 (i.e., as-
suming that there is a single triangulation of a 2–gon), the sequence pFn qně1 starts like
this:
1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, . . .
For n ą 1 and any such triangulation, consider the triangle containing the edge which
connects the vertices 1 and 2: This triangle “cuts” the pn ` 1q–gon in two parts, namely
a pk ` 1q–gon and an pn ´ k ` 1q–gon, where 1 ď k ď n ´ 1. This leads to the recur-
sion
n´1
ÿ
Fn “ Fk ¨ Fn´k for n ą 1,
k“1
pF pzqq2 “ F pzq ´ z.
we obtain:
1 r ´1 z
Jzn K F pzq “ z pz ´ z2 q´n
n
1 r ´1 z ´n
“ z z p1 ´ zq´n
n
1 r n´1 z
“ z p1 ´ zq´n
nˆ ˙
1 ´n
“ p´1qn´1
n n´1
ˆ ˙
1 2n ´ 2
“ .
n n´1
Clearly, Fn “ Cn´1 , i.e., we get (once again) the well–known Catalan–numbers.
CHAPTER 2
The theory of species provides a general framework for labelled and unlabelled
combinatorial objects and their generating functions. One can tackle this theory
very abstractly; however here we would like to approach the theme by means
of a couple of illustrative examples. (Much more material than treated here can
be found in the textbook [1].)
3 1 6 5 2 9 10 7 12 11 8 4
Here one could have the idea to “forget the numbering”. One could imagine
this graphically for the above example like this:
2 ` 3 `1` 2 ` 4
A simple “modification” would be the species sets1 of nonempty sets (of atoms); the
generating function is then
z 1
gfsets1 pzq “ z ` z2 ` z3 ` ¨ ¨ ¨ “ “ ´ 1.
1´z 1´z
R EMARK 2.3.3. Since the atoms here are unlabelled (i.e., not distinguishable), their
ordering does not matter: Therefore the species sets (of atoms) is equal to the species
sequences (of atoms). This however depends on the indistinguishability of the el-
ements (atoms): sequences of distinguishable elements are of course not the same
combinatorial objects as (multi-)sets of these elements!
D EFINITION 2.3.4. A tree is a simply connected graph without circles.
A rooted tree is a tree with a special vertex, which is referred to as root. A rooted
tree always has a “natural orientation” (“from the root away”), it thus emerges as a
directed graph: Vertices with outgoing degree 0 are called leaves, all other vertices
are called inner vertices.
Every tree is of course planar (is embeddable in the plane such that two arbitrary em-
bedded edges intersect in at most one vertex): A rooted tree in a fixed embedding in
the plane, where the order of the subtrees is relevant, is called planar oder ordered
rooted tree.
The size of a tree (in the sense of Definition 2.2.1) can be given by the number of its
vertices: Trees are thus a further example of species.
For instance, there are the following 2 planar rooted trees with 3 vertices:
With 4 vertices there are the following 5 planar rooted trees, although the sec-
ond and the third of them would be viewed as the same “normal” (i.e., “not-
planar”) rooted tree:
R EMARK 2.3.5 (Here always: rooted trees!). Especially for computer sciences algo-
rithms (e.g., search algorithms) and data structures which can be “composed” of trees
are interesting objects. Since these trees often have a “distinguished vertex” (a root), we
will in the following always consider rooted trees (if not explicitly speficied otherwise)
but for simplicity talk about trees.
As we have seen in the introductory Example 2.1.2, the combinatorial construc-
tion “form finite sequences of objects” transfers to the algebraic construction
“insert the generating function of the objects into the geometric series”. We
shall examine this “transfer”
combinatorial construction Ñ algebraic construction
12 2. SPECIES
in further examples:
D EFINITION 2.3.6 (Disjoint union). Given two species A and B , we can consider
their disjoint union C “ A Y
9 B ; the corresponding size function is
#
# pxq x P A,
#C pxq :“ A
#B pxq x P B .
řn
` ˘ k“0 m k ¨ wn´k
H, (Here: n “ 3.)
m 0 ¨ w3 “ 1 ¨ 2 “ 2
` ˘
H,
` ˘
t˝u, m 1 ¨ w2 “ 1 ¨ 1 “ 1
` ˘
t˝, ˝u, m 2 ¨ w1 “ 1 ¨ 1 “ 1
` ˘
t˝, ˝, ˝u, m 3 ¨ w0 “ 1 ¨ 1 “ 1
A special case is the k–fold cartesian product of a species A whose objects are sequences
of length k of objects from A:
A ˆ A ˆ ¨ ˆ A : gf Ak “ pgf A qk .
Ak “ loooooooomoooooooon
k times
If k “ 0, we obtain the single-element species A0 “ tǫu which only contains the empty
sequence ǫ “ pq with }ǫ} “ 0.
For instance, we can view the combinatorial objects “compositions of a number n with
exactly k parts” as species K (with size function n): This apparently corresponds to
ˆ ˙k
k
` ˘k z
K “ sets1 , thus gf K “ gfsets1 “ .
1´z
2.3. UNLABELLED SPECIES AND THE ENUMERATION OF TREES 13
D EFINITION 2.3.8 (Sequences). Given a species A which does not contain any object
of size 0, we can consider arbitrarily long finite sequences of objects from A, i.e., the
species A‹ ď
A‹ :“ Ak
kě0
According to the previous examples this transfers to generating functions as follows:
ÿ 1
gf A‹ “ pgf A qk “ .
1 ´ gf A
kě0
We have already seen a concrete examples for this in Example 2.1.2:
compositions “ sets‹1 .
R EMARK 2.3.9. Why did we have to assume here that A contains no objects of size 0?
Because otherwise A‹ would contain infinitely many objects of size 0, in contradiction
to our assumption for species:
} pq } “ } pǫq } “ } pǫ, ǫq } “ } pǫ, ǫ, ǫq } “ ¨ ¨ ¨ “ 0.
A PPLICATION 2.3.10. If one removes the root from a nonempty rooted tree W, a
(possibly empty) “list” of connected components arises which themselves consist of
nonempty rooten trees: If the rooted tree W is
‚ ordered, then one should also view the “list” as ordered (i.e., as sequence) —
then from the (ordered) sequence one can of course again uniquely reconstruct
the original (ordered) rooted tree W,
‚ unordered, then one should also view the “list” as unordered (i.e., as multi-
set) — then from the (unordered) multiset one can of course again uniquely
reconstruct the original (unordered) rooted tree.
The following graphic illustrates this simple thought:
´ ¯
n
1
, 2 , 3 ,...,
n
planar rooted tree
1 2 3
! )
2 n
, , 1
,..., 3
W
“normal” rooted tree
If we denote the number of planar rooted trees with n vertices with f n , then we
obtain by definition the generating function of the species of planar rooted trees
with at least one vertex:
ÿ
f pzq “ f n ¨ zn “ z ` z2 ` 2z3 ` 5z4 ` ¨ ¨ ¨
ně1
From the Application 2.3.10 we first get
planarrootedtrees 1 “ atom ˆ planarrootedtrees ‹1
and hence then the following functional equation:
´ ¯ z
f “ z ¨ f0 ` f1 ` f2 ` ¨¨¨ “ ðñ f 2 ´ f ` z “ 0.
1´ f
14 2. SPECIES
This quadratic equation for f has two solutions of which the following has the
“right” (positive) coefficients:
?
1 ´ 1 ´ 4z ÿ 1 ˆ2n˙
f pzq “ “ zn`1 . (2.4)
2 ně0
n ` 1 n
If an identity for analytic functions also makes sense for the corresponding for-
mal power series, then it is automatically also an identity for formal power
seriesn.
Conversely: If an identity for formal power series also makes sense for the cor-
responding analytic functions (this means that there exists a nontrivial common
radius of convergence for all involved series), then it is automatically also an
identity for analytic functions.
E XAMPLE 2.3.11 (Sets). Let A be a species which does not contain any object of size 0.
If we in sequences of A “forget about the order”, we obtain the species C “ sets pAq
(as said: actually multisets — elements can occur repeatedly): the size function for an
object
tA1 , A2 , . . . , Ak u P C
is given (just as in sequences) by
k
ÿ
#C ptA1 , A2 , . . . , Ak uq “ #A pAi q .
i“1
{ } { }
{ }
B:
{ }
{ } { }
A pB q:
A:
{ }
However, unfortunately this construction does not translate to the composition of the
respective generating functions! Foir instance, we have:
sets psets1 q “ numberpartitions,
but
ź 18
1
z ‰ .
1 ´ 1´z 1 ´ zi
i“1
(The problem is apparently: the construction as such is not unique — in the above visu-
alization we could also have distributed the objects from B differently among the atoms
of the objects from A. We would need to, say, first construct all possible “insertions
of objects in atoms”, and then figure out which of those are actually different!)
We would like to determine the number of complete binary ordered trees with
n inner vertices: This is the same as the number of ordered binary trees with n
vertices. The following graphic shows the objects of this species B with size
(number of vertices) n “ 3:
Hint: Find an equation for the generating function and use Langrange’s inversion formula.
Exercise 2: How many ways are there to (properly) parenthesize n pairwise non–commuting elements
of a monoid? And how does this number change if the n elements are pairwise commuting?
For example, consider 6 non–commuting elements x1 , x2 , . . . , x6 . Two different ways to parenthesize
them properly would be
ppx2 x5 qppx1 px4 x6 qqx3 qq and ppx3 px1 px4 x6 qqqpx5 x2 qq.
However, these would be equivalent for commuting elements.
Hint: Translate parentheses to labelled binary trees: The outermost pair of parentheses corresponds
to the root, and the elements of the monoid correspond to the leaves.
It is noticeable that the generating functions for nonempty planar rooted trees
(2.4) and ordered binary trees (2.5) are equal up to factor z, and the numbers of
these objects with fixed size (number of vertices) is essentially a Catalan number
` ˘
1 2k
Ck “ k`1 k : This can also be shown using bijections.
4
u
3 d
u u u u
2 d d d d
u u u
1 d d d
u
0 d
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
2.4.1.1. Planar trees Ø Dyck paths. From the graphic it becomes apparent
that a Dyck path can be coded as as word of length 2n over the alphabet tu, du,
where
‚ the number of u’s is equal to the number of d’s, i.e. n,
‚ and in each “initial part” of the word at least as many u’s occur as d’s.
For the above graphic this “code word” is:
uududuududuudduddd.
This immediately translates to a coding for ordered trees (plant louse coding).
Each tree on n ` 1 vertices has exactly n edges: We imagine a “plant louse”
which sequentially walks along the edges of the ordered rooted tree (always
starting with the left subtree).
ud
u
ud
d
udu
ud
ud
ud
ud
2.4.1.2. Ordered binary trees Ø Planar trees. The corresponding rotation corre-
spondence is best made clear in a graphic:
rotate by 45˝ : new root: new edges:
E XAMPLE 2.6.1. The labelled species atom again only consists of a single element, for
which there is clearly just one labeling: “˝ is labelled by 1”. The generating function is
hence again
z
egf atom pzq “ “ z.
1!
E XAMPLE 2.6.2. A further species is sets (of atoms): In the labelled objects the atoms
are distinguishable but, since in a set the order does not matter, there is exactly one
labelled set with n elements:
sets “ tH, t1u , t1, 2u , . . . u ,
the (exponential) generating function thus is
z z2
egf sets pzq “ 1 ` ` ` ¨ ¨ ¨ “ exp pzq .
1! 2!
The (exponential) generating function of sets1 then is
z z2 z3
` ` ` ¨ ¨ ¨ “ exp pzq ´ 1.
1! 2! 3!
E XAMPLE 2.6.3. Since the atoms here are labelled, the (labelled) species sets (of
atoms) is not anymore equal to the (labelled) species sequences (of atoms): The lat-
ter is simply einfach equal to the (labelled) species permutations with the generating
function
1! 2! 1
egfpermutations pzq “ 1 ` z ` z2 ` ¨ ¨ ¨ “ .
1! 2! 1´z
20 2. SPECIES
The species cycles of the cyclic permutations (i.e. of the permutations whose unique
cycle decomposition consists of a single cycle) apparently has the generating function
0! 1! 2! 1
egfcycles pzq “ z ` z2 ` z3 ` ¨ ¨ ¨ “ log .
1! 2! 3! 1´z
Just as with the unlabelled species there is also here a “transfer”
combinatorial construction Ñ algebraic construction.
E XAMPLE 2.6.4 (Disjoint union). Given two labelled species A and B , we can (just
as with the unlabelled species) consider their disjoint union C “ A Y
9 B , and also here
one has for the generating functions:
egf C “ egf A ` egf B .
E XAMPLE 2.6.5 (Product). Given two labelled species A and B , we define their prod-
uct C “ A ‹ B as follows: We first consider the cartesian product A ˆ B . From
`m`n ˘ a pair of
objects pA, Bq P A ˆ B (where }A} “ m and }B} “ n) we construct n ddifferent
labelled objects according to the following rule:
‚ Consider all decompositions of the set rm ` ns in two disjoint parts
rm ` ns “ X Y
9 Y with |X| “ m, |Y| “ n,
‚ For every such decomposition replace the numbers from A resp. B according
to the order-preserving bijections2
rms Ñ X
rns Ñ Y
We illustrate this by an example:
2 3 1
pA, Bq “ 1 2
}A}`}B}“3`2“5:
2 3 4
tt1,2,3u ,t4,5uu Ñ 1 5
2 4 3
tt1,2,4u ,t3,5uu Ñ 1 5
4 5 1
tt3,4,5u ,t1,2uu Ñ 3 2
A ‹ A ‹ ¨ ¨ ¨ ‹ A : egf Ak “ pegf A qk .
Ak :“ loooooooomoooooooon
k times
For instance:
‚ setsk “ ordered decomposition in k blocks, which one could also in-
terpret as “functions to rks”.
‚ sets1 k are thus to be viewed as “surjective functions to rks”.
E XAMPLE 2.6.6 (Sequences). Also for a labelled species A which contains no object
of size 0 we can consider arbitrary finite sequences of objects from A, i.e. the species A‹
ď
A‹ :“ Ak .
kě0
For instance, we obtain for the species of surjections τ : rns Ñ rks (with }τ} “ n)
surjections “ sets1 ‹ .
It immediately follows: The exponential generating function of the numbers3
n
ÿ
Sn,k ¨ k!
k“0
(where Sn,k denotes the Stirling numbers of the second kind) thus is
1 1
“ .
1 ´ pe z ´ 1q 2 ´ ez
E XAMPLE 2.6.7 (Sets). If in an object from A‹ we “forget about the order”, we obtain
the species
C “ sets pAq .
3These are the numbers of all surjections from rns Ñ rks, where k “ 0, 1, . . . , n is arbitrary.
22 2. SPECIES
pegf B pzqqk
.
k!
‚ Take an A–object A of size #A pAq “ k and imagine the atoms of A “blown
up”: In the labelled atoms now insert in the “canonical order” the objects
B j1 , ¨ ¨ ¨ , B jk (i.e., in the A–atom with number 1 one inserts B j1 , in the A–atom
with number 2 one inserts B j2 , etc.).
‚ The size of the thus produced object from C P C :“ A pB q is then of course
given by
#C pCq “ #B pB1 q ` ¨ ¨ ¨ ` #B pBk q .
2.6. LABELLED SPECIES 23
Here the construction indeed translates to the composition of the generating functions:
ÿ ÿ 1 ÿ z}Bj1 }`¨¨¨`}Bjk }
egf C pzq “ ` ˘
k! !
kě0 APA
) }B j1 } ` ¨ ¨ ¨ ` }B jk } !
}A}“k B j1 ,...,B j PB k
k
ÿ pegf B pzqqk
“ ak “ egf A pegf B pzqq .
k!
kě0
Exercise 3: Develop a theory for weighted generating functions (for labelled and unlabelled species).
I.e., let A be some species with weight function ω which assigns to every object A P A some element
in a ring R (for instance, R “ Z rys, the ring of polynomials in y with coefficients in Z). So the
generating function to be considered is
ÿ
z} A} ¨ ω pAq .
A PA
How should we define the weight function for sums, products and composition of species, so that the
corresponding assertions for generating functions remain valid?
Exercise 4: Let f pm, nq be the number of all paths from p0, 0q to pm, nq in N ˆ N, where each single
steps is either p1, 0q (step to the right) or p0, 1q (step upwards) or p1, 1q (diagonal step upwards).
Use the language of species to show that
ÿ 1
f pm, nq x m yn “ .
1´x´y´x¨y
m,ně0
E XAMPLE 2.6.10 (Cayley’s theorem). Let T denote the species of nonempty unorderd
labelled rooted trees. With the same consideration as in Section 2.5 we obtain the same
“species equation” (2.9) also for the labelled rooted trees, thus
T “ t˝u ˆ sets pT q .
This equation here translates according to Example 2.6.7 resp. 2.6.9 to
egf T pzq “ z ¨ eegfT pzq ,
hence egf T pzq is the compositional inverse of z ¨ e´z . We can very easily determine
their coefficients using the Lagrange inversion formula (Corollary 1.2.5: (2.7) for k “
1):
1 r ´1 z ` ˘´n
Jzn K gf T pzq “ z z ¨ e´z
n
1 r ´1 z ´n nz
“ z z e
n
1 r n´1 z nz
“ z e
n
1 nn´1
“
n pn ´ 1q!
for n ą 0. This means, the number of labelled unordered rooted trees on n ą 0 vertices
is
tn “ nn´1 .
24 2. SPECIES
If we “forget” about the root, we then obtain Cayley’s formula for the number of all
labelled trees (not rooted trees!):
# plabelled treeson n verticesq “ nn´2 . (2.10)
Exercise 5: Determine the number of all unlabelled ordered binary rooted trees with n vertices and
k leaves.
Hint: Consider the generating function in 2 variables z and y, where every rooted tree W with
n vertices and k leaves is assigned ω pWq :“ zn yk . The following picture shows these trees for
n “ 1, 2, 3, 4:
n“1:
n“2:
n“3:
n“4:
Exercise 6: Determine the number of all labelled unordered rooted trees with n vertices and k leaves.
Hint: Consider the exponential generating function in 2 variables z and y (as in Exercise 5) and use
Lagrange’s inversion formula.
and it is clear that the root has number 2: With this the rooted tree can be
uniquely reconstructed.
Exercise 7: Prove Cayley’s formula (the number of labelled trees on n vertices equals nn´2 ) as
follows: Take a labelled tree on n vertices and tag two vertices S and E. View S and E as the
starting point and ending point of the unique path p connecting S and E in the tree. Now orient
all edges belonging to p “from S to E”, and all edges not belonging to p “towards p”. Now travel
along p from S to E and write down the labels of the vertices: Whenever a new maximal label is
encountered, close a cycle (by inserting an oriented edge from the vertex before this new maximum
to the start of the “current cycle”) and start a new cycle. Interpret the resulting directed graph as a
function rns Ñ rns (i.e., a directed edge from a to b indicates that the function maps a to b).
Exercise 8: Show that the number of all graphs on n vertices, m edges and k components equals
the coefficient of un αm βk {n! in
˜ ¸β
ÿ n
pn2 q u
p1 ` αq .
n!
n ě0
Hint: Find a connection between the generating function of all labelled graphs (weight ω pGq :“
u|V pGq| α|EpGq| ) and the generating function of all connected labelled graphs.
Exercise 9: Show that the number of labelled unicyclic graphs (i.e., connected graphs with exactly
one cycle) on n vertices equals
n ˆ ˙
1ÿ n
j! nn´1´ j .
2 j
j “3
Hint: For a bijective proof consider the species W of labelled trees, where the vertex with the largest
label (i.e.: n, if the tree has n vertices) is tagged as the root, but this label is erased, and the root
does not contribute to the size of the tree (i.e: if t has n vertices (including the root), then we have
}t}W “ n ´ 1): Obviously, T “ gfW .
Now consider functions f : rls Ñ rl ` js. Visualize such function f as a directed graph with vertex
set rl ` js and directed edges px, f pxqq. The following graphic illlustrates this for the case l “ 20,
j “ 2 and
p f pnqqln“1 “ p2, 6, 2, 5, 2, 7, 8, 5, 4, 4, 7, 21, 21, 18, 16, 18, 22, 15, 15, 15q :
9
1 3 12
4 10
2 13 21
6 5 14
7 8 16 18
19
11
20 15 17 22
The components of this graph are rooted trees or unicyclic graphs; all vertices in rl ` js z rls appear
as roots of corresponding trees.
Exercise 11: The derivative A1 of a labelled species A is defined as follows: Objects of species A1
with size n ´ 1 are objects of A with size n, whose atoms are numbered from 1 to n ´ 1 (not from
1 to n), such that there is one atom without a label. A typical element of Sequences1 is
p3, 1, 2, 5, ˝, 4q ,
where ˝ indicates the unlabelled atom.
Show: The generating function of A1 is precisely the derivative of the generating function of A.
Moreover, show:
(1) pA ` B q1 “ A1 ` B 1.
(2) pA ‹ B q1 “ `A1 ‹ B `
˘ B 1 ‹ A.
1
(3) pA ˝ B q “ A1 ˝ B ‹ B 1.
These equations are to be understood as size–preserving bijections.
Exercise 13: Let A be the (labelled) species of (unordered) rooted trees, U the (labelled) species of
trees (without root) and F the (labelled) species of rooted forests. Show the following equations:
(1) A1 “ F ‹ Sequences pAq,
(2) U 2 “ F ‹ A1,
` ˘2 ` ˘2
(3) A2 “ A1 ` A1 ‹ Sequences pAq.
n“4
k“2
n“5
k“2
3 2 4
S3 xǫ,p1 2q,p3 4q,p2 4qp1 3qp5 6qy
5 6
1 2 1 3
Further, let R “ rks be a given set of k colours. We are interested in the number of all
nonequivalent mappings f : rns Ñ R, where f 1 „ f 2 if f 1 ˝ g “ f 2 for any g P G.
With other words: We search for the number of all equivalence classes of mappings
under the equivalence relation
f 1 „ f 2 :ðñ Dg P G : f 1 ˝ g “ f 2 .
We call these equivalence classes of “colourings” patterns.
Let each element r P R be assigned a weight śnω prq (usually a monomial, e.g. simply
r
z ), the weight of a function be ω p f q :“ i“1 ω p f piqq. It is clear: f 1 „ f 2 ùñ
ω p f 1 q “ ω p f 2 q; the weight of a pattern shall be ω p f q, where f is any representative
of the pattern (i.e. of the equivalence class).
As always we are interested in the generating function of all patterns, i.e. in the
sum ÿ
ω pfq.
pattern
E XAMPLE 2.7.4. We consider the example of the cube with ω preadq “ x, ω pblueq “ 1:
The weight of a pattern is hence x#pred facesq . The generating function is then:
x6 ` x5 ` 2x4 ` 2x3 ` 2x2 ` x ` 1.
L EMMA 2.7.5 (Burnside’s Lemma, Theorem of Cauchy–Frobenius). Consider the
problem from Definition 2.7.3. Let ψα pgq be the number of all functions f with weight
ω p f q “ α and f ˝ g “ f . Then the number of all patterns with weight α is equal to
1 ÿ
ψα pgq .
|G|
gPG
I.e., for the elements of a fixed right coset of G f one always gets the same func-
tions. Hence we have:
ˇ ˇ ˇ ˇ ˇ ˇ ˇ ˇ
ˇ ˇ ˇ ˇ ˇ ˇ ˇ ˇ
|G| “ ˇ f ˇ ¨ G f ùñ η p f q “ G f “ |G| { ˇ f ˇ .
2.7. PÓLYA’S THEOREM 29
T HEOREM 2.7.9 (Pólya’s theorem). Let G Ď Sn . The generating function for all
patterns f , f : rns Ñ R, is
˜ ¸
ÿ ÿ ÿ ÿ
ω pmq “ PG ω prq , ω prq2 , ω prq3 , . . .
m pattern rPR rPR rPR
R EMARK 2.7.10. If the set R of colours contains only one element which is simply
assigned the weight z, then there is of courses only one pattern, and the generating
2.8. A GENERALISATION 31
function of this single pattern is simply zn , where n is the size of the object — and
indeed also ´ ¯
PG z, z2 , . . . “ zn
holds.
E XAMPLE 2.7.11. For the cube we indeed obtain:
ˆ ´ ¯ ´ ¯2 ´ ¯2
1
6 x4 ` 1 px ` 1q2 ` 8 x3 ` 1 ` 3 x2 ` 1 px ` 1q2 `
24
´ ¯3 ˙
6 x ` 1 ` px ` 1q “ x6 ` x5 ` 2x4 ` 2x3 ` 2x2 ` x ` 1.
2 6
E XAMPLE 2.7.12. For a chain of length n which is coloured with k colours, we assign
the weight xi to the colour i and obtain:
1ÿ ´ ¯n{d
ϕ pdq x1d ` x2d ` ¨ ¨ ¨ ` xkd .
n
d|n
Exercise 15: How many different necklaces of n pearls in k colours are there? (This should be
understood “as in real life”, where rotations and reflections of necklaces are considered equal; in
contrast to the presentation in the lecture course.)
2.8. A generalisation
We can apparently generalise Definition 2.7.3, by considering also on the set of
colours R a permutation group H which acts on R: How many patterns are there
which under the action of H are yet different? I.e, we now consider two patterns
f 1 , f 2 as equivalent if there is a g P G (G is the perrmutation group on the atoms)
and a h P H (H is permutation group on the colours) such that
f 1 „ f 2 :ðñ f 1 ˝ g “ h ˝ f 2
E XAMPLE 2.8.1. In Example 2.7.1 (colourings of the faces of the cube with colours
R “ tblue, redu), let the group be H “ SR : I.e., there are now only 6 possibilities for
“essentially different” patterns:
32 2. SPECIES
The same object, viewed as ordered rooted tree, would however aber only be fixed by
the identity (compare also with Remark 2.6.8).
D EFINITION 2.9.3 (Cycle index series). Let F be a labelled species. For a permuta-
tion σ P Sn we define:
fixF pσq :“ tF P Fn : σ pFq “ Fu .
|fixF pσq| is hence the number of objects F in F for which σ belongs to the symmetry
group of F.
The cycle type of a permutation σ we denote by
c pσq :“ pc1 , c2 , . . . , cn q .
(I.e., ci is the number of cycles of length i in σ.)
Further, for a sequence x :“ px1 , x2 , . . . q of infinitely many variables we define:
xcpσq :“ x1c1 ¨ x2c2 ¨ ¨ ¨ xncn .
The cycle index series of a species F is then the formal power series (in infinitely
many variables)
ÿ xcpσq
ZF pxq “
}F}!
σ,F
σpFq“F
ÿ 1 ÿ
“ |fixF pσq| ¨ xcpσq (2.13)
ně0
n! σPS n
ÿ 1 ÿ ÿ
“ xcpσq
n!
ně0 σPSn FPFn :
σpFq“F
ÿ 1 ÿ ÿ
“ xcpσq . (2.14)
n!
ně0 FPFn σPSn :
σpFq“F
R EMARK 2.9.4. Unordered labelled rooted trees differ from ordered labelled rooted trees
“only” by the property that the first have a nontrivial symmetry group (i.e..: not only
consisting of the identity), compare with Example 2.9.2. With other words: For the
species of ordered labelled rooted trees
W1 „ W2 :ðñ “W1 and W2 are identical as unordered objects”
is clearly an equivalence relation whose equivalence classes can be naturally identi-
fied with the unordered labelled rooted trees. “From the view of the unordered trees”
this equivalence relation is given by its symmetry group G (which consists of “rear-
rangements” of subtrees):
W1 „ W2 :ðñ D σ P G : σ pW1 q “ W2 .
This can be generalized as follows: Let Fp be a labelled specied which contains only
the trivial symmetry group (i.e., all n! relabelings of an object Fp P Fp are viewed as
different): ´ ¯
sym p Fp “ tidu for all Fp P Fp .
F
34 2. SPECIES
Let G be a “larger” (casually spoken) symmetry group for the objects from Fp ; i.e.,
tidu Ă G “ GFp Ă Sn for all Fp P Fpn Ď Fp .
Then one can consider the species F of the equivalence classes of Fp with respect to
this symmetry group, whose objects F P F are determined by the prpoerty that for any
two representatives F1 , F2 P F Ď Fp
F1 „ F2 :ðñ D σ P G : σ pF1 q “ F2 (2.15)
holds. For a Fp P Fpn let
´ ¯ ! ´ ¯ )
Sn Fp :“ σ Fp : σ P Sn (all “relabelings” of F),
p
´ ¯ ´ ¯ ´ ¯
´ ¯ τ2 Fp τ3 Fp τm Fp
Fp
r
F “ Sn Fp “
F“GF p F
pq τ3 pFq“pτ3 GF τ3´1 qpτ3 p Fpqq
Casually spoken: Generally labelled species F appear “in the middle between the cor-
responding unlabelled species Fr and the corresponding totally unsymmetric species
Fp ”; expressed by the corresponding symmetry groups:
´ ¯ ´ ¯
symFp Fp “ tidu Ă GF “ symF pFq Ă Sn “ symFr Fr .
E XAMPLE 2.9.5. To make this completely concrete, we consider again the unordered
labelled rooted tree from Example 2.9.2:
p
“Fix numbering of F”:
1 2
Fp 1 Fp 2
3
3
GF “ tid, p3q p12qu » Z2 .
´ ¯
Consider S6 Fp , i.e., all reorderings of F:
p
σ“p123q τ“p1qp23q id
2 3 1 3 1 2
1 2 1 2 1 2
´ ¯ 1 2 3
S6 Fp “ 3 3 2 3 3 1 2 3 1
1 2 1 2 1 2
1 2 3
3 3 3
σ pFq τ pFq F
L EMMA 2.9.6. With the notation introduced in Remark 2.9.4 there holds: The cycle
indikator series of a labelled species F is the sum over the cycle indices of all objects of
the corresponding unlabelled species Fr . This means:
ÿ
PGF pxq “ ZF pxq . (2.16)
r Fr
FP
36 2. SPECIES
Suppose i and π piq “ j belong to cycles of different length in σ: W.l.o.g. let the σ–cycle
of i be the shorter one (otherwise consider π ´1 instaed of π); let k be its length:
σ σ σ
i ÝÝÝÑ σ piq ÝÝÝÑ ¨ ¨ ¨ ÝÝÝÑ σk piq “ i
§ § § §
§
πđ
§
πđ πđ
§ §
πđ
σ σ σ
j ÝÝÝÑ σ pjq ÝÝÝÑ ¨ ¨ ¨ ÝÝÝÑ σk pjq ‰ j
Hence in π always only those numbers are permuted which in σ belong to cycles of
equal length! Hoe many possibilities are there? We consider the ck cycles of length k
of σ:
` ˘
pi . . . q pi . . . q
loomoon loomoon looomooon
1 2 . . . i c k
. . . .
k k k
Let i be the smallest number in these ck k–cycles: The permutation π can map i to ck ¨ k
numbers: Thereby, at the same time, also the images
´ ¯
k´1
π piq , π pσ piqq , . . . , π σ piq
are already determined! For the next smaller number j there are thus only pck ´ 1q ¨ k
possibilities leftn, etc.: It is clear that in this manner all permutations π, which are
fixed by σ, are correctly counted; the number of these permutations thus is in total
c 1 ! ¨ 1 c1 ¨ c 2 ! ¨ 2 c2 ¨ ¨ ¨
On the other side there are exactly
n!
1 c1 ¨ c 1 ! ¨ 2 c2 ¨ c 2 ! ¨ ¨ ¨
permutations σ P Sn of cycle type pc1 , c2 , . . . q (see (2.19)). Thus we obtain:
ÿ 1 ÿ ˇ ˇ
Zpermutations “ ˇfixpermutations pσqˇ ¨ xcpσq
ně0
n!
σPSn
ÿ 1 ÿ n! ¨ pc1 ! ¨ 1c1 ¨ ¨ ¨ cn ! ¨ ncn q c1
“ ¨ x1 ¨ ¨ ¨ xncn
ně0
n! 1 ¨ c1 ! ¨ ¨ ¨ c n ! ¨ n
c 1 c n
pc1 ,...,cn q
ÿ c
ź 1
“ x11 ¨ x2c2 ¨ ¨ ¨ “ .
1 ´ xi
pc1 ,c2 ,... q iě1
The following theorems are concerned with product and composition of species
with symmetries: We actually would first need to define the symmetries for
these constructions — but how these symmetries are to be meant8, will become
clear in the proofs.
8The interpretations being dealt with are anyway the most obvious ones.
2.9. CYCLE INDEX SERIES 39
τ1 τ2
B
τl τ3
This cycle is “completed” by the maps (“relabeling of the original object B”)
τ1 , τ2 , . . . , τl to cycles of the “full” permutation σ. Here it should be observed
that
τ :“ τ1 ˝ τ2 ¨ ¨ ¨ τl (2.23)
gives a permutation which fixes B. For instance, in the following schematic
graphic the cycle length is l “ 4 and the permutation (in cycle notation) is
τ “ p12q p3q:
τ1
1 4
2 5
3 6
τ4 τ2
12 9
11 8
10 7
τ3
From the graphic it is immediately clear: If the cycle decomposition of τ is given
by
τ “ ξ1 ˝ ξ2 ¨ ¨ ¨ ξs ,
then the corresponding cycles of the “composed permutation” σ have the lengths
pℓ pξ 1 q ¨ l, ℓ pξ 2 q ¨ l, . . . q.
Further is follows: For a fixed permutation τ which fixes B, there are
p}B}!ql´1
“completions” pτ1 , τ2 , . . . , τl q. For the first pl ´ 1q τi ’s can be chosen arbitrarily;
τl is then uniquely determined by (2.23).
In total we can thus continue the calculation in the following way:
` ˘ ` ˘
ÿ 1 ÿ ÿ xl , x2l , . . . cpτ1 q ÿ xl , x2l , . . . cpτm q
ZA p B q “ 1 1
¨¨¨ m m
k! }B1 }! }Bm }!
kě0 APAk B1 Pζ 1 Bm Pζ m
ρpAq“A ℓpζ 1 q“l1 ℓpζ m q“lm
ρ“ζ 1 ˝ζ 2 ...ζ m τ1 pB1 q“B1 τm pBm q“Bm
“ ZA pZB px1 , x2 , . . . q , ZB px2 , x4 , . . . q , . . . q .
42 2. SPECIES
Exercise 16: Determine the cycle index series of the species Fixfree of permutations without fixed
points.
Hint: Show the relation Sets ¨ Fixfree “ Permutations.
Exercise 17: Determine the cycle index series for the species “set partitions”.
b
e
c
(2) P “ 2rns (Power set of rns) ordered by set inclusion: This poset is said to be
the Boolean algebra Bn . E.g., for n “ 3:
t1,2,3u
321
231 312
213 132
123
6 10 15
2 3 5
tt1u,t2u,t3uu
(6) P “ Vn,q , the family of subspaces of a finite vectorspace pGF pqqqn over a finite
field GF pqq, ordered by “subspace–relation” (that is set inclusion).
D EFINITION 3.1.3. Given two posets P, Q with order relation ďP (for P) and ďQ (for
Q). A map φ : P Ñ Q is said to be order-preserving, if, for all x, y P P, we have:
x ďP y ùñ φ pxq ďQ φ pyq .
The posets P and Q are said to be isomorphic if there exists a bijective map φ : P Ñ Q
such that φ and φ´1 both are order-preserving (such a map φ is said to be an order
isomorphism), that is
x ďP y ðñ φ pxq ďQ φ pyq .
E XAMPLE 3.1.4. Let P “ t0, 1un be the set of all t0, 1u–vectors of length n, ordered by
pv1 , . . . , vn q ď pw1 , . . . , wn q :ðñ vi ď wi for i “ 1, 2, . . . , n.
The interpretation of these vectors as characteristic functions of subsets induces an
order-preserving bijection P Ñ Bn on the Boolean algebra Bn (see example 3.1.2).
Exercise 20: Let P be some finite poset and f : P Ñ P an order–preserving bijection. Show that
f ´1 is also order–preserving.
Show that this is not true in general for infinite posets.
b
Q induziert
P e e
c
Q schwach
f f e f
The set of all order ideals of P, ordered by inclusion, is itself a poset, which is denoted
by J pPq in this text: Such a poset always has a minimum 0̂ (the empty set) and always
has a maximum 1̂ (P itself).
P ROPOSITION 3.1.9. Let P be a poset, I P J pPq be a order ideal of P, and x a maximal
element in I.
Then also Iz txu is an order ideal of P.
In case that P is finite, there is a bijection between order ideals I and antichains
A:
I ÞÑ A :“ tx : x maximal in Iu
A ÞÑ I “ xAy :“ ty : Dx P A : y ď xu
D EFINITION 3.1.10 (Linear extension). An order-preserving bijection from a finite
poset P with |P| “ n to rns is said to be a linear extension of P.
P ROPOSITION 3.1.11. The number of linear extensions of a finite poset P equals the
number of maximal chains in J pPq.
In case ρ pxq “ i, we say: x has rank i. For such a poset P, let pi be the number of
elements of rank i: Then the polynomial
n
ÿ
F pP, qq :“ p i qi
i“0
is said to be the rank generating function of P.
Exercise 22: Consider the “zigzag–poset” Zn with elements x1 , x2 , . . . , x n and cover relations
x2i´1 Ì x2i for i ě 1, 2i ď n and x2i Í x2i`1 for i ě 1, 2i ` 1 ď n
a) How many order ideals are there in Zn ?
b) Let Wn pqq be the rank generating function of the lattice of order ideals J pZn q of Zn . For instance,
W0 pqq “ 1, W1 pqq “ 1 ` q, W2 pqq “ 1 ` q ` q2 , W3 pqq “ 1 ` 2q ` q2 ` q3 . Show:
8
ÿ 1 ` p1 ` qqz ´ q2 z3
W pq, zq :“ Wn pqq zn “ .
n “0
1 ´ p1 ` q ` q2 qz2 ` q2 z4
Using Hasse diagrams, the direct sum and the ordinal sum can be illustrated as
follows:
P
P`Q P‘Q
Q
3.2. CONSTRUCTION OF POSETS 49
We obviously have
P ` Q » Q ` P,
however, in general,
P ‘ Q fi Q ‘ P.
Both, the direct sum and the ordinal sum, are associative.
An antichain with n elements is obviously isomorphic to the n-fold direct sum
n ¨ ˝ :“ ˝loooooomoooooon
` ˝ ` ¨ ¨ ¨ ˝,
n times
and a chain with n elements is isomorphic to the n-fold ordinal sum
˝loooooomoooooon
‘˝‘¨¨¨˝.
n times
Using Hasse diagrams, the direct sum and the ordinal sum can be illustrated as
follows:
P
PˆQ PbQ
Q
Again, we have
P ˆ Q » Q ˆ P,
but, in general,
P b Q fi Q b P.
The direct product is associative.
Exercise 23: Let P, Q be graded posets, let r and s be the maximal ranks of P and Q, respectively,
and let F pP, qq and F pQ, qq be the corresponding rank generating functions. Show:
a) If r “ s (otherwise maximal chains would be of different lengths), then F pP ` Q, qq “ F pP, qq `
F pQ, qq.
b) F pP ‘ Q, qq “ F pP, qq ` qr`1 F pQ, qq.
c) F pP ˆ Q, qq “ F pP, qq ¨ F pQ, qq.
` ˘
d) F pP b Q, qq “ F P, qs`1 ¨ F pQ, qq.
50 3. POSETS
ř
Exercise 25: Let P be a finite poset and define GP pq, tq :“ I q| I | tmp I q , where the summation range
is the set of all order ideals I of P, and where m pIq denotes the number of maximal elements of I.
(For instance: GP pq, 1q is the rank generating function of J pPq.)
a) Let Q be a poset with n elements. Show:
` ` ˘˘
GPbQ pq, tq “ GP qn , q´n ¨ GQ pq, tq ´ 1 ,
where P b Q denotes the ordinal product.
b) Let P be a poset with p elements. Show:
ˆ ˙
q´1
GP q, “ qp.
q
3.3. Lattices
D EFINITION 3.3.1 (Lattice). Let P be a poset, and let x, y P P: An element z with
z ě x and z ě y
is said to be an upper bound of x and y; z is said to be the least upper bound or
supremum of x and y, if, for every upper bound w of x and y, we have
w ě z.
If there exists a least upper bound z of x and y, then it is of course unique; we write
x _ y and say “x sup y”.
Analogously, we define the greatest lower bound or infimum (if it exists): We write
x ^ y and say “x inf y”.
A poset in which each pair px, yq of elements has a least upper bound as well as a
greatest lower bound is said to be a lattice.
A finite lattice has a (unique) minimum 0̂ and a (unique) maximum 1̂.
An element x in a lattice that covers 0̂ is said to be an atom.
An element x in a lattice that is covered by 1̂ is said to be a coatom.
E XAMPLE 3.3.2. Examples of finite lattices are as follows:
(1) rns with the usual linear order.
(2) the power set 2rns , ordered by set inclusion.
(3) Sn with the weak order.
3.3. LATTICES 51
P ROOF. First we show the following: The relation defined above is indeed a
partial order.
P ROOF. It suffices to show that supremum x _ y exists for all x, y P P. For this
purpose, it makes sense to define
ľ
x _ y :“ z.
zěx,y
The set over which the infimum is taken is not empty (as 1̂ ě x, y) and also
finite, which guarantees that x _ y is well-defined. It also has the property of
a supremum: Ź As a, b ě x, y ùñ x, y ď a ^ b, it can be proven (induction)
that x, y ď zěx,y z, and, for every upper bound a ě x, y, we obviously have
Ź
a ě zěx,y z.
R EMARK 3.3.8. Let P and Q be lattices. Then the direct product P ˆ Q is also a lattice,
and we have
px1 , y1 q _ px2 , y2 q “ px1 _ x2 , y1 _ y2 q ,
px1 , y1 q ^ px2 , y2 q “ px1 ^ x2 , y1 ^ y2 q .
Similarly, Q P is a lattice with
p f _ gq pxq “ f pxq _ g pxq ,
p f ^ gq pxq “ f pxq ^ g pxq .
Exercise 26: Let L be a finite lattice. Show that the following three conditions are equivalent for all
x, y P L:
(a) L is graded (i.e.: all maximal chains have in L the same length), and for the rank function rk of
L there holds
rk pxq ` rk pyq ě rk px ^ yq ` rk px _ yq .
(b) If y covers the element x ^ y, then x _ y covers the element x.
(c) If x and y both cover element x ^ y, then x _ y covers both elements x and y.
(A lattice L obeying one of these conditions is called semimodular.)
Hint: Employ an indirect proof for (c) ùñ (a): For the first assertion in (a), if there are intervals
which are not graded, then we may choose an interval ru, vs among them which is minimal with respect
to set–inclusion (i.e., every sub–interval is graded). Then there are two elements x1 , x2 P ru, vs, which
both cover u, and the length of all maximal chains in rx i , vs is ℓi , such that ℓ1 ‰ ℓ2 . Now apply (b)
or (c) to x1 , x2 .
For the second assertion in (a), take a pair x, y P L with
rk pxq ` rk pyq ă rk px ^ yq ` rk px _ yq , (3.4)
such that the length of the interval rx ^ y, x _ ys is minimal, and under all such pairs, also rk pxq `
rk pyq is minimal. Since it is impossible that both x and y cover x ^ y (why?), w.l.o.g. there is an
element x1 with x ^ y ă x1 ă x. Show that X “ x, Y “ x1 _ y is a pair such that rk pXq ` rk pYq ă
rk pX ^ Yq ` rk pX _ Yq, but where the length of the interval rX ^ Y, X _ Ys is less than the length
of rx ^ y, x _ ys.
Exercise 27: Let L be a finite semimodular lattice. Show that the following two conditions are
equivalent:
a) For all elements x, y, z P L with z P rx, ys (i.e., x ď y) there is an element u P rx, ys, such that
z ^ u “ x and z _ u “ y (u is a “complement” of z in the interval rx, ys).
b) L is atomic, i.e.: Every element can be represented as the supremum of atoms.
(A finite semimodular lattice obeying one of these conditions is called geometric.)
54 3. POSETS
Exercise 28: Let G be a (labelled) graph on n vertices. A partition of the vertex–set V pGq is called
connected if every block of the partition corresponds to a connected induced subgraph of G. The set
of all connected partitions is a subposet of the poset of partitions of V pGq, and thus a poset itself.
(If G is the complete graph, the poset of connected partitions of V pGq is the same as the poset of
all partitions of V pGq.)
Show that the poset of connected partitions of G is a geometric lattice.
Exercise 29: A lattice L is called modular if it is graded and for all x, y P L there holds:
rk pxq ` rk pyq “ rk px ^ yq ` rk px _ yq . (3.5)
(In particular, the lattice L pVq of subspaces of a finite vector space is modular.)
Show: A finite lattice L is modular if and only if for all x, y, z P L with x ď z there holds:
x _ py ^ zq “ px _ yq ^ z. (3.6)
Hint: Show that (3.6) implies the Diamond Property: The mappings
ψ : ry, y _ zs Ñ ry ^ z, zs , ψ pxq “ x ^ z
ϕ : ry ^ z, zs Ñ ry, y _ zs , ϕ pxq “ x _ y
are order preserving bijections with φ ˝ ψ “ id, see the following picture:
y_z
ϕ
y z
y^z
Exercise 30: Show: The lattice Πn of all partitions of an n–element set is not modular.
“ a _ ppa _ bq ^ cq Ðabsorption
R EMARK 3.3.11. Let L1 and L2 be two distributive lattices. Then also the direct prod-
uct L1 ˆ L2 is a distributive lattice.
E XAMPLE 3.3.12. Examples for distributive lattices are:
(1) Linear order
(2) Boolean lattice Bn : Here we have A _ B “ A Y B and A ^ B “ A X B; the
laws of distributivity are the set theoretic identities
A X pB Y Cq “ pA X Bq Y pA X Cq
A Y pB X Cq “ pA Y Bq X pA Y Cq
(3) Divisor lattice Tn “ td P N : d | nu: Indeed, let n “ p1k1 ¨ ¨ ¨ prkr be the unique
prime factor decomposition of n, then
Tn » r0, k1 s ˆ ¨ ¨ ¨ r0, kr s .
(4) For an arbitrary poset P, the lattice of order ideals J pPq is always distributive.
On the other hand, simple examples of lattices that are not distributive are the following
(given by their Hasse diagram):
1̂ 1̂
c
M5 a b c N5 b
a
0̂ 0̂
a^pb_cq‰pa^bq_pa^cq a_pb^cq‰pa_bq^pa_cq
ta,b,cu
b a b a tbu tau
T HEOREM 3.3.17. Let L be a finite distributive lattice. Then L has a rank function rk;
it is given by
rk pxq “ |I pxq| ,
where I pxq is the ideal defined in (3.7) (of the poset of supremum-irreducible elements
of L). In particular, it follows that the rank rk pLq of L is equal to the number of
supremum-irreducible elements of L.
P ROOF. By Birkhoff’s Theorem (Theorem 3.3.13), we have
L » J pPq ,
where P is the poset of supremum-irreducible elements in L. Consequently, it
suffices to show that all maximal chains in J pPq have the same length.
Let
0̂ “ H “ I0 Ì I1 Ì ¨ ¨ ¨ Ì Ir “ 1̂ “ P
be a maximal chain in P. We claim: |Ik | “ |Ik´1 | ` 1. Indeed, if S :“ Ik zIk´1 “
tx1 , . . . , xm u, then, using S ‰ H, there exists a maximal element xi (with respect
to the order of P) in S. But then Ik z txi u is also an order ideal of P (see Corol-
lary 3.1.9), and, since Ik Í Ik´1 , we can conclude Ik´1 “ Ik z txi u. Hence all
maximal chains in J pPq have the same length |P|.
AzD CzD
D
BzD
Evidently, I pPq is a vector space over C under the addition and scalar multiplication.
As for the multiplication (convolution), it satisfies the law of associativity:
ÿ
pp f ‹ gq ‹ hq px, yq “ p f ‹ gq px, zq h pz, yq
zPP
˜ ¸
ÿ ÿ
“ f px, uq g pu, zq h pz, yq
zPP uPP
ÿÿ
“ f px, uq g pu, zq h pz, yq
uPP zPP
ÿ
“ f px, uq pg ‹ hq pu, yq
uPP
“ p f ‹ pg ‹ hqq px, yq .
Furthermore, we have the laws of distributivity
f ‹ pg ` hq “ f ‹ g ` f ‹ h,
p f ` gq ‹ h “ f ‹ h ` g ‹ h
(as can be shown easily). Finally, with
#
1 : x “ y,
δ P I pPq : δ px, yq “ ,
0 otherwise
we obviously have
δ ‹ f “ f ‹ δ “ f,
that is, I pPq is a (non-commutative) algebra with unit δ: We refer to it as the inci-
dence algebra of P.
3.4. INCIDENCE ALGEBRA AND MÖBIUS INVERSION 59
E XAMPLE 3.4.2. For P “ rns, I pPq is equal to the algebra of n ˆ n upper triangular
matrices over C: Indeed, each value f px, yq of an element f P I pPq can be interpreted
as the entry in row x and columns y of a corresponding triangular matrix M, and,
conversely, each such matrix defines uniquely a function f P I pPq. It can also be seen
easily that the multiplication in I pPq corresponds to matrix multiplication.
T HEOREM 3.4.3. Let P be a locally finite poset. An element f P I pPq has an inverse
element f ´1 (i.e.: D f ´1 P I pPq with f ‹ f ´1 “ f ´1 ‹ f “ δ) if and only if f px, xq ‰
0 @x P P. If an inverse exists, then it is uniquely determined.
and, using a well-known trick (telescoping series) for the geometric series, we
have
l
ÿ
pδ ´ pζ ´ δqq ‹ pζ ´ δqk “ δ ´ loooomoooon
pζ ´ δql`1 “ δ,
k“0
“0
D EFINITION 3.4.7. The zeta function ζ of a locally finite poset P is according to The-
orem 3.4.3 invertible: µ :“ ζ ´1 is said to be the Möbiusfunction of P. By definition,
we have:
ÿ ÿ
ζ px, zq ¨µ pz, yq “
loomoon µ pz, yq “ 0 for all x ă y, (3.11)
xďzďy xďzďy
“1
ÿ ÿ
µ px, zq ¨ loomoon
ζ pz, yq “ µ px, zq “ 0 for all x ă y (3.12)
xďzďy xďzďy
“1
and µ px, xq ” 1.
T HEOREM 3.4.8 (Möbius inversion). Let P be a locally finite poset such that every
principal order ideal is finite. Let f , g : P Ñ C be functions on P with values in C.
Then the following assertions are equivalent:
ÿ
g pxq “ f pyq @x P P (3.13)
yďx
ÿ
f pxq “ g pyq µ py, xq @x P P. (3.14)
yďx
ÿ ÿ ´n¯
g pnq “ f pdq ðñ f pnq “ µ g pdq
d
d|n d|n
(the interval rd, ns in this poset is isomorphic to r1, n{ds, that is µ pd, nq “ µ p1, n{dq “:
µ pn{dq).
62 3. POSETS
P ROOF OF THE M ÖBIUS INVERSION . For the proof, we only need to perform
two simple calculations:
ÿ ÿ ÿ
g pyq ¨ µ py, xq “ f pzq ¨ µ py, xq Ðby (3.13)
yďx yďx zďy
ÿ ÿ
“ µ py, xq f pzq ¨ ζ pz, yq Ðdefinition ζ
yďx zďy
ÿ ÿ
“ f pzq ζ pz, yq ¨ µ py, xq
zPP zďyďx
ÿ
“ f pzq δ pz, xq
zPP
“ f pxq ,
and, conversely:
ÿ ÿ ÿ
f pxq “ g pyq ¨ µ py, xq Ðby (3.14)
xďs xďs yďx
ÿ ÿ
“ ζ px, sq g pyq ¨ µ py, xq Ðdefinition ζ
xďs yďx
ÿ ÿ
“ g pyq µ py, xq ¨ ζ px, sq
yPP yďxďs
ÿ
“ g pyq δ py, sq
yPP
“ g psq .
R EMARK 3.4.10. From an abstract point of view, the key observation is that the equa-
tions (3.13) and (3.14) are equivalent to the “symmetric” equations
ÿ
g pxq “ f pyq ¨ ζ py, xq @x P P,
yPP
ÿ
f pxq “ g pyq ¨ µ py, xq @x P P,
yPP
and that I pPq acts on the vector space C P of all functions P Ñ C as algebra of linear
mappings (from the right) with
ÿ
p f ξq pxq :“ f pyq ξ py, xq für f P C P , ξ P I pPq .
yPP
E XAMPLE 3.4.12. Let S1 , S2 , . . . , Sn be sets. Consider the poset that consists of all in-
tersections of selections of these sets, ordered by inclusion. Here the empty intersection
is defined as č
Si :“ S1 Y S2 Y ¨ ¨ ¨ Y Sn “ 1̂.
iPH
We aim at deriving a formula for |S1 Y S2 Y ¨ ¨ ¨ Y Sn | and
Ť define for this purpose func-
tions g, f : P Ñ C as g pTq :“ |T| and f pTq :“ |Tz p T1 ăT T 1 q| (in words: f pTq is
the number of elements in T that are not contained in any T 1 P P with T 1 Ř T).
First we observe that in general we have the following (think for a moment!):
ÿ ` ˘
g pTq “ f T1 .
T 1 ĎT
Using Möbius inversion, we obtain
ÿ ` ˘ ` ˘
f pTq “ g T1 ¨ µ T1, T .
T 1 ĎT
In particular, we have
` ˘ ÿ ` ˘ ` ˘ ÿ ` ˘ˇ ˇ
0 “ f 1̂ “ g T 1 ¨ µ T 1 , 1̂ “ µ T 1 , 1̂ ˇT 1 ˇ ,
T 1 Ď1̂ T 1 Ď1̂
we leads us to the following formulation of the inclusion-exclusion principle:
ˇ ˇ ÿ ` ˘ˇ ˇ
|S1 Y S2 Y ¨ ¨ ¨ Y Sn | “ ˇ1̂ˇ “ ´ µ T 1 , 1̂ ˇT 1 ˇ .
T 1 Ř1̂
In the introductory example (see the beginning of Abschnitt 3.4) we then have
the following:
` ˘
1̂ µ 1̂, 1̂ “ 1
` ˘
A B C µ ., 1̂ “ ´1
` ˘
0̂ µ 0̂, 1̂ “ 2
The Möbius inversion implies the “discrete analogue of the fundamental theorem of
calculus”
#
ÿn
g p1q for n “ 1,
g pnq “ f piq @n ě 1 ðñ f pnq “
i“1
g pnq ´ g pn ´ 1q for n ą 1.
From an abstract point of view, this means that the summation operator
n
ÿ
pΣ cq pnq :“ c piq
i“1
and the difference operator
#
c p1q n“1
p∆ cq pnq :“
c pnq ´ c pn ´ 1q ną1
are inverse operators on the vector space CN of all functions c : N Ñ C.
T HEOREM 3.4.14. Let P and Q be locally finite posets. For the direct product P ˆ Q,
we have
µPˆQ ppx1 , y1 q , px2 , y2 qq “ µP px1 , y1 q ¨ µQ px2 , y2 q .
P ROOF. Recall that in the direct product P ˆ Q we have
px1 , y1 q ď px2 , y2 q :ðñ x1 ď x2 and y1 ď y2 .
We check whether the given function has the defining property of the Möbius
function:
ÿ
µP px3 , x2 q ¨ µQ py3 , y2 q “
px1 ,y1 qďpx3 ,y3 qďpx2 ,y2 q
ÿ ÿ
µP px3 , x2 q ¨ µQ py3 , y2 q “ δP px1 , x2 q ¨ δQ py1 , y2 q
x1 ďx3 ďx2 y1 ďy3 ďy2
“ δPˆQ ppx1 , y1 q , px2 , y2 qq .
E XAMPLE 3.4.15 (Classical Möbius function). The divisor lattice Tn :“ td : d | nu
(order relation “ divisor relation) for n with prime factor decomposition n “ p1k1 ¨ p2k2 ¨
pkmm is isomorphic to the product of linear posets
Tn » r0, k1 s ˆ ¨ ¨ ¨ ˆ r0, km s .
By Theorem 3.4.14, we obtain:
´ ¯
β β
µTn px, yq “ µTn p1α1 ¨ ¨ ¨ pαmm , p1 1 ¨ ¨ ¨ pmm
“ µr0,k1 s pα1 , β 1 q ¨ ¨ ¨ µr0,km s pαm , β m q
# řm
p´1q i“1 pβi ´αi q if pβ i ´ αi q P t0, 1u @i,
“
0 otherwise.
Since we have the following isomorphy
”yı
rx, ys “ 1,
x
3.4. INCIDENCE ALGEBRA AND MÖBIUS INVERSION 65
E XAMPLE 3.4.17 (Boolean lattice). The power set 2rns of rns can be interpreted as
family of characteristic functions (see Example 3.1.4):
A Ď rns Ø f A : rns Ñ t0, 1u , f A pxq “ rx P As 1.
1Iversons notation.
66 3. POSETS
C OROLLARY 3.4.22. Let L be a finite lattice with at least 2 elements, and let 1̂ ‰ a P L.
Then we have ÿ ` ˘
µ x, 1̂ “ 0. (3.20)
x:x^a“0̂
It follows
` that we ˘can compute the Möbius function of Πn in general, if we can deter-
mine µ 0̂Πm , 1̂Πm for all m P N. We use Corollary 3.4.22. Let a “ 1 2 . . . pm ´ 1q |m.
If x ^ a “ 0̂Πm , then we either have x “ 0̂Πm or x “ m i| premainder in blocks with 1 elementq
for an i “ 1, 2, . . . , m ´ 1. According to (3.20), we have
` ˘ ÿ ` ˘
µm :“ µ 0̂Πm , 1̂Πm “ ´ µ x, 1̂Πm “ ´ pm ´ 1q µm´1 .
x“m i|premainderq
On the other hand, we can also simply expand the product in order to obtain
an expansion in terms of the basis L:
ÿ ÿ ÿ ` ˘
1̂ L ´ x` x^y´`¨¨¨ “ µ x, 1̂ L x .
xPX x,yPX looooooomooooooon
xPL
δ1̂
L
(See also the paper of Rota and Stanley [6] regarding Möbiusfunctions.)
E XAMPLE 3.4.28 (Möbiusfunction of Vn pqq “ GF pqqn ). Let V,”!W)ď Vn pqqı be two
subspaces with V ď W. Then the interval rV, Ws is isomorphic to ~0 , W{V , where
´! ) dim pV{Wq
the latter interval only depends on the dimension ¯ “ dim pVq ´ dim pWq.
Therefore, it suffices to compute µn :“ µ ~0 , Vn pqq . For this purpose, we use the
dualised version of Corollary 3.4.22:
ÿ ` ˘
µ 0̂, x “ 0. (3.22)
x:x_a“1̂
µn “ p´1qn qp 2q .
n
As a small application, we compute the number of spanning sets in Vn pqq. Let f pWq
be the number of subsets that span W, and let g pWq be the number of non-empty
subspaces of W. Evidently, we have:
ÿ
g pWq “ f pTq .
TďW
in particular
ÿ n „
ÿ ´ ¯
n
p´1qn´k qp q 2q k ´ 1 .
n´k
f pVn pqqq “ µ pT, Vn pqqq g pTq “ 2
k q
T k“0
Exercise 31: Prove the “NBC–Theorem” (“Non–broken circuit theorem”) of G.–C. Rota: Let L be
geometric lattice. We assume that the atoms Žof L are labelled (with natural numbers 1, 2, . . . ). A
set B of atoms is called independent, if rk p Bq “ |B|, otherwise it is called dependent. A set C
of atoms is called a circuit, if C is a minimal dependent set. A broken circuit is a set corresponding
to a circuit from which its largest atom (with respect to the labeling of atoms) was removed. A
non–broken circuit is a set B of atoms which does not contain a broken circuit. Then Rota’s Theorem
states: ´ ¯
` ˘ ł
µ 0̂, x “ p´1qrkp xq ¨ # non–broken circuits B with B“x .
Exercise 32: Show: The Möbius function µ px, yq of a semimodular lattice is alternating, i.e.
Moreover, show that the Möbius function of a geometric lattice is strictly alternating, i.e.
p´1qlength of r x,ys µ px, yq ą 0.
Hint: Use the following formula for the Möbius function of a lattice:
ÿ ` ˘
µ 0̂, x “ 0 for all a P L. (3.23)
x:x _ a“1̂
Show that if a is an atom, then there follows:
` ˘ ÿ ` ˘
µ 0̂, 1̂ “ ´ µ 0̂, x . (3.24)
x coatom
xğa
CHAPTER 4
Asymptotic enumeration
Now we may ask: What is the “order of magnitude” of the number Bn , for
n large? The analogous question “what is the order of magnitude of n!?” is
answered by Stirling’s Formula
? ´ n ¯n
n! „ 2πn for n Ñ 8. (4.1)
e
We want to find similar “asymptotic” results also in other cases.
z2 “ o pzq pz Ñ 0q
1
“ o p1q pn Ñ 8q
n ´ n ¯n
?
n! “ 2πn p1 ` o p1qq pn Ñ 8q
e
o p f pzq ¨ g pzqq “ f pzq ¨ o pg pzqq pz Ñ ζq
The meaning of the last “equation” is (cf. Remark 4.1.2): The class of functions φ, for
which
φ pzq
lim “0
zÑζ f pzq ¨ g pzq
holds, is contained in the class of functions f ¨ ψ, for which
ψ pzq
lim “0
zÑζ g pzq
φpzq
holds (simply set ψ pzq :“ f pzq ).
4.1. LANDAU’S NOTATION 75
Exercise 33: Let f , g : N Ñ R` be two functions (from the natural numbers to the nonnegative
reals). Which of the following rules is valid for n Ñ 8 (under which preconditions)?
O p f pnqq ` O pg pnqq “ O p f pnq ` g pnqq O p f pnqq ´ O pg pnqq “ O p f pnq ´ g pnqq
O p f pnqq
O p f pnqq ¨ O pg pnqq “ O p f pnq ¨ g pnqq “ O p f pnq {g pnqq
O pg pnqq
´ ¯
O p f pnqqOpgpnqq “ O f pnq gpnq exp pO p f pnqqq “ O pexp p f pnqqq
b ˆb ˙
O p f pnqq “ O f pnq g pO p f pnqqq “ O pg f pnqq
Exercise 34: Same question as in the preceding exercise, where O p.q is replaced by o p.q.
Exercise 35: Let f 1 , f 2 , g1 , g2 be functions N Ñ C, such that f 1 pnq „ f 2 pnq and g1 pnq „ g2 pnq
for n Ñ 8. Which of the following rules are valid for n Ñ 8 (under which preconditions)?
f 1 pnq ` g1 pnq „ f 2 pnq ` g2 pnq f 1 pnq ´ g1 pnq „ f 2 pnq ´ g2 pnq
f 1 pnq f pnq
f 1 pnq ¨ g1 pnq „ f 2 pnq ¨ g2 pnq „ 2
g1 pnq g2 pnq
f 1 pnq g 1 pnq „ f 2 pnq g2 pnq exp p f 1 pnqq „ exp p f 2 pnqqq
b b
f 1 pnq „ f 2 pnq g1 p f 1 pnqq „ g2 p f 2 pnqq
log p f 1 pnqq „ log p f 2 pnqqq .
E XAMPLE 4.1.7. Here are some typical generating functions that we want to examine
in the following:
ř n
‚ Derangements: Dn zn! “ e´z ¨ 1´z 1
,
ř ?
n 1´ 1´4z
‚ Catalan numbers: Cn z “ 2z ,
ř zn e z ´1
‚ Bell–numbers: Bn n! “ e .
76 4. ASYMPTOTIC ENUMERATION
The asymptotics for the Catalan numbers shows the typical feature: If a gen-
erating function f has singularities, let ρ be a singularity of minimal absolute
value. Then the coefficients of f show the following asymptotic behaviour:
ˆ ˙n
n 1
Jz K f pzq „ ¨ Θ pnq
|ρ|
a
with limnÑ8 n Θ pnq “ 1; the type of the singularity determines Θ pnq. This
kind of asymptotic analysis is called singularity analysis (cf. Abschnitt 4.3). If
there are no (or no “simple”) singularities, then the saddle point method proves to
be useful in many cases (cf. Abschnitt 4.4).
We will quickly recapitulate some basic concepts from complex analysis, which
we shall need in the following:
D EFINITION 4.2.1. A domain G is an open connected subset of the complex num-
bers C: G Ď C. A function f : G Ñ C is called analytic at z0 P G if there is a
neighbourhood U Ď G of z0 in which the series expansion
ÿ
f pzq “ cn pz ´ z0 qn
ně0
where the sum runs over all poles s of h, and ns denotes the winding number1 of Γ
with respect to s.
ř
C OROLLARY 4.2.9 (Cauchy’s integral formula). Let f pzq “ ně0 f n zn be analytic
in a disk centered at 0, and let Γ be a contour in the interior of this disk, which winds
around 0 exactly once (in positive orientation) (i.e., the winding number with respect
to 0 is 1), then we have
ż
n 1 f pzq
Jz K f pzq “ f n “ dz. (4.5)
2π i Γ zn`1
ř
D EFINITION 4.2.10. Let f pzq “ ně0 cn pz ´ z0 qn be an analytic function with ra-
dius of convergence R; let ζ be a point in the boundary
BDpz0 , Rq “ tz P C : |z ´ z0 | “ Ru
of the disk of convergence. If there is a neighbourhood U of ζ and an analytic function
fˆ on U, such that ˇ
ˆf ˇˇ “ f |UXDpz0 ,Rq ,
UXDpz0 ,Rq
Exercise 36: Let f , g be complex functions which are analytic on some given domain. Which of the
following rules are valid (under which preconditions)?
Sing p f ˘ gq Ď Sing p f q Y Sing pgq Sing p f ¨ gq Ď Sing p f q Y Sing pgq
Sing p f {gq Ď Sing p f q Y Sing pgq Y Null pgq Sing p f ˝ gq Ď Sing pgq Y gp´1q pSing p f qq
´a ¯
Sing f Ď Sing p f q Y Null p f q Sing plog f q Ď Sing p f q Y Null p f q
´ ¯ ` ` ˘˘
Sing f p´1q Ď f pSing p f qq Y f Null f 1
Here, Sing p f q denotes the set of singular points of f , and Null p f q denotes the set of zeroes of f .
1Loosely speaking, the winding number enumerates how often the contour Γ “winds
around” s: In the examples we consider here, the winding number will always be equal to
1.
4.2. RECAPITULATION: ELEMENTS OF COMPLEX ANALYSIS 79
ř
T HEOREM 4.2.12. Let f pzq “ f n ¨ zn be analytic at 0 with radius of convergence
R ă 8. Then there exists a singularity ζ in the boundary of the disk of convergence
(i.e., |ζ| “ R).
We may parametrize the circle with center 0 and radius R0 as follows: z pθq “
R0 ¨ eiθ and dθ
d
z pθq “ iz pθq. Then we obtain for this contour:
ż 2π ` ˘
1 f R 0 ¨ ei θ
fn “ dθ.
2π 0 R0n ¨ en¨iθ
1
| fn| ď ¨ max | f pzq| , (4.6)
R0n |z|“R0
` ˘ ř
i.e., f n “ O R´n 0 pn Ñ 8q. Let q with R ă q ă R 0 , then f pzq “ f n zn
2
converges for |z| ă q according to the well–known root test , a contradiction
(since q ą R).
ř
T HEOREM 4.2.13 (Pringsheim). Let f pzq “ f n zn be analytic at 0 with radius of
convergence R. Assume that f n ě 0 for all n ě 0. Then the (real) number R is a
singularity of f .
2From | f | C
a q ?
n
n ď R0n we get for |z| ă q ă R0 : n
| f n zn | ă R0 C, in particular
a q
lim supnÑ8 n | f n zn | ă 1 (since R0 ă 1).
80 4. ASYMPTOTIC ENUMERATION
2h r
z0
Equivalently, we have:
ÿ 2r
f pz0 ` zq “ gm ¨ zm for |z| ď 2h “
m
3
and therefore
ˆ ˙m
m 1 d
gm “ Jz K f pz ` z0 q “ f pz ` z0 q |z“0
m! dz
1 pmq
“ f pz0 q
m!
1 ÿ
“ f n ¨ n ¨ pn ´ 1q ¨ ¨ ¨ pn ´ m ` 1q z0n´m
m! n
ÿ ˆn˙
“ f n ¨ z0n´m .
něm
m
So we have ˜ ¸
ÿ ÿ ˆn˙
pz0 ` 2hq “
f looomooon f n ¨ z0n´m p2hqm .
m něm
m
“R`h
The right–hand side is a converging double sum (since R ` h lies in the inte-
rior of DpR, rq) and has only non–negative summands: Therefore it is absolutely
convergent, and we may reorder the summands arbitrarily.
ÿ ÿ n ˆ ˙
n n´m
f pR ` hq “ fn z0 p2hqm
n m“0
m
ÿ ÿ
“ f n pz0 ` 2hqn “ f n pR ` hqn ,
n n
a contradiction, since R ` h ą R.
Stated otherwise:
an “ K n ¨ Θ pnq ,
where Θ is subexponential, i.e.
a
lim sup n
|Θ pnq| “ 1.
nÑ8
Typical examples for subexponential sequences are
? 2 π
1, n3 , n´3{2 , log n, p´1qn , e n
, elog n , n2 sin n , . . .
2
82 4. ASYMPTOTIC ENUMERATION
ř
T HEOREM 4.3.4. Let f pzq “ n f n ¨ zn be analytic at 0, let R be the absolute value of
a dominant singularity. Then we have:
ˆ ˙n
1
fn ’ .
R
P ROOF. This follows immediately from řCauchy–Hadamard’s formula for the ra-
dius of convergence R of a power series n f n zn :
1
R“ a . (4.10)
lim supnÑ8 n
| fn|
(For the absolute value of a dominant singularity clearly equals the radius of
convergence.)
Exercise 37: Let p pnq be the number of (integer) partitions of n. We know that
8
ÿ 8
ź 1
p pnq zn “ .
n “0 i “1
1 ´ zi
What are the (dominant) singular points of this generating function? What does this imply for the
asymptotic behaviour of p pnq for n Ñ 8?
R EMARK 4.3.8. The logarithmus is a covering map which is not analytic on any
domain containing 0, so 0 is not an isolated singularity of the logarithm (and thus of
?
the root function n z “ z1{n “ exp n1 log z).
P ROPOSITION 4.3.9. Let f pzq be analytic on |z| ă R with 0 ď R ă 8. Then we have
for arbitrary r P p0, Rq:
max|z|“r | f pzq|
| fn| ď ,
rn
i.e.
max|z|“r | f pzq|
| f n | ď inf .
0ărăR rn
If f pzq only has non–negative coefficients, then we have
f prq f prq
fn ď n
ùñ f n ď inf .
r 0ărăR r n
T HEOREM 4.3.14. Let f pzq be meromorphic for |z| ď R with poles at a1 , . . . , am ; let
e1 , . . . , em be the orders of these poles. Let f be analytic at 0 and for all z with |z| “ R.
Then we have:
ÿm
` ´n ˘
fn “ a´n
i polynomial i pnq `O R
loooooooomoooooooon
i“1
of degree ei ´1
Here Sa1 pzq is the “singular part” at a1 , while Ha1 pzq is analytic at a1 . This
“subtraction of singular parts” can, of course, be repeated for all the remaining
poles, whence we obtain
˜ ¸
ÿm ÿm
f pzq “ Sai pzq ` f pzq ´ Sai pzq ,
i“1 i“1
loooomoooon looooooooooomooooooooooon
“:Spzq “:Hpzq
E XAMPLE 4.3.16. The (exponential) generating function of cycles‹ (i.e., sequences pcyclesq)
is
1
f pzq “ 1
.
1 ´ log 1´z
Obviously there is one pole at z “ 1 and another at
1 1 e´1 1
log “ 1 ðñ “ e ðñ z “ “ 1 ´ » 0.632121,
1´z 1´z e e
hence the dominant singularity is ζ “ 1 ´ e´1 . Again, by de L’Hospital’s rule, we get
z´ζ 1
lim 1
“ lim ´1 “ ζ ´ 1 “ ´e´1 ,
zÑζ 1 ´ log 1´z zÑζ
1´z
i.e.,
e´1
f pzq “ ´ ` ˘ ` H pzq .
z ´ 1 ´ e´1
H pzq is analytic in the (closed!) disk with radius R “ 1 ´ ǫ (ǫ ą 0), therefore we
obtain the asympotics
ˆ ˙n`1
1 1 ` ´n ˘ 1
fn “ ` O p1 ´ ǫq » 1.58198n`1 ` O p1.0001n q .
e 1´e ´1 e
E XAMPLE 4.3.17. The number Cn,k of surjections rns Ñ rks can be expressed in terms
of the Stirling–numbers of the second kind:
Cn,k “ k! ¨ Sn,k .
We want to investigate the following question: If we assume uniform probability
distribution on the set of all surjections rns Ñ rks (n fixed, k arbitray), what is the ex-
pected value of the cardinality of the surjection’s image? Our approach is to consider
the generating function in two variables
ÿ ÿ zn k ÿ k z 1
f pt, zq :“ Cn,k ¨ ¨t “ t pe ´ 1qk “ .
n! 1 ´ t pez ´ 1q
kě0 ně0 kě0
Then we have ˜ ¸
ˇ n
d f ˇˇ ÿ ÿ zn ez ´ 1
“ k ¨ C ¨ “ .
dt ˇt“1 ně0 n,k
n! p2 ´ e z q2
k“0
The poles (order: 2) of this function are log 2 ` k ¨ 2π i, k P Z, and the Laurent–series
expansion is
ez ´ 1 1 1
“ ` ` H pzq ,
p2 ´ ez q2
4 pz ´ log 2q 2 4 pz ´ log 2q
´ ¯
where H pzq “ ´ 487 1
` 48 pz ´ log 2q ` O pz ´ log 2q2 is analytic in the (closed!)
b
disk of radius R “ 4π 2 ` plog 2q2 ´ ǫ (ǫ ą 0; we may choose R » 6.3213). From
4.3. SINGULARITY ANALYSIS 87
For the second term on the right–hand side we have according to (4.13)
O pr pnqq
“ O pr pnqq ,
g pnq ` O pr pnqq
1
since also |g pnq ` O pr pnqq| ą ǫ{2 for n “large enough”, so gpnq`Oprpnqq “ O p1q.
Of course we have O pr pnqq ` O pr pnqq “ O pr pnqq, so also (4.16) is proved.
4.3.4. Several dominant singularities.
1
E XAMPLE 4.3.19 (Cancellation). Consider the sum of rational functions 1`z2
and
1
1´z3
:
1 1 2 ` z2 ´ z3 2 ´ z2 ` z3 ` z4 ` z8 ` z9 ´ z10
` “ “ .
1 ` z2 1 ´ z3 p1 ` z2 q p1 ´ z3 q 1 ´ z12
The coefficients of zn vanish for n ” 1, 5, 6, 7, 11 pmod 12q.
E XAMPLE 4.3.20 (Nonperiodic fluctuations). Consider the rational function
1 6 11 84 3 779 4 2574 5
6
“ 1 ` z ` z2 ´ z ´ z ´ z `¨¨¨
1´ 5z ` z
2 5 25 125 625 3125
with partial fraction decomposition
1 3i 1 3i
1 2` 8 ¯ 2´ 8 ¯
“ ´ ` ´ .
1 ´ 56 z ` z2 1 ´ 53 ´ 45i z 1 ´ 53 ` 45i z
The roots of the denominators have absolute value 1 and are inverses of one another, i.e.
3 4i 3 4i
` “ eiθ and ´ “ e´iθ .
5 5 5 5
So the coefficient of zn may be written as
einθ ` e´inθ 3 einθ ´ e´inθ 3
` ¨ “ cos pn ¨ θq ` sin pn ¨ θq
2 4 2i 4
cos θ
“ cos pn ¨ θq ` sin pn ¨ θq
sin θ
sin ppn ` 1q ¨ θq
“ .
sin θ
4.3.5. “Standardized” singularity analysis. We shall now present the “(quasi–
) mechanical algorithm” of (standard) singularity analysis for determining the
asymptotics of the coefficients of some generating function f pzq, as developped
by Philippe Flajolet and Andrew Odlyzko. It consists of the following steps:
Normalization of the singularity: If ρ is the “relevant singularity” of f pzq,
consider fˆ pzq :“ f pρ ¨ zq: then 1 is the “relevant singularity” of fˆ pzq.
Asymptotics for “standard functions”: For “standard functions” of the
form ˆ ˙β
´α 1 1
p1 ´ zq log
z 1´z
there is a general result describing the asymptotics of the coefficients.
90 4. ASYMPTOTIC ENUMERATION
p
The following formula gives evaluations of the Ψ–function at rational numbers z “ q ;
cf. [2, section 1.7.3, equation (29)]:
ˆ ˙ ˆ ˙
p π πp
Ψ “ ´γ ´ log q ´ cot
q 2 q
q
21
ÿ ˆ ˙ ˆ ˆ ˙˙
pn n
` cos 2π log 2 ´ 2 cos 2π , (4.24)
q q
n“1
q
where the prime ( 2 1)
in the upper bound of the summation means that for evev q the
q
last summand (with index n “ 2 ) is multiplied by 21 . (Here, γ denotes the Euler–
Mascheroni constant.) For example:
ˆ ˙ ´π ¯ 1
1 π
Ψ “ ´γ ´ log 2 ´ cot ` cos π log p2 ´ 2 cos pπqq
2 2 2 2
“ ´γ ´ 2 log 2, (4.25)
`π˘
since cot 2 “ 0, cos π “ ´1 and log 4 “ 2 log 2.
R EMARK 4.3.23. For the following Theorem it is useful is to recall the polar coordi-
nates for complex numbers: Let z P C with ℜ pzq “ x and ℑ pzq “ y, i.e., z “ x ` iy.
Then we may express z in polar coordinates pr, φq, i.e. z “ r ¨ pcos φ ` i sin φq with
r “ |z| and φ “ φ0 ` 2πk for k P Z, where φ0 denotes the main value (in the interval
´π ă arg pzq ď π)of the argument arg pzq. (for r “ 0 we may choose the angle φ
arbitrarily.) We have:
ex`iy “ ex ¨ pcos y ` i sin yq , (4.26)
log px ` iyq “ log |x ` iy| ` i arg px ` iyq . (4.27)
In particular, this means that the (complex) logarithmus (the inverse function of the ex-
ponential function) is a multi–valued function, which might be “forced” to be single–
valued by constraining its range (for instance, by considering only the main values of
the argument).
Moreover, we have the following connection of the trigonometricn functions with the
exponential function:
eiθ ` e´iθ
cos θ “ ,
2
eiθ ´ e´iθ
sin θ “ .
2i
T HEOREM 4.3.24 (Presentation with Hankel–contour integral). For the Gamma–
function we have ż
1 1
“ p´tq´s ¨ e´t dt, (4.28)
Γ psq 2π i H
where H “ H´ ` H˝ ` H` is a contour which cosists of the following parts
‚ H´ “ tw ´ i : w ě 0u, (
‚ H˝ “ ´eiφ : ´ π2 ď φ ď π2 ,
‚ H` “ tw ` i : w ě 0u
See the following graphical illustration:
92 4. ASYMPTOTIC ENUMERATION
D
ρ
D1 ρ
In the first (straight) part of this new contour D1 we have arg p´tq “ ´π, i.e.,
p´tqz´1 “ epz´1qplog t´iπq ,
and in the last (straight) part we have arg p´tq “ π, i.e.,
p´tqz´1 “ epz´1qplog t`iπq .
(For both parts we have log t P R.) Parametrizing the small circle as ´t pθq “
δeiθ we obtain:
ż żδ żρ
´iπpz´1q z´1 ´t
p´tq z´1 ´t
e dt “ e t e ` eiπpz´1q tz´1 e´t
D ρ δ
żπ ´ ¯z´1 iθ
`i δ ei θ eδe δeiθ dθ “
´π
żρ żπ
iθ
´ 2i sin pπzq t e dt ` iδ
z´1 ´t z
eizθ`δe dθ.
δ ´π
This equality holds for all positive real numbers δ ď ρ: for δ Ñ 0 we have
şπ iθ şπ
δz Ñ 0 (since ℜ pzq ą 0) and ´π eizθ`δe dθ Ñ ´π eizθ dθ (since the integrand
converges uniformly); so we obtain:
ż żρ
p´tq e dt “ ´2i sin pπzq tz´1 e´t dt.
z´1 ´t
D 0
This holds for all positive real numbers ρ: H is obviously the “limiting case of
the contour” D for ρ Ñ 8, so there holds
ż ż8
z´1 ´t
p´tq e dt “ ´2i sin pπzq tz´1 e´t dt,
H 0
i.e., ż
1
Γ pzq “ ´ p´tqz´1 e´t dt. (4.30)
2i sin pπzq H
The contour H is not close to 0, therefore we do not need the assumption ℜ pzq ą
0 any more4: The integral in (4.30) is a (single–valued) analytic function, which
equals Γ pzq for ℜ pzq ą 0; so Hankel’s formula (4.30) for the Gamma–function
holds for all z P CzZ. Replacing z by 1 ´ z and applying Euler’s Reflection
Formula (4.20), we obtain the assertion.
4.3.5.2. Asymptotics for “standard functions”. According to the binomial the-
orem, for the “standard function” p1 ´ zq´α (where α P C is arbitrary) we obtain
ˆ ˙ ˆ ˙
n ´α n ´α n`α´1
Jz K p1 ´ zq “ p´1q “ .
n n
If α P t0, ´1, ´2, . . . u, then the asymptotics of these coefficients is trivial (since
in that case p1 ´ zq´α is a polynomial, and the coefficients are all 0 for n ą ´α).
Otherwise, we can write
ˆ ˙
n`α´1 α ¨ pα ` 1q ¨ ¨ ¨ pα ` n ´ 1q Γ pn ` αq
“ “ .
n n! Γ pαq Γ pn ` 1q
Hence, we could determine the asymptotics by means of Stirling’s formula for
the gamma function (4.23)
ˆ ˆ ˙˙
? 1
Γ ps ` 1q “ s e
s ´s
2πs 1 ` O .
s
However, it is more convenient to apply Cauchy’s integral formula (4.5):
T HEOREM 4.3.25 (Asymptotics for coefficients of standard functions, 1). Let α P
Cz t0, ´1, ´2, . . . u. We consider the “standard function”
f pzq “ p1 ´ zq´α .
Then, asymptotically as n Ñ 8, we have
˜ ¸
ÿ8
nα´1 ek pαq
Jzn K f pzq „ 1` , (4.31)
Γ pαq nk
k“1
4An analytic continuation for the logarithm exists only in the slitted complex plane!
94 4. ASYMPTOTIC ENUMERATION
Kurve CR,n
R
1{n
0 1 1{n
For all n ą ℜ p´αq, the integral becomes arbitrarily small along the “circle
piece” of C R,n if we let R Ñ 85, so that “only” the “Hankel-like piece” Hn “
H´ ˝ `
n ` H n ` H n of C R,n plays a role, which consists of the subpieces
i (
‚ H´n “!w ´ n : w ě 1 , )
iφ
‚ H˝n “ 1 ´ en : ´ π2 ď φ ď π2 ,
‚ H` “ w ` i : w ě 1(.
n n
Performing the substitution,
t 1
z “ 1 ` , dz “ dt
n n
we obtain ż ˆ ˙
n nα´1
´α ´α t ´n´1
Jz K p1 ´ zq “ p´tq 1` dt, (4.34)
2i π H n
where H is the Hankel contour from Theorem 4.3.24.
ˇ ˇ
5The reason is that the modulus of the integrand can be estimated from above as ˇˇ p1´zq ˇ
´α
ˇď
ˇ z n`1 ˇ
Rℜp´αq´n´1 , and hence the integral is less than or equal to 2πRℜp´αq´n .
4.3. SINGULARITY ANALYSIS 95
` t
˘
For |t| ă n the power series expansion of log 1 ` n converges, hence we ob-
tain
ˆ ˙
t ´n´1
1` “ e´pn`1q logp1`t{nq
n
ˆ ˆ ˙˙
t t2 t3
“ exp ´ pn ` 1q ´ ` ´`¨¨¨
n 2n2 3n3
ˆ ˙
t2 ´ 2t 3t2 ´ 2t3
“ exp ´t ` ` ` ¨ ¨ ¨ Ðsort n´k
2n 6n2
ˆ 2 ˙
t ´ 2t 3t2 ´ 2t3
“ e ¨ exp
´t
` `¨¨¨
2n 6n2
ˆ ˙
t2 ´ 2t 3t4 ´ 20t3 ` 24t2
“e ´t
1` ` `¨¨¨ . (4.35)
2n 24n2
That is, for |t| ď R (R ă 8 arbitrary but fixed) the integrand of (4.34) converges
for n Ñ 8 uniformly to p´tq´α e´t . “Asymptotically”, we therefore have
ˆ ˙ ˆ ˆ ˙˙
t ´n´1 1
1` “e ´t
1`O for n Ñ 8.
n n
“Formally”, by substitution in (4.34), we obtain
ż ˆ ˆ ˙˙
´α nα´1 ´α ´t 1
n
Jz K p1 ´ zq “ p´tq e 1`O dt
2i π H n
ˆ ˆ ˙˙
nα´1 1
“ 1`O .
Γ pαq n
That was admittedly somewhat bold, and it needs justification. To this end, we
split the Hankel contour H
‚ into the part with ℜ ptq ě log2 n: we shall show that this part is “negli-
gible”;
‚ and into the part with ℜ ptq ă log2 n: inbthat case, for n large enough
2
log n ?
we have |t|n „ n ă 1 (because |t| ď ℜ ptq 2
` 1 and log z ă z for
all z ą 0) and the expansion (4.35) “works”.
To see that the first part is negligible, some finer work is required for the esti-
mation of the integral. The subparts of the integrals over H` and H´ can be
dealt with “in once” because, to start with, we have
ˇż ˇ ż
ˇ 8 p´t ˘ i q α ˇ 8
|t ` i|α
ˇ ˇ
ˇ dt ˇ ď ˇ ˇn`1 dt.
ˇ log2 n `1 ` t˘i ˘n`1 ˇ log2 n ˇ1 ` t`i ˇ
n n
From
zα “ eα log z “ epℜpαq`iℑpαqqplog|z|`i arg zq
“ eℜpαq log|z|´ℑpαq arg z ¨ eipℑpαq log|z|`ℜpαq arg zq
it follows that
|zα | “ eℜpαq log|z|´ℑpαq arg z ,
96 4. ASYMPTOTIC ENUMERATION
Furthermore, we have ˇ ˇ ˆ ˙
ˇ1 ` t ` i ˇ ě 1 ` t ,
ˇ ˇ
ˇ n ˇ n
hence we may continue the estimation of the integral (we put α1 “ ℜ pαq):
ˇż ˇ ż8
ˇ 8 p´t ˘ iqα ˇ tα
1
ˇ ˇ
ˇ dtˇ ď C ¨ ˇ ˇn`1 dt.
ˇ log2 n `1 ` t˘i ˘n`1 ˇ log2 n ˇ1 ` t ˇ
n n
Now we claim: for fixed k, we may choose ann such that for t ě log2 n we have
1
tα 1
` ˘ ď . (4.37)
t n`1 n k t2
1` n
For, if we take the logarithm on both sides of this claim, then we obtain the
equivalent inequality
ˆ ˙
` 1 ˘ t
α ` 2 log t ´ pn ` 1q log 1 ` ď ´k log n (4.38)
n
(since the real logarithm is a monotone increasing function), and if we substi-
tute for t Ñ log2 n, then it is obvious that the following inequality holds for n
large: ˜ ¸
` 1 ˘ 2
log n
α ` 2 log log2 n ´ pn ` 1q log 1 ` ď ´k log n
n
looooooooooooooomooooooooooooooon
ěp1´ǫq log2 n by (4.36)
(since log2 n grows faster than log log2 n and log n). If we consider the deriva-
tive of (4.38), then we see that
2 ` α1 n ` 1 1 2 ` α1 n ` 1
´ “ ´ ď 0,
t n 1 ` t{n t n`t
np2`α1 q
for t ě n´α1 ´1 , and since the right-hand side of this inequality converges to
np2`α1 q
2 ` α1 , for n large enough we have also log2 n ě n´α1 ´1 . Thus, the inequality
(4.37) is established, and we may now complete our estimation:
ˇż ˇ ż8
ˇ 8 p´t ˘ i q α ˇ 1 C
ˇ ˇ
ˇ ` ˘ dt ˇ ď C ¨ dt “ .
ˇ log2 n 1 ` t˘i n`1 ˇ k 2
log2 n n t nk log2 n
n
By (4.35), for large n the integral (4.34) equals a sum of terms of the form
ż
ck p´tq´α e´t tk dt,
H :ℜptqďlog2 n
“except for the above error term”. We would of course like to argue that these
integrals “essentially” are equal to
ż
2π i 2π i ¨ pα ´ kq ¨ pα ´ k ` 1q ¨ ¨ ¨ pα ´ 1q
p´tq´α`k e´t dt “ “
H Γ pα ´ kq Γ pαq
(by Hankel’s formula (4.28) and the functional equation (4.21), respectively). In
order to see this, we argue that the “cut-off remaining integral” (using the same
manipulations as in the proof of Theorem 4.3.24) leads to the real integral
ż8
2i sin pπαq e´t p´tq´α dt,
log2 n
whose integrand clearly becomes very small if t becomes large:
ˇ ´t ˇ
ˇe p´tq´α ˇ ď e´t tℜpαq ď 1
n k t2
98 4. ASYMPTOTIC ENUMERATION
(since limtÑ8 e´t tℜpαq`2 “ 0). Hence, we can bound the error as before.
If we put everything together, we arrive at the claim.
E XAMPLE 4.3.26 (Catalan numbers). Applying (4.32) to the generating function of
the Catalan numbers, ?
ÿ 1´ 1´4¨z
n
Cn´1 ¨ z “ ,
2
ně1
we have α “ ´ 12 . Setting y “ 4 ¨ z, we obtain:
4n n
Cn´1 “ ´ Jy K p1 ´ yq´α
2
´3{2
ˆ ´ ¯˙
n´1 n? 3 25 105 ´4
»4 ¨ ¨ 1` ` ` `O n .
π 8n 128n2 1024n3
R EMARK 4.3.27. For the coefficients ek in (4.31) we have
2k
ÿ
ek “ λk,l pα ´ 1q ¨ pα ´ 2q ¨ ¨ ¨ pα ´ lq
l“k
q y
with λk,l “ vk t l et ¨ p1 ` vtq´1´1{v .
T HEOREM 4.3.28 (Asymptotics for coefficients of standard functions, 2). Let α P
Cz t0, ´1, ´2, . . . u. We consider the “standard function”
ˆ ˙β
´α 1 1
f pzq “ p1 ´ zq log . (4.39)
z 1´z
Then, asymptotically as n Ñ 8, we have
˜ ¸
8
ÿ
nα´1 Ck pαq
Jzn K f pzq „ plog nq β 1` , (4.40)
Γ pαq logk n
k“1
` β˘ ˇ
dk 1 ˇ
where Ck “ k
Γ pαq ˇ .
dsk Γpsq s“α
R EMARK 4.3.29. The factor 1z in (4.39) is “just” there, in order that f pzq is a power
series, also in the case where β R Z (the logarithm is analytic in a neighbourhood of
z “ 1):
ˆ ˙β ˆ ˆ ˙˙
1 1 1 1
log “ exp β log log
z 1´z z 1´z
ˆ ˆ ˙˙
z z2 z3
“ exp β log 1 ` ` ` ` ¨ ¨ ¨
2 3 4
ˆ ˆ ˙˙
z 5z2 z3
“ exp β ` ` `¨¨¨
2 24 8
R EMARK 4.3.30. The formulae in Theorems 4.3.25 and 4.3.28 actually remain valid
also in the case α P t0, ´1, ´2, . . . u if we interpret them as limit cases, i.e., if we set
1 1
“ lim “ 0.
Γ pαq zÑα Γ pzq
4.3. SINGULARITY ANALYSIS 99
In (4.40), the limit should be taken after cancelling the Gamma–factors Γ pαq, so only
the first term vanishes: For instance, we have
˜ ¸
z 1 2 ¨ γ 1
Jzn K “´ ` `O .
log p1 ´ zq´1 n ¨ log2 n n ¨ log3 n n ¨ log4 n
S KETCH OF PROOF. In analogy to the proof of Theorem 4.3.25, we choose a
“Hankel-type” contour in Cauchy’s integral formula. After the substitution
z Ñ 1 ` nt , we obtain an integral over the “ordinary” Hankel contour H with
integrand
ˆ ˙ˆ ˙ ˆ ˙ ´ n ¯¯β ˆ ˙
1 t t ´n´1 1 t ´α ´ t ´β´n´1
f 1` 1` “ ´ log ´ 1`
n n n n n t n
¨ ˛
logp´tq β
1 ´ log n
„ nα´1 p´tq´α plog nq β ˝ ‚ e´t .
1 ´ ´tn
R r γ2
φ γ1 1 θ
n
0 1 0 1
γ4
γ3
T HEOREM 4.3.33. Let α, β P R (not C!) be arbitrary, and let f pzq be a function that
is analytic on a ∆-domain. If f in the intersection of a neighbourhood of 1 and the
∆-domain satisfies
‚ the property that, if
˜ ˆ ˙β ¸
´α 1
f pzq “ O p1 ´ zq log ,
1´z
then ´ ¯
n α´1 β
Jz K f pzq “ O n plog nq ; (4.41)
‚ the property that, if
˜ ˆ ˙β ¸
´α 1
f pzq “ o p1 ´ zq log ,
1´z
then ´ ¯
Jzn K f pzq “ o nα´1 plog nq β . (4.42)
7Philippe Flajolet, in line with his cultural background, used to call these domains
“Camembert domains”.
4.3. SINGULARITY ANALYSIS 101
for a contour γ in the ∆-domain, which encircles the origin once in positive
direction. Concretely, we choose (see the right part of the above figure) the
contour γ “ γ1 Y γ2 Y γ3 Y γ4 which consists of the following parts:
" *
1
γ1 “ z P C : |z ´ 1| “ , |arg pz ´ 1q| ě θ ,
n
" *
1
γ2 “ z P C : ď |z ´ 1| , |z| ď r, |arg pz ´ 1q| “ θ ,
n
γ3 “ tz P C : |z| “ r, |arg pz ´ 1q| ě θu ,
" *
1
γ4 “ z P C : ď |z ´ 1| , |z| ď r, |arg pz ´ 1q| “ ´θ ;
n
where 1 ă r ă R and φ ă θ ă π2 . We treat each of the individual corresponding
parts of the integral separately:
ż
j 1 f pzq
fn “ dz,
2π i γ j zn`1
j “ 1, 2, 3, 4. Then we have:
ˆ ˙ ˜ˆ ˙ ¸
1 1 ´α ´ ´ ¯¯
f n1 “ O ¨O plog nq β “ O nα´1 log n β ,
n n
´ ¯
since f “ O nα plog nq β along γ1 by assumption, the length of the integration
´ ¯
path is O n1 , and z´n´1 “ O p1q along γ1 .
Furthermore, for r1 :“ n ¨ r we have
ˇ ˇ ż r1 ˇ ˇ´α ˇ ˇ β ˇˇ iθ t ˇˇ´n´1
ˇ 2ˇ 1 ˇtˇ ˇ n
ˇ ˇ ˇlog ˇ ˇ1 `ˇ e ˇ
ˇ fn ˇ ď dt Ðlog nt ďlog n
2πn 1 ˇ n ˇ t ˇ n ˇ
ż8 ˇ ˇ´n´1
1 α´1 β ´α ˇ
ˇ eiθ t ˇˇ ˇ ˇ
ˇ eiθ t ˇ
´ iθ ¯
e t
ď n plog nq t ˇ1 ` ˇ dt Ð ˇ1` n ˇě1` ℜ n
2π 1 n
ż8 ˇ ˇ
1 α´1 β ´α ˇ
ˇ t cos θ ˇˇ´n´1
ď n plog nq t ˇ1 ` dt.
2π 1 n ˇ
Since
ż8 ˇ ˇ ż8
ˇ t cos θ ˇˇ´n´1
lim t ´α ˇ
ˇ1 ` n ˇ dt “ t´α e´t cos θ dt ă 8, Ð because 0ăθă π2
nÑ8 1 1
it follows that ´ ¯
f n2 “ O nα´1 plog nq β ,
and completely analogously this also holds for f n4 .
` ˘
Along γ3 , the function f pzq is bounded, but z´n “ O r´n ; consequently f n3 is
“exponentially small”.
The estimation with o (instead of O) works in an analogous fashion. (It is how-
ever slightly more complicated: the straight pieces must be “decomposed” in a
similar way as in the proof of Theorem 4.3.25. . . .)
102 4. ASYMPTOTIC ENUMERATION
T “ ` `
T T T
This decomposition can be directly translated into the following functional equation for
the (ordinary) generating function T pzq:
´ ¯
T pzq “ z 1 ` T pzq ` T 2 pzq .
Its solution is
a
1´z´ p1 ` zq p1 ´ 3zq
T pzq “ .
2z
´1 0 1
3
we obtain
c´ ¯
1 ´ 1´u
3 ´ 1´u
1` 3 u
T puq “
2 1´u
b 3
2 ` u ´ 3 13 u p4 ´ uq
“ Ð2`u“2p1´uq`3u
2 p1 ´ uq
? ? b` ˘
3u ´ 3 u 2 1 ´ u4
“ 1` Ð2 binom. th. “trivial”
2 p1 ´ uq
? ?
3u ´ 2 3 u p1 ` O puqq
“ 1` p1 ` O puqq
2
? ? 3u ´ ¯
“ 1´ 3 u` ` O u3{2 .
2
Altogether, this leads to the expansion (now again in z):
? ? 3 ´ ¯
T pzq “ 1 ´ 3 1 ´ 3z ` p1 ´ 3zq ` O p1 ´ 3zq3{2 .
2
Altogether, we obtain
c ´ ¯
n n 3 n ´ 52
Tn “ Jz K T pzq “ 3 `O 3 n
4πn3
for the coefficient Tn of zn in T pzq.
E XAMPLE 4.3.35 (2-regular graphs). The labelled species of 2-regular graphs can be
interpreted as sets of cycles of length at least 3 (“ cycles in the sense of graph theory),
where we identify cycles if they differ only by a reflection. Hence, the corresponding
104 4. ASYMPTOTIC ENUMERATION
We basically proceed as in the previous example and apply Theorem 4.3.25 to the stan-
dard functions
This gives
3 1 3 5 ´ 5¯
f n “ e´ 4 ? ´ e´ 4 ? ` O n´ 2 .
πn 8 πn3
Hence, the number of 2-regular graphs is asymptotically equal to n! ¨ f n .
3 7
6 1 4
2 5
n!
There are exactly n´1 different such “rounds” with n children (since cyclic permuta-
tions of the children that sit around in a circle are identified); the exponential generating
function of such “single rounds” is therefore
z3 z4
z2 ` ` `¨¨¨ ,
2 3
4.4. SADDLE POINT METHOD 105
and “children’s rounds” are simply sets of such “single rounds”. Consequently, the
corresponding generating function equals
ˆ ˙
2 z3 z4
f pzq “ exp z ` ` ` ¨ ¨ ¨
2 3
ˆ ˙
1
“ exp z log
1´z
1
“ p1 ´ zq´z “ p1 ´ zq´z`1
1´z
1
“ exp pp1 ´ zq log p1 ´ zqq .
1´z
This function is analytic in CzRě1 , and it has a singularity at ζ “ 1. We expand the
exponential function about z0 “ 0:
plogp1´zqp1´zqq2
1 ` log p1 ´ zq p1 ´ zq ` 2! `¨¨¨
f pzq “
1´z
1 ÿ plog p1 ´ zqqn p1 ´ zqn´1
“ ` log p1 ´ zq `
1´z n!
ně2
ˆ ˙
1 1 1 2 1
“ ´ log ` ¨ p1 ´ zq ¨ log
1´z 1´z 2 1´z
ˆ ˆ ˙˙
1
` O p1 ´ zq2 ¨ log3 (see [3, p. 397]).
1´z
1 1
The part 1´z ´ log 1´z contributes 1 ´ n1 to the coefficient. After application of Re-
mark 4.3.30, we obtain an estimation for the rest from (4.41) (with α “ ´1 and β “ 2,
respectively α “ ´2 and β “ 3). Thus, altogether we have
ˆ ˆ ˙˙ ´ ¯
1 ´2 1 ´3 2
f n “ 1 ´ ´ n ¨ plog nq ¨ 1 ` O ` O n ¨ log n
n log n
1 ´ ¯
´2 ´2
“ 1 ´ ´ n ¨ log n ` O n .
n
4.4. Saddle point method
Singularity analysis turned out to be extremely useful for the determination of
asymptotics. However, what do we do if the generating function
z2
‚ either has no singularities at all, such as e.g. ez (sets), ez` 2 (involutions)
or ee ´1 (Bell numbers),
z
4.4.1. Heuristics: contour through saddle point. In such cases, the follow-
ing heuristic argument may help. Starting point is again Cauchy’s integral for-
mula: ż
n 1 f pzq
Jz K f pzq “ dz. (4.44)
2π i C zn`1
106 4. ASYMPTOTIC ENUMERATION
The idea is to deform the contour C around 0 with winding number 1 inside the
domain G where the function f is analytic such that
f pzq
‚ the integrand zn`1 is so small on a part C1 of the contour C such that the
ş f pzq
part 2π1 i C1 zn`1 dz of the integral is negligible in the asymptotic consid-
erations,
‚ so that the asymptotically relevant part of the integral is the one coming
from the “rest” C2 of the contour C (where moreover this rest should
be easily computable, at least asymptotically).
This sounds somewhat unclear and vague, so that one cannot imagine how
this is going to work in practice. As a matter of fact, this vague idea can be
successfully carried through in so many cases such that it acquired a name:
saddle point method.8
So, to begin with, we have to estimate the modulus of the integrand in (4.44),
ˇ ˇ
ˇ f pzq ˇ
ˇ ˇ
ˇ zn`1 ˇ .
This real valued function appears graphically as a “mountain landscape” over
the complex plane; according to the above vague requirement, the path de-
scribed by the contour C should “almost always be close to the plane, and only
once lead through a mountain pass”. (See the following figure.)
f pzq
So, let us consider the function g pzq :“ zn`1
. By Taylor’s theorem, we have
1 ´ ¯
1 2 2 3
g pzq “ g pz0 q ` pz ´ z0 q g pzo q ` pz ´ z0 q g pz0 q ` O pz ´ z0 q for z Ñ z0 .
2
We write |g| in a neighbourhood of z0 :
Case 1: g pz0 q “ 0. In this case, the function |g| has a global minimum in z0 , see
the following figure:
8Recall that a trick that works at least three times becomes a method . . .
4.4. SADDLE POINT METHOD 107
g1 pz q
Case 2: g pz0 q ‰ 0 and g1 pz0 q ‰ 0. Then we write gpz 0q “ λeiφ in polar coordi-
0
nates, and we see that
ˇ ´ ¯ˇ
ˇ ˇ
|g pzq| “ |g pz0 q| ¨ ˇ1 ` r ¨ λ ¨ eipφ`θq ` O r2 ˇ .
For small r along a line segment through z0 with angle θ, we have in first ap-
proximation:
‚ for θ “ ´φ ˘ π, the function |g pzq| grows, from a minimum „ |g pz0 q| p1 ´ rλq
to a maximum „ |g pz0 q| p1 ` rλq; ? ` ˘
‚ for θ “ ´φ ˘ π2 , we have |1 ˘ rλi| “ 1 ` r2 λ2 “ 1 ` O r2 “ 1 ` o prq,
so that the function g pzq is “essentially constant”.
g2 pz q
Case 3: g pz0 q ‰ 0 and g1 pz0 q “ 0. Let 21 gpz 0q “ λeiφ in polar coordinates. Then
0
ˇ ´ ¯ˇ
ˇ 2 ip2θ`φq ˇ
|g pzq| “ |g pz0 q| ˇ1 ` λr e ` O r3 ˇ .
“:In0 “:In1
2 ?
5n
´n δ2
„e “e . ´ 2
´ ´ 1
¯¯
Altogether we obtain In1 “ O exp ´ 21 ¨ n 5 .
Now we consider In0 :
żδ żδ ´ 2 ¯
npeiθ ´1´iθ q n ´ θ2 `Opθ 3 q ř
e dθ “ e dθ Ðe z “ zn
n n!
´δ ´δ
ˆ ˙
żδ 2 ´ 51
´n θ2 `O n
´ 65
“ e dθ Ð|θ 3 |ďδ3 “n
´δ
żδ ´ ´ ¯¯
θ2 1
“ e´n 2 ¨ 1 ` O n´ 5 dθ
´δ
´ ´ ¯¯ ż δ θ2
´ 15
“ 1`O n e´n 2 ¨ dθ.
´δ
We make small a computation on the side:
ż8 2
ż8
e´ 2 ¨pδ `2δt`t q dt Ðθ“t`δ
n 2 2
´n¨ θ2
e dθ “
δ 0
2
ż8
´ nδ2
ďe e´nδt dt
ˆ0 ˙ˇ
nδ2 1 ´nδt ˇˇ8
“e ´ 2
´ e ˇ
nδ 0
ˆ 1
˙
1 ´ nδ 2
2
“ e 2 “O e ´n 5
. Ðδ“n´ 5
nδ
Hence, we have
żδ ż8 ´ ¯
2
´n¨ θ2 ´n¨ θ2
2
´n1{5 ?n
e dθ “ e dθ ` O e Ðs“θ 2
´δ ´8
c ż ´ ¯
2 8 ´s2 ´n1{5
“ e ds `O e
n looooomooooon
´8
?
π
c ´ 1¯
2π
“ ` O n´ 5 . ÐOpe´x q“Op 1x q
n
If we put all this together, then we obtain the following variant of Stirling’s formula:
1 ´ n ¯´n 1 ´ ´ ¯¯
“ ? 1 ` O n´1{5 .
n! e 2πn
110 4. ASYMPTOTIC ENUMERATION
E XAMPLE 4.4.3. Let un be the number of ordered set partitions; these are like ordinary
(set) partitions except that the order of the elements in each block is relevant. In the
language of species, these are
sets ppermutationens 1 q ,
from which we directly read off the (exponential) generating function:
ÿ zn ˆ ˙
z
un “ exp .
ně0
n! 1 ´ z
Altogether, we have
1 1 ´ ¯
z1 “ 1 ´ ? ` ` O n´3{2 .
n 2 pn ` 1q
Again, we approach the situation “generously”. We choose the “asymptotic saddle
point”
1
r “ 1´ ?
n
on the real axis, which is at the same time the radius of the circle about 0 that we choose
as integration contour (see the preceding example, in particular (4.45)). Then Cauchy’s
integral formula yields
ż ˆ ˙
un 1 2π r ei θ
“ exp ´ niθ r´n dθ.
n! 2π 0 1 ´ r ei θ
As in the preceding example, we split this integral “skilfully”:
ż 2π ż 2π´δ żδ ż 2π´δ
“ “ ` ,
0 ´δ ´δon loδomoon
loomo
“:In0 “:In1
r ei θ
“ r ¨ eiθ ` r2 ¨ e2iθ ` r3 ¨ e3iθ ` ¨ ¨ ¨ Ðeiθ “1`iθ´ θ2 `Opθ3 q
2
1 ´ re i θ
¨ ˛ ¨ ˛
˚ ‹ ˚ ‹
r ` r2 ` r3 ` ¨ ¨ ¨‚` iθ ˚
“ 1 ˝loooooooomoooooooon ˝ r ` 2r 2
` 3r 3
` ¨ ¨
loooooooooomoooooooooon ¨‹
‚
r ? r ?
1´r “ n´1 p1´rq2
“n´ n
¨ ˛
θ2 ˚ ‹
˚ 2 2 2 3 ‹
´ r ` 2 r ` 3 r ` ¨ ¨ ¨‹
˚looooooooooomooooooooooon
2 ˝ ‚
rp1`rq ?
“ n´3n`2n3{2
p1´rq3
¨ ˛
˚ ‹
´ ¯˚ ‹
3 ˚ 3 2 3 3 ‹
r ` 2 r ` 3 r ` ¨ ¨ ¨ ‹ pn Ñ 8q .
` O θ ˚ looooooooooomooooooooooon
˚ ‹
˝ rp1`4r`r2q ˆ ˙ ‚
p1´rq4
“O 1
p1´rq4
“O p n 2 q
112 4. ASYMPTOTIC ENUMERATION
` 3˘ ` 2˘ ´ ¯
3
For the part In0 ,
we have |θ| ď δ3 “ n´21{10 , hence O θ O n “ O n ´1{10 .
Consequently, we get
ˆ ˙ ˆ ˙
żδ rp1`rq θ3
1´r `iθ
r r
1 ´n ´θ 2 `O
In0 r n “ e p1´rq2 2p1´rq3 p1´rq4
dθ
2π ´δ
żδ ´? ´ ¯ ´ ¯¯
1 ?
“ exp n ´ 1 ´ iθ n ´ θ 2 n3{2 ` O pnq ` O n´1{10 dθ
2π ´δ
żδ ´? ´ ¯¯
1
“ exp n ´ 1 ´ θ 2 n3{2 ` O n´1{10 dθ,
2π ´δ
where the last simplification follows from
? ´ ¯ ´ ¯
´2{10 2 ´4{10
θ n“O n , θ O pnq “ O n .
we eventually arrive at
?
´
e2 ´n´ 12 ¯¯
In0 “ ? 3{4 1 ` O n´1{10 .
2 πn
4.4. SADDLE POINT METHOD 113
For the remaining part of the integral, which we denoted by In1 , we collect the following
z 1
facts: we have (of course) e 1´z “ e 1´z ´1 and (equally obviously) 1´z 1 1´z
“ p1´zqp1´zq ;
from this, it follows immediately that
ˇ z ˇ 1 ˇˇ ˆ
1 ´ r e ´i θ
˙ˇ
ˇ 1
ˆ
1 ´ r cos θ
˙
ˇ 1´z ˇ
ˇe ˇ “ ˇˇexp ˇ “ exp .
e 1 ´ 2r cos θ ` r2 ˇ e 1 ´ 2r cos θ ` r2
1´rz
The derivative of the function 1´2rz`r 2 (z P r´1, 1s, 0 ă r ă 1), with respect to z
rp1´r q
2
is 2 ą 0. Consequently, it is (strictly) monotone increasing. However, this
p1´2rz`r2 q
implies that, for z “ reiθ and δ ď θ ď 2π ´ δ, we have
ˇ z ˇ ˆ ˙
ˇ 1´z ˇ 1 ´ r cos δ
ˇe ˇ ď e exp
´1
,
1 ´ 2r cos δ ` r2
and ˜ ¸
1 ´ r cos n´7{10 ´? ¯
1{10
exp “ exp n´n ` O p1q .
1 ´ 2r cos n´7{10 ` r2
(This is seen — after the substitutions r “ 1 ´ ?1n and n “ z´10 — from the Laurent
series expansion
` 5 ˘ ` ˘ ´ ¯
z ´ 1 cos z7 ` 1 1 1
“ 5 ´ ` O z3 ,
z ´ 2z ` 2 pz ´ 1q cos pz q ` 2
10 5 5 7 z z
and the “backward substitution” z “ n´1{10 .) Furthermore,
2π ´n ?n´n1{10 Op1q ´ ? 1{10 ¯
In1 ď r e e “ O e2 n´n ,
2π
with the expansion of r´n as before. This may not look very spectacular: however,
asymptotically, In1 is “much smaller” than In0 :
ˇ 1ˇ ˇ ˇ
ˇ In ˇ
ˇ ˇ ď C ˇˇn3{4 e´n1{10 ˇˇ ,
ˇ I0 ˇ
n
hence ´ ¯
In1 “ o In0 pn Ñ 8q .
R EMARK 4.4.4 (Longer asymptotic expansion). We may take more terms in the
Taylor series expansion in θ — for ez (see Example 4.4.2) we get for instance
ż8 ż8 ´ 2 3 4 θ 5 `O θ 6
¯
npeiθ ´1´iθ q n ´ θ2 ´i θ6 ` θ24 `i 120 p q
e dθ “ e dθ Ðez “1`z`Opz2 q!
´8 ´8
ż8 ˆ ´ ¯˙
2
´n θ2 θ3 θ4 θ5
“ e 1 ´ ni ` n ` ni ` nO θ 6 dθ
´8 6 24 120
? ż
2π n 8 4 ´θ2 n
niO p. . . q `
“ ? ´ loooomoooon θ e 2 dθ
n 24 ´8
“0, see below.
ˆ ż8 ˙
6 ´θ 2 n2
`O n θ e dθ .
´8
114 4. ASYMPTOTIC ENUMERATION
ř 4.4.2. nHayman’s theorem. Our typical situation is the following: let f pzq “
ně0 f n z be analytic about z “ 0; we want to express the coefficient f n by using
Cauchy’s integral formula, where we integrate along a circle z pθq “ reiθ :
ż ` ˘
r´n π f reiθ
fn “ dθ.
2π ´π eniθ
The subsequent steps follow the general pattern:
‚ show that the “tails” of the integral (i.e., θ R p´δ, δq) are “asymptotically
negligible in comparison to the main part”;
‚ show that, for the main part of the integral (i.e., θ P p´δ, δq), there exists
2
an approximation of the integrand of the type e´θ ;
‚ show that the tails of these “new, approximative integrals” (i.e., the
integration domains θ ă ´δ and θ ą δ) are again asymptotically negli-
gible;
see the following schematic illustration.
ignore tails
I:
´π ´δ δ π
2
approximate by e´t
II: in the centre
´δ δ
add tails
III:
´8 `8
We now make situations in which this “general pattern” works more precise.
For r ě 0 in the circle of convergence with f prq ą 0 and small θ, we have the
series expansion in θ
´ ´ ¯¯ 8
ÿ piθq p
log f reiθ “ log p f prqq ` α p prq . (4.46)
p!
p“1
4.4. SADDLE POINT METHOD 115
Of particular interest for our purposes are the first two terms, which, using the
notation h prq :“ log p f prqq, we write
f 1 prq
a prq :“ α1 prq “ rh1 prq “ r (4.47)
f prq
b prq :“ α2 prq “ r2 h2 prq ` rh1 prq
˜ ˆ 1 ˙2 ¸
f 1 prq f 2 prq f prq
“r ` r2 ´ . (4.48)
f prq f prq f prq
In practice, the first thing to do is to find the function δ prq: with the notation
from (4.46), by H2 and H3 we should for r Ñ R on the one hand have
α3 prq2
lim “ 0.
rÑR α2 prq3
For the “saddle point Ansatz”, we need
´ ¯
´1{2 ´1{3
α2 prq “ o pδ prqq and δ prq “ o α3 prq , (4.49)
2
1) are Hayman-admissible functions. Not Hayman-admissible is for example ez
(becomes “too large” about π and violates H3).
116 4. ASYMPTOTIC ENUMERATION
ř
T HEOREM 4.4.7 (Hayman). Let f pzq “ něo f n zn be Hayman-admissible. For n P
N, let ζ “ ζ pnq be the largest (and, thus, unique) solution of the equation
f 1 pzq
a pzq “ z “n
f pzq
f pζq
fn „ a ,
ζ n 2πb pζq
d 2 d
where b pzq “ z2 dz 2 log f pzq ` z dz log f pzq.
where we set δ “ δ prq. The part In1 can be immediately approximated by as-
sumption H3:
˜ ¸
ˇ ˇ 1 f prq
ˇ 1ˇ
ˇ In ˇ ď 2π ¨ o a .
2π b prq
The assumptions H2 and H3 both hold for θ “ δ prq. Together, they yield
a
b prqeiδprqaprq´δprq
2
bprq{2
Ñ 0 ùñ δ prq2 b prq Ñ 8 for r Ñ R.
4.4. SADDLE POINT METHOD 117
b ?
bprq
We perform the substitution t “ 2 θ (i.e., dθ “ ? 2
dt), and we obtain for
bprq
r Ñ R:
¨ b ˛
ż δ bprq d b ´ ¯
f prq ˝ 2 2 ´t `it bprq paprq´nq
2 2
In0 “ ¨ b e dt ` o b prq´1{2 ‚
2π ´δ
bprq
2
b prq
¨ b ˛
ż δ bprq b
f prq ˝ 2 ´t `it bprq paprq´nq
2 2
“ a b e dt ` o p1q‚ Ðform a square
π 2b prq ´δ
bprq
2
˜ż ´ b ¯2 ¸
8 paprq´nq2
f prq ´ t´i 2bprq1
paprq´nq ´ 2bprq ?
“ a e dt ` o p1q Ðδ bprq{2Ñ8
π 2b prq ´8
¨ ˛
˚ b
ż 8´i 1 paprq´nq ‹
f prq ˚ paprq´nq2
˚ ´ 2bprq 2bprq 2
‹
‹
“ a ˚e b e ´s
ds `o p1q ‹.
˚
π 2b prq ˝ ´8´i 2bprq paprq´nq
1 ‹
looooooooooooooomooooooooooooooon ‚
?
π
Thus, we established (4.50), and the claim of the theorem follows immediately:
for, by H1, we have limrÑR a prq “ `8, and hence
lim ζ pnq “ R;
nÑ8
we simply put r “ ζ pnq.
The convenient feature of the concept of Hayman-admissibility is that there
hold the following “closure properties”: “Large classes of functions are Hayman-
admissible”.
T HEOREM 4.4.8. Let f pzq and g pzq be Hayman-admissible functions, and let p pzq be
a polynomial with real coefficients. Then:
‚ f pzq ¨ g pzq and e f pzq are also Hayman-admissible;
‚ f pzq ` p pzq is Hayman-admissible; if the leading coefficient of p pzq is posi-
tive, then also f pzq p pzq and p p f pzqq are Hayman-admissible;
‚ If almost all Taylor coefficients of e ppzq are positive, then also e ppzq is Hayman-
admissible.
S KETCH . In the first place, one has to find the function δ prq: to this end, one
uses (4.49).
We give some examples in order to illustrate the concept of Hayman-admissibility.
E XAMPLE 4.4.9. For f pzq “ ez , we have the expansion
´ ¯ θ2
i “ er`riθ´r 2 `Oprθ q
3
f re θ
2
“ f prq ¨ eriθ´r 2 `Oprθ q .
θ 3
We have
´ ¯ 1 1
log f reiθ “ e2iθ r2 ` e3iθ r3
2 3
ˆ 2 3 ˙ ´ ¯ ´ ¯ 2 ´ ¯
r r 3 θ
“ ` ` r ` r iθ ´ 2r ` 3r
2 3 2
` O θ3
2 3 loooomoooon looooomooooon 2
aprq bprq
10The name comes from the idea to pull oneself out of one’s boots, similar to Münch-
hausen’s famous story where he pulls himself out of a swamp. This is somehow how the
method works.
120 4. ASYMPTOTIC ENUMERATION
or, equivalently,
´ ¯
n2{3 p1 ` 3ρq ` n1{3 2ρ ` 3ρ2 ` ρ2 ` ρ3 “ 0. (4.53)
Now we suppose that ρ is of larger asymptotic order than a constant. In that case, the
term 3n2{3 ρ is asymptotically larger than all other terms on the left-hand side of (4.53).
This is a contradiction because then the left-hand side cannot be zero. Hence, ρ must be
asymptotically constant. In that case, the asymptotically largest terms on the left-hand
side of (4.53) are n2{3 p1 ` 3ρq. They must (asymptotically) “cancel” each other, by
which we mean that ρ must be chosen such that they actually are of smaller order. We
infer that ρ “ ´ 13 ` o p1q.
We need to continue since, by our findings so far, the term ζ 3 (which appears in the
exponential in (4.51)) is still not approximated to an order that would tend to zero. So,
we let ζ “ n1{3 ´ 31 ` ρ, with ρ of smaller (asymptotic) order than a constant.11 If we
substitute this Ansatz in (4.53), then we obtain
1 1 2
3n2{3 ρ ` 3n1{3 ρ2 ´ n1{3 ` ρ3 ´ ρ ` “ 0. (4.54)
3 3 27
By arguments
´ that¯are analogous to the previous ones, here we find that ρ must be of
the order O n´1{3 . In that case, the asymptotically largest terms in (4.54) are 3n2{3 ρ
and ´ 13 n1{3 . They
´ must
¯ (asymptotically) cancel each other, and therefore we see that
ρ “ 91 n´1{3 ` o n´1{3 .
This is still not good enough. Another round of “bootstrap” reveals that
1 1 ´1{3 2 ´ ¯
ζ “ n1{3 ´ ` n ´ n´2{3 ` o n´2{3 .
3 9 81
If this is substituted in (4.51), we obtain
1 1 2{3 ´ 1 n1{3 ` 1
e 3 n` 2 n 6 18
f n „ n! ? n 1
, as n Ñ 8.
6πn 3 ` 2
Here we used that
` ˘n
ζ n “ nn{3 ζn´1{3
ˆ ´ ¯˙n
n{3 1 ´1{3 1 ´2{3 2 ´1 ´1
“n 1´ n ` n ´ n `o n
3 9 81
ˆ ˆ ´ ¯˙˙
n{3 1 ´1{3 1 ´2{3 2 ´1 ´1
“ n exp n log 1 ´ n ` n ´ n `o n
3 9 81
´ ¯
“ nn{3 exp ´ 13 n2{3 ` 18
1 1{3
n ` o p1q , as n Ñ 8.
11We allow ourselves the liberty of using the same symbol again although the ρ here has
nothing to do with the earlier ρ.
4.4. SADDLE POINT METHOD 121
Equality for a z “ reiθ can only hold if all powers zn “ r n eniθ (interpreted as
vectors in R2 ) point in the same direction. So let two arbitrary coefficients be
non-zero, say f n1 , f n2 ‰ 0. Then we must have
en1iθ “ en2iθ ðñ |n1 ´ n2 | θ “ 2kπ,
k
and thus θ “ 2π |n , where we may of course choose k ă |n1 ´ n2 |. Let
1 ´n2 |
p k
q “ |n1 ´n2 | be the reduced fraction. Then we have
q | n1 ´ n2 ,
or, equivalently, n1 ” n2 pmod qq.
θ2 ´ A ¯ ´ ¯
´ α2 pζq ` nα2B pζq ` O n ¨ θ 3 q,
2
for α3A pζq and α3B pζq are O p1q for n Ñ 8 since ζ does not depend on n. Now let
θ0 :“ n´2{5 . Then, for θ0 ď θ ď θ1 , we have (of course)
2 ´ ¯ ´ ¯
´ θ2 pα2A pζq`nα2B pζqq`Opn¨θ 3 q ´n1{5 Opn´1{5 q ´n1{5
e “O e ¨e “O e ,
´ 1{5
¯
for n Ñ 8. Hence also J pθ1 q ´ J pθ0 q “ O e´n for n Ñ 8, that is, exponen-
tially small. The desired coefficient is therefore´asymptotically ¯ J pθ0 q: for θ ă θ0
2
we have however (of course) θ2 ¨ α2A pζq “ O n´4{5 pn Ñ 8q — that is, the
desired coefficient is (after “appending” the asymptotically negligible integral
ends) asymptotically equal to
ż
A pζq B pζq n 8 ´ θ 2 nαB pζq A pζq B pζq n
e 2 2 dθ “ b .
2πζ N ´8 ζ N B
2πn ¨ α2 pζq
` ˘
The coefficient α2B is the third Taylor coefficient in the expansion of log B reiθ ,
hence
d2 ´ ¯ˇˇ B1 prq 2
ˆ 1 ˙2
i 2 B prq 2 B prq
α2 prq “ ´ 2 log B re ˇ
B θ ˇ
“r `r ´r ,
dθ θ“0 B prq B prq B prq
and we have
ˇ ˆ 1 ˙2
d2 ˇ
ˇ B 2 pζq B pζq λ α2B pζq
ξ“ plog B psq ´ λ log sq ˇ “ ´ ` “ .
ds2 s“ζ B pζq B pζq ζ2 ζ2
This establishes our claim.
E XAMPLE 4.4.17. Let f pzq be given by the equation
f pzq “ z ¨ e f pzq ðñ f pzq e´ f pzq “ z.
I.e., f is the compositional inverse of z ¨ e´z . By Lagrange inversion, we get
1 r ´1 z ` ´z ˘´n
fn “ z ze
n
1 r ´1 z enz
“ z
n zn
r
1 n´1 z
“ z pez qn
n
1
“ Jzn K pez qn ¨ z.
n
124 4. ASYMPTOTIC ENUMERATION
The summands s pn, kq in this sum “optically” show a similar behaviour as the function
2
e´t that strongly decreases away from zero, see the following plot for the summands
s p100, kq:
Summands for the 100-th Motzkin-Number
1.2 × 1044
1.0 × 1044
8.0 × 1043
6.0 × 1043
4.0 × 1043
2.0 × 1043
10 20 30 40 50
For the determination of the place of the “peak” of these summands, we look for the k,
for which s pn, k ` 1q {s pn, kq “ 1. This leads to the equation
k ` 1 pn ´ 2kq pn ´ 2k ´ 1q
“ 1,
k`2 pk ` 1q2
respectively equivalently
pn ´ 2kq pn ´ 2k ´ 1q “ pk ` 1q pk ` 2q .
126 4. ASYMPTOTIC ENUMERATION
Again, we are not necessarily interested in the “exact” location of the peak. With the
“approximate” Ansatz k „ λn, we are “approximatively” led to the quadratic equation
p1 ´ 2λq2 “ λ2 ,
with obvious solutions λ1 “ 1 and λ2 “ 31 . Therefore, we set k0 “ n3 , and approximate
s pn, n{3 ´ k1 q with the help of Stirling’s formula (4.1), in logarithmic form:
ˆ ˙
1 1
log n! „ ´n ` n ` log n ` plog 2πq .
2 2
We have
n!
s pn, kq “
pk ` 1q! pkq! pn ´ 2kq!
˜ˆ ˙
1
“ exp n` log n ` 1 ´ log 2π
2
ˆ ˙
1
´ n ´ 2k ` log pn ´ 2kq
2
ˆ ˙ ˆ ˙ ¸
1 3
´ k` logpkq ´ k ` logpk ` 1q .
2 2
Now we substitute k “ n{3 ´ k1 and use
´ ˆ ˙
n¯ 3z
log ´z ` “ log n ´ log 3 ` log 1 ´
3 n
(for z “ k1 and z “ k1 ´ 1) together with
z2 ´ ¯
log p1 ´ zq “ z ` ` O z3
2
to obtain (after some computation):
˜
9k21 ´ 3k1 ` 3
s pn, kq “ exp ´ 2 logpnq ` n logp3q ` 2
n
˜ ¸¸
5 logp3q k31
´ logp2πq ` `O
2 n2
5
˜ ˜ ¸¸
3n` 2 2´ 1 p9k2 ´3k1 `3q k31
“ e n 1 1 ` O .
2πn2 n2
(The expression
´ ¯ becomes so simple because many terms can be “subsumed” in the error
k31
term O n2
.) computations ok, fix this??? Consequently, we have
ÿ 3n` 25 ˆ ˆ 3 ˙˙
´ n1 p9k2 ´3k`3q k
Mn “ e 2
2
e 1`O . (4.57)
2πn n2
k
We split this sum as follows:
ÿ ÿ
Mn “ p¨ ¨ ¨ q ` p¨ ¨ ¨ q .
|k|ďn3{5 |k|ąn3{5
4.5. ASYMPTOTICS OF COMBINATORIAL SUMS 127
´Q U ¯
For |k| ą n3{5 (i.e., k “ ˘ n3{5 ` k1 with k1 ą 0), the summand for |k| “ n3{5
dominates the other summands, provided n is large enough:
´ ¯
2 ?
´ n1 9pn3{5 `k1 q ˘3pn3{5 `k1 q`3 2
“ e´9pn nq `Opn´2{5 q
3{5´1{2 `k
e 1{ .
Consequently, we can bound that part of the sum by
ˇ ˇ
ˇ ˇ ´ ¯
ˇ ÿ ˇ
ˇ p¨ ¨ ¨ q ˇ “ O n ¨ 3n n´2 e´9n1{5 .
ˇ ˇ
ˇ|k|ąn3{5 ˇ
This is exponentially small in comparison with 3n {n3{2 . The other part of the sum is
ÿ 3n` 25 2´9´ ?k ¯2 ´ ´ ¯¯
e n ¨ 1`O n ´1{5
.
3{5
2πn2
|k|ďn
Exercise 39: Let Dn,k be the number of permutations of rns, whose disjoint cycle decomposition
does not contain any cycle of length ď k. (So Dn,1 is the number of fixed–point–free permutations
of rns.)
(1) Show:
z 2 z k
ÿ Dn,k e´z´ 2 ´¨¨¨´ k
zn “ .
n! 1´z
n ě0
(2) For k fixed, what is the asymptotic behaviour of Dn,k for n Ñ 8?
Exercise 40: Let wn,k be the number of possibilities of paying an amount of n Euro using 1–Euro–
coins, 2–Euro–coins and 5–Euro–notes, where exactly k coins or notes are used (again, the order of
the coins and notes is irrelevant).
ř
(1) Determine the generating function n,kě0 wn,k zn tk .
ř
(2) If we assume that all possibilities which are enumerated by wn “ k wn,k have the same
probability: What is the asymptotic behaviour of the expected value for the number of
coins and notes which are used to pay an amount of n Euro, for n Ñ 8?
128 4. ASYMPTOTIC ENUMERATION
Hint: The fact that the dominant singularity can not be calculated explicitly is not an obstacle. One
has to continue to calculate with the dominant singularity symbolically.
Exercise 42: Let p pn, kq be the number of all (integer) partitions of n with at most k summands.
Show:
ÿ 1
p pn, kq zn “ .
n ě0
p1 ´ zqp1 ´ z2 q ¨ ¨ ¨ p1 ´ zk q
Determine the asymptotic behaviour of ppn, kq for fixed k and n Ñ 8.
1
Exercise 44: The fraction Γ pz q
is an entire function with zeroes 0, ´1, ´2, . . . . Show Weierstraß’
product representation:
8 ´
1 ź s ¯ ´s{n
“ seγs 1` e ,
Γ psq n
n “1
where γ denotes Eulers constant
γ “ lim pHn ´ log nq
nÑ8
and
n
ÿ 1
Hn “
j
j “1
denotes the n–th harmonic number.
Exercise 45: Show the reflection formula for the gamma function:
π
Γ pzq Γ p1 ´ zq “ . (4.58)
sin πz
Hint: Use the product representation of the sine:
ź8 ˆ ˙
z2
sin z “ z 1´ 2 2 . (4.59)
n “1
n π
Exercise 46: Show the duplication formula for the gamma function:
ˆ ˙
1 1
Γ pzq Γ z ` “ 21´2z π 2 Γ p2zq
2
and its generalisation
źm ˆ ˙
j 1 m´1
Γ z` “ m 2 ´mz p2πq 2 Γ pmzq .
m
j “0
4.5. ASYMPTOTICS OF COMBINATORIAL SUMS 129
Exercise 48: A Motzkin path is a lattice path, where every step is of the form p1, 0q, p1, 1q, p1, 1q
(horizontal, up– and down–steps), which starts at the origin, returns to the x–axis and never goes
below the x–axis. Let Mn be the number of all Motzkin paths of length (i.e., number of steps) n.
Show that the generating function for Motzkin paths is given by
?
ÿ 1 ´ z ´ 1 ´ 2z ´ 3z2
n
Mn z “ .
n ě0
2z2
Derive an explizit formula for Mn . Use the generating function to determine the asymptotic behaviour
of Mn for n Ñ 8.
Exercise 49: A Schröder path is a lattice path consisting of steps p2, 0q, p1, 1q and p1, ´1q (i.e.,
double horizontal, upward and downward steps) which starts at the origin, returns to the x–axis
but never falls below the x–axis. If we assume that all Schröder paths of length n have the same
probability: What is the asymptotics of the expected value of the number of steps for a Schröder
path of length n for n Ñ 8?
Exercise 50: Consider the number of cycles in the disjoint cycle decomposition of permutations of
rns on average: What is the asymptotics for this average for n Ñ 8?
řn
Exercise 51: Let Hn “ j “1 j´1 be the n–th harmonic number. Show that
ÿ 1 1
Hn z n “ log ,
n ě0
1´z 1´z
and use this result together with singularity analysis to obtain an asymptotic expansion of Hn for
n Ñ 8.
Exercise 52: Let un be the number of permutations of rns, which only have cycles of odd length in
their decomposition into disjoint cycles. Determine the asymptotic behaviour of un for n Ñ 8.
Hint: Observe that the generating function is analytic in a “Double–Delta–Domain” (i.e., in a disk
with two “dents” at the two singularities), so we have two contributions according to the Transfer
Theorem.
Hint: Compute the generating function of these sums, i.e., multiply the above expression by zn and
sum over all n ě 0; apply the binomial theorem for simplifying the double sum thus obtained.
Exercise 54: Denote by In the number of all involutions (an involution is a self–inverse bijection) on
rns.
(1) Show
ÿ zn z2
In “ ez` 2 .
n!
n ě0
(2) Use the saddle point method to determine the asymptotic behaviour of In for n Ñ 8.
Exercise 55: The exponential generating function of the Bell–numbers Bn (Bn is the number of all
partitions of rns) is
ÿ zn
“ ee ´ 1 .
z
Bn
n ě0
n!
Use the saddle point method to determine the asymptotic behaviour of Bn for n Ñ 8.
Exercise 57: The saddle point method can also be used for the asymptotics of the Motzkin numbers
Mn (see exercise 48) for n Ñ 8:
Show r z ´ ¯n r z ´ ¯n
Mn “ z 0 z ` 1 ` z ´1 ´ z 2 z ` 1 ` z ´1
and obtain a complex contour integral for Mn , which can be dealt with using the saddle point method.
[1] F. Bergeron, G. Labelle, and P. Leroux. Combinatorial Species and tree–like Structures, vol-
ume 67 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, 1998.
[2] A. Erdélyi, editor. Higher Transcendental Functions, volume 1. McGraw–Hill, 1953.
[3] P. Flajolet and R. Sedgewick. Analytic Combinatorics. Cambridge University Press, 2009.
[4] M. Fulmek and C. Krattenthaler. Diskrete Mathematik. Vorlesungsskriptum.
[5] H. Heuser. Lehrbuch der Analysis. B.T. Teubner, 12. edition, 1998.
[6] G. C. Rota and R. Stanley. Theory of Möbiusfunction. Z. Wahrsch. Theorie, 2:340–368, 1964.
[7] R. Stanley. Enumerative Combinatorics, volume 1. Wadsworth & Brooks/Cole, 1986.
133