Sucesiones completas

En esta entrada abordaré un tema interesante, y es el de la posibilidad de que los elementos de una sucesión de números naturales, finita o infinita, generen otros números, a partir de sumas, con o sin repetición.

Ya traté ese tema en algunas entradas antiguas, como

https://hojaynumeros.blogspot.com/2010/02/frobenius-y-los-mcnuggets.html

“Frobenius y los mcnuggets”, en la que trataba el problema de las monedas necesarias para alcanzar una cantidad y el número de Frobenius.

https://hojaynumeros.blogspot.com/2012/11/descomposicion-de-un-numero-segun-una.html

En esta entrada se trataba de generar un número a partir de una lista de ellos.

En ambos casos las sucesiones eran finitas, y al sumar los elementos se podían repetir los sumandos.

También se relacionaba algo con lo que pretendo presentar ahora, la entrada sobre algoritmos voraces o codiciosos:

https://hojaynumeros.blogspot.com/2024/09/algoritmos-codiciosos-para-sumas.html

Sucesión completa

Desarrollaré aquí el tipo de sucesión más eficaz para generar todos los números naturales, el de sucesión completa.

Una sucesión de números naturales se llama completa cuando cualquier otro número natural se puede escribir como suma de elementos de esa sucesión sin repetir ninguno.

El ejemplo básico más popular es el de las potencias de 2, que da lugar al sistema binario de representación de los números. En efecto, la sucesión 1, 2, 4, 8, 16, … permite representar cualquier número como suma de algunos términos, sin repetición. La forma de conseguirlo es similar a la de los algoritmos codiciosos que traté en la entrada referenciada más arriba.

Vemos un ejemplo: representar el número 53:

Buscamos el término de la sucesión más cercano a 53 y menor que él, que es el 32, tomamos nota y restamos: 53-32=21. Reiteramos: 21-16=5, 5-4=1, con lo que obtenemos 53=32+16+4+1, que da lugar a la representación binaria 110101.

Este ejemplo nos abre camino para caracterizar las sucesiones completas. Podemos dar tres condiciones para que una sucesión sea de este tipo:

·    El primer término ha de ser 1, a₁=1

·   Aunque no es imprescindible, pero sí algo operativo, la sucesión debe presentar un orden creciente.

·    Si llamamos s_n a la suma de los n primeros términos de la sucesión, se debe cumplir: s_n-1 ≥ a_n-1 para todo n ≥ 1

La primera condición es imprescindible para poder generar todos los números y que los algoritmos tengan parada. Por ejemplo, los números pares no son completos, pues nunca pueden generar un impar.

La segunda se incluye para que tengan sentido los algoritmos.

La tercera es fundamental, porque garantiza que se pueda avanzar en descubrir sumandos hasta llegar, si es necesario al 1.

Seguimos el proceso del 53 visto más arriba. Dado el número N, deberemos encontrar el mayor elemento de la sucesión que sea menor que N, en el ejemplo 32. La diferencia 21 se ha de poder tratar del mismo modo, con lo que s_n-1 ha de ser suficiente para ello, es decir, s_n ≥ a_n-1. Piénsese en el caso 63, en el que 63-32=31, y debe poder ser generado por los elementos anteriores, es decir por s_n-1, y de ahí la condición s_n ≥ a_n-1.

Si no se cumple esta condición, la sucesión no es completa. Por ejemplo, la sucesión de potencias de cuatro, 1, 4, 16, 64, …no la cumple, y, por ejemplo, no puede generar el número 70, porque 70-64=6, y 6 no se puede formar con 1 y 4. Aquí la suma 1+4+16 no es suficiente.

La sucesión de potencias de 2 la cumple, porque s_n=2ⁿ-1, y da lugar a la representación binaria. Esto se logra porque no se admiten repeticiones en los sumandos. También es minimal, pues si eliminamos un término de ella, habrá números naturales que no se puedan representar. Por ejemplo, si elimino el 4, no podré representar el 6. Otra característica es que la representación es única, algo que no se contempla en las definiciones.

Existe una condición suficiente que puede sustituir a la tercera, y es que 2a_k ≥ a_k+1 para k≥1. El ejemplo de las potencias de 2 la cumple, pero falta ver su suficiencia. Lo razonamos con el ejemplo del N=53. La primera operación fue buscar el término menor que 53 más cercano a él, y fue el 32. Restamos y nos resultó 21 para proseguir. Lo vemos en general:

