OFFSET
1,1
MATHEMATICA
s[p_] := Solve[1 < r < s < 2r && p == r s + s^2 && GCD[r, s] == 1, {r, s}, Integers];
Reap[For[p = 1, p <= 2000, p++, sp = s[p]; If[sp =!= {}, Print[r /. sp[[1]] ]; Sow[r /. sp[[1]]]]]][[2, 1]] (* Jean-François Alcover, Mar 06 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, May 04 2005
EXTENSIONS
Extended by Ray Chandler, May 09 2005
STATUS
approved