Números y cifras

Números y cifras

Edición Octubre-2013

Colección [Link]
© Antonio Roldán Martínez
[Link]

1
P RESENT AC IÓ N

El sistema de numeración no es una cuestión importante en la teoría

matemática, pero da lugar a temas divertidos, curiosos, que no hay que
tomar demasiado en serio, pues todos sus resultados tendrán siempre
un ámbito limitado. No obstante, entretienen y son fuente de desarrollo
de habilidades matemáticas. Un reto importante es el de intentar
trasladar las propiedades a otros sistemas de numeración.

Como curiosidades trataremos los temas, pero en ocasiones nos

derivarán a cuestiones mucho más profundas. En esta edición
incorporamos temas nuevos en la sección de Primos y múltiplos.
Algunos introducen temas estadísticos, como el relativo a Pandigitales,
cromos y Benford. El de “Como en casita…” es una extensión de otro
anterior, al que ha englobado. Similar a este es el del mayor divisor
propio.

Como advertiremos en todos los documentos de esta colección, el

material presentado no contiene desarrollos sistemáticos, ni pretende
ser un manual teórico. En cada tema se incluirán cuestiones curiosas o
relacionadas con las hojas de cálculo, con la única pretensión de
explicar algunos conceptos de forma amena.

2
T AB L A D E C ON TE N ID O

Presentación............................................................................................. 2
Operaciones con cifras ............................................................................ 4
La suma divide la concatenación ............................................................ 4
Algoritmo derivado de un problema ........................................................ 7
Función “Dígitos” ...................................................................................10
Número más la suma de sus cifras ........................................................12
Igualdades curiosas ................................................................................15
Casi narcisistas .....................................................................................15
Cuadrado del simétrico..........................................................................16
Suma pandigital ....................................................................................17
Aprovechando las cifras ........................................................................17
Las primeras, doble de las segundas ....................................................18
Números automórficos ..........................................................................20
El fósil de un número.............................................................................23
Con primos y múltiplos ...........................................................................25
Primos con cifras consecutivas..............................................................25
Primos reversibles (Primo-Omirp)..........................................................27
Múltiplos de 11 ......................................................................................30
Damos vueltas a primos y al 18 .............................................................31
Los múltiplos acunan.............................................................................36
Como en casita en ningún sitio..............................................................39
Mayor divisor propio con la misma suma de cifras .................................46
Pandigitales, cromos y Benford .............................................................49
Soluciones ...............................................................................................56
Igualdades curiosas ..............................................................................58
Con primos y múltiplos ..........................................................................64
Apéndice ..................................................................................................67

3
O PER AC IO NES CON CIFR AS

L A SUMA DI VI DE L A CO NCAT EN ACI Ó N

El blog Números de Claudio Meller nos presentó el día 2 una

interesante propuesta:

La suma divide la concatenación

1+2 divide a 12 - 12/3=4

4+5 divide a 45 - 45/9=5
16+17 divide a 1617 - 1617/33=49
49+50 divide a 4950 - 4950/99=50

¿Cuáles son los siguientes números consecutivos tal que la suma de

ellos divide a la concatenación de los mismos?

Aunque desde este blog le enviamos un comentario con posibles

soluciones, parecía interesante aprovechar esta cuestión para recorrer
un razonamiento mixto (hoja de cálculo y Álgebra) en esa búsqueda.
Usamos el proceso Exploración – Conjetura – Demostración de la conjetura –
Complementos, que siempre hemos recomendado en los procesos de
investigación en el aula de Matemáticas.

Exploración

Al tratar de números consecutivos y dos operaciones sencillas, era

atractivo organizar una búsqueda con una hoja de cálculo. Bastaba
crear una tabla similar a la siguiente:

Número Consecutivo Concatenación Suma Cociente

…
164 165 164165 329 498,98
165 166 165166 331 498,99
166 167 166167 333 499
167 168 167168 335 499,01
168 169 168169 337 499,02
…
y esperar a que aparecieran números enteros en el cociente.

