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