This really is Vieta Jumping, I ought to have put that first line. I put about five pictures of the relevant hyperbolas at the end of this post.
Actual reference: In 1907, A. Hurwitz wrote Über eine Aufgabe der unbestimmten Analysis, in Archiv der Mathematik und Physik, pages 185-196. I give a summary and a little discussion at https://mathoverflow.net/questions/84927/conjecture-on-markov-hurwitz-diophantine-equation
We will be carefully examining the part $x,y > 0$ in the graph of
$$ x^2 - q x y + y^2 + x + y = 0, $$ especially what happens as we get near the origin. Note that for $q=3$ we get the integer point $(2,2),$ for $q=4$ we get the lattice point $(1,1),$ but for $q=5$ we seem to miss any integer points. It turns out that is all that needs to be proved, that for $q \geq 5$ we get no lattice points near the origin. Technique written out very well in Hurwitz (1907).
In the fourth image, we see how the Hurwitz fundamental solutions lie between $y = \frac{qx-1}{2}$ and $x = \frac{qy-1}{2}$ along the curve
$ x^2 - q x y + y^2 + x + y = 0, $ and that both $x,y < 1$ there. So that proves it.
Part of a careful proof is to confirm that the Hurwitz mappings below, also called Vieta Jumping, do not magically take us from the hyperbola branch in the first quadrant to the other branch, which passes through the origin and contains other integer solutions as well. It suffices to show that the largest value of $y$ on that branch is smaller than $1.$ This occurs, by implicit differentiation, where $2x-qy+1 = 0,$ or $$y = \frac{1}{q}(2x+1),$$ with $q \geq 3$ no matter what. When $x < 0,$ $2x+1 < 1,$ finally $ \frac{1}{q}(2x+1) < 1.$ That is, while there really is a solution at $(0,0),$ there are no integer solutions with $xy < 0.$
Hurwitz: given a solution with $x,y > 0$ integers, we get a solution that usually shrinks $x+y$ with either
$$ (x,y) \mapsto (qx-y-1,x), $$ or
$$ (x,y) \mapsto (y,qy-x-1). $$
I should add that Hurwitz did not bother about the order of the variables, he passed very quickly over the case of two variables. I wrote the "jumps" this way so as to keep on the same half-branch ($y \geq x$ or $x \geq y$) of the hyperbola in the first quadrant, so as to fit in with some presentations I found of Vieta jumping; I will leave them written this way, but it hardly matters.
Since $x,y>0,$ we know that both $qx-y-1,$ $qy-x-1$ are positive as well; there are no solutions with $xy<0,$ and when one of $x,y$ is $0,$ the other is $0$ or $-1.$
We keep applying these transformations to get a sequence of solutions, until we reach a fundamental solution, at the minimum of $x+y,$ i.e. $qx-y-1 \geq y$ and $qy-x-1 \geq x,$
or
$$ y \leq \frac{qx-1}{2} \; \; \mbox{AND} \; \; x \leq \frac{qy-1}{2}, $$
which are the extra pair of lines I drew in the final two pictures. We do not need to pay much attention to whether $x$ or $y$ is larger, really.
I added in a picture with my $q=9,$ the entire arc is even more firmly within the strict standard unit square, that is $x,y < 1.$
The $y$ coordinate at $x = (2y+1)/ q$ satisfies
$$ (q^2 - 4)y^2 - (2q+4)y - (q+1) = 0. $$
For $q=3,4$ the positive root exceeds $1.$ However, for $q \geq 5,$ the positive root is smaller than one. For $q=5,$ $y = \frac{14 + 10 \sqrt 7}{42} \approx 0.963,$ so true for $q=5.$
Then
$$ y < \frac{1 + \sqrt{q+2}}{q-2}, $$
and that inequality finishes the proof for $q \geq 6.$
Wednesday, September 2. I found a prettier way to to the final bound on the positive value of $y,$ or, rather, show that it is strictly below $1.$ We have, where the Hurwitz bounding line $x = (2y+1)/ q$ meets the hyperbola in the first quadrant,
$$ (q^2 - 4)y^2 - (2q+4)y - (q+1) = 0. $$ Subtract $1$ from both sides, we have
$$ (q^2 - 4)y^2 - (2q+4)y - (q+2) = -1. $$ Next divide through by $(q+2),$ for
$$ (q - 2)y^2 - 2y - 1 = \frac{-1}{q+2} < 0. $$
For the left hand side to be negative with $y$ positive, the quadratic formula says we need strict
$$ \color{blue}{ y < \frac{1 + \sqrt{q-1}}{q-2}}. $$
For $q=5$ this says $y<1,$ exactly what we need. For $q > 5$ the restriction is stronger still.