4
La concatenación se programó con fórmulas del tipo 10^N*a+a+1,
siendo N el número de cifras de a. De esta forma fueron apareciendo
las soluciones 1, 4, 16, 49, 166, 499, 1666, 4999,
…([Link]

Conjetura

A la vista de los resultados, parecía que las soluciones eran de dos

tipos:

A1=5*10n-1 y A2=(5*10n-2)/3

Y que los cocientes siempre estaban comprendidos entre 5*10n-2 y

5*10n+1

¿Sería siempre así?

Demostración

El cociente estudiado entre concatenación y suma se puede

representar por la expresión

Donde N es el número de cifras de a. Este cociente siempre está

cercano al número 5*10N-1. Precisemos más:

En efecto:

Luego los cocientes no llegarán a 6, 51, 501, 5001, 50001,…

Por otra parte

5
Esto hace que los cocientes enteros puedan ser también del tipo 48,
49, 498, 499,…

Así que tenemos tres posibilidades, aunque la primera no ha aparecido

en la experimentación. Seguro que se puede lograr una acotación más
fina (ver comentario de Robert Israel en A173712)

K1=5*10N-1, K2=5*10N-1-1, K1=5*10N-1-2

Primer caso:

Nos lleva, al despejar la incógnita a a la expresión: a=5*10 N-1-1

que nos da las soluciones 4, 49, 499, 4999, 49999,…

Segundo caso

Despejando a tendremos

Que siempre da un resultado entero, porque 5*10N-1 es congruente

módulo 3 con 2 (¿por qué?) y nos devuelve las soluciones 1, 16, 166,
1666, 16666,…

Tercer caso

Dejamos como ejercicio ver que no puede dar solución entera.

6
AL G O RI T MO DERI VADO D E UN PRO BL EMA

2758620689655172413793103448 * 3=8275862068965517241379310344

¿Qué tiene de particular este resultado?

Hace días leí en un libro de problemas el siguiente:

Encontrar un número entero que termine en 6, y que si esa cifra 6

se mueve hasta situarse delante del resto de las cifras del
número, el resultado equivalga a multiplicar ese número por 4:
6abc..de=4*abc..de6

Un caso similar es el número 205128, que si movemos el 8 a la

primera posición 820512 el resultado equivale a cuatro veces el
primitivo: 205128*4=820512

¿Cuál es la forma más rápida de resolver este tipo de problemas?

Intenté analizar el problema propuesto por la parte izquierda, y aunque

llegué a alguna solución, vi que resultaba mucho más eficiente trabajar
por las unidades, después las decenas, etc. En efecto, si las unidades
son 6, al multiplicar por 4 han de resultar 4 unidades. Luego el número
termina en 4. Por un razonamiento similar, las decenas han de valer 8
(4*4+2=18), las centenas…seguí así hasta encontrar la solución.
¿Puedes encontrarla tú?

Este razonamiento se puede convertir en un algoritmo, en el que dada

la cifra de las unidades (la que ha de moverse) y el número a
7
multiplicar te devuelva el resultado, si es que existe. El problema es
que hay que darle dos condiciones de parada:

(a) Debe aparecer la cifra buscada

(b) El algoritmo se detiene cuando la cifra multiplicadora es igual a
cero.
Si lo implementamos en hoja de cálculo el límite es el número de filas o
columnas.
Si te animas a encontrar un procedimiento de resolución e incluso a
convertirlo en algoritmo, intenta conseguir resultados tan
espectaculares como el que encabeza esta entrada.

Aquí tienes otro resultado del algoritmo

Equivale a encontrar que

1304347826086956521739 * 7=9130434782608695652173
En las soluciones puedes ver otros casos con más cifras

¿Existirán datos que produzcan algoritmos sin parada?

Hasta donde hemos experimentado, siempre se llega a la cifra primitiva
y a un multiplicador cero que detiene el algoritmo.

Planteamiento algebraico

Una de las características más elegantes de las Matemáticas es la

concurrencia de resultados procedentes de métodos muy distintos.

8
Basta recordar, por ejemplo, las demostraciones que existen del
teorema de Pitágoras procedentes de planteamientos muy variados.

Intentaremos un planteamiento del problema de tipo más algebraico.

En lugar del 6 consideraremos cualquier cifra a entre 2 y 9 (el 0 y el 1
dan casos triviales) y llamaremos b al número formado por el resto de
cifras. Así, el número que debemos buscar se puede expresar como
10b+a. Llamemos N al factor por el que hay que multiplicar ese
número. Si se cumplen las condiciones del problema se podrá escribir

N(10b+a) =10xa+b, siendo x el número de cifras de b

Despejando: (10N-1)b=(10x-N)a

Serán soluciones del problema (no todas) las procedentes de la

condición de que 10N-1 sea divisor de 10x-N, que estará formado por
varios 99…9 seguidos de la cifra 10-N.

Podemos ir probando los valores de N, con lo que 10N-1 irá teniendo el

valor de 9, 19, 29,…89, y deberán ser divisores de 10x-N. Es fácil
programarlo en una hoja de cálculo. En la siguiente tabla se han
descubierto tres posibilidades (si siguiéramos encontraríamos más)

El primer 1 corresponde a N=4 y expresa que 99996 es divisible entre

39. Si multiplicamos su cociente 2564 por las distintas cifras, nos
resultarán valores de b, y por tanto soluciones del problema: 102564*4
= 410256; 128205*4 = 512820; 153846*4 = 615384 (Esta es la solución
para el enunciado de arriba). Intenta encontrar más.

El segundo 1 nos da otras soluciones con más cifras:

99999999996/39=2564102564, que al multiplicar por a y añadirle una
cifra nos devuelve más soluciones, por ejemplo:

128205128205*4 = 512820512820
9
Prueba a obtener soluciones del último 1, que corresponde a N=8

Este método no agota las soluciones, pues 10N-1 puede tener factores
comunes con a que alteren las condiciones, pero como diversión ya
está bien con lo estudiado.

F UNCI Ó N “ DÍ G I TO S”

Si escribimos la serie de números 1, 2, 3,….N-1, N, ¿cuántos

dígitos hemos escrito en el sistema decimal de numeración?

Esta cuestión se puede expresar de forma inversa mediante un

problema:

Para numerar las páginas de un libro hemos tenido que escribir

702 dígitos ¿Cuántas páginas tiene el libro?

Es interesante estudiar la relación entre un número natural N y los

dígitos empleados en escribir desde 1 hasta N. ¿Cómo se te ocurre
abordar esta cuestión? Damos tres pistas:

(a) Recuento simple

Deberemos contar un dígito por cada número escrito, otro por cada
número a partir de 10, otro a partir de 100, etc. Esto nos daría, para el
número de tres cifras del ejemplo, la expresión

N+N-9+N-99 = 702; 3N=810; N=270

luego el libro tiene 270 páginas.

(b) Truco de Hoja de Cálculo

A partir de la resolución anterior, ¿podríamos construir una función tal

que dado un número N nos devolviera el número de dígitos empleados
10
en la sucesión 1...N?

Aprovechamos un truco. En las hojas de cálculo una igualdad o

desigualdad verdadera posee el valor 1 y la falsa 0. Podríamos
entonces construir esta función:

D(N)=N+(N-9)*(N>9)+(N-99)*(N>99)+(N-999)*(N>999)+(N-
9999)*(N>9999)

que nos devolvería el valor deseado para cada entero positivo menor
que 100000.

Así se puede construir una tabla para esta función en Hoja de Cálculo

(c) Uso en el aula

Esta función definida en Z puede usarse en las clases de Matemáticas,

como ejemplo de

* Función definida entre números enteros

* Definición por intervalos
* Ejemplo de linealidad a trozos

Este tipo de ejemplos ayuda a extender el concepto de función, que a

veces se queda tan solo en funciones reales, continuas y de definición
simple.

11
¿Te atreverías con la definición de la función inversa de esta? Es
evidente que su dominio no contendría a todos los números naturales.

Lo anterior puede sugerir otras cuestiones similares. Por ejemplo ésta:

Si escribimos en sistema de numeración binario todos los números

naturales comprendidos entre 1 y 2 n, ¿cuántos “unos” hemos escrito?
Se debe expresar mediante una expresión dependiente de n.

La siguiente imagen puede sugerirte la solución:

NÚMERO MÁS L A SUM A D E SU S CI F RAS

Esta propuesta parte de un problema incluido en el libro “Concurso

intercentros de Matemáticas”, de Joaquín Hernández y Juan Jesús
Donaire.

“Encuentra razonadamente un número positivo “n” tal que la

suma de “n” y la suma de sus cifras, resulte ser 379”

No es muy difícil encontrar la solución, 365, aunque el razonamiento

debe ser cuidadoso.

12
¿Sería posible crear un algoritmo que resolviera el problema para
cualquier otro número de tres cifras distinto del 379? Por ejemplo, para
832 la solución sería n=821.

El problema radica en que para algunos datos, como 717, existen dos
soluciones: 696 y 705, ya que 696+6+9+6 = 717 y 705+7+0+5=717.
Para más complicación, existen datos que no producen ninguna
solución, como 222.

Proponemos algunos estudios sobre este problema:

(1) Encontrar un algoritmo que resuelva la cuestión para números de

tres cifras (quizás deba tener dos ramas). La siguiente imagen recoge
uno de ellos.

(2) ¿Existirá alguna caracterización para aquellos números que

admitan, como 717 o 218, dos soluciones?

(3) ¿Se podrá encontrar, igualmente, alguna condición que cumplan

los números que no producen soluciones, como 198 o 266?

Puede ayudar el estudio de las igualdades del tipo 101X+11Y+2Z = N,

y también la construcción de tablas con hoja de cálculo como la que
sigue.

13
Número Cifras Suma

190 1 9 0 200
191 1 9 1 202
192 1 9 2 204
193 1 9 3 206
194 1 9 4 208
195 1 9 5 210
196 1 9 6 212
197 1 9 7 214
198 1 9 8 216
199 1 9 9 218
200 2 0 0 202

14
I GU ALD ADES C UR IOS AS

CASI NARCI SI ST AS

Esta propuesta está basada en un ejemplo incluido por Eugenio

Manuel Fernández Aguilar en su blog Ciencia del Siglo XXI - Mirando
con la mente.

882+332=8833

¿Existirán otros números de cuatro cifras con esta propiedad?

Observa que es una situación parecida a la de los números narcisistas,

tales como 153 = 13+53+33 ó 1033 = 81+80+83+83

Si deseas averiguarlo, implementa este código, que te sirve para Excel

y para [Link] Calc

(Entre paréntesis los comentarios)

Sub busqueda
dim i,j,k,l (Cifras del número)
dim a,b,c (a y b están formados por dos cifras)

for i=0 to 9
for j=0 to 9
for k=0 to 9
for l=0 to 9

a=10*i+j (Formamos un número con las dos primeras cifras)

b=10*k+l (Formamos otro con las dos últimas)
c=100*a+b (Formamos el número total)

if a^2+b^2=c then (Si se cumple, tenemos la solución)

msgbox(a) (Se comunican las soluciones)

msgbox(b)
msgbox(c)
end if

next l
next k
next j
next i

End Sub

15
La respuesta es que existe otro número de cuatro cifras con la misma
propiedad ¿Cuál?

Otras propuestas

(a) Ya que tienes el código adecuado para resolver la cuestión, si

efectúas en él algunos cambios puedes encontrar números narcisistas
de cuatro cifras. Uno es 1634=1 4+64+34+44. ¿Cuáles son los otros dos?

(b) Con otro pequeño cambio puedes encontrar números de cifras

abcd que cumplen que abcd=a3+(bc)3+d3. Hay dos

(c) Con otro cambio más puedes encontrar números de cifras abcd
que cumplen que abcd = (cd)2-(ab)2. Sólo existe una solución.

CUADRADO DEL SI MÉT RI CO

Claudi Alsina, en su libro “Vitaminas matemáticas”, señala como una

propiedad del número 12 la siguiente: 122 = 144 y 212 = 441, es decir,
que el cuadrado de su número simétrico en cifras coincide con el
simétrico de su cuadrado.

Esta propiedad la poseen otras parejas de números, en concreto hay,

si la hoja de cálculo no falla, las siguientes:

Dos parejas de dos cifras: 12 y 21, 13 y 31

Cinco parejas de tres cifras, desde 102 con 201 hasta 311 y 113

Dieciocho de cuatro cifras, desde 1002-2001 hasta 3111-1113

Cuarenta y una parejas de cinco cifras…

(a) Una cuestión sencilla: ¿Qué cifras no pueden figurar entre las
componentes de esos números? ¿Cuál es la causa?

(b) Otra algo más compleja: De las cifras que pueden figurar, ¿qué
combinaciones de ellas habría que desechar?

16
(c) Y más difícil, porque hay que contar bastante: ¿Por qué aparecen
estos números de parejas?: 2 de dos cifras, 5 de tres cifras, 18 de
cuatro y 41 de cinco…

SUMA PANDI G I T AL

En la entrada del blog “Espejo lúdico” ([Link]

[Link]/) de fecha 18 de agosto de 2008, se presentó la
siguiente propuesta:

APRO VEC HANDO L AS CI F RA S

Buscar números tales que entre su cuadrado y su cubo se utilicen

todas las cifras (del 0 al 9) y una sola vez cada una.

Podríamos darle la vuelta a esta propuesta, y en lugar de aconsejar

que no se use la hoja de cálculo, que era lo recomendado, promover
su uso, y de manera más fuerte, exigiendo que sea la propia hoja, sin
ayuda nuestra, quien encuentre la solución. Evidentemente, en ese
caso el objetivo es algorítmico, y no los razonamientos matemáticos
que pedía el Espejo Lúdico.

¿Te atreves a crear una "trampa automática" en la que caigan los

números que cumplan la condición exigida?

Para conseguirlo puedes plantear las siguientes operaciones de hoja

de cálculo

(1) El cuadrado y el cubo del número a probar se descomponen en

17
cifras, una por celda (zona verde de la imagen). Es el primer problema
a resolver.

(2) Se construye una tabla con las cifras del 0 al 9 y se cuenta el

número de veces que cada una aparece tanto en el desarrollo en cifras
del cuadrado como del cubo. (zona amarilla)

(3) La celda de "Se cumple" o "No se cumple" examina los contadores,

y si todos presentan el valor 1 (¿cómo se averigua eso en una sola
operación?) da por válido el número.

(4) Se va probando, de forma manual o automática (mediante un bucle

con ayuda de macros) en un rango de búsqueda, y se espera a que
aparezca el número probado como válido. Esto ocurre muy pronto.

¿Te atreves a construir algo similar?

L AS PRI M ERAS, DO BL E DE L AS SE G UNDAS

El siguiente problema consiste en encontrar los dos únicos números de

seis cifras que son iguales a un cuadrado menos uno, y en los que la
última mitad (los tres últimos dígitos tomados como un número de tres
cifras) es el doble que la primera.(Propuesta del blog “Números”)

Es decir que se cumple: abcdef = n 2-1 y abc = 2 x def

Los dígitos abcdef no tienen que ser todos diferentes.

Si las tres primeras son el doble de las tres segundas las soluciones
son: 190095, 446223 y 806403 y si es al revés: 112224 y 444888

Lo bueno de este problema es que se puede abordar con distintas

técnicas:

18
Algebraica

Es la que ofrece el autor del blog, que en esencia consiste en lo

siguiente:

Podemos llamar x al número formado por las tres cifras inferiores, con
lo que el resto del número sería 2000x (o bien al revés) y pplanteamos
que 2001x = n(n+2), ya que todo cuadrado menos 1 equivale al
producto de dos enteros cuya diferencia es 2 (en el caso simétrico
sería 1002x=n(n+2) y deberemos intentar descomponer 2001 en
factores y ver cuáles de ellos se pueden completar a un producto de
dos factores que se diferencien en dos unidades. Los factores de 2001
son: 2001*1; 667*3; 87*23 y 29*69 y se deberán completar
multiplicando por un número de tres cifras hasta conseguir el producto
del tipo n(n+2)

(Ver Soluciones)

Hoja de cálculo sin macros

Se forma una columna con todos los múltiplos de 2001 (o de 1002) que
tengan seis cifras (supongamos que es la D) y en la columna paralela
siguiente (la E) se inserta una fórmula similar a la
siguiente:

=SI(D6+1=ENTERO(RAIZ(D6+1))^2;"SI";"")

Que viene a expresar que si D6+1 es cuadrado

perfecto (igual al cuadrado de la parte entera de
la raíz) se escribirá un SI, y en caso contrario se
dejará en blanco. Al rellenar esa fórmula
observaremos que aparece un SI en las soluciones 190095, 446223 y
806403. Cambia a 1002 y obtendrás las otras.

Hoja de cálculo con macro

Si se intenta mediante Basic, el código de macro adecuado sería:

Sub buscar
v=7

for i=1 to 999

a=1002*i (o bien 2001)
19
if a+1=int(sqr(a+1))^2 then (se prueba si es cuadrado perfecto)
v=v+1
[Link](0).GetCellByPosition(4,v).value=a (En
OpenOffice Calc)
[Link](1).Cells(v,5).Value = " " (En Excel)
end if
next i
End sub

Sería esta buena ocasión para iniciarte en la programación de macros.

Puedes consultar la Guía correspondiente.

([Link]

Variantes

¿Existen soluciones similares con cuatro cifras? Sí: 4623 si las dos
primeras son el doble de las segundas, y 1224 y 4488 en el caso
contrario. Intenta demostrarlo o encontrarlas con la hoja de cálculo.

¿Se podrían estudiar cuestiones similares con n2-4 o con n2 – 9? ¿Y

con n2+1?

(Ver Soluciones)

NÚMERO S AUT O MÓ RF I CO S

Los números de la primera columna de la siguiente tabla son

automórficos. Si los estudias adivinarás pronto qué propiedad tienen
para recibir este nombre.

1 1
5 25
6 36
25 625
76 5776
376 141376
625 390625
9376 87909376
90625 8212890625
109376 11963109376

20
Efectivamente, interviene su cuadrado en la propiedad que está
patente en la tabla.

¿Cómo podríamos encontrarlos con una hoja de cálculo? Para

construir la tabla que se incluye se han usado macros, pero se puede
prescindir de ellas. Puedes crear una tabla de números consecutivos y
después aplicarles una condición. ¿Cuál?

Esta tabla es complementaria de la anterior. ¿Qué relación tiene con

ella?

1 0 0
5 4 20
6 5 30
25 24 600
76 75 5700
376 375 141000
625 624 390000
9376 9375 87900000
90625 90624 8212800000
109376 109375 11963000000

En ella tienes contenido el procedimiento de búsqueda.

Notas
Después de publicar esta entrada, se recabó más información sobre
este tipo de números, mucha de ella interesante, que se añade en
forma de notas. En la sección de Soluciones se completan algunos
razonamientos.

(1) Salvo los casos triviales de 0 y 1, todos los números automórficos

terminan en 6 ó 5.

Basta ver que si a es automórfico, a2–a = (a-1)a es múltiplo de 10, y

que esto se reduce a cuatro casos. Dos de ellos son triviales, el 0 y el
1, y los otros terminan en 5 o en 6.

(2) Si un número es automórfico, no sólo coincide en sus cifras con las

últimas de su cuadrado, sino también con las del cubo y todas las
21
demás potencias. Por ejemplo, las potencias de 76 terminan todas en
76

76
5776
438976
33362176
2535525376
192699928576
14645194571776

(3) Si a un número automórfico de varias cifras se le suprime la

primera, sigue siendo automórfico.

Esto nos lleva a que un automórfico cualquiera, como 109376, contiene

en sí todos los automórficos que comparten la última cifra con él: 6, 76, 376,
9376 y 109376.

(4) En cada nivel de cifras sólo pueden existir dos números

automórficos, uno terminado en 5 y otro en 6.

En efecto, de una cifra existen dos, 5 y 6, de dos cifras otros dos, 25 y

76, y se puede demostrar que si a cada uno de ellos se le añade una
determinada cifra, se convierten en automórficos de una cifra más (ver
Soluciones).

Si se suma los dos automórficos del mismo número de cifras, resulta

siempre un número de la forma 10 k+1:

5+6=11
25+76=101
625+376=1001

(5) Si m es automórfico de k cifras, entonces 3m 2 – 2m3 (mod 102k)

también es automórfico. Esta propiedad permite generar otro
automórfico con doble de cifras.

Así, a partir del automórfico 109376 se genera

22
3*1093762- 2*1093763 (mod 1012) = -2616918212890624 (mod 10 12) =
7383081787109376 (mod 10 12) = 081787109376

No se incluye demostración a causa de su longitud.

EL F Ó SI L DE UN NÚMERO

(Problema propuesto en la Fase provincial de Alicante de la XIX

Olimpiada Matemática, 2008)

Dado un número natural N, se multiplican todas sus cifras. Se repite el

proceso con el resultado obtenido, hasta obtener un número de una
cifra únicamente; a ese número se le llama el fósil de N. Por ejemplo,
el fósil de 327 es 8. Hallar el mayor número natural, con todas sus
cifras distintas, cuyo fósil sea impar.

Nosotros le daremos unas vueltas a la idea de “fósil” de un número.

(1) ¿Tienen fósil todos los números naturales?

Te lo puedes plantar en dos pasos:

23
 (a) El algoritmo de multiplicar todas las cifras produce una sucesión
estrictamente decreciente y llega a términos de una cifra.

(b) Sólo los números de una cifra son invariantes en el proceso.

(2) Construye un algoritmo de hoja de cálculo tal que dado un número

natural, encuentre su fósil. Puedes restringirlo sólo a números de tres o
cuatro cifras, pero ten en cuenta que si disminuye el número de cifras
no pueden aparecer ceros, que arruinarían el cálculo. En el algoritmo
de la imagen, cuando disminuye el número de cifras aparece la unidad,
para no desvirtuar el producto.

(3) ¿Obtendríamos otro tipo de fósil si sumáramos las cifras en lugar de

multiplicarlas?

(4) Se pueden aplicar estas ideas al aula si se restringe el estudio a tres

cifras, por ejemplo. Se podrían formar grupos e intentar que cada uno,
con calculadora u hoja de cálculo lograra todos los fósiles posibles
entre 0 y 9, y después se discutieran algunos casos:

¿Cuándo el fósil resulta ser cero? ¿Qué crees que hay más, fósiles
pares o impares? ¿Por qué siempre se desemboca en una cifra?...

24
C ON PR IMOS Y MÚLT IPLOS

PRI MO S CO N CI F RAS CO N SECUT I V AS

(Sobre una propuesta del blog “Números”-

[Link]

a) Encontrar todos los primos cuyos dígitos son consecutivos y

están ordenados de menor a mayor (yo encontré cinco), y de
mayor a menor.

Esta cuestión es muy interesante para la construcción de un algoritmo

sobre ella. La idea es comenzar con la cifra de las unidades, que sólo
puede ser 1, 3, 7 ó 9, e ir adosando a esa cifra por la izquierda cifras
decrecientes en una unidad, analizando en cada paso si es primo o no.

Por ejemplo:

7 es primo. Le adoso un 6.

67 es primo. Le adoso un 5.

567 no es primo. Lo salto, pero le adoso un 4.

4567 es primo. Añado un 3.

Así se van añadiendo hasta llegar a n cifras.

Se puede plantear en el Basic de la hoja de cálculo, con estos

resultados:

25
Primos con cifras ascendentes:
Código

Sub buscaprimos(n) (n es una cifra igual a 1,3,7 ó 9)

dim m,p,n

m=n:p=n (Variables que albergan las cifras)

for i=1 to n (Se recorren n cifras para agotar las posibilidades)

p=p-1 (Cifras crecientes, van disminuyendo a la izquierda)
m=10^i*p+m (Se forma el número con cifras crecientes)
if esprimo(m) then msgbox(m) (Si es primo, se comunica)
next i

End Sub

Si se dan a n los valores 1,3,7 y 9, resultan las soluciones 23, 67, 89,
4567 y 23456789

Para los descendentes sólo hay que corregir un detalle:

Sub buscaprimos(n) (n es una cifra igual a 1...9 )

dim m,p,n

m=n:p=n

for i=1 to n
p=p-1 (Cifras decrecientes)
m=10*m+p (Se forma el número con cifras decrecientes)

if esprimo(m) then msgbox(m) (Si es primo, se comunica)

next i

End Sub

Con este otro código, dando valores a n entre 1 y 9, resultan dos

soluciones: 43 y 76543.

Esta cuestión no es importante, pero resalta la simplicidad que puede

tener un algoritmo que después resulta bastante potente.

26
PRI MO S R EVER SI BL ES ( PRI MO - O MI RP)

Es muy popular la definición de los pares de números primo-omirp, o

primos reversibles, que son aquellos en los que uno se forma
invirtiendo las cifras del otro (en base de numeración decimal) y que
ambos son primos, como los pares 199 y 991, 7589 y 9857. Se suelen
excluir los capicúas.

No vamos a insistir en el concepto, que incluso se recoge en la

Wikipedia, sino en la posibilidad de encontrarlos con Hoja de Cálculo.

Para ello necesitamos las dos funciones del Apéndice,

INVERTIR_CIFRAS y ESCAPICUA. Además, deberemos contar con la
función ESPRIMO, uno de cuyos posibles códigos figura también en el
Apéndice.

Se pueden encontrar así (excluyendo capicúas) 4 parejas de dos cifras

(13 – 31, 17 – 71, 37 – 73, 79 – 97), 14 parejas de tres cifras, desde
107-701 hasta 991-199, y 102 de cuatro cifras. Se pueden ordenar bien
los cálculos usando las mencionadas funciones.

Aquellas personas no interesadas en la programación de macros

pueden descargar el archivo de texto [Link], que contiene una lista
27
con los números “omirp” inferiores a 200.000, de la dirección
[Link] Una vez descargado es fácil
trasladar los datos a una hoja de cálculo. Sobre esta lista doble se
pueden plantear cuestiones bastante atractivas:

Cuestiones

(a) Muchos números “omirp” restados con su pareja producen

cuadrados perfectos. ¿Qué se puede escribir en las celdas de la
derecha de las listas para que se detecte esta propiedad?

Llama la atención la repetición de los resultados 30, 60 y 78 en las

raíces cuadradas, así como 738, todos múltiplos de 6. ¿Será siempre
así? ¿Será siempre divisible entre 6 la raíz cuadrada de la diferencia
entre dos “omirp”?

Como curiosidad, si siguiéramos escribiendo números hasta un millón,

una raíz cuadrada que aparece varias veces es 666.

746.203 – 302.647 = 443.556 = 666 2

767.323 – 323.767 = 443.556 = 666 2

785.143 – 341.587 = 443.556 = 666 2

797.353 – 353.797 = 443.556 = 666 2

(b) Si hemos buscado diferencias que sean cuadrados perfectos, nada

nos impide intentarlo con triangulares. Un criterio fácilmente
programable en hoja de cálculo es que un número N es triangular si
8n+1 es cuadrado perfecto. ¿Existirán diferencias triangulares?

28
(c) Un poco más difícil: ¿Existirán diferencias que sean cubos
perfectos?

Notas

(1) Una lista más amplia de los “omirps” que producen cuadrados
perfectos al restar es
37, 73, 1237, 3019, 7321, 9103, 104801, 105601, 106501, 108401, 111211, 112111,
120121, 121021, 137831, 138731, 144541, 145441, 150151, 151051, 161561, 165161,
167861, 168761, 171271, 172171, 180181, 181081, 185681, 186581, 189337, 194891,
198491, 302647, 305603, 306503… ([Link]

La distribución de apariciones de cuadrados perfectos al restar las

parejas primo-omipr (contando una diferencia por pareja) inferiores a
200000 es la siguiente:

Diferencia Frecuencia
36 1
900 10
3600 5
544644 1

(2) Ya que hablamos de parejas, podríamos definir como “hijo” al

mayor número primo que divide a la diferencia primo-omirp. Por
ejemplo, para números menores que 10000 el hijo más crecidito es 449
y el más pequeño 3. Como las diferencias son pares, eso significa que
nunca serán potencias de 2.

Si ampliamos la búsqueda hasta 200000 nos aparecen dos hijos algo

monstruosos:

996001 y 100699 tienen como hijo mayor a 49739

997001 y 100799 tienen como hijo mayor a 49789

29
(3) Otros omirps al restarse producen cubos perfectos. Los primeros
son
1523, 3251, 7529, 9257, 154747, 165857, 171467, 174767, 312509, 322519, 373669,
747451, 758561, 764171, 767471, 905213, 915223, 966373, 000033, 1020233,
1077733, 1078733, 1083833, 1099933, 1165643, 1173743, 1175743
([Link]

MÚL T I PL O S DE 1 1

En un tomo de la colección “La tortuga de Aquiles” hemos encontrado

el siguiente problema:

Determinar todos los números N de tres cifras que tengan la propiedad

de ser divisibles por 11 y que N/11 sea igual a la suma de los
cuadrados de los dígitos de N

Es un problema complicado, por lo que, con un poco de humor,

recorreremos varias opciones de resolución según el ánimo que nos
dejen los primeros intentos:

Método directo: Sea N=100a+10b+c, luego se cumplirá que

N=100a+10b+c = 11(a2+b2+c2)

(1) Después de simplificar esto un poco, intentar eliminar alguna

variable y seguir un desarrollo tremebundo de tipo algebraico, se
desemboca en dos discriminantes de ecuaciones de segundo grado
que han de ser cuadrados perfectos, y ¡oh maravilla!, descubrimos las
dos soluciones. No es recomendable.

(2) Con ayudita: El mismo método anterior desemboca mejor en esos

discriminantes si consideramos que los múltiplos de 11 de tres cifras
sólo pueden tener estas dos expresiones:

a(a+b)b si a+b<10 o (a+1)(a+b-10)b si a+b>10 (Nos referimos a

expresión decimal y no a un producto)

30
Seguimos un método similar al anterior, pero ya tenemos eliminada
una variable. Se desemboca básicamente en dos expresiones, cada
una con una sola solución.

(3) Sin Álgebra: Si lo anterior nos asusta, podemos emprender una

búsqueda (buena para un cálculo mental mientras paseamos. Así lo
resolvimos hace días).

¿En qué terminarán esos números? Expresemos como 10a+b el

número N/11 y sólo buscaremos entre 10 y 90.

Si b=9,8,7,6 es fácil ver que no hay solución, porque 10a+b ha de ser

mayor que el cuadrado de b, y si a+b<10, incluso mayor que el doble
de su cuadrado. Con pocas pruebas se desecha esta posibilidad.

Si b=5,4,3 la condición anterior se cumple con más facilidad, por lo

que hay que ir con más cuidado. Serán mejores candidatos aquellos en
los que a+b>=10 ¿Por qué? Quizás encontremos alguno terminado en
5,4 ó 3.

Si b=2,1,0, caliente, caliente…

(4) Con hoja de cálculo: Una búsqueda sistemática se puede organizar

creando una columna con los números que van desde 10 hasta 90,
multiplicándolos por 11 en otra columna paralela. Después se
descomponen estos últimos en sus tres cifras ¿cómo? y finalmente se
calcula la suma de sus cuadrados y se comparan con la primera
columna.

Apúntate a un método o dos y encuentra las dos soluciones que

existen.

DAMO S VUEL T AS A P RI MO S Y AL 1 8

Hace unos días Honorio, un seguidor de este blog, nos envió la

siguiente conjetura: “Entre dos números primos consecutivos cuyos
dígitos sumen lo mismo, como mínimo hay una diferencia de 18 entre
ambos”.

31
Me causó sorpresa y aunque el tema de primos y cifras no es de los
que más me entusiasman me puse a pensar en ella. Pronto descubrí
que esta propiedad no la tienen por ser primos, sino por ser impares (el
2 no entra en la conjetura porque no coincide su suma con el
consecutivo). Lo podemos demostrar:

La diferencia entre dos números impares distintos que presenten la

misma suma de cifras es siempre un múltiplo (no nulo) de 18.

En efecto, si tienen la misma suma de cifras ambos presentarán el

mismo resto módulo 9 (recuerda el criterio de divisibilidad entre 9),
luego su diferencia es múltiplo de 9. Pero como ambos son impares, su
diferencia es par, luego también es múltiplo de 18, no nulo, porque
ambos números son distintos. Luego el valor mínimo de la diferencia
es 18, y todas las demás, múltiplos de dicho número.

Esta propiedad abre un abanico de posibilidades: los primos pueden

ser consecutivos o no. La diferencia suele ser 18 pero puede ser
mayor. Podíamos dar algunas vueltecitas al tema:

V1) Primos consecutivos con la misma suma de cifras y diferencia

18

Si disponemos de las funciones PRIMPROX (próximo primo),

ESPRIMO y SUMACIFRAS, ya tenemos las condiciones de búsqueda.
Lo hemos realizado con el resultado de
523, 1069, 1259, 1759, 1913, 2503, 3803, 4159, 4373, 4423, 4463, 4603, 4703, 4733,
5059, 5209. 6229. 6529, 6619, 7159, 7433, 7459, 8191, 9109, 9749, 9949, 10691,
10753, 12619, 12763, 12923, 13763, 14033, 14303, 14369, 15859, 15973

523 (Sólo se escribe el primer número primo de cada par)

1069
1259
1759 Con nuestro Buscador de naturales puedes reproducirla
1913
2503 planteando las condiciones
3803
4159
4373
ES PRIMO(N)
4423 ES SUMACIF(N)=SUMACIF(PRIMPROX(N))
4463
4603
4703
4733
ES PRIMPROX(N)=18+N
5059
5209

32
Se exige que N sea primo, que tenga la misma suma de cifras que el
siguiente primo y que su diferencia sea 18. Si deseas ver el par
completo añade EVALUAR PRIMPROX(N)

Siempre que encontramos una secuencia la comprobamos en OEIS

para ver si está publicada, y en este caso no lo está, por lo que la
hemos propuesto con el número A209875 [Link] Hoy
la nombraré como V1

V2) Primos con la misma suma de cifras que se diferencian en 18

Parece la misma cuestión, pero es que no exigimos que sean

consecutivos. Por ejemplo, el 2 y el 11 presentan la misma suma y se
diferencian en 9. Para buscarlos bastará ver que p sea primo y p+18
también, y que tengan la misma suma de cifras. Como las condiciones
son menos restrictivas, es normal que aparezcan muchos más.

El resultado es este:
5, 13, 19, 23, 29, 43, 53, 79, 109, 113, 139, 149, 163, 173, 179, 223, 233, 239, 263,
313, 349, 379, 439, 443, 449, 491, 503, 523, 569, 613, 643, 659,…

Se puede reproducir con el Buscador con las siguientes condiciones:

PRIMO
ES PRIMO(N+18)
ES SUMACIF(N)=SUMACIF(N+18)

5 En la imagen tienes el resultado. También aquí puedes

13
19 ver el par completo con EVALUAR N+18
23
29
43 Esta sucesión incluye a la V1. No estaba publicada en
53
79 OEIS, por lo que la hemos incluido con el número
109
113 A209663 [Link]
139
149
163
173
179
223 Si la nombramos como V2, ya tenemos que V1 V2.
233
239
263

33
V3) Primos consecutivos con la misma suma de cifras

En este caso las diferencias entre ellos serán múltiplos de 18.

El resultado es muy parecido al de V1 y está publicado en OEIS hace

tiempo
523, 1069, 1259, 1759, 1913, 2503, 3803, 4159, 4373, 4423, 4463, 4603, 4703, 4733,
5059, 5209, 6229, 6529, 6619, 7159, 7433, 7459, 8191, 9109, 9749, 9949, 10691,
10753, 12619, 12763, 12923, 13763, 14033, 14107, 14303,…
[Link]

Puedes reproducirla en el Buscador de Naturales con

PRIMO
ES SUMACIF(N)=SUMACIF(PRIMPROX(N))

El primer par con diferencia 36 es (14107,14143). El primero con

diferencia 54 es (35617, 35671) y el primero con 72 (31397, 31469)

Es claro que V1 es un subconjunto de V3, porque 14107 o 35617

pertenecen a V3 y no a V1

Estos pares de consecutivos se pueden ampliar a tripletes: tres

números primos consecutivos con la misma suma de dígitos

Los primeros que hemos encontrado son:

22193 22229 22247

25373 25391 25409
69539 69557 69593
107509 107563 107581
111373 111409 111427
167917 167953 167971
200807 200843 200861
202291 202309 202327
208591 208609 208627
217253 217271 217307
221873 221891 221909
34
236573 236609 236627
238573 238591 238627
250073 250091 250109
250307 250343 250361
274591 274609 274627
290539 290557 290593

Estos tripletes tampoco figuraban en OEIS. Ya es de prever que los

hemos incorporado (ver A209396)

Me he puesto a buscar conjuntos de primos consecutivos con la misma

suma de cifras. Después de encontrar este me he cansado. Si alguien
quiere seguir…

1442173, 1442191, 1442209, 1442227

(Claudio Meller, en la entrada que enlazamos al final, presenta estos

cuatro, aunque referidos a igual promedio: 8508473, 8508491,
8508509, 8508527. También nos ha indicado dónde se pueden
consultar los primeros elementos de los pares, tripletes y demás
conjuntos de primos consecutivos con la misma suma. Los puedes
encontrar en [Link] Gracias, Claudio)

V4) Otra vuelta más.

Si dos números presentan la misma suma de cifras también coinciden

en el valor de su raíz digital, que es el número entre 0 y 8 que resulta
si sumamos sus cifras, y después volvemos a sumar las cifras de esa
suma y reiteramos hasta obtener un número menor que 9. Es fácil
razonar que ese número es el resto de dividir el número primitivo entre
9.

El inverso no es cierto: si se da la misma raíz digital las sumas de

cifras no han de ser iguales, sino congruentes módulo 9.

Pues bien, si sólo exigimos que dos números primos consecutivos

tengan la misma raíz digital nos resulta otra sucesión más amplia que
la V1 y la V3, que también se ha publicado en OEIS

35
523, 1069, 1259, 1381, 1759, 1913, 2161, 2503, 2861, 3803, 3889, 4159, 4373, 4423,
4463, 4603, 4703, 4733, 5059, 5209, 5483, 6011, 6229, 6451, 6529, 6581, 6619, 7159,
7351, 7393, 7433, 7459, 7621, 7883, 8191, 8761, 9109, 9293, 9551, 9749, 9949,…
([Link]

Aquí se puede razonar también que las diferencias han de ser

múltiplos de 18. Inténtalo.

V5) Aún quedan vueltas que dar, pero lejos de mí producir mareos
irreversibles. Las presento con breves referencias:

V51) Tener la misma suma de dígitos es una condición fuerte, pero es

