Solución. Haciendo $x=y$, tenemos que $2xf(x)=f(x^2)+x^2-f(0)$ y, como el miembro de la derecha es par, el de la izquierda también debe serlo, de donde deducimos que $f$ es impar. Si consideramos la función impar $g(x)=f(x)-x$ y sustituimos $f(x)=g(x)+x$ en la ecuación original, llegamos a que $g$ cumple la ecuación funcional $g((x-y)^2)=g(x)^2+2xg(x)-2xg(y)$. Haciendo $x=0$ en esta última ecuación, tenemos que $g(y^2)=g(0)^2$ para todo $y\in\mathbb{R}$, es decir, $g$ es constante en el intervalo $[0,\infty)$ y, como es impar, es constante en todo $\mathbb{R}$, pongamos $g(x)=a$ para cierto $a\in\mathbb{R}$ y, en sustituyendo en la ecuación funcional de $g$, se cumple que $a=a^2$ luego $a=0$ ó $a=1$. Deshaciendo ahora el cambio de función, las únicas posibles soluciones para $f$ son $f(x)=x$ y $f(x)=x+1$ y es fácil comprobar que ambas cumplen la ecuación del enunciado.