OFFSET
1,1
LINKS
Abhiram R Devesh, Table of n, a(n) for n = 1..675
EXAMPLE
a(1): p = 419; (2*p)+1 = 839
Prime Pairs of the form (p,p**3+2): (419, 73560061) and (839, 590589721)
a(2): p = 761; (2*p)+1 = 1523
Prime Pairs (761, 440711083) and (1523, 3532642669)
PROG
(Python)
p1=2
n=2
count=0
while p1>2:
....## Generate the pair
....cc=[]
....cc.append(p1)
....for i in range(1, n):
........cc.append((2**(i)*p1+((2**i)-1)))
....## chain entries cubed + 2
....cc2=[(c*c*c)+2 for c in cc]
....## check if cc is a Sophie Germain Pair or not
....## pf.isp_list returns True or false for a given list of numbers
....## if they are prime or not
....##
....pcc=pf.isp_list(cc)
....pcc2=pf.isp_list(cc2)
....## Number of primes for cc
....npcc=pcc.count(True)
....## Number of primes for cc2
....npcc2=pcc2.count(True)
....if npcc==n and npcc2==n:
........print "For length ", n, " the series is : ", cc, " and ", cc2
....p1=pf.nextp(p1)
CROSSREFS
KEYWORD
nonn,hard
AUTHOR
Abhiram R Devesh, May 10 2014
STATUS
approved