más exigente pedir que sean los mismos dígitos, aunque en distinto
orden, los que tengan dos primos consecutivos. Puedes verlos en
[Link] y se llaman pares de Ormiston. Los tienes
completos en [Link] También existen tripletes de
Ormiston.

V52) Y otra vuelta: Claudio Meller, de forma casi simultánea a nosotros

ha tratado el tema, pero con promedios
(ver[Link]
[Link])

Bueno, bueno, ya vale de dar vueltas. Si encontráis temas similares los

incorporo como extensión.

L O S MÚL T I PL O S ACUNAN

Toma el número 231. Ve encontrando sus múltiplos hasta que uno de

ellos contenga 231 entre sus cifras, sin ser la primera ni la última.
Tardarás un poco, porque esto no ocurre hasta 231*533=123123.

Elige otro, como el 2011. Necesitarás multiplicarlo por 5973 para

obtener 12011703.

Llamaremos, con un poco de humor, “cuna” a ese múltiplo que acoge

las cifras del número, pero rodeadas de otras que le den calor.

36
¿Cómo encontrar la cuna de cualquier número entero? Aquí
entraríamos en la programación de macros para las hojas de cálculo.
Lo intentamos.

Mediante potencias de 10
Todo número de cinco cifras (por ejemplo) abcde es igual a
a*104+b*103+c*102+d*10+e. Si otro número más pequeño contiene
algunas de esas cifras en el mismo orden, por ejemplo bcd, su
expresión será b*102+c*10+d. Luego deberemos intentar localizar bcd
en abcde mediante esas expresiones.

