OFFSET
1,1
LINKS
Richard Fischer, Fermatquotienten von 2 bis 1052, Dec 19 2019.
Wikipedia, Wieferich prime
MATHEMATICA
Select[Range[3*10^6], PrimeQ[#] && PowerMod[31, # - 1, #^2] == 1 &] (* Amiram Eldar, May 05 2021 *)
PROG
(PARI) forprime(p=2, 1e8, if(Mod(31, p^2)^(p-1)==1, print1(p", ")))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Seiichi Manyama, Jan 16 2020
STATUS
approved