OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
J. W. L. Glaisher, On the representations of a number as the sum of two, four, six, eight, ten, and twelve squares, Quart. J. Math. 38 (1907), 1-62 (see p. 4 and p. 8).
FORMULA
a(n) = sigma(n)*delta(n).
Multiplicative with a(2^e) = 2^(e+1) - 1 and a(p^e) = (e+1)*(p^(e+1)-1)/(p-1) for p > 2. - Amiram Eldar, Nov 01 2022
From Amiram Eldar, Dec 04 2023: (Start)
Dirichlet g.f.: (4^s - 3*2^s + 2)/(4^s - 2) * (zeta(s)*zeta(s-1))^2/zeta(2*s-1).
Sum_{k=1..n} a(k) ~ (Pi^4/(168*zeta(3))) * n^2 * (log(n) + 2*gamma - 1/2 + 22*log(2)/21 + 2*zeta'(2)/zeta(2) - 2*zeta'(3)/zeta(3)), where gamma is Euler's constant (A001620). (End)
MATHEMATICA
Array[DivisorSum[#, 1 &, OddQ] DivisorSigma[1, #] &, 69] (* Michael De Vlieger, Nov 22 2019 *)
f[p_, e_] := (e+1)*(p^(e+1)-1)/(p-1); f[2, e_] := 2^(e+1) - 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 01 2022 *)
CROSSREFS
KEYWORD
nonn,mult,easy
AUTHOR
Omar E. Pol, Jul 14 2019
STATUS
approved