Necesitaremos las funciones numcifras y cifra

(ver[Link]
[Link])

Un posible código para ver si b es cuna de a sería:

Function escuna(b,a) as boolean Ver si las cifras de b acunan a las de a
Dim m,n,q
Dim e as boolean

m=numcifras(a)
n=numcifras(b)
e=false
for q=n-1 to m+1 step -1 Se procura no leer ni la primera cifra ni la última
c=0
for k=1 to m
c=c+cifra(b,q-k+1)*10^(m-k) Se construye un número con m cifras
next k
if a=c then e=true El número construido es a. Lo hemos encontrado
next q
escuna=e
End function

Mediante la traducción a texto

Si tanto a como b los traducimos a formato de texto (string) es más
fácil localizar dentro de b. Aquí tienes el código:
Function escuna(b,a) as boolean Ver si las cifras de b acunan a las de a
dim s$,t$
dim m,n

s$=str$(b) Convierte b en un texto

n=len(s$)-1
s$=right(s$,n) Le suprime el primer carácter en blanco
37
t$=str$(a)
m=len(t$)-1
t$=right(t$,m) Hace lo mismo con a
r=instr(b$,a$) Averigua si las cifras de a están incluidas en las de b
if r>1 then escuna=true else escuna=false
end function

Una vez que tienes definida la función escuna te bastará con buscar
resultados en tablas de hoja de cálculo o creando bucles en Basic.
Aquí tienes algunos resultados:

