Fuss–Catalan number

From HandWiki

In combinatorial mathematics and statistics, the Fuss–Catalan numbers are numbers of the form

[math]\displaystyle{ A_m(p,r)\equiv\frac{r}{mp+r}\binom{mp+r}{m} = \frac{r}{m!}\prod_{i=1}^{m-1}(mp+r-i) = r\frac{\Gamma(mp+r)}{\Gamma(1+m)\Gamma(m(p-1)+r+1)}. }[/math]

They are named after N. I. Fuss and Eugène Charles Catalan.

In some publications this equation is sometimes referred to as Two-parameter Fuss–Catalan numbers or Raney numbers. The implication is the single-parameter Fuss-Catalan numbers are when [math]\displaystyle{ \,r=1\, }[/math] and [math]\displaystyle{ \,p=2\, }[/math].


The Fuss-Catalan represents the number of legal permutations or allowed ways of arranging a number of articles, that is restricted in some way. This means that they are related to the Binomial Coefficient. The key difference between Fuss-Catalan and the Binomial Coefficient is that there are no "illegal" arrangement permutations within Binomial Coefficient, but there are within Fuss-Catalan. An example of legal and illegal permutations can be better demonstrated by a specific problem such as balanced brackets (see Dyck language).

A general problem is to count the number of balanced brackets (or legal permutations) that a string of m open and m closed brackets forms (total of 2m brackets). By legally arranged, the following rules apply:

  • For the sequence as a whole, the number of open brackets must equal the number of closed brackets
  • Working along the sequence, the number of open brackets must be greater than the number of closed brackets

