Fuss–Catalan number

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

A m ( p , r ) r m p + r ( m p + r m ) = r m ! i = 1 m 1 ( m p + r i ) = r Γ ( m p + r ) Γ ( 1 + m ) Γ ( m ( p 1 ) + r + 1 ) . {\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)}}.}

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 r = 1 {\displaystyle \,r=1\,} and p = 2 {\displaystyle \,p=2\,} .

Uses

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 ( 2 m m ) {\displaystyle {\tbinom {2m}{m}}} or numerically ( 6 3 ) {\displaystyle {\tbinom {6}{3}}} = 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 1 2 m + 1 ( 2 m m ) {\displaystyle {\tfrac {1}{2m+1}}{\tbinom {2m}{m}}} or 1 4 ( 6 3 ) {\displaystyle {\tfrac {1}{4}}{\tbinom {6}{3}}} =5. By simple subtraction, there are m m + 1 ( 2 m m ) {\displaystyle {\tfrac {m}{m+1}}{\tbinom {2m}{m}}} or 3 4 ( 6 3 ) {\displaystyle {\tfrac {3}{4}}{\tbinom {6}{3}}} =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 ( 2 m 2 m 1 ) {\displaystyle {\tbinom {2m-2}{m-1}}} or ( 4 2 ) {\displaystyle {\tbinom {4}{2}}} =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
A 0 ( p , r ) = 1 {\displaystyle A_{0}(p,r)=1} A 0 ( p , p ) = A 1 ( p , 1 ) = 1 {\displaystyle A_{0}(p,p)=A_{1}(p,1)=1}
A 1 ( p , r ) = r {\displaystyle A_{1}(p,r)=r} A 1 ( p , p ) = A 2 ( p , 1 ) = p {\displaystyle A_{1}(p,p)=A_{2}(p,1)=p}
A 2 ( p , r ) = r ( 2 p + r 1 ) / 2 {\displaystyle A_{2}(p,r)=r(2p+r-1)/2} A 2 ( p , p ) = A 3 ( p , 1 ) = p ( 3 p 1 ) / 2 {\displaystyle A_{2}(p,p)=A_{3}(p,1)=p(3p-1)/2}
A 3 ( p , r ) = r ( 3 p + r 1 ) ( 3 p + r 2 ) / 6 {\displaystyle A_{3}(p,r)=r(3p+r-1)(3p+r-2)/6} A 3 ( p , p ) = A 4 ( p , 1 ) = p ( 4 p 1 ) ( 4 p 2 ) / 6 {\displaystyle A_{3}(p,p)=A_{4}(p,1)=p(4p-1)(4p-2)/6}
A 4 ( p , r ) = r ( 4 p + r 1 ) ( 4 p + r 2 ) ( 4 p + r 3 ) / 24 {\displaystyle A_{4}(p,r)=r(4p+r-1)(4p+r-2)(4p+r-3)/24} A 4 ( p , p ) = A 5 ( p , 1 ) = p ( 5 p 1 ) ( 5 p 2 ) ( 5 p 3 ) / 24 {\displaystyle A_{4}(p,p)=A_{5}(p,1)=p(5p-1)(5p-2)(5p-3)/24}


If p = 0 {\displaystyle p=0} , we recover the Binomial coefficients A m ( 0 , r ) = ( r m ) {\displaystyle A_{m}(0,r){=}{\binom {r}{m}}}

A m ( 0 , 1 ) = 1 , 1 {\displaystyle A_{m}(0,1)=1,1} ;
A m ( 0 , 2 ) = 1 , 2 , 1 {\displaystyle A_{m}(0,2)=1,2,1} ;
A m ( 0 , 3 ) = 1 , 3 , 3 , 1 {\displaystyle A_{m}(0,3)=1,3,3,1} ;
A m ( 0 , 4 ) = 1 , 4 , 6 , 4 , 1 {\displaystyle A_{m}(0,4)=1,4,6,4,1} .

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

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

Examples

For subindex m 0 {\displaystyle m\geq 0} the numbers are:

Examples with p = 2 {\displaystyle p=2} :

