Solovay–Strassen primality test

The Solovay–Strassen primality test, developed by Robert M. Solovay and Volker Strassen, is a probabilistic test to determine if a number is composite or probably prime. It has been largely superseded by the Baillie-PSW primality test and the Miller–Rabin primality test, but has great historical importance in showing the practical feasibility of the RSA cryptosystem. The Solovay–Strassen test is essentially a Euler-Jacobi pseudoprime test.

Concepts

Euler proved^[1] that for any prime number p and any integer a,

a^{(p-1)/2}\equiv \left({\frac {a}{p}}\right){\pmod {p}}

where $\left({\tfrac {a}{p}}\right)$ is the Legendre symbol. The Jacobi symbol is a generalisation of the Legendre symbol to $\left({\tfrac {a}{n}}\right)$ , where n can be any odd integer. The Jacobi symbol can be computed in time O((log n)²) using Jacobi's generalization of law of quadratic reciprocity.

Given an odd number n we can contemplate whether or not the congruence

a^{(n-1)/2}\equiv \left({\frac {a}{n}}\right){\pmod {n}}

holds for various values of the "base" a, given that a is relatively prime to n. If n is prime then this congruence is true for all a. So if we pick values of a at random and test the congruence, then as soon as we find an a which doesn't fit the congruence we know that n is not prime (but this does not tell us a nontrivial factorization of n). This base a is called an Euler witness for n; it is a witness for the compositeness of n. The base a is called an Euler liar for n if the congruence is true while n is composite.

For every composite odd n at least half of all bases

a\in (\mathbb {Z} /n\mathbb {Z} )^{*}

are (Euler) witnesses:^[2] this contrasts with the Fermat primality test, for which the proportion of witnesses may be much smaller. Therefore, there are no (odd) composite n without many witnesses, unlike the case of Carmichael numbers for Fermat's test.

Example

Suppose we wish to determine if n = 221 is prime. We write (n−1)/2=110.

We randomly select an a (smaller than n): 47. Using an efficient method for raising a number to a power (mod n) such as binary exponentiation, we compute:

a^(n−1)/2 mod n = 47¹¹⁰ mod 221 = −1 mod 221
$({\tfrac {a}{n}})$ mod n = $({\tfrac {47}{221}})$ mod 221 = −1 mod 221.

This gives that, either 221 is prime, or 47 is an Euler liar for 221. We try another random a, this time choosing a = 2:

a^(n−1)/2 mod n = 2¹¹⁰ mod 221 = 30 mod 221
$({\tfrac {a}{n}})$ mod n = $({\tfrac {2}{221}})$ mod 221 = −1 mod 221.

Hence 2 is an Euler witness for the compositeness of 221, and 47 was in fact an Euler liar. Note that this tells us nothing about the prime factors of 221, which are actually 13 and 17.

Algorithm and running time

The algorithm can be written in pseudocode as follows:

Inputs: n, a value to test for primality; k, a parameter that determines the accuracy of the test
Output: composite if n is composite, otherwise probably prime

repeat k times:
   choose a randomly in the range [2,n − 1]
    $x\gets \left({\tfrac {a}{n}}\right)$ 
   if  $x = 0$  or  $a^{(n-1)/2}\not \equiv x{\pmod {n}}$  then 
      return composite
return probably prime

Using fast algorithms for modular exponentiation, the running time of this algorithm is O(k·log³ n), where k is the number of different values of a we test.

Accuracy of the test

It is possible for the algorithm to return an incorrect answer. If the input n is indeed prime, then the output will always correctly be probably prime. However, if the input n is composite then it is possible for the output to be incorrectly probably prime. The number n is then called a Euler-Jacobi pseudoprime.

When n is odd and composite, at least half of all a with gcd(a,n) = 1 are Euler witnesses. We can prove this as follows: let {a₁, a₂, ..., a_m} be the Euler liars and a an Euler witness. Then, for i = 1,2,...,m:

(a\cdot a_{i})^{(n-1)/2}=a^{(n-1)/2}\cdot a_{i}^{(n-1)/2}=a^{(n-1)/2}\cdot \left({\frac {a_{i}}{n}}\right)\not \equiv \left({\frac {a}{n}}\right)\left({\frac {a_{i}}{n}}\right){\pmod {n}}.

