OFFSET
1,3
COMMENTS
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
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..16384
M. Baake and R. V. Moody, Similarity submodules and root systems in four dimensions, Canad. J. Math. (1999), 51 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.
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