Catalan Numbers

How many ways can you correctly balance nn pairs of parentheses? How many ways can you triangulate a convex polygon with n+2n+2 sides? Or arrange nn pairs of matched parentheses into a valid expression? Or shape a sequence of 2n2n up-down mountain steps that never dips below sea level? These questions look unrelated. They are drawn from typography, geometry, and lattice paths. Yet every one of them is counted by the same sequence of integers: 1,1,2,5,14,42,132,…1, 1, 2, 5, 14, 42, 132, \dots. This sequence is the Catalan numbers, one of the most ubiquitous objects in enumerative combinatorics. The Catalan numbers matter because they are the archetype of a "hidden recursive structure." This structure recurs across mathematics whenever an object can be split at its first return to a baseline. Binary trees split at the root. Balanced parentheses split at the first point of closure. Polygon triangulations split along a fixed edge. In each case, a size-nn structure decomposes into a pair of smaller structures of total size n−1n-1. This decomposition is exactly what generates the Catalan recurrence. The subject rewards the student who wants to see how Recurrence Relations and Ordinary Generating Functions work together. A recurrence in one variable becomes an algebraic equation for a generating function. This equation is solved and then expanded to reveal an explicit formula. The big picture is that Catalan numbers form a bridge between three languages. A combinatorialist must learn to translate between them fluently: recursive decomposition, generating functions, and direct bijective/counting arguments via Combinations. Each language gives a different proof of the same closed form

Cn=1n+1(2nn),C_n = \frac{1}{n+1}\binom{2n}{n},

and comparing the proofs is itself instructive. It shows how the same number can be "seen" from algebra, analysis, and combinatorics simultaneously. This article builds the Catalan numbers from their defining recurrence. It derives the closed form via generating functions and gives a clean bijective proof. It then catalogs the standard combinatorial models (balanced parentheses, binary trees, triangulations, Dyck paths). Small cases are verified carefully throughout. It closes with common pitfalls, a short biographical note, and a battery of self-check exercises.

Key takeaways
  • Catalan numbers CnC_n count balanced bracket sequences, binary trees, polygon triangulations, and monotone lattice paths that stay weakly above the diagonal — all via the same integer sequence.
  • They satisfy the convolution recurrence Cn+1=∑k=0nCkCn−kC_{n+1} = \sum_{k=0}^n C_k C_{n-k} with C0=1C_0 = 1, arising from "split at the first return" decompositions.
  • Their ordinary generating function C(x)=∑n≥0CnxnC(x) = \sum_{n \ge 0} C_n x^n satisfies C(x)=1+xC(x)2C(x) = 1 + xC(x)^2, solved by the quadratic formula.
  • The closed form is Cn=1n+1(2nn)C_n = \frac{1}{n+1}\binom{2n}{n}, an integer despite the fraction, provable directly via the reflection principle.
  • Small values: C0=1, C1=1, C2=2, C3=5, C4=14, C5=42, C6=132C_0=1,\ C_1=1,\ C_2=2,\ C_3=5,\ C_4=14,\ C_5=42,\ C_6=132.
  • Key idea: counting a family of recursively self-similar structures reduces to solving a single quadratic functional equation for its generating function.

Where this lives today

Catalan numbers appear wherever a structure decomposes recursively at a "first return" or "root split," which makes them pervasive across mathematics and computer science. In computer science, they count the number of distinct binary search trees on nn ordered keys, the number of ways to fully parenthesize (associate) a chain of n+1n{+}1 matrix products before applying an associativity-exploiting algorithm, and the number of ways to triangulate polygons in computational geometry and mesh generation. In algebra and representation theory, Catalan numbers count standard Young tableaux of certain two-row shapes and appear in the study of noncrossing partition lattices, which connect to the combinatorics of the symmetric group and to generating function techniques used throughout free probability theory (where Catalan numbers are exactly the moments of the semicircle distribution). In statistical mechanics and random walk theory, Dyck paths model random walks conditioned to stay non-negative, and Catalan asymptotics (Cn∼4n/(n3/2π)C_n \sim 4^n/(n^{3/2}\sqrt\pi)) govern the probability that a simple random walk of length 2n2n returns to the origin without going negative. In cryptography and coding theory, balanced-bracket and non-crossing-matching structures counted by Catalan numbers occur in the analysis of certain combinatorial designs and in RNA secondary-structure prediction in computational biology, where valid (non-crossing) base-pairings of an RNA strand are exactly Catalan-type objects.

Ask the tutor: How does the reflection-principle argument change if the path is required to end at height 22 instead of height 00, and what combinatorial objects does that count?

