|
|
A074932
|
|
Row sums of unsigned triangle A075513.
|
|
3
|
|
|
1, 3, 18, 170, 2200, 36232, 725200, 17095248, 463936896, 14246942336, 488428297984, 18491942300416, 766293946203136, 34498781924766720, 1676731077272217600, 87501958444207351808, 4880017252828686155776
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = sum(|A075513(n, m)|, m=0..n-1) = sum(binomial(n-1, m)*(m+1)^(n-1), m=0..n-1), n>=1.
E.g.f.: log(G(x)) where G(x) = (1 + exp(2*x*G(x)))/2 is the e.g.f. of A201595. - Paul D. Hanna, Aug 02 2012
a(n) ~ n!/(sqrt(2*Pi*(1+LambertW(exp(-1))))*n^(3/2)*LambertW(exp(-1))^n). - Vaclav Kotesovec, Dec 04 2012
|
|
EXAMPLE
|
E.g.f.: A(x) = x + 3*x^2/2! + 18*x^3/3! + 170*x^4/4! + 2200*x^5/5! +...
where exp(A(x)) = 1 + x + 4*x^2/2! + 28*x^3/3! + 288*x^4/4! + 3936*x^5/5! + 67328*x^6/6! +...+ A201595(n)*x^n/n! +...
|
|
MATHEMATICA
|
Rest[CoefficientList[Series[Log[x-LambertW[-x*Exp[x]]]-Log[2*x], {x, 0, 20}], x]* Range[0, 20]!] (* Vaclav Kotesovec, Dec 04 2012 *)
|
|
PROG
|
(PARI) {a(n)=sum(k=0, n-1, binomial(n-1, k)*(k+1)^(n-1))} \\ Paul D. Hanna, Aug 02 2012
(PARI) {a(n) = my(A); if( n<0, 0, A = O(x); for(k=1, n, A = log( (1 + exp( 2*x * exp(A))) / 2 )); n! * polcoeff(A, n))}; /* Michael Somos, Apr 10 2018 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|