The property of the solutions of Pell-Fermat equation Write
$$ 
(1+\sqrt{2})^n=a_n+b_n\sqrt{2}
$$
for $n,a_n,b_n$ are positive integers.
For each prime number $p$,
I would like to ask  if $p|b_n$ for some $n$？
Thanks.

Yes - The $a_n$ and $b_n$ satisfy the recurrences $a_{n+1}=a_n+2b_n$ and $b_{n+1} = a_n+b_n$, for $n\in\mathbb{N}$. This can be used to extend the definition of the $a_n$ and $b_n$ to all $n$, including zero and negative. Now consider these sequences $\pmod{p}$. Since there are only $p^2$ possible (ordered) pairs of remainders $(a_n \pmod{p},b_n \pmod{p})$, the sequence of pairs must repeat in a loop of at most $p^2$. Since, with the extended definition, we have that $b_0 \equiv 0 \pmod{p}$, we will eventually have $b_k \equiv 0 \pmod{p}$ for some $k$ with $0 < k \leq p^2$, and then for infinitely many $b_n$.
