OFFSET
1,1
PROG
(PARI) isp(n) = for(b=2, oo, if(Mod(b, n^2)^(n-1)==1, return(ispseudoprime(b))))
forprime(p=1, 200, if(!isp(p), print1(p, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Jul 12 2022
STATUS
approved