OFFSET
1,1
COMMENTS
The entries 600 and 780 occur in two solutions each; the others entries uniquely.
MAPLE
a := []; for n from 2 to 1000 do for m from 2 to n do if (numer(simplify(1/n^2+1/m^2))=1) then if type(sqrt(denom(simplify(1/n^2+1/m^2))), integer)=true then print(n, m, ifactor(n), ifactor(m), ifactor(denom(simplify(1/n^2+1/m^2))), sqrt(denom(simplify(1/n^2+1/m^2)))); a := [op(a), n]; fi fi od od; print(a);
CROSSREFS
KEYWORD
nonn
AUTHOR
Len Smiley, Dec 01 2001
EXTENSIONS
Offset corrected by Sean A. Irvine, Sep 09 2023
STATUS
approved