In my previous article, I talked about the RSA cryptosystem which is widely used on the Internet for secure data transmission. The power and security of the RSA cryptosystem derives from the fact that the factoring problem is “hard.” That is, it is believed that the full decryption of an RSA ciphertext is infeasible because no efficient classical algorithm currently exists for factoring large numbers. However, in 1994 Peter Shor showed that a quantum computer could be used to factor a number in polynomial time, thus effectively breaking RSA.
It may be tempting to use the speed of a quantum computer to simply check all possible divisors in parallel. In this case, we would be performing a classical algorithm on a quantum computer, making use only of the increased speed of the quantum machine. Unfortunately, this is not going to work. In a way, it is possible for a quantum computer to try all possible divisors. However, due to the nature of quantum computing, when measuring the outcome of the computations, you will get a random possible divisor, which is almost certainly not the one you want.
How, then, can we use a quantum computer to solve the factoring problem? The key to a fast and accurate quantum factoring algorithm is to make use of the structure of the factoring problem itself. Instead of looking for factors directly, we must use some mathematical property of factoring. Fortunately, the factoring problem has plenty of special properties from which to choose. For example, given a positive integer, even if we do not know its prime factorization we do know that it has exactly one factorization. This fact does not help us solve the factorization problem, but it does give us hope that the problem has other nice mathematical properties that will.
The property we will use is the ability to reduce the prime factorization problem into a problem of order (or period) finding. Let’s start by looking at an example. Consider the sequence of numbers
2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, …
Now, let’s look at this same sequence of powers of two, but taken “mod 15.” In other words, we will create a new sequence of numbers consisting of the remainders when each power of two is divided by 15. This gives us the new sequence
2, 4, 8, 1, 2, 4, 8, 1, 2, 4, …
We see that taking the powers of two mod 15 gives us a periodic sequence whose period (or order) is four. For another example, consider the same powers of two, but taken mod 21. In this case, we have the new sequence
2, 4, 8, 16, 11, 1, 2, 4, 8, 16, …
Here, we get a periodic sequence whose period is six.
In the 1760s, Euler discovered a beautiful pattern to this period finding problem. Let be the product of two prime numbers, and . Consider the sequence
…
If is not divisible by or , then the above sequence will repeat with some period that evenly divides .
In our examples above, we have . In the first example, we have which has the prime factors and . Then, , which is divisible by the period of 4. In the second example, we have which has the prime factors and . Then, , which is divisible by the period of 6.
But how does this help us solve the factorization problem? If we can find the period of the sequence
…
then we learn something about the prime factors of . In particular, we learn a divisor of . Of course, we’d rather learn the factors and themselves, but this, at least, represents progress. If we determine several random divisors of by trying different random values of , then we can multiply those divisors together to obtain itself. Once we know , we can then determine and .
However, there’s still a problem with applying our observations to the factorization problem. Even though the sequence
…
will eventually start repeating itself, the number of steps before it repeats could be almost as large as , which in the RSA cryptosystem is a very large number! This issue is why finding the period of the sequence does not appear to lead to a classical factoring algorithm. However, with the help of quantum mechanics, we can define a quantum algorithm that works in a reasonable amount of time.
Shor’s algorithm is composed of two parts. The first part turns the factoring problem into the period finding problem, and can be computed on a classical computer. The second part (step 2 below) finds the period using the quantum Fourier transform and is responsible for the quantum speedup of the algorithm. We begin by briefly describing all five steps. After that, we will focus on the quantum part of the algorithm (i.e. step 2). To factor a large integer (which, without loss of generality, we may assume is odd), we use Shor’s algorithm:
1. Choose a random positive integer . Compute gcd, which may be done in polynomial time using the Euclidean algorithm. If gcd, then we have found a non-trivial factor of , and we are done. If, on the other hand, gcd, proceed to step 2.
2. Use a quantum computer to determine the unknown period of the sequence
…
3. If is an odd integer, then return to step 1. The probability that is odd is , where is the number of distinct prime factors of . If is even, then proceed to step 4.
4. Since is even,
If , then go to step 1. If , then proceed to step 5. It can be shown that the probability that is less than , where denotes the number of distinct prime factors of .
5. Compute gcd using the Euclidean algorithm. Since , it can be shown that is a non-trivial factor of . Exit with the answer .
Thus, the task of factoring an odd positive integer reduces to the problem of finding the period of a function/sequence. Shor’s period-finding algorithm (step 2 above) relies heavily on the ability of a quantum computer to be in many states simultaneously (a superposition of states). To compute the period of a function , we evaluate the function at all points simultaneously.
Unfortunately, quantum mechanics does not allow us to access all of this information directly. Instead, a measurement must be taken which will yield only one of the possible values (destroying all others). Because of this issue, we must transform the superposition to another state that will return the correct period with high probability. This is achieved by using the quantum Fourier transform. The main components of Shor’s period-finding algorithm are as follows:
1. Create a superposition of states. This can be done by applying Hadamard gates to all qubits in the input register.
2. Implement the function as a quantum transform. To achieve this, Shor used repeated squaring for his modular exponentiation transform.
3. Perform a quantum Fourier transform.
After these transformations, a measurement will yield an approximation to the period .
Now let’s look at an example of how can be factored using Shor’s algorithm.
Step 1. Choose a random positive integer , say . Since gcd, proceed to step 2 to find the period of the function given by
.
Step 2. Run Shor’s period-finding algorithm on a quantum computer to find (with high probability) that the period .
Step 3. Since is even, we proceed to step 4.
Step 4. Since
proceed to step 5.
Step 5. With the Euclidean algorithm, compute
gcd= gcd= gcd
We have succeeded in using Shor’s algorithm to find a non-trivial factor of , namely .
Sources:
Shor, P. W. (1997). Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer. SIAM Journal on Computing, 26(5), 1484-26. doi:http://dx.doi.org/10.1137/S0097539795293172
Lomonaco, Jr, S. J. (2000) Shor’s Quantum Factoring Algorithm. arXiv:quant-ph/0010034
Thanks very much for this. I found it very helpful for getting started on thinking about what one might actually do with a quantum computer.
I think there is an error, just a typo, in one of your equations. Where it says:
(m^(P/2)-1)(m^(P/2)-1) = m^P – 1
I think it should say:
(m^(P/2)-1)(m^(P/2)+1) = m^P – 1
I would be 100% sure about this, except that the topic is number theory, which often pulls the rug out from under my physicist brain. For example, the number 1 has entirely too many square roots in number theory, where the modulo n part is sometimes sort of an afterthought.
I think most people will just read around that error and either not notice, or else just know what you meant, although like I said, I am still not 100% sure. If it’s actually correct as is, you would help people like me by explaining how.
Again, thanks much,
Bill Peria
Hello, Bill.
I see the same typo in the equation.
I think you got this a little wrong. To begin with, Bill Peria is correct. You factored m^p – 1 incorrectly.
The next error is in step 4. You don’t want the = sign there, you want the ≡ congruence symbol.
m^P – 1 ≡ 0 Mod N
But if the next relation
m^(P/2) + 1 ≡ 0 Mod N
is true, stop. You found a factor of N. Why would you go back to step 1 at this point?