24 abril 2019

Termux : Soluciones enteras a la ecuación diofántica k = x³+y³+z³

Hace poco leí un interesante artículo en el  blog "La ciencia de la mula francis "  relacionado con la última solución encontrada a las ecuaciones diofánticas k = x³+y³+z³  para números enteros.

Una ecuación tan sencilla requirió tres semanas de un superordenador para encontrar la primera solución conocida de la ecuación 33=k = x³+y³+z³  ...

Si queremos encontrar soluciones enteras a esta ecuación podemos asumir las siguientes premisas:
  • x, y, z números enteros positivos
  • Opción 1:  k = x^3 + y^3 + z^3  (válida para valores pequeños de x,y,z)
  • Opción 2:  k = x^3 - y^3 + z^3 
  • Opción 3: k = x^3 - y^3 - z^3
Hay números que no pueden ser representados por esta ecuación : 4, 5, 13, 14, 22, 23 ,31 ,33, etc. Son números que se pueden expresar como 9*n+4 ó 9*n+5.

Y por otro lado hay números que pueden ser representados por infinitas combinaciones de x,y,z. Por ejemplo, el número 1 puede expresarse mediante la fórmula:
  1= (1 + 9*m^3)^3 + (9*m^4)^3 - (9*m^4 + 3*m)^3

Para comprobar la dificultad de encontrar soluciones a este tipo de ecuaciones desarrollé un pequeño escript utilizando la Opcion 3 en PARI-GP para números menores a 1000.

Después de un dia de cálculo todavía no he encontrado soluciones para muchos números pequeños ( 7, 12, 24, 26, etc).

Por lo visto no es tan sencilla esta ecuación !!!




23 abril 2019

Termux : Factorización de números enteros 2^n+1 (2<=n<=500)

Al igual que en el post anterior, se calcularon las factorizaciones de los números enteros con formato 2^n+1 en el rango 2<=n<=500.

El programa  se ejecutó en una tablet Praga utilizando  un núcleo de los cuatro disponibles.
El tiempo necesario fué de cuatro días, mayor que en el caso anterior ya que es más dificil la factorización de estos números.

El script tiene varias líneas de comandos similares al siguiente:

gmpfac -is -e100,100000,10000 -c2,250+   (en este caso para el número 2^250+1)

20 abril 2019

Termux : Factorización de números enteros 2^n-1 (2<=n<=500)

Se utilizó el emulador  termux en una tablet Praga para factorizar  números enteros con formato  2^n-1 en el rango 2 <= n <= 500 .

Se creó un fichero script  mediante el sistema  PARI-GP para crear comandos de ejecución utilizando la aplicación gmpfac12 (Conrad Curry) (un solo núcleo). Cada línea tiene el formato:

   gmpfac -is -e100,100000,10000 -c2,250-   (en este caso para el número 2^250-1)

El proceso se demoró 3 días en ejecutarse.
Se puede observar que hay líneas con factores que tiene el prefijo (C). En este caso el programa no pudo factorizar dicho número a pesar de usar el algoritmo de factorización ecm con 100  curvas.

La salida del script se enrutó a un fichero texto el cual fué procesado por un programa creado en el lenguaje  free pascal para obtener el formato final:

Termux: ¿ La suma de varios números primos consecutivos es un número primo ? II

La entrada anterior finalizaba con una pregunta: ¿Cual será el máximo de secuencias consecutivas? Pues la respuesta parece ser que no hay ...