Method, software and apparatus for computing discrete logarithms modulo a prime

Information

  • Patent Grant
  • 12143490
  • Patent Number
    12,143,490
  • Date Filed
    Wednesday, February 26, 2020
    5 years ago
  • Date Issued
    Tuesday, November 12, 2024
    3 months ago
  • Inventors
    • Coraluppi; Giorgio (Pittsburgh, PA, US)
  • Original Assignees
  • Examiners
    • Abedin; Shanto
    Agents
    • Schwartz; Ansel M.
Abstract
A decoding apparatus having a non-transient memory in which is stored an electromagnetic signal representative of data which were encrypted relying on the difficulty of computing discrete logarithms. The decoding apparatus has a computer in communication with the memory that decodes the encrypted data in the memory by computing the data's discrete logarithm. The decoding apparatus has a display on which the decoded encrypted data are displayed by the computer. A method for decoding.
Description
I. FIELD OF THE INVENTION

The present invention considers the exponential congruence

a0x≡y0(mod p)  (1)

where p is prime and a0 is a primitive root modulo p. Since a0 is primitive, x and y0 are in a one-to-one correspondence for integer values in the range 1≤x, y0≤p−1 [3]. Let G denote the set of integers {1, 2, . . . , p−1} and let |G| denote their number. Given p and a0 and given y0 in G, it is desired to find x modulo p−1. The integer x is usually referred to as the discrete logarithm of y0 in base a0 modulo p. (As used herein, references to the “present invention” or “invention” relate to exemplary embodiments and not necessarily to every embodiment encompassed by the appended claims.)


BACKGROUND OF THE INVENTION

Pohlig and Hellman discussed the significance of this problem for cryptographic systems [3]. It was concluded by Pohlig and Hellman that, if p−1 has only small prime factors, x can be computed in a time of the order of log2 p. However, if p−1 has a large prime factor p′, the search for x requires a time of the order p′·log p and may be untractable. As an illustration, Pohlig and Hellman presented two large primes of the form p=2·p′+1, where p′ is also prime and where

p′=213·5·7·11·13·17·19·23·29·31·37·41·43·47·53·59+1  (2)

or

p′=2121·52·72·112·13·17·19·23·29·31·37·41·43·47·53·59+1  (3)


In general, let p=2·p′+1, where p′ is prime and

p′−1=2ε0·q1ε1·q2ε2· . . . ·qiεi· . . . ·qhεh,  (4)

where ε0≥1 and, for 1≤i≤h, qi denotes an odd prime and εi>0. Also, for 1≤i<h, 2<qi<qi+1.


NOTE 1: Pohlig and Hellman observed that q1≠3. In fact p=2·p′+1=2·(p′−1)+3. Since p is prime, it must be gcd (3, p′−1)=1.


NOTE 2: Let X denote the set of elements of G which are relatively prime to p−1 and let A denote the set of primitive roots modulo p. Then |X|=|A|=φ(p−1) , where φ(n) denotes the Euler totient function.


NOTE 3: The elements of X form a commutative (abelian) group under the operation of multiplication modulo p−1. An integer m≥1 has a primitive root if and only if m=1, 2, 4, pd or 2·pd, where p is prime number and a is a positive integer [1, p. 211]. When X is cyclic, there exist integers ρ which are primitive roots of X modulo p−1. When primitive roots of X exist, let Y denote the set of elements of X which are primitive roots of X modulo p−1.


NOTE 4: Section VIII below shows that, when p′−1 can be described as in (4), X is cyclic only if ε0<3


BRIEF SUMMARY OF THE INVENTION

The present invention introduces an algorithm which, when p=2·p′+1, p′ is prime and p′−1 contains only small prime factors, produces the solution of (1) in a time of the order of loglog p·log2p.


The present invention pertains to a decoding apparatus. The decoding apparatus comprises a non-transient memory in which is stored an electromagnetic signal representative of data which were encrypted relying on the difficulty of computing discrete logarithms. The decoding apparatus comprises a computer in communication with the memory that decodes the encrypted data in the memory by computing the data's discrete logarithm. The decoding apparatus comprises a display on which the decoded encrypted data are displayed by the computer.


The present invention pertains to a method for processing an electromagnetic signal representative of encrypted data which were produced relying on the difficulty of computer discrete logarithms, comprising a first computer. The method comprises the steps of storing the encrypted data in a non-transient memory of a second computer. There is the step of performing with the second computer in communication with the memory the computer-generated steps of decoding the encrypted data in the memory by computing the data's discrete logarithms, and displaying on a display the decoded data.


The present invention pertains to a computer program stored in a non-transient memory for decoding an electromagnetic signal which is encrypted relying on the difficulty of computing discrete logarithms. The program has the computer-generated steps of storing the encrypted data in a non-transient memory. There is the step of decoding the encrypted data in the memory by computing the data's discrete logarithms. There is the step of displaying on a display the decoded data.


The present invention pertains to a method for reducing the complexity of an exponential congruence, preferably for decoding, which is defined modulo p, where p=2·p′+1, p′ is also a prime and p′−1 contains only factors which are smaller than 100,000. The method comprises the steps of executing with a computer a sequence of reversible transformations supported by a non-transient memory in such a way that the exponential congruence modulo p is restated as a problem involving new relationships modulo p and a concurrent independent congruence modulo p−1. There is the step of reporting the restated problem on a display.


