OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Mario Hampejs and László Tóth, On the subgroups of finite Abelian groups of rank three, Annales Univ. Sci. Budapest., Sect. Comp. 39 (2013) 111-124.
FORMULA
For a prime p: a(p) = 2*(p^2+p+2). - Vladeta Jovovic, Oct 22 2001
Multiplicative with a(p^e) = Sum_{j=0..2*e} (e - floor((j - 1)/2))*(2*j - floor((j - 1)/2))*p^(2*e - j). - Amiram Eldar, Nov 29 2022
MAPLE
A064803 := proc(n)
local a, f, nu, p, j ;
a := 1 ;
for f in ifactors(n)[2] do
nu := op(2, f) ;
p := op(1, f) ;
add( (nu-floor((j-1)/2))*(2*j-floor((j-1)/2))*p^(2*nu-j), j=0..2*nu) ;
a := a*% ;
end do:
a ;
end proc: # R. J. Mathar, May 11 2013
MATHEMATICA
f[p_, e_] := Sum[(e - Floor[(j - 1)/2])*(2*j - Floor[(j - 1)/2])*p^(2*e - j), {j, 0, 2*e}]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 29 2022 *)
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Oct 21 2001
EXTENSIONS
More terms from Laszlo Toth, May 11 2013
STATUS
approved