Número Factor Cuna

1 210 210
2 60 120
3 44 132
4 35 140
5 30 150
6 27 162

7 25 175

8 23 184

9 22 198

10 110 1100

11 192 2112

12 94 1128

13 87 1131

14 82 1148

15 77 1155

Observa que se necesitan factores grandes. No es tan simple

reproducir las cifras en un múltiplo.

Entre los primeros 100 números el más tardío es el 27, que necesita
ser multiplicado por 471 para que se reproduzcan las cifras 27:
27*471=12717.

38
Entre los primeros 500 es el número 271 el que necesita un factor
mayor: 271*4691=1271261.

El menos exigente es el 91, que sólo con el factor 21 ya reproduce las

cifras 91*21=1911.

Los primeros valores los tienes en [Link]

¿Se podrán acunar todos los enteros? La intuición nos dice que es
probable, ya que el factor se puede hacer tender a infinito y la
probabilidad de que no aparezca la pauta deseada se haría cada vez
más pequeña. Supongamos que buscamos la pauta 253. En un
número de n cifras se podrían leer n-2 pautas consecutivas. La
probabilidad de encontrar 253 en una pauta es de 0,001, luego la de no
encontrarla será de 0,999, y la de no encontrarla en n-2 pautas de
0,999n-2 y esa expresión tiende a cero si n tiende a infinito.

Bueno, no hay que tomarse el párrafo anterior en serio. Sólo es una

idea vaga y aproximada, pero que quizás apunte en buena dirección.

C O MO E N C A S IT A E N N I N G Ú N S IT IO

En el apartado anterior estudiamos la coincidencia de cifras entre un

número y sus múltiplos. Ahora realizaremos un trabajo similar con
funciones.

A partir de un número natural se pueden definir muchas funciones de

variable entera. Sólo algunas de ellas tienen la propiedad de que, en
valores particulares, sus cifras son una subcadena (substring) de las
propias cifras del número elegido expresado en base decimal.

Por ejemplo, si sumamos las cifras de un número en el sistema

decimal el resultado constituye una función de ese número. Pues bien,
en algunos casos, la expresión de la suma de sus cifras está incluida
en el conjunto ordenado de las cifras del número. Es lo que ocurre con
el 2711, cuya suma de cifras, 11, es una subcadena de 2711. Son
muchos los números que tienen esa propiedad. Aquí tienes los primos
que la cumplen:
39
2, 3, 5, 7, 109, 139, 149, 179, 199, 911, 919, 1009, 1063, 1109, 1163, 1181, 1327,
1381, 1409, 1427, 1481, 1609, 1627, 1663, 1709, 1811, 2099, 2137, 2399,
2699, 2711,…

[Link]

y en esta tienes todos los casos, primos o compuestos

[Link]

Puedes comprobar cualquiera de ellos.

Otros presentan una propiedad similar con una función tan sofisticada
como la indicatriz de Euler

([Link]
[Link])
1, 1320, 1640, 1768, 1996, 2640, 3960, 13200, 16400, 19984, 19996, 26400, 39600,
132000, …

Los puedes consultar en [Link] y además leer las

interesantes propiedades que tienen.

En estos números sus cifras son como la gran casa que acoge a una
función concreta. Hay más casos:
15, 25, 125, 1537, 3977, 11371, 38117, 110317, 117197, 123679, 143323, 146137,
179297, 197513, 316619, 390913, 397139, 399797, 485357, 779917, 797191,
990919…

contienen las cifras de su mayor divisor propio

([Link]

[Link] con el radical en los no libres de cuadrados.

Y existen más curiosidades:

[Link] [Link]
[Link]

40
Propiedades presentadas por nosotros

Aportamos algo más: buscaremos funciones que en ciertos números

estén como “en casita” dentro de sus cifras.

Suma de partes alícuotas

Existen números compuestos (en los primos esto carece de interés) en

los que la suma de sus partes alícuotas (divisores de N menores que
N) tienen sus cifras incluidas como cadenas en las suyas propias. Son
estos:
6, 28, 121, 437, 496, 611, 1331, 1397, 8128, 10201, 14641, 27019, 40301, 40991,
41347, 41917, 45743, 47873, 49901, 51101, 67997, 76459, 97637, 99101, 99553,
99779, 120353, 133307, 133961, 134179, 153091, 161051, 165101, 165743, 166171,
182525, 186503, 190987, 193121, 357101, 357307, 359573, 360397, 418153, 464353,
924611…

Los hemos publicado en [Link] Recuerda que sólo

consideramos los compuestos.

Entre ellos están los números perfectos. Razona el porqué. Todos los
demás es claro que son deficientes e impares. Como el
razonamiento es un poco largo, lo dejamos para el final (ver
Complemento abajo) y así, si no te apetece leerlo, no te estorba para
ver los siguientes.

Un código PARI para obtenerlos puede ser

indigit(a,b)={ u=Vec(Str(a));v=Vec(Str(b));indi=0;la=#u;lb=#v;i=1;while(i<=la-
lb+1&&indi==0,d=0;for(x=1,lb,if(v[x]==u[i+x-1],d+=1));indi=(d==lb)
;i+=1);return(indi)}
{ for(i=4,10^7,if(indigit(i,sigma(i,1)-i)&&isprime(i)==0,print(i)))}

Se presentan casos espectaculares, como 161051, cuya suma de

partes alícuotas es 16105, y 182525 con suma 82525

Suma de factores primos con repetición

Hemos experimentado con nuestra querida función SOPFR (logaritmo

entero: [Link]
[Link], suma de factores primos con repetición). Como en el caso de

41
los números primos la propiedad es trivial (¿por qué?), hemos buscado
la propiedad sólo para compuestos y los primeros son estos:
4, 18, 144, 150, 168, 175, 198, 220, 230, 242, 246, 255, 322, 366, 444, 624, 1166,
1243, 1323, 1330, 1331, 1462, 1480, 1530, 1992, 2187, 2230, 2240, 2406, 2436, 2625,
2650, 2673, 2730, 2744, 2808, 2925, 3024, 3125, 3182, 3264, 3286, 3366, 3388, 3420,
3484, 3591…

Un caso notable es el de 1330 y 1331, ambos con el mismo valor de

SOPFR. En efecto, 1330=2*5*7*19, con suma 33 y 1331=11*11*11 con
igual suma.

Una subsucesión de este ejemplo la tienes en [Link]

Suma de factores primos sin repetición

Con la función afín a la anterior SOPF (suma de los factores primos sin
repetir) también existen números que poseen la propiedad

Son estos
25, 32, 54, 98, 125, 126, 128, 140, 196, 230, 243, 246, 255, 256, 315, 322, 348, 366,
392, 512, 520, 576, 625, 810, 828, 896, 1024, 1029, 1060, 1080, 1152, 1166…

Es fácil comprender que tienen menos interés porque en la potencias

de primos resulta más fácil su cumplimiento. Los hemos incorporado a
OEIS: [Link]

Se obtienen con el código PARI

indigit(a,b)={u=Vec(Str(a));v=Vec(Str(b));indi=0;la=#u;lb=#v;i=1;while(i<=la-
lb+1&&indi==0,d=0;for(x=1,lb,if(v[x]==u[i+x-1],d+=1));indi=(d==lb)
;i+=1);return(indi)}
sopf(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) }
{ for(i=2,10^5,if(indigit(i,sopf(i))&&isprime(i)==0,print(i)))}

Destaca el caso de 1243, cuyos divisores primos son 111 y 13, y su

suma sopf(1243)=111+13=124, que sólo se diferencia del número en
un 3.

42
 Función TAU

La función TAU cuenta el número de divisores de un número. También

ella puede ser una subcadena. Lo es en muchos ejemplos, por lo que
es menos interesante

Aquí tienes un listado de los primeros:

2, 14, 23, 29, 34, 46, 63, 68, 74, 76, 78, 88, 94, 116, 126, 127, 128, 134, 138, 141, 142,
143, 145, 146, 164, 168, 180, 182, 184, 186, 189, 194, 196, 211, 214, 216, 223, 227,
229, 233, 236, 238, 239, 241, 247, 248, 249, 251, 254, 257, 258, 261, 263, 268, 269,
271, 274, 277, 281, 282

PARI para obtenerlos:

Indigit(a,b)={u=Vec(Str(a));v=Vec(Str(b));indi=0;la=#u;lb=#v;i=1;while(i<=la-
lb+1&&indi==0,d=0;for(x=1,lb,if(v[x]==u[i+x-1],d+=1));indi=(d==lb)
;i+=1);return(indi)}
{ for(i=1,1000,if(Indigit(i,sigma(i,0)),print(i)))}

Otros casos de menos interés.

Los ejemplos que presentaremos a continuación como simples

curiosidades provienen de listados mayores, en los que abundan los
casos triviales. Con eso se aumentan excesivamente los números y se
llega a hacer aburrido. Hemos preferido presentar los destacados.

Con divisores de cierto tipo

En casi todos los casos aparecen demasiadas soluciones, con muchos

casos triviales que le quitan interés. Destacamos algunos:

 11371 contiene a su mayor divisor impar propio, 137

 Si a 22742 le suprimes las cifras extremas se convierte en su mayor
divisor par propio, 274.
 Igual le ocurre a 31218 con su mayor divisor cuadrado 121. También
36250 se convertiría en 625.
 Si a 11300 le suprimos las cifras extremas se convierte en su suma de
divisores cuadrados: 130=100+25+4+1
 5145 contiene a la suma de sus divisores triangulares:
145=105+21+15+3+1

43
 Por último, 1584 contiene a la suma de sus divisores que pertenecen a
la sucesión de Fibonacci: 158=144+8+3+2+1

Esto hay que tomarlo como un pasatiempo sin mayor interés.

Otros ejemplos

La parte cuadrada de un número o la parte libre

([Link]
[Link])

Con la parte cuadrada tienes dos ejemplos sencillos pero no triviales:

9225 contiene a su parte cuadrada 225, porque 9225=15^2*41, y
10625=5^4*17 contiene a su parte cuadrada 625. Intenta otros
ejemplos.

Con la parte libre podemos destacar como menos triviales 2835, que
contiene a su parte libre 35 (Comprueba que es así) o el número
2772=2^2-3^2*7*11 que contiene a 77=7*11

Bigomega cuenta los divisores primos con repetición. Con bigomega

destacamos estos:

