OFFSET
1,2
COMMENTS
Multiplicative with a(2^p) = 3, a(p^e) = (e+1)*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
Research Group Michael Baake, Preprints & Recent Articles: Algebra, Combinatorics and Number Theory
Michael Baake and Robert V. Moody, Similarity submodules and root systems in four dimensions, arXiv:math/9904028 [math.MG].
Michael Baake and Robert V. Moody, Similarity submodules and root systems in four dimensions, Canad. J. Math. (1999), 51 1258-1276.
FORMULA
Baake and Moody give Dirichlet generating function.
MATHEMATICA
Clear[ a, f ]; a[ {p_, r_} ] := If[ p == 2, 3, (r + 1)*p^r + (2*(1 - (r + 1)*p^r + r*p^(r + 1)))/(p - 1)^2 ]; f[ m_Integer ] := f[ m ] = Times @@ a /@ FactorInteger[ m ]; (* f[ m ] is number of similar sublattices of Z^4 of index m^2 *)
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)); a045771 = prod(i=1, #f~, fp(f[i, 1], f[i, 2])); if (n % 2, a045771, 3*a045771); } \\ Michel Marcus, Mar 03 2014
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Michael Baake (michael.baake(AT)uni-tuebingen.de)
STATUS
approved