OFFSET
1,2
COMMENTS
Of course m = n^2 + 1 is the sum of two squares, by definition. Here there should be just one other way to write m as a different sum of two squares.
Let p and q be primes of the form 1+4k. Then n^2+1 must be pq or 2pq. - T. D. Noe, May 27 2008
LINKS
EXAMPLE
E.g., 111^2 + 1 = 21^2 + 109^2 only.
MATHEMATICA
ok[1] = True; ok[n_] := Length[ {ToRules[ Reduce[ 1 < x <= y && n^2 + 1 == x^2 + y^2, {x, y}, Integers] ] } ] == 1; Select[ Range[136], ok] (* Jean-François Alcover, Feb 16 2012 *)
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Patrick De Geest, Sep 15 1999
EXTENSIONS
Better definition from T. D. Noe, May 27 2008
STATUS
approved