OFFSET
1,2
FORMULA
E.g.f.: Sum_{k>=1} (exp(2*x^k) - 1) / (exp(2*x^k) + 1).
a(n) = n! * Sum_{d|n} A155585(d) / d!.
a(n) = n! * Sum_{d|n, d odd} 2^(d + 1) * (2^(d + 1) - 1) * Bernoulli(d + 1) / (d + 1)!.
MATHEMATICA
nmax = 21; CoefficientList[Series[Sum[Tanh[x^k], {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]! // Rest
A155585[n_] := Sum[StirlingS2[n, k] (-2)^(n - k) k!, {k, 0, n}]; a[n_] := n! DivisorSum[n, A155585[#]/#! &]; Table[a[n], {n, 1, 21}]
Table[n! DivisorSum[n, 2^(# + 1) (2^(# + 1) - 1) BernoulliB[# + 1]/(# + 1)! &, OddQ[#] &], {n, 1, 21}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 16 2019
STATUS
approved