Probable prime

Not to be confused with Provable prime.

In number theory, a probable prime (PRP) is an integer that satisfies a specific condition that is satisfied by all prime numbers, but which is not satisfied by most composite numbers. Different types of probable primes have different specific conditions. While there may be probable primes that are composite (called pseudoprimes), the condition is generally chosen in order to make such exceptions rare.

Fermat's test for compositeness, which is based on Fermat's little theorem, works as follows: given an integer n, choose some integer a that is not a multiple of n; (typically, we choose a in the range 1 < a < n − 1). Calculate a^{n− 1} modulo n. If the result is not 1, then n is composite. If the result is 1, then n is likely to be prime; n is then called a probable prime to basea. A weak probable prime to basea is an integer that is a probable prime to base a, but which is not a strong probable prime to base a (see below).

For a fixed base a, it is unusual for a composite number to be a probable prime (that is, a pseudoprime) to that base. For example, up to 25 × 10⁹, there are 11,408,012,595 odd composite numbers, but only 21,853 pseudoprimes base 2. ^{[1] : 1005} The number of odd primes in the same interval is 1,091,987,404.

Properties

Probable primality is a basis for efficient primality testing algorithms, which find application in cryptography. These algorithms are usually probabilistic in nature. The idea is that while there are composite probable primes to base a for any fixed a, we may hope there exists some fixed P<1 such that for any given composite n, if we choose a at random, then the probability that n is pseudoprime to base a is at most P. If we repeat this test k times, choosing a new a each time, the probability of n being pseudoprime to all the as tested is hence at most P^k, and as this decreases exponentially, only moderate k is required to make this probability negligibly small (compared to, for example, the probability of computer hardware error).

This is unfortunately false for weak probable primes, because there exist Carmichael numbers; but it is true for more refined notions of probable primality, such as strong probable primes (P = 1/4, Miller–Rabin algorithm), or Euler probable primes (P = 1/2, Solovay–Strassen algorithm).

Even when a deterministic primality proof is required, a useful first step is to test for probable primality. This can quickly eliminate (with certainty) most composites.

A PRP test is sometimes combined with a table of small pseudoprimes to quickly establish the primality of a given number smaller than some threshold.

Variations

An Euler probable prime to basea is an integer that is indicated prime by the somewhat stronger theorem that for any prime p, a^(p−1)/2 equals ${\displaystyle ({\tfrac {a}{p}})}$ modulo p, where ${\displaystyle ({\tfrac {a}{p}})}$ is the Jacobi symbol. An Euler probable prime which is composite is called an Euler–Jacobi pseudoprime to base a. The smallest Euler-Jacobi pseudoprime to base 2 is 561. ^{[1] : 1004} There are 11347 Euler-Jacobi pseudoprimes base 2 that are less than 25·10⁹. ^{[1] : 1005}

This test may be improved by using the fact that the only square roots of 1 modulo a prime are 1 and −1. Write n = d · 2^s + 1, where d is odd. The number n is a strong probable prime (SPRP) to basea if:

${\displaystyle a^{d}\equiv 1{\pmod {n}},\;}$

or

${\displaystyle a^{d\cdot 2^{r}}\equiv -1{\pmod {n}}{\text{ for some }}0\leq r\leq s-1.\,}$

A composite strong probable prime to base a is called a strong pseudoprime to base a. Every strong probable prime to base a is also an Euler probable prime to the same base, but not vice versa.

The smallest strong pseudoprime base 2 is 2047. ^{[1] : 1004} There are 4842 strong pseudoprimes base 2 that are less than 25·10⁹. ^{[1] : 1005}

There are also Lucas probable primes, which are based on Lucas sequences. A Lucas probable prime test can be used alone. The Baillie–PSW primality test combines a Lucas test with a strong probable prime test.

Example of testing for a strong probably prime

To test whether 97 is a strong probable prime base 2:

• Step 1: Find ${\displaystyle d}$ and ${\displaystyle s}$ for which ${\displaystyle 96=d\cdot 2^{s}}$, where ${\displaystyle d}$ is odd
  • Beginning with ${\displaystyle s=0}$, ${\displaystyle d}$ would be ${\displaystyle 96}$
  • Increasing ${\displaystyle s}$, we see that ${\displaystyle d=3}$ and ${\displaystyle s=5}$, since ${\displaystyle 96=3\cdot 2^{5}}$
• Step 2: Choose ${\displaystyle a}$, ${\displaystyle 1<a<97-1}$. We will choose ${\displaystyle a=2}$.
• Step 3: Calculate ${\displaystyle a^{d}{\bmod {n}}}$, i.e. ${\displaystyle 2^{3}{\bmod {9}}7}$. Since it isn't congruent to ${\displaystyle 1}$, we continue to test the next condition
• Step 4: Calculate ${\displaystyle 2^{3\cdot 2^{r}}{\bmod {9}}7}$ for ${\displaystyle 0\leq r<s}$. If it is congruent to ${\displaystyle 96}$, ${\displaystyle 97}$ is probably prime. Otherwise, ${\displaystyle 97}$ is definitely composite
  • ${\displaystyle r=0:2^{3}\equiv 8{\pmod {97}}}$
  • ${\displaystyle r=1:2^{6}\equiv 64{\pmod {97}}}$
  • ${\displaystyle r=2:2^{12}\equiv 22{\pmod {97}}}$
  • ${\displaystyle r=3:2^{24}\equiv 96{\pmod {97}}}$
• Therefore, ${\displaystyle 97}$ is a strong probable prime base 2 (and is therefore a probable prime base 2).

See also

• Provable prime
• Baillie–PSW primality test
• Euler–Jacobi pseudoprime
• Lucas pseudoprime
• Miller–Rabin primality test
• Perrin primality test
• Carmichael number

External links

• The prime glossary – Probable prime
• The PRP Top 10000 (the largest known probable primes)

References

1. Carl Pomerance; John L. Selfridge; Samuel S. Wagstaff, Jr. (July 1980). "The pseudoprimes to 25·10⁹" (PDF). Mathematics of Computation. 35 (151): 1003–1026. doi: 10.1090/S0025-5718-1980-0572872-7 . JSTOR   2006210.