OFFSET
1,4
LINKS
FORMULA
MATHEMATICA
f[p_, e_] := Module[{q = NextPrime[p]}, If[e == 1, q - p - 1, q^e - (p + 1)*q^(e - 1) + p*q^(e - 2)]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Dec 11 2024 *)
PROG
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Dec 11 2024
STATUS
approved