|
|
A319025
|
|
Primes p such that W_p == 2 (mod p), where W_p = A007619(n) and p = prime(n).
|
|
0
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
These are the members of René Gy's set W_2 (cf. Gy, 2018).
The sequence is complete to 2*10^13, with the higher terms coming from a list of primes with small Wilson quotients in the article by Costa, Gerbicz, and Harvey. - John Blythe Dobson, Jan 05 2021
|
|
LINKS
|
|
|
PROG
|
(PARI) forprime(p=1, , if(Mod(((p-1)!+1)/p, p)==2, print1(p, ", ")))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,hard,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|