Meromorphic Coefficient Asymptotics
Anchor (Master): Flajolet-Sedgewick 2009 *Analytic Combinatorics* (Cambridge) Ch. IV-V (rational and meromorphic asymptotics, applications to the supercritical sequence and to surjections, compositions, denumerants, and the analytic Perron-Frobenius reading); Titchmarsh 1939 *The Theory of Functions* 2e (Oxford) Ch. 7 for the meromorphic-function and Mittag-Leffler background; Odlyzko 1995 'Asymptotic enumeration methods' in *Handbook of Combinatorics* (Elsevier) §3-5 for the residue and subtraction-of-singularities technique
Intuition Beginner
The previous unit handed you a way to write a generating function by recipe: describe what your objects are made of, and a fixed dictionary spells out the series. That series packages every count at once, but it does not yet tell you the one thing you most often want — how fast do the counts grow? This unit reads the growth rate straight off the formula, without expanding a single extra term.
The key change in viewpoint is to stop treating the variable as a bookkeeping marker and start treating it as a real (in fact complex) number you can plug in. A generating function is then a genuine function with a graph, and that graph has trouble spots: places where it blows up to infinity, or otherwise stops being smooth. These trouble spots are called singularities. For the composition series from last unit, the trouble spot is at , where the bottom hits zero and the value shoots off.
Here is the headline fact. The trouble spot nearest the origin sets the growth rate of the counts. If the nearest one sits at distance from zero, then the count of size- objects grows like — one over that distance, raised to the power . For compositions, , so the counts grow like , matching the you found by hand.
Visual Beginner
Picture the generating function as a landscape over a flat map. The flat map is the plane of input values, with the origin in the middle. The height of the landscape is how big the function is at each input. Most of the map is gently rolling, but at a singularity a thin spike rises to infinity — a pole, like a flagpole punched through the surface.
Now grow a circle outward from the origin. As long as the circle stays in the rolling part of the map, the series behaves; the moment the circle's edge reaches the first spike, the series stops converging. So the radius of the largest safe circle is exactly the distance to the nearest spike. That radius is the number , and the counts grow like .
| where the nearest spike sits | how the counts grow |
|---|---|
| at distance (compositions, ) | like |
| at distance on the positive axis | like |
| farther out (bigger ) | slower growth |
| closer in (smaller ) | faster growth |
A second spike farther away barely matters: its contribution is a smaller growth rate, swamped by the nearest one. So the nearest spike rules, and a useful bonus is that for honest counting series — whose counts are never negative — the nearest spike always sits on the positive axis, straight out to the right, where it is easy to find.
Worked example Beginner
Let us read the growth rate of surjections — the ways to split labelled tasks into an ordered list of nonempty groups (an ordered set partition). The number of these is called , and the labelled symbolic method gives the generating function , with the count attached to the peg . We will not need that derivation; we only need to find the nearest spike.
Step 1. Find where the function blows up. The value shoots to infinity exactly when the bottom is zero, that is when . On the positive axis this happens at , which is about . That is a spike.
Step 2. Check it is the nearest one. Any other solution of needs the same size for , so it sits at least as far out as (the others are off the axis and farther from the origin). So the nearest spike is at distance .
Step 3. Read the growth rate. The counts (per the peg) grow like . Since , the growth base is about .
Step 4. Sanity-check the size. The first few surjection numbers are , , , , , . Dividing each by the one before gives ratios drifting up toward roughly — and indeed the full count per object is about , whose successive ratios are , growing because of the . The base is the part the nearest spike pinned down.
What this tells us: locating one spike — solving — handed us the exponential growth base of the surjection numbers with no recurrence and no table.
Check your understanding Beginner
Formal definition Intermediate+
A generating function with positive radius of convergence is treated as an analytic function of a complex variable on its disc of convergence, following the formal-power-series substrate of 40.08.01 now read analytically. Two analytic facts anchor everything below; both are imported, not reproved here.
Cauchy's coefficient formula. If is analytic in a disc containing and is a positively oriented simple loop around inside the disc of analyticity, then
$$
f_n = [z^n]f = \frac{1}{2\pi i}\oint_\lambda \frac{f(z)}{z^{n+1}},dz,
$$
the Cauchy integral formula of 06.01.02 applied to extract the -th coefficient. The radius of convergence of the series equals the modulus of a singularity of nearest the origin: is analytic on and has at least one singularity on .
A point on the circle is a singularity of if has no analytic continuation to any neighbourhood of . The function is meromorphic on a domain (in the sense of 06.01.05) if it is analytic there except at isolated poles: points near which with analytic at , ; the integer is the order of the pole, and is its residue. A pole of order is simple. A rational function ( polynomials, , ) is meromorphic on with poles exactly at the zeros of .
The dominant singularities of are those of smallest modulus, all lying on ; the subdominant singularities lie strictly farther out. The notation , , , , and (radius / dominant-singularity modulus) is recorded in _meta/NOTATION.md.
Exponential Growth Formula. Writing for the radius of convergence, , so with subexponential (). The first principle of coefficient asymptotics is this statement: the location of the dominant singularity fixes the exponential growth factor . The second principle, the subject of 40.08.04, is that the nature (type) of the dominant singularity fixes the subexponential factor .
Counterexamples to common slips Intermediate+
"The nearest singularity can always be read off the positive real axis." Only for series with nonnegative coefficients — there Pringsheim's theorem guarantees a singularity at itself. A series such as has coefficients of mixed sign, and its dominant singularities are the non-real points ; nothing sits on the positive axis. Combinatorial OGFs, having nonnegative counts, are the favourable case.
"A pole of order gives growth ." The exponential factor is regardless of order, but the subexponential factor depends on the order: a pole of order at contributes . A simple pole () gives a pure ; a double pole gives .
"All singularities contribute equally." Only the dominant ones contribute to the leading asymptotic. A subdominant pole at contributes , exponentially smaller than , so it appears only in lower-order correction terms obtained by subtracting the dominant singular part.
"There is always a unique dominant singularity." Several singularities can share the modulus . When two simple poles sit at , their residues combine into an oscillating factor -style term; more generally poles equally spaced on produce coefficients that are nonzero only on a residue class mod (periodicity). Pringsheim guarantees one of these is the positive real point, but it need not be alone.
Key theorem with proof Intermediate+
The signature theorem is the residue expansion of meromorphic coefficients: it converts coefficient extraction into a finite sum of residues by deforming the Cauchy contour outward past the poles, and the dominant pole's residue is the leading asymptotic term [Flajolet-Sedgewick Ch. IV.5].
Theorem (meromorphic coefficient asymptotics). Let be analytic at and meromorphic in the closed disc , with no pole on , and let its poles in be . Then $$ [z^n]f = -\sum_{k=1}^{p}\mathrm{Res}\Big(\frac{f(z)}{z^{n+1}};, z = s_k\Big) + \frac{1}{2\pi i}\oint_{|z|=R}\frac{f(z)}{z^{n+1}},dz, $$ and the remainder integral is . In particular, if the dominant pole is the unique pole of smallest modulus and has order with singular part , then $$ [z^n]f \sim \frac{c_{-m}}{(m-1)!},(-1)^m,\rho^{-m}\binom{n+m-1}{m-1}\rho^{-n} \sim \frac{c_{-m}(-1)^m}{(m-1)!,\rho^{m}},n^{m-1}\rho^{-n}. $$
Proof. By Cauchy's coefficient formula, for any with (a circle inside the disc of analyticity at the origin). The integrand is meromorphic in with poles at (from , with residue — but that pole is inside the small circle and already accounted as the integral itself) and at . Deform the contour from outward to . By the residue theorem of 06.01.03, the two circular integrals differ by times the sum of residues of at the poles crossed, namely :
$$
\frac{1}{2\pi i}\oint_{|z|=R}\frac{f(z)}{z^{n+1}},dz - \frac{1}{2\pi i}\oint_{|z|=r}\frac{f(z)}{z^{n+1}},dz = \sum_{k=1}^{p}\mathrm{Res}\Big(\frac{f(z)}{z^{n+1}};, s_k\Big).
$$
Rearranging and using the inner integral gives the stated identity. For the remainder, is continuous hence bounded on the compact circle , say ; then , which is .
For the asymptotic form, compute the dominant residue. Near write ; then has the residue at equal to the coefficient of in plus lower-order pole contributions. Expanding in powers of and extracting the coefficient gives $$ \mathrm{Res}\Big(\frac{f}{z^{n+1}}; \rho\Big) = c_{-m},\rho^{-n-1}\binom{-n-1}{m-1}\rho^{-(m-1)} = c_{-m}\frac{(-1)^{m-1}}{\rho^{n+m}}\binom{n+m-1}{m-1}, $$ using . Since and the subdominant residues are each with , the dominant term is , and gives the claimed form.
Corollary (simple dominant pole). If is a simple dominant pole with near (equivalently residue of at equal to ), then . For compositions the pole is at , simple, with and local form , , recovering exactly for .
Bridge. This theorem is the foundational reason rational and meromorphic enumeration are governed by a finite amount of local data: the entire infinite sequence is captured, up to an exponentially small error, by the residues of at finitely many poles, and the dominant pole's order and residue alone fix the leading asymptotic. This is exactly the analytic completion of the partial-fraction reading of 40.01.04, where a rational was expanded as and its coefficients read as ; the residue at each reciprocal root is precisely the partial-fraction coefficient, so the residue sum is the partial-fraction sum, and the dominant root with largest modulus is the reciprocal of the nearest pole, the central insight tying spectral radius to growth rate. The simple-pole corollary builds toward the supercritical-sequence schema previewed in 40.08.01, where acquires a simple pole at the root of , and it appears again in 40.08.04 as the rational/meromorphic base case of the singularity-analysis transfer, where algebraic and logarithmic singularities replace poles. Putting these together, the residue expansion generalises the dominant-eigenvalue asymptotics of the transfer matrix of 40.01.04: the largest eigenvalue is the reciprocal of the nearest pole of , so meromorphic coefficient asymptotics and Perron-Frobenius are the same statement read analytically.
Exercises Intermediate+
Advanced results Master
Beyond the single dominant pole, the meromorphic theory extends through the systematic subtraction of singular parts, the treatment of several singularities on the circle of convergence, the supercritical-sequence schema as the canonical simple-pole source, the denumerant / multiple-pole-at-roots-of-unity case, and the analytic recasting of Perron-Frobenius.
Theorem 1 (full singular expansion by subtraction of singularities). Let be meromorphic in with poles ordered by increasing modulus, of order and singular part . Then is analytic in , so its coefficients are , and $$ [z^n]f = \sum_{j=1}^{p}[z^n]S_j(z) + O(R^{-n}), \qquad [z^n]\frac{c_{-k}}{(s_j - z)^k} = c_{-k}\binom{n+k-1}{k-1}s_j^{-n-k}, $$ giving an asymptotic expansion to any prescribed exponential precision by taking past successively more poles. Each pole contributes a polynomial-in- times , and ordering by orders the terms by decreasing exponential weight [Flajolet-Sedgewick Ch. IV.5, Odlyzko 1995 §4].
Theorem 2 (several dominant singularities and periodicity). If has poles of equal smallest modulus , located at with on the unit circle, the leading coefficient asymptotics is the sum of the dominant residues, . When the are roots of unity (the typical case for combinatorial OGFs with a periodicity / gcd structure in the allowed sizes), the factor is periodic in , so is a periodic modulation of ; the period is the order of the group generated by the . A function has all its dominant singularities at the -th roots of , the structural source of period- coefficients [Flajolet-Sedgewick Ch. IV.6].
Theorem 3 (supercritical sequence, analytic form). If with analytic at , , radius of convergence , and , then has a dominant singularity at the unique with , which is a simple pole because (positive coefficients). By the simple-pole corollary,
$$
[z^n]F(z) = [z^n]\frac{1}{1 - A(z)} \sim \frac{1}{\sigma A'(\sigma)},\sigma^{-n},
$$
the residue of at producing the clean growth previewed combinatorially in 40.08.01. The number of -components in a random size- object is then asymptotically Gaussian with mean and variance linear in [Flajolet-Sedgewick Ch. V.2].
Theorem 4 (denumerants and poles at roots of unity). The denumerant counting representations in nonnegative integers has OGF , a rational function all of whose poles lie on the unit circle (at roots of unity), with the pole at of order dominant. Hence , a quasi-polynomial in of degree : the order- pole at gives the polynomial growth , while the lower-order poles at the other roots of unity contribute the periodic (quasi-polynomial) fluctuations. This is the Sylvester-Schur denumerant, the prototype of Ehrhart quasi-polynomials [Flajolet-Sedgewick Ch. IV.6].
Theorem 5 (analytic Perron-Frobenius). For a strongly connected weighted digraph with nonnegative transfer matrix , the walk-generating function of 40.01.04 is rational with denominator , the eigenvalues of . Its dominant pole is , the reciprocal of the Perron-Frobenius eigenvalue , which is simple by Perron-Frobenius; therefore with built from the left/right Perron eigenvectors. The meromorphic residue at is exactly the Perron projector, so coefficient asymptotics and spectral dominance coincide [Flajolet-Sedgewick Ch. V.5].
Synthesis. Meromorphic coefficient asymptotics is the foundational reason a generating function's finite local data — the location, order, and residue of its dominant pole — determines the infinite asymptotic profile of its coefficients. The central insight is that Cauchy's coefficient formula turns coefficient extraction into a contour integral, and pushing that contour outward past the poles, by the residue theorem of 06.01.03, replaces the integral with a finite residue sum plus an exponentially small remainder; the dominant pole's residue is the leading term, and this is exactly the analytic completion of the partial-fraction reading of 40.01.04, where each residue at a reciprocal characteristic root is the partial-fraction coefficient. Pringsheim's theorem is the foundational reason combinatorics is the favourable case: nonnegative counts force the dominant singularity onto the positive real axis, so the growth base is found by solving a real equation — for compositions, for supercritical sequences, for surjections. This is dual to the spectral reading of the transfer matrix: the nearest pole is the reciprocal Perron-Frobenius eigenvalue, so meromorphic asymptotics and Perron-Frobenius dominance are one statement, and the residue at the dominant pole is the Perron projector — a reading that generalises the dominant-eigenvalue growth of 40.01.04.
Putting these together, the first principle (singularity location fixes the exponential factor ) is complete for the rational and meromorphic world, where every singularity is a pole and the subexponential factor is the polynomial of the pole order; the bridge to 40.08.04 is that the second principle — singularity type fixes the subexponential factor — generalises this to the algebraic and logarithmic singularities (, ) where the factor becomes , , or , and the residue calculus is replaced by the Hankel-contour transfer theorems. The recipe-to-series dictionary of 40.08.01 and the series-to-growth reading of this unit are one analytic-combinatorics pipeline.
Full proof set Master
Proposition 1 (Cauchy's coefficient formula and the radius as nearest singularity). If is analytic on , then for , ; and equals the distance from to the nearest singularity of .
Proof. On the series converges uniformly, so it may be integrated term by term against : , using . This is the Cauchy integral formula for the -th derivative, 06.01.02. For the radius statement: is analytic on by hypothesis, so the nearest singularity is at distance ; conversely if were analytic on an open disc of radius about , its Taylor series there would converge on and agree with , contradicting that is the radius of convergence (Cauchy-Hadamard). Hence a singularity sits on and none nearer.
Proposition 2 (exponential growth formula). With the radius of convergence of , , so with .
Proof. This is Cauchy-Hadamard: the radius of convergence of a power series is . Indeed if then eventually, so converges by the root test; and if the terms do not tend to . Thus . Setting gives , the subexponential normalisation.
Proposition 3 (residue expansion with remainder). Let be analytic at , meromorphic in with poles in and none on . Then , with the last term .
Proof. The function is meromorphic in with poles among . By the residue theorem 06.01.03 applied to the circle (which encloses all of them and meets none), . The residue at : near , is analytic, so and the coefficient of is , giving . Rearranging, . The bound on the boundary integral is , and by continuity on the compact circle, so the term is .
Proposition 4 (order- pole asymptotic). If the dominant pole of is unique of smallest modulus and has order with leading Laurent coefficient , then .
Proof. By Proposition 3 with chosen between and the next-smallest , with (no other pole is dominant). Near , write ; the residue of at is the coefficient of in the Taylor expansion of about , namely . The lower-order singular terms , , contribute residues of order , subleading. Thus ; using , the leading term is , and the remainder with is exponentially smaller.
Connections Master
This unit is the analytic completion of the rational-generating-function and transfer-matrix theory of
40.01.04: there the closed form was obtained algebraically by partial fractions, and here each is the reciprocal of a pole , the partial-fraction coefficient is the residue of at that pole, and the dominant root of largest modulus is the reciprocal of the nearest pole . The transfer-matrix denominator has its dominant root at the reciprocal Perron-Frobenius eigenvalue, so the analytic Perron-Frobenius theorem of the Advanced results is the meromorphic reading of that unit's spectral asymptotics.The generating functions this unit dissects are produced by the symbolic method of the prerequisite
40.08.01: the supercritical-sequence schema with its simple pole at , the composition OGF , and the partition / denumerant products all arise from the dictionary there, and this unit supplies the residue calculus that reads their growth rates. The labelled half40.08.02supplies the surjection EGF and the set-partition / permutation EGFs whose meromorphic poles this unit's method locates.The first/second-principle split previewed here is completed in the co-produced singularity-analysis unit
40.08.04: when the dominant singularity is algebraic (, from tree recursions and the implicit-function schema) or logarithmic rather than a pole, the residue calculus fails and is replaced by the Hankel-contour transfer theorems, which read the subexponential factor off the singular exponent. This unit is the rational/meromorphic base case ( a positive integer, factor ) of that general scale; the co-produced saddle-point unit40.08.06handles the remaining case of entire generating functions with no finite singularity, where growth is read from the saddle of rather than from a pole.The complex-analytic substrate is imported wholesale from the Riemann-surfaces chapter: Cauchy's integral formula
06.01.02supplies the coefficient integral, the residue theorem06.01.03supplies the contour-deformation residue sum, and the theory of meromorphic functions06.01.05supplies the pole/Laurent/residue vocabulary. This unit is one of the principal combinatorial consumers of that machinery, and the denumerant quasi-polynomials of the Advanced results connect forward to the Ehrhart theory of lattice-point counting in40.03.02, where poles at roots of unity reappear as the periodic part of the Ehrhart quasi-polynomial.
Historical & philosophical context Master
The use of complex analysis to extract the growth of combinatorial sequences begins with the circle method and the saddle-point estimates of Hardy and Ramanujan's 1918 asymptotic for the partition function , and with Darboux's late-nineteenth-century method of reading coefficient asymptotics from the singularities of a generating function on its circle of convergence. The clean statement that the radius of convergence is the distance to the nearest singularity, and that nonnegative coefficients force a singularity onto the positive axis, is Alfred Pringsheim's theorem (1894), foundational for the favourable combinatorial case [Titchmarsh 1939].
The systematic residue / partial-fraction treatment of rational and meromorphic coefficient asymptotics, with the dominant pole giving and the subtraction-of-singularities method exposing lower-order terms, was consolidated as the entry point to analytic combinatorics by Philippe Flajolet and Robert Sedgewick, whose Analytic Combinatorics (2009) [Flajolet-Sedgewick Ch. IV] organises it as the rational/meromorphic precursor to the algebraic-singularity transfer theory; Andrew Odlyzko's 1995 Handbook of Combinatorics survey [Odlyzko 1995] gives the parallel systematic account with explicit error bounds. The recognition that the dominant pole of a transfer-matrix generating function is the reciprocal Perron-Frobenius eigenvalue, tying meromorphic asymptotics to the spectral theory of nonnegative matrices, runs back to the statistical-mechanics transfer-matrix tradition of Kramers and Wannier (1941) and was made combinatorially explicit in Stanley's enumerative-combinatorics framework and in Flajolet-Sedgewick's Chapter V.
Bibliography Master
@book{flajoletsedgewick2009,
author = {Flajolet, Philippe and Sedgewick, Robert},
title = {Analytic Combinatorics},
publisher = {Cambridge University Press},
year = {2009}
}
@book{titchmarsh1939functions,
author = {Titchmarsh, E. C.},
title = {The Theory of Functions},
edition = {2},
publisher = {Oxford University Press},
year = {1939}
}
@incollection{odlyzko1995asymptotic,
author = {Odlyzko, Andrew M.},
title = {Asymptotic enumeration methods},
booktitle = {Handbook of Combinatorics},
editor = {Graham, R. L. and Gr\"{o}tschel, M. and Lov\'{a}sz, L.},
volume = {2},
publisher = {Elsevier and MIT Press},
year = {1995},
pages = {1063--1229}
}
@book{henrici1974applied,
author = {Henrici, Peter},
title = {Applied and Computational Complex Analysis, Volume 1},
publisher = {John Wiley \& Sons},
year = {1974}
}
@article{hardyramanujan1918,
author = {Hardy, G. H. and Ramanujan, S.},
title = {Asymptotic formulae in combinatory analysis},
journal = {Proceedings of the London Mathematical Society},
volume = {17},
year = {1918},
pages = {75--115}
}
@article{pringsheim1894,
author = {Pringsheim, Alfred},
title = {\"{U}ber Functionen, welche in gewissen Punkten endliche Differentialquotienten jeder endlichen Ordnung, aber keine Taylor'sche Reihenentwickelung besitzen},
journal = {Mathematische Annalen},
volume = {44},
year = {1894},
pages = {41--56}
}
@article{kramerswannier1941,
author = {Kramers, H. A. and Wannier, G. H.},
title = {Statistics of the two-dimensional ferromagnet. Part I},
journal = {Physical Review},
volume = {60},
year = {1941},
pages = {252--262}
}