Find expressions for $a_{n+1}$ and $b_{n+1}$ as linear combinations of $a_n$ and $b_n$ with coefficients independent of n. With some of your comments, I see $a_{n+1} +ib_{n+1} = (a_n+ib_n)(3+i) = 3a_n + ia_n + i3b_n-b_n$. So the imaginary parts have to be equal which means that $b_{n+1} = a_n +3b_n$ and the real parts have to be equal so $a_{n+1} = 3a_n - b_n$, right? So that's that question I believe
Show that for $n\geq 1, a_n \equiv3\pmod 5$ and $b_n \equiv 1\pmod 5$ Here we know that if n = 1, $(3+i) = a_1 + ib_1$ so $a_1 = 3$ and $b_1 = 1$ which means that $a_1 \equiv 3(mod 5)$ and $b_1 \equiv 1 (mod 5)$. We can use these as our base case for $a_n$ and $b_n$ and see that $a_{n+1} = 3*3-1 = 8 \equiv 3(mod 5)$ and $b_{n+1} = 3*1+3 = 6 \equiv 1(mod 5)$. [Thank you to J.W. Tanner for your help with this]
Hint: $a_{n+1}+i b_{n+1} = (a_n+ib_n)(3+i)$.