The present invention pertains to a method for decoding. The method comprises the steps of selecting with a computer primitives of sub-groups of a group stored in a non-transient memory, where the group is defined modulo φ(p−1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product. There is the step of reporting the exponents on a display.





BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING

In the accompanying drawings, the preferred embodiment of the invention and preferred methods of practicing the invention are illustrated in which:



FIG. 1 is a block diagram of the apparatus of the claimed invention.



FIG. 2 is a representation of ρ1x1·ρ2x2 (mod 70) using orthogonal primitives.





DETAILED DESCRIPTION OF THE INVENTION

Referring now to the drawings wherein like reference numerals refer to similar or identical parts throughout the several views, and more specifically to FIG. 1 thereof, there is shown a decoding apparatus 10. The decoding apparatus 10 comprises a non-transient memory 14 in which is stored an electromagnetic signal representative of data which were encrypted relying on the difficulty of computing discrete logarithms. The decoding apparatus 10 comprises a computer 12 in communication with the memory 14 that decodes the encrypted data in the memory 14 by computing the data's discrete logarithm. The decoding apparatus 10 comprises a display 18 on which the decoded encrypted data are displayed by the computer 12.


The computer 12 may reduce the complexity of an exponential congruence which is defined modulo p, where p=2·p′+1, p′ is also a prime and p′−1 contains only factors which are smaller than 100,000, and executes a sequence of reversible transformations supported by the non-transient memory 14 in such a way that the exponential congruence modulo p is restated as a problem involving new relationships modulo p and a concurrent independent congruence modulo p−1. The computer 12 may select primitives of sub-groups of a group stored in the non-transient memory 14, where the group is defined modulo φ(p−1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product.


The present invention pertains to a method for processing an electromagnetic signal representative of encrypted data which were produced relying on the difficulty of computing discrete logarithms. The method comprises the steps of producing the electromagnetic signal by a first computer 12. There is the step of providing the signal to a second computer 22 through an input 20 of the second computer 22. The input 20 can be a keyboard in communication with the second computer 22 or a memory port, such as a USB port that receives a flash drive or a CD reader that receives a CD with the signal; or the input 20 can be a network interface card in communication with the second computer 22 having a network port which is in communication with a network 24 over which the signal is transmitted from the first computer 12. The second computer 22 obtains the signal from the network 24 through the input 20 of the second computer 22. There is the step of storing the encrypted data in a non-transient memory 14 of a second computer 22. There is the step of performing with the second computer 22 in communication with the memory 14 the computer-generated steps of decoding the encrypted data in the memory 14 by computing the data's discrete logarithms, and displaying on a display 18 the decoded data.


The performing step may include the steps of reducing the complexity of an exponential congruence which is defined modulo p, where p=2·p′+1, p′ is also a prime and p′−1 contains only factors which are smaller than 100,000. There may be the step of executing with the computer a sequence of reversible transformations supported by a non-transient memory 14 in such a way that the exponential congruence modulo p is restated as a problem involving new relationships modulo p and a concurrent independent congruence modulo p−1. There may be the step of reporting the restated problem on a display 18. The performing step may include the step of selecting with the computer primitives of sub-groups of a group stored in the non-transient memory 14, where the group is defined modulo φ(p−1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product.


The present invention pertains to a computer program 16 stored in a non-transient memory 14 for decoding an electromagnetic signal which is encrypted relying on the difficulty of computing discrete logarithms. The program has the computer-generated steps of storing the encrypted data in a non-transient memory 14. There is the step of decoding the encrypted data in the memory 14 by computing the data's discrete logarithms. There is the step of displaying on a display 18 the decoded data.


The decoding step may include the steps of reducing the complexity of an exponential congruence which is defined modulo p, where p=2·p′+1, p′ is also a prime and p′−1 contains only factors which are smaller than 100,000. There may be the step of executing with the computer a sequence of reversible transformations supported by a non-transient memory 14 in such a way that the exponential congruence modulo p is restated as a problem involving new relationships modulo p and a concurrent independent congruence modulo p−1.


The decoding step may include the steps of selecting with the computer primitives of sub-groups of a group stored in the non-transient memory 14, where the group is defined modulo φ(p −1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product.


The present invention pertains to a method for reducing the complexity of an exponential congruence, preferably for decoding, which is defined modulo p, where p=2·p′+1, p′ is also a prime and p′−1 contains only factors which are smaller than 100,000. The method comprises the steps of executing with a computer a sequence of reversible transformations supported by a non-transient memory 14 in such a way that the exponential congruence modulo p is restated as a problem involving new relationships modulo p and a concurrent independent congruence modulo p−1. There is the step of reporting the restated problem on a display 18.


The present invention pertains to an apparatus 10 for reducing the complexity of an exponential congruence, preferably for decoding, which is defined modulo p, where p=2·p′+1, p′ is also a prime and p′−1 contains only factors which are smaller than 100,000. The apparatus 10 comprises a non-transient memory 14. The apparatus 10 comprises a computer in communication with the non-transient memory 14 which executes a sequence of reversible transformations supported by the non-transient memory 14 in such a way that the exponential congruence modulo p is restated as a problem involving new relationships modulo p and a concurrent independent congruence modulo p−1. The apparatus 10 comprises a display 18 on which the restated problem is reported.


The present invention pertains to a computer program 16 stored in a non-transient memory 14 for reducing the complexity of an exponential congruence, preferably for decoding, which is defined modulo p, where p=2·p′+1, p′ is also a prime and p′−1 contains only factors which are smaller than 100,000. The program comprises the computer generated steps of executing a sequence of reversible transformations supported by a non-transient memory 14 in such a way that the exponential congruence modulo p is restated as a problem involving new relationships modulo p and a concurrent independent congruence modulo p−1. There is the step of reporting the restated problem on a display 18.


The present invention pertains to a method for decoding. The method comprises the steps of selecting with a computer primitives of sub-groups of a group stored in a non-transient memory 14, where the group is defined modulo φ(p−1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product. There is the step of reporting the exponents on a display 18.


The present invention pertains to a computer program 16 stored in a non-transient memory 14 for decoding. The program comprises the computer generated steps of selecting primitives of sub-groups of a group stored in a non-transient memory 14, where the group is defined modulo φ(p−1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product. There is the step of reporting the exponents on a display 18.


The present invention pertains to an apparatus 10 for decoding. The apparatus 10 comprises a non-transient memory 14. The apparatus 10 comprises a computer in communication with the memory 14 which selects primitives of sub-groups of a group stored in the non-transient memory 14, where the group is defined modulo φ(p−1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product. The apparatus 10 comprises a display 18 in communication with the computer on which the exponents are reported.


In the operation of the invention, the following is a description of the solution of (1).


II. THE CASE WHEN ε0=1. A RESTATEMENT

1) Step One. Definition of “Superprimitives”.


In general in (1) a0 is not a primitive root of X modulo p−1. It is convenient to restate (1) in such a way that on the LHS of (1) a0 be replaced by a primitive of X modulo p−1.


If ρ denotes a primitive of X modulo p−1, consider the process of raising both sides of (1) by ρl. As l increases, a0ρl modulo p traces an orbit of primitives modulo p.


If p is large, and if p=2·p′+1, where p′ is also prime, approximately half of the elements of G are elements of this orbit [2, p. 269].


For some integer {tilde over (l)}, a0ρl is also a primitive of X modulo p−1. In this case, define









{






a
0

ρ

i
_





a

(

mod

p

)








y
0

ρ

i
¯





y

(

mod

p

)






.





(
5
)








Then

ax≡y(mod p).  (6)


An integer which is a primitive root of p and also a primitive of X modulo p−1 will be referred to as a superprimitive of p.


Table I shows the superprimitives of a set of small primes (ε0≤2).


Table II shows some relevant variables for such primes.










TABLE I





p
Superprimitives of p
















11
7


23
7, 17, 19


47
5, 11, 15, 19, 33, 43


59
11, 31, 37, 39, 43, 47, 55


107
5, 21, 31, 45, 51, 55, 65, 67, 71, 73, 103


167
5, 13, 15, 35, 39, 43, 45, 53, 55, 67, 71, 73, 79, 91, 101, 103, 105,



117, 125, 129, 135, 139, 143, 145, 149, 155, 159, 163


263
29, 57, 67, 85, 87, 97, 115, 119, 127, 130, 139, 141, 161, 171, 185,



197, 213, 219, 227, 229, 237, 241, 247, 251, 255, 257, 259


347
5, 7, 17, 19, 45, 63, 65, 69, 79, 91, 97, 101, 103, 111, 123, 125, 141,



145, 153, 155, 165, 171, 175, 191, 193, 203, 215, 217, 221, 223, 231,



239, 245, 247, 283, 301, 307, 335, 343,


359
7, 21, 35, 53, 63, 71, 97, 103, 105, 109, 113, 119, 137, 143, 157, 159,



163, 167, 175, 189, 197, 209, 211, 213, 223, 251, 257, 263, 265, 269,



271, 277, 291, 293, 299, 309, 311, 313, 315, 319, 327, 329, 339,



341, 343, 349, 353, 355


383
33, 35, 47, 53, 61, 83, 89, 91, 95, 99, 105, 111, 123, 127, 131, 141,



145, 151, 157, 167, 179, 181, 183, 187, 233, 247, 249, 253, 267, 285,



297, 307, 315, 337, 355, 359, 365, 367, 369, 379


479
13, 19, 37, 39, 41, 43, 47, 53, 57, 59, 65, 79, 95, 117, 119, 123, 129,



143, 149, 159, 171, 177, 179, 185, 191, 205, 209, 223, 227, 235, 237,



265, 281, 285, 295, 317, 325, 333, 351, 353, 357, 369, 379, 387, 391,



395, 423, 429, 433, 447, 449, 451, 461, 463, 467, 469, 473, 475


503
19, 29, 37, 53, 55, 57, 71, 87, 107, 109, 111, 127, 133, 137, 139, 159,



163, 165, 167, 191, 193, 203, 213, 215, 239, 269, 277, 295, 305, 307,



313, 321, 327, 333, 341, 347, 349, 371, 381, 385, 399, 409, 417, 419,



437, 453, 457, 461, 467, 471, 475, 479, 481, 485, 487, 489, 495, 499


587
5, 11, 13, 19, 23, 41, 45, 85, 99, 103, 105, 111, 117, 125, 127, 131,



139, 157, 171, 173, 183, 187, 207, 213, 215, 221, 227, 231, 241, 245,



251, 259, 265, 263, 273, 275, 291, 295, 321, 323, 325, 327, 335, 337,



341, 365, 367, 369, 373, 391, 399, 403, 405, 415, 427, 435, 467,



473, 475, 483, 487, 497, 523, 539, 541, 557, 559, 583




















TABLE II





p
|A| = |X|
|Y|
|A ∩ Y|
|A|/(A ∩ Y)



















11
4
2
1
4.0000


23
10
4
3
3.3333


47
22
10
6
3.6666


59
28
12
7
4.0000


107
52
24
11
4.7273


167
82
40
28
2.9286


263
130
48
26
5.0000


347
172
84
39
4.4103


359
178
88
48
3.7083


383
190
72
40
4.7500


479
238
96
58
4.1034


503
250
100
58
4.3103


587
292
144
68
4.2941










NOTE 1: If ε0≤2, X is cyclic and there exist an integer ρ which is a primitive root of X modulo p−1. If p is large, to determine ρ it is sufficient to select any random integer and to verify that a) ρ is an element of X, which means that it is relatively prime to p−1, and b) ρ is an element of Y, which means that it is relatively prime to φ(p−1)=p′−1. The process of producing p should not be long, because, if p is large, the probability that two integers be prime to one another is 6/π2 [2, p. 269]. Thus, the probability that an integer chosen at random be prime to p−1 and p′−1 is approximately (6/π2)2 or 1/2.7055.


NOTE 2: The ratio |A|/|(A∩Y)| is relevant because it is related to the number of trials which should be expected when is employed in the search for a.


NOTE 3: The ratio |A|/|(A∩Y)| may grow when p increases. As an example, when p=6466463=2·p′+1 and p′=2·5·7·11·13·17·19+1, |A|/|A ∩Y|=7.7931.


NOTE 4: Comparing the data for p2=6466463 and p1=587, observe that, when p1 is replaced by p2, the ratio |A|/|(A∩Y)| is multiplied by a factor of 7.7931/4.2941, which is less than 2, while 6466463 is greater than 5872.


2) Step Two


In general in (6) y is not a primitive root modulo p. It is convenient to restate (6) in such a way that the RHS of (6) be a primitive modulo p. This can be accomplished by multiplying both sides of (6) by a, a sufficient number of times until the desired condition is satisfied. If after {tilde over (r)} iterations this condition is satisfied, let









{





b



a

r
~


·

y

(

mod

p

)








s


x
+


r
~

(

mod


(

p
-
1

)


)







.





(
7
)








Then

as≡b(mod p).  (8)


After this restatement the search for x is conducted in a smaller, more structured environment. Since b is a primitive modulo p and a is a primitive of X modulo p−1, s is relative prime to p−1 and can be represented as follows









{





s



a
t

(

mod



(

p
-
1

)


)








a

a





b

(

mod


p

)





,





(
9
)








where t denotes an integer and 0≤t<φ(p−1).


3) Step Three


Consider the process of raising the second of (9) to au modulo p. Let d denote the least positive residue modulo p of the corresponding RHS of (9). As u increases, the integer d describes an orbit of primitives modulo p. It is desired that d be also a primitive of X modulo p−1. If, after ũ operations this condition is satisfied, define

aaν≡d(mod p)  (10)

where









{





v


t
+


u
~

(

mod



φ

(

p
-
1

)


)








d



b

a

u
~



(

mod


p

)





.





(
11
)







Consider the integer ddw, where w denotes an integer. Since d is a primitive modulo p and a primitive of X modulo p−1, when w varies ddw traces an orbit which contains all the primitives modulo p, including a.


Therefore, there does exist an integer w such that

a≡ddw(mod p).  (12)


4) Conclusion


The exponential congruence (1) is referred to as a “one-way” transaction, meaning that, when x is known, it is easy to compute a0x modulo p, while, when y0 is known, the computation of x may be untractable. The restatement introduced by this section produces the congruences (10) and (12), which have similar structure and comparable complexity.


In order to determine the relationship between ν and w, raise (12) to aν modulo p. It will be

aaν≡daν·dw(mod p)  (13)

whence, by (10),

aν·dw≡1(mod p−1)  (14)

or

aν≡d−w(mod p−1).  (15)

As a conclusion: ν and w are exponents of known superprimitives of p, a and d, respectively. The integers aν and dw are related in a congruence which is defined modulo p−1.


NOTE 1: In general, given (1), the integers a and d which result from the proposed restatements are not unique.


NOTE 2: In principle, it would be possible to explore the case when a is a superprimitive of p and p−1. As an example, 19 is a superprimitive of 47 and 23. However, not all primes have superprimitives modulo p and modulo p−1.


III. ORTHOGONAL PRIMITIVES when ε0=1

Refer to (15). Let V denote the set of integers ν and w (1≤ν, w≤p′−1) which are candidate solutions of (10), (12), and (15) and let |V| denote their number.


A) It is desired to represent V as the direct product of distinct subsets of V, each one associated with one of the factors (qiε1 or 2ε0) of p′−1.


B) It is desired to partition and process independently the corresponding sets of candidate solutions.


To reach these aims:


A) The number of significant candidate elements associated with each of such sets is φ(qiε1) . Then the total number of candidate elements, say |V|, would be










|
V
|

=


φ

(


p


-
1

)

=


2


ε
0

-
1


·




i
=
1

h




φ

(

q
i

ε
i


)

.








(
16
)







The candidate elements associated to φ(qiε1), say ρlν1, are relatively prime to qi and can be represented as the elements of a cyclic group having ρi as its generator.


Notice that, thus far, nothing has been stated concerning the divisibility of ρi by qj when i≠j.


B) Consider the case when ν and w are represented as the direct product of their component subgroups. In the case when ε0=1, ε0−1=0. In order to process independently the cyclic subsets of V, consider the case when the primitive of the cycle i is defined as follows:










ρ
i

=


1
+


λ
i

·


φ

(

p
-
1

)

/

q
i

ε
i





=


σ
i

+


μ
i

·

q
i

ε
i









(
17
)








where σi denotes any primitive modulo qi and (λi, μi) denotes a pair of integers. Given σi, the pair (λi, μi) can be any one of the solution pairs of the following:

σi−1+μi·qiε1i·φ(p−1)/qiε1.  (18)

Given any solution pair ({tilde over (λ)}i, {tilde over (μ)}i), its substitution into (17) produces ρi modulo φ(p−1). After this restatement, ρi is relatively prime to φ(p−1).


Consider the case when p′−1 has a structure of the form (2) or (3), that is

    • a) 5 is the smallest odd prime divisor of p′−1, and
    • b) each divisor qi is the smallest odd prime greater than qi−1.


