OFFSET
1,1
COMMENTS
Except for 3, all terms end in 9.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 379 is a term because 379, 2*379-1 = 757 and (2*379-1)^2+(2*379)^2 = 1147613 are prime.
MAPLE
filter:= proc(p) isprime(p) and isprime(2*p-1) and isprime(8*p^2-4*p+1) end proc:
select(filter, [3, seq(i, i=9..50000, 10)]);
PROG
(Python)
from sympy import isprime, primerange
def ok(p): return isprime(2*p-1) and isprime((2*p-1)**2 + (2*p)**2)
def aupto(limit): return list(filter(ok, primerange(2, limit+1)))
print(aupto(44450)) # Michael S. Branicky, Aug 20 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Aug 20 2021
STATUS
approved