There is no real need to look for substitutions- it's pretty simple as is. This is an "exact" equation.
\frac{dy}{dx}= \frac{1- 2y- 4x}{1+ y+ 2x}
(1+ y+ 2x)dy= (1- 2y- 4x)dx
(1+ y+ 2x)dy+ (4x+ 2y- 1)dx= 0
Now we see that (1+ y+ 2x)_x= 2= (4x+ 2y- 1)_y so this is an "exact differential"- there exist a function F(x,y) such that dF= F_xdx+ F_ydy= (4x+2y- 1)dx+ (1+ y+ 2x)dx.
From F_x= 4x+ 2y- 1 so that F= 2x^2+ 2xy- g(y).
(Since F_x is a differentiation with respect to x, treating y as a constant, when we "reverse" that the "constant of integration" may be a function of y.)
From that, F_y= 2x- g'(y)= 1+ y+ 2x so that -g'(y)= 1+ y and g'(y)= -1- y. g(y)= -y- \frac{1}{2}y^2+ C.
So F(x,y)= 2x^2+ 2xy- (-y- \frac{1}{2}y^2+ C)= 2x^2+ 2xy+ y+ \frac{1}{2}y^2+ c.
Saying that "dF= 0" gives F(x, y)= 2x^2+ 2xy- (-y- \frac{1}{2}y^2+ C)= 2x^2+ 2xy+ y+ \frac{1}{2}y^2+ c equals a constant. Including the "c" in that constant,
2x^2+ 2xy+ y+ \frac{1}{2}y^2= C.