Teoría de NúmerosContenidos

Sucesiones recurrentes lineales: Fibonacci, Lucas y más

Una sucesión recurrente lineal de orden satisface . Su comportamiento modular es periódico, su forma cerrada explícita vía raíces del polinomio característico, y aparece en divisibilidad, combinatoria y geometría.

DificultadRegional
Etiquetasrecurrenciasfibonaccilucasbinetpisanodiofanticas
Requisitospolinomioscongruencias
AutorAdrián García Bouzas
Actualizado2026-06-03

La sucesión de Fibonacci (0,1,1,2,3,5,8,13,21,0, 1, 1, 2, 3, 5, 8, 13, 21, \ldots) es probablemente la sucesión más famosa de las matemáticas. Aparece en la naturaleza (espirales de girasoles, disposición de hojas), en la geometría (razón áurea), en la combinatoria (número de formas de cubrir un tablero con dominós) y, por supuesto, en la teoría de números.

Pero Fibonacci no es especial por sí mismo: es el representante más simple de una familia mucho más amplia — las recurrencias lineales homogéneas — cuya teoría es elegante y completa. Toda sucesión de esta familia tiene una forma cerrada explícita, un comportamiento modular periódico, y satisface una familia de identidades algebraicas que la hacen extraordinariamente útil en problemas olímpicos.

Definición

Una recurrencia lineal homogénea de orden kk es una sucesión {an}n0\{a_n\}_{n \geq 0} que satisface, para ciertos coeficientes c1,,ckc_1, \ldots, c_k con ck0c_k \neq 0:

an+k=c1an+k1+c2an+k2++ckan,n0.a_{n+k} = c_1 a_{n+k-1} + c_2 a_{n+k-2} + \cdots + c_k a_n, \qquad n \geq 0.

La sucesión queda completamente determinada por las condiciones iniciales a0,a1,,ak1a_0, a_1, \ldots, a_{k-1}.

Las sucesiones clásicas

Fibonacci. F0=0F_0 = 0, F1=1F_1 = 1, Fn+2=Fn+1+FnF_{n+2} = F_{n+1} + F_n. 0,1,1,2,3,5,8,13,21,34,55,89,144,0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, \ldots

Lucas. L0=2L_0 = 2, L1=1L_1 = 1, Ln+2=Ln+1+LnL_{n+2} = L_{n+1} + L_n. 2,1,3,4,7,11,18,29,47,76,2, 1, 3, 4, 7, 11, 18, 29, 47, 76, \ldots

Ambas tienen la misma recurrencia pero distintas condiciones iniciales. Son las dos sucesiones de Fibonacci «fundamentales».

Pell. P0=0P_0 = 0, P1=1P_1 = 1, Pn+2=2Pn+1+PnP_{n+2} = 2P_{n+1} + P_n. Aparece en la ecuación de Pell x22y2=±1x^2 - 2y^2 = \pm 1: las soluciones son (Pn+1,Pn)(P_{n+1}, P_n).

Sucesión de Lucas generalizada. Para parámetros P,QP, Q: U0=0U_0 = 0, U1=1U_1 = 1, Un+2=PUn+1QUnU_{n+2} = P U_{n+1} - Q U_n. Fibonacci es P=1P = 1, Q=1Q = -1; Pell es P=2P = 2, Q=1Q = -1.

Teorema

(Solución general via polinomio característico) El polinomio característico de la recurrencia an+k=c1an+k1++ckana_{n+k} = c_1 a_{n+k-1} + \cdots + c_k a_n es

χ(x)=xkc1xk1c2xk2ck.\chi(x) = x^k - c_1 x^{k-1} - c_2 x^{k-2} - \cdots - c_k.

Si χ\chi tiene kk raíces distintas r1,r2,,rkr_1, r_2, \ldots, r_k (complejas), la solución general es:

an=A1r1n+A2r2n++Akrkn,a_n = A_1 r_1^n + A_2 r_2^n + \cdots + A_k r_k^n,

donde A1,,AkA_1, \ldots, A_k se determinan de las condiciones iniciales. Si una raíz rjr_j tiene multiplicidad mm, su contribución es (B0+B1n++Bm1nm1)rjn(B_0 + B_1 n + \cdots + B_{m-1} n^{m-1}) r_j^n.

Demostración

(Idea) El conjunto de soluciones de la recurrencia forma un espacio vectorial de dimensión kk (pues los datos iniciales a0,,ak1a_0, \ldots, a_{k-1} determinan la sucesión). Las funciones nrjnn \mapsto r_j^n con rjr_j raíz del polinomio característico son soluciones independientes (verificable por sustitución), y si las kk raíces son distintas, son linealmente independientes (por la no-singularidad de la matriz de Vandermonde de los rjr_j). Por tanto forman una base del espacio. \blacksquare

