Dados dos números enteros relativamente primos $r,s,t$ Estoy buscando una solución completa ( es decir, parametrización entera o similar) para la ecuación diofantina $$ rX_1^2+sY_1^2+tZ_1^2=rX_2^2+sY_2^2+tZ_2^2. \tag{$ \bigstar $} $$ Sé que puedo juntar términos en un lado y aplicar la parametrización Bachet-Cauchy-Bezout en términos lineales pero esperaba que esta rueda ya estuviera resuelta.
EDICIÓN/ACLARACIÓN:
Utilizando la solución lineal completa clásica de Bachet-Cauchy-Bezout, podemos resolver la ecuación $$ r(x_1^2-x_2^2)+s(y_1^2-y_2^2)+t(z_1^2-z_2^2)=0, $$ para obtener tres ecuaciones lineales de la forma $$ sc_1+t(-b_1)+(x_2-x_1)(x_2+x_1)=0, $$ donde $a_1,b_1,c_1$ son números enteros arbitrarios. Aplicando de nuevo la solución completa a cada una de las tres nuevas ecuaciones se obtienen tres soluciones de la forma \begin{align} x_1 &= \frac{(sb_2-ta_2)^2+(sc_1-tb_1)}{2(sb_2-ta_2)} &&\text{and}& x_2 &= \frac{(sb_2-ta_2)^2-(sc_1-tb_1)}{2(sb_2-ta_2)}, \end{align} donde el $a_i,b_i,c_i$ son números enteros arbitrarios. Una vez igualados y resueltos, el resultado sería ostensiblemente la parametrización entera completa de $(\bigstar)$ .