Los números de Perrin se pueden expresar como sumas de los tres términos iniciales.
Los primeros catorce números primos de Perrin son
Historia
En 1876, la secuencia y su ecuación fueron mencionadas inicialmente por Édouard Lucas , quien señaló que el índice n divide al término P(n) si n es primo. [5] En 1899, Raoul Perrin preguntó si había algún contraejemplo de esta propiedad. [6] El primer P(n) divisible por el índice compuesto n no fue encontrado hasta 1982 por William Adams y Daniel Shanks . [7] Presentaron una investigación detallada de la secuencia, con una secuela que apareció cuatro años después. [8]
Propiedades
Un microcosmos de Perrin: el algoritmo de tiempo de escape se aplica al mapa de Newton de toda la función de Perrin F(z) alrededor del punto crítico z = 1,225432 con un ancho de ventana gráfica de 0,05320. Las cuencas de atracción que emanan del centro corresponden al número infinito de raíces reales (izquierda) y complejas (derecha) F(z) = 0 .
La secuencia de Perrin también satisface la relación de recurrencia. A partir de esto y de la recurrencia que la define, se pueden crear un número infinito de relaciones adicionales, por ejemplo
La forma polar es con Dado que la fórmula se reduce al primer o segundo término sucesivamente para n grandes positivos o negativos , y los números con subíndices negativos oscilan. Siempre que α se calcule con suficiente precisión, estas fórmulas se pueden utilizar para calcular los números de Perrin para n grande .
Expandir la identidad proporciona la importante regla de duplicación del índice mediante la cual se vinculan las partes directa e inversa de la secuencia.
El índice primo p divide a P(p)
Si la ecuación característica de la secuencia se escribe como entonces los coeficientes se pueden expresar en términos de raíces con las fórmulas de Vieta :
El teorema fundamental de los polinomios simétricos establece que cada polinomio simétrico en las raíces complejas de monic se puede representar como otra función polinómica en los coeficientes enteros de
Reorganice en sumandos monomios simétricos , permutando los exponentes i, j, k:
Sustituya potencias por números primos y raíces complejas por números enteros y calcule representaciones en términos de para todas las funciones polinomiales simétricas. Por ejemplo, is y la suma de potencias se pueden expresar en los coeficientes con el esquema recursivo de Newton . Se deduce que la identidad tiene términos enteros y ambos lados son divisibles por números primos.
Para demostrar que un primo se divide en secuencia inversa, se debe reflejar la ecuación característica . Las raíces son entonces los coeficientes y se aplica el mismo razonamiento.
Prueba de primalidad de Perrin
Consulta 1484. La curiosa proposición de origen chino que es objeto de la consulta 1401 [10] proporcionaría, de ser cierta, un criterio más práctico que el teorema de Wilson para verificar si un número dado m es primo o no; bastaría con calcular los residuos con respecto a m de términos sucesivos de la secuencia de recurrencia u n = 3u n−1 − 2u n−2 con valores iniciales u 0 = −1, u 1 = 0. [11]
He encontrado otra secuencia de recurrencia que parece poseer la misma propiedad; es aquel cuyo término general es v n = v n−2 + v n−3 con valores iniciales v 0 = 3, v 1 = 0, v 2 = 2. Es fácil demostrar que v n es divisible por n , si n es primo; He comprobado, hasta valores bastante altos de n, que en el caso contrario no lo es; pero sería interesante saber si esto es realmente así, especialmente porque la secuencia v n da números que crecen mucho menos rápidamente que la secuencia u n (para n = 17, por ejemplo, se encuentra u n = 131070, v n = 119 ), lo que conduce a cálculos más simples cuando n es un número grande. La misma prueba, aplicable a una de las sucesiones, se aplicará indudablemente a la otra, si la propiedad enunciada es cierta para ambas: sólo es cuestión de descubrirla. [12]
La secuencia de Perrin tiene la propiedad de Fermat : si p es primo , P(p) ≡ P(1) ≡ 0 (mod p) . Sin embargo, lo contrario no es cierto: algún compuesto n aún puede dividir P(n) . Un número con esta propiedad se llama pseudoprimo de Perrin .
Malo y Jarden consideraron la cuestión de la existencia de pseudoprimos de Perrin, [13] pero no se conoció ninguno hasta que Adams y Shanks encontraron el más pequeño, 271441 = 521 2 (el número P(271441) tiene 33150 dígitos decimales). [14] Jon Grantham demostró más tarde que hay infinitos pseudoprimos de Perrin. [15]
Los diecisiete pseudoprimos de Perrin por debajo de 10 9 son
Adams y Shanks notaron que los números primos también satisfacen la congruencia P(−p) ≡ P(−1) ≡ −1 (mod p) . Los compuestos en los que se mantienen ambas propiedades se denominan pseudoprimos de Perrin restringidos. Sólo hay nueve números de este tipo por debajo de 10 9 . [17] [18] [19]
Si bien los pseudoprimos de Perrin son raros, se superponen con los pseudoprimos de Fermat . De los diecisiete números anteriores, cuatro también son fermatianos de base 2. Por el contrario, los pseudoprimos de Lucas están anticorrelacionados. [20] Presumiblemente, la combinación de las pruebas de Perrin y Lucas debería hacer una prueba de primalidad tan sólida como la confiable prueba BPSW que no tiene pseudoprimos conocidos, aunque a un costo computacional más alto.
Pseudocódigo
La prueba de primalidad de Adams y Shanks O (log n) Perrin de 1982. [21]
Dos matrices de enteros u(3) y v(3) se inicializan en los términos más bajos de la secuencia de Perrin, con índices positivos t = 0, 1, 2 en u( ) e índices negativos t = 0,−1,−2 en v( ).
El bucle principal de doble y suma , originalmente diseñado para ejecutarse en una calculadora de bolsillo HP-41C , calcula P(n) mod n y el inverso P(−n) mod n al costo de seis elevaciones al cuadrado modulares para cada bit de n. .
Los subíndices de los números de Perrin se duplican usando la identidad P(2t) = P 2 (t) − 2P(−t) . Las brechas resultantes entre P(±2t) y P(±2t ± 2) se cierran aplicando la relación definitoria P(t) = P(t − 2) + P(t − 3) .
Los valores iniciales son int u(0):= 3, u(1):= 0, u(2):= 2 son int v(0):= 3, v(1):=−1, v(2) := 1Pruebe el número positivo impar n entrada int n set int h:= bit más significativo de n para k:= h − 1 hasta 0 Duplica los índices de los seis números de Perrin. para yo = 0, 1, 2 temp:= u(i)^2 − 2v(i) (mod n) v(i):= v(i)^2 − 2u(i) (mod n) u(i):= temperatura fin paraCopie P(2t + 2) y P(−2t − 2) en los extremos de la matriz y utilícelos en la siguiente declaración if. tu(3):= tu(2) v(3):=v(2) Sobrescribir P(2t ± 2) con P(2t ± 1) temperatura:= u(2) − u(1) u(2):= u(0) + temperatura u(0):= temperatura Sobrescribir P(−2t ± 2) con P(−2t ± 1) temperatura:= v(0) − v(1) v(0):= v(2) + temperatura v(2):= temperatura si n tiene el bit k establecido , aumente los índices de ambos triples de Perrin en 1. para i = 0, 1, 2 tu(yo):= tu(yo + 1) v(yo):= v(yo + 1) fin para endiffin paraResultado imprimir v(2), v(1), v(0) imprimir u(0), u(1), u(2)
Si P(−n) = −1 y P(n) = 0 entonces n es un primo probable , es decir: realmente primo o un pseudoprimo de Perrin restringido.
Shanks y cols. observaron que para todos los pseudoprimos restringidos que encontraron, el estado final de los seis registros anteriores (la "firma" de n ) es igual al estado inicial 1, −1,3, 3,0,2. [22] Lo mismo ocurre con ≈ 1/6 de todos los primos, por lo que los dos conjuntos no se pueden distinguir basándose únicamente en esta prueba. [23] Para esos casos, recomiendan utilizar también la secuencia hermana de Narayana-Lucas con relación de recurrencia A(n) = A(n − 1) + A(n − 3) y valores iniciales
Se aplica la misma regla de duplicación y las fórmulas para llenar los espacios son
temperatura:= u(0) + u(1)u(0):= u(2) − temperaturau(2):= temperatura temperatura:= v(2) + v(1)v(2):= v(0) − temperaturav(0):= temperatura
Aquí, n es un primo probable si A(−n) = 0 y A(n) = 1 .
Kurtz et al. no encontró superposición entre los pseudoprimos impares para las dos secuencias por debajo de 50∙10 9 y supuso que 2,277,740,968,903 = 1067179 ∙ 2134357 es el número compuesto más pequeño que pasa ambas pruebas. [24]
Si también se utiliza la recurrencia de Pell-Lucas de tercer orden A(n) = 2A(n − 1) + A(n − 3) , este límite se elevará a 4.057.052.731.496.380.171 = 1424263447 ∙ 2848526893. [25]
Además, las raíces de la congruencia de la regla de duplicación distintas de −1 o 3 exponen números compuestos, como las raíces cuadradas no triviales de 1 en la prueba de Miller-Rabin . [26] Esto reduce el número de pseudoprimos restringidos para cada secuencia en aproximadamente un tercio y es especialmente eficiente en la detección de números de Carmichael . [27]
El pseudoprimo de Perrin restringido menos fuerte es 46672291 y los dos límites anteriores se expanden sucesivamente a 173,536,465,910,671 y 79,720,990,309,209,574,421. [28]
^ Ninguno de los 2402549 pseudoprimos de Lucas-Selfridge por debajo de 10 15 enumerados por Dana Jacobsen (2020) es también un pseudoprimo de Perrin.
^ Adams y Shanks (1982, pág.265, 269-270) harvtxt error: no target: CITEREFAdamsShanks1982 (help)
^ Adams y Shanks (1982, p. 275) , Kurtz, Shanks y Williams (1986, p. 694) . Esto fue confirmado posteriormente para n < 10· 14 por Steven Arno (1991). harvtxt error: no target: CITEREFAdamsShanks1982 (help) harvtxt error: no target: CITEREFKurtzShanksWilliams1986 (help)
^ La firma proporciona información discriminatoria sobre los dos tipos restantes de números primos. Por ejemplo, el pseudoprimo de tipo Q más pequeño 50,972,694,899,204,437,633 calculado por Holger Stephan (2019) está expuesto por las condiciones de firma 14a y 14c en Adams & Shanks (1982, p. 257) . harvtxt error: no target: CITEREFAdamsShanks1982 (help)
^ Kurtz, Shanks y Williams (1986, pág. 697) harvtxt error: no target: CITEREFKurtzShanksWilliams1986 (help)
^ Esteban (2019)
^ Adams y Shanks (1982, págs. 280-283) harvtxt error: no target: CITEREFAdamsShanks1982 (help)
^ La implementación AC/C++ de la prueba Perrin extendida se puede encontrar en la subsección final de una versión anterior de este artículo.
^ Esteban (2019)
Referencias
Lucas, E. (1878). "Teoría de las funciones numéricas simples y periódicas". Revista Estadounidense de Matemáticas (en francés). 1 (3). Prensa de la Universidad Johns Hopkins: 229–231. doi : 10.2307/2369311 . JSTOR 2369311.
Füredi, Zoltán (1987). "El número de conjuntos independientes máximos en gráficos conectados". Revista de teoría de grafos . 11 (4): 463–470. doi :10.1002/jgt.3190110403.
Grantham, Jon (2010). "Hay infinitos pseudoprimos de Perrin". Revista de teoría de números . 130 (5): 1117-1128. arXiv : 1903.06825 . doi :10.1016/j.jnt.2009.11.008.
Jacobsen, Dana (2020). "Tablas y estadísticas de pseudoprimes". ntheory.org . Consultado el 7 de marzo de 2024 . #LPSP Lucas-Selfridge
Stephan, Holger (2020). "Millones de pseudoprimos de Perrin, incluidos algunos gigantes". arXiv : 2002.03756v1 [matemáticas.NA].
Stephan, Holger (2019), Pseudoprimos de Perrin. Conjuntos de datos , Berlín: Instituto Weierstrass , doi : 10.20347/WIAS.DATA.4
enlaces externos
Jacobsen, Dana (2016). "Pruebas de primalidad de Perrin".
Wright, Colin (2015). "Encontrar pseudoprimos de Perrin".
"Secuencia de Perrin". MathPages.com .
"Lucas y Perrin Pseudoprimes". MathPages.com .
Holzbaur, Christian (1997). "Perrin Pseudoprimos".