login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A167293
Long legs of Pythagorean triangles that are not divisible by any other long leg of a Pythagorean triangle.
0
4, 15, 21, 35, 55, 77, 91, 99, 117, 143, 153, 171, 187, 209, 221, 247, 253, 299, 323, 325, 377, 391, 403, 425, 437, 475, 493, 527, 551, 575, 589, 621, 629, 667, 697, 703, 713, 725, 775, 779, 783, 817, 837, 851, 899, 925, 943, 957, 989, 999, 1023, 1025, 1073
OFFSET
1,1
COMMENTS
All long legs of Pythagorean triangles (A009023) are multiples of these values, so these values can be thought of as "primes" of the sequence of long legs.
PROG
(PARI) llp = vector(60); np = 1; llp[np] = 4;
notdiv(k) = for(j=1, np, if(k%llp[j], 1, return(0))); return(1);
isLongLeg(n) = local(b); b=0; for(k=1, n-1, if(issquare(k^2+n^2), b=1)); return(b);
for(k=4, 1175, if(notdiv(k), if(isLongLeg(k), np+=1; llp[np]=k)))
for(n=1, 60, print1(llp[n], ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Gerald McGarvey, Nov 01 2009
EXTENSIONS
Comments and PARI program corrected by Gerald McGarvey, Nov 03 2009
STATUS
approved