A365013 E.g.f. satisfies A(x) = exp( x*A(x)/(1 - x * A(x)^3) ).

1, 1, 5, 58, 1061, 26536, 843457, 32553424, 1478813513, 77304347776, 4571222616701, 301696674682624, 21985118975444077, 1753288356936334336, 151887264799071753785, 14203597499192539334656, 1426051485043745729079953, 153000280727938469281693696

Offset

0,3

Links

Michael De Vlieger, Table of n, a(n) for n = 0..341

Formula

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

Mathematica

Array[#!*Sum[ (3 # - 2 k + 1)^(k - 1)*Binomial[# - 1, # - k]/k!, {k, 0, #}] &, 18, 0] (* Michael De Vlieger, Aug 18 2023 *)

Prog

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

Crossrefs

Cf. A052873, A365012.

Sequence in context: A256218 A151424 A097631 * A130768 A195947 A371329

Adjacent sequences: A365010 A365011 A365012 * A365014 A365015 A365016

Keyword

nonn

Author

Seiichi Manyama, Aug 15 2023

Status

approved