Because the following holds:

\left({\frac {a}{n}}\right)\left({\frac {a_{i}}{n}}\right)=\left({\frac {a\cdot a_{i}}{n}}\right),

now we know that

(a\cdot a_{i})^{(n-1)/2}\not \equiv \left({\frac {a\cdot a_{i}}{n}}\right){\pmod {n}}.

This gives that each a_i gives a number a·a_i, which is also an Euler witness. So each Euler liar gives an Euler witness and so the number of Euler witnesses is larger or equal to the number of Euler liars. Therefore, when n is composite, at least half of all a with gcd(a,n) = 1 is an Euler witness.

Hence, the probability of failure is at most 2^−k (compare this with the probability of failure for the Miller-Rabin primality test, which is at most 4^−k).

For purposes of cryptography the more bases a we test, i.e. if we pick a sufficiently large value of k, the better the accuracy of test. Hence the chance of the algorithm failing in this way is so small that the (pseudo) prime is used in practice in cryptographic applications, but for applications for which it is important to have a prime, a test like ECPP or Pocklington^[3] should be used which proves primality.

Average-case behaviour

The bound 1/2 on the error probability of a single round of the Solovay–Strassen test holds for any input n, but those numbers n for which the bound is (approximately) attained are extremely rare. On the average, the error probability of the algorithm is significantly smaller: it is less than

2^{-k}\exp \left(-(1+o(1)){\frac {\log x\,\log \log \log x}{\log \log x}}\right)

for k rounds of the test, applied to uniformly random n ≤ x.^[4]^[5] The same bound also applies to the related problem of what is the conditional probability of n being composite for a random number n ≤ x which has been declared prime in k rounds of the test.

Complexity

The Solovay–Strassen algorithm shows that the decision problem COMPOSITE is in the complexity class RP.^[6]

References

↑ Euler's criterion
↑ PlanetMath
↑ Pocklington test on Mathworld
↑ P. Erdős; C. Pomerance (1986). "On the number of false witnesses for a composite number". Mathematics of Computation. American Mathematical Society. 46 (173): 259–279. JSTOR 2008231. doi:10.2307/2008231.
↑ I. Damgård; P. Landrock; C. Pomerance (1993). "Average case error estimates for the strong probable prime test". Mathematics of Computation. American Mathematical Society. 61 (203): 177–194. JSTOR 2152945. doi:10.2307/2152945.
↑ R. Motwani; P. Raghavan (1995). Randomized Algorithms. Cambridge University Press. pp. 417–423. ISBN 0-521-47465-5.

External links

Solovay-Strassen Implementation of the Solovay–Strassen primality test in Maple

Number-theoretic algorithms
Primality tests	AKS test APR test Baillie–PSW Elliptic curve Pocklington Fermat Lucas Lucas–Lehmer Lucas–Lehmer–Riesel Proth's theorem Pépin's Quadratic Frobenius test Solovay–Strassen Miller–Rabin
Prime-generating	Sieve of Atkin Sieve of Eratosthenes Sieve of Sundaram Wheel factorization
Integer factorization	Continued fraction (CFRAC) Dixon's Lenstra elliptic curve (ECM) Euler's Pollard's rho p − 1 p + 1 Quadratic sieve (QS) General number field sieve (GNFS) Special number field sieve (SNFS) Rational sieve Fermat's Shanks' square forms Trial division Shor's
Multiplication	Ancient Egyptian Long Karatsuba Toom–Cook Schönhage–Strassen Fürer's
Discrete logarithm	Baby-step giant-step Pollard rho Pollard kangaroo Pohlig–Hellman Index calculus Function field sieve
Greatest common divisor	Binary Euclidean Extended Euclidean Lehmer's
Modular square root	Cipolla Pocklington's Tonelli–Shanks
Other algorithms	Chakravala Cornacchia LLL Integer square root Modular exponentiation Schoof's
Italics indicate that algorithm is for numbers of special forms

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.