Ejemplo

La fórmula de Binet

Ejemplo 1. Derivar la fórmula de Binet para Fibonacci.

El polinomio característico de Fn+2=Fn+1+FnF_{n+2} = F_{n+1} + F_n es x2x1x^2 - x - 1. Las raíces son:

φ=1+521.618,φ^=1520.618.\varphi = \frac{1 + \sqrt 5}{2} \approx 1.618, \qquad \hat\varphi = \frac{1 - \sqrt 5}{2} \approx -0.618.

La solución general es Fn=Aφn+Bφ^nF_n = A\varphi^n + B\hat\varphi^n. De las condiciones iniciales:

  • F0=0F_0 = 0: A+B=0A + B = 0, así B=AB = -A.
  • F1=1F_1 = 1: AφAφ^=A(φφ^)=A5=1A\varphi - A\hat\varphi = A(\varphi - \hat\varphi) = A\sqrt{5} = 1, así A=1/5A = 1/\sqrt{5}.

Fn=φnφ^n5.\boxed{F_n = \frac{\varphi^n - \hat\varphi^n}{\sqrt 5}.}

Consecuencias:

  • FnF_n es el entero más cercano a φn/5\varphi^n / \sqrt 5 (pues φ^<1|\hat\varphi| < 1, así φ^n/5<1/2|\hat\varphi^n/\sqrt 5| < 1/2).
  • Fnφn/5F_n \sim \varphi^n / \sqrt 5 exponencialmente.

Ejemplo 2. Los números de Lucas y su fórmula.

Para LnL_n con la misma recurrencia: Ln=φn+φ^nL_n = \varphi^n + \hat\varphi^n.

Verificación: L0=1+1=2L_0 = 1 + 1 = 2 ✓, L1=φ+φ^=1L_1 = \varphi + \hat\varphi = 1 ✓.

La relación entre Fibonacci y Lucas: Ln=Fn1+Fn+1L_n = F_{n-1} + F_{n+1} y Fn=(Ln1+Ln+1)/5F_n = (L_{n-1} + L_{n+1})/5.


Identidades clásicas

Ejemplo 3. Demostrar la identidad de Cassini: Fn1Fn+1Fn2=(1)nF_{n-1} F_{n+1} - F_n^2 = (-1)^n.

Por Binet:

Fn1Fn+1Fn2=(φn1φ^n1)(φn+1φ^n+1)(φnφ^n)25.F_{n-1} F_{n+1} - F_n^2 = \frac{(\varphi^{n-1} - \hat\varphi^{n-1})(\varphi^{n+1} - \hat\varphi^{n+1}) - (\varphi^n - \hat\varphi^n)^2}{5}.

Expandiendo el numerador:

φ2nφn1φ^n+1φ^n1φn+1+φ^2nφ2n+2(φφ^)nφ^2n\varphi^{2n} - \varphi^{n-1}\hat\varphi^{n+1} - \hat\varphi^{n-1}\varphi^{n+1} + \hat\varphi^{2n} - \varphi^{2n} + 2(\varphi\hat\varphi)^n - \hat\varphi^{2n}

=φn1φ^n+1φ^n1φn+1+2(φφ^)n= -\varphi^{n-1}\hat\varphi^{n+1} - \hat\varphi^{n-1}\varphi^{n+1} + 2(\varphi\hat\varphi)^n

=(φφ^)n1(φ^2φ2+2φφ^)= (\varphi\hat\varphi)^{n-1}(-\hat\varphi^2 - \varphi^2 + 2\varphi\hat\varphi)

=(φφ^)n1((φφ^)2)=(1)n1(5)=5(1)n.= (\varphi\hat\varphi)^{n-1} \cdot (-(\varphi - \hat\varphi)^2) = (-1)^{n-1} \cdot (-5) = 5 (-1)^n.

Dividiendo por 55: Fn1Fn+1Fn2=(1)nF_{n-1} F_{n+1} - F_n^2 = (-1)^n. \square

(Alternativa: inducción o matrices 2×22 \times 2, ambas más elegantes.)


Ejemplo 4. Demostrar gcd(Fm,Fn)=Fgcd(m,n)\gcd(F_m, F_n) = F_{\gcd(m, n)}.

La clave es la identidad de adición: Fm+n=FmFn+1+Fm1FnF_{m+n} = F_m F_{n+1} + F_{m-1} F_n.