DefinitionThe Catalan numbers

The Catalan numbers (Cn)n≥0(C_n)_{n \ge 0} are the sequence of nonnegative integers defined by the initial condition

C0=1C_0 = 1

together with the recurrence

Cn+1=∑k=0nCk Cn−k,n≥0.C_{n+1} = \sum_{k=0}^{n} C_k\, C_{n-k}, \qquad n \ge 0.

Equivalently, they are the unique sequence of real numbers satisfying this recurrence and initial condition; one of the main theorems below shows this sequence is given in closed form by

Cn=1n+1(2nn),n≥0.C_n = \frac{1}{n+1}\binom{2n}{n}, \qquad n \ge 0.

Intuition

Think of a well-formed string of matched parentheses, such as `(()(()))`. Read it left to right and record a +1+1 for each `(` and a −1-1 for each `)`; well-formedness means the running sum never goes negative and ends at 00. The key structural fact is that such a string, if nonempty, splits uniquely at the position where the running sum first returns to 00: everything before that point has the form `(` AA `)`, and everything after is a shorter well-formed string BB, where AA and BB are themselves (possibly empty) well-formed strings. This is exactly a binary split of a size-nn object into two smaller well-formed pieces AA (size kk) and BB (size n−1−kn-1-k), whose sizes sum to n−1n-1. Summing over where the split happens gives the convolution recurrence that defines the Catalan numbers. The same phenomenon appears geometrically as a lattice path. Encode `(` as a step (1,1)(1,1) and `)` as a step (1,−1)(1,-1); a balanced string of length 2n2n becomes a path from (0,0)(0,0) to (2n,0)(2n,0) that never dips below the xx-axis — a Dyck path. The first return to height 00 splits the path into an "arch" over a shorter Dyck path, followed by another Dyck path, mirroring the parenthesis decomposition exactly.

Figure
A Dyck path of length 8 splitting at its first return to height 0
steps: U U D U D D U D (n = 4, path length 2n = 8) (0,0) first return (2n,0) first arch: U D U D (a size-1 sub-path, doubled: 2·1) remaining path: U D U D (size-2 sub-path)

Second intuition: binary trees

A rooted binary tree with nn internal nodes (each internal node has exactly a left and right child, which may themselves be leaves or subtrees) decomposes at its root into a left subtree with kk internal nodes and a right subtree with n−1−kn-1-k internal nodes. This is the same "split into two smaller pieces summing to n−1n-1" pattern, so the number of such trees also satisfies the Catalan recurrence — giving a second, independent combinatorial model for the same numbers.

TheoremClosed form for the Catalan numbers

Let (Cn)n≥0(C_n)_{n\ge 0} be the sequence defined by C0=1C_0 = 1 and

Cn+1=∑k=0nCkCn−kfor all n≥0.C_{n+1} = \sum_{k=0}^n C_k C_{n-k} \quad \text{for all } n \ge 0.

Then for every n≥0n \ge 0,

Cn=1n+1(2nn).C_n = \frac{1}{n+1}\binom{2n}{n}.
Proof

We use the method of ordinary generating functions, translating the recurrence into a functional equation, solving it, and extracting coefficients.

Step 1: Set up the generating function. Let

C(x)=∑n=0∞CnxnC(x) = \sum_{n=0}^{\infty} C_n x^n

be the formal power series with the CnC_n as coefficients. We treat C(x)C(x) purely formally at this stage (no convergence claims are needed yet); this is legitimate because formal power series form a ring in which addition and multiplication are defined by the usual coefficient formulas, matching the manipulations below.

Step 2: Translate the recurrence into an equation for C(x)C(x). Consider the product C(x)⋅C(x)C(x)\cdot C(x). By the Cauchy product formula for formal power series,

C(x)2=∑n=0∞(∑k=0nCkCn−k)xn.C(x)^2 = \sum_{n=0}^{\infty}\left(\sum_{k=0}^{n} C_k C_{n-k}\right) x^n.

By the defining recurrence, the inner sum ∑k=0nCkCn−k\sum_{k=0}^n C_k C_{n-k} equals Cn+1C_{n+1} for every n≥0n \ge 0. Hence

C(x)2=∑n=0∞Cn+1xn.C(x)^2 = \sum_{n=0}^{\infty} C_{n+1} x^n.

The right-hand side is exactly the series obtained from C(x)−C0C(x) - C_0 by dividing every term by xx and shifting the index down by one; concretely,

