OFFSET
1,1
COMMENTS
LINKS
Felix Huber, Table of n, a(n) for n = 1..10001
Eric Weisstein's World of Mathematics, Pythagorean Triple
EXAMPLE
28 is in the sequence because 21^2 + 28^2 = 35^2 and 21 + 28 = 7^2.
MAPLE
# Calculates the first 10001 terms
A380073:=proc(M)
local i, m, p, q, r, v, w, L, F;
L:=[];
m:=M^2+2*M+2;
for p from 2 to M do
for q to p-1 do
if gcd(p, q)=1 and (is(p, even) or is(q, even)) then
r:=1;
for i in ifactors(p^2-q^2+2*p*q)[2] do
if is(i[2], odd) then
r:=r*i[1]
fi
od;
w:=r*(p^2+q^2);
if w<=m then
v:=r*max(p^2-q^2, 2*p*q);
L:=[op(L), seq([i^2*w, i^2*v], i=1..floor(sqrt(m/w)))]
fi
fi
od
od;
F:=[];
for i in sort(L) do
F:=[op(F), i[2]]
od;
return op(F)
end proc;
A380073(4330);
CROSSREFS
KEYWORD
nonn,new
AUTHOR
Felix Huber, Jan 18 2025
STATUS
approved