El grueso de este artículo es una colaboración enviada por fede a gaussianos (arroba) gmail (punto) com.

Introducción

Una fracción continua es una expresión del tipo:

x=a_0+\cfrac{1}{a_1+\cfrac{1}{a_2+\cfrac{1}{a_3+\cfrac{1}{\ddots}}}}

donde a_0 es un número entero y los demás a_i son enteros positivos.

La representación de un número real de este tipo en fracción continua tiene varias propiedades que hacen que dicha representación sea más interesante que la representación decimal habitual:

  • La representación en fracción continua de un número es finita si y solo si ese número es racional.
  • La representación en fracción continua de un racional simple es generalmente corta.
  • La representación en fracción continua de un racional es única siempre que no acabe en 1.
  • Los términos de una fracción continua se repetirán si y solo si representa a un irracional cuadrático, es decir, si es solución de una ecuación cuadrática con coeficientes enteros. Por ejemplo, la fracción continua [1, 1, 1, \ldots ] representa al número áureo y [1, 2, 2, 2, \ldots ] a \sqrt{2}.
  • El truncado de la representación en fracción continua de un número x da una aproximación racional que es, en cierto sentido, la mejor posible.

Todo número real puede representarse como fracción continua, pero en este artículo vamos a centrarnos en la representación continua de ciertos números racionales.

Fracción continua de un número racional


Si a_0, a_1, \ldots, a_n son los cocientes parciales que obtenemos aplicando el algoritmo de Euclides a una fracción irreducible \textstyle{\frac{p}{q}}, con p > q, todos los a_i serán enteros positivos.

Escribimos:

\cfrac{p}{q}= \cfrac{N[a_0, a_1, \ldots, a_n]}{D[a_0, a_1, \ldots, a_n]}=a_0+\cfrac{1}{a_1+\cfrac{1}{a_2+\cfrac{1}{\ddots+\cfrac{1}{a_n}}}} = [a_0, a_1, \ldots, a_n]

de forma que N[a_0,\ldots a_n] = p es el numerador de la fracción irreducible cuyo desarrollo en fracción continua tiene los cocientes parciales a_0,\ldots a_n y D[a_0,\ldots a_n] = q es el denominador de esa fracción.

En la teoría elemental de las fracciones continuas se demuestra fácilmente (por inducción) que

  • N[a_0]=a_0 y N[a_0,a_1]=a_0a_1+1
  • N[a_0,\ldots a_n] = a_nN[a_0,\ldots a_{n-1}]+N[a_0,\ldots a_{n-2}]
  • D[a_0]= 1 y D[a_0,a_1,\ldots a_n]=N[a_1,\ldots a_n]

Aunque en el resto del artículo no se va a utilizar esto es interesante resaltar que N[a_0,\ldots,a_n] también se puede expresar como un determinante:

N[a_0,\ldots,a_n]=\begin{vmatrix} a_0 & -1 & 0 & 0 & \ldots & 0 & 0 \\ 1 & a_1 & -1 & 0 & \ldots & 0 & 0 \\ 0 & 1 & a_2 & -1 & \ldots & 0 & 0 \\ & & \ldots & & \ldots \\ & & \ldots & & \ldots \\ 0 & 0 & 0 & 0 & \ldots & a_{n-1} & -1 \\ 0 & 0 & 0 & 0 & \ldots & 1 & a_n \end{vmatrix}

Interpretación combinatoria

En este artículo (pdf, 488kb) Benjamin, Quinn y Su dan la siguiente interpretación combinatoria de las cantidades N[a_0,\ldots a_n], que hace visualmente inmediatas algunas relaciones entre ellas.

Supongamos una regla de longitud n+1, dividida en n+1 casillas C_i \ (i=0,\ldots,n) de longitud 1.
Asignamos a cada casilla C_i un peso máximo a_i (que será siempre un entero positivo).

Formamos configuraciones con dos tipos de piezas, unas simples, de longitud 1, y otras dobles, de longitud 2, colocándolas sobre la regla de longitud n+1 de la siguiente forma:

  • Cubrimos toda la regla con una fila de piezas simples y dobles, colocadas arbitrariamente.
  • A continuación apilamos un número arbitrario (cero o más) de piezas simples siempre encima de piezas simples ya colocadas y de forma que el número total de piezas situadas sobre la casilla C_i no supere el peso máximo a_i de esa casilla.

La figura adjunta es una configuración que cumple la condición anterior con n+1=9 y con los pesos máximos a_i anotados bajo cada casilla: a_0=5, a_1=4, \ldots ,a_8=6.

Si C[a_0, \ldots,a_n] es el número de las distintas configuraciones de piezas que se pueden formar cumpliendo las condiciones anteriores, resulta que C[a_0, \ldots ,a_n] = N[a_0, \ldots ,a_n] , es decir, el número de configuraciones es el numerador de la fracción cuya fracción continua tiene los cocientes parciales a_i.

Usando por ejemplo esta calculadora obtenemos que el numero de posibles configuraciones distintas para los pesos máximos indicados en la figura es 225995.

Tenemos que C[a_0] = a_0 y contando las configuraciones que no tienen y que tienen una pieza doble cubriendo las 2 últimas casillas se ve enseguida que C[a_0,a_1] = a_0a_1 +1 , y C[a_0, \ldots,a_n] = a_nC[a_0, \ldots,a_{n-1}]+C[a_0, \ldots,a_{n-2}].

Como C[a_0,\ldots,a_n] y N[a_0,\ldots,a_n] cumplen la misma relación de recurrencia con los mismos valores iniciales resulta que C[a_0,\ldots,a_n] = N[a_0,\ldots,a_n].

La regla de Euler

La anterior interpretación combinatoria no es más que una reformulación de la regla de Euler que dice que N[a_0,\ldots,a_n] es la suma de determinados productos de los cocientes parciales a_i, que son:

  • El producto de todos los cocientes parciales.
  • Los productos obtenidos suprimiendo de todas la formas posibles 1 pareja de cocientes parciales de indices consecutivos.
  • Los productos obtenidos suprimiendo de todas la formas posibles 2 parejas (disjuntas) de cocientes parciales de indices consecutivos.
  • Etc.

Por otro lado la regla de Euler se obtiene inmediatamente de la interpretación combinatoria.

Consecuencias

La interpretación combinatoria permite demostrar sin usar inducción, estableciendo biyecciones entre conjuntos, identidades relativas a las cantidades N[a_0,\ldots,a_n] .

Por ejemplo, resulta visualmente claro que N[a_0,\ldots ,a_n] = N[a_n,\ldots,a_0], pues hay tantas posibles configuraciones si los mismos pesos están situados de izquierda a derecha que de derecha a izquierda.

También es inmediato, considerando las configuraciones que tienen o no una pieza doble cubriendo las casillas de índices k y k+1, que:

N[a_0,\ldots,a_k,\ldots,a_n] = N[a_0,\ldots,a_k] N[a_{k+1},\ldots,a_n]+N[a_0,\ldots,a_{k-1}] N[a_{k+2},\ldots,a_n]

En el artículo referenciado anteriormente se dan demostraciones biyectivas de otras identidades, por ejemplo la que da la diferencia entre dos convergentes consecutivos, que en la notación que hemos usado, es equivalente a:

N[a_0 \ldots,a_n] N[a_1, \ldots,a_{n-1}] - N[a_0 \ldots,a_{n-1}] N[a_1, \ldots,a_n]=(-1)^{n-1}
Print Friendly, PDF & Email
5 1 vote
Article Rating

¿Te ha gustado la entrada? Puedes invitarme a un café, Gauss te lo agradecerá 😉