OFFSET
0,6
FORMULA
a(0) = 0; a(n) = binomial(n,4) + (1/n) * Sum_{k=1..n-1} binomial(n,k) * binomial(n-k,4) * k * a(k).
a(n) ~ (n-1)! / (4*LambertW(3^(1/4)/2^(5/4)))^n. - Vaclav Kotesovec, Aug 08 2021
a(n) = n! * Sum_{k=1..floor(n/4)} k^(n-4*k-1)/(24^k * (n-4*k)!). - Seiichi Manyama, Dec 14 2023
MATHEMATICA
nmax = 24; CoefficientList[Series[-Log[1 - x^4 Exp[x]/4!], {x, 0, nmax}], x] Range[0, nmax]!
a[0] = 0; a[n_] := a[n] = Binomial[n, 4] + (1/n) Sum[Binomial[n, k] Binomial[n - k, 4] k a[k], {k, 1, n - 1}]; Table[a[n], {n, 0, 24}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 01 2021
STATUS
approved