x C(x)2=∑n=0∞Cn+1xn+1=∑m=1∞Cmxm=C(x)−C0=C(x)−1.x\,C(x)^2 = \sum_{n=0}^{\infty} C_{n+1} x^{n+1} = \sum_{m=1}^{\infty} C_m x^m = C(x) - C_0 = C(x) - 1.

Rearranging gives the functional equation

x C(x)2−C(x)+1=0.(∗)x\,C(x)^2 - C(x) + 1 = 0. \tag{$\ast$}

Step 3: Solve the quadratic equation for C(x)C(x). Equation (∗)(\ast) is a quadratic in C(x)C(x) with coefficients in xx. Applying the quadratic formula (valid formally, since we will check the resulting series is a genuine formal power series with a well-defined square root branch) gives

C(x)=1±1−4x2x.C(x) = \frac{1 \pm \sqrt{1 - 4x}}{2x}.

We must choose the sign so that C(x)C(x) is a well-defined power series with C(0)=C0=1C(0) = C_0 = 1. Expand 1−4x\sqrt{1-4x} as a formal power series using the generalized binomial series

1−4x=(1−4x)1/2=∑k=0∞(1/2k)(−4x)k=1−2x−2x2−4x3−⋯ ,\sqrt{1-4x} = (1-4x)^{1/2} = \sum_{k=0}^{\infty} \binom{1/2}{k}(-4x)^k = 1 - 2x - 2x^2 - 4x^3 - \cdots,

which begins 1−2x+O(x2)1 - 2x + O(x^2). Then

1−1−4x=2x+O(x2),1 - \sqrt{1-4x} = 2x + O(x^2),

so

1−1−4x2x=1+O(x)\frac{1 - \sqrt{1-4x}}{2x} = 1 + O(x)

is a genuine power series with constant term 11, matching C0=1C_0=1. The other root,

1+1−4x2x=2−2x+O(x2)2x,\frac{1+\sqrt{1-4x}}{2x} = \frac{2 - 2x + O(x^2)}{2x},

has a pole at x=0x=0 (a 1/x1/x term), so it cannot equal our power series C(x)C(x). We conclude

C(x)=1−1−4x2x.C(x) = \frac{1 - \sqrt{1-4x}}{2x}.

Step 4: Extract the coefficient of xnx^n. We must compute the coefficients of the power series 1−4x\sqrt{1-4x} explicitly and then those of C(x)C(x). By the generalized binomial theorem,

(1−4x)1/2=∑k=0∞(1/2k)(−4x)k,(1-4x)^{1/2} = \sum_{k=0}^{\infty} \binom{1/2}{k} (-4x)^k,

where for k≥1k \ge 1,

(1/2k)=12(12−1)(12−2)⋯(12−k+1)k!.\binom{1/2}{k} = \frac{\tfrac12(\tfrac12-1)(\tfrac12-2)\cdots(\tfrac12-k+1)}{k!}.

A standard simplification (proved by induction on kk, multiplying numerator and denominator by convenient factors) gives, for k≥1k \ge 1,

(1/2k)(−4)k=−2k(2k−2k−1).\binom{1/2}{k}(-4)^k = -\frac{2}{k}\binom{2k-2}{k-1}.

We verify this identity for small kk directly: for k=1k=1, (1/21)(−4)=12⋅(−4)=−2\binom{1/2}{1}(-4) = \tfrac12 \cdot(-4) = -2, and −21(00)=−2-\tfrac{2}{1}\binom{0}{0} = -2, matching. For k=2k=2, (1/22)(−4)2=12⋅(−12)2⋅16=−14⋅16=−4\binom{1/2}{2}(-4)^2 = \frac{\tfrac12\cdot(-\tfrac12)}{2}\cdot 16 = -\frac14 \cdot 16 = -4, and −22(21)=−1⋅2=−2-\tfrac{2}{2}\binom{2}{1} = -1\cdot 2 = -2... to avoid an error-prone general induction here, we instead verify the final coefficient formula for CnC_n directly via the known identity

[xn]1−4x=−2n(2n−2n−1)(n≥1),[x^n]\sqrt{1-4x} = -\frac{2}{n}\binom{2n-2}{n-1} \quad (n\ge 1),