A m ( 2 , 1 ) = 1 , 1 , 2 , 5 , 14 , 42 , 132 , 429 , 1430 , 4862 , {\displaystyle A_{m}(2,1)=1,1,2,5,14,42,132,429,1430,4862,\ldots } OEIS: A000108, known as the Catalan Numbers
A m ( 2 , 2 ) = 1 , 2 , 5 , 14 , 42 , 132 , 429 , 1430 , 4862 , 16796 , = A m + 1 ( 2 , 1 ) {\displaystyle A_{m}(2,2)=1,2,5,14,42,132,429,1430,4862,16796,\ldots =A_{m+1}(2,1)}
A m ( 2 , 3 ) = 1 , 3 , 9 , 28 , 90 , 297 , 1001 , 3432 , 11934 , 41990 , {\displaystyle A_{m}(2,3)=1,3,9,28,90,297,1001,3432,11934,41990,\ldots } OEIS: A000245
A m ( 2 , 4 ) = 1 , 4 , 14 , 48 , 165 , 572 , 2002 , 7072 , 25194 , 90440 , {\displaystyle A_{m}(2,4)=1,4,14,48,165,572,2002,7072,25194,90440,\ldots } OEIS: A002057

Examples with p = 3 {\displaystyle p=3} :

A m ( 3 , 1 ) = 1 , 1 , 3 , 12 , 55 , 273 , 1428 , 7752 , 43263 , 246675 , {\displaystyle A_{m}(3,1)=1,1,3,12,55,273,1428,7752,43263,246675,\ldots } OEIS: A001764
A m ( 3 , 2 ) = 1 , 2 , 7 , 30 , 143 , 728 , 3876 , 21318 , 120175 , 690690 , {\displaystyle A_{m}(3,2)=1,2,7,30,143,728,3876,21318,120175,690690,\ldots } OEIS: A006013
A m ( 3 , 3 ) = 1 , 3 , 12 , 55 , 273 , 1428 , 7752 , 43263 , 246675 , 1430715 , = A m + 1 ( 3 , 1 ) {\displaystyle A_{m}(3,3)=1,3,12,55,273,1428,7752,43263,246675,1430715,\ldots =A_{m+1}(3,1)}
A m ( 3 , 4 ) = 1 , 4 , 18 , 88 , 455 , 2448 , 13566 , 76912 , 444015 , 2601300 , {\displaystyle A_{m}(3,4)=1,4,18,88,455,2448,13566,76912,444015,2601300,\ldots } OEIS: A006629

Examples with p = 4 {\displaystyle p=4} :

A m ( 4 , 1 ) = 1 , 1 , 4 , 22 , 140 , 969 , 7084 , 53820 , 420732 , 3362260 , {\displaystyle A_{m}(4,1)=1,1,4,22,140,969,7084,53820,420732,3362260,\ldots } OEIS: A002293
A m ( 4 , 2 ) = 1 , 2 , 9 , 52 , 340 , 2394 , 17710 , 135720 , 1068012 , 8579560 , {\displaystyle A_{m}(4,2)=1,2,9,52,340,2394,17710,135720,1068012,8579560,\ldots } OEIS: A069271
A m ( 4 , 3 ) = 1 , 3 , 15 , 91 , 612 , 4389 , 32890 , 254475 , 2017356 , 16301164 , {\displaystyle A_{m}(4,3)=1,3,15,91,612,4389,32890,254475,2017356,16301164,\ldots } OEIS: A006632
A m ( 4 , 4 ) = 1 , 4 , 22 , 140 , 969 , 7084 , 53820 , 420732 , 3362260 , 27343888 , = A m + 1 ( 4 , 1 ) {\displaystyle A_{m}(4,4)=1,4,22,140,969,7084,53820,420732,3362260,27343888,\ldots =A_{m+1}(4,1)}

Algebra

Recurrence

A m ( p , r ) = A m ( p , r 1 ) + A m 1 ( p , p + r 1 ) {\displaystyle A_{m}(p,r)=A_{m}(p,r-1)+A_{m-1}(p,p+r-1)} equation (1)

This means in particular that from

A m ( p , 0 ) = 0 {\displaystyle A_{m}(p,0)=0} equation (2)

and

A 0 ( p , r ) = 1 {\displaystyle A_{0}(p,r)=1} equation (3)

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

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

A m ( p , s + r ) = k = 0 m A k ( p , r ) A m k ( p , s ) {\displaystyle A_{m}(p,s+r)=\sum _{k=0}^{m}A_{k}(p,r)A_{m-k}(p,s)} 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:

B p , r ( z ) := m = 0 A m ( p , r ) z m {\displaystyle B_{p,r}(z):=\sum _{m=0}^{\infty }A_{m}(p,r)z^{m}} equation (5).

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

A m ( p , p ) = A m + 1 ( p , 1 ) {\displaystyle A_{m}(p,p)=A_{m+1}(p,1)} 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

