OFFSET
1,2
LINKS
Ray Chandler, Table of n, a(n) for n = 1..10000
Saltire Software, Pythagorean Triples
Eric Weisstein's World of Mathematics, Pythagorean Triple.
FORMULA
Arrange all primitive Pythagorean triples a, b, c by value of hypotenuse c, then by long leg b; for n-th value of c, sequence gives radius of largest inscribed circle, (a+b-c)/2.
MAPLE
with(NumberTheory):
f := proc(N)
local k, l, m, x, y, z;
l := [];
for m from 2 while m^2 + 1 <= N do
for k from 1 to m - 1 do
if igcd(m, k) > 1 then next end if;
if (m - k) mod 2 = 0 then next end if;
z := m^2 + k^2;
if z > N then break end if;
x := m^2 - k^2; y := 2*m*k; l := [op(l), [z, max(x, y), min(x, y), m, k]]
end do
end do;
sort(l, (u, v) -> u[1] < v[1] or (u[1] = v[1] and u[2] < v[2]));
end proc:
local a, i, k, l, m, r;
l := f(N); a := [];
for i to nops(l) do
m := l[i][4]; k := l[i][5]; r := m*k - k^2; a := [op(a), r]
end do;
a
end proc:
A014498(500); # Felix Huber, Jul 07 2026
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Ralph Peterson (ralphp(AT)library.nrl.navy.mil)
EXTENSIONS
More terms from Asher Auel May 05 2000
Extended by Ray Chandler, Mar 09 2004
STATUS
approved
