OFFSET
0,2
FORMULA
Expansion of q^(-1/3) * c(q) * (c(q)^3 + b(q)^3 / 3) in powers of q where b(), c() are cubic AGM functions. - Michael Somos, Aug 22 2007
If b(3*n) = 0, b(3*n+1) = a(n), b(3*n+2) = A092343(n), then b(n) is multiplicative with b(3^e) = 0^e, b(p^e) = (p^(3*e+3) - 1) / (p^3 - 1) otherwise. - Michael Somos, Aug 22 2007
Sum_{k=0..n} a(k) ~ (20*zeta(4)/3) * n^4. - Amiram Eldar, Dec 12 2023
EXAMPLE
q + 73*q^4 + 344*q^7 + 1134*q^10 + 2198*q^13 + 4681*q^16 + ...
MATHEMATICA
DivisorSigma[3, 3*Range[0, 40]+1] (* Harvey P. Dale, Apr 22 2019 *)
PROG
(PARI) {a(n) = if(n<0, 0, sigma(3*n+1, 3))} /* Michael Somos, Aug 22 2007 */
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Mar 20 2004
STATUS
approved