Pandigitales inesperados

Hay varias acepciones de la palabra pandigital. Aquí nos referiremos a los pandigitales completos sin repetición, es decir, que consideraremos pandigitales a aquellos números de 10 cifras en las que figuren todas del 0 al 9 en cualquier orden y sin repetir ninguna. No se admite cero inicial. Excel lo suele suprimir, pero habrá que tener cuidado en otros lenguajes de programación.

Con esta definición es sencillo calcular por Combinatoria cuantos pandigitales de este tipo puede haber. Inténtalo.

Función para detectar pandigitales

Para esta detección se crea un vector de 10 elementos, ci(10), y en cada uno se cuentan las apariciones de una cifra determinada del número. Para que sea pandigital, todos los elementos han de ser iguales a 1. En cualquier otro caso, no lo será.

Función pandigital

Public Function pandigital(a) As Boolean

 

'han de estar todas las cifras sin repetición y sin cero inicial

 

Dim ci(10)’ Contenedor de cifras

Dim i

Dim t As Boolean

 

t = True ‘En principio suponemos que es pandigital

n = numcifras(a) ‘Se puede conseguir también sumando 1 a su logaritmo decimal

If n <> 10 Then pandigital = False: Exit Function ‘No tiene diez cifras

For i = 0 To 9: ci(i) = 0: Next I ‘Contadores a cero

For i = 1 To n

ci(cifra(a, i)) = ci(cifra(a, i)) + 1’ Se cargan las cifras encontradas

If ci(cifra(a, i)) > 1 Then pandigital = False: Exit Function ‘Si hay repetidas, no es pandigital

Next i

s = 0

For i = 0 To 9

If ci(i) = 0 Then t = False ‘Si falta una cifra, no es pandigital

Next i

pandigital = t ‘Devuelve VERDADERO o FALSO

End Function

 

Resultados inesperados

 

La idea de hoy es buscar pandigitales en distintas circunstancias u operaciones. Predominarán productos y potencias, como ejemplos favorables a encontrar pandigitales, pero sólo serán ejemplos de uso para facilitar a los lectores la búsqueda en otros casos similares. Se pueden abordar muchos.

 

Vemos algunos ejemplos con nuestro buscador:

 

Multiplicar un número por otro relacionado

 

A partir de este ejemplo hay que tener en cuenta que los pandigitales están comprendidos entre 10^9=1000000000 y 10^10=10000000000. Sería inútil buscar fuera de este intervalo.

 

Multiplicamos N por N+1

 

En este primer ejemplo, podemos comenzar la búsqueda por RAIZ(10^9)-1, para que N(N+1) entre en ese intervalo. La salida puede ser RAIZ(10^10)+1. Esto no es demasiado importante. Se puede ampliar algo el rango. Nuestro buscador nos da 52 soluciones. Estas son las primeras:

No esperaba tantos resultados. Se ve que el producto de consecutivos tiene características favorables. He probado con N(N+2) y resultan 34. Con N(N+3), 99. Como era de esperar, es algo que resulta aleatorio en la práctica. Con otras diferencias se obtienen números similares.

 

Un número multiplicado por su doble produce resultados parecidos.

 

Uso de potencias

 

Cuadrado de N

 

Si lo reducimos a su cuadrado, es decir un número por sí mismo, está publicado en https://oeis.org/A054038 con otros planteamientos, pero coincide mi buscador para pandigitales sin repetición, que resultan ser 87.

 

Estos son los primeros:

 

Podíamos operar N con N², a ver qué ocurre.

Con suma aparecen bastantes resultados. Uno es:

40508+40508^2=1640938572

Con cubos no he obtenido ningún pandigital.

Otras potencias

Para las potencias es preferible usar el lenguaje PARI, porque maneja las cifras con más soltura que Excel.

 

El criterio para saber si un número es pandigital puede ser este:

 

vecsort(digits(n))==[0..9] &&n>=10^9

 

Se interpreta fácilmente: se ordenan los dígitos de n en un vector y ha de coincidir con [0..9]. Además, no puede tener cero inicial, ha de ser mayor o igual que 10⁹.

 

Con estas líneas de PARI podemos investigar potencias. Están adaptadas a cuadrados para probar, escribiendo k=2. Cambiando k se pueden investigar otras potencias.

