There are various kinds of experience you can gain that would lead you more naturally to this solution. Roughly in order, they are as follows:
The first kind is having experience with the relationship between solving linear recurrences and taking powers of a matrix by diagonalizing it. The keyword here is companion matrix; the notion of the adjacency matrix of a graph is also relevant. The canonical example here is the relationship between the Fibonacci numbers and the Lucas numbers on the one hand and taking powers of the matrix $\left[ \begin{array}{cc} 1 & 1 \\\ 1 & 0 \end{array} \right]$ on the other, which is the adjacency matrix of a certain particularly simple graph. Slightly more generally, if $A$ is any matrix with integer coefficients, then the sequence
$$a_n = \text{tr}(A^n)$$
is a sequence of integers satisfying a linear recurrence whose coefficients are determined by the characteristic polynomial of $A$. But using linear algebra, we can be more specific: in fact $a_n$ must have the form
$$a_n = \sum_{i=1}^k \lambda_i^n$$
where $\lambda_i$ are the eigenvalues of the matrix $A$. In the Fibonacci example, $a_n$ is the Lucas numbers, and the eigenvalues are $\frac{1 \pm \sqrt{5}}{2}$; in particular, the Lucas numbers satisfy
$$L_n = \left( \frac{1 + \sqrt{5}}{2} \right)^n + \left( \frac{1 - \sqrt{5}}{2} \right)^n.$$
The second kind is having experience with some basic Galois theory and algebraic number theory. Here the key notion is understanding what a Galois conjugate of an algebraic number is, and in particular why the sum of all the Galois conjugates of an algebraic number is rational and why the sum of all the Galois conjugates of an algebraic integer is an integer. If you understand this you will immediately understand why, for example,
$$(a + \sqrt{b})^n + (a - \sqrt{b})^n$$
is always an integer. (There are more elementary ways to see this, but Galois theory / algebraic number theory is the appropriate context in which this observation should be placed.)
The third kind is having seen the notion of a Pisot-Vijarayaghavan number. These are real algebraic integers $\alpha > 1$ all of whose Galois conjugates have absolute value less than $1$. Consequently, the sequence
$$a_n = \sum_{i=1}^k \alpha_k^n$$
(where $\alpha_k$ runs over all of the Galois conjugates of $\alpha$) is dominated by $\alpha^n$, and in fact the remainder term vanishes so quickly that $a_n$ is necessarily the closest integer to $\alpha^n$.
The golden ratio $\frac{1 + \sqrt{5}}{2}$, as well as $\frac{7 + \sqrt{37}}{2}$, are both examples of Pisot-Vijarayaghavan numbers. It follows that the Lucas number $L_n$ is the closest integer to $\left( \frac{1 + \sqrt{5}}{2} \right)^n$, and if you have ever seen this result before or the corresponding result for the Fibonacci numbers then you should instantly think of using linear recurrences to solve this problem.