N Bigomega(N) Factorización
11264 11 [2,10][11,1]
11520 11 [2,8][3,2][5,1]
26112 11 [2,9][3,1][17,1]
38912 12 [2,11][19,1]
49152 15 [2,14][3,1]

Se ha adjuntado la factorización (cada primo con su exponente) para

que los compruebes.

Adjuntamos ahora la demostración anunciada:

Complemento

Si un número aloja a una función suya como substring, la relación entre

sus valores está limitada por unas desigualdades fáciles de obtener. Si
ambos números son iguales, en el caso de las partes alícuotas
resultaría un número perfecto. Dejamos ese caso.

44
Si los números no son iguales, sino que la relación de substring es
estricta, las cifras alojadas pueden ser las primeras, como cuando
2187 aloja a su función SOPFR 21, estar en el interior rodeadas por
otras cifras, como 1331 con sopfr(1331)=33, o bien al final, como
ocurre con la función phi: phi(1768)=768.

Veamos los tres casos, en los que llamamos A al número total y B a las
cifras alojadas. Su cociente A/B es el que vamos a acotar.

(a) Al principio

Si las cifras están al principio, A=B*10^k+C, siendo C un número de k

cifras. El cociente pedido sería: A/B=10^k+C/b, luego A/B>=10^k. En el
más desfavorable de los casos A sería más de diez veces mayor que B

(b) En el interior

Entonces A=D*10^m+B*10^n+C, siendo C un número de n cifras. Así

quedaría

A/B=D*10^m/B+10^n+C/B>=10^n

También, en el caso más desfavorable, A/B>=10

(c) Al final

En ese caso A=C*10^h+B, con B<C*10^h, luego A/B ha de ser mayor

que 2. Por ejemplo, el caso más desfavorable con tres cifras sería
1999/999=2,001001

¿Qué sacamos de todo esto? Pues que en el caso de las partes

alícuotas el número ha de ser deficiente (si no es perfecto), pues su
abundancia B/A<1/2. Ahora bien, no puede ser par, porque en estos
casos el mayor divisor propio M de N es N/2, con lo que tendríamos
que la suma de partes alicuotas sería mayor que N/2 y por tanto la
abundancia sería mayor que 1/2 en contra de lo demostrado mediante
cifras:

Los elementos de la sucesión, o son perfectos o son deficientes

impares.

45
MAYO R DI VI SO R PRO PI O CO N L A MI S MA SU MA DE
CI F RAS

A partir de una cuestión simple (que no sencilla) desarrollaremos varias

técnicas y conceptos, ejercicio que nos agrada mucho en este blog.

¿Qué números poseen la misma suma de cifras que su mayor divisor propio?

Uno de ellos es el 12673 que se descompone como 12673=19*23*29, luego

su mayor divisor propio es 23*29=667 y, en efecto la suma de las cifras de
12673 es 19 y la de 667 también es 19.

No hay que irse tan lejos: el mayor divisor de 18 es 9 y también se cumple

esta propiedad. Puede que hayas pensado que la cumplen todos los múltiplos
de 9 y no es así, porque 45 tiene como mayor divisor 15 y ahí no coinciden las
sumas, y tampoco en 63. Sin embargo sí se cumple en 18, 27, 36, 54,…

Esto se complica, porque tienen la propiedad números no múltiplos de 9 y

entre los que sí lo son, unos la cumplen y otros no. Reflexionemos:

Relación entre un número y su mayor divisor propio

Si B es el mayor divisor propio de A se tendrá que A/B será el menor divisor
de A (mayor que 1 pues si no B no sería propio), pero por uno de los más
elegantes teoremas sobre divisores, ese cociente A/B ha de ser primo. Por
tanto

Si B es el mayor divisor propio de A se cumple A=Bp siendo p primo

(igual o menor que B)

Condición de igualdad de sumas

Si dos números presentan la misma suma de cifras es que son congruentes
módulo 9, como bien se sabe en Aritmética Modular (recuerda el criterio de
divisibilidad por 9).

