OFFSET
1,2
FORMULA
E.g.f.: Sum_{i>=1} Sum_{j>=1} (-1)^(j + 1) * x^(i*(2*j - 1)) / (2*j - 1).
a(n) = (n - 1)! * Sum_{d|n, n/d odd} (-1)^((n/d - 1)/2) * d.
MATHEMATICA
nmax = 21; CoefficientList[Series[Sum[ArcTan[x^k], {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]! // Rest
Table[(n - 1)! DivisorSum[n, (-1)^((n/# - 1)/2) # &, OddQ[n/#] &], {n, 1, 21}]
PROG
(PARI) a(n) = (n-1)!*sumdiv(n, d, if (n/d % 2, (-1)^((n/d - 1)/2)*d)); \\ Michel Marcus, Dec 17 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 16 2019
STATUS
approved