OFFSET
1,1
COMMENTS
Numbers of the form m^2 + n^2 where 1 <= n < m, m and n are coprime with one odd and one even, and 2*m^2 + m*n + n^2 and 2*m^2 + 3*m*n - n^2 are prime.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 25 is a term because 25 is the hypotenuse of the primitive Pythagorean triangle with sides 24, 7, 25 having perimeter 24 + 7 + 25 = 56 and inradius (24 + 7 - 25)/2 = 3 and both 56 - 3 = 53 and 56 + 3 = 59 are prime.
a(14) = 265 is the first term that corresponds to more than one primitive Pythagorean triangle: sides (23, 264, 265) and (247, 96, 265).
MAPLE
N:= 10000: # for terms <= N
R:= NULL:
for m from 1 while m^2 < N do
for n from 1 + (m mod 2) to m by 2 while m^2 + n^2 <= N do
if igcd(m, n) > 1 then next fi;
if isprime(2*m^2 + m*n + n^2) and isprime(2*m^2 + 3*m*n - n^2) then
R:= R, m^2+n^2;
fi od od:
sort(convert({R}, list));
CROSSREFS
KEYWORD
nonn
AUTHOR
Will Gosnell and Robert Israel, Dec 01 2025
STATUS
approved
