OFFSET
1,1
MATHEMATICA
f[p_, e_] := If[EvenQ[e], (p^(e + 1) - 1)/(p - 1), 1]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[300], PrimeQ[s[#]] &] (* Amiram Eldar, Feb 25 2022 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 24 2022
STATUS
approved