OFFSET
1,1
COMMENTS
All terms == 1 (mod 3).
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 85 is a term because (3*85-3)^2 + 85^2 = 70729 and (3*85)^2 + 86^2 = 72421 are both primes.
MAPLE
filter:= t -> isprime((3*t-3)^2 + t^2) and isprime((3*t)^2 + (t+1)^2):
select(filter, [$1..10000]);
MATHEMATICA
Select[Range[1, 5000, 3], PrimeQ[2*#*(5*#-9) + 9] && PrimeQ[2*#*(5*#+1) + 1] &] (* Paolo Xausa, Jan 15 2026 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Will Gosnell and Robert Israel, Jan 13 2026
STATUS
approved
