He explorado un poco por ordenador y mis resultados son los siguientes, presentados sin pruebas.
Por cada n>3 podemos generar una secuencia infinita de soluciones.
Definir una secuencia por a_{k}=(n-2)a_{k-1}-a_{k-2} con a_1=1 y a_2=n-1 . Es bastante sencillo demostrar por inducción que a_k^2 = a_{k-1}a_{k+1} + n . De esto se deduce que cualquier par de elementos consecutivos de la secuencia, x=a_k , y=a_{k+1} satisface las ecuaciones.
Así que para n=4 se obtiene 1, 3, 5, 7, 9, 11, 13,... .
Para n=5 se obtiene 1, 4, 11, 20, 76, 199, 521,... .
Para n=6 se obtiene 1, 5, 19, 71, 265, 989,... .
Para n=7 se obtiene 1, 6, 29, 139, 666,... .
etc.
La misma construcción funciona también para n<0 .
Así que para n=-1 se obtiene 1, -2, 5, -13, 34, ... .
Para n=-2 se obtiene 1, -3, 8, -21, 85, ... .
etc.
Para n=0 tienes las infinitas soluciones: x=k , y=k para cualquier k .
Para n=1 tienes las infinitas soluciones: x=k , y=1 para cualquier k .
Para n=3 hay una fórmula recursiva diferente para una secuencia infinita de soluciones.
(x_n,y_n)= (6,3) , (69,39) , (753,426) , (8214,4647) , (89601,50691) , ...
Estos satisfacen x_k=11x_{k-1}-x_{k-2} y y_k=11y_{k-1}-y_{k-2} .
n=2 es el caso más difícil. Sólo he encontrado las siguientes soluciones: (1,1) , (2,2) , (1262,94) y (3602578,58594) .