pandigital(n)=vecsort(digits(n))==[0..9]&&n>=10^9

k=2

a=truncate(10^(9/k))

b=truncate(10^(10/k))

for(i=a,b,m=i^k;if(pandigital(m),print(i," , ",m)))

Este es un recorte de pantalla para k=2 en la web de PARI:

Observamos que coincide con el resultado de Excel.

No he encontrado pandigitales en los exponentes del 3 al 10, pero nos ha servido para iniciar el uso de PARI.

Otras búsquedas

Número por su simétrico en cifras

Aquí sí tendremos éxito. En Excel uso mi función CIFRAINVER y en PARI mi función REVERSE. Elijo este porque lo tengo desarrollado más brevemente. Podríamos usar lo siguiente:

pandigital(n)=vecsort(digits(n))==[0..9]&&n>=10^9

reverse(n)=eval(concat(Vecrev(Str(n))))

a=2

b=1000

for(i=a,b,m=i*reverse(i);if(pandigital(m),print(i," , ",m)))

Definimos pandigital y reverse y después multiplicamos el número por su reverso:

Si multiplicas cualquier número de la primera columna por su simétrico, obtendrás el pandigital de la segunda.

Otros ejemplos sin desarrollar

Un número multiplicado por su doble:

Por ejemplo. 22887*45774=1047629538

Cuadrado de N menos N

Así, 44028^2-44028=1938420756

Todos terminarán en 0, 2 o 6.

Y hasta aquí llega el tema. Quedan muchas posibilidades, e incluso eligiendo las primeras cifras de un número real, pero aquí sólo se estudian los enteros.

 

Un número como diferencia entre potencias (2)

En la anterior entrada sobre este tema se usó la descomposición en factores primos de un número. Esto supone declarar los vectores primo() y expo() de forma global y puede resultar un poco técnico. En esta segunda entrada consideraré otras alternativas.

Si no te quieres complicar con factores primos, puedes usar esta función, que te da las primeras soluciones y resulta rápida:

Function dospoten$(n)

 

Dim i, p, q, a, b, c, x

Dim s$

 

s$ = ""

x = n

For p = 2 To 12 ‘Se buscan potencias de exponentes del 2 al 12. ‘Esto se puede cambiar

a = Int(x ^ (1 / p)) + 1 ‘Posible raíz p-nésima

For i = a To 10 * a

b = i ^ p – n ‘Se halla la diferencia con N

For q = 2 To 20

c = Int(b ^ (1 / q) + 0.000001) ‘Se identifica como potencia

If c ^ q = b Then ‘Si todo va bien, se publica

s$ = s$ + "# " + ajusta(i) + "^" + ajusta(p) + "-" + ajusta(c) + "^" + ajusta(q) ‘Si no se cuenta con la función AJUSTA, se puede usar Str$

End If

Next q

Next i

Next p

If s = "" Then s = "NO"

dospoten = s$ ‘Publica la lista de diferencias de potencias

End Function

Con alguna limitación, es útil, y yo la uso en mis cálculos diarios. Puede que no llegue a alguna solución, pero suele funcionar bien.

Unos ejemplos:

60=8^2-2^2# 16^2-14^2# 4^3-2^2# 4^4-14^2# 2^6-2^2# 2^8-14^2

112=# 11^2-3^2# 12^2-2^5# 16^2-12^2# 29^2-27^2# 29^2-9^3# 29^2-3^6# 8^3-20^2# 4^4-12^2# 2^7-4^2# 2^7-2^4# 2^8-12^2# 2^9-20^2# 2^11-44^2

13# 7^2-6^2# 16^2-3^5# 17^3-70^2# 4^4-3^5# 2^8-3^5

Uso de la herramienta CARTESIUS

Mi máquina de construir productos cartesianos condicionados también nos puede ayudar a encontrar diferencias de potencias dentro de un nivel elemental. En ella se puede exigir que las columnas a combinar sean potencias, lo que facilita su uso en esta cuestión. La puedes descargar desde

https://www.hojamat.es/sindecimales/combinatoria/herramientas/herrcomb.htm#cartesius

El problema es que la actualizo a veces, pero no la versión publicada, pero aquí valdrá.

La programación de esta búsqueda sería esta, concretada para el número 28:

xtotal=2

xt=1..100

xt=etiq(potencia)

