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”).

A036078
E.g.f.: exp((exp(p*x)-p-1)/p+exp(x)) for p=8.
0
1, 2, 13, 127, 1508, 20859, 332557, 6019108, 121462267, 2692076295, 64846340130, 1684713690917, 46916754353013, 1393010598959594, 43889040801834505, 1461369418905803027, 51243270154712083052
OFFSET
0,2
REFERENCES
T. S. Motzkin, Sorting numbers for cylinders and other classification numbers, in Combinatorics, Proc. Symp. Pure Math. 19, AMS, 1971, pp. 167-176.
T. S. Motzkin, Sorting numbers ...: for a link to an annotated scanned version of this paper see A000262.
FORMULA
a(n) ~ exp(exp(p*r)/p + exp(r) - 1 - 1/p - n) * (n/r)^(n + 1/2) / sqrt((1 + p*r)*exp(p*r) + (1 + r)*exp(r)), where r = LambertW(p*n)/p - 1/(1 + p/LambertW(p*n) + n^(1 - 1/p) * (1 + LambertW(p*n)) * (p/LambertW(p*n))^(2 - 1/p)) for p=8. - Vaclav Kotesovec, Jul 03 2022
a(n) ~ (8*n/LambertW(8*n))^n * exp(n/LambertW(8*n) + (8*n/LambertW(8*n))^(1/8) - n - 9/8) / sqrt(1 + LambertW(8*n)). - Vaclav Kotesovec, Jul 10 2022
MATHEMATICA
mx = 16; p = 8; Range[0, mx]! CoefficientList[ Series[ Exp[ (Exp[p*x] - p - 1)/p + Exp[x]], {x, 0, mx}], x] (* Robert G. Wilson v, Dec 12 2012 *)
Table[Sum[Binomial[n, k] * 8^k * BellB[k, 1/8] * BellB[n-k], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Jun 29 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Edited by N. J. A. Sloane, Jul 11 2008 at the suggestion of Franklin T. Adams-Watters.
STATUS
approved