OFFSET
1,1
COMMENTS
If k is a term then all the positive multiples of k are terms of A359565.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
q[n_] := Max[Tally[EulerPhi[Divisors[n]]][[;; , 2]]] > 2; primQ[n_] := q[n] && AllTrue[Divisors[n], # == n || !q[#] &]; Select[Range[17000], primQ]
PROG
(PARI) is1(k) = vecmax(matreduce(apply(x->eulerphi(x), divisors(k)))[, 2]) > 2;
is(k) = fordiv(k, d, if(is1(d), return(d==k))); return(0);
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jan 06 2023
STATUS
approved