Solve[x == y + y^2, y, MaxExtraConditions -> Automatic]