OFFSET
0,2
LINKS
Paolo Xausa, Table of n, a(n) for n = 0..445
FORMULA
a(n) = (n + 3)! * ((n + 5) * H(n + 3) - 3*n - 9), where H(n) is the n-th harmonic number.
a(n) = A396981(n, 12).
MAPLE
prec := 20: egf := (1 - 23*z - (30 - 6*z)*log(1-z)) /(1 - z)^5:
ser := series(egf, z, prec+1): seq(n!*coeff(ser, z, n), n = 0..19);
# Alternative:
a := n -> (n + 3)! * ((n + 5) * harmonic(n + 3) - 3*n - 9):
MATHEMATICA
A397013[n_] := (n + 3)!*((n + 5)*HarmonicNumber[n + 3] - 3*n - 9);
Array[A397013, 20, 0] (* Paolo Xausa, Jun 14 2026 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Jun 14 2026
STATUS
approved