ES x1-x2=28

 

Declara que se combinarán dos columnas (X1 y X2)

En ellas se escribirán 100 datos (pueden ser más a costa del tiempo de proceso)

Cada columna contendrá las primeras potencias.

Se deberá cumplir que X1 menos X2 sea 28 (esto es un ejemplo)

Al pulsar sobre el botón Iniciar se construirán las columnas de potencias:

En la siguiente página Producto obtendremos las primeras diferencias de potencias que equivalen a 28.

 

Esto es un simple complemento, para comprobar otros resultados y para conocer una posibilidad alternativa. No tiene más importancia.

Esta sería la solución para 28 con las anteriores herramientas:

# 6^2-2^3# 8^2-6^2# 4^3-6^2# 8^3-22^2# 37^3-225^2# 37^3-15^4# 2^5-2^2# 2^6-6^2# 2^7-10^2# 2^9-22^2

Las soluciones de Cartesius están contenidas en esta lista.

Uso del Buscador de Naturales

Mi herramienta Buscador de Naturales posee la función POTENCIA, tanto de forma booleana, como True o False como actuando sobre un parámetro. Esto nos permite plantear una búsqueda de diferencia de potencias de forma sencilla. Basta condicionar N a que sea potencia y después volverla a usar pero con parámetro N+K.

Lo vemos más claramente con el número 24 como ejemplo:

Las condiciones han sido:

POTENCIA: Busca los números que son potencias

ES POTENCIA(N+24): La partícula ES anuncia que se exigirá una condición complementaria mediante expresiones o funciones. Aquí es que N+24 sea también potencia.

EVALUAR: Indica qué se escribirá en la segunda columna. Aquí hemos pedido que sea N, una coma, y N+24

Puedes descargar esta herramienta desde

https://www.hojamat.es/sindecimales/divisibilidad/herramientas/herrdiv.htm#buscador

No resulta muy lento. Le hemos exigido tratar el 28 hasta un tope de 20000 y no ha tardado demasiado:

Coincide con anteriores planteamientos:

28 =2^5-2^2=6^2-2^3=2^6-6^2=2^7-10^2=2^9-22^2

Se han creado dos columnas que después se pueden copiar en otra parte de Excel y añadirles detalles:

Despedida simplista

¿No queremos seguir con algoritmos? Pues podemos construir una tabla de doble entrada en la que figuren las potencias ordenadas tanto en filas como en columnas, y usar las búsquedas de Excel para localizar el número deseado. Se podía comenzar el tema por aquí, y el trabajo pesado que da puede justificar el pensar en algoritmos.

Aquí se pueden localizar algunos de los ejemplos que se han usado. Los he destacado en rojo:

Un número como diferencia entre potencias (1)

Antecedentes

En este blog se han estudiado varios tipos de diferencias de potencias. Por ejemplo, ya hemos visto aquí las diferencias de cuadrados, como puedes buscar en

https://hojaynumeros.blogspot.com/

También se han tratado las de cubos:

https://hojaynumeros.blogspot.com/2024/09/diferencias-de-cubos-enteros-positivos.html

O bien de cubos menos cuadrados y también el caso en el que los dos exponentes son iguales:

https://hojaynumeros.blogspot.com/2025/01/regresos-13-diferencia-de-potencias.html

Por último, algunas curiosidades sobre diferencias de potencias, como, por ejemplo. Acercamientos entre potencias

https://hojaynumeros.blogspot.com/2018/10/acercamiento-entre-potencias.html

