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