editing
approved
editing
approved
a(20)=145 corrected by Philippe Guglielmetti, Dec 14 2016
a(43)=401 inserted by Michael Somos, Dec 21 2016
approved
editing
editing
approved
(PARI) v=vector(M=10^4); for(a=1, M, v[a] = []; fordiv(2*a, x, if(x<(y=2*a/x) && issquare(x^2+y^2, &z) && 1==gcd([x, y, z]), v[a] = concat(z, v[a])))); concat(v) /* Michael Somos, Dec 21 2016 */
approved
editing
editing
approved
5, 13, 17, 25, 41, 29, 37, 61, 65, 85, 53, 113, 65, 101, 145, 73, 85, 89, 181, 145, 221, 97, 125, 109, 197, 265, 149, 313, 257, 173, 137, 365, 185, 157, 325, 421, 229, 169, 481, 205, 185, 193, 401, 269, 545, 293
approved
editing
editing
approved
v=vector(M=10^4); for(a=1, M, v[a] = []; fordiv(2*a, x, if(x<(y=2*a/x) && issquare(x^2+y^2, &z) && 1==gcd([x, y, z]), v[a] = concat(z, v[a])))); concat(v) /* Michael Somos, Dec 21 2016 */
proposed
editing
editing
proposed
corrected 20th term of the sequence, which is 145, not 142. the triple is (24, 143, 145), with area 1716, which is between 1710 (19, 180, 181) and 2310 (21, 220, 221). BTW there is no primitive triple with hypotenuse=142 (A020882)