login
A323705
Odd primes p such that ((p+1)/2)^(p-1) == 1 (mod p^2).
0
1897121, 52368101, 126233057
OFFSET
1,1
COMMENTS
Primes p such that the Fermat quotient of p base 2 (A007663) is congruent to -1 modulo p. In other words, 2^(p-1) == 1-p (mod p^2). - Max Alekseyev, May 09 2026
The corresponding sequence when (p+1) in the congruence is replaced with (p-1) is A125854.
a(4) > 4262937421 if it exists.
MATHEMATICA
Select[Prime[Range[718*10^4]], PowerMod[(#+1)/2, #-1, #^2]==1&] (* Harvey P. Dale, May 14 2021 *)
PROG
(PARI) forprime(p=3, , my(x=(p+1)/2); if(Mod(x, p^2)^(p-1)==1, print1(p, ", ")))
CROSSREFS
KEYWORD
nonn,hard,bref,more
AUTHOR
Felix Fröhlich, Jan 24 2019
STATUS
approved