OFFSET
1,2
LINKS
MATHEMATICA
f[p_, e_] := (e + 1)*(pn = NextPrime[p])^e*(pn - 1)/(pn^(e + 1) - 1); a[1] = 1; a[n_] := Numerator[Times @@ f @@@ FactorInteger[n]]; Array[a, 100] (* Amiram Eldar, Oct 31 2021 *)
PROG
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Antti Karttunen, Oct 31 2021
STATUS
approved