Proving a Pellian connection in the divisibility condition $(a^2+b^2+1) \mid 2(2ab+1)$ I'm trying to prove that all integer solutions $a > b \ge 0$ to the divisibility condition in the title, namely
$$(a^2+b^2+1) \mid 2(2ab+1),$$
are given by
$$(a,b)=(1,0),(4,1),(15,4),(56,15),(209,56),\dots$$
where $a_n$ is in http://oeis.org/A001353, and $b_{n+1}=a_n$. This sequence is related to the Pell equation $x^2-3y^2=1$, in that the sequence $0,1,4,15,56,\dots$ give the $y$ values in the solutions $(x,y)$.
Related MSE posts of possible interest include Proving there are an infinite number of pairs of positive integers $(m,n)$ such that $\frac{m+1}{n}+\frac{n+1}{m}$ is a positive integer and Conjecture on integer solutions to the equation $ (ab + 1) \mid (a^{2}+b^{2})$. Vieta jumping doesn't seem to apply (since it's best used to derive a contradiction, or prove a constant value). Any help/pointers would be appreciated.
 A: At a quick glance, Vieta jumping still applies. Note that your consecutive solution pairs share a common term e.g. $(1,4) \rightarrow (4, 15) $. This strongly hints that we do apply Vieta's, and in fact the standard Vieta's jumping argument leads to your conclusion.

For $ a > b \geq 0$, we have $ a^2 + b^2 + 1 > 2ab + 1$. Hence, $ \frac{ 2(2ab+1) } { a^2 + b^2 + 1} < 2 $, which tells us that divisibility is satisfied only when 
$$ a^2 + b^2 + 1 = 2 (2ab+1). $$
Now, observe that if $(a,b)$ is a solution, then by Vieta's jumping, so is $(4a-b, b)$ and $(b, 4b-a)$. 

You can verify that this leads to the same solution set of Pell's equation, because they both satisfy the same underlying recurrence of consecutive terms, namely $ y_{n+2} = 4 y_{n+1} - y_n$.
A: Your link does not like me. Solutions for hard to find such equations. So easy.
the equation:  $\frac{m+1}{n}+\frac{n+1}{m}=a$
Can be solved using the equation Pell:
$p^2-(a^2-4)s^2=1$
Solutions have the form:
$n=2(p-(a+2)s)s$
$m=-2(p+(a+2)s)s$
Solutions have the form:
$n=\frac{2p(p+(a-2)s)}{a-2}$
$m=\frac{2p(p-(a-2)s)}{a-2}$
If there is a solution of the following equation Pell: $p^2-(a^2-4)s^2=4$
We use its solutions and the formula is:
$n=\frac{p-(a-2)s+2}{2(a-2)}$
$m=\frac{p+(a-2)s+2}{2(a-2)}$
For your case.
To solve the equation in general form: $\frac{X^2+aX+Y^2+bY+c}{XY}=j$
In the case where a square:
$t=\sqrt{(b+a)^2+4c(j-2)}$
We use the solutions of Pell's equation:
$p^2-(j^2-4)s^2=1$
Solutions can be written.
$X=\frac{(b+a\pm{t})}{2(j-2)}p^2+(t\mp{(b-a)})ps-\frac{(b(3j-2)+a(6-j)\mp{(j+2)t})}{2}s^2$
$Y=\frac{(b+a\pm{t})}{2(j-2)}p^2+(-t\mp{(b-a)})ps-\frac{(b(6-j)+a(3j-2)\mp{(j+2)t})}{2}s^2$
It must be remembered that the numbers $p,s$ may be any mark.
