OFFSET
0,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..300
FORMULA
E.g.f. A(x) satisfies A(x) = 1/(exp(-x*A(x)) - (x*A(x))^3).
a(n) = n! * Sum_{k=0..floor(n/3)} (n+k+1)^(n-3*k-1) * binomial(n+k+1,k)/(n-3*k)!.
MATHEMATICA
Table[n!*Sum[(n+k+1)^(n-3*k-1)*Binomial[n+k+1, k]/(n-3*k)!, {k, 0, Floor[n/3]}], {n, 0, 18}] (* Vincenzo Librandi, Jan 30 2026 *)
PROG
(PARI) a(n) = n!*sum(k=0, n\3, (n+k+1)^(n-3*k-1)*binomial(n+k+1, k)/(n-3*k)!);
(Magma) [Factorial(n) * &+[(n+k+1)^(n-3*k-1)* Binomial(n+k+1, k) / Factorial(n-3*k) : k in [0..Floor(n/3)]] : n in [0..25] ]; // Vincenzo Librandi, Jan 30 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 30 2026
STATUS
approved
