OFFSET
1,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
M. A. Basoco, On the Fourier developments of a certain class of theta quotients, Bull. Amer. Math. Soc. 49 (1943), 299-306.
FORMULA
Multiplicative with a(2^e) = 4^e, a(p^e) = p^e * (p^(e+1) - 1) / (p - 1) if p>2.
G.f.: Sum_{k>0} k^2 * (x^k + x^(3*k)) / (1 - x^(2*k))^2. [see Basoco (1943) bottom page 305]
G.f.: Sum_{k>0} k^2 * (3 - (-1)^k)/4 * x^k / (1 - x^k)^2.
G.f.: Sum_{k>0 odd} k * (x^k + x^(2*k)) / (1 - x^k)^3.
a(n) = n * A002131(n). a(2*n) = 4 * a(n).
a(n) = A007331(n) - 4 * Sum_{k>0} A002131(k) * A002131(n-k). [see Basoco (1943) page 305 equation (9)]
Sum_{k=1..n} a(k) ~ c * n^3, where c = Pi^2/24 = 0.411233... (A222171). - Amiram Eldar, Nov 30 2022
EXAMPLE
G.f. = x + 4*x^2 + 12*x^3 + 16*x^4 + 30*x^5 + 48*x^6 + 56*x^7 + 64*x^8 + ...
a(6) = 48 = 6 * (2 + 6). a(9) = 117 = 9 * (1 + 3 + 9). a(10) = 120 = 10 * (2 + 10).
MATHEMATICA
PROG
(PARI) {a(n) = if( n<1, 0, n * sumdiv(n, d, d * (n/d % 2)))};
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Michael Somos, Aug 19 2013
STATUS
approved