DSolve[-((-(x*f[x^2 + a*y[x]^2]) + y[x])/(x + a*f[x^2 + a*y[x]^2]*y[x])) + Derivative[1][y][x] == 0,y[x],x]