OFFSET
1,1
COMMENTS
It is not known if there are infinitely many Sophie Germain pairs with this property.
The sequence is infinite under Dickson's conjecture. Aside from a(1) = 5, all terms are 29 or 179 mod 210. - Charles R Greathouse IV, Feb 05 2014
LINKS
Abhiram R Devesh and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 135 terms from Devesh)
Eric Weisstein's World of Mathematics, Sophie Germain Prime
Wikipedia, Sophie Germain prime
EXAMPLE
a(1): p = 5; (2*p)+1 = 11
Prime triples (5,7,13);(11,13,19)
a(2): p = 29; (2*p)+1=59
Prime triples (29,31,37);(59,61,67)
MATHEMATICA
sgpQ[n_]:=Module[{sg=2n+1}, AllTrue[Flatten[{sg+{0, 2, 8}, n+{2, 8}}], PrimeQ]]; Select[Prime[ Range[ 300000]], sgpQ] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Feb 02 2016 *)
PROG
(Python)
from sympy import isprime
print([5] + [n for m in range(29, 10**8, 210) for n in (m, m+150) if isprime(n) and isprime(n+2) and isprime(n+8) and isprime(2*n+1) and isprime(2*n+3) and isprime(2*n+9)]) # David Radcliffe, Aug 07 2025
(PARI) is(n)=isprime(n) && isprime(n+2) && isprime(n+8) && isprime(2*n+1) && isprime(2*n+3) && isprime(2*n+9) \\ Charles R Greathouse IV, Feb 05 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Abhiram R Devesh, Feb 04 2014
STATUS
approved
