Suppose $m$ and $n$ are integers (without loss of generality, nonnegative integers) satisfying the relations $m|n^2-1$ and $n|m^2-1$. What can be said about $m$ and $n$?
One family of solutions can be constructed as follows: let $k$ be an arbitrary integer, and consider the sequence $a^k_n$ defined recursively by $$ a^k_0=0\\ a^k_1=1\\ a^k_j=ka^k_{j-1}-a^k_{j-2} $$
This sequence satisfies the "Cassini-like" identity $$ a^k_{j+1}a^k_{j-1}=(a^k_j)^2-1 $$ which means that we can take $(m,n)=(a^k_j,a^k_{j+1})$ for any $k,j$.
Note that this family includes the obvious solutions $n=m+1$, which are consecutive values of the sequence $a^2_j$. It also includes the slightly-less-obvious solutions $n=m^2-1$, as $a^k_2=k$ and $a^k_3=k^2-1$.
A computer search shows that every pair $(m,n)$ satisfying these relations with $m,n<10^5$ is of this form.
Is every pair which satisfies these relations of this form?
(This question arose out of an early attempt to solve this problem, before I found some stronger divisibility conditions which led to a complete solution.)