Inspired by the egregious problem in IMO 1988, I simulated the integer solutions to the equation
$$ (ab + 1) \mid (a^{2} + b^{2}) \tag{*}$$
for $1 \leq a, b \leq 3000$ and conjectured that every solution arises as an adjacent pair of numbers in the following sequence
$$ a_{0} = 0, \quad a_{1} = m, \quad \text{and} \quad a_{n+2} = ma_{n+1} - a_{n} $$
for some $m \in \Bbb{Z}$. Indeed, you can check that any pair $(a, b) = (a_{n}, a_{n+1})$ gives rise to a solution to $\text{(*)}$ with
$$ \frac{a_{n+1}^{2} + a_{n}^{2}}{a_{n+1}a_{n} + 1} = m^{2}. $$
I was unable to prove this conjecture as I'm ham-handed at number theory. Can you help me prove or disprove this?