OFFSET
1,1
MATHEMATICA
Select[Range[2, 4100000], Function[b, AllTrue[{2, 3, 5, 7, 11, 13}, PowerMod[b, (# - 1), #^2] == 1 &]]] (* Michael De Vlieger, Dec 10 2020 *)
PROG
(PARI) is(n) = forprime(p=1, 13, if(Mod(n, p^2)^(p-1)!=1, return(0))); 1
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Dec 08 2020
STATUS
approved