OFFSET
1,1
COMMENTS
Solutions to 2*m^2 + m*n + n^2 = k^2 where m and n are coprime and one of m and n is even.
LINKS
Robert Israel, Table of n, a(n) for n = 1..2500
EXAMPLE
a(3) = 23 is a term because with m = 16 and n = 1, the primitive right triangle with sides a = m^2 - n^2 = 255, b = 2*m*n = 32 and c = m^2 + n^2 = 257 has perimeter a+b+c = 544, inradius (a+b-c)/2 = 15 and perimeter - inradius = 544 - 15 = 529 = 23^2.
MAPLE
N:= 1000: # for terms <= N
R:= {}:
for m from 1 while 2*m^2 < N^2 do
for n from 1+(m mod 2) to m-1 by 2 do
v:= 2*m^2 + m*n + n^2; if v > N^2 then break fi;
if igcd(m, n) = 1 and issqr(v) then
a:= m^2 - n^2; b:= 2*m*n; c:= m^2 + n^2; y:= sqrt(v);
R:= R union {y}
fi
od od:
sort(convert(R, list));
CROSSREFS
KEYWORD
nonn
AUTHOR
Will Gosnell and Robert Israel, Nov 20 2025
STATUS
approved
