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