OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 13 is a term because 13 and 13^2 + 2*13 + 4 = 199 are prime.
MAPLE
select(p -> isprime(p^2+2*p+4), [seq(ithprime(i), i=1..1000)]);
PROG
(Python)
from sympy import isprime
for p in range (2, 3700):
if isprime(p) and isprime(p**2 + 2*p + 4):
print (p, end=", ") # Karl-Heinz Hofmann, Jan 11 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Jan 10 2022
STATUS
approved