OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
23 is in the sequence because 23, 2*23+1 = 47 and 2*23+13 = 59 are all prime.
MATHEMATICA
Select[Prime[Range[10000]], PrimeQ[2 # + 1] && PrimeQ[2 # + 13] &] (* Vincenzo Librandi, Feb 15 2014 *)
Select[Prime[Range[1000]], AllTrue[2#+{1, 13}, PrimeQ]&] (* Harvey P. Dale, Jun 27 2023 *)
PROG
(PARI) s=[]; forprime(p=2, 10000, if(isprime(2*p+1) && isprime(2*p+13), s=concat(s, p))); s
(Magma) [p: p in PrimesUpTo(9200) | IsPrime(2*p+1) and IsPrime(2*p+13)]; // Vincenzo Librandi, Feb 15 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Feb 13 2014
STATUS
approved