Esto implica que FmFknF_m \mid F_{kn} para todo kk entero (por inducción). En particular, Fgcd(m,n)FmF_{\gcd(m,n)} \mid F_m y Fgcd(m,n)FnF_{\gcd(m,n)} \mid F_n.

Por otro lado, si dFmd \mid F_m y dFnd \mid F_n, de la identidad de adición dFmnd \mid F_{m-n} (por linealidad). El argumento del algoritmo de Euclides aplicado a (Fm,Fn)(F_m, F_n) con la operación FmFmnF_m \mapsto F_{m-n} reproduce el algoritmo de Euclides en los índices, dando gcd(Fm,Fn)=Fgcd(m,n)\gcd(F_m, F_n) = F_{\gcd(m, n)}. \square

Corolario: FmFn    mnF_m \mid F_n \iff m \mid n.


Periodicidad modular y período de Pisano

Ejemplo 5. Demostrar que {Fnmodm}\{F_n \bmod m\} es periódica para todo m1m \geq 1.

Consideramos los pares (Fnmodm,Fn+1modm)(F_n \bmod m, F_{n+1} \bmod m). Hay a lo sumo m2m^2 pares posibles. Por el principio del palomar, dos índices i<ji < j dan el mismo par: (Fi,Fi+1)(Fj,Fj+1)(modm)(F_i, F_{i+1}) \equiv (F_j, F_{j+1}) \pmod m.

Como la recurrencia es invertible (también vale hacia atrás: Fn=Fn+2Fn+1F_n = F_{n+2} - F_{n+1}), la sucesión es periódica en ambas direcciones, y el período es jij - i. \square

El período de {Fnmodm}\{F_n \bmod m\} se llama período de Pisano y se denota π(m)\pi(m).

Valores conocidos: π(2)=3\pi(2) = 3, π(3)=8\pi(3) = 8, π(4)=6\pi(4) = 6, π(5)=20\pi(5) = 20, π(7)=16\pi(7) = 16, π(8)=12\pi(8) = 12, π(10)=60\pi(10) = 60.


Ejemplo 6. Calcular F10100mod7F_{10^{100}} \bmod 7.

π(7)=16\pi(7) = 16. Necesitamos 10100mod1610^{100} \bmod 16.

1010(mod16)10 \equiv 10 \pmod{16}, 102=1004(mod16)10^2 = 100 \equiv 4 \pmod{16}, 10416010^4 \equiv 16 \equiv 0... wait: 102=100=616+410^2 = 100 = 6 \cdot 16 + 4, así 102410^2 \equiv 4. 10416010^4 \equiv 16 \equiv 0? No: 42=160(mod16)4^2 = 16 \equiv 0 \pmod{16}. Entonces 10100=(102)50450=(42)2516250(mod16)10^{100} = (10^2)^{50} \equiv 4^{50} = (4^2)^{25} \equiv 16^{25} \equiv 0 \pmod{16}.

Entonces F10100F00(mod7)F_{10^{100}} \equiv F_0 \equiv 0 \pmod 7.

Verificación: 161010016 \mid 10^{100} pues 102=10010^2 = 100 y 10100=(100)5010^{100} = (100)^{50}, y 1004(mod16)100 \equiv 4 \pmod{16}, 42=1604^2 = 16 \equiv 0, 450=4225=(42)250(mod16)4^{50} = 4^{2 \cdot 25} = (4^2)^{25} \equiv 0 \pmod{16}. ✓


Ejemplo 7. Demostrar que FnF_n es divisible por 55 si y solo si n0(mod5)n \equiv 0 \pmod 5.

π(5)=20\pi(5) = 20. Calculamos los primeros 2020 valores módulo 55:

F0,F1,,F190,1,1,2,3,0,3,3,1,4,0,4,4,3,2,0,2,2,4,1(mod5).F_0, F_1, \ldots, F_{19} \equiv 0, 1, 1, 2, 3, 0, 3, 3, 1, 4, 0, 4, 4, 3, 2, 0, 2, 2, 4, 1 \pmod 5.

Los ceros aparecen en posiciones 0,5,10,150, 5, 10, 15, exactamente los múltiplos de 55. Como el período es 2020 y 5205 \mid 20, el patrón se repite. \square

(Esto también se sigue de gcd(Fm,Fn)=Fgcd(m,n)\gcd(F_m, F_n) = F_{\gcd(m,n)}: 5Fn    F5Fn    5n5 \mid F_n \iff F_5 \mid F_n \iff 5 \mid n.)


Recurrencias en problemas olímpicos

Ejemplo 8. Demostrar que Fn2+Fn+12=F2n+1F_n^2 + F_{n+1}^2 = F_{2n+1}.

