OFFSET
1,2
COMMENTS
LINKS
EXAMPLE
MATHEMATICA
q[k_] := Module[{e = FactorInteger[k][[;; , 2]]}, AllTrue[e, # > 2 &] && Divisible[Times @@ e, Times @@ DivisorSigma[0, e]]]; Select[Range[140000], # == 1 || q[#] &]
PROG
(PARI) isok(k) = if(k == 1, 1, my(e = factor(k)[, 2]); vecmin(e) > 2 && !(vecprod(e) % vecprod(apply(x -> numdiv(x), e))));
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 14 2025
STATUS
approved
