Step 1: Express the equation in differential form
Rewriting:
\[
(2y^3 - 4xy + y) dx - (2y^2 + 1) dy = 0.
\]
Step 2: Check for exactness
Computing partial derivatives,
\[
M(x, y) = 2y^3 - 4xy + y, \quad N(x, y) = -(2y^2 + 1).
\]
Since \( \frac{\partial M}{\partial y} = \frac{\partial N}{\partial x} \), the equation is exact.
Step 3: Solve using integration
\[
\int (2y^3 - 4xy + y) dx = F(x, y).
\]
\[
\int (2y^2 + 1) dy = G(x, y).
\]
Solving, we get:
\[
4xy^2 + 2x = y^4 + y^2 + c.
\]
Thus, the correct answer is \( \boxed{4xy^2 + 2x = y^4 + y^2 + c} \).