Probaremos por inducción. Caso base n=0n = 0: 02+12=1=F10^2 + 1^2 = 1 = F_1 ✓.

Paso inductivo: suponemos Fn2+Fn+12=F2n+1F_n^2 + F_{n+1}^2 = F_{2n+1} y probamos para n+1n + 1:

Fn+12+Fn+22=Fn+12+(Fn+1+Fn)2=2Fn+12+2FnFn+1+Fn2F_{n+1}^2 + F_{n+2}^2 = F_{n+1}^2 + (F_{n+1} + F_n)^2 = 2F_{n+1}^2 + 2F_n F_{n+1} + F_n^2

=Fn+1(2Fn+1+2Fn)+Fn2Fn+12+2Fn+12= F_{n+1}(2F_{n+1} + 2F_n) + F_n^2 - F_{n+1}^2 + 2F_{n+1}^2

Mejor usar la identidad de adición: Fm+n=FmFn+1+Fm1FnF_{m+n} = F_m F_{n+1} + F_{m-1} F_n.

Con m=n+1m = n + 1, n=n+1n = n + 1: F2n+2=Fn+12+FnFn+1F_{2n+2} = F_{n+1}^2 + F_n F_{n+1} y F2n+3=Fn+2Fn+1+Fn+1Fn+...F_{2n+3} = F_{n+2}F_{n+1} + F_{n+1}F_n + ...

Más directo: la identidad F2n+1=Fn+12+Fn2F_{2n+1} = F_{n+1}^2 + F_n^2 se puede probar con matrices:

(Fn+1FnFnFn1)=(1110)n.\begin{pmatrix} F_{n+1} & F_n \\ F_n & F_{n-1} \end{pmatrix} = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}^n.

Elevando al cuadrado (potencia 2n2n) y comparando entradas: F2n+1=Fn+12+Fn2F_{2n+1} = F_{n+1}^2 + F_n^2. \square

Identidades de adición y de duplicación

Recopilamos las identidades más útiles:

  • Adición: Fm+n=FmFn+1+Fm1FnF_{m+n} = F_m F_{n+1} + F_{m-1} F_n.
  • Duplicación: F2n=Fn(2Fn+1Fn)=FnLnF_{2n} = F_n(2F_{n+1} - F_n) = F_n L_n.
  • Cassini: Fn1Fn+1Fn2=(1)nF_{n-1}F_{n+1} - F_n^2 = (-1)^n.
  • Catalan generalizada: Fn+rFnrFn2=(1)nrFr2F_{n+r}F_{n-r} - F_n^2 = (-1)^{n-r} F_r^2.
  • Suma: k=1nFk=Fn+21\sum_{k=1}^n F_k = F_{n+2} - 1.
  • Suma de cuadrados: k=1nFk2=FnFn+1\sum_{k=1}^n F_k^2 = F_n F_{n+1}.
  • Zeckendorf: todo entero positivo se escribe de forma única como suma de Fibonacci no consecutivos.
Aplicaciones

Teorema de Zeckendorf. Todo entero positivo nn se escribe de manera única como n=Fi1+Fi2++Firn = F_{i_1} + F_{i_2} + \cdots + F_{i_r} con ijij+12i_j - i_{j+1} \geq 2 (Fibonacci no consecutivos) y ir2i_r \geq 2.

Dominós y combinatoria. El número de formas de cubrir un tablero de 1×n1 \times n con piezas de 1×11 \times 1 y 1×21 \times 2 es Fn+1F_{n+1} (recurrencia inmediata: la primera pieza ocupa 11 o 22 casillas).

Ecuación de Pell. La sucesión de Pell PnP_n da todas las soluciones (Pn+1,Pn)(P_{n+1}, P_n) de x22y2=±1x^2 - 2y^2 = \pm 1. La recurrencia garantiza que todas las soluciones se obtienen así.

Observación

Por qué las recurrencias lineales son ubicuas. Una sucesión recurrente lineal es esencialmente la misma cosa que un sistema lineal discreto. La «forma cerrada» via raíces del polinomio característico es el análogo discreto de la solución de una EDO lineal con coeficientes constantes. Esta analogía permite transferir toda la teoría de EDOs al mundo discreto: estabilidad, resonancia, comportamiento asintótico.

Módulo primo. Para pp primo, la sucesión {Fnmodp}\{F_n \bmod p\} tiene período π(p)\pi(p) que divide a 2(p1)2(p - 1) o 2(p+1)2(p + 1) según p±1(mod5)p \equiv \pm 1 \pmod 5 o p±2(mod5)p \equiv \pm 2 \pmod 5. Este es un resultado profundo conectado con la teoría de cuerpos de clase.