OFFSET
1,4
COMMENTS
a(n) = A000688(n) for n < 32.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(p^e) = A001400(e).
Sum_{k=1..n} a(k) ~ Pi^6*zeta(3)*n/540 + Pi^2*zeta(1/2)*zeta(3/2)*sqrt(n)/6 + zeta(1/3)*zeta(2/3)*zeta(4/3)*n^(1/3) + zeta(1/4)*zeta(1/2)*zeta(3/4)*n^(1/4).
PROG
(PARI) for(n=1, 200, print1(direuler(p=2, n, 1/(1 - X)/(1 - X^2)/(1 - X^3)/(1 - X^4))[n], ", "))
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Vaclav Kotesovec, Mar 14 2023
STATUS
approved