OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 100000: # for triples with area <= N
R:= NULL:
for n from 1 while (2*n+1)*(n+1)*n <= N do
for m from n+1 by 2 while (m^2 - n^2)*m*n <= N do
if igcd(m, n) = 1 then
R:= R, [m^2-n^2, 2*m*n, m^2+n^2, (m^2-n^2)*m*n]
fi
od od:
R:= sort(R, (s, t) -> s[4] < t[4] or (s[4] = t[4] and s[3] < t[3])):
seq(min(t[1], t[2]), t=R); # Robert Israel, Dec 30 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Aug 17 2006
EXTENSIONS
a(43)=40 inserted by Ray Chandler, Jun 26 2017
STATUS
approved