Given $x_1$ and $x_2$ are roots of $P(x)=x^2-6x+1$ prove using Binomial expansion that $x_1^n+x_2^n$ is not divisible by 5 if $n\in \Bbb N\setminus\{0\}$.
Source: list of problems for the preparation for math contests.
Notice: this problem is different from other SE question, as it requires a specific proof that was suggested by @DeepSea in a comment in that post, but not developed.
My attempt: By solving the roots of $P(x)$ it is easy to see that $$x_1=3+2\sqrt{2}\ \ \text{and}\ \ \ x_2=3-2\sqrt{2}$$ and, using the Binomial expansion, $$x_1^n+x_2^n=(x_1+x_2)^n-\sum_{k=1}^{n-1}{{n-1}\choose {k}} x_1^{n-k}x_2^k=6^n-\sum_{k=1}^{n-1} {{n-1}\choose {k}}x_1^{n-k}x_2^k$$
It is easy to see that $6^n\equiv 1\pmod{5}$ for all possible $n$. Therefore the problem is showing that for $$A(n)=\sum_{k=1}^{n-1}{{n-1}\choose {k}}x_1^{n-k}x_2^k,\ \ A(n)\not\equiv 1\pmod{5},\ \forall n\in \Bbb N\setminus\{0\}.$$
My problem is on how to proceed with the proof of this last statement. Hints and answers are welcomed.
I think you misinterpreted the suggestion by @DeepSea. He probably meant expanding $(3\pm2\sqrt{2})^n$ instead of $(x_1+x_2)^n$.
So if we follow this and binomial expand $(3\pm 2\sqrt{2})^n$, we have $$ x_1^n+x_2^n=\sum_{k=0}^n\binom{n}{k}3^{n-k} (2\sqrt{2})^k \underbrace{[1+(-1)^k]}_{=\begin{cases}2&k\text{ even}\\0&k\text{ odd}\end{cases}} $$ So we cancel the odd $k$ terms and are left with $$ x_1^n+x_2^n=2\sum_{k=0}^{\lfloor n/2\rfloor}\binom{n}{2k}3^{n-2k} (2\sqrt{2})^{2k} \equiv 2\cdot3^n\sum_{k=0}^{\lfloor n/2\rfloor}\binom{n}{2k}3^{-k}\pmod{5}. $$ Hence we want to prove $$ \sum_{k=0}^{\lfloor n/2\rfloor}\binom{n}{2k}(-3)^k\not\equiv 0\pmod{5} $$ which I don't see a way without
Cheating
So we look at an equivalent problem with roots $1\pm\sqrt{-3}\pmod{5}$, which if we take out a common factor of $2$ (doesn't change mod 5 nonzero) becomes a problem with roots $y_\pm=\frac12(1\pm\sqrt{-3})$ (which are of course the roots of $y^2-y+1$, what you get from reducing coefficients of $P$ mod $5$). Now these roots are sixth roots of unity, so $y_+^n+y_-^n$ only take values $\pm 1,\pm 2$.
This leaves the question: why not reduce $P$ mod 5 right at the start and be done with it?