OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..170
FORMULA
a(n) ~ n! * (1 + 1/n + 3/n^2 + 12/n^3 + 66/n^4 + 450/n^5 + 3679/n^6 + 35260/n^7 + 388511/n^8 + 4844584/n^9 + 67502450/n^10), for next coefficients see A248871. - Vaclav Kotesovec, Mar 14 2015
G.f.: Product_{n>=1} 1/(1-x^n)^(n!). - Vaclav Kotesovec, Aug 04 2015
MAPLE
EulerTrans := proc(p) local b; b := proc(n) option remember; local d, j;
`if`(n=0, 1, add(add(d*p(d), d=numtheory[divisors](j)) *b(n-j), j=1..n)/n) end end:
MATHEMATICA
EulerTrans[p_] := Module[{b}, b[n_] := b[n] = Module[{d, j}, If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[j]}]*b[n-j], {j, 1, n}]/n]]; b]; A107895 = EulerTrans[Factorial]; Table[A107895[n], {n, 0, 22}] (* Jean-François Alcover, Feb 25 2014, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Thomas Wieder, May 26 2005
STATUS
approved