OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
q[n_] := !UnsameQ @@ EulerPhi[Divisors[n]]; primQ[n_] := q[n] && AllTrue[Divisors[n], # == n || !q[#] &]; Select[Range[1, 30000, 2], primQ]
PROG
(PARI) is1(k) = k>1 && k%2 && numdiv(k) > #Set(apply(x->eulerphi(x), divisors(k)));
is(k) = fordiv(k, d, if(is1(d), return(d==k))); return(0);
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jan 06 2023
STATUS
approved