OFFSET
1,1
COMMENTS
Prime p>2 is in this sequence iff (p-1)/2 is in A015950. - Max Alekseyev, Dec 26 2017
LINKS
Max Alekseyev, Table of n, a(n) for n = 1..177
MATHEMATICA
Select[Prime[Range[1000]], Mod[2^# + 2, # - 1] == 0 &]
PROG
(PARI)
N=10^9;
default(primelimit, N);
forprime(p=2, N, if (-2==Mod(2, p-1)^p, print1(p, ", ")));
/* Joerg Arndt, Feb 06 2013 */
(Python)
from sympy import primerange
A211349_list = [p for p in primerange(1, 10**6) if p == 2 or pow(2, p, p-1) == p-3] # Chai Wah Wu, Mar 25 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Philip A. Hoskins, Feb 06 2013
EXTENSIONS
a(19)-a(47) from Giovanni Resta, Feb 10 2013
a(48)-a(177) from Max Alekseyev, Jan 06 2018
STATUS
approved