My question is in bold below. Here is a summary of my working on this problem which is question 6 from the IMO 1988 paper:
For $\frac{a^2+b^2}{ab+1}=n$, where $a,b,n \in Z^+$
let $n=x^2$
then for all $x \in Z^+, (a,b)$ have infinite solutions such that
$\hspace{2cm}$$a_m = a_0,a_1,a_2,... = 0, x,x^3,x^5-x,x^7-2x^3,x^9-3x^5+x,...$
$\hspace{2cm}$and $b=a_{m-1}, a_{m+1}$ for all $a_m$.
e.g. when $x=2$
valid solutions include:
$\hspace{2cm}$$(2,0)$ and $(2,8)$
$\hspace{2cm}$$(8,2)$ and $(8,30)$
$\hspace{2cm}$$(30,8)$ and $(30, 112)$ etc.
when $x=3,$
$\hspace{2cm}$$(3,0)$ and $(3,27)$
$\hspace{2cm}$$(27,3)$ and $(27,240)$
$\hspace{2cm}$$(240,27)$ and $(240,2133)$ etc.
An observation is that $a_m$ is a list of infinite solutions. For all $x$ (for all positive integers), $n$ must be a square. Another way to state this would be to say $n=x$, so then $a_m= \sqrt{x},x^\frac{3}{2},...$ and so on, $\therefore$ not integer solutions.
I am hoping I can either prove $a_m$ is a list of all solutions, therefore n must be a square integer, or to find the $m^{th}$ term for $a_m$.
My main question is where does this sequence for $a$ come from? It is easily reproducible, and looks like some sort of binomial expansion, but an expansion of what? And could I possibly prove this is a list of all solutions?
Some more terms of $a_m$ for curiosity:
$x^{11}-4x^7+3x^3$
$x^{13}-5x^9+6x^5-x$
$x^{15}-6x^{11}+10x^7-4x^{3}$... a clear pattern emerges.
I have briefly read that this is connected to the geometry of conics. If anybody has a good reference to more reading on this it would be greatly appreciated.
Also is there any links with Vieta jumping and $a_m$?