OFFSET
1,1
MATHEMATICA
onpQ[n_]:=Module[{pr=PrimitiveRoot[n, 1]}, OddQ[pr]&&!PrimeQ[pr]]; Select[ Prime[ Range[7100]], onpQ] (* Harvey P. Dale, Jun 22 2020 *)
PROG
(PARI) forprime(n=1, 100000, if((-1)^(lift(znprimroot(n))*(1-isprime(lift(znprimroot(n)))))==-1, print1(n, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, May 09 2002
STATUS
approved