OFFSET
1,1
MATHEMATICA
s = {}; p = 2; Do[q = NextPrime[p]; If[BitAnd[p, q] == 1, AppendTo[s, p]]; p = q, {10^5}]; s (* Amiram Eldar, Apr 16 2020 *)
Select[ NextPrime[ 2^Range[82], -1], BitAnd[#, NextPrime@ #] == 1 &] (* Giovanni Resta, Apr 16 2020 *)
PROG
(PARI) isok(p) = isprime(p) && (bitand(p, nextprime(p+1)) == 1);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Marcus, Apr 16 2020
EXTENSIONS
a(9)-a(10) from Amiram Eldar, Apr 16 2020
a(11)-a(19) from Giovanni Resta, Apr 16 2020
STATUS
approved