OFFSET
0,1
COMMENTS
LINKS
PROG
(PARI) for(n=0, 20, b=2; goodwief=0; while(goodwief==0, badwief=0; forprime(p=1, 10^n, if(Mod(b, p^2)^(p-1)==1, badwief++; break({1}))); if(badwief==0, forprime(p=10^n, 10^(n+1), if(Mod(b, p^2)^(p-1)==1, print1(b, ", "); goodwief++; break({1})))); b++))
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Felix Fröhlich, Apr 02 2015
STATUS
approved