OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
q[k_] := DivisorSigma[-1, k^2] > 2 && AllTrue[Divisors[k], DivisorSigma[-1, #^2] <= 2 || # == k &]; Select[Range[3200], q]
PROG
(PARI) is1(k) = {my(f = factor(k)); prod(i = 1, #f~, f[i, 2] *= 2); sigma(f, -1) > 2; }
isok(k) = if(!is1(k), 0, fordiv(k, d, if(d < k && is1(d), return(0))); 1);
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Mar 05 2025
STATUS
approved
