OFFSET
1
PROG
(PARI) a(n) = my(p=prime(n)); for(b=2, oo, if(Mod(b, p^2)^(p-1)==1, return(ispseudoprime(b))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Jul 12 2022
STATUS
approved