[go: up one dir, main page]

login
A036074
Expansion of e.g.f. exp((exp(p*x) - p - 1)/p + exp(x)) for p=4.
8
1, 2, 9, 55, 412, 3619, 36333, 408888, 5080907, 68914023, 1011165446, 15935379409, 268125052373, 4792458452162, 90605469012877, 1805135197261131, 37775862401203916, 827992670793489263
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.
LINKS
FORMULA
a(n) = sum(sum(binomial(m,i)*sum(binomial(i,j)*(1/4)^j*(3*j+i)^n,j,0,i)*(-5/4)^(m-i),i,0,m)/m!,m,1,n), n > 0. - Vladimir Kruchinin, Sep 14 2010
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=4. - Vaclav Kotesovec, Jul 03 2022
a(n) ~ (4*n/LambertW(4*n))^n * exp(n/LambertW(4*n) + (4*n/LambertW(4*n))^(1/4) - n - 5/4) / sqrt(1 + LambertW(4*n)). - Vaclav Kotesovec, Jul 10 2022
MATHEMATICA
mx = 16; p = 4; 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] * 4^k * BellB[k, 1/4] * BellB[n-k], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Jun 29 2022 *)
PROG
(Maxima) a(n):=sum(sum(binomial(m, i)*sum(binomial(i, j)*(1/4)^j*(3*j+i)^n, j, 0, i)*(-5/4)^(m-i), i, 0, m)/m!, m, 1, n); /* Vladimir Kruchinin, Sep 14 2010 */
CROSSREFS
KEYWORD
nonn
EXTENSIONS
Edited by N. J. A. Sloane, Jul 11 2008 at the suggestion of Franklin T. Adams-Watters
STATUS
approved