Solución. La desigualdad entre las medias aritmética y geométrica nos dice que
\begin{align*}
\frac{3^{x^2-x-y}+3^{y^2-y-z}+3^{z^2-z-x}}{3}&\geq\sqrt[3]{3^{x^2-x-y}\cdot 3^{y^2-y-z}\cdot 3^{z^2-z-x}}\\
&=3^{\frac{x^2-2x+y^2-2y+z^2-2z}{3}}=3^{\frac{(x-1)^2+(y-1)^2+(z-1)^2-3}{3}}\\
&=\tfrac{1}{3}\cdot 3^{\frac{(x-1)^2+(y-1)^2+(z-1)^2}{3}}\geq \frac{1}{3}.
\end{align*}
Por lo tanto, $3^{x^2-x-y}+3^{y^2-y-z}+3^{z^2-z-x}\geq 1$ para todo $x,y,z\in\mathbb{R}$ y, si la igualdad se alcanza, tiene que ser $(x-1)^2+(y-1)^2+(z-1)^2=0$, es decir, $x=y=z=1$. Como $x=y=z=1$ verifica la ecuación del enunciado, deducimos que esta es la única solución.
Nota. Las exponenciales pueden ocultar la aplicación de la desigualdad entre las medias aritmética-geométrica, pero una solución similar se tiene aplicando la desigualdade de Jensen a la función convexa $f(t)=3^t$. ¿Sabrías escribir los detalles?