OFFSET
1,1
MATHEMATICA
q[n_] := Module[{s = DivisorSigma[1, n], d = DivisorSigma[0, n]}, Divisible[s*d, s - d]]; Select[Range[2, 100], q] (* Amiram Eldar, Apr 07 2024 *)
PROG
(PARI) is(n) = {my(f = factor(n), s = sigma(f), d = numdiv(f)); n > 1 && !((s*d) % (s-d)); } \\ Amiram Eldar, Apr 07 2024
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Ctibor O. Zizka, Nov 01 2008
STATUS
approved