OFFSET
1,1
COMMENTS
See comments in A174825
c is the prime hypotenuse c, i. e. c of a primitive Pythagorean triple: a^2 + b^2 = c^2
REFERENCES
W. W. R. Ball, H. S. M. Coxeter: Mathematical Recreations and Essays, New York: Dover, 1987
L. E. Dickson, "Rational Right Triangles", ch. 4 in History of the Theory of numbers, vol. II, Dover Publications 2005
W. Sierpinski, Pythagorean Triangles, Mineola, NY, Dover Publications, Inc, 2003
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 29 with 29^2 = 20^2 + 21^2 and 29 and 29||20||21 = 292021 prime.
a(2) = 409 with 409^2 = 120^2 + 391^2 and 409 and 409||120||391 = 409120391 prime.
a(3) = 461 with 461^2 = 380^2 + 261^2 and 461 and 461||380||261 = 461380261 prime.
a(4) = 661 with 661^2 = 300^2 + 589^2 and 661 and 661||300||589 = 661300589 prime.
a(5) = 929 with 929^2 = 920^2 + 129^2 and 929 and 929||920||129 = 929920129 prime.
MAPLE
N:= 20000: # for terms <= N
S:= {}:
for x from 2 while x^2 < N do
for y from 3 - (x mod 2) to x-1 by 2 while x^2 + y^2 <= N do
if igcd(x, y) > 1 then next fi;
c:= x^2 + y^2;
if not isprime(c) then next fi;
a:= 2*x*y;
b:= x^2-y^2;
if isprime((c*10^(1+ilog10(a))+a)*10^(1+ilog10(b))+b) then
S:= S union {c}
fi
od od:
sort(convert(S, list)); # Robert Israel, Mar 06 2026
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Ulrich Krug (leuchtfeuer37(AT)gmx.de), Apr 01 2010
EXTENSIONS
More terms from Zak Seidov, Apr 04 2010
Edited by Robert Israel, Mar 07 2026
STATUS
approved