En mi reciente publicación “Números y potencias” están contenidas varias cuestiones interesantes (descargable desde

https://www.hojamat.es/publicaciones/numpot.pdf)

Diferencias entre potencias de distinto exponente

No he abordado nunca (creo que es así, porque a mi edad la memoria falla), la diferencia entre dos potencias de distinto exponente. Aprovecharé el tema para presentar varias alternativas.

Función ESPOTENCIA

Hay varias formas de identificar una potencia perfecta. La más directa es la de descomponer el número en factores primos y encontrar el MCD de sus exponentes. Si este es mayor que 1, estaremos ante una potencia. Por ejemplo, 576=2⁶*3². EL MCD de 6 y 2 es 2, luego 576 es un cuadrado. Esta técnica nos da el mayor exponente si el número es multipotencia, como el 64, que lo catalogará como sexta potencia. Para la cuestión actual lo que nos interesa es que sea verdaderamente una potencia.

Para sacar los factores primos disponemos de nuestra función SACAPRIMOS, que la puedes encontrar en

https://hojaynumeros.blogspot.com/2022/06/numeros-con-cifras-crecientes-o_02105820646.html

Esta función te devuelve los factores primos en el vector primo() y sus exponentes en expo(), además de indicar el número de esos factores en NUMOMEGA.

Esas variables se han de declarar de forma global en VBasic:

Global primo(50), expo(50)

Global numomega

 

Esta es la función ESPOTENCIA:

Public Function espotencia(n)

Dim i, j, s, p

 

If n = 1 Or n = 0 Then espotencia = 1: Exit Function

p = n

j = sacaprimos(p) ‘Explicado en el párrafo anterior

s = expo(1) ‘El vector expo está también explicado

If j > 1 Then

For i = 2 To j

s = mcd(s, expo(i)) ‘Calcula el MCD

Next i

End If

If s = 1 Then s = 0 ‘Si MCD es 1, no es potencia

espotencia = s ‘Devuelve el exponente o un cero

End Function

 

Búsqueda de diferencias con la función ESPOTENCIA

Haré uso de mi función ESPOTENCIA y de la correspondiente ISPOWER del lenguaje PARI. Esta última me permitirá llegar a números grandes con más facilidad.

La idea para buscar diferencias de potencias que coincidan con N será la de recorrer todos los números K hasta un tope prefijado, y estudiar si K y K+N son potencias. En caso afirmativo, habremos dado con una solución. Se les halla la raíz exacta y se incorporan al resultado.

La función quedaría así:

Function difepote$(n, tope)

Dim b, m, p, q, r, t

Dim s$

 

s = "" ‘Contenedor de soluciones

m = 0 ‘Contador de soluciones

For b = 1 To tope – n ‘El tope es parámetro prefijado

q = espotencia(b) ‘Daría cero o un exponente

p = espotencia(b + n) ‘Daría cero o un exponente

If p > 0 And q > 0 Then ‘Si no dan cero, es que son potencias

r = raiz_exacta(b + n, p)

t = raiz_exacta(b, q)

m = m + 1 ‘Nueva solución

s = s + "=" + ajusta(r) + "^" + ajusta(p) + "-" + ajusta(t) + "^" + ajusta(q)

End If

Next b

If s = "" Then s = "NO" Else s = ajusta(m) + ": " + s

difepote = s

End Function

La función RAIZ_EXACTA también se basa en los vectores PRIMO y EXPO. Para quien tenga interés, adjunto su código sin comentarios:

Function raiz_exacta(n, k)

Dim p, r, i

 

If n = 1 Or n = 0 Then raiz_exacta = n: Exit Function

p = espotencia(n)

If p Mod k = 0 Then

r = 1

For i = 1 To numomega

r = r * primo(i) ^ (expo(i) / k)

Next i

End If

If r = 1 Then r = 0

raiz_exacta = r

End Function

Un ejemplo del uso de esta función DIFEPOTE lo vemos con el número 28 y tope 200000:

Difepote(28;200000)=7:=2^5-2^2=6^2-2^3=2^6-6^2=2^7-10^2=2^9-22^2=37^3-15^4=2^17-362^2

Nos indica que ha encontrado siete soluciones y escribe sus desarrollos.

No es un algoritmo rápido. En el caso del 28 sí devuelve el resultado en pocos segundos, pero con un tope mayor puede tardar.

Elijo ahora el número 120 para destacar dos resultados que pueden resultar:

120=11^2-1^1=2^7-2^3=13^2-7^2=17^2-13^2=31^2-29^2

Una de las bases puede ser 1 (en ese caso tendrá exponente 1) y las dos bases pueden ser iguales, como en 2^7-2^3. Es normal que esto ocurra.

Versión en PARI

Para lograr más rapidez y llegar a números más grandes, he traducido la función DIFEPOTE al lenguaje PARI. No necesita explicación, salvo que en él es preferible presentar los resultados como lista y no como cadena de texto.

Una novedad en este lenguaje es que ISPOWER, además de devolver el exponente, guarda la base en una variable. Así, ispower(343,,&k) devolverá 3, porque es un cubo, y guardará la base 7 en la variable k. Tiene el inconveniente de que no maneja bien al caso del 1, pues lo considera siempre potencia de exponente cero. Esto se puede ignorar o intentar corregirlo en el código. He elegido esta última posibilidad, aunque el código se alarga demasiado. Es este:

difepote(n,tope)={my(b,b1,s=List(),r1,r2);p=ispower(n+1,,&r1);if(p>0,listput(~s,r1);listput(~s,p);listput(~s,1);listput(~s,1);listput(~s,"#"));for(b=2,tope-n,q=ispower(b,,&r2);b1=b+n;p=ispower(b1,,&r1);if(p<>0&&q<>0,listput(~s,r1);listput(~s,p);listput(~s,r2);listput(~s,q);listput(~s,"#")));s}

print(difepote(132608,2000000))

Aquí se trocea en líneas, pero hay que escribirlo sin saltos de línea. Este es el resultado en la web de PARI:

He elegido este ejemplo porque en él se destaca la mayor potencia de cálculo de PARI, ya que Excel tardaría más. Una solución más restringida, con un tope menor, es

4: =372^2-76^2=52^3-20^3=387^2-131^2=408^2-184^2

Si llegamos a un tope de 2000000, tarda sus buenos minutos en Excel, y al final devuelve:

10:=372^2-76^2=52^3-20^3=387^2-131^2=408^2-184^2=522^2-374^2=582^2-454^2=648^2-536^2=933^2-859^2=1068^2-1004^2=1212^2-34^4

El mismo resultado, pero con más lentitud.

Es una preferencia personal el pasar o no a PARI en esta cuestión.

Seguiré con el tema en la siguiente entrada.

Polidivisibles

Un número se llama polidivisible (aquí se limitará el estudio a base 10) cuando al recorrer sus cifras de izquierda a derecha, las dos primeras forman un múltiplo de 2, las tres primeras, de 3, las cuatro de 4, y así sucesivamente. Por ejemplo, 126006 es polidivisible, porque 12=2*6, 126=3*42, 1260=4*315, 12600=5*2520, 126006=6*21001. Se supone que no se han escrito ceros a la izquierda, o lo que es lo mismo, que la primera cifra es no nula.

Puedes comprobar la afirmación de Wikipedia de que 381654729 es polidivisible.

Si se entiende bien el troceado de cifras, no es difícil crear una función que determine si un número es polidivisible. Propongo esta para VBasic de Excel, fácilmente traducible a otros lenguajes:

Function polidivisible(n) As Boolean ‘Devuelve verdadero o falso

Dim m, i, t

Dim es

 

 

m = numcifras(n)’Función contenida en mi blog “Números y hoja de cálculo”. Es fácil de copiar.

If m < 2 Then polidivisible = False: Exit Function ‘Caso de una cifra

i = 2 ‘Número de cifras primeras a considerar

es = True ‘Suponemos que sí es polidivisible

While i <= m And es ‘Recorremos las primeras cifras

t = Int(n / potencia(10, m - i)) ‘Trozo de cifras

If t / i <> t \ i Then es = False ‘Ha de ser múltiplo de su número de cifras

i = i + 1

Wend

polidivisible = es

End Function

Esta función te devuelve VERDADERO si el número es polidivisible. Con ella y un buscador podemos encontrar los primeros números polidivisibles. En la lista se descubren los comprendidos entre 100 y 200:

102, 105, 108, 120, 123, 126, 129, 141, 144, 147, 162, 165, 168, 180, 183, 186, 189

Los puedes comprobar en https://oeis.org/A144688. Allí se les llama “magic”. Su definición sugiere que si un número es polidivisible, también lo son los trozos de cifras que nos han servido para la definición. Eso ocurre con 20445, que 20, 204 y 2044 también son polidivisibles. Más adelante estudiaremos el proceso contrario, extender un polidivisible a más cifras.

 

Versión en PARI

Excel no es útil para números enteros de muchas cifras, ya que pasa automáticamente al formato científico. Por ello es conveniente el uso de una función en PARI. La de arriba se traduce fácilmente a esta otra:

k=381654729

polidivisible(n)={my(m=#digits(n),i=2,es=1,t);while(es==1&&i<=m,t=truncate(n/10^(m-i));if(t%i<>0,es=0);i+=1);es}

print(polidivisible(k))

Escribimos un valor de k en la primera línea y nos devolverá un 1 si k es polidivisible o un 0 si no lo es. Lo he probado en la web de PARI:

Nos devuelve un 1 porque 381654729 sí es polidivisible.

Añadiendo un bucle podremos buscar polidivisibles en un rango. En la imagen figuran los primeros a partir de 2000000:

Mentalmente se puede comprobar alguno de ellos.

En https://en.wikipedia.org/wiki/Polydivisible_number puedes encontrar un procedimiento sencillo para extender la definición a cualquier base b.

Extensión de un número polidivisible

Si N lo es, podemos intentar añadirle otra cifra más y probar si también es polidivisible. Si N posee k cifras, el siguiente polidivisible estará entre 10*N y 10*N+9 y deberá ser múltiplo de k+1. Si k no es mayor que 10, siempre existirá un múltiplo en ese rango. En valores superiores no se puede garantizar la extensión.

Uso PARI para garantizar un buen número de cifras:

polidivisible(n)={my(m=#digits(n),i=2,es=1,t);while(es==1&&i<=m,t=truncate(n/10^(m-i));if(t%i<>0,es=0);i+=1);es}

extenpoli(n)={my(g=#digits(n),nn=0,i=0,es=0);while(i<=9&&es==0,nn=i+n*10;if(polidivisible(nn),es=1);i+=1);nn*(es==1)}

print(extenpoli(126006))

En primer lugar, vuelvo a definir la función polidivisible. Después, extenpoli, que, como vemos, recorre una cifra adicional (i+n*10) para encontrar un nuevo polidivisible. Si no lo encuentra, devuelve un cero.

En el caso de 126006 nos devuelve un múltiplo de 7. Lo vemos en la web de PARI:

Sería 1260063 el siguiente polidivisible, divisible entre 7. Si proseguimos, el siguiente resulta ser 12600632, divisible entre 8. Seguirían 126006327, 1260063270, 12600632704, 126006327048. Aquí se detiene la secuencia, porque al añadir otra cifra no se logra un múltiplo de 13.

Todo esto funciona si el primer número es polidivisible con seguridad. Si no, llegaríamos a resultados erróneos. Estas secuencias de extensiones no tienen que ser únicas. Un polidivisible puede dar lugar a dos o más cuando se le añade una cifra, ya que en un rango de 10 puede haber dos o más múltiplos del número de cifras.

Secuencia de extensiones

Si usamos una lista, es posible lograr que PARI nos devuelva el conjunto de extensiones (o uno de ellos). Bastará añadir otra función nueva que descubra extensiones mientras sea posible. Todo el conjunto quedaría así:

polidivisible(n)={my(m=#digits(n),i=2,es=1,t);while(es==1&&i<=m,t=truncate(n/10^(m-i));if(t%i<>0,es=0);i+=1);es}

extenpoli(n)={my(g=#digits(n),nn=0,i=0,es=0);while(i<=9&&es==0,nn=i+n*10;if(polidivisible(nn),es=1);i+=1);nn*(es==1)}

secuenpoli(n)={my(m=n,s=List());while(m<>0,listput(~s,m);m=extenpoli(m));s}

print(secuenpoli(126006))

 

Nos daría el mismo resultado para 126006 en forma de lista:

 

List([126006, 1260063, 12600632, 126006327, 1260063270, 12600632704, 126006327048])

 

Otros ejemplos:

 

N=62765

List([62765, 627654, 6276543, 62765432, 627654321, 6276543210, 62765432103, 627654321036, 6276543210366, 62765432103668])

 

N=180402

List([180402, 1804026, 18040264, 180402642, 1804026420, 18040264209, 180402642096, 1804026420963, 18040264209632])

N=747

List([747, 7472, 74720, 747204, 7472045, 74720456, 747204561, 7472045610, 74720456107, 747204561072])

Ejemplo de los primeros párrafos:

N=381654729

List([381654729, 3816547290, 38165472906, 381654729060, 3816547290608])

Con esto terminamos la introducción a este tipo de números. Quedan muchas curiosidades, que puedes encontrar en las dos webs enlazadas más arriba. Aquí nos quedamos con las que son fácilmente tratables con VBasic y PARI.