Under these conditions, all the odd prime divisors of φ(p′−1), with the exception of 3, are also divisors of φ(p−1). It is possible to select σi in such a way that ρi is not a multiple of 3. In this case, ρiν1 is relatively prime with φ(p−1) and φ(p′−1).


Thus, when p′−1 has the structure of (8) and (10) and ν is relatively prime with φ(p−1) and 3, it is possible to represent ν and w as follows









{




υ





i
=
1

h



ρ
j

υ
i


(

mod



φ

(

p
-
1

)


)









w





l
=
1

h



ρ
i

w
i


(

mod



φ

(

p
-
1

)


)




,








(
19
)








where νi and wi denote integers defined modulo φ(qiε1).


It will be









{







φ

(

p
-
1

)


q
i

ε
i



·

ρ
j






φ

(

p
-
1

)


q
i

ε
i





(

mod



φ

(

p
-
1

)


)



for


i


j










φ

(

p
-
1

)


q
i

ε
i



·

ρ
j







φ

(

p
-
1

)


q
i

ε
i



·


σ
j

(

mod



φ

(

p
-
1

)


)




for


i


=

j
.









(
20
)







The congruences (20) define the orthogonality between ρi and ρj, for i≠j, and validate the definition of ρi offered by (17).


Notice that the definitions (17) imply that










ρ
i

φ

(

q
l

ε
l


)




1



(

mod



φ

(

p
-
1

)


)

.






(
21
)







In fact,










ρ
i

φ

(

q
i

ε
i


)


=



(


σ
i

+


μ
i

·

q
i

ε
i




)


φ

(

q
i

ε
i


)




1
+


χ
i

·


q
i

ε
i


(

mod



q
i

ε
i



)








(
22
)








and also, for all positive integers n,










ρ
i
n

=



(

1
+


λ
i

·


φ

(

p
-
1

)

/

q
i

ε
i





)

n

=

1
+



ψ
i

·

φ

(

p
-
1

)




q
i

ε
i









(
23
)








for some χi and Ψi integers. Combining (22) and (23), (21) follows.


Refer to Section I of the Appendix.


IV. THE RELATIONSHIP BETWEEN νi AND wi MODULO φ(qiεi) WHEN ε0=1

Using orthogonal primitives (17), consider raising (15) to







φ

(

p
-
1

)


q
i

ε
i







modulo p−1.


It will be











(

a


φ

(

p
-
1

)


q
i

ε
i




)


σ
i

v
i







(

d


φ

(

p
-
1

)


q
i

ε
i




)


-

σ
i

w
i







(

mod



(

p
-
1

)


)

.






(
24
)







Let









{






α
i




a


φ

(

p
-
1

)


q
i

ε
i





(

mod


(

p
-
1

)


)








δ
i




d


φ

(

p
-
1

)


q
i

ε
i




(

mod

(

p
-
1

)

)






.





(
25
)







Then










α
i

σ
i

v
i







δ
i

-

σ
i

w
i




(

mod

(

p
-
1

)

)

.





(
26
)







This congruence establishes a relationship between νi and wi which does not depend on any of the values of νj and wj, for i≠j . In fact, given {tilde over (v)}i, the condition (26) defines {tilde over (w)}i modulo φ(qiε1).


NOTE 1: In (26), let wi, m denote the value of wi when νi≡0(mod φ(qiε1)). Then










α
i

-
1






δ
i

σ
i

w

i
,
m




(


mod


p

-
1

)

.





(
27
)








Likewise, let νi, m denote the value of νi when wi≡0(mod φ(qiε1)). Then










α
i

σ
i

v

i
,
m








δ
i

-
1


(


mod


p

-
1

)

.





(
28
)








NOTE 2: Consider the case when all the νj's are congruent to zero modulo φ(qjεj). In this case, from (15),










a

-
1






d







j
=
1

h



ρ
j

w

j
,
m





(


mod


p

-
1

)

.





(
29
)








Let









W
=




j
=
1

h



ρ
j

w

j
,
m



.






(
30
)








Then

a−1≡dW(mod p−1).  (31)


Likewise, consider the case when all the wj's are congruent to zero modulo φ(qj249j). In this case










a





j
=
1


h


ρ
j

υ

j
,
m









d

-
1


(


mod


p

-
1

)

.





(
32
)







Let









U
=




j
=
1

h



ρ
j

v

j
,
m



.






(
33
)







Then

aU≡d−1(mod p−1).  (34)

NOTE 3: Consider the case when all the νj's other than νi are congruent to zero modulo φ(qjεj) .


Assume {tilde over (ν)}i≢0 (mod φ(qiεi)) and νj≡0 (mod φ(qjεj)) for i≠j . In this case, using (31), it will be










a

-

ρ
ι


υ
~

i








d


ρ
ι


w
~

i


·

(

W
/

ρ
ι

w

ι
,
m




)



(


mod


p

-
1

)

.





(
35
)








By (31)










d

W
·

ρ
ι

v
i







d

W
/

ρ
ι

w

ι
,
m





·



d

ρ
ι

w
i



(


mod


p

-
1

)

.






(
36
)








Then

ρiwl, m·ρivi≡ρiwi(mod φ(p−1)),  (37)

or

wi, m+vi−wi≡0(mod φ(qiεi)).  (38)

NOTE 4: In a similar fashion, consider the case when all the wj other then wi are congruent to zero modulo φ(qjεj).


Assume wj≢0 (mod qiεi) for j≢i. In this case, using (32), it will be












a

ρ
ι

υ
i



·
U

/

ρ
i

V

ι
,
m






d

-

ρ
ι

w
i








(
39
)








By (32)

ρivl≡ρivl, m·ρiwl  (40)

or

vi≡wi+vi, m(mod φ(qjεj)).  (41)


Compare (38) and (41):

wi, m=vi, m(mod φ(qjεj)).  (42)

As a conclusion:


There exists a linear relationship between vi and wi modulo φ(qjεj).


V. THE TRIAD

Given a pair (a, d), after the determination of the corresponding relationships (vi, wi), it is interesting to explore comparable relationships for a pair (d, g), where g is a superprimitive modulo p.


Let (vi(A, D), wi(A, D)) denote the pair (vi, wi) associated with the pair (a, d).


Likewise, let (vi(D, G), wi(D, G)) denote the pair (vi, wi) associated with the pair (d, g).


Then, by (26),









{







α
i

σ
i


v
i

(

A
,
D

)



·

δ
i

σ
i


w
i

(

A
,
D

)






1


(


mod


p

-
1

)










δ
i

σ
i


v
i

(

D
,
G

)



·

γ
i

σ
i


w
i

(

D
,
G

)






1



(


mod


p

-
1

)






,





(
43
)








where









γ
i





g


φ

(

p
-
1

)


q
i

ε
i




(

mod

(

p
-
1

)

)

.





(
44
)







Raising the first of (43) to σivi(D, G) and the second to σiwi(A, D) modulo p−1 yields









{






α
i

σ
i



v
i

(

A
,
D

)

+


v
i

(

D
,
G

)




·

δ
i

σ
i



w
i

(

A
,
D

)

+


v
i

(

D
,
G

)







1


(

mod

(

p
-
1

)

)










δ
i

σ
i



w
i

(

A
,
D

)

+


v
i

(

D
,
G

)




·

γ
i

σ
i



w
i

(

A
,
D

)

+


w
i

(

D
,
G

)







1


(

mod

(

p
-
1

)

)










(
45
)








whence









α
i

σ
i



v
i

(

A
,
D

)

+


v
i

(

D
,
G

)








γ
i

σ
i



w
i

(

A
,
D

)

+


w
i

(

D
,
G

)




(

mod

(

p
-
1

)

)

.





(
46
)








This congruence establishes a relationship between νi and wi which does not depend on any of the values of νj and wj, for i≠j. However, this relationship does not identify the value of νi which is consistent with (6).


NOTE 1: In general, in (27) wi, m depends on the pair (a, d). Therefore, it will be necessary to distinguish wi, m (A, D) from wi, m (D, G). The same observation applies to vi, m (A, D) and vi, m (D, G).


NOTE 2: a and d are primitives modulo p−1. Therefore, they are relatively prime with φ(p−1). When a or d are raised to a divisor of φ(p−1), such as φ(p−1)/qiε1, they produce primitives modulo φ(p−1) for the sets {σiv1} and {σiw1}, respectively.


VI. INVERTIBLE SUPERPRIMITIVES

1) Introduction


Thus far, given a pair (a, d), the relationships between the sets {νi} and {wi} have been based on the condition (14).


It should be emphasized the fact that (14) is defined modulo (p−1) and is only a necessary condition for the solution of (10) and (12). It does not produce the solution of (10) and (12) because no adequate link has been produced between the domain modulo p−1 and the domain modulo p.


To this end, this Section introduces the notion of “Invertible Superprimitive”.


2) Definition


A superprimitive of p is defined as invertible if its inverse modulo p is also a superprimitive of p. In general, only some of the superprimitives are invertible. Table III shows the invertible superprimitives of the set of primes which are included in Tables I and II. This congruence establishes a relationship between νi and wi which does not depend on any of the values of νj and wj, for i≠j.














TABLE III









Number of





Invertible
Number of
Invertible



p
Superprimitives of p
Superprimitives
Superprimitives
Superprimitives
Ratio




















11
7

1
0
0.


23
7, 17, 19
17, 19
3
2
0.666667


47
5, 11, 15, 19, 33, 43
5, 19
6
2
0.333333


59
11, 31, 37, 39, 43, 47, 55
11, 43
7
2
0.285714


107
5, 21, 31, 45, 51, 55, 65, 67, 71, 73, 103
21, 51
11
2
0.181818


167
5, 13, 15, 35, 39, 43, 45, 53, 55, 67, 71, 73, 79,
5, 35, 43, 67, 101, 105,
28
10
0.357143



91, 101, 103, 105, 117, 125, 129, 135, 139, 143,
125, 129, 145, 163






145, 149, 155, 159, 163






263
29, 57, 67, 85, 87, 97, 115, 119, 127, 130, 139,
29, 97, 115, 127, 141,
26
12
0.461538



141, 161, 171, 185, 197, 213, 219, 227, 229,
197, 219, 241, 247,






237, 241, 247, 251, 255, 257, 259
251, 257, 259





347
5, 7, 17, 19, 45, 63, 65, 69, 79, 91, 97, 101, 103,
17, 69, 79, 103,
39
8
0.205128



111, 123, 125, 141, 145, 153, 155, 165, 171,
123, 171, 245, 283






175, 191, 193, 203, 215, 217, 221, 223, 231,







239, 245, 247, 283, 301, 307, 335, 343,






359
7, 21, 35, 53, 63, 71, 97, 103, 105, 109, 113,
157, 197, 209, 213,
48
16
0.333333



119, 137, 143, 157, 159, 163, 167, 175, 189,
223, 257, 269, 271,






197, 209, 211, 213, 223, 251, 257, 263, 265,
277, 293, 299, 339,






269, 271, 277, 291, 293, 299, 309, 311, 313,
341, 343, 353, 355






315, 319, 327, 329, 339, 341, 343, 349,







353, 355






383
33, 35, 47, 53, 61, 83, 89, 91, 95, 99, 105, 111,
359, 367
40
2
0.05



123, 127, 131, 141, 145, 151, 157, 167, 179,







181, 183, 187, 233, 247, 249, 253, 267, 285,







297, 307, 315, 337, 355, 359, 365, 367,







369, 379






479
13, 19, 37, 39, 41, 43, 47, 53, 57, 59, 65, 79, 95,
19, 47, 53, 177, 235,
58
12
0.206897



117, 119, 123, 129, 143, 149, 159, 171, 177,
265, 325, 353, 433,






179, 185, 191, 205, 209, 223, 227, 235, 237,
449, 451, 463






265, 281, 285, 295, 317, 325, 333, 351, 353,







357, 369, 379, 387, 391, 395, 423, 429, 433,







447, 449, 451, 461, 463, 467, 469, 473, 475






503
19, 29, 37, 53, 55, 57, 71, 87, 107, 109, 111,
19, 53, 133, 193, 213,
58
14
0.241379



127, 133, 137, 139, 159, 163, 165, 167, 191,
295, 305, 307, 409,






193, 203, 213, 215, 239, 269, 277, 295, 305,
417, 467, 475, 485,






307, 313, 321, 327, 333, 341, 347, 349, 371,
489






381, 385, 399, 409, 417, 419, 437, 453, 457,







461, 467, 471, 475, 479, 481, 485, 487, 489,







495, 499






587
5, 11, 13, 19, 23, 41, 45, 85, 99, 103, 105, 111,
11, 85, 111, 117, 215,
69
16
0.231884



117, 125, 127, 131, 139, 157, 171, 173, 183,
221, 241, 275, 291,






187, 207, 213, 215, 221, 227, 231, 241, 245,
321, 341, 415, 427,






251, 259, 265, 263, 273, 275, 291, 295, 321,
435, 475, 523






323, 325, 327, 335, 337, 341, 365, 367, 369,







373, 391, 399, 403, 405, 415, 427, 435, 467,







473, 475, 483, 487, 497, 523, 539, 541, 557,







559, 583










Consider the case when a denotes an invertible superprimitive, and let g denote its inverse modulo p. Then, for some integers ν and w, the conditions (10) and (12) take the following forms:

aaν≡a−1(mod p)  (47)

and

g−1≡ggw(mod p).  (48)

Therefore,

aaν+1≡1(mod p)  (49)

whence

aν+1≡0(mod p−1)  (50)

or

a2·ν≡1(mod(p−1))  (51)


Similarly,

g2·w≡1(mod(p−1)) .  (52)

Then

2·ν≡0(mod φ(p−1))  (53)

and

w≡0(mod φ(p−1)).  (54)

NOTE 1: In (53) and (54), v and w cannot be represented as in (19), because, by (17), ρi must be relative prime to φ(p−1). Thus (19) is applicable only if









{





g

c


d

(

v
,

φ

(

p
-
1

)


)


=
1







g

c


d

(

w
,

φ

(

p
-
1

)


)


=
1.








(
55
)








NOTE 2: Compare (56) with (46). The congruence (46) establishes a necessary condition which relates vi(A, G) and wi(A, G). For g≡a−1(mod p), (46) constrains the values of v(A, G) and w(A, G) to be congruent to each other modulo φ(p−1) and also to be congruent to zero modulo φ(p−1).


These properties will be fundamental in the search for a solution of the problem.


APPENDIX A
NOTES ON ORTHOGONAL PRIMITIVES

I. AN EXAMPLE FOR p−1 CONTAINING ONLY SMALL PRIMES


Let

ax≡b(mod 71),  (A.1)

where a and b are primitive roots modulo 71.


Then x is an element of the set X, containing all the integers which are relatively prime to p−1=70=2·5·7.


Let













p
-
1

=


2
·

q
1

·

q
2








=


2
·
5
·
7.








(

A
.2

)







The order of X is φ(70)=φ(5)·φ(7)=24. The exponent of X is e(X)=1 cm(4, 6)=12. Then X can be described as the direct product of a cyclic subgroup of order 2 and a cyclic subgroup of order 12 as follows:

X=C1(2)×C2(12).  (A.3)


Also, the elements of X can be represented by using orthogonal primitives. In this case, given a selection of σ1(mod 7) and σ2(mod 5), ρ1(mod 70) and ρ2(mod 70) can be computed by letting










ρ
1

=


1
+


λ
1

·


p
-
1


q
1




=


σ
1

+


μ
1

·

q
1








(

A
.

4

)








and









ρ
2

=


1
+


λ
2

·


p
-
1


q
2




=


σ
2

+


μ
2

·


q
2

.








(

A
.

5

)







For σ1≡3(mod 5) and σ2≡5(mod 7) , it will be









{





ρ
1



4

3


(

mod


70

)









ρ
2



6

1



(

mod


70

)

.










(

A
.6

)








Then

x≡43x1·61x2(mod70).  (A.7)


Given a pair (a, b), to solve (A.1), observe that, by (A.3) and (A.4),









{







7

0

7

·

ρ
1







7

0

7

·

σ
1

·

(

mod


7

0

)











7

0

7

·

ρ
2






7

0

7




(

mod


70

)

.










(

A
.8

)







Therefore, raising (A.1) to 10 (modulo 71) yields











(

a
10

)


5

x
1







b

1

0


(

mod


71

)

.





(

A
.9

)







Likewise, raising (A.1) to 14 (modulo 71) yields











(

a
14

)


3

x
2







b
14

(

mod


71

)

.





(

A
.10

)








Then (A.7) produces x.


Therefore, in the example, x2 and x1 can be determined independently of each other.

FIG. 1 shows the elements of X as intersections of vertical and horizontal straight lines through 61x2(mod 70) and 43x1(mod 70), respectively.


It is apparent that the elements on a vertical line (constant x2) are congruent to one another modulo 14=2·7. Likewise, the elements on a horizontal line are congruent to one another modulo 2·5=10.


Also, each elements of X is the product of its horizontal and vertical components.


Different selections of the primitives σ1 and σ2 would cause appropriate permutations of the vertical and horizontal lines, respectively.


APPENDIX B
THE ORDER OF σ032 4·ODD+1 MODULO 2ε0

When σ0=4·ODD+1, the order of σ0 modulo 2ε0 equals 2ε0−2:










σ
0

2


ε
0

-
2





1



(

mod



2

ε
0



)

.






(

B
.1

)







Consider the case when σ0=4·ODD+1. Then









{





σ
0





4
·

(


2
·
k

+
1

)


+

1


(

mod



2

ε
0



)












1
+

2
2

+


8
·
k



(

mod



2

ε
0



)










σ
0
2




1
+

2
3

+

16
·


k
1

(

mod



2

ε
0



)










σ
0
4




1
+

2
4

+

32
·


k
2

(

mod



2

ε
0



)










……………………………







σ
0

2


ε
0

-
3






1
+


2


ε
0

-
1




(

mod



2

ε
0



)










σ
0

2


ε
0

-
2






1



(

mod



2

ε
0



)


.










(

B
.2

)








(k, k1, k2 integers).


Notice that the integer σ0,0≡−1+2ε0−1 cannot be produced as a power of σ0.


By using the methods described herein, the encrypted data in the memory is decoded in a time of an order of log log p·log2p by computing the data's discrete logarithm. This speed is important, which only the operation of the second computer performing the second computer generated steps can achieve, because by having this speed for factoring, the signal representative of an encrypted message can be effectively decrypted and deciphered in real time so any threat to property or individuals can be quickly acted upon to eliminate the threat before it occurs and actual damage to property or injury to individuals is prevented or mitigated. In other words, for an encrypted message to be effectively understood, it must by decrypted fast enough that any threat identified in the signal can be stopped. The present invention with the use of the second computer allows for this capability. Here, it is inherent that to save lives if required, the second computer is required.


There may be the step of obtaining the electromagnetic signal representative of a message from a telecommunications network, or a data network or an Internet or a non-transient memory. Law enforcement departments, such as Homeland Security, the FBI, the CIA, NSA, state and local Police or the Military have the well-known capability of obtaining or intercepting messages sent encrypted by a first computer operated by a potential terrorist or criminal as an electromagnetic signal, such as by smart phone or computer or internet, or stored in the memory of a smart phone or computer, or a flash drive. The encrypted electromagnetic signal can be extracted from such messages or memories and operated upon by the techniques described herein to decrypt the encrypted messages and read them to determine whether there is any violation of law or threat to property or individuals. Of course, the intended recipient of the encrypted message by the first computer has the key so the recipient can decrypt the encrypted message the recipient has received and understand it. It is the object of this invention, and the problem this invention solves, to allow a recipient of the encrypted message who does not have the key to read it, to determine what the key is by the techniques described here to compute the data's discrete logarithm, and then using the determined key, decrypting the encrypted message, reviewing what the decrypted message says, and acting as necessary to protect property damage or bodily injury or any type of crime, as deemed appropriate.


Although the invention has been described in detail in the foregoing embodiments for the purpose of illustration, it is to be understood that such detail is solely for that purpose and that variations can be made therein by those skilled in the art without departing from the spirit and scope of the invention except as it may be described by the following claims.


REFERENCES, all of which are incorporated by reference herein.



  • [1] T. M. Apostol, Introduction to Analytic Number Theory, New York, NY: Springer-Verlag, 1976.

  • [2] G. H. Hardy, E. M. Wright, An Introduction to the Theory of Numbers, Oxford, UK: Clarendon Press, 1979.

  • [3] S. C. Pohlig, M. E. Hellman, “An Improved Algorithm for Computing Logarithms over GF(p) and its Cryptographic Significance”, IEEE Trans, Inform. Theory, Vol IT-24, pp. 106-110, 1978.


Claims
  • 1. A method for reading an encrypted message as an electromagnetic signal representative of encrypted data comprising the steps of: entering the encrypted message as the electromagnetic signal into a computer;decrypting by the computer the encrypted message as the electromagnetic signal in a time of the order of log log p x log squared p by selecting with the computer primitives of sub-groups of a group stored in a non-transient memory, where the group is defined modulo ϕ(p−1) in such a way that an exponent of any one primitive is independent on an exponent of any other primitive, thus reducing the complexity of a search for such exponents to a number of operations of the order of a sum of such exponents as opposed to their product to obtain a decrypted message as a decrypted electromagnetic signal;determining if there is any threat in the decrypted message as the decrypted electromagnetic signal; andacting on the decrypted message as the decrypted electromagnetic signal with the computer so any threats to property or individuals can be quickly acted upon to eliminate the threat before it occurs and actual damage to property or injury to individuals is prevented or mitigated.
CROSS-REFERENCE TO RELATED APPLICATIONS

This is a continuation-in-part of U.S. patent application Ser. No. 15/875,737 filed Jan. 19, 2018, now U.S. Pat. No. 10,579,337, which is a continuation-in-part of U.S. patent application Ser. No. 14/886,404 filed Oct. 19, 2015, which claims priority from U.S. provisional application Ser. No. 62/181,322 filed Jun. 18, 2015, all of which are incorporated by reference herein.

US Referenced Citations (10)
Number Name Date Kind
6285761 Patel Sep 2001 B1
6480605 Uchiyama Nov 2002 B1
7031468 Hoffstein Apr 2006 B2
7286665 Wang Oct 2007 B1
20020007457 Neff Jan 2002 A1
20020062330 Paar May 2002 A1
20020186837 Hopkins Dec 2002 A1
20030120944 Kim Jun 2003 A1
20070118746 Lauter May 2007 A1
20100046741 Isogai Feb 2010 A1
Non-Patent Literature Citations (3)
Entry
Extremely fast method for modular exponentiation with modulus and exponent of several million digits, downloaded from http//stackoverlow.com/ questions; Stack Exchange Inc. (Year: 2015).
Ryan Williams et al, Finding orthogonal vectors in discrete structures, Society for Industrial and Applied Mathematics. (Year: 2014).
Machine translation of Japanese patent application No. JP H0548599 A, M. Natsume et al. (Year: 1991).
Provisional Applications (1)
Number Date Country
62181322 Jun 2015 US
Continuation in Parts (2)
Number Date Country
Parent 15875737 Jan 2018 US
Child 16802018 US
Parent 14886404 Oct 2015 US
Child 15875737 US