Vamos a ver que esto implica que $a=b=c=0$. En efecto, si alguno de los tres números $a,b,c\in\mathbb{Z}$ es no nulo, entonces podríamos haber comenzado suponiendo que la solución $(a,b,c)$ es tal que la suma $a^2+b^2+c^2\gt 0$ es lo más pequeña posible (de entre todas las soluciones no nulas habrá una que cumpla esto), pero entonces $(x,y,z)$ es una solución con \[x^2+y^2+z^2=\frac{1}{9}(a^2+b^2+c^2)\lt a^2+b^2+c^2,\] en contradicción con el hecho de que $a^2+b^2+c^2$ es mínimo. Esto es lo que se llama técnica del descenso infinito o principio de minimalidad.
Nota. También se podría haber eliminado $y$ obteniendo la igualdad \[7(25x+31y)-31(3x+7y)=82x,\] y el razonamiento a partir de aquí es similar.
Ahora bien para $n\leq -2$, la función decrece desde el límite $\frac{4}{3}$ hasta $f(-2)=1$, lo que nos dice que para $n$ en este intervalo sólo tenemos la solución $(m,n)=(1,-2)$. Ahora comprobamos algunos valores: \begin{align*} f(-1)&=\frac{8}{7},&f(0)&\text{ no definido},&f(1)&=-8,\\ f(2)&=5,&f(3)&=\frac{8}{3},&f(4)&=\frac{17}{8}. \end{align*} Como $f(5)=\frac{104}{55}\lt 2$, deducimos que la función decrece desde este valor hasta el límite $\frac{4}{3}\gt 1$, luego no hay soluciones enteras para $n\geq 5$.
Hemos probado que $(1,-2)$, $(-8,1)$ y $(5,2)$ son las únicas soluciones.
Nota. Si sólo buscamos las soluciones positivas, hay otro truco que merece la pena comentar. Si $m=1$ ó $n=1$, entonces $\frac{1}{m}+\frac{1}{n}+\frac{1}{mn^2}\gt 1\gt\frac{3}{4}$, luego supondremos que $m\geq 2$ y $n\geq 2$ en lo que sigue. Comenzamos probando valores: