Teoria Algebraica de Numeros - Carppintero, Tesis
Teoria Algebraica de Numeros - Carppintero, Tesis
Teoria Algebraica de Numeros - Carppintero, Tesis
Facultad de Ciencias
9 de julio de 2023
Índice general
Introducción 4
1
1.3. Alternativa al Teorema fundamental de la Aritmética . . . . . 24
2.1. Preliminares . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.1.3. Q-isomorfismos . . . . . . . . . . . . . . . . . . . . . . 31
2.1.4. Retículos en Rn . . . . . . . . . . . . . . . . . . . . . . 32
2.4. Unidades . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3. Aplicaciones a criptografía 50
2
3.1.1. RSA . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.2.1. Pseudoprimos . . . . . . . . . . . . . . . . . . . . . . . 59
Conclusiones 64
Bibliografía 65
3
Introducción
4
y, por supuesto, no puede esperarse que se resuelvan todos ellos sólo usando
métodos algebraicos.
5
Capítulo 1
6
N es un número primo distinto de los n dados, lo que lleva a contradicción.
CQD
hay por lo menos N primos, cada uno en uno de los intervalos (2n−1 , 2n ],
n ∈ {1, 2, . . . , N }.
7
Por hipótesis de inducción, como b < a, tenemos las descomposiciones
b = p1 · · · ps , a = p · p1 · · · ps
y a no tiene otra descomposición en la que aparezca p.
Sea otra descomposición prima de a en la que no aparece p:
a = q1 · · · qr
Tenemos que q1 ̸= p. Por definición de p, p < q1 , pues q1 ∈ D. Dado c =
q2 · · · qr , se define
a0 = a − pc = p(b − c) = (q1 − p)c
Tenemos que 1 ≤ a0 < a, y los divisores (b − c), (q1 − p) y c son menores
que a. Por hipótesis de inducción, ellos y a0 tienen una única factorización
por primos, y p aparece en la de a0 , por lo que aparece en la de c o en la de
(q1 − p) .
No puede aparecer en la de (q1 − p), pues eso requeriría que p|q1 , y es impo-
sible por ser primos distintos. Pero hemos supuesto que tampoco está en la
descomposición de c. Suponer que hay dos descomposiciones de a distintas
nos lleva a contradicción. CQD
8
Demostración: Sin perder generalidad, suponemos que a es positivo.
Entonces, es posible usar el método de inducción. Para a = 1, es evidente.
Asumimos que es cierto para a = n, entonces:
p
p
X p
(n + 1) = ni = np + pnp−1 + · · · + pn + 1
i=0
i
p p!
Para 0 < i < p, = tiene numerador divisible por p y deno-
i i!(p − i)!
minador
no divisible por p. Por el Teorema Fundamental de la Aritmética,
p
es divisible por p para dichos i. Por tanto
i
(n + 1)p ≡ np + 1 ≡ n + 1 (mód p)
9
Demostración: Los factores propios que dividen a 2p−1 q se pueden se-
parar según sean o no divisibles por q:
La suma de los factores no divisibles es
1 + 2 + 22 + · · · + 2p−1 = 2p − 1 = q
La de los divisibles es
Demostración:
Por ser x el único factor impar de la expresión, se tiene que 2n+1 −1 (también
x
impar) divide a x. Consideramos entonces y = n+1 y dividiento en la
2 −1
n+1
expresión anterior por 2 − 1, se obtiene la identidad
10
Teorema: Denotando Mn al n-ésimo número de Mersenne, 2n − 1, se
tiene que dado p primo y q divisor no trivial de Mp , q ≡ 1 (mód p).
n
Se denota al n-ésimo número de Fermat como Fn = 22 + 1. Esta
expresión da como resultado un número primo cuando n toma valores entre
0 y 4. Euler probó que F5 es divisible por 641 y se sospecha que solo hay un
número finito de primos de Fermat.
Por otro lado, la factorización de los números de Fermat es bastante sencilla,
como nos muestra el siguiente resultado:
11
n
Demostración: Sabemos que 22 ≡ −1 (mód p) y p es impar. Por tanto,
n+1
22 ≡1 (mód p)
(n − 1)! ≡ −1 (mód n)
Demostración
12
⇒ La demostración es trivial para n = 2. Sea n = p con p > 2 primo.
Si consideramos los enteros a tales que 1 < a < p − 1, todos tienen inver-
so multiplicativo módulo p. Como su unicidad es obvia, si cada inverso es
distinto del propio elemento, entonces se cancelarían todos los términos de
(p − 1)! (excepto el propio n − 1 ≡ −1) módulo p.
Supongamos que no se cumple, es decir, a2 ≡ 1 (mód p). Entonces
p|(a + 1)(a − 1), y por primalidad se sigue que a ≡ ±1 (mód p), esto es, a es
1 ó p − 1. CQD
para cualquier natural a. Parece funcionar entonces como test escoger al azar
a con 1 < a < p y comprobar si la congruencia se verifica. Sin embargo, exis-
ten números compuestos tales que para toda base coprima con ellos se verifica
la congruencia del Pequeño Teorema de Fermat. Estos son conocidos como
números de Carmichael.
El hecho de que haya infinitos de ellos vuelve al test descrito anteriormente
13
poco válido como test de primalidad.
14
Así, obtenemos las soluciones primitivas de la ecuación de Pitágoras,
x = 2mn, z = m2 + n2 , y = m2 − n2
n = x2 + y 2
T2 + 1 ≡ 0 (mód p)
15
Demostración: Para p=2 es evidente. Sea p = 4n + 1 para algún n > 0.
Por el teorema de Wilson
(p − 1)! = (p − 1)(p − 2) · · · 1 ≡ −1 (mód p)
Considerando
4n = p − 1 ≡ −1 (mód p)
4n − 1 = p − 2 ≡ −2 (mód p)
.
.
.
2n + 1 = p − 2n ≡ −2n (mód p)
En consecuencia,
(p − 1)! = (4n)! = 1 · 2 · · · (2n − 1)(2n)(2n + 1) · · · (4n) ≡
1 · 2 · · · (2n − 1)(2n)(−2n) · · · (−1) ≡ ((2n)!)2 (−1)2n ≡ −1 (mód p)
16
En cuanto a las soluciones para n cualquiera, se tiene que n es suma de
dos cuadrados si y solo si cada uno de sus factores primos congruentes con 3
módulo 4 aparece en la factorización de n con exponente par.
Demostración:
Definimos
los conjuntos
p − 1
A = a2 | 0 ≤ a ≤
2
2 p−1
B = −b − 1| 0 ≤ b ≤
2
17
Análogamente, los elementos de B módulo p también son distintos entre
sí. Como en A y en B hay p+1
2
elementos módulo p, por solo existir p enteros
módulo p entonces existen a2 ∈ A, −b2 − 1 ∈ B tales que
a2 + b2 + 1 ≡ 0 (mód p) CQD
x+y x−y
entonces x e y tienen la misma paridad y n = + es suma
2 2
de dos cuadrados.
Aplicando esta identidad a (a2 + b2 ) y a (c2 + d2 ) (ambas expresiones pares),
a2 + b 2 c2 + d 2
entonces tanto como son sumas de dos cuadrados. Como
2 2
m
conclusión inmediata, p es suma de cuatro cuadrados.
2
18
la parte de la izquierda es kpm2 . Como módulo m ax ≡ bw y cz ≡ dy,
ax − bw − cz + dy (y su cuadrado) es divisible por m2 . Se puede razonar de
forma análoga para ay + bz − cw − dz y az − by + cx − dw.
aw + bx + cy + dz ≡ w2 + x2 + y 2 + z 2 ≡ 0 (mód m), por lo que su cuadrado
es divisible por m2 .
Todo sumando de la expresión de la derecha, y en particular la expresión
2
entera,
2 m . Concluimos que kp es suma de cuatro cuadrados:
es divisible por
aw + bx + cy + dz
, etc. CQD
m
√
Cuando se elige un entero d > 1 libre de cuadrados, se tiene que Z[ d]
es un anillo con infinitas
√ unidades y en el que se puede definir una norma N
dada por N (x + y d) = x2 − dy 2 . Esto nos da una pista para la solución de
otra ecuación:
x2 − dy 2 = 1
∆ = b2 − 4ac
19
Tenemos el siguiente teorema:
Sean
b+ρ
r = aαδ + cβδ +
2
s = aδ + bδγ + cγ 2
2
20
El discriminante de (C) es (2r + ρ)2 − 4sn =∆ (D),por lo cual tenemos
∆−ρ
una solución de (B) dada por z = r, r2 + ρr − = sn
4
21
Demostración: Es evidente si p es factor de a, así que suponemos que a y
p son coprimos.
Por ser coprimos, ap−1 ≡ 1 (mód p) y como las únicas raíces de 1 módulo p
son ±1, a(p−1)/2 ≡ ±1 (mód p)
Deducimos
entonces que j es par y a es residuo cuadrático módulo p, es
a
decir, =1 CQD
p
Si al menos uno de los dos es congruente con 1 módulo 4, los símbolos son
iguales.
22
Este teorema junto con las propiedades vistas y el siguiente resultado
nos permitirá verificar de manera eficiente si un entero cualquiera es residuo
cuadrático de un primo cualquiera.
2
Concluimos que f (p) = por ser G distinto de 0, y f (p) es 1 si y solo
p
si p es congruente con ±1 módulo 8 CQD.
23
los de Legendre. Tomemos como ejemplo:
1003 201 3 67 401 401
= = = (−1)2 =
401 401 401 401 3 67
2 66 2 3 11 67 67 1 1
= = (−1) = = =1
3 67 67 67 67 3 11 3 11
Se sabe que en los anillos que poseen una norma se puede definir un
algoritmo de división, y como consecuencia se prueba la existencia de un
TFA. Tal norma es relativamente intuitiva cuando se consideran algunos
anillos de enteros algebraicos sobre Z (es decir, anillos formados por las raíces
de polinomios mónicos con coeficientes enteros), como ya se ha visto. Sin
embargo, no todos esos anillos son euclídeos, y tampoco es posible obtener
factorizaciones en primos de sus elementos.
Desde
√ aquí hasta el fin del capítulo se va a trabajar en un cuerpo
K = Q( d) y en su anillo de enteros sobre Z, para algún d libre de cuadrados.
Denotamos
√ este anillo OK Consideramos entonces la norma y traza de α =
x+y d
N (α) = x2 − dy 2
T (α) = 2x
24
√
Caractericemos
√ el anillo de enteros sobre Q( √d). Si d ≡ 2 o 3 (mód 4)
éste es Z[ d], y si d ≡ 1 (mód 4), este es Z[(1 + d)/2]
√
Demostración: Por ser Q( d) extensión de grado 2, tenemos que en
particular todos los elementos de su anillo de enteros son solución de una
ecuación mónica de coeficientes enteros y de √grado 2, x2 + bx + c = 0.
−b ± b2 − 4c
Consideramos sus posibles soluciones . Podemos considerar
2
que el término de la raíz es distinto de 0 para buscar elementos fuera de Z.
Consideramos 2 casos:
Los elementos generadores de los anillo de enteros en cada caso son los
buscados.
√ √
Dado α = x+y d, definimos el conjugado de α, α∗ = x−y d. Dado I un
ideal, se define el ideal I ∗ como el conjunto de conjugados de los elementos
de I. Se tiene que el conjugado de la suma y producto de ideales es la suma
y producto de los ideales conjugados, respectivamente.
25
Añadimos como notación, para α1 , . . . , αk ∈ OK ,
(α1 , . . . , αk ) = (α1 ) + · · · + (αk )
Demostración:
Como el producto de ideales está contenido en su intersección, I|J implica
J=IK y J está contenido en I.
Sean J, I ideales, J ⊆ I. Entonces JI ∗ ⊆ II ∗ = (N (I)). Para todo elemento
a ∈ JI ∗ , N(I) divide a a. Podemos considerar
1
K= JI ∗
N (I)
ideal contenido en OK . Se tiene que
1 N (I)
IK = I(JI ∗ ) = J =J
N (I) N (I)
Por definición, I|J.
26
propio no nulo de OK se puede escribir como producto de ideales primos,
siendo la factorización única salvo el orden.
27
Capítulo 2
Fundamentos algebraicos de la
teoría de números
2.1. Preliminares
28
2.1.1. Ideales fraccionarios
29
Clausura entera de A en R: Dicho subanillo de los enteros de R
sobre A.
30
Sean B un anillo y A un subanillo de B tal que B es A-módulo libre de
rango finito n. Para un sistema (x1 , . . . , xn ) ∈ B n , se define el discriminante
de (x1 , . . . , xn ) al elemento de A dado por D(x1 , . . . , xn ) = det(T r(xi xj ))
Proposición:
Pn Si (y1 , . . . , yn ) es otro sistema de elementos de B tal que
yi = i=1 aij xj , aij ∈ A, entonces D(y1 , . . . , yn ) = det(aij )2 D(x1 , . . . , xn )
2.1.3. Q-isomorfismos
31
ción sobre x, siendo esta σi (x) = xi , y siendo {xi }i=1,...,n el conjunto de raíces
del polinomio irreducible de x, que sabemos también que son distintas:
2.1.4. Retículos en Rn
32
Se tiene que µ(Pe ) es independiente Pn de la base elegida. En efecto, sea f =
(f1 , . . . , fn ) otra base de H. fi = j=1 aij ej , aij ∈ Z. Sabemos que el deter-
minante de la matriz (aij ) de cambio de base es una unidad en Z por lo que,
considerando que µ(Pf ) = |det(aij )|µ(Pe ), y la igualdad es evidente al ver
que las únicas unidades de Z son ±1.
33
c) Todo submódulo de M es de tipo finito
Demostración:
c ⇒ b)
S Sea (Mn )n≥0 una sucesión creciente de submódulos de M. Enton-
ces E = n≥0 Mn es un submódulo de M. Por c), E admite un sistema gene-
rador finito, {x1 , . . . , xq }. Para cada i ∈ {1, . . . , q}, ∃n(i) tal que xi ∈ Mn(i) .
Sea n0 = max{n(1), . . . , n(q)}. Claramente, E ⊆ Mn0 y E = Mn0 . Se con-
cluye que para n ≥ n0 , Mn = Mn0 y la sucesión es estacionaria a partir de
n0 .
Demostración:
34
de elección, existe f : S → S tal que f (x) > x. Por ser S no vacío, existe
t0 ∈ S. Definimos una sucesión (tn )n≥0 tomando tn+1 = f (tn ), y resulta ser
una sucesión creciente y no estacionaria, lo que contradice b). CQD
Demostración:
35
Proposición: Dados A un anillo noetheriano íntegramente cerrado, K su
cuerpo de fracciones de característica 0, L una extensión finita de K y A’ la
clausura entera de A en L, se tiene que A’ es un A-módulo de tipo finito y
un anillo noetheriano.
Este resultado se puede aplicar a los números enteros, siendo Z anillo íntegra-
mente cerrado, y obtenemos que todo anillo de enteros sobre una extensión
de Q es noetheriano.
xn + an−1 xn−1 + · · · + a1 x1 + a0 = 0 ai ∈ A
36
Con los siguientes teoremas se intenta obtener una descomposición en pro-
ducto de ideales primos de los ideales fraccionarios de los anillos de Dedekind.
37
b ⊈ Aa por ser n mínimo. Entonces ∃b ∈ b tal que b ∈
/ Aa.
Como mb ⊆ Aa, mb ⊆ Aa y mba ⊆ A. Por definición, ba−1 ∈ m′ . Sin
−1
38
Veamos ahora la unicidad. Supongamos p∈P pnp (b) = p∈P pmp (b) ,
Q Q
esto es, p∈P pnp (b)−mp (b) = A. Se separan los n(p) y los m(p) no nulos nega-
Q
tivos y positivos para obtener
Como última parte de la sección, veamos que dados dos ideales enteros
no nulos de A, N (ab) = N (a)N (b).
39
card(A/am) = card(A/a)card(a/am). Veamos que card(a/am) = card(A/m).
Como a/am es un A-módulo anulado por m, es un espacio vectorial sobre A/m
cuyos subespacios son sus A-submódulos de la forma q/am, donde q es un
ideal tal que am ⊆ q ⊆ a. Por la descomposición en ideales maximales de am
y a, no existen ideales comprendidos estrictamente entre ambos. Se deduce
que no existen subespacios vectoriales propios de a/am sobre A/m, por lo
que la dimensión de a/am como espacio vectorial sobre A/m es 1. Por tanto,
card(a/am) = card(A/m). CQD
40
Demostración: Como consecuencia directa del primer corolario, hay una
cantidad finita de enteros que se pueden corresponder con la norma de un
ideal b entero. Será suficiente ver que el conjunto de ideales enteros con norma
igual alguno de esos enteros es también finito. Sea b tal que card(A/b) = q
(con A el anillo de enteros de K), se sigue que q ∈ b; A/b es un grupo de orden
q, y q ∈ A. Se deduce fácilmente que en tal grupo cociente q se identifica con
0, por lo que q ∈ b.
Obtenemos entonces que b contiene a Aq, y por finitud de A/Aq (la norma de
todo ideal de A es finita), tenemos un número finito de ideales que contienen
a Aq. CQD
2.4. Unidades
41
de las unidades de K A* es isomorfo a Zr × G, siendo G un grupo cíclico
formado por las raíces de la unidad contenidas en K.
42
un único generador mayor que 1 denominado unidad fundamental de K
Demostración:
√
√ Consideramos K = Q[ d] con d ≥ 2 sin factores cuadrados y x = a +
b d una unidad √ de K. √ Sabemos que N(x) es una unidad en Z, entonces
N (x) = (a + b d)(a − b d)√= ±1. De esta igualdad se puede deducir que
x, −x, x−1 , −x−1 son ±a ± b d. Si x ̸= ±1, solo uno de ellos es mayor√que 1,
y las unidades mayores que 1 de K son las unidades de la forma a + b d con
a, b > 0. La existencia de un único generador es evidente por la isomorfía a
Z.
43
a
Representando la clase de equivalencia de (a, s) podemos identificar
s
S −1 A con na o
∈ K | a ∈ A, s ∈ S
s
y es un subanillo de K que contiene a A
a
Veamos la inclusión opuesta. Sea x ∈ b′ , x = . Entonces sx ∈ b′ , a ∈ b′
s
1
y a ∈ b ∩ A. Por tanto, x = · a ∈ A (b ∩ A). Deducimos b′ ⊆ A′ (b′ ∩ A).
′ ′ ′
s
La inyectividad de φ se deduce de la existencia de la aplicación θ(b) = A′ b,
que cumple θ ◦ φ = idA′ .
44
np o
Veamos que pA′ = | p ∈ p, s ∈ S .
′
P asi
Sea x ∈ pA . x = pi , pi ∈ p, ai ∈ A, s ∈ S, suma finita. Se reduce a
si
P bi
x= pi , bi ∈ p, s el denominador común de los si . Como p es ideal de
s
p P p
A, x = , p = bi pi ∈ p. Como S ∩ p = ∅, entonces 1 ∈
/ p y 1 ̸= .
s s
a b ab
Sigue ver que p′ es primo: sean , ∈ A′ tales que ∈ pA′ . Se tiene
s t st
ab p
que = , p ∈ p, u ∈ S y se deduce abu = pst ∈ p. Como p ∩ S = ∅, u ∈ /p
st u
a
y por ser p ideal primo, o bien a o bien b son elementos de p, de donde o
s
b
bien son elementos de pA′ = p′
t
45
la modificamos dividiendo por sn para obtener
bn + an−1 bn−1 + · · · + a0 = 0
n n−1
b an−1 b a0
+ + ··· + n = 0
s s s s
b
y es entero sobre S −1 A. De forma similar, si xs es un elemento de S −1 R
s
entero sobre S −1 A, su ecuación de dependencia entera es
x n c x n−1 c0
n−1
+ + ··· + =0
s tn−1 s t0
para alguna colección finita de ci ∈ A, ti ∈ S. Si se multiplica la expresión
anterior por (t0 · t1 . . . tn−1 )n , se obtiene una ecuación de dependencia entera
x · t0 . . . tn−1
de con coeficientes en A, por lo que este es entero sobre A
s
x 1 x · t0 · t1 . . . tn−1
y pertenece a B. Para finalizar, = es un
s t0 · t1 . . . tn−1 s
elemento de S −1 B y se concluye que S −1 B es la clausura entera de S −1 A en
S −1 R. CQD
46
ideales son de la forma buscada. Además p es primo por ser generador de un
ideal primo. CQD
⇒ Se tiene que p ⊆ D, pB ⊆ DB = D
47
B es un A-módulo de tipo finito y ambos A/p y B/Bi son cuerpos, se puede
considerar B/Bi como un espacio vectorial de dimensión finita sobre A/p,
denotando por fi a esta dimensión.
Cabe destacar que Bp ∩ A = p y B/Bp es también un espacio vectorial de
dimensión finita sobre A/p.
48
ducimos que pB ′ = qi=1 (B ′ Bi )ei .
Q
Los B ′ Bi son ideales primos no nulos de B’, pues Bi ∩ S = ∅ por ser
Bi ∩ A = p.
[B ′ /B ′ Bi : A′ /A′ p] = fi
Qq
Proposición: B/Bp es isomorfo a i=1 B/Bei i
49
Capítulo 3
Aplicaciones a criptografía
50
je. La seguridad se basa en parte de la dificultad de que el adversario pueda
llegar a obtener la clave privada.
51
Utilizaremos funciones de una vía que posean una trampa, es decir, cuya
inversa se pueda computar fácilmente al conocer alguna información previa,
sin la cual se comportaría como una función de una vía. Nos centraremos en
ejemplos de algunos esquemas cuya seguridad está basada en distintos pro-
blemas computacionalmente difíciles de resolver, en particular los problemas
de factorización, logaritmo discreto y el problema de la mochila. La trampa
en este caso sería, trivialmente, el haber construido el problema a partir de
una información conocida.
3.1.1. RSA
52
El resultado que nos garantiza que P ′ = P es una generalización del
Pequeño Teorema de Fermat
53
desincentiva recurrir a tipos de primos especiales, como por ejemplo los de
Mersenne. Sería válido como método para generar un primo aleatorio generar
un entero aleatorio m y luego tomar el menor primo p tal que m ≥ p. Un
procedimiento similar podría utilizarse para obtener el exponente de cifrado
e, comprobando la coprimalidad con φ(N ) en lugar de la primalidad.
54
no se ha descubierto un método sencillo para la operación inversa. Algunos
algoritmos para la resolución de logaritmos discretos son
55
Se eligen un primo p y g una raíz primitiva módulo p.
5. B devuelve a A M eA eB (mód p)
56
7. Por último, B descifra el mensaje, M (mód p) = M eB dB (mód p)
57
y en caso contrario ϵi = 0, y se considera que el volumen disponible se reduce
correspondientemente.
Esta diferencia de dificultad nos da una idea para un criptosistema de cla-
ve pública: el criptosistema knapsack (mochila en inglés) o sistema Merkle-
Hellman.
Pk−1
1. Se eligen una k-tupla supercreciente {vk−1 . . . v0 }, m ≥ i=0 vi y a un
entero coprimo con m y menor que m.
2. Se calculan b = a−1 (mód m) y la k-tupla {wi } dada por wi = avi
(mód m)
3. La clave pública es la k-tupla {wk−1 . . . w0 }
4. La clave privada es la tupla (b, m, {vi }k−1
i=0 )
58
3.2. Tests de primalidad en criptografía
3.2.1. Pseudoprimos
59
Si n es número de Carmichael entonces es producto de al menos 3
números primos.
60
2. Existe al menos un a para el que no se cumple (*).
El siguiente criterio está basado en que para un número primo p las únicas
raíces de la unidad módulo p son ±1.
Dada una base cualquiera a coprimo con n se tendría que an−1 ≡ 1 (mód n)
y se procedería a tomar las posibles raíces de esa potencia, las potencias
(n − 1)/2, (n − 1)/4, . . . , (n − 1)/2s , siendo n − 1 = 2s t con t impar. El
test en sí se realiza partiendo de at y elevando sucesivamente al cuadrado,
2 s
a2t , a2 t , . . . , a2 t . El test termina cuando se obtiene en uno de los cuadrados
s−1
−1, y el número n se declara compuesto si se llega a la potencia a2 t y esta
no es −1. Este es el llamado test de Miller-Rabin.
at ≡ 1 (mód n)
r
∃r tal que 0 ≤ r < s y a2 t ≡ −1 (mód n)
61
3.2.4. Algoritmo AKS
62
La comprobación de la primalidad de los r en el primer bucle está justi-
ficada por el siguiente lema que acota los números con las condiciones que
exigimos:
63
Conclusiones
64
Bibliografía
65
[10] E. Kranakis, Primality and Cryptography (John Wiley & Sons, 1986)
66