Deje $\xi_n$ ser la solución positiva de la ecuación $$x=\frac{1}{x}+n\tag{1}\label{1}$$ La secuencia de $\{a_n\}_{n=1}^{\infty}$ se define como: $$a_n=\xi_n(\xi_n-1)\cdots(\xi_n-n+1)\tag{2}\label{2}$$ Ahora, dado que (1) implica $\xi_n^2=1+n \xi_n$, es seguro decir que el $a_n$ puede ser en última instancia, se expresa como una función lineal de la $\xi_n$.
En otras palabras, si ampliamos $a_n$ en términos de $\xi_n$, todos los de la $\xi^m,m>1$ términos puede ser eliminado y reemplazado con bajos grados. Y si repetimos este proceso, por fin tendría: $$a_n=b_n\xi_n+c_n$$ El problema es que estoy teniendo problemas en encontrar una fórmula general para $b_n$$c_n$. Traté de ampliar el polinomio (2) y la simplificación de acuerdo a (1), pero se convirtió en demasiado complejo. También traté de calcular el $b_n,c_n$ algunos $n$ con el fin de encontrar un patrón, pero no tuvo suerte. Aquí está el resultado de $n=1,...,10$: $$ \begin{array}{c|lcr} n & \xi_n & b_n & c_n \\ \hline 1 & \frac{\sqrt{5}+1}{2} & 1 & 0 \\ 2 & \sqrt{2}+1 & 1 & 1 \\ 3 & \frac{\sqrt{13}+3}{2} & 3 & 0 \\ 4 & \sqrt{5}+2 & 8 & 4 \\ 5 & \frac{\sqrt{29}+5}{2} & 35 & 0 \\ 6 & \sqrt{10}+3 & 162 & 54 \\ 7 & \frac{\sqrt{53}+7}{2} & 1001 & 0 \\ 8 & \sqrt{17}+4 & 6656 & 1664 \\ 9 & \frac{\sqrt{85}+9}{2} & 53865 & 0 \\ 10 & \sqrt{26}+5 & 467500 & 93500 \\ \end{array}$$ Bueno, al parecer, $c_n=0$ para valores impares de $n$. Pero estoy atrapado aquí por ahora y no se puede ir más allá. Todas las ideas o sugerencias se agradece.
Editar: Se me olvidó mencionar que si $c_n\ne 0$ ($n$ es incluso) a continuación, se asemeja $b_n=\frac{n}{2} c_n$. Aunque no tengo ni idea de por qué?