Más diversión con los números
El número 12345678987654321=111111111x111111111
se menciona en el PO. Y, lo creas o no, $1234567898765432111$ ¡¡¡es un número primo también!!!
Hay que creerlo. Lo siguiente $maxima$ La salida lo demostrará utilizando El certificado de primalidad de Pratt .
(%i1) m:1234567898765432111;a:29;power_mod(a,m-1,m);factor(m-1); makelist(power_mod(a,(m-1)/p,m),p,map(first,ifactors(m-1))); (%o1) 1234567898765432111 (%o2) 29 (%o3) 1 (%o4) 2*5*11*11223344534231201 (%o5) [1234567898765432110, 972745681039223016, 1223153431857342200, 1089307852892054661]
%i1
es la línea de entrada.
%o1
es la salida que muestra el modul $m$ , el número que comprobamos si es primo.
%o2
es una línea de salida que muestra el número $a$ exponemos.
%o3
es la salida de $a^{m-1} \pmod{m}$ . Debe ser $1$ si $m$ es un primo. Pero podría ser $1$ incluso si $M$ es compuesto.
%o4
es la factorización de $m-1$ y
%o5
son los números $a^{\frac{m-1}{p}}$ para todos los factores primos de $m-1$
Los factores $2$ , $5$ y $11$ son primos pero no es obvio que $11223344534231201$ es un número primo. Así que añadimos una prueba de que $11223344534231201$ también es un primo:
(%i1) m:11223344534231201;a:3;power_mod(a,m-1,m);factor(m-1); makelist(power_mod(a,(m-1)/p,m),p,map(first,ifactors(m-1))); (%o1) 11223344534231201 (%o2) 3 (%o3) 1 (%o4) 2^5*5^2*7*73*101*109^2*137*167 (%o5) [11223344534231200, 7251767246727978,1687389182360412, 5679768249246961, 2181793601204580, 9078160829860754, 8735272021960592, 1320423471360269]
Porque es fácil comprobar que $2,5,7,73,109,137,167$ son primos estamos acabados.
( ambos $a=29$ y $a=3$ puede ser sustituido por el número más grande pero más divertido $1111111$ en los certificados)
Pero, ¿cómo funciona este certificado de primalidad?
Lo siguiente es bien conocido:
- todas las clases de residuos $a$ que son relativamente primordiales para el módulo $m$ constituyen un grupo multiplicativo $\pmod{m}$
- Si $a^{m-1} \pmod{m}$ no es igual a $1$ entonces $m$ no puede ser un primo porque esto contradice El pequeño teorema de Fermat .
- $\text{ord}_{m-1}(a)$ es la menor potencia $e$ tal que $a^{e} \equiv 1 \pmod{m}$ es un divisor de cada $e$ . Especialmente de $\phi(m)$ .
- si $a^{m-1} \equiv 1 \pmod{m}$ entonces $\text{ord}_{m-1}(a)$ es un divisor de $m-1$ . Si es un divisor propio de $m-1$ debe ser un divisor de $\frac{m-1}{p}$ para un primer $p$ dividiendo $m$ . Entonces $a^\frac{m-1}{p} \equiv 1 \pmod{m}$ para tal $p$
- todas las clases de residuos $a,a^2,...,a^{\text{ord}(a)}$ son diferentes entre sí
Así que si $\text{ord}(a)=m-1$ y $\phi(m)=m-1$ y por lo tanto $m$ es un primo.
@ronno anotaría este certificado como
1234567898765432111,29,1234567898765432111-1=2*5*11*11223344534231201
11223344534231201,3,11223344534231201-1=2^5*5^2*7*73*101*109^2*137*167
167,...
137,...
109,...
101,...
0 votos
Gracias; avísame si desarrollas un argumento más preciso.
0 votos
Mi intento: $\frac{123456789098765432111 - 11}{100} = 1234567890987654321$ no sé a dónde ir desde allí... espero que ayude, de alguna manera
0 votos
Relevante: certificado de primacía .
0 votos
Tenga en cuenta que $123456789009876543211$ también es primo, y $1111111111\cdot 1111111111=1234567900987654321$
0 votos
@dtldarek: Sí. No es algo que haya estudiado mucho, aunque soy consciente de (la existencia de) el algoritmo AKS relacionado. Sería interesante que esto llevara a una certificación que un humano pudiera resolver a lápiz en una o dos páginas.
2 votos
@abiessu: Como traté de indicar, me gustaría algo conceptual y no horriblemente de fuerza bruta (porque conociendo a Charles Weibel, apuesto a que tiene algo más conceptual en mente). Aunque no soy un teórico de los números, soy un matemático y estaría dispuesto a tolerar una buena dosis de sofisticación teórica de los números en una respuesta. Me interesaría que la pequeña observación que ha hecho no fuera una pista falsa que ha lanzado por gusto.
0 votos
Tal vez intente algo usando la base $2$ ? Si las tres últimas cifras de un número en base $10$ son $111$ entonces sus tres últimos dígitos en base $2$ también lo son. (Estoy perplejo...)
0 votos
Esto parece un gran candidato para los clásicos certificados del "pequeño teorema de Fermat", probando la primalidad de $p$ mediante el factoring $p-1$ y luego comprobar las potencias de los primos pequeños a los distintos divisores de $p-1$ , $\bmod p$ . El volumen 4 de The Art Of Computer Programming repasa el método con un par de ilustraciones; eso sería lo primero que miraría. (Pero mis copias están en casa y no puedo conseguirlas fácilmente ahora mismo. :/ )
0 votos
@StevenStadnicki Lamentablemente, no tengo el volumen 4, pero tu comentario es en espíritu muy similar a la respuesta de ronno más abajo. Sí, tienes razón -y metodológicamente no veo nada malo en ese tipo de enfoque-, pero me sigo preguntando por qué Weibel cree que es una "pregunta divertida", si no hay una "solución divertida" particularmente por ahí. (¿Pero tal vez no la haya?)
0 votos
@user43208 Corrección menor a mi comentario - me refería al capítulo 4, en el vol. 2; el café aún no ha hecho efecto. :-) Estoy de acuerdo, sin embargo; me parece que hay una versión de este enfoque que es mucho más "bonita" en sus cálculos que la bruta-forja que se ofrece a continuación, pero no estoy seguro de lo que sería de antemano.
0 votos
11111111110^2 - 1000000^2 + 11 = 123456789098765432111 , por si sirve de algo.