Triángulo trinomial

Al igual que construíamos el triángulo de Pascal para binomiales sumando cada dos consecutivos, podemos definir el triángulo trinomial si cada coeficiente es la suma de los tres situados en la fila anterior y sobre él. Los coeficientes de los extremos serán suma de dos o de uno, como si a izquierda y derecha del triángulo existieran ceros. Así quedaría con una hoja de cálculo.

 

A diferencia de los binomiales, estos coeficientes del triángulo forman filas de longitud impar, 2n+1, por lo que se puede hablar siempre de términos centrales, que en la imagen son 1, 1, 3, 7 y 19.

Este proceso equivale a la siguiente recurrencia:

Como curiosidad, John D. Cook relaciona estos números con los pasos del rey en el ajedrez para llegar a las distintas posiciones

 (ver https://www.johndcook.com/blog/2025/05/16/trinomial-coefficients-and-kings/)

Este triángulo posee diversas propiedades, tal como ocurría con el binomial. Por ejemplo, la primera diagonal de la izquierda contiene los números naturales, 1, 2, 3, 4, … y la segunda los triangulares, 1, 3, 6, 10, 15, …Igualmente, si en los binomiales la suma por filas era 2ⁿ, aquí es 3ⁿ. Por ejemplo, 1+3+6+7+6+3+1=27=3³

Estos números se pueden interpretar también como los coeficientes de la potencia (1+x+x²)ⁿ

Por ejemplo, para n=4 nos queda, en wxMaxima:

Los coeficientes coinciden con los de la fila 4 del triángulo.

Para distinguir estos coeficientes de los binomiales, se suele añadir un 2 a la derecha de su símbolo. El resto queda igual porque vemos que dependen de solo dos parámetros. Así que, a partir de ahora, los representaremos así:

En ese símbolo k es un entero que cumple -n ≤ k ≤ n, por lo que puede tomar valores negativos. Lo podemos interpretar como la diferencia entre el lugar central de la fila y el del término definido.

Cálculo directo

Según la Wikipedia, existe una fórmula para el cálculo directo:

Usa un sumatorio algo complicado, por lo que es más directo y formativo el seguir la recurrencia definida.

Por recurrencia

El esquema de recurrencia ya presentado se puede traducir a una función directa, en la que dados n y k nos devuelva el coeficiente correspondiente en el triángulo trinomial. Para ello, traducimos las filas de hoja de cálculo a dos vectores u(i) y v(i) que comiencen y terminen con ceros, y con ellos ir construyendo la recurrencia hasta llegar al coeficiente deseado:

Function trinomial2(n, k)

Dim u(50), v(50), i, j, h, t

 

 

For i = 1 To 2 * n + 3: u(i) = 0: v(i) = 0:  Next I ‘Se rellenan los vectores con ceros

u(n + 3) = 1 ‘Primer coeficiente

For i = 1 To n ‘Se trabaja hasta la fila n

For j = 3 To 2 * n + 3 ‘Longitud de la fila contando con ceros iniciales

v(j) = u(j - 1) + u(j) + u(j + 1) ‘Recurrencia

Next j

For h = 3 To 2 * n + 3: u(h) = v(h: Next h ‘Se copia v(i) en u(i)

For h = 1 To 2 * n + 1: v(h) = 0: Next h ‘Y se rellena con ceros

Next i

t = u(n + k + 3) ‘Se localiza el coeficiente pedido

trinomial2 = t

End Function

 

Aunque parece algo complejo, funciona con rapidez. En la imagen se recogen los coeficientes correspondientes a n=5 y n=6:

Coeficientes centrales

Los números centrales de cada fila poseen importancia propia. Por eso, en OEIS, se les dedican muchas referencias y descripciones de propiedades. Aquí sólo se darán algunos detalles, remitiendo para un estudio completo a la página https://oeis.org/A002426.

Su cálculo no es difícil. En la página mencionada se usa la definición como coeficientes de  (1+x+x²)ⁿ, pero eligiendo el coeficiente número n. Se propone este código PARI de fácil comprensión, extendido aquí con un bucle FOR:

a(n) = if( n<0, 0, polcoeff( (1 + x + x^2)^n, n))

for(k=0,20,print1(a(k),", "))

Su resultado, en la página oficial de PARI es:

Coincide con lo publicado.

Con nuestra función trinomial2(n,k) basta dar a k el valor 0. Así se ha efectuado para crear esta columna de Excel:

Una interpretación interesante de estos números es la de David Callan, que los identifica con el número de sucesiones de n elementos tomados del conjunto (1, 2, 3, … n), crecientes en sentido amplio, en las que ningún elemento se repite más de dos veces. Propone como ejemplo si n = 3, las sucesiones son 112, 113, 122, 123, 133, 223, 233, es decir, son siete como indica su coeficiente central.

He añadido la condición FMAX a mi herramienta Cartesius (https://www.hojamat.es/sindecimales/combinatoria/herramientas/herrcomb.htm#cartesius) y, aunque aún no está publicada, la he usado para comprobar esta propiedad. He concretado las condiciones siguientes:

xtotal=5

xt=1..5

FMAX<3

Creciente

Con ellas se combinan en arreglos crecientes los números del 1 al 5, con la condición de que la frecuencia máxima sea 2, tal como plantea David Callan, y se han obtenido 51 posibilidades, como corresponde al coeficiente central de orden 5. Este recorte contiene los primeros y su total:

 

 

 

Coeficientes trinomiales

Los coeficientes trinomiales son una extensión de los binomiales, definidos estos como

Y que son los coeficientes del binomio de Newton:

De forma análoga se pueden definir los coeficientes trinomiales, propios de un trinomio similar al anterior.

En realidad, los binomiales se podrían expresar de otra forma:

A partir de aquí podemos explicar mejor la definición de coeficientes trinomiales:

También, al igual que con los binomiales, los podemos identificar con los coeficientes de la potencia de un trinomio:

Las demostraciones son similares a las del caso binomial.

Esta fórmula produce siempre resultados enteros positivos, debido al teorema de que un producto de k números consecutivos decrecientes es siempre múltiplo de los k primeros números enteros positivos. Es fácil recordar que coincide con el número de permutaciones con repetición en las que se fija el número de cada una.

Para entender esto mejor podemos descomponer la definición en tres bloques. Lo vemos con un ejemplo:

A la vista de este esquema de tres cocientes, observamos que el tercer número de la definición no influye en el resultado, luego podemos usar sólo n₁ y n₂. Incluso estos son intercambiables, los que nos lleva a que el coeficiente trinomial se puede expresar como

En nuestro ejemplo, usando la función COMBINAT de Excel y Calc, nos quedaría:

COMBINAT(8;2)*COMBINAT(6;3)=560

COMBINAT(8;3)*COMBINAT(5;2)=560

 

Cálculo de un coeficiente trinomial en la práctica

Si no se desea depender de la función COMBINAT, se puede calcular directamente. El problema con el cálculo práctico de un coeficiente de este tipo es que los factoriales alcanzan números muy grandes. Si usamos calculadora o bien hoja de cálculo, los resultados pasan a coma flotante y se pierde el carácter entero de estos coeficientes. Por ello, es conveniente usar el desarrollo del párrafo anterior: dividir los números consecutivos a partir de m entre el primer factorial, los siguientes entre el segundo, y parar el cálculo, porque el tercer factorial produce un uno. En VBasic y otros lenguajes se puede construir una función similar a la siguiente:

Function trinomial(m, x, y, z)

Dim t, k, v

If x + y + z <> m Then trinomial = 0: Exit Function’ Datos no válidos

t = m: k = x: v = 1

While k > 0

v = v * t / k: t = t - 1: k = k – 1 ‘Divide factores consecutivos entre el primer coeficiente. Es el primer combinatorio de la fórmula.

Wend

k = y

While k > 0

v = v * t / k: t = t - 1: k = k – 1’ Igual, en el segundo coeficiente para calcular el segundo combinatorio.

Wend

trinomial = v

End Function

Presenta un pequeño riesgo de redondeo inadecuado, pero suele funcionar bastante bien. Devuelve un cero si los coeficientes de abajo no suman el de arriba. No analiza si los datos son enteros.

Es muy sencillo el traslado a cualquier otro lenguaje de programación, por ejemplo, en PARI, con un rango de datos y resultados mucho mayor:

trinomial(m, x, y, z)=my(t=m,k=x,v=1);while(k>0,v=v*t/k;t-=1;k-=1);k=y;while(k>0,v=v*t/k;t-=1;k-=1);v*(x+y+z==m)

Tal como se esperaba, los valores de estos coeficientes pueden ser grandes. Aquí tenemos un ejemplo:

Al ser la primitiva definición simétrica respecto a los índices de abajo, el número total de posibles coeficientes trinomiales para un índice superior m dado, es el de particiones de m en tres conjuntos. Por ejemplo, para m=5 las particiones son (5,0,0), (4,1,0), (3,2,0), (3,1,1), (2,2,1), que se corresponden con los posibles valores de los coeficientes, que son 1, 5, 10, 20 y 30. Se comprende que, en el desarrollo de (a+b+c)⁵, se tendrán que repetir, como vemos en el siguiente desarrollo:

Con nuestra herramienta Cartesius se puede encontrar el número de estos coeficientes distintos.(Ver https://www.hojamat.es/sindecimales/combinatoria/herramientas/herrcomb.htm#cartesius)

Por ejemplo, para m=7 plantearíamos:

xtotal=3

xt=0..7

creciente

suma=7

 

Se interpreta que deseamos encontrar todos los tipos de sumas con resultado 7 y sin tener en cuenta el orden (creciente)

Al resultado se le han añadido manualmente los valores de los coeficientes:

Estos resultados se corresponden con el desarrollo de una potencia séptima de un trinomio:

Como era de esperar, los coeficientes aparecen muy repetidos.

Tetraedro de Pascal

Con la expresión COMBINAT(m;n₁)*COMBINAT(m-n₁;n₂) se puede construir un triángulo de coeficientes para un m dado. En una hoja de cálculo se puede visualizar bien. En la imagen figuran los correspondientes a m=7, vistos anteriormente:

Observamos que los lados del triángulo son filas del triángulo de Pascal. Si ahora apilamos triángulos según los valores de m, nos resultará un tetraedro, al que también se le llama de Pascal. En la siguiente página de Wikipedia se explica el proceso:

https://en.wikipedia.org/wiki/Pascal%27s_pyramid

La siguiente imagen pertenece a dicha página:

(Fuente: Wikipedia https://en.wikipedia.org/wiki/Pascal%27s_pyramid)

Para quienes deseéis profundizar en el tema. La lectura de esta página es recomendable.

Recurrencia

Es fácil la comprobación de esta recurrencia para calcular coeficientes trinomiales:

En realidad, este es el fundamento del tetraedro de Pascal.

 

Funciones definidas para tipos de números (3)

 Esta entrada cierra el estudio de funciones definidas para ciertos tipos de números, así como de sus funciones inversas. En esta tercera se estudiarán las definidas por recurrencia.

Funciones sobre recurrencias

 

Son aquellas que se definen mediante los términos iniciales, y además, expresan a(n+1) en función de los anteriores términos. Por ejemplo, la sucesión de Fibonacci se define como a(1)=1, a(2)=1, y para el resto, a(n+1)=a(n)+a(n-1). Si esta última relación sólo usa el término a(n) se llama de primer orden, y así, se pueden plantear de segundo, tercero o más órdenes.

 

Recurrencias de primer orden

 

Cada función de este tipo deberá estar definida por dos parámetros, a(1) y otro según la definición, que puede ser DIFERENCIA en las progresiones aritméticas, RAZÓN en las geométricas, o bien cualquier otra relación entre a(n+1) y a(n). No existen ejemplos populares concretos para definir funciones.

 

Recurrencias de segundo orden

 

Existen muchas recurrencias de segundo orden, como la ya citada de Fibonacci. Las más simples son las llamadas Horadam, y se caracterizan porque la relación de recurrencia es lineal. Vemos las más importantes.

 

Caso particular: sucesiones Horadam

 

Son aquellas que están definidas por cuatro parámetros, a1, a2, c1 y c2, de la forma a(1) = a1, a(2) =a2, a(n) = a(n-2)*c1+a(n*1)*c2. Así, tenemos muchos ejemplos concretos de este tipo de sucesiones recurrentes de segundo grado:

 

Horadam(1,1,1,1): Números de Fibonacci: 1, 1, 2, 3, 5

Horadam(2,1,1,1): Números de Lucas: 2, 1, 3, 4, 7, 11, 18, 29, 47, …

Horadam(1,2,1,2): Potencias de 2:1, 2, 4, 8, 16, 32, …

Horadam(0,1,2,1): Números de Pell: 0, 1, 2, 5, 12, 29, 70, …

Horadam(1,1,2,1): Números de Pell-Lucas: 1, 1, 3, 7, 17, 41, 99, …

Horadam(2,3,3,-2): Sucesión de Pisot (2ⁿ+1): 2, 3, 5, 9, 17, 33, …

Horadam(0,1,1,2): Sucesión de Jacobsthal: 2, 3, 5, 9, 17, 33, …

 

Son muy sencillas de programar, y su búsqueda suele ser rápida. Para encontrar un término general hay que usar los cuatro parámetros citados más el número de orden n.

  

En lenguaje PARI

 

horadam(n,a1, a2, c1, c2) = {my(i,f,t,v,u); if(n==1,f=a1); if(n==2,f=a2);if(n>2,v=a1;t=a2;for(i=3,n,u=c1*t+c2*v;v=t;t=u);f=u);f}

 

En Vbasic

 

Function horadam(n, a1, a2, c1, c2)

Dim i, f, v, u

If n = 1 Then horadam = a1: Exit Function

If n = 2 Then horadam = a2: Exit Function

v = a1: t = a2

For i = 3 To n

u = c1 * t + c2 * v: v = t: t = u

Next i

horadam = u

End Function

 

Si se dispone de estas dos funciones, podemos definir las funciones FIBONACCI(N), PELL(N), LUCAS(N) y demás, igualando su definición a la de Horadam. Por ejemplo, LUCAS(N) se definirá así:

 

Lucas(n)=horadam(n, 2,1,1,1)

 

Así tendremos que lucas(10)=horadam(10,2,1,1,1)=76, como puedes ver en cualquier descripción de estos números. Imagina que puedes definir así todas estas sucesiones.

El autor lleva años usando FIBONACCI(N) en sus cálculos.

 

La función ORDEN se puede construir de similar a la que usamos para Primos o capicúas. Puede ser esta:

 

Function eshoradam(n, a1, a2, c1, c2) As Long

Dim v, u, k, t As Long

Dim es As Long

 

If n = a1 Then eshoradam = 1: Exit Function

If n = a2 Then eshoradam = 2: Exit Function

v = a1: t = a2: k = 2

While t <= n

If t = n Then

es = k

ElseIf t > n Then

es = 0

End If

u = c1 * t + c2 * v: v = t: t = u: k = k + 1

Wend

eshoradam = es

End Function

 

Nos devuelve el orden k si es un término de la sucesión, o un cero si no lo es.

 

El autor usa en su hoja de cálculos diarios esta versión, que le permite averiguar si un número dado es de uno de los tipos populares de esta sección:

 

Function estipohoradam$(n)

Dim s$

Dim es As Long

s = ""

es = eshoradam(n, 1, 1, 1, 1): If es <> 0 Then s = s + " Fibonacci " + " & " + ajusta(es)

es = eshoradam(n, 2, 1, 1, 1): If es <> 0 Then s = s + " Lucas " + " & " + ajusta(es)

es = eshoradam(n, 0, 1, 2, 1): If es <> 0 Then s = s + " Pell " + " & " + ajusta(es)

es = eshoradam(n, 1, 1, 2, 1): If es <> 0 Then s = s + " Pell-Lucas " + " & " + ajusta(es)

es = eshoradam(n, 2, 3, 3, -2): If es <> 0 Then s = s + " Pisot " + " & " + ajusta(es)

es = eshoradam(n, 0, 1, 1, 2): If es <> 0 Then s = s + " Jacobsthal " + " & " + ajusta(es)

es = eshoradam(n, 2, 1, 1, 2): If es <> 0 Then s = s + " Jacobsthal-Lucas " + " & " + ajusta(es)

 

If s = "" Then s = "No es de ningún tipo"

estipohoradam = s

End Function

 

Así, por ejemplo, estipohoradam(76) nos devuelve Lucas  & 10, y estipohoradam(81), No es de ningún tipo.

 

Podemos seguir recorriendo tipos de sucesiones recurrentes, como las de Perrin y otras, pero no suelen demandar nuestros trabajos unas funciones directas como las que aquí se han presentado.

 

Otras recurrencias

 

El contenido de esta entrada nos ayuda si necesitamos definir una función para sucesiones que no sean de Horadam, pero el procedimiento para encontrar funciones de este tipo es el mismo, salvo la fórmula de recurrencia.