I already know that $b_{n+1} = a_n +3b_n$ and $a_{n+1} = 3a_n - b_n$. So
$a_{n+2} = 3(3a_n-b_n)-(3b_n+a_n) = 9a_n-3b_n-3b_n-a_n = 8a_n-6b_n$ and
$b_{n+2} = 8b_n+6a_n$.
So we can rewrite the whole thing as
$8a_n-6b_n+r(3a_n-b_n)+sa_n = 8a_n-6b_n+r(3b_n+a_n)+sb_n$ which, in turn is:
$(-4r+s)a_n = (4r+s)b_n$.
The original problem states that $(3+i)^n = a_n+ib_n$ so I tried using n=1 so a = 3 and b = 1
so then I have $-12r+3s=4r+s$ so
$2s=16r$.
From here I'm stuck. Do I just try random options? Honestly, I don't feel like I've done it right so far.