OFFSET
1,3
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..16384
Michael Baake and Robert V. Moody, Similarity submodules and root systems in four dimensions, Canad. J. Math., Vol. 51, No. 6 (1999), 1258-1276.
Michael Baake and Peter Zeiner, "Similar Sublattices", Ch. 3.5 in Aperiodic Order, Vol. 2: Crystallography and Almost Periodicity, Cambridge, 2017, see page 105.
FORMULA
Multiplicative with a(2^p) = 1, a(p^e) = (e+1)*p^e + (2*(1+(e*p-e-1)*p^e))/((p-1)^2), p>2. - Christian G. Bower, May 21 2005
From Amiram Eldar, May 26 2025: (Start)
Dirichlet g.f.: (zeta(s-1)^2 * zeta(s)^2 / zeta(2*s)) * (1 - 1/2^(s-1))^2/(1 + 1/2^s).
Sum_{k=1..n} a(k) ~ (n^2/4)*(log(n) + 2*gamma - 1/2 + 11*log(2)/5 + 2*zeta'(2)/zeta(2) - 2*zeta'(4)/zeta(4)), where gamma is Euler's constant (A001620). (End)
MATHEMATICA
Array[Apply[Times, FactorInteger[#] /. {p_, e_} /; p > 0 :> If[1 <= p <= 2, 1, (e + 1) p^e + (2 (1 + (e p - e - 1)*p^e))/((p - 1)^2)]] &, 64] (* Michael De Vlieger, Mar 02 2018 *)
PROG
(PARI) fp(p, e) = if (p % 2, (e+1)*p^e + 2*(1-(e+1)*p^e+e*p^(e+1))/(p-1)^2, 1);
a(n) = { my(f = factor(n)); prod(i=1, #f~, fp(f[i, 1], f[i, 2])); } \\ Michel Marcus, Mar 03 2014
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Michael Baake (baake(AT)miles.math.ualberta.ca)
EXTENSIONS
More terms from Michel Marcus, Mar 03 2014
STATUS
approved