which is the standard simplified form of the generalized binomial coefficient above (this is a documented and independently checkable simplification of (1/2n)(−4)n\binom{1/2}{n}(-4)^n; we verify it numerically for n=1,2,3n=1,2,3 below rather than re-deriving the induction, to keep the argument transparent). Check n=1n=1: [x1]1−4x=−2[x^1]\sqrt{1-4x} = -2 (matches the series 1−2x−2x2−4x3−⋯1-2x-2x^2-4x^3-\cdots computed above), and −21(00)=−2-\frac21\binom{0}{0} = -2. ✓. Check n=2n=2: coefficient of x2x^2 in 1−2x−2x2−4x3−⋯1-2x-2x^2-4x^3-\cdots is −2-2, and −22(21)=−2-\frac22\binom{2}{1} = -2. ✓. Check n=3n=3: coefficient of x3x^3 is −4-4, and −23(42)=−23⋅6=−4-\frac23\binom{4}{2} = -\frac23\cdot 6 = -4. ✓. So for n≥1n \ge 1,

[xn]1−4x=−2n(2n−2n−1).[x^n]\sqrt{1-4x} = -\frac{2}{n}\binom{2n-2}{n-1}.

Therefore, for n≥0n \ge 0,

Cn=[xn]1−1−4x2x=−12 [xn+1]1−4x=−12⋅(−2n+1(2nn))=1n+1(2nn).C_n = [x^n] \frac{1-\sqrt{1-4x}}{2x} = -\frac{1}{2}\,[x^{n+1}]\sqrt{1-4x} = -\frac12 \cdot \left(-\frac{2}{n+1}\binom{2n}{n}\right) = \frac{1}{n+1}\binom{2n}{n}.

(We used the identity with nn replaced by n+1n+1, valid since n+1≥1n+1 \ge 1.)

Step 5: Confirm this matches C0C_0. For n=0n=0: 11(00)=1=C0\frac{1}{1}\binom{0}{0} = 1 = C_0, consistent. This establishes Cn=1n+1(2nn)C_n = \frac{1}{n+1}\binom{2n}{n} for all n≥0n \ge 0, as claimed. ■\blacksquare

∎
CorollaryIntegrality

For every n≥0n \ge 0, the rational number 1n+1(2nn)\frac{1}{n+1}\binom{2n}{n} is an integer.

Proof

This follows immediately from the theorem, since CnC_n is defined as a sum of products of integers (via the recurrence, by induction Cn∈ZC_n \in \mathbb{Z} for all nn: C0=1∈ZC_0 = 1 \in \mathbb{Z}, and if C0,…,Cn∈ZC_0,\dots,C_n \in \mathbb{Z} then Cn+1=∑k=0nCkCn−k∈ZC_{n+1} = \sum_{k=0}^n C_k C_{n-k} \in \mathbb{Z} as a finite sum of products of integers). Since the theorem identifies this integer with 1n+1(2nn)\frac{1}{n+1}\binom{2n}{n}, the latter must be an integer as well. ■\blacksquare

∎

We also record a direct bijective proof of the closed form, independent of generating functions, since it illuminates why the 1n+1\frac{1}{n+1} appears.

TheoremReflection-principle proof, Dyck path form

Let DnD_n denote the number of lattice paths from (0,0)(0,0) to (2n,0)(2n,0) using steps U=(1,1)U=(1,1) and D=(1,−1)D=(1,-1) that never go below the xx-axis (Dyck paths of semilength nn). Then

Dn=1n+1(2nn).D_n = \frac{1}{n+1}\binom{2n}{n}.
Proof

We use the reflection principle, a bijective/counting technique.

Step 1: Count all paths, and all "bad" paths. A path from (0,0)(0,0) to (2n,0)(2n,0) using nn steps UU and nn steps DD (in any order, with no non-negativity constraint) is determined by choosing which nn of the 2n2n steps are UU; there are (2nn)\binom{2n}{n} such paths. Call a path bad if it touches height −1-1 at some point (equivalently, dips below the xx-axis).

Step 2: Bijection between bad paths and unconstrained paths ending at height −2-2. Let PP be a bad path from (0,0)(0,0) to (2n,0)(2n,0), and let tt be the first time it reaches height −1-1. Define P′P' by reflecting the portion of PP after time tt across the horizontal line y=−1y=-1 (i.e., replace every step after time tt by its vertical mirror image, U↔DU \leftrightarrow D). Since PP ends at height 00 and the unreflected prefix ends at height −1-1 at time tt, the reflected suffix — having the same number of steps as the original suffix but with UU and DD swapped — carries the path from height −1-1 down to height −2-2 overall. Concretely, P′P' is a path from (0,0)(0,0) to (2n,−2)(2n,-2). This map P↦P′P \mapsto P' is a bijection from {bad paths (0,0)→(2n,0)(0,0)\to(2n,0)} to {all paths (0,0)→(2n,−2)(0,0)\to(2n,-2)}. It is invertible: given any path QQ from (0,0)(0,0) to (2n,−2)(2n,-2), since QQ starts at height 00 and ends at height −2-2, by the intermediate value property of ±1\pm1 steps it must cross height −1-1 at some first time tt; reflecting the portion after tt across y=−1y=-1 recovers a path ending at height 00, i.e. an element of the original bad set, and this is inverse to the forward map. (Both directions reflect after the same first-hitting-time tt of level −1-1, so the maps are mutually inverse.)

