OFFSET
1,2
COMMENTS
MATHEMATICA
f[p_, e_] := (e*(p - 1)/p + 1)*p^e; pil[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[500], Divide @@ DivisorSigma[{1, 0}, pil[#]] == # &] (* Amiram Eldar, May 17 2025 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ctibor O. Zizka, May 17 2025
STATUS
approved
