Define $c_0 = 2$, $c_1 = 4$, and $c_{n + 2} = 4c_{n + 1} + c_n$ for $n \geqslant 0$, then$$
c_n = (2 + \sqrt{5})^n + (2 - \sqrt{5})^n. \quad \forall n \in \mathbb{N}
$$
Note that$$
\begin{pmatrix}c_{n + 2}\\c_{n + 1}\end{pmatrix} = \begin{pmatrix}4 & 1\\ 1 & 0\end{pmatrix} \begin{pmatrix}c_{n + 1}\\c_n\end{pmatrix}.
$$
Denote $\displaystyle A = \begin{pmatrix}4 & 1 \\ 1 & 0\end{pmatrix}$, then$$
\begin{pmatrix}c_n\\c_{n - 1}\end{pmatrix} = A^{n - 1} \begin{pmatrix}c_1\\c_0\end{pmatrix}.
$$
(Tedious computation starts.) Because$$
2016 = 2^{10} + 2^9 + 2^8 + 2^7 + 2^6 + 2^5,
$$
and$$
A^2 = \begin{pmatrix}17 & 4 \\ 4 & 1\end{pmatrix}, A^{2^2} = \begin{pmatrix}305 & 72 \\ 72 & 17\end{pmatrix}, A^{2^3} = \begin{pmatrix}98209 & 23184 \\ 23184 & 5473\end{pmatrix}, \cdots
$$
then the explicit answer can be computed, but it is too hard to continue computing manually.