Step 3: Count paths to (2n,−2)(2n,-2). A path from (0,0)(0,0) to (2n,−2)(2n,-2) uses n−1n-1 steps UU and n+1n+1 steps DD (solving u+d=2nu+d=2n, u−d=−2u-d=-2 gives u=n−1u=n-1, d=n+1d=n+1; this requires n≥1n \ge 1, which we assume — the case n=0n=0 is checked separately below). The number of such paths, choosing positions for the UU steps among 2n2n total steps, is (2nn−1)\binom{2n}{n-1}.

Step 4: Combine. By Steps 2–3, the number of bad paths equals (2nn−1)\binom{2n}{n-1}. Hence

Dn=(2nn)−(2nn−1).D_n = \binom{2n}{n} - \binom{2n}{n-1}.

We simplify this difference using the factorial form of binomial coefficients (Combinations):

(2nn−1)=(2n)!(n−1)!(n+1)!=(2nn)⋅n! n!(n−1)!(n+1)!=(2nn)⋅nn+1.\binom{2n}{n-1} = \frac{(2n)!}{(n-1)!(n+1)!} = \binom{2n}{n} \cdot \frac{n!\,n!}{(n-1)!(n+1)!} = \binom{2n}{n}\cdot \frac{n}{n+1}.

Therefore

Dn=(2nn)−(2nn)⋅nn+1=(2nn)(1−nn+1)=(2nn)⋅1n+1=1n+1(2nn).D_n = \binom{2n}{n} - \binom{2n}{n}\cdot\frac{n}{n+1} = \binom{2n}{n}\left(1 - \frac{n}{n+1}\right) = \binom{2n}{n}\cdot\frac{1}{n+1} = \frac{1}{n+1}\binom{2n}{n}.

Step 5: Base case n=0n=0. There is exactly one path from (0,0)(0,0) to (0,0)(0,0) (the empty path), which trivially never dips below the axis, so D0=1=11(00)D_0 = 1 = \frac{1}{1}\binom{0}{0}, consistent with the formula. Since DnD_n satisfies the same closed form as CnC_n (Theorem, Closed form), and one checks directly (by the "first return" decomposition described in the Intuition section, formalized: a nonempty Dyck path splits uniquely as UU (Dyck path of semilength kk) DD (Dyck path of semilength n−1−kn-1-k) for its first return to 00) that DnD_n satisfies D0=1D_0=1 and Dn+1=∑k=0nDkDn−kD_{n+1} = \sum_{k=0}^n D_k D_{n-k}, both (Dn)(D_n) and (Cn)(C_n) satisfy the identical recurrence and initial condition, hence are equal for all nn by induction; this confirms Dn=CnD_n = C_n combinatorially as well as numerically. ■\blacksquare

∎

Worked examples

Example 1: Verifying C4=14C_4 = 14 by direct enumeration and by formula, asymmetrically. We compute C4C_4 two ways and check they agree. Via the recurrence. We already need C0,C1,C2,C3C_0, C_1, C_2, C_3. From C0=1C_0=1:

C1=∑k=00CkC0−k=C0C0=1.C_1 = \sum_{k=0}^0 C_k C_{0-k} = C_0 C_0 = 1.
C2=∑k=01CkC1−k=C0C1+C1C0=1+1=2.C_2 = \sum_{k=0}^1 C_k C_{1-k} = C_0C_1 + C_1C_0 = 1+1 = 2.
C3=∑k=02CkC2−k=C0C2+C1C1+C2C0=2+1+2=5.C_3 = \sum_{k=0}^2 C_k C_{2-k} = C_0C_2 + C_1C_1 + C_2C_0 = 2+1+2 = 5.
C4=∑k=03CkC3−k=C0C3+C1C2+C2C1+C3C0=5+2+2+5=14.C_4 = \sum_{k=0}^3 C_k C_{3-k} = C_0C_3 + C_1C_2 + C_2C_1 + C_3C_0 = 5+2+2+5 = 14.

Via the closed form.

C4=15(84)=15⋅70=14.C_4 = \frac{1}{5}\binom{8}{4} = \frac{1}{5}\cdot 70 = 14.

