A092342 a(n) = sigma_3(3n+1).

1, 73, 344, 1134, 2198, 4681, 6860, 11988, 15751, 25112, 29792, 44226, 50654, 73710, 79508, 109512, 117993, 160454, 167832, 219510, 226982, 299593, 300764, 390096, 389018, 500780, 493040, 620298, 619164, 779220, 756112, 934416, 912674, 1149823, 1092728

Offset

0,2

Links

Table of n, a(n) for n=0..34.

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

a(n) = A000731(n) + 81*A033690(n-1). - 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

Trisection of A001158.

Cf. A000731, A013662, A033690, A045823, A091986, A092341, A092343.

Sequence in context: A005108 A142810 A142279 * A170916 A200908 A031418

Adjacent sequences: A092339 A092340 A092341 * A092343 A092344 A092345

Keyword

nonn,easy

Author

N. J. A. Sloane, Mar 20 2004

Status

approved