OFFSET
1,4
FORMULA
MATHEMATICA
nmax = 39; CoefficientList[Series[Sum[x^(2 k)/(1 - x^k - x^(2 k) - x^(3 k)), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
PROG
(PARI) f(n) = ([0, 1, 0; 0, 0, 1; 1, 1, 1]^n)[1, 3]; \\ A000073
a(n) = sumdiv(n, d, f(d)); \\ Michel Marcus, Sep 20 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Sep 19 2022
STATUS
approved