A106500 Values r associated with A106499.
2, 3, 3, 4, 5, 4, 5, 6, 5, 7, 5, 7, 8, 7, 6, 9, 7, 8, 9, 7, 10, 7, 8, 11, 9, 10, 11, 9, 12, 8, 11, 13, 11, 9, 13, 11, 14, 9, 10, 13, 11, 12, 15, 13, 11, 14, 15, 10, 13, 16, 11, 12, 13, 11, 14, 17, 15, 13, 16, 11, 17, 18, 13, 17, 13, 16, 19, 17, 12, 15, 13, 19, 14, 17, 20, 15, 13, 16, 19
Offset: 1
Keywords
Programs
-
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 *)
Extensions
Extended by Ray Chandler, May 09 2005