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
35 is in the sequence because 21^2 + 28^2 = 35^2 and 21 + 28 = 7^2.
206125 is twice in the sequence because 31525^2 + 203700^2 = 206125^2 and 31525 + 203700 = 485^2 as well as 94588^2 + 183141^2 = 206125^2 and 94588 + 183141 = 527^2.
MAPLE
# Calculates the first 10001 terms
A380072:=proc(M)
local i, m, p, q, r, w, L;
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
L:=[op(L), seq(i^2*w, i=1..floor(sqrt(m/w)))]
fi
fi
od
od;
return op(sort(L))
end proc;
A380072(4330);
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Huber, Jan 18 2025
STATUS
approved