N estará entre a_k y a_k+1, pero a_k+1 estará entre a_k y 2a_k, luego tendremos: a_k ≤ N ≤ a_k+1 ≤ 2a_k. Restando a_k se verificará:

0 ≤ N - a_k ≤ a_k

Por tanto, la siguiente diferencia (en el ejemplo, el 21) será menor que a_k, y se podrá seguir el algoritmo con un término menor, hasta llegar a 1.

Ejemplo reciente

En una entrada anterior vimos la sucesión “catering perezoso”, formada por los términos de fórmula

Sus primeros términos son 1, 2, 4, 7, 11, 16, 22, 29, 37, 46, 56, 67, 79, 92, 106, 121, 137, 154, 172, 191, 211, 232, 254, 277, 301, 326, 352, 379, …https://oeis.org/A000124

Es evidente que cumple 2a_k ≥ a_k+1.

2a_k - a_k+1 = 2(k(k+1)/2+1)-(k+1)(k+2)/2-1 = (k+1)(k-(k+2)/2)+1 =(k+1)(k/2-1)+1, y esta diferencia es positiva o cero para k>1

Podemos afirmar que esta sucesión es completa, como se afirma en la página en la que está publicada. Según esto, cualquier número natural es suma de algunos de sus términos tomados sin repetición. He efectuado sencillas comprobaciones con mi hoja de cálculo partlista

 (Ver https://www.hojamat.es/sindecimales/aritmetica/herramientas/herrarit.htm#reprenum)

He elegido el número 62 al azar, y deseo generarlo con los términos de esta sucesión 1, 2, 4, 7, 11, 16, 22, 29, 37. He concretado en la hoja que no se repitan sumandos, y me ha devuelto estos resultados:

No es de extrañar que resulten siete soluciones, porque en estos procesos no tiene que existir solución única.

 
Sucesión de Fibonacci

Esta sucesión es claramente completa. Basta recordar que es creciente y en ella a_k+1 = a_k+a_k-1, luego a_k+1 < a_k+a_k = 2a_k, que es condición suficiente para la completitud.

Existe una descomposición basada en esta propiedad, y es la de Zeckendorf, que consiste en ir restando a un número N el mayor número de Fibonacci posible. Es interesante, y acudo a ella con frecuencia. Puedes profundizar en el tema leyendo mi entrada https://hojaynumeros.blogspot.com/2020/09/representacion-de-zeckendorf.html y descubrir en qué consiste la multiplicación de Fibonacci.

Por ejemplo, esta es la representación de Zeckendorf del número 75:

75 = F(10)+F(7)+F(5)+F(3) = 55+13+5+2

No es el único desarrollo. Con nuestra herramienta Cartesius (https://www.hojamat.es/sindecimales/combinatoria/herramientas/herrcomb.htm#cartesius) se puede conseguir otro. Estas son las condiciones:

xrango=5 

xt=1..55    

xt=filtro(fibonacci)    

suma=75  

creciente 

no repite  

Y este el resultado:

Nos devuelve el ya conocido, 2+5+13+55, que es el desarrollo minimal, pero nos ofrece otro: 75=34+21+13+5+2

Números primos con la unidad

Si a la sucesión de números primos le adjunto a(1) = 1, resulta una sucesión completa. Es creciente, comienza en 1 y cumple que 2p(k) ≥ p(k+1). Esta afirmación se basa en el Postulado de Bertrand, que afirma, entre otras formulaciones, esta desigualdad. Por ejemplo, con la hoja partlista, podemos generar el número 35 con esta sucesión, y resulta:

Es evidente que esta descomposición da lugar a muchas soluciones distintas.

Se pueden encontrar más ejemplos de sucesiones completas, pero con estos se comprende bien el concepto. Otros que he encontrado son demasiado particulares y no son muy interesantes.

 

El catering perezoso

El título es la traducción literal de la sucesión “Lazy Caterer's Sequence” publicada en https://oeis.org/A000124 con el título Central polygonal numbers, que se refiere al máximo número de porciones que se puede obtener en una tarta si la sometemos a n cortes. Su referencia a los números de Hogben, ya tratados en este blog

(https://hojaynumeros.blogspot.com/2021/11/numeros-de-hogben.html)

 me ha llevado a estudiarla, sólo en los aspectos que trata normalmente este blog.

Fórmula de la sucesión

Si solo nos interesan los números máximos de piezas resultantes, deberemos suponer que cada corte posee intersección con todos los anteriores, evitando paralelismos, por lo que se cumplirá que c(n+1)=c(n)+n. Como c(1)=2, ya que el primer corte produce dos piezas, es fácil ver que la sucesión será: 2, 2+2,2+2+3, 2+2+3+4 ,… pero en OEIS se considera siempre c(0), que aquí sería 1 (ningún corte), por lo que resultaría al final 1, 2, 4, 7, 11, 16, …, que es la publicada.

El hecho de que la suma de los primeros números naturales sea un número triangular, n(n+1)/2, nos da directamente la fórmula de estos números:

En la imagen vemos la generación de c(4), porque de los tres cortes existentes, el cuarto produce 4 nuevos, y el resultado es 4*5/2+1=11

Si la fórmula de c(n) es n(n+1)/2+1, la de c(n-1) será, desarrollando (n²-n+2)/2, y podemos identificar el numerador con un número de Hogben aumentado en una unidad. Esto relaciona las dos sucesiones, pues si añadimos una unidad a la Hogben y dividimos entre 2 nos resulta la que estamos estudiando.

Si emparejamos índices, se comprueba esta relación:

Hogben:             1, 1, 3, 7, 13, 21, 31, 43,  57, 73, 91

Lazy Caterer's:       1, 2, 4,  7,  11, 16,  22,  29, 37, 46

(13+1)/2=7, (21+1)/2=11,  (31+1)/2=16, …

Los números de Hogben, según es fácil de demostrar, se generan con la recurrencia h(n+1)=h(n)+2n, ajustando bien el índice inicial. Es otra forma de razonar la relación entre las dos sucesiones.

Si repasamos mi estudio de estos números, podemos descubrir que los c(n) son la cuarta parte de la suma de dos de h(n)  diferenciados en dos órdenes:

2=(1+7)/4, 4=(13+3)/4, 7=(21+7)/4, 11=(31+13)/4, …

Los primeros términos de la sucesión c(n) son:

1, 2, 4, 7, 11, 16, 22, 29, 37, 46, 56, 67, 79, 92, 106, 121, 137, 154, 172, 191, 211, 232, 254, 277, 301, 326, 352, 379, …

https://oeis.org/A000124

Basta aplicar la formula a los primeros números naturales para conseguirla.

Algunas curiosidades

En la página citada se proponen varias propiedades y curiosidades, pero sin desarrollar. Estudiaremos algunas.

Todo proceso en el que el término enésimo produzca el siguiente al sumarle n puede interpretarse como número del tipo C(n).

Un ejemplo que se propone que el de contar los términos de la expresión (x+y)*(x²+y²)*(x³+y³)*...*(xⁿ+yⁿ).

For n >= 1 a(n) is the number of terms in the expansion of (x+y)*(x^2+y^2)*(x^3+y^3)*...*(x^n+y^n). - Yuval Dekel (dekelyuval(AT)hotmail.com), Jul 28 2003

La justificación es la siguiente: Los primeros términos la cumplen, porque c(0)=1, c(1)=2 y c(2)=4, como fácilmente se comprueba. Los siguientes deberán incrementar los términos en n elementos, porque, por ejemplo, la expresión segunda x³+x²y+xy²+y³ se convertiría en x³*(x³+x²y+xy²+y³)+y³*(x³+x²y+xy²+y³), aparentemente ocho términos, pero se pierde uno al sumar los dos términos del tipo x³y³, con lo que se suman sólo tres términos, quedando 7, que es c(3). Lo vemos con wxMaxima:

         ratsimp((x+y)*(x^2+y^2)*(x^3+y^3));

     y^6+x*y^5+x^2*y^4+2*x^3*y^3+x^4*y^2+x^5*y+x^6

Siete términos, c(3)=7

Se razonaría de igual forma para grado 4. Lo comprobamos:

         ratsimp((x+y)*(x^2+y^2)*(x^3+y^3)*(x^4+y^4));

y^10+x*y^9+x^2*y^8+2*x^3*y^7+2*x^4*y^6+2*x^5*y^5+2*x^6*y^4+2*x^7*y^3+x^8*y^2+x^9*y+x^10

Once términos, c(4)=11

Otra forma de ver las curiosidades es la relación con técnicas de elegir dos elementos en un conjunto. La propuesta de nuestro compatriota Arregui aprovecha ese detalle:

Number of interval subsets of {1, 2, 3, ..., n} (cf. A002662). - Jose Luis Arregui (arregui(AT)unizar.es), Jun 27 2006

Para fijar un intervalo en el conjunto {1, 2, 3, ..., n} bastará elegir un inicio y un final, pero existen n(n-1)/2 formas de efectuar la elección, y añadimos el intervalo vacío, obtendremos los términos C(n). Por ejemplo, {1} sólo posee el intervalo vacío (), {1, 2} los intervalos (), (1,2), (2,3), (1,3), que son c(2)=4,  {1, 2, 3} tendría estos: intervalos (), (1,2), (2,3), (1,3) (1,4), (2,4), (3,4), en total c(3)=7.

Cualquier otro proceso de este tipo, ajustando los índices, se puede representar con los números c(n)

Termino con dos propiedades triviales:

Numbers m such that 8m - 7 is a square. - Bruce J. Nicholson, Jul 24 2017

En efecto, si T(n) es triangular, 8T(n)+1 es un cuadrado, según una popular propiedad, luego 8C(n)-7=8T(n)+8-7=8T(n)+1, que es un cuadrado.

a(n) is the sum of the first three entries of row n of Pascal's triangle. - Daniel T. Martin, Apr 13 2022

En efecto, los tres primeros términos son 1, n y n(n-1)/2 y su suma 1+(n²-n+2n)/2=1+n(n+1)/2=c(n). Por ejemplo, para n=4 queda 1+4+6=11=c(4)

Recurrencia lineal

Las sucesiones dependientes de cuadrados suelen presentar esta recurrencia: a(n+3) = 3*a(n+2) - 3*a(n+1) + a(n). En nuestro caso los términos iniciales son 1, 2, 4. Con mi herramienta de recurrencias se puede comprobar.

(https://www.hojamat.es/sindecimales/aritmetica/herramientas/herrarit.htm#recurre2)

Rellenamos las condiciones iniciales:

Pedimos ver la sucesión:

Como era de esperar, obtenemos C(n)

Como los números de Hogben también dependen de un cuadrado, se obtendrán con una recurrencia idéntica, pero iniciando en 1, 1, 3:

 

Pedimos ver la sucesión y obtenemos H(n):

En nuestra entrada de blog sobre estos números se explica el fundamento de que esta recurrencia sea válida para expresiones de segundo grado.

 

Por una unidad, no son de Fibonacci (2)

En la entrada anterior repasamos las recurrencias lineales no homogéneas, alterando la definición de la sucesión de Fibonacci, ya sea sumando o restando una unidad y cambiando las condiciones iniciales.

En esta segunda entrada restaremos una unidad a la definición clásica, y fijaremos a(1)=1 y a(2)=k, número que variaremos para ver qué ocurre. La recurrencia será a(n)=a(n-1)+a(n-2)-1. Tomaremos k=4 para comenzar. Como en la anterior entrada, generaremos la sucesión directamente en columnas de hoja de cálculo:

Procedemos a encontrar la solución homogénea, y vemos la diferencia entre ambas. Se puede consultar la entrada anterior para seguir este proceso.

Pantalla de la hoja recurre_lineal:

Copio la sucesión de la recurrencia homogénea junto a la obtenida para la no homogénea. Creo una tabla en la que al leer de izquierda a derecha desembocaremos en la expresión definitiva:

En primer lugar figuran los índices de las sucesiones. Aquí comenzamos con 1. Después emparejamos las soluciones homogéneas con las que no lo son, y se puede descubrir que su diferencia es 1-F(0), con lo que quiero expresar que es el número de Fibonacci correspondiente al índice de su fila. En la siguiente columna situamos la solución para la homogénea, deducida de mi hoja recurre_lineal, y, por último, sumo 1-F(0). Al final (ver celda naranja) quedamos con 3F(-1)+1. Por cuestiones de índices, en OEIS usan FIBONACCI(N)+1. Al final da igual donde comencemos.

Se puede comprobar con PARI:

Observa el código tan simple que figura arriba.

Estos resultados coinciden con los publicados en https://oeis.org/A187893

Lo interesante del proceso que hemos seguido es que es válido para cualquier otro valor de K, y la sucesión vendría dada por (k-1)*FIBONACCI(N)+1

Por ejemplo, para k=3 obtenemos 1, 3, 3, 5, 7, 11, 17, 27, 43, 69, 111, 179, 289, 467, 755, 1221, … que sigue la expresión 2*FIBONACCI(N)+1.

Para k=6 obtenemos 1, 6, 6, 11, 16, 26, 41, 66, 106, 171, 276, 446, … definida por 5*FIBONACCI(N)+1:

for(n=0,20,print1(5*fibonacci(n)+1,", "))

1, 6, 6, 11, 16, 26, 41, 66, 106, 171, 276, 446, 721, 1166, 1886, 3051, 4936, 7986, 12921, 20906, 33826, …

La mayoría de casos para k están sin publicar, como era esperable.

 

Casos más simples

Si las complicaciones de la recurrencia, otra forma de añadir una unidad a la sucesión de Fibonacci es realizarlo una vez esta formada. Nos resultan dos sucesiones sencillas:

Fibonacci más una unidad

Usaré PARI, que nos está viniendo muy bien en este tema:

for(n=0,20,print1(fibonacci(n)+1,", "))

1, 2, 2, 3, 4, 6, 9, 14, 22, 35, 56, 90, 145, 234, 378, 611, 988, 1598, 2585, 4182, 6766, …

Está publicada en https://oeis.org/A001611

En esa página destacan que tres elementos pueden formar triángulo, porque cada uno es menor que la suma de los dos anteriores y mayor que su diferencia. También se afirma que los únicos números primos presentes son el 2 y el 3. Podemos añadir algo más:

Entre los primeros términos sólo hay tres cuadrados: 1, 4 y 9.

for(n=0,2000,a=fibonacci(n)+1;if(issquare(a),print1(a,", ")))

1, 4, 9, …

Aparecen cuatro triangulares (observa el código, que se basa en que 8T+1 es cuadrado):

? for(n=0,2000,a=fibonacci(n)+1;if(issquare(8*a+1),print1(a,", ")))

1, 3, 6, 378, …

También contamos con oblongos (tipo N(N+1)):

? for(n=0,2000,a=fibonacci(n)+1;if(issquare(4*a+1),print1(a,", ")))

2, 2, 6, 56, 90, …

 

Fibonacci menos una unidad

Con la experiencia adquirida todo es más fácil, y me limitaré a resultados:

? for(n=1,30,a=fibonacci(n)-1;print1(a,", "))

0, 0, 1, 2, 4, 7, 12, 20, 33, 54, 88, 143, 232, 376, 609, 986, 1596, 2583, 4180, 6764, 10945, 17710, 28656, 46367, 75024, 121392, 196417, 317810, 514228, 832039, …

Publicada en https://oeis.org/A000071, donde se desarrollan múltiples propiedades.

Los únicos primos presentes entre los primeros son 2 y 7:

? for(n=1,3000,a=fibonacci(n)-1;if(isprime(a),print1(a,", ")))

2, 7, …

Cuadrados sólo aparecen 0, 1 y 4

Triangulares 0, 1 y 1596

Oblongos 0, 2, 12 y 20.

Hasta aquí llega el recorrido de este tema. El resto se lo dejo a quien me lea.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Por una unidad, no son de Fibonacci (1)

La sucesión de Fibonacci, definida como a(1) = 1, a(2) = 1, a(n) = a(n-1) + a(n-2), da lugar a otras muchas sucesiones que mantienen básicamente la misma fórmula de recurrencia o bien cambian las condiciones iniciales. Hoy nos vamos a dedicar a estudiar aquellas sucesiones definidas mediante a(n)=a(n-1) + a(n-2) ± 1. Respecto a las condiciones iniciales, elegiré aquellas que estén publicadas o alguna otra que sea afín, porque el objetivo no es descubrir ada nnuevo, sino practicar con las herramientas que suelo usar y explicar lo que se descubra.

Una unidad más

En primer lugar, estudiaré a(n) = a(n-1) + a(n-2) + 1, con a(1) = 0 y a(2) = 2. (Ver https://oeis.org/A001610)

Con una hoja de cálculo se construyen perfectamente las sucesiones definidas por recurrencia. Escribimos a(1), en la celda de abajo a(2), y en las siguientes la fórmula de recurrencia. No hay más dificultad. Es conveniente situar a su izquierda una columna de índices. Por ejemplo, en este caso obtendríamos:

Coinciden los términos con los publicados en https://oeis.org/A001610

0, 2, 3, 6, 10, 17, 28, 46, 75, 122, 198, 321, 520, 842, 1363, 2206, 3570, 5777, 9348, 15126, 24475, 39602, 64078, 103681, 167760, …

Recurrencia lineal no homogénea

Si deseamos profundizar algo más sobre esta sucesión es conveniente tratarla como recurrencia no homogénea. En estos casos se suele resolver la homogénea y después tratar aparte el sumando no lineal, que aquí sería el 1.

La parte homogénea se resuelve mediante la ecuación característica. La base teórica la tienes en cualquier texto o página web, como

https://es.wikipedia.org/wiki/Relaci%C3%B3n_de_recurrencia

Poseo una herramienta descargable para resolver estas ecuaciones con las limitaciones propias de una hoja de cálculo

 (https://www.hojamat.es/sindecimales/aritmetica/herramientas/herrarit.htm#recurre2)

Le suministro los datos de la parte homogénea del problema, que aquí sería la de números de Fibonacci con las condiciones iniciales 0 y 2.

Pulso el botón Resolver para ver la ecuación característica:

La hoja de cálculo no es simbólica, por lo que hay que saber interpretar estos coeficientes en función del número áureo o de la raíz cuadrada de 5:

Aparecen dos soluciones reales, que dan lugar a la expresión de más abajo:

Sería:

Esta fórmula daría lugar a la parte homogénea, que nos da la sucesión que presenta la herramienta usada:

La parte no homogénea es un 1, y según la teoría, habría que someterla a las mismas condiciones que la homogénea, y después encontrar el coeficiente adecuado.

En nuestro caso esta es la situación:

Valores de la homogénea:   0  2  2  4  6  10  16  26  42

Valores encontrados:           0  2  3  6 10  17  28  46 75

Diferencias:                          0  0  1  2  4  7  12  20 33

Las diferencias equivalen a los números de Fibonacci clásicos menos una unidad, y las soluciones de la homogénea son el doble de los mismos de Fibonacci con un índice una unidad menor. Con esta consideración nos liberamos de cálculos algebraicos.

Así, la solución buscada por recurrencia se puede expresar de esta forma, que da el valor directamente:

A(n)=2*F(n-1)+F(n)-1=F(n+1)+F(n-1)-1

Podríamos expresarlo en función de la raíz de 5 en una fórmula demasiado larga. Es preferible expresarla usando la función FIBONACCI(N). Así se puede encontrar con PARI:

 for(i=1,40,a=2*fibonacci(i-1)+fibonacci(i)-1;print1(a,", "))

La web de PARI (https://pari.math.u-bordeaux.fr/gpwasm.html), nos ofrece las mismas soluciones que las publicadas en OEIS, que en las fórmulas incluyen F(n+2)+F(n)-1, porque suelen usar el 0 como primer índice. En la programación en PARI,  G. C. Greubel usa una expresión similar a la obtenida aquí.

0, 2, 3, 6, 10, 17, 28, 46, 75, 122, 198, 321, 520, 842, 1363, 2206, 3570, 5777, 9348, 15126, 24475, 39602, 64078, 103681, 167760, 271442, 439203, 710646, 1149850, 1860497, 3010348, 4870846, 7881195, 12752042, 20633238, 33385281, 54018520, 87403802, 141422323, 228826126, …

Se descubren varios números primos en el listado. Podemos modificar nuestra búsqueda para encontrarlos:

for(i=1,100,a=2*fibonacci(i-1)+fibonacci(i)-1;if(isprime(a),print1(a,", ")))

Estos serían los primeros números primos:

2, 3, 17, 103681, 10749957121, 115561578124838522881

No he encontrado cuadrados, pero si se encuentran números que son una potencia más una unidad:

10, 17, 28, 122, 842, 5777, 39602, 271442, 1860497, 12752042, 87403802, …

Este ha sido un estudio que abre caminos a otras técnicas matemáticas. Es sólo un inicio para quienes deseen profundizar. En la siguiente entrada se completará con algún caso más.