OFFSET
1,1
PROG
(PARI) g(n)= { for(x=1, n, for(y=x, n, z=sqrt(x^2+y^2); s=x*y*z/(z^2+x*y); if(s==floor(s)&&x%37, print1(floor(x)", ") ) ) ) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, Nov 17 2006
STATUS
approved