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/(1 + (1 - exp((x*A(x))^4))/(x*A(x))^3).
a(n) = (1/(n+1)) * Sum_{k=0..floor(n/4)} (2*n-4*k)! * Stirling2(n-3*k,n-4*k)/(n-3*k)!.
MATHEMATICA
Table[(1/(n+1))*Sum[(2*n-4*k)!*Abs[StirlingS2[n-3*k, n-4*k]/(n-3*k)!], {k, 0, Floor[n/4]}], {n, 0, 21}] (* Vincenzo Librandi, Jan 24 2026 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(serreverse(x+(1-exp(x^4))/x^2)/x))
(Magma) [ 1/(n+1)* &+[ Factorial(2*n-4*k)*StirlingSecond(n-3*k, n-4*k)/Factorial(n-3*k) : k in [0..Floor(n/4)] ] : n in [0..18] ]; // Vincenzo Librandi, Jan 24 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 22 2026
STATUS
approved
