OFFSET
1,2
FORMULA
MATHEMATICA
nmax = 48; CoefficientList[Series[Sum[x^(2^k) /(1 - x^(2^k))^4, {k, 0, Floor[Log[2, nmax]] + 1}], {x, 0, nmax}], x] // Rest
a[n_] := If[EvenQ[n], a[n/2] + n (n + 1) (n + 2)/6, n (n + 1) (n + 2)/6]; Table[a[n], {n, 1, 48}]
Table[(1/6) DivisorSum[n, Boole[IntegerQ[Log[2, n/#]]] # (# + 1) (# + 2) &], {n, 1, 48}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Oct 08 2020
STATUS
approved