Let me provide a method, for part a, applied only for finding the impulse response $h(t)$ of an LTI system characterised by an LCCDE of the form $ \sum_{k=0}^{N}{ a_k {{d^k y(t)}\over {dt^k}}} = \sum_{k=0}^{M}{ b_k {{d^k x(t)}\over {dt^k}}}$ by using the classical time domain approach. A method which is generally ignored.
Note that the input $x(t) = \delta(t)$ is formally a problematic signal (function) and therefore mathematically oriented classical books on ODEs, do generally avoid any discussion of such functions and their solutions, unless the scope of the book specifically includes generalised functions (distributions), usually to be used in some engineering or physical fields.
The procedure: Consider an LTI system which is causal with initial rest conditions. We say that when the input to this system is $x(t)=\delta(t)$ an impulse, then its output $y(t)=h(t)$ is the impulse-response of the system. We find this solution $h(t)$ in two steps, by breaking the LCCDE into two parts as inspired by a serial (cascade) implementation of two LTI systems corresponding to the right and left sides of the LCCDE.
I/O relationship of Part-I is given by the LCCDE: $$\sum_{k=0}^{N}{ a_k {{d^k y(t)}\over {dt^k}}} = x(t)$$ which represents the first stage with its input $x(t)=\delta(t)$ and we denote its solution (stage1-output) as $h_0(t)$ which is actually the impulse response of Part-I.
And the I/O relationship of the Part-II is given by the equation: $$ y(t) = \sum_{k=0}^{M}{b_k{{d^k x(t)}\over {dt^k}}} $$ which requires nothing but simple summation of its input, $x(t) = h_0(t)$ ,and its derivatives to compute the output as $$h(t) = \sum_{k=0}^{M}{b_k {{d^k h_0(t)}\over {dt^k}}}$$ Therefore we need to find $h_0(t)$ of the Part-I to simply compute the impulse response $h(t)$ of the complete system.
To find the solution of part-I: $\sum_{k=0}^{N}{ a_k {{d^k y(t)}\over {dt^k}}} = x(t)$ whose solution $y(t)$ is what we call $h_0(t)$ , we we make the following observation that $y(t)$ is composed of two parts as $$y(t) = y_{h}(t) + y_{p}(t)$$ where $y_{h}(t)$ is the homogeneous solution corresponding to $x(t)=0$ for all t, and $y_{p}(t)$ is the particular solution corresponding to input $x(t)$ for $ t > 0$, in particular as what remains as input when time goes to infinity.
Now, we can see that for the input $x(t)=\delta(t)$ the particular solution $y_p(t)$ for $t>0$ is zero as the input is zero for $t > 0$. Then the output simplifies to $h_0(t) = y(t) = y_h(t)$
Now, we have to find the homogeneous part of $y(t)$ as the solution of the equation $$\sum_{k=0}^{N}{ a_k {{d^k y(t)}\over {dt^k}}} = 0$$
following the usual procedure by finding the roots of its characteristic polynomial in complex-s: $$p(s) = \sum_{k=0}^{N}{a_k s^k}$$
which is obtained when we simply insert $y_h(t) = Ke^{st}$ as a solution canditade to the homogeneous differential equation $\sum_{k=0}^{N}{ a_k {{d^k y(t)}\over {dt^k}}} = 0$. The solution is then hypothesised to be $$y_h(t) = \sum_{k=0}^{N}{ A_k e^{st}}$$ for the case when all the roots are distinct. If they are not, you know what to do right ;)
The solution therefore now requires to determine the N unknown coefficients $A_k$ from N-inital conditions, at time t=0 for the framework of this method, possesed by the output $y(t)$ and its derivates up to order N-1.
Remember that this system was defined to be under inital rest with complete zero inital conditions prior to the application of any excitation. Hence if we were to solve the homogeneous solution $y_h(t)$ for an ordinary (formally valid) input signal $x(t)$ to the system, the homogeneous solution will be identically zero, as the initial conditions would all be. (all the coefficients $A_k$ of the exponential terms turn out to be zero). But the magic of the generalised input $\delta(t)$ is that it will set one of those inital conditions to non-zero, thereby, enabling a non-zero homogeneous response to exist, which will become the solution of the part-I as well.
It can be shown that the inital conditions of the part-I due to the excitation $\delta(t)$ will be set as $y(t)=0 ~,~ y(t)' = 0 ~,~ y(t)''=0 ~,~... ~,~ y(t)^{N-1}=1/a_N$ all at $t=0$.
From this information you can proceed to solve the algebraic equaitons as usual and find $h_0(t)$ from which you can compute the impulse response $h(t)$ for all $t > 0$. What happens for all $t <0 $? As we stated, this system was LTI and causal with initial rest condition. For such a system if the input $x(t)$ is zero for all $t<0$ then that output is also zero for all $t<0$ which is achieved by multiplying $h_0(t)$ via $u(t)$ to represent the solution of Part-I for all t.
It is instructive to apply this procedure to your example:
$$y'' + 6 y' + 8y = 2x$$
Lets first solve Part-I: $$y'' + 6 y' + 8y = x$$
to find $h_0(t)$ when the input $x(t)$ = $\delta(t)$. solve the characteristic equation: $s^2 + 6s + 8 = 0$ which yields $s=-4$ and $s=-2$ as solutions. A stable system, isn't it? The proposed homogeneous solution is therefore $$y_h(t) = A_1 e^{-2t} + A_2 e^{-4t}$$
For which we need 2 inital conditions to find $A_1$ and $A_2$.
These initial conditions are set by the impulse $\delta(t)$ to be: $y(0)=0$ and $y(0)'=1$. Using these we find them to be $A_1=1/2$ and $A_2=-1/2$ and this yields $$ h_0(t) = {1 \over 2} [ e^{-2t} - e^{-4t}] $$ for all $t > 0$. Also note that this system was LTI and causal with initial rest and input $x(t)=\delta(t)$ is zero for all $t<0$ implies that output is also zero for all $t<0$ which is simplified by using a unit step $u(t)$ function to format the output as: $$ h_0(t) = {1 \over 2} [ e^{-2t} - e^{-4t}]u(t)$$ for all t.
Finally apply this $h_0(t)$ as an input to Part-II whose output is simply $$y(t)= 2x(t)$$ to get the impulse response of the overall system as:
$$h(t) = [ e^{-2t} - e^{-4t}]u(t)$$
