(Python)
from itertools import count, islice
from sympy import isprime, nextprime
from sympy.abc import x, y
from sympy.solvers.diophantine.diophantine import diop_quadratic
def A114370_gen(): # generator of terms
for l in count(1):
c = []
for m in range(1, 10):
k = m*(10**l-1)//9<<1
for a, b in diop_quadratic((x-y)*(x+y-1)-k):
if isprime(b) and a == nextprime(b):
c.append(b)
yield from sorted(c)
|