As an numeric example how many combinations can 3 pairs of brackets be legally arranged? From the Binomial interpretation there are [math]\displaystyle{ \tbinom {2m}{m} }[/math] or numerically [math]\displaystyle{ \tbinom 63 }[/math] = 20 ways of arranging 3 open and 3 closed brackets. However, there are fewer legal combinations than these when all of the above restrictions apply. Evaluating these by hand, there are 5 legal combinations, namely: ()()(); (())(); ()(()); (()()); ((())). This corresponds to the Fuss-Catalan formula when p=2, r=1 which is the Catalan number formula [math]\displaystyle{ \tfrac{1}{2m+1}\tbinom{2m}{m} }[/math] or [math]\displaystyle{ \tfrac{1}{4}\tbinom63 }[/math]=5. By simple subtraction, there are [math]\displaystyle{ \tfrac{m}{m+1}\tbinom{2m}{m} }[/math] or [math]\displaystyle{ \tfrac34\tbinom63 }[/math] =15 illegal combinations. To further illustrate the subtlety of the problem, if one were to persist with solving the problem just using the Binomial formula, it would be realised that the 2 rules imply that the sequence must start with an open bracket and finish with a closed bracket. This implies that there are [math]\displaystyle{ \tbinom{2m-2}{m-1} }[/math] or [math]\displaystyle{ \tbinom42 }[/math]=6 combinations. This is inconsistent with the above answer of 5, and the missing combination is: ())((), which is illegal and would complete the binomial interpretation.

Whilst the above is a concrete example Catalan numbers, similar problems can be evaluated using Fuss-Catalan formula:

  • Computer Stack: ways of arranging and completing a computer stack of instructions, each time step 1 instruction is processed and p new instructions arrive randomly. If at the beginning of the sequence there are r instructions outstanding.
  • Betting: ways of losing all money when betting. A player has a total stake pot that allows them to make r bets, and plays a game of chance that pays p times the bet stake.
  • Tries: Calculating the number of order m tries on n nodes.[1]

Special Cases

Below is listed a few formulae, along with a few notable special cases

General Form Special Case
[math]\displaystyle{ A_{0}(p,r) = 1 }[/math] [math]\displaystyle{ A_{0}(p,p) = A_{1}(p,1) = 1 }[/math]
[math]\displaystyle{ A_{1}(p,r) = r }[/math] [math]\displaystyle{ A_{1}(p,p) = A_{2}(p,1) = p }[/math]
[math]\displaystyle{ A_{2}(p,r) = r(2p+r-1)/2 }[/math] [math]\displaystyle{ A_{2}(p,p) = A_{3}(p,1) = p(3p-1)/2 }[/math]
[math]\displaystyle{ A_{3}(p,r) = r(3p+r-1)(3p+r-2)/6 }[/math] [math]\displaystyle{ A_{3}(p,p) = A_{4}(p,1) = p(4p-1)(4p-2)/6 }[/math]
[math]\displaystyle{ A_{4}(p,r) = r(4p+r-1)(4p+r-2)(4p+r-3)/24 }[/math] [math]\displaystyle{ A_{4}(p,p) = A_{5}(p,1) = p(5p-1)(5p-2)(5p-3)/24 }[/math]

If [math]\displaystyle{ p=0 }[/math], we recover the Binomial coefficients [math]\displaystyle{ A_m(0,r){{=}}\binom{r}{m} }[/math]

[math]\displaystyle{ A_m(0,1) = 1,1 }[/math];
[math]\displaystyle{ A_m(0,2) = 1,2,1 }[/math];
[math]\displaystyle{ A_m(0,3) = 1,3,3,1 }[/math];
[math]\displaystyle{ A_m(0,4) = 1,4,6,4,1 }[/math].

If [math]\displaystyle{ p=1 }[/math], Pascal's Triangle appears, read along diagonals:

[math]\displaystyle{ A_m(1,1) = 1,1,1,1,1,1,1,1,1,1,\ldots }[/math];
[math]\displaystyle{ A_m(1,2) = 1,2,3,4,5,6,7,8,9,10,\ldots }[/math];
[math]\displaystyle{ A_m(1,3) = 1,3,6,10,15,21,28,35,45,55,\ldots }[/math];
[math]\displaystyle{ A_m(1,4) = 1,4,10,20,35,56,84,120,165,220,\ldots }[/math];
[math]\displaystyle{ A_m(1,5) = 1,5,15,35,70,126,210,330,495,715,\ldots }[/math];
[math]\displaystyle{ A_m(1,6) = 1,6,21,56,126,252,462,792,1287,2002,\ldots }[/math].


For subindex [math]\displaystyle{ m\ge 0 }[/math] the numbers are:

Examples with [math]\displaystyle{ p=2 }[/math]:

[math]\displaystyle{ A_m(2,1) = 1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862,\ldots }[/math] OEISA000108, known as the Catalan Numbers
[math]\displaystyle{ A_m(2,2) = 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, \ldots = A_{m+1}(2,1) }[/math]
[math]\displaystyle{ A_m(2,3) = 1, 3, 9, 28, 90, 297, 1001, 3432, 11934, 41990, \ldots }[/math] OEISA000245
[math]\displaystyle{ A_m(2,4) = 1, 4, 14, 48, 165, 572, 2002, 7072, 25194, 90440,\ldots }[/math] OEISA002057

Examples with [math]\displaystyle{ p=3 }[/math]:

[math]\displaystyle{ A_m(3,1) = 1, 1, 3, 12, 55, 273, 1428, 7752, 43263, 246675, \ldots }[/math] OEISA001764
[math]\displaystyle{ A_m(3,2) = 1, 2, 7, 30, 143, 728, 3876, 21318, 120175, 690690,\ldots }[/math] OEISA006013
[math]\displaystyle{ A_m(3,3) = 1, 3, 12, 55, 273, 1428, 7752, 43263, 246675, 1430715,\ldots = A_{m+1}(3,1) }[/math]
[math]\displaystyle{ A_m(3,4) = 1, 4, 18, 88, 455, 2448, 13566, 76912, 444015, 2601300,\ldots }[/math] OEISA006629

Examples with [math]\displaystyle{ p=4 }[/math]:

[math]\displaystyle{ A_m(4,1) = 1, 1, 4, 22, 140, 969, 7084, 53820, 420732, 3362260,\ldots }[/math] OEISA002293
[math]\displaystyle{ A_m(4,2) = 1, 2, 9, 52, 340, 2394, 17710, 135720, 1068012, 8579560,\ldots }[/math] OEISA069271
[math]\displaystyle{ A_m(4,3) = 1, 3, 15, 91, 612, 4389, 32890, 254475, 2017356, 16301164, \ldots }[/math] OEISA006632
[math]\displaystyle{ A_m(4,4) = 1, 4, 22, 140, 969, 7084, 53820, 420732, 3362260, 27343888,\ldots = A_{m+1}(4,1) }[/math]



[math]\displaystyle{ A_m(p,r)=A_m(p,r-1)+A_{m-1}(p,p+r-1) }[/math] equation (1)

This means in particular that from

[math]\displaystyle{ A_m(p,0)=0 }[/math] equation (2)


[math]\displaystyle{ A_0(p,r)=1 }[/math] equation (3)

one can generate all other Fuss–Catalan numbers if p is an integer.

Riordan (see references) obtains a convolution type of recurrence:

[math]\displaystyle{ A_m(p,s+r) = \sum_{k=0}^m A_k(p,r)A_{m-k}(p,s) }[/math] equation(4)

Generating Function

Paraphrasing the Densities of the Raney distributions [2] paper, let the ordinary generating function with respect to the index m be defined as follows:

[math]\displaystyle{ B_{p,r}(z):=\sum_{m=0}^\infty A_m(p,r)z^m }[/math] equation (5).

Looking at equations (1) and (2), when r=1 it follows that

[math]\displaystyle{ A_m(p,p)=A_{m+1}(p,1) }[/math] equation (6).

Also note this result can be derived by similar substitutions into the other formulas representation, such as the Gamma ratio representation at the top of this article. Using (6) and substituting into (5) an equivalent representation expressed as a generating function can be formulated as

[math]\displaystyle{ B_{p,1}(z) = 1+zB_{p,p}(z) }[/math].

Finally, extending this result by using Lambert's equivalence

[math]\displaystyle{ B_{p,1}(z)^r=B_{p,r}(z) }[/math].

The following result can be derived for the ordinary generating function for all the Fuss-Catalan sequences.

[math]\displaystyle{ B_{p,r}(z) = [1+zB_{p,r}(z)^{p/r}]^r }[/math].

Recursion Representation

Recursion forms of this are as follows: The most obvious form is:

[math]\displaystyle{ A_m(p,r) = \frac{m-1}{m}\frac{\binom{mp+r-1}{m-1}}{\binom{(m-1)p+r-1}{m-2}} A_{m-1}(p,r) }[/math]

Also, a less obvious form is

[math]\displaystyle{ A_m(p,r) = \frac{(m-1)p+r}{m}\frac{\binom{mp+r-1}{p-1}}{\binom{m(p-1)+r}{p-1}} A_{m-1}(p,r) }[/math]

Alternate Representations

In some problems it is easier to use different formula configurations or variations. Below are a two examples using just the binomial function:

[math]\displaystyle{ A_m(p,r)\equiv\frac{r}{mp+r}\binom{mp+r}{m} = \frac{r}{m(p-1)+r}\binom{mp+r-1}{m} = \frac{r}{m}\binom{mp+r-1}{m-1} }[/math]

These variants can be converted into a product, Gamma or Factorial representations too.

See also


  1. Clark, David (1996). "Compact Tries". Compact Pat Trees (Thesis). p. 34.
  2. Mlotkowski, Wojciech; Penson, Karol A.; Zyczkowski, Karol (2013). "Densities of the Raney distributions". Documenta Mathematica 18: 1593–1596. doi:10.4171/dm/437. Bibcode2012arXiv1211.7259M.