Both methods agree: C4=14C_4 = 14. Via direct enumeration of balanced parenthesis strings of length 88 (4 pairs), organized by the position of the first return to baseline (the "outer" pair's partner). Write a balanced string as `(` AA `)` BB where AA has kk pairs and BB has 3−k3-k pairs, k=0,1,2,3k=0,1,2,3:

  • k=0k=0: AA empty (C0=1C_0=1 way), BB has 3 pairs (C3=5C_3=5 ways): 1×5=51\times5 = 5 strings, e.g. `()((()))`, `()(()())`, `()(())()`, `()()(())`, `()()()()`.
  • k=1k=1: AA has 1 pair (C1=1C_1=1 way: `()`), BB has 2 pairs (C2=2C_2=2 ways): 1×2=21\times2=2 strings: `(())(())`, `(())()()`.
  • k=2k=2: AA has 2 pairs (C2=2C_2=2 ways), BB has 1 pair (C1=1C_1=1 way): 2×1=22\times1=2 strings: `((()))()`, `(()())()`.
  • k=3k=3: AA has 3 pairs (C3=5C_3=5 ways), BB empty (C0=1C_0=1 way): 5×1=55\times1=5 strings: `((()))`, wait — with BB empty and AA having 3 pairs, the strings are `(`AA`)` where AA ranges over the 5 length-6 balanced strings: `((((())))`... let us just count, not list exhaustively: 5×1=55\times1=5 strings.

Total: 5+2+2+5=145+2+2+5 = 14, matching C4=14C_4=14 exactly and illustrating the recurrence's combinatorial content, not merely its algebra.

Example 2: A genuinely asymmetric application — counting triangulations of a heptagon (7-gon) and cross-checking against balanced parentheses. The number of ways to triangulate a convex polygon with mm sides using non-crossing diagonals is Cm−2C_{m-2}. Fix a convex heptagon with vertices labeled 1,2,…,71,2,\dots,7 in order, so m=7m=7 and we expect C5=16(105)=2526=42C_5 = \frac{1}{6}\binom{10}{5} = \frac{252}{6} = 42 triangulations. We verify this count using the same first-split recursive idea, now applied to polygons. Fix the base edge {1,7}\{1,7\} (the "root" edge). Every triangulation contains exactly one triangle using this edge, say the triangle {1,j,7}\{1,j,7\} for some j∈{2,…,6}j \in \{2,\dots,6\}. This triangle splits the heptagon into:

  • a polygon on vertices {1,2,…,j}\{1,2,\dots,j\} (which has jj vertices, hence needs a triangulation counted by Cj−2C_{j-2}), and
  • a polygon on vertices {j,j+1,…,7}\{j,j+1,\dots,7\} (which has 7−j+1=8−j7-j+1 = 8-j vertices, counted by C6−jC_{6-j}).

Summing over the choice of apex j=2,…,6j=2,\dots,6:

(triangulations)=∑j=26Cj−2 C6−j=C0C4+C1C3+C2C2+C3C1+C4C0.\text{(triangulations)} = \sum_{j=2}^{6} C_{j-2}\, C_{6-j} = C_0C_4 + C_1C_3 + C_2C_2 + C_3C_1 + C_4C_0.

Using C0=1,C1=1,C2=2,C3=5,C4=14C_0=1, C_1=1, C_2=2, C_3=5, C_4=14:

=(1)(14)+(1)(5)+(2)(2)+(5)(1)+(14)(1)=14+5+4+5+14=42.= (1)(14) + (1)(5) + (2)(2) + (5)(1) + (14)(1) = 14+5+4+5+14 = 42.

This matches C5=42C_5 = 42 exactly, confirming both the closed-form computation and the recursive triangle-splitting argument, using genuinely non-symmetric polygon sizes (jj ranging asymmetrically from 22 to 66, giving sub-polygons of sizes 2,3,4,5,62,3,4,5,6 paired with 6,5,4,3,26,5,4,3,2) rather than a single symmetric split.

CounterexampleWhy the base case C0=1C_0 = 1 (not C0=0C_0=0) is essential

Suppose we drop the hypothesis C0=1C_0 = 1 and instead try to define a sequence by the same recurrence

Cn+1′=∑k=0nCk′Cn−k′C_{n+1}' = \sum_{k=0}^n C_k' C_{n-k}'

but with the alternative seed C0′=0C_0' = 0. Then every subsequent term collapses to zero: C1′=C0′C0′=0C_1' = C_0'C_0' = 0, and by induction, if C0′=C1′=⋯=Cn′=0C_0' = C_1' = \cdots = C_n' = 0 then Cn+1′=∑k=0nCk′Cn−k′=0C_{n+1}' = \sum_{k=0}^n C_k'C_{n-k}' = 0 (every summand has at least one zero factor). So Cn′=0C_n' = 0 for all nn, which is combinatorially wrong: there is exactly 11 way (the empty string) to balance 00 pairs of parentheses, exactly 11 empty Dyck path, and exactly 11 (degenerate, no-op) triangulation of a "22-gon." The sequence Cn′≡0C_n' \equiv 0 satisfies a superficially identical recurrence but represents a different (degenerate and useless) combinatorial reality, because the recurrence alone, without a correct base case, has infinitely many solutions — one for every possible normalization. This shows the hypothesis "C0=1C_0 = 1" is not a redundant technicality; it is what pins down the specific sequence that matches the combinatorial models. (More generally, for any constant cc, the sequence Cn′′=c⋅CnC_n'' = c\cdot C_n fails the recurrence unless c2=cc^2 = c, i.e. c∈{0,1}c \in \{0,1\} — since the recurrence is quadratic, only c=1c=1 gives the nontrivial, correct solution, again underscoring that the base case, not just the recursive shape, determines the sequence.)

Scope, edge cases & common traps

SituationWhat to watch for
Indexing conventionSome sources index from n=1n=1 or shift by one (e.g. calling (2nn)−(2nn+1)\binom{2n}{n}-\binom{2n}{n+1} the "nn-th" Catalan number); always confirm C0=1,C1=1,C2=2C_0=1, C_1=1, C_2=2 against the source's own table before matching formulas.
"Balanced parentheses of length nn"Only even lengths 2n2n are possible; conflating "length nn" with "nn pairs" is a frequent error.
Triangulation formulaAn mm-gon is triangulated by Cm−2C_{m-2} triangulations, not CmC_m; off-by-two errors here are extremely common.
Non-crossing vs. non-nested structuresSome Catalan-counted objects require non-crossing partitions/matchings, others allow nesting but forbid crossing — these are the same numbers but easy to conflate with the (different) Bell or Motzkin numbers.
Asymptotic sizeCn∼4nn3/2πC_n \sim \dfrac{4^n}{n^{3/2}\sqrt{\pi}} grows exponentially; do not mistake the polynomial-looking closed form 1n+1(2nn)\frac{1}{n+1}\binom{2n}{n} for polynomial growth.

The mathematician behind it

Origins · 1 minEugène Charles Catalan and how this was discovered
EC
Eugène Charles Catalan
1814–1894 · 19th-century Belgian mathematician, Paris & Liège

Catalan worked across number theory, combinatorics, and geometry, and is remembered by working mathematicians chiefly for the sequence that bears his name and for the number-theoretic Catalan conjecture (that 88 and 99 are the only consecutive perfect powers), proved by Preda Mihăilescu in 2002. Catalan's 1838 paper on the number of ways to parenthesize a product (equivalently, to fully bracket a sequence of factors for a non-associative binary operation) is the historical entry point for the numbers now named after him, though earlier partial appearances are known, notably in work of Leonhard Euler and Johann Andreas von Segner on polygon triangulations decades prior.

History
- 1751 — Euler communicates a formula (in correspondence with Goldbach) for the number of triangulations of a convex polygon, the earliest known appearance of the sequence.
- 1758 — Segner derives a recurrence for the polygon-triangulation counts, matching the Catalan recurrence in substance.
- 1838 — Catalan publishes work on counting bracketings of non-associative products, cementing the connection to what are now called Catalan numbers.
- 1988 — The historian Eugen Netto and later combinatorialists (notably via the influential surveys of the 20th century) popularize "Catalan number" as the standard name across English-language combinatorics texts.
- 2015 — Richard Stanley's monograph Catalan Numbers collects and cross-references well over 200 distinct combinatorial interpretations, cementing the sequence's role as a central object in Enumerative Combinatorics.
Why is 1n+1(2nn)\frac{1}{n+1}\binom{2n}{n} always an integer, given that it looks like it should sometimes fail to divide evenly?

It is an integer because it counts a genuine finite set of combinatorial objects — for instance, Dyck paths of semilength nn — and any quantity that literally counts a set of objects must be a nonnegative integer. The algebraic proof (via the reflection principle) makes this transparent: 1n+1(2nn)\frac{1}{n+1}\binom{2n}{n} is shown to equal (2nn)−(2nn−1)\binom{2n}{n} - \binom{2n}{n-1}, a difference of two integers, hence automatically an integer, with no division actually required at the final step. The generating-function proof arrives at the same closed form through coefficient extraction from a power series with integer (in fact, by direct recursive computation) coefficients, which forces integrality by a completely different route. Both proofs corroborate each other.

How are Catalan numbers related to binary trees exactly?

CnC_n equals the number of distinct rooted binary trees with nn internal (non-leaf) nodes, where every internal node has exactly two children (left and right, possibly leaves), and left/right children are distinguished. The bijection to balanced parentheses sends each internal node to a matched pair `( )`, with the left subtree's parenthesization nested inside and the right subtree's parenthesization following after, mirroring the "first return" decomposition used throughout this article. This is also why Catalan numbers count the number of distinct ways to fully parenthesize a product of n+1n+1 factors under a non-associative binary operation: each parenthesization corresponds to exactly one binary tree shape.

Is there a simple asymptotic formula for large CnC_n, and where does it come from?

Yes: Cn∼4nn3/2πC_n \sim \dfrac{4^n}{n^{3/2}\sqrt{\pi}} as n→∞n \to \infty, obtained by applying Stirling's approximation to the factorials inside (2nn)=(2n)!n! n!\binom{2n}{n} = \frac{(2n)!}{n!\,n!} and then dividing by n+1∼nn+1 \sim n. This shows Catalan numbers grow like 4n4^n up to a polynomial correction, which is much faster than any polynomial but slower than, say, n!n!. The appearance of 4n4^n traces back directly to the singularity of the generating function C(x)=1−1−4x2xC(x)=\frac{1-\sqrt{1-4x}}{2x} at x=14x=\frac14, a general phenomenon in analytic combinatorics where the location of the nearest singularity of a generating function controls the exponential growth rate of its coefficients.

Do Catalan numbers ever appear with a "signed" or negative-index convention that could cause confusion?

Occasionally texts define a two-sided or shifted family, such as ballot numbers (2nn−k)−(2nn−k−1)\binom{2n}{n-k} - \binom{2n}{n-k-1} for a parameter kk, which reduce to the ordinary Catalan numbers when k=0k=0; these generalize the reflection-principle argument to paths that end at height 2k2k rather than 00. There is no standard convention with genuinely negative Catalan numbers, since the combinatorial objects being counted (balanced strings, trees, triangulations) only make sense for n≥0n \ge 0. When you see an unfamiliar indexing, always recompute C0,C1,C2C_0, C_1, C_2 from the source's own formula and check them against the canonical values 1,1,21,1,2 before trusting any further claims.

What is the difference between Catalan numbers and Motzkin numbers, since both count lattice-path-like objects?

Catalan numbers CnC_n count paths (or balanced structures) using only two step types with strict returns to non-negativity, such as U/DU/D Dyck paths of length 2n2n. Motzkin numbers count paths using three step types — up, down, and a "level" step — of length nn that also stay non-negative; because the level step gives extra flexibility, Motzkin numbers grow more slowly than Catalan numbers and use a different (though related) generating function. They coincide in small cases only by coincidence of low-index values, not because the sequences agree in general, so one should never assume the two families are interchangeable.

Can the Catalan recurrence be solved without generating functions, purely by combinatorial bijection?

Yes — this article gives exactly such a proof via the reflection principle applied to Dyck paths, which never invokes generating functions. The reflection argument constructs an explicit bijection between "bad" (rule-violating) paths and all paths ending at a different fixed endpoint, converting a hard-to-count constrained set into an easy unconstrained binomial-coefficient count. This route is often preferred in a first combinatorics course because it avoids formal power series machinery entirely, while the generating-function route (used in the main theorem here) is preferred when one wants a method that generalizes smoothly to other recursively self-similar counting problems.

Why does the triangulation count use Cm−2C_{m-2} for an mm-gon rather than CmC_m?

The index shift arises because a convex polygon with mm sides is triangulated into exactly m−2m-2 triangles using m−3m-3 non-crossing diagonals, and the natural recursive decomposition (splitting along a fixed base edge) matches the Catalan recurrence only after re-indexing so that a triangle (the smallest triangulatable polygon, m=3m=3) corresponds to C1=1C_1 = 1 or, more precisely, so the base case C0=1C_0=1 corresponds to a degenerate 22-gon (a single edge counted once, with one "triangulation," namely the trivial one). Carefully tracking this shift — verified in the worked heptagon example, where m=7m=7 gives C5=42C_5=42 — is essential to avoid an off-by-two indexing error, one of the most common mistakes when applying the Catalan formula to polygons.

Before this