A364980 E.g.f. satisfies A(x) = 1 + x*A(x)*exp(x*A(x)^2).

1, 1, 4, 33, 412, 6945, 147846, 3807601, 115151464, 4001162913, 157096369450, 6878742553881, 332361857826780, 17566215943990753, 1008161606338206334, 62440146891413434305, 4151012174991960338896, 294834882756167048975553

Offset

0,3

Links

Table of n, a(n) for n=0..17.

Formula

a(n) = n! * Sum_{k=0..n} k^(n-k) * binomial(2*n-k+1,k)/( (2*n-k+1)*(n-k)! ).

a(n) ~ sqrt((1 + r*s^2)/(6 + 4*r*s^2)) * n^(n-1) / (exp(n) * r^(n + 1/2)), where r = 0.2190923703746024362724546703711998154573791458000... and s = 1.747404632046819382844696016554403302840973484745... are real roots of the system of equations 1 + exp(r*s^2)*r*s = s, 2*r*s^2*(s-1) = 1. - Vaclav Kotesovec, Nov 18 2023

Mathematica

Join[{1}, Table[n! * Sum[k^(n-k) * Binomial[2*n-k+1, k] / ((2*n-k+1)*(n-k)!), {k, 0, n}], {n, 1, 20}]] (* Vaclav Kotesovec, Nov 18 2023 *)

Prog

(PARI) a(n) = n!*sum(k=0, n, k^(n-k)*binomial(2*n-k+1, k)/((2*n-k+1)*(n-k)!));

Crossrefs

Cf. A006153, A161633, A161635, A364981.

Sequence in context: A213641 A343673 A343686 * A360234 A111534 A162655

Adjacent sequences: A364977 A364978 A364979 * A364981 A364982 A364983

Keyword

nonn

Author

Seiichi Manyama, Aug 15 2023

Status

approved