The convolution is defined as
\begin{align*}
x(t) * h(t)
&=\int_{-\infty}^{\infty}x(\tau)\,h(t-\tau)\,d\tau\\
&=\int_{-\infty}^{\infty}[\delta(\tau-3) - 2\delta(\tau-4)]\,[\delta(t-\tau-1) + \delta(t-\tau-3)]\,d\tau\\
&=\int_{-\infty}^{\infty}\delta(\tau-3)\,\delta(t-\tau-1)\,d\tau\\
&+\int_{-\infty}^{\infty}\delta(\tau-3)\,\delta(t-\tau-3)\,d\tau\\
&-2\int_{-\infty}^{\infty}\delta(\tau-4)\,\delta(t-\tau-1)\,d\tau\\
&-2\int_{-\infty}^{\infty}\delta(\tau-4)\,\delta(t-\tau-3)\,d\tau.
\end{align*}
Now the way the Dirac delta function works, it "picks out" values under the integral sign. The fact that the other function is a Dirac delta function doesn't change this. So we obtain:
\begin{align*}
x(t) * h(t)
&=\delta(t-3-1)+\delta(t-3-3)-2\delta(t-4-1)-2\delta(t-7-3)\\
&=\delta(t-4)+\delta(t-6)-2\delta(t-5)-2\delta(t-10).
\end{align*}