OFFSET
1,1
COMMENTS
Tested up to n = 4.09*10^8.
LINKS
Eric Weisstein's World of Mathematics, Pythagorean Quadruple
EXAMPLE
If n = 6336 then n' = 23808, n'' = 103936 and sqrt(n^2 + n'^2 + n''^2) = 106816.
MAPLE
with(numtheory): P:= proc(q) local a1, a2, n, p;
for n from 2 to q do a1:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
a2:=a1*add(op(2, p)/op(1, p), p=ifactors(a1)[2]);
if type(sqrt(n^2+a1^2+a2^2), integer) then print(n);
fi; od; end: P(10^10);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, Oct 25 2013
EXTENSIONS
a(16)-a(18) from Giovanni Resta, Oct 25 2013
a(19) from Ray Chandler, Dec 22 2016
a(20) from Ray Chandler, Dec 31 2016
a(21) from Ray Chandler, Jan 05 2017
a(22) from Ray Chandler, Jan 09 2017
STATUS
approved