login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A152059
a(n) is the number of ways 2n-1 seats can be occupied by at most n people for n>=1, with a(0)=1.
3
1, 2, 13, 136, 1961, 36046, 805597, 21204548, 642451441, 22021483546, 842527453421, 35591363004352, 1645373927307673, 82625931422081126, 4478815087922020861, 260648364396903639676, 16208855884741850686817
OFFSET
0,2
COMMENTS
Let A(x) be the e.g.f. of this sequence, and B(x) be the e.g.f. of A082545, then B(x)/A(x) = C(x) where C(x) = 1 + x*C(x)^2 is the Catalan function (A000108). This follows from the fact that this sequence and A082545 form adjacent semi-diagonals of table A088699. - Paul D. Hanna, Aug 16 2022
FORMULA
a(n) = Sum_{k=0..n} k! * C(2*n-1,k) * C(n,k).
Central terms of triangle A086885 (after initial term).
a(n) = n! * [x^n] exp(x/(1 - x))/(1 - x)^n. - Ilya Gutkovskiy, Oct 02 2017
a(n) ~ 2^(2*n - 1/2) * n^n / exp(n-1). - Vaclav Kotesovec, Oct 02 2017
a(n) = n! * pollaguerre(n, n-1, -1). - Seiichi Manyama, May 01 2021
From Paul D. Hanna, Aug 16 2022: (Start)
E.g.f.: exp( (1-2*x - sqrt(1-4*x))/(2*x) ) / ((sqrt(1-4*x) - (1-4*x))/(2*x)), derived from the e.g.f for A082545 given by Mark van Hoeij.
E.g.f.: exp(C(x) - 1) / (2 - C(x)), where C(x) = (1 - sqrt(1-4*x))/(2*x) is the Catalan function (A000108). (End)
MATHEMATICA
Table[(-1)^n * HypergeometricU[-n, n, -1], {n, 0, 20}] (* Vaclav Kotesovec, Oct 02 2017 *)
PROG
(PARI) a(n)=sum(k=0, n, k!*binomial(2*n-1, k)*binomial(n, k))
(PARI) a(n) = n!*pollaguerre(n, n-1, -1); \\ Seiichi Manyama, May 01 2021
(Magma) [Factorial(n)*Evaluate(LaguerrePolynomial(n, n-1), -1): n in [0..40]]; // G. C. Greubel, Aug 11 2022
(SageMath) [factorial(n)*gen_laguerre(n, n-1, -1) for n in (0..40)] # G. C. Greubel, Aug 11 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 22 2008
STATUS
approved