1,1

444,1333,1405; 444,5467,5485; 444,12317,12325; 444,49283,49285;

429, 460, 629; 429, 700, 821; 429,10220,10229; 429,92020,92021;

420, 851, 949; 420,1189,1261; 420, 1739, 1789; 420,4891,4909; 420,11021,11029; 420,44099,44101;

385,552,673; 385,1488,1537; 385,2952,2977; 385,74112,74113;

372,925,997; 372,3835,3853; 372,8645,8653; 372,34595,34597;

348,805,877; 348,3355,3373; 348,7565,7573; 348,30275,30277;

336,377,505; 336,527,625; 336,3127,3145; 336,28223,28225;

315,572,653; 315,988,1037; 315,1972,1997; 315,49612,49613;

300,589,661; 300,2491,2509; 300,5621,5629; 300,22499,22501;

280,351,449; 280,759,809; 280,1209,1241; 280,19599,19601;

276,493,565; 276,2107,2125; 276,4757,4765; 276,19043,19045;

228,325,397; 228,1435,1453; 228,3245,3253; 228,12995,12997;

204,253,325; 204,1147,1165; 204,2597,2605; 204,10403,10405;...

aa=aaa=aaaa=1; s=""; For[a=1, a<=10^3, For[b=a+1, ((b+1)^2-b^2)<=a^2, c=(a^2+b^2)^0.5; If[c==Round[c]&&GCD[a, b]==1, If[a==aaaa, s=s<>ToString[a]<>", "]; aaaa=aaa; aaa=aa; aa=a; ]; b++ ]; a++ ]; s

Sequence in context: A159382 A090486 A145635 this_sequence A108877 A163393 A154518

Adjacent sequences: A138602 A138603 A138604 this_sequence A138606 A138607 A138608

nonn

Vladimir Orlovsky (4vladimir(AT)gmail.com), May 13 2008