OFFSET
1,1
MATHEMATICA
f[p_, e_] := (p^(e+1)-2*p+1) * (p^(e+1)-1)/((p-1)^2 * (p^e - 1)); s[n_] := (Times @@ f @@@ FactorInteger[n])/n; Select[s /@ Range[2, 10^5], IntegerQ] (* Amiram Eldar, Sep 19 2022 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Yasutoshi Kohmoto Sep 30 2006
EXTENSIONS
More terms from R. J. Mathar, Oct 01 2006
a(12)-a(14) from Amiram Eldar, Sep 19 2022
STATUS
approved