close
login
A390015
E.g.f. A(x) satisfies A(x) = exp( x * (1-x^3)^2 * A(x) ).
4
1, 1, 3, 16, 77, 576, 5287, 57184, 791289, 12866464, 244736171, 5315199264, 128583465397, 3426391314256, 99322564531023, 3107106602767696, 104266181082952433, 3733781580291963456, 142103965349570443219, 5727647815719442313536, 243730370521404874282221
OFFSET
0,3
LINKS
FORMULA
a(n) = n! * Sum_{k=0..floor(n/3)} (-1)^k * (n-3*k+1)^(n-3*k-1) * binomial(2*(n-3*k),k)/(n-3*k)!.
E.g.f.: exp( -LambertW(-x * (1-x^3)^2) ).
MATHEMATICA
a[n_]:=n!*Sum[(-1)^k*(n-3*k+1)^(n-3*k-1)*Binomial[2*(n-3*k), k]/(n-3*k)!, {k, 0, Floor[n/3]}]; Table[a[n], {n, 0, 25}] (* Vincenzo Librandi, Oct 25 2025 *)
PROG
(PARI) a(n) = n!*sum(k=0, n\3, (-1)^k*(n-3*k+1)^(n-3*k-1)*binomial(2*(n-3*k), k)/(n-3*k)!);
(Magma) a := func< n | Factorial(n) * &+[ (-1)^k * (n-3*k + 1)^(n-3*k - 1) * Binomial(2*(n-3*k), k) / Factorial(n-3*k) : k in [0..Floor(n/3)]] >;
[a(n) : n in [0..25]]; // Vincenzo Librandi, Oct 25 2025
CROSSREFS
Cf. A389987.
Sequence in context: A005386 A053572 A329806 * A390216 A371350 A309915
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Oct 22 2025
STATUS
approved