OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..384
FORMULA
E.g.f.: log(1 + Sum_{k>=1} k^4*x^k/k!).
a(0) = 0; a(n) = n^4 - (1/n)*Sum_{k=1..n-1} binomial(n,k)*(n - k)^4*k*a(k).
MAPLE
a:=series(log(1 + exp(x)*x*(1 + 7*x + 6*x^2 + x^3)), x=0, 21): seq(n!*coeff(a, x, n), n=0..20); # Paolo P. Lava, Mar 26 2019
MATHEMATICA
nmax = 20; CoefficientList[Series[Log[1 + Exp[x] x (1 + 7 x + 6 x^2 + x^3)], {x, 0, nmax}], x] Range[0, nmax]!
a[n_] := a[n] = n^4 - Sum[Binomial[n, k] (n - k)^4 k a[k], {k, 1, n - 1}]/n; a[0] = 0; Table[a[n], {n, 0, 20}]
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Feb 07 2019
STATUS
approved