OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..775
MATHEMATICA
q[k_] := DivisorSigma[-1, k^2] > 2 && AllTrue[Divisors[k], DivisorSigma[-1, #*k] <= 2 || # == k &]; Select[Range[10^6], q]
PROG
(PARI) isok(k) = fordiv(k, d, if(d < k && sigma(d*k, -1) > 2, return(0))); sigma(k^2, -1) > 2;
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 06 2025
STATUS
approved