Por tanto tendremos, con la notación anterior que AB (mod 9, es decir BBp
(mod 9 y por tanto Bp-B=B(p-1) deberá ser múltiplo de 9. Ten cuidado, que
esta condición es necesaria pero no suficiente.

Esto nos lleva a tres posibilidades: (a) B es múltiplo de 9 (b) B es múltiplo de 3

pero no de 9 (c) B no es múltiplo de 3

46
(a) Si B es múltiplo de 9, el número primo p sólo puede ser 2 o 3. Si fuera
mayor no se cumpliría, porque entonces B no sería el mayor divisor, ya que el
menor sería 3 y por tanto el mayor sería A/3. Esto divide a los múltiplos de 9
en dos clases:

(a1) Los números en los que un múltiplo de 9 está multiplicado por 2 o 3 (y por
otros), sí pueden cumplir la igualdad de sumas. De hecho son estos:

18, 27, 36, 54, 72, 81, 90, 108, 126, 135, 144, 162, 180, 198, 216, 234, 243,…

No sé si echas de menos alguno. No estará el 288, a pesar de cumplir el

contener el factor 2, pero es que sus cifras suman 18 y las de su mayor divisor
144 sólo 9.

Aquí tienes la trampa: dijimos que la condición era necesaria pero no

suficiente. Por eso hemos usado la palabra “pueden ser”. Lo que sí ocurrirá es
que las sumas sean congruentes módulo 9 (razónalo)

(a2) Aquellos que tienen la forma 9p con p producto de primos mayores que 3.
Estos no tienen que cumplir la igualdad de sumas. Por ejemplo 873=9*97. Sus
cifras suman 18. Su mayor divisor es 873/3=291 y sus cifras suman 12.

(b) Para que B(p-1) sea múltiplo de 9 también puede ocurrir que B sea
múltiplo de 3 pero no de 9, luego el otro factor 3 debe pertenecer a p-1, de
donde se deduce que p tiene la forma de 3k+1 con k>0, luego p será un primo
mayor que 3 y entonces B no será el mayor divisor de A sino A/3. No se
puede dar este caso.

(c) Si B no es múltiplo de 9, lo será p-1. Así que si A no es múltiplo de 9,

tampoco lo será B y si se cumple la igualdad de suma de cifras, ha de ser
divisible entre un número primo del tipo 9k+1 con k>0. Más exigente aún:
como p-1 es par, p será del tipo 18k+1

Efectivamente, a continuación presentamos los números que cumplen la

propiedad que estamos exigiendo y que no son múltiplos de 9. En todos ellos
figura un factor primo del tipo 18k+1. En la factorización el primer número de
cada corchete es el factor primo y el segundo su exponente.

361 [19,2]
551 [19,1][29,1]
703 [19,1][37,1]
1007 [19,1][53,1]
1273 [19,1][67,1]
1691 [19,1][89,1]
47
 1843 [19,1][97,1]
2033 [19,1][107,1]
2071 [19,1][109,1]
2183 [37,1][59,1]
2413 [19,1][127,1]
2603 [19,1][137,1]
2641 [19,1][139,1]
2701 [37,1][73,1]
2831 [19,1][149,1]
2923 [37,1][79,1]
3071 [37,1][83,1]
3173 [19,1][167,1]
3293 [37,1][89,1]
3743 [19,1][197,1]
3781 [19,1][199,1]

No creas que todos son semiprimos. Hay otros que no lo son: 18791 está en
la lista y se descompone como 19*23*43.

Recuerda también que la condición no es suficiente aquí tampoco.

Hemos publicado esta lista en OEIS con la referencia [Link]

El código PARI que engendra esta secuencia lo tienes a continuación, aunque

en este blog la primera búsqueda se efectúa con hoja de cálculo y después se
comprueba con PARI, Wmaxima o Wiris cuando es posible.

(PARI) digsum(n)={local (d,p); d=0; p=n; while(p,d+=p%10;p=floor(p/10));

return(d)}
largdiv(n)=if(n==1, 1, n/factor(n)[1, 1]) \\ Charles R Greathouse IV, Jun 15 2011
{ k=0; for (n=2, 10^5, if(digsum(n)==digsum(largdiv(n))&&n%9>0,
k=k+1;write("[Link]",k,", ",n))); }

Sí, esta era una cuestión menor, pero nos ha divertido estudiarla. Se puede
aprender mucho con este tipo de planteamientos.

48
PANDI G I T AL ES, CRO MO S Y BE NF O RD

Hace unas semanas conocí esta conjetura:

El número 168 es el mayor N que cumple que la potencia 2^N no contiene

todas las cifras del 0 al 9

[Link]
1/#comment-316640

Es decir, a partir de 2^169 todas las potencias de 2 son pandigitales en

sentido amplio, pues contienen todas las cifras, pero repetidas (usualmente se
exige que los pandigitales presenten cada cifra una sola vez).

2^168 = 374144419156711147060143317175368453031918731001856

(le falta la cifra 2)

2^169 = 748288838313422294120286634350736906063837462003712

2^170 = 1496577676626844588240573268701473812127674924007424

2^171 = 2993155353253689176481146537402947624255349848014848

2^172 = 5986310706507378352962293074805895248510699696029696

2^173 = 11972621413014756705924586149611790497021399392059392

(todos contienen las cifras 0 al 9)

Esta conjetura también está publicada en [Link]

Comienzo de los pandigitales

Nos podíamos preguntar qué ocurre con las demás bases y sus potencias.
Hemos trabajado un poco con la hoja de cálculo y llegado a esta tabla, en la
que figuran las siguientes bases (no múltiplos de 10, que serían casi triviales)
y los exponentes hasta donde llega la carencia de alguna de las cifras en sus
potencias

49
Esta tabla “huele” a inverso de un logaritmo. En efecto, si en lugar del tope en
el que se acaban las potencias no pandigitales (con repetición) nos fijamos en
las cifras de esas potencias llegamos a una cierta uniformidad, especialmente
en las primeras:

Para que una potencia alcance un número de cifras se deberá cumplir de

forma aproximada esta igualdad:

B es la base dada, T el tope no pandigital y C el número de cifras a partir del

cual están representadas todas las posibles. Si tomamos este número de
cifras en un promedio de 50, por ejemplo, nos daría una aproximación del
tope:

El logaritmo es decimal, evidentemente. Si aplicáramos esta fórmula

obtendríamos:

Resulta coherente con los cálculos, luego lo importante es el número de cifras.

El tope es una consecuencia de ellas.

Esto no funciona como algo aleatorio

Este problema, si tuviera una base aleatoria se parecería al de completar una
colección de cromos. Aquí la colección completa sería el conjunto
{0,1,2,3,4,5,6,7,8,9} y los “cromos” se incorporarían uno a uno a la colección.
Cuando aparezcan todos la colección estará completa, pero se habrán
producido repeticiones.

En este blog estudiamos dichas colecciones de sobre en sobre, lo que no es

este caso.

[Link]

50
 [Link]

En otras direcciones puedes consultar una fórmula sencilla para cuando se

incorporan a la colección los cromos de uno en uno, caso más parecido al que
nos ocupa.

[Link]
coleccion-de-cromos/

[Link]

En ellas puedes estudiar una fórmula que te da el total de cromos T que has
de comprar para completar una colección de N

En el caso de diez cifras T=29,29

En nuestro ejemplo hemos necesitado más, unos 50. Es claro. Estamos

comparando como mera diversión dos conceptos diferentes:

 Los cromos aparecen de forma aleatoria y las cifras de las potencias

constituyen un cálculo exacto, determinista.
 En los cromos cada vez que sale uno ya lo tenemos definitivo y aquí en cada
potencia hay que volver a empezar. Esto, en parte, justifica la discrepancia
entre 29 y 50.
 Aquí existe una relación clara de causalidad entre las cifras de 2N y las de
2N+1

Acabamos de afirmar que estudiamos un fenómeno determinista, pero si la

distribución de cifras fuera muy uniforme, sus resultados se acercarían a los
aleatorios. Cuidado: no confundas aleatorio con uniforme. Sólo afirmamos que
los resultados serían más parecidos.

¿Cómo se comportan las potencias respecto a la frecuencia de las distintas

cifras? ¿Qué grado de uniformidad presentan?

Estadísticas de cifras
Según lo anterior, hemos de tener cuidado en considerar como casi uniforme
la aparición de cifras en las potencias de un número. Si así aparecieran,
deberíamos encontrar más similitud entre lo que se espera de un fenómeno
aleatorio y este que nos ocupa.

51
La uniformidad

Para estudiar de forma empírica la distribución de cifras en las potencias

hemos elegido las bases entre 2 y 9, a cada una la hemos elevado a todos los
exponentes comprendidos entre 1 y 50, pues son los cálculos antecedentes
de la región en la que desaparecen los resultados que no son pandigitales.
Para ello nos ha sido útil nuestra calculadora STCALCU para hoja de cálculo
(que presentaremos próximamente). Hemos obtenido este resultado:

(1) Las desviaciones típicas mayores se corresponden con los divisores del
10, 2 y 5. Después baja algo en los números no coprimos con 10: 4, 6 y 8. Por
último, son más homogéneas en los coprimos, 3, 7 y 9. Esto tiene cierto
sentido, pero no seguiremos por ahí.

(2) La distribución por cifras presentan un máximo en la cifra 1 (como en la

Ley de Benford) de 11,2%, muy alejado del 8,7% de la cifra 0. Además, las
cifras impares aparecen más que las pares. Esto lo afirmamos
descriptivamente, pues una prueba chi-cuadrado no da significación.

(3) Casi todas las cifras presentan un máximo en la base igual a ellas. Las
hemos destacado en rojo. Llama la atención el 14% del 5 y del 6. A eso no es
ajeno el que en esas cifras coincidan las terminaciones de sus potencias
sucesivas: 5, 25, 125, 625, … y 6, 36, 216,…Te puedes divertir intentando
analizar otros casos.

Resumiendo, no aparece una uniformidad clara en los resultados, que más

bien parecen sesgados hacia el 1 y los impares. ¿Se mantendrá esta
tendencia para potencias mayores?

Hemos acumulado los resultados desde exponente 1 al 200, para ver cómo
evoluciona la distribución de cifras, llegando a esto:

52
Aquí el panorama cambia algo: se percibe más uniformidad, aunque el 1 es la
cifra que presenta mayor frecuencia. Por tanto debemos pensar que en las
primeras potencias las cifras aparecen con frecuencias más alejadas del 10%
y que eso es lo que produce que se tenga que llegar a unas 50 cifras para
llegar a completar el carácter pandigital

La herencia

Otra pregunta sería pertinente: estas desviaciones de la uniformidad ¿se

mantienen de cierta forma entre unas potencias y las posteriores dentro de
una misma base? Si una cifra presenta una frecuencia en 2N sería interesante
saber cómo se comporta en 2N+1. Pues bien, aquí tampoco se ve relación
clara y significativa entre las frecuencias de un exponente con el siguiente.
Tomamos como ejemplo la base 5 haciendo trampa, porque podía esperase
que la cifra 5 y la 0 se mantuvieran en sus frecuencias al crecer N.

Basta ver los máximos y mínimos para darnos cuenta de lo alejada de la

uniformidad que está la distribución de cifras. Respecto a la herencia, si

53
recorres los porcentajes correspondientes a cada cifra sí se percibe una cierta
constancia en la tendencia. No es importante. Le hemos aplicado la prueba
Chi-cuadrado y no nos da una diferencia significativa respecto a la
homogeneidad máxima.

Así que, por si acaso, no uses potencias para extraer números

psudoaleatorios, que te puedes llevar sorpresas.

Nuestra Ley de Benford

Y ya puestos, ¿cómo se comportan las primeras cifras de cada potencia?

Recuerda que según la Ley de Benford (en la Red tienes muchas referencias
a ella, por ejemplo en
[Link] se podría
esperar un 30% para el 1, un 17% para el 2, 12% para el 3 y así
disminuyendo para el resto, como se ve en la gráfica incluida en la página
recomendada.

Lo intentamos: elevaremos las distintas bases de 2 a 9 (podían ser otras) a

todos los exponentes comprendidos entre 1 y 250 y recogeremos las
estadísticas de la primera cifra.

Son estas:

Esto quiere decir que respecto a la Ley de Benford las potencias se

comportan admirablemente. Hemos comparado nuestras frecuencias con la
fórmula de Benford LOG((d+1)/d) y nos ha resultado:

54
 Potencias Benford

30,2% 30%

17,4% 18%

12,3% 12%

10,0% 10%

7,8% 8%

6,8% 7%

6,0% 6%

5,3% 5%

4,3% 5%

No necesita comentario. El comportamiento de las estadísticas globales viene

dado más por las cifras intermedias que por la primera, que sigue la
distribución esperada. A partir de aquí puedes emprender un estudio del que
sólo hemos esbozado el principio.

55
S OLUCIO NES

O PERACI Ó NE S CO N CI F R AS
Suma pandigital

La solución es el número 69. Puedes utilizar y cambiar a tu gusto el

buscador de números de este tipo alojado en la dirección de Internet
[Link]

Algoritmo derivado de un problema

La solución es 153846, porque 153846*4=615384

Otros resultados relevantes son:

1016949152542372881355932203389830508474576271186440677966 * 6=
6101694915254237288135593220338983050847457627118644067796

163265306122448979591836734693877551020408 * 5 =
816326530612244897959183673469387755102040

Función “Dígitos”

La solución es n*2n-1+1

Número más la suma de sus cifras

El algoritmo que se puede considerar en primer lugar es el de dividir el

número entre 101 para calcular X, hallar el resto y dividirlo entre 11
para hallar Y. Reiterando se calcula Z.

Esta idea tiene algún inconveniente:

(a) Al dividir entre 101, o posteriormente entre 11, se debe considerar

también una unidad menos, pues es muy probable que también sea
solución si Y y Z son cercanos a 9.

56
(b) El resto de dividir el número entre 101 no debe sobrepasar al
número 117, pues entonces Y o Z sobrepasarían el valor de 9.

Según esto, el algoritmo debería ser:

Divido N entre 101 y le llamo X al cociente y H al resto.

Tanto para el valor X como para X-1

Si H>117 paro el cálculo

Si H<=117 lo divido entre 11 con cociente Y y resto M

Tanto para Y como para Y-1

Si M es par y menor que 20, Z será su mitad.

Tienen solución los números del tipo 101X+11Y+Z con Y>0 y Z par,
como es evidente. También tienen solución los terminados en 1, 3, 5 y
7, porque detrayendo 11 de 11Y se pueden convertir en 12, 14, 16 y 18
respectivamente, que al ser pares pueden ofrecer una Z entera.

No tienen solución los números de la forma: 101X+11Y+9 con Y>11,

porque si le añadimos 11 ya el 9 se convertiría en 20, lo que produciría
una Z=10.

Si Y=0, no tienen solución los de la forma 101X+K con K igual a 7 ó 9,

porque aunque detraigamos un 101 de X, o se forma una Y>=10 o una
Z no entera o mayor que 9.

Tienen dos soluciones 101X+K con K par entre 0 y 16, pues si

separamos un 101 de 101X, al ser 101=99Y+2, la nueva última cifra
será par y producirá un Z entero.

57
I G UAL DADES CURI O SAS

Casi narcisistas

Si ejecutamos el código obtendremos otra solución, que es 1233 =

122+332

(a) Los demás números narcisistas de cuatro cifras son

8208=84+24+04+84 y 9474=94+44+74+44.

(b) En el otro grupo propuesto las soluciones son 4160=4 3+163+03 y

4161=43+163+13
(No se cuentan soluciones triviales formadas por 1 y 0)

(c) La solución única es 3468 = 682-342

Cuadrado del simétrico

Si excluimos los capicúas o palindrómicos y los terminados en cero,

para no alterar el número de cifras, las parejas obtenidas son (sólo se
escribe el primer número):

12, 13

102, 103, 112, 113, 122

1002, 1003, 1011, 1012, 1013, 1021, 1022, 1031, 1102, 1103, 1112,
1113, 1121, 1122, 1202, 1212, 2012, 2022

10002, 10003, 10011, 10012, … (hasta 41)

(1) No figurarán en estos números las cifras del 4 al 9, pues en sus

productos por sí mismas producen cifras de arrastre que rompen la
simetría.

(2) Tampoco puede figurar la combinación 23, pues también produce

arrastre 2*3+3*2=12

58
(3) Para contar los casos hay que organizar un diagrama en árbol e
incluir como cifras 0, 1, 2 y 3, evitando la combinación 23.

Las primeras, doble de las segundas

Orientación algebraica

En el caso 2001*1 = 3*23*39 no es posible igualar a n(n+2), pues se

necesitaría un número de cuatro cifras.

Si es 667*3 completamos a 667*3*223 = 667*669 = 446223 = (668-

1)(668+1) = 6682 -1, que constituye la primera solución.

En el caso 87*23 completamos a 87*23*5*19 = 435*437 = 190095 =

(436-1)(436+1) = 4362-1

Por último, 29*69 se completa con 29*31*69*13 = 899*897 = 806403 =

(898+1)(898-1) = 8982-1

El caso contrario, si las segundas son el doble de las primeras, habrá

que analizar el número 1002 = 2*3*167 = 1002*1=167*6 = 501*2 =
334*3

Para abreviar, sólo se darán las soluciones válidas:

167*2*6*56 = 334*336 = 112224 = 3352-1

334*2*3*222= 668*666 = 444888 = 6692-1

Variantes

Para n2-4 existen estas soluciones de seis cifras:

442221 = 6652-4; 760380 = 7822-4; 110220 = 3322-4; 448896 = 6702-4

Para n2-9 existen esta otra solución: 480240 = 693 2-9;

Números automórficos

La propiedad que comparten los números de la primera tabla es que

sus cifras coinciden con las últimas de sus cuadrados (en el sistema de
numeración decimal).

59
La condición que han de cumplir en la hoja de cálculo ha de ser similar
a esta: =SI(RESIDUO(E77-D77;100)=0;”Sí”;"No"), en la que 100 se va
sustituyendo por 1000, 10000, etc. según el número de cifras.

Es decir, que a2-a = a(a-1) = N.10, que es lo que expresa la segunda

tabla.

Para encontrarlos con [Link] Calc puedes usar este código.

Cambia el 10000 por un número mayor si deseas encontrar más.
Sub busquedas
Dim x,fila,i,j as long

fila=10

for i=1 to 10000

x=10^(int(log(i)/log(10))+1)
[Link](0).GetCellByPosition(2,2).value=i
j=(i*i-i)
if int(j/x)=j/x then
fila=fila+1
[Link](0).GetCellByPosition(5,fila).value=i
end if
next i

End Sub

Este código dará como soluciones falsas algunas potencias de 10,

debido al uso de logaritmos decimales.

Notas

(1) En efecto, sólo hay cuatro casos para a-1 y a

(1a) a-1 es múltiplo de 10k y a congruente con 1 módulo 10 k. En este

caso sólo obtendremos el caso trivial 1.
(1b) a es múltiplo de 10k , lo que nos lleva al caso trivial 0.
(1c) a-1 es múltiplo de 5 y a lo es de 2. En este caso a termina en la
cifra 6
(1d) a es múltiplo de 5 y a-1 de 2. El número termina en 5.

60
(2) Si a2≡a (mod 10k), multiplicando por a resulta a 3≡a2 (mod 10k) y
reiterando a≡a2≡a3≡a4 …

(3) Sea a=m.10k-1+n la descomposición del número automórfico de k

cifras en la suma de la primera m y del resto n. Se deberá cumplir que
(m.10k-1+n)2 - (m.10k-1+n) ≡ 0 (mod 10k). Si desarrollamos se tendrá
que la expresión

m2.102(k-1)+n2+2mn10k-1 - m.10k-1-n ≡ 0 (mod 10k) ≡ 0 (mod 10k-1)

lo que nos lleva a que n 2-n ≡ 0 (mod 10k-1) y que por tanto, n es
automórfico con una fila menos.

(4) Si la situación es la inversa de la anterior, hay que encontrar un

valor de m que haga que la expresión anterior sea automórfica para k
cifras, lo que obliga a que 2mn-m+(n2-n)/ 10k-1 termine en cero, es
decir, que se puede despejar m de la ecuación en congruencias

m(2n-1)+(n2-n)/ 10k-1≡0 (mod 10)

En los casos tratados, 2n-1 es congruente con 1 o con 9, lo que

garantiza la solución única.

Por ejemplo, generemos automórficos a partir del 6:

n=6; k=1 Ecuación 11m+3≡0 (mod 10) Solución m=7

n=76; k=2 Ecuación 151m+57≡0 (mod 10) Solución m=3
n=376; k=3 Ecuación 751m+141≡0 (mod 10) Solución m=9
y así se puede seguir hasta la generación de infinitos automórficos.

Para n=5 quedaría este proceso:

n=5; k=1 Ecuación 9m+2≡0 (mod 10) Solución m=2

n=25; k=2 Ecuación 49m+6≡0 (mod 10) Solución m=3
…

El fósil de un número

Solución del problema

Puesto que buscamos un fósil impar, ninguna de las cifras del número
que buscamos puede ser par.

61
Como todas sus cifras deben ser diferentes, para ser lo más grande
posible deberíamos usar todas las cifras impares, es decir, el 1, el 3, el
5, el 7 y el 9. Sin embargo, el producto de todos estos números es el
3*5*7*9 = 945, y tiene una cifra par, de forma que al final tendría un
fósil par.

Por tanto, tenemos que usar cuatro cifras. Quitar la más pequeña, el 1,
no ayuda, pues el producto seguiría siendo el mismo, de forma que
probamos a quitar el 3, así que el producto 5*7*9 = 315, y 3*5 = 15,
proporciona un fósil impar, el 5.

Ya sabemos que debemos usar los números 1, 5, 7 y 9, y para que sea

lo mayor posible pondremos los mayores en las posiciones más
significativas (más a la izquierda). El número buscado es, entonces, el
9751.

(1) Todos los números naturales tienen fósil. Lo veremos en dos pasos:

(a) Todo número natural es mayor que el producto de sus cifras. En

efecto, sea N=a0*10n + a1*10n-1….+an-1*10+an > a0*10n > a0* a1*…* an-1*
an

Por tanto, la operación de multiplicar las cifras produce una sucesión

estrictamente decreciente, pero al tener como cota el 0 y ser finito el
conjunto posible de valores del producto, éste tiene que llegar
necesariamente a un número de una cifra.

(b) Los números de una cifra son invariantes

Es un hecho evidente, por lo que completa el razonamiento.

(2) Un algoritmo para calcular el fósil ha de contener una fase que

descomponga el número en sus cifras

9 8 0 7

9807 807 7 7

1 9 72 72 504

Esta operación se corresponde con las dos primeras filas. En la

primera se calculan los cocientes de la segunda fila entre 1000, 100…y
62
en esta segunda los residuos del número anterior respecto a los
mismos 1000, 100…De esta forma la primera fila contendrá las cifras
del número.

Si el algoritmo está preparado para números de cuatro cifras, por

ejemplo, y escribimos un número de tres, aparecerían ceros al principio
de la primera fila y alterarían el producto final. Para evitarlo en la
tercera fila se van multiplicando las cifras siempre que no sean ceros
iniciales.

Desde 1 hasta el final sólo se van multiplicando cifras si son reales. Por
ejemplo, la primera cifra se multiplica si el número es mayor que 1000,
la segunda si es mayor que 100, y así hasta el final.

Este bloque se repite varias veces, tomando como entrada de cada

una la salida de la otra, hasta llegar a una cifra. En l siguiente imagen
se puede estudiar el proceso:

63
(3) Si sumáramos en lugar de multiplicar, la sucesión que lleva hasta el
fósil también sería estrictamente decreciente:

N=a0*10n + a1*10n-1….+an-1*10+an > a0+ a1+…+an

y por tanto existiría un límite que sería el fósil.

CO N PRI MO S Y MÚL T I PL O S

Primos reversibles (Primo-Omirp)

Cuestión 1: Tiene varias soluciones. Una sencilla para que aparezcan

destacados los cuadrados perfectos es crear una columna nueva con
las raíces cuadradas de las diferencias (supongamos que comienza en
G5) y después crear otra columna nueva con la fórmula
=SI(G5=ENTERO(G5);G5;"")

13 31 4,24264069
17 71 7,34846923
37 73 6 6
79 97 4,24264069
107 701 24,3721152
113 311 14,0712473
149 941 28,1424946
157 751 24,3721152

El que sea múltiplo de 6 se justifica porque la diferencia entre dos

“omirps” es múltiplo de 18, ya que tiene que ser par como diferencia
entre impares, y múltiplo de 9 por ser diferencia entre dos números
simétricos (su diferencia está formada con múltiplos de 9). Si esa
diferencia es además cuadrado perfecto, deberá contener otro factor 2
que formará un 36.

Cuestión 2: La respuesta es afirmativa. Basta someter a las diferencias

a un criterio similar al siguiente:

64
=SI(8*G5+1=ENTERO(RAÍZ(8*G5+1))^2;G5;"")

Entre las diferencias triangulares más frecuentes están 630, 990, 4950
y 16290. La más pequeña es 36.

Cuestión 3: Aunque parecía difícil, también aparecen cubos perfectos,

como en

3251 – 1523 = 1728 = 123

747451 – 154747 = 592704 = 843

758561 – 165857 = 592704 = 843

764171 – 171467 = 592704 = 843

767471 – 174767 = 592704 = 843

Se encuentran con el criterio

=SI(G5=ENTERO(G5^(1/3))^3;G5;"")

Múltiplos de 11

(1) El primer método no es muy recomendable, salvo que se desee

sufrir.

(2) Sea un número de dos cifras a y b entre 10 y 90. Si lo multiplicamos

por 11 obtendremos un número de tres cifras como el considerado en
el problema. Si la suma a+b es menor que 10, ese número estará
compuesto por las cifras a, a+b y b, y en caso contrario por a+1,a+b-10
y b. Esto nos permite plantear dos ecuaciones:

(a) 10a+b = a2+(a+b)2+b2

Desarrollando como una ecuación de segundo grado en “a” (por

ejemplo), y exigiendo que el discriminante esa igual a un cuadrado
perfecto M2 llegamos a la condición

b(3b+8)=25 - M2, que sólo tiene solución si b=0. Esto nos da el valor 5
para M y sustituyendo en la primera ecuación, que a=5. Luego el
numero de dos cifras sería el 50, y el de tres cifras pedido el 550, que
constituye nuestra primera solución.

65
(b) 10a+b = (a+1)2+(a+b-10)2+b2

Procediendo de la misma forma se llega a la condición

b(14-3b)=6+ M2, con solución b=3 M=3 y a=7. El número de dos cifras
sería el 73, y su múltiplo 803, que sería la segunda solución.

Por tanto, las dos soluciones son 803 y 550.

(3) Esta búsqueda está explicada en gran parte y no necesita más

aclaraciones.

(4) Para implementar un algoritmo que resuelva el problema hay que

vencer la dificultad de descomponer el número en tres cifras. Puedes
usar para ello la función CIFRA contenida en el Apéndice..

Una vez resuelto este paso se pueden organizar siete columnas en la

hoja de cálculo. En la primera se situará el número de dos cifras, en la
siguiente su producto por 11, que se descompondrá en tres cifras en
las siguientes columnas. Se calcula la suma de sus cuadrados en la
penúltima cifra y si usa la función condicional SI para que en la última
columna se nos responda si es solución o no.

En la siguiente tabla puedes ver unos resultados:

N N*11 Centenas Decenas Unidades Suma

cuad.
69 759 7 5 9 155
70 770 7 7 0 98
71 781 7 8 1 114
72 792 7 9 2 134
73 803 8 0 3 73 Sí
74 814 8 1 4 81
75 825 8 2 5 93
76 836 8 3 6 109

66
A PÉ ND ICE

FUNCIÓN CIFRA
Argumento: Un número entero m

Valor: Devuelve la cifra de orden n del número

Código en Basic

Public function cifra(m,n)

dim v

v=int(m/10^(n-1)) Divide de forma entera el número entre 10^(n-1)

v=v-10*int(v/10) Le suprime las cifras de la izquierda
cifra=v Se recoge la función Cifra
end function

FUNC IÓN ESCA PI CUA

Argumento: Un número entero

Valor: Devuelve 1 si es capicúa y 0 si no lo es.

Código en Basic

Public function escapicua(n)

dim l%,i%,c%
dim auxi$
auxi = str$(n) (Convierte el numero en texto)
l = len(auxi)
auxi=mid(auxi,2,l-1)
l = l-1
if l<2 then
escapicua = 0
else

67
c=1
for i=1 to int(l/2) (Comprueba cifra a cifra si es igual a su simétrica)
if mid(auxi,i,1)<>mid(auxi,l-i+1,1) then c = 0 (Hay un fallo)
next i
escapicua = c
end if
end function

F U N C I Ó N IN V ER T IR _ C IFR AS

Argumento: Un número entero positivo

Valor: Devuelve el mismo número con sus cifras invertidas

Código en Basic

Public function invertir_cifras(n)

dim l%,i%
dim auxi$,auxi2$,c$
dim m

auxi = str$(n) (Convierte el número en texto)

l = len(auxi)
auxi=mid(auxi,2,l-1) (Le quita un espacio en blanco)
l = l-1 (En todas estas líneas cambia el orden de las cifras)
auxi2=""
for i=1 to l
c=mid(auxi,i,1)
auxi2=c+auxi2
next i
m=val(auxi2) (Convierte de nuevo el texto en número)
invertir_cifras=m
End function

68