Textbook answer: $a_{k}=-\frac{1}{2} + \frac{3}{2}3^{k}$
Here is what I did so far:
$G(x) = \frac{1}{(1-x)(1-3x)} = \left ( \sum_{k=0}^{\infty}x^{k} \right )\left ( \sum_{k=0}^{\infty}(3x)^{k} \right )$
$= 3^{k}\left ( \sum_{k=0}^{\infty}x^{k} \right )\left ( \sum_{k=0}^{\infty}x^{k} \right )=3^{k}\left ( \sum_{k=0}^{\infty}x^{k} \right )^{2}$
$= 3^{k}(1+x^{2}+x^{3}+\cdots)$
$= 3^{k}(1+x+x^{2}+x^{3}+\cdots)-x$ $= 3^{k}\frac{1}{1-x}-x$
But how did they get the $x=-\frac{1}{2}$ and $\frac{1}{1-x}=\frac{3}{2}$?
Many of the steps above are careless errors and irrecoverably wrong. Here is how to do it:
You have: $$ G(x)= \left(\sum x^k\right)\left(\sum 3^k x^k\right) $$
So you realise this is a convolution and work out the coefficient of $x^k$:
$$ G(x) = \sum_{k=0} \left(\sum_{\ell=0}^k x^{k-\ell}3^\ell x^\ell\right) $$
$$ G(x) = \sum_{k=0} x^k \left(\sum_{\ell=0}^k 3^\ell \right) $$
Now we can simplify the inner sum as we know $\sum_{\ell=0}^k 3^\ell = \frac{3^{k+1}}{2}-\frac 12$ and the result follows