OFFSET
1,2
LINKS
T. D. Noe, Table of n, a(n) for n=1..10000
FORMULA
Dirichlet g.f.: zeta^2(x)zeta^2(x-1).
Multiplicative with a(2^e) = (e-1) 2^(e+2) + e + 5, a(p^e) = ((1+e)p^(e+3) - (3+e)(p^(e+2)-p+1) + 2)/(p-1)^3, p > 2. - Mitch Harris, Jun 27 2005 [corrected by Amiram Eldar, Oct 16 2022 and Sep 12 2023]
Sum_{k=1..n} a(k) ~ Pi^2 * n^2 * (2*Pi^2 * log(n) + (4*gamma - 1)*Pi^2 + 24*zeta'(2)) / 144, where gamma is the Euler-Mascheroni constant A001620 and Zeta'(2) = A073002. Equivalently, Sum_{k=1..n} a(k) ~ Pi^4 * n^2 * (2*log(n) - 1 + 8*gamma - 48*log(A) + 4*log(2*Pi)) / 144, where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Jan 28 2019
MATHEMATICA
f[p_, e_] := ((e + 1)*p^(e + 3) - (e + 3)*(p^(e + 2) - p + 1) + 2)/(p - 1)^3; f[2, e_] := (e - 1)*2^(e + 2) + e + 5; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Oct 16 2022 *)
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
STATUS
approved