Su $\,a_i\,$ cumplen la misma ley de adición que los números de Fibonacci (con la misma facilidad de demostración que en esta respuesta mediante multiplicación matricial). Por lo tanto, la breve demostración que di en esta respuesta muestra que esta secuencia también es una secuencia de divisibilidad fuerte es decir $\, (a_m,a_n) = a_{\large (m,n)}\,$ que produce inmediatamente el resultado buscado.
Observación $\ $ En respuesta a los comentarios, a continuación se ofrecen más detalles. Como en el primer enlace tenemos
$$ \begin{bmatrix}a_2 &a_1\\ a_1 & a_0\end{bmatrix} = \begin{bmatrix}4 &1\\ 1 & 0\end{bmatrix},\quad \begin{bmatrix} a_{n+2} &\!\!\! a_{n+1}\\ a_{n+1} & \!\!\!a_n\end{bmatrix} = \begin{bmatrix} a_{n+1} &\!\! a_{n}\\ a_{n} & \!\!\!\!a_{n-1}\end{bmatrix} \begin{bmatrix}4 &1\\ 1 & 0\end{bmatrix}$$
Así pues, deducimos por inducción
$$ A_n := \begin{bmatrix} a_{n+1} &\!\! a_{n}\\ a_{n} &\!\!\! a_{n-1}\end{bmatrix} = \begin{bmatrix}4 &1\\ 1 & 0\end{bmatrix}^n\! =\, A_1^n $$
Por lo tanto, deducimos que $\,A_{m+n} = A_1^{m+n} = A_1^m A_1^n = A_m A_n,\ $ es decir
$$ \begin{align} \begin{bmatrix} a_{m+n+1} &\!\! a_{m+n}\\ a_{m+n} &\!\!\!\! a_{m+n-1}\end{bmatrix} &= \begin{bmatrix} a_{m+1} &\!\! a_{m}\\ \!a_{m} & \!\!\!\!a_{m-1}\end{bmatrix} \begin{bmatrix} a_{n+1} &\!\! a_{n}\\ a_{n} & \!\!\!\!a_{n-1}\end{bmatrix}\\[.5em] &= \begin{bmatrix}a_{m+1}a_{n+1}+a_m a_n &\! a_{m+1}a_n+a_m a_{n-1}\\ a_m a_{n+1}+a_{m-1}a_n &\! a_m a_n + a_{m-1} a_{n-1} \end{bmatrix}\\ \end{align} $$
Se obtiene la ley de adición $\ a_{m+n} =\, a_{m+1} a_n +a_m a_{n-1}.\ $ Por ejemplo
$$ \begin{align} \begin{bmatrix}a_8 & a_7\\ a_7 & a_6\end{bmatrix} &= \begin{bmatrix}a_4 & a_3\\ a_3 & a_2\end{bmatrix} \begin{bmatrix}a_5 & a_4\\ a_4 & a_3\end{bmatrix}\\[.4em] &= \begin{bmatrix}72 & \!17\\ 17 &\! 4\end{bmatrix} \begin{bmatrix}305 &\!\!\!\!72\\ 72 &\!\!\!\!\! 17\end{bmatrix} = \begin{bmatrix}23184 &\!\!\! 5473\\ 5473 &\!\!\! 1292\end{bmatrix} \end{align} $$
En relación con $\,f_n = \dfrac{x^n-y^n}{x-y},\,$ que satisface $\,f_{n+2} = (x\!+\!y) f_{n+1}-xy\, f_n,\,$ una prueba similar a la anterior demuestra que cumple la ley de adición $\, f_{m+n} = f_{m+1} f_n - xy\, f_m f_{n-1},\ $ es decir
$$ \dfrac{x^{m+n}\!-y^{m+n}}{x-y}\,=\, \dfrac{x^{m+1}\!-y^{m+1}}{x-y}\,\dfrac{x^{n}\!-y^{n}}{x-y} - xy\, \dfrac{x^{m}-y^{m}}{x-y}\,\dfrac{x^{n-1}-y^{n-1}}{x-y}$$
Para ayudar a disipar dudas en los comentarios, aquí%2F(x-y)%20%3D%20-xy(x%5E(m)-y%5E(m))%2F(x-y)(x%5E(n-1)-y%5E(n-1))%2F(x-y)%2B(x%5E(m%2B1)-y%5E(m%2B1))%2F(x-y)(x%5E(n)-y%5E(n))%2F(x-y)%3F) es una verificación Alfa de la ecuación anterior.