B p , 1 ( z ) = 1 + z B p , p ( z ) {\displaystyle B_{p,1}(z)=1+zB_{p,p}(z)} .

Finally, extending this result by using Lambert's equivalence

B p , 1 ( z ) r = B p , r ( z ) {\displaystyle B_{p,1}(z)^{r}=B_{p,r}(z)} .

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

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

Recursion Representation

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

A m ( p , r ) = m 1 m ( m p + r 1 m 1 ) ( ( m 1 ) p + r 1 m 2 ) A m 1 ( p , r ) {\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)}

Also, a less obvious form is

A m ( p , r ) = ( m 1 ) p + r m ( m p + r 1 p 1 ) ( m ( p 1 ) + r p 1 ) A m 1 ( p , r ) {\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)}

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:

A m ( p , r ) r m p + r ( m p + r m ) = r m ( p 1 ) + r ( m p + r 1 m ) = r m ( m p + r 1 m 1 ) {\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}}}

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

See also

References

  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. arXiv:1211.7259. Bibcode:2012arXiv1211.7259M. doi:10.4171/dm/437. S2CID 17493895.
  • Fuss, N. I. (1791). "Solutio quaestionis, quot modis polygonum n laterum in polygona m laterum, per diagonales resolvi queat". Nova Acta Academiae Scientiarum Imperialis Petropolitanae. 9: 243–251.
  • Riordan, J. (1968). Combinatorial identities. Wiley. ISBN 978-0471722755.
  • Bisch, Dietmar; Jones, Vaughan (1997). "Algebras associated to intermediate subfactors". Inventiones Mathematicae. 128 (1): 89–157. Bibcode:1997InMat.128...89J. doi:10.1007/s002220050137. S2CID 119372640.
  • Przytycki, Jozef H.; Sikora, Adam S. (2000). "Polygon dissections and Euler, Fuss, Kirkman, and Cayley Numbers". Journal of Combinatorial Theory, Series A. 92: 68–76. arXiv:math/9811086. doi:10.1006/jcta.1999.3042. S2CID 15724561.
  • Aval, Jean-Christophe (2008). "Multivariate Fuss-Catalan numbers". Discrete Mathematics. 20 (308): 4660–4669. arXiv:0711.0906. doi:10.1016/j.disc.2007.08.100. S2CID 509695.
  • Alexeev, N.; Götze, F; Tikhomirov, A. (2010). "Asymptotic distribution of singular values of powers of random matrices". Lithuanian Mathematical Journal. 50 (2): 121–132. arXiv:1012.2743. doi:10.1007/s10986-010-9074-4. S2CID 3799186.
  • Mlotkowski, Wojciech (2010). "Fuss-Catalan Numbers in noncommutative probability". Documenta Mathematica. 15: 939–955. doi:10.4171/dm/318. S2CID 8083529.
  • Penson, Karol A.; Zyczkowski, Karol (2011). "Product of Ginibre matrices: Fuss-Catalan and Raney distributions". Physical Review E. 83 (6): 061118. arXiv:1103.3453. Bibcode:2011PhRvE..83f1118P. doi:10.1103/PhysRevE.83.061118. PMID 21797313. S2CID 15289531.
  • Gordon, Ian G.; Griffeth, Stephen (2012). "Catalan numbers for complex reflection groups". American Journal of Mathematics. 134 (6): 1491–1502. arXiv:0912.1578. doi:10.1353/ajm.2012.0047. S2CID 2654664.
  • Mlotkowski, W.; Penson, K. A. (2015). "A Fuss-type family of positive definite sequences". arXiv:1507.07312 [math.PR].
  • He, Tian-Xiao; Shapiro, Louis W. (2017). "Fuss-Catalan matrices, their weighted sums, and stabilizer subgroups of the Riordan group". Linear Algebra and Its Applications. 532: 25–42. doi:10.1016/j.laa.2017.06.025.
  • v
  • t
  • e
Classes of natural numbers
Of the form a × 2b ± 1
Other polynomial numbers
Recursively defined numbers
Possessing a specific set of other numbers
Expressible via specific sums
2-dimensional
centered
non-centered
3-dimensional
centered
non-centered
pyramidal
4-dimensional
non-centered
Combinatorial numbers
Divisor functions
Prime omega functions
Euler's totient function
Aliquot sequences
Primorial
Numeral system-dependent numbers
Arithmetic functions
and dynamics
Digit sum
Digit product
Coding-related
Other
P-adic numbers-related
Digit-composition related
Digit-permutation related
Divisor-related
Other
Generated via a sieve
  • Mathematics portal