|
|
A266829
|
|
Primes p such that a prime q < p exists with p^(q-1) == 1 (mod q^2) and q^(p-1) == 1 (mod p^2), i.e., primes that are the larger member of a double Wieferich prime pair.
|
|
5
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
There are no further terms less than 10^6 (cf. Ernvall, Metsänkylä, 1997, p. 1360).
There are no further terms p less than 2^(1/3)*10^10 with p*q <= 10^15 and p and q both odd. (cf. Logan, Mossinghoff, results 4.2.). - Felix Fröhlich, May 29 2016 [Corrected. Felix Fröhlich, Aug 05 2018]
|
|
LINKS
|
|
|
MATHEMATICA
|
fQ[p_] := Block[{q = 2}, While[q < p && (PowerMod[p, q - 1, q^2] != 1 || PowerMod[q, p - 1, p^2] != 1), q = NextPrime@ q]; If[q < p, True, False]]; p = 3; lst = {}; While[p < 1000000, If[fQ@ p, AppendTo[lst, p]]; p = NextPrime@ p]; lst (* Robert G. Wilson v, Jan 04 2016 *)
|
|
PROG
|
(PARI) forprime(p=3, , forprime(q=2, p-1, if(Mod(p, q^2)^(q-1)==1 && Mod(q, p^2)^(p-1)==1, print1(p, ", "); break({1}))))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,hard,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|