Con la ecuación
$$x^3 + y^3 = x^2 + y^2 + 42xy \tag{1}\label{eq1A}$$
ya ha manejado los casos de $x = 0$ o $y = 0$ . Para los casos no nulos, tener
$$\gcd(x,y) = d, \; x = de, \; y = df, \; \gcd(e, f) = 1 \tag{2}\label{eq2A}$$
Entonces \eqref {eq1A} se convierte en
$$\begin{equation}\begin{aligned} (de)^3 + (df)^3 & = (de)^2 + (df)^2 + 42(de)(df) \\ d^3(e^3 + f^3) & = d^2(e^2 + f^2 + 42ef) \\ d(e + f)(e^2 - ef + f^2) & = e^2 + f^2 + 42ef \end{aligned}\end{equation}\tag{3}\label{eq3A}$$
Esto demuestra que $e^2 - ef + f^2$ debe dividir el lado derecho, así que usando esto y $e^2 + f^2 \equiv ef \pmod{e^2 - ef + f^2}$ , da
$$e^2 + f^2 + 42ef \equiv ef + 42ef \equiv 43ef \pmod{e^2 - ef + f^2} \tag{4}\label{eq4A}$$
Debido a $\gcd(e,f) = 1 \implies \gcd(ef, e^2 - ef + f^2) = 1$ entonces
$$e^2 - ef + f^2 \mid 43 \tag{5}\label{eq5A}$$
Desde $d$ es el $\gcd$ Será positivo. Si ambos $e$ y $f$ son negativos, entonces no hay solución (ya que el lado izquierdo de \eqref {eq1A} será negativo y el lado derecho será positivo). Por lo tanto, o bien $e$ y $f$ son ambos positivos o uno es negativo siendo el otro positivo.
Cualquiera de los dos casos da $e^2 - ef + f^2 \gt 0$ por lo que, desde \eqref {eq5A}, esto significa $e^2 - ef + f^2 = 1$ o $e^2 - ef + f^2 = 43$ . Por ser igual a $1$ Debemos tener $e$ y $f$ siendo positivo y, dado que $e^2 - 2ef + f^2 = (e - f)^2 = 1 - ef \ge 0 \implies ef = 1$ entonces $e = f = 1$ es la única solución posible. Esto da en \eqref {eq3A}
$$d(1 + 1)(1 - 1 + 1) = 1 + 1 + 42 \implies d = 22 \tag{6}\label{eq6A}$$
Esto muestra $(x, y) = (22, 22)$ es una solución, ya que wimi La pregunta del Sr. G. de la Cruz es la siguiente comentario estados. A continuación, consideremos el otro caso de
$$e^2 - ef + f^2 = 43 \implies e^2 - (f)e + (f^2 - 43) = 0 \tag{7}\label{eq7A}$$
Tratamiento de $f$ como si fuera constante, entonces se resuelve para $e$ utilizando la fórmula cuadrática, da un discriminante de $f^2 - 4(f^2 - 43) = 172 - 3f^2$ que debe ser $\ge 0$ así que $\lvert f \rvert \le 7$ y es un cuadrado perfecto. Probando las distintas posibilidades de números enteros se observa que $f \in \{\pm 1, \pm 6, \pm 7\}$ . Utilizando estos valores para comprobar $e = \frac{f \pm \sqrt{172 - 3f^2}}{2}$ entonces $d$ de \eqref {eq3A} y $(x,y)$ de \eqref {eq2A} (que dejaré que tú hagas), da como resultado las restantes soluciones no nulas de $(x, y)$ son $\{(-6, 1)$ , $(1, -6)$ , $(1, 7)$ , $(7, 1)\}$ , ya que Christian Blatter La pregunta del Sr. G. de la Cruz es la siguiente comentario indica.