|
%I #23 Jun 15 2022 11:46:34
%S 1,2,36,1200,58800,3810240,307359360,29682132480,3339239904000,
%T 428906814336000,61934143990118400,9931984545324441600,
%U 1751339941492209868800,336796142594655744000000,70149825129001153536000000,15732267448930658699673600000
%N Central Lah numbers: a(n) = A105278(2*n,n) = A008297(2*n,n).
%C a(n) is the number of Lah partitions of a set of size 2n with n blocks.
%F a(n) = binomial(2n-1,n-1)*(2n)!/n! (for n>0).
%F D-finite with recurrence (n+1)*a(n+1) = 4*(2n+1)^2*a(n) - delta(n,0).
%F a(n) ~ 2^(4*n)*n^n*exp(-n)/sqrt(2*n*Pi).
%F a(n)*a(n+2) - a(n+1)^2 is >= 0 and is a multiple of 2^(n+3) for all nonnegative n.
%F a(n) == 0 (mod 10) for n>3.
%F E.g.f.: 1/2 + K(16x)/Pi, where K(z) is the complete elliptic integral of the first kind, which can also be written as a Legendre function of the second kind.
%F a(n) = Catalan(n)*C(2*n-1,n)*(n+1)!. - _Peter Luschny_, Oct 07 2014
%F a(n) = A125558(n)*(n+1)! = A090181(2*n,n)*(n+1)!. - _Peter Luschny_, Oct 07 2014
%F a(n) = (2/n)*(Gamma(2*n)^2/Gamma(n)^3) for n>0. - _Peter Luschny_, Oct 17 2014
%p A187535:= n -> if n=0 then 1 else binomial(2*n-1,n-1)*(2*n)!/n! fi;
%p seq(A187535(n),n=0..12);
%t a[n_]:=If[n==0,1,Binomial[2n-1,n-1](2n)!/n!]
%t Table[a[n],{n,0,12}]
%t (* Alternative: *)
%t a[n_] := Binomial[2*n, n] FactorialPower[2*n - 1, n];
%t Table[a[n], {n, 0, 15}] (* _Peter Luschny_, Jun 15 2022 *)
%o (Maxima) a(n) := if n=0 then 1 else binomial(2*n-1,n-1)*(2*n)!/n!;
%o makelist(a(n),n,0,12);
%o (Sage)
%o [catalan_number(n)*binomial(2*n-1,n)*factorial(n+1) for n in range(15)] # _Peter Luschny_, Oct 07 2014
%Y Cf. A008297, A111596, A066667, A187536, A187538, A187539, A187540, A187542 - A187548, A090181, A125558.
%K nonn,easy
%O 0,2
%A _Emanuele Munarini_, Mar 11 2011
|
