Mersenne prime

From Wikipedia, the free encyclopedia

In mathematics, a Mersenne number is a number that is one less than a power of two,

Mn = 2n − 1.

A Mersenne prime is a Mersenne number that is a prime number. As of August 2007, only 44 Mersenne primes are known; the largest known prime number (232,582,657−1) is a Mersenne prime and in modern times the largest known prime has nearly always been a Mersenne prime.[1] Like several previous Mersenne primes, it was discovered by a distributed computing project on the Internet, known as the Great Internet Mersenne Prime Search (GIMPS).

Contents

[edit] About Mersenne primes

Many fundamental questions about Mersenne primes remain unresolved. It is not even known whether there is a largest Mersenne prime, which would mean that the set of Mersenne primes is finite. The Lenstra-Pomerance-Wagstaff conjecture asserts that, on the contrary, there are infinitely many Mersenne primes and predicts their order of growth. It is also not known whether infinitely many Mersenne numbers with prime exponents are composite, although this would follow from widely believed conjectures about prime numbers, for example, the infinitude of Sophie Germain primes.

A basic theorem about Mersenne numbers states that in order for Mn to be a Mersenne prime, the exponent n itself must be a prime number. This rules out primality for numbers such as M4 = 24−1 = 15: since the exponent 4=2×2 is composite, the theorem predicts that 15 is also composite; indeed, 15 = 3×5. The three smallest Mersenne primes are

M2 = 3, M3 = 7, M5 = 31.

While it is true that only Mersenne numbers Mp, where p = 2, 3, 5, … could be prime, it may nevertheless turn out that Mp is not prime even for a prime exponent p. The smallest counterexample is the Mersenne number

M11 = 211 − 1 = 2047 = 23 × 89,

which is not a Mersenne prime, even though 11 is a prime number. The lack of an obvious rule to determine whether a given Mersenne number is prime makes the search for Mersenne primes an interesting task, which becomes difficult very soon, since Mersenne numbers grow very fast. The Lucas–Lehmer test for Mersenne numbers is an efficient primality test that greatly aids this task. Search for the largest known prime has somewhat of a cult following. Consequently, a lot of computer power has been expended searching for new Mersenne primes, much of which is now done using distributed computing.

Mersenne primes are used in pseudorandom number generators such as Mersenne Twister and Park–Miller RNG.

[edit] Searching for Mersenne primes

The identity

2^{ab}-1=(2^a-1)\cdot \left(1+2^a+2^{2a}+2^{3a}+\dots+2^{(b-1)a}\right)

shows that Mn can be prime only if n itself is prime—that is the primeness of n is necessary but not sufficient for Mn to be prime—which simplifies the search for Mersenne primes considerably. (This follows very simply from the Mersenne property of the sequence of numbers of the form xn − yn. This states that xa − ya | xb − yb if and only if a|b.) The converse statement, namely that Mn is necessarily prime if n is prime, is false. The smallest counterexample is 211−1 = 23×89, a composite number.

Fast algorithms for finding Mersenne primes are available, and the largest known prime numbers as of 2007 are Mersenne primes.

The first four Mersenne primes M2 = 3, M3 = 7, M5 = 31 and M7 = 127 were known in antiquity. The fifth, M13 = 8191, was discovered anonymously before 1461; the next two (M17 and M19) were found by Cataldi in 1588. After nearly two centuries, M31 was verified to be prime by Euler in 1772. The next (in historical, not numerical order) was M127, found by Lucas in 1876, then M61 by Pervushin in 1883. Two more (M89 and M107) were found early in the 20th century, by Powers in 1911 and 1914, respectively.

The best method presently known for testing the primality of Mersenne numbers is based on the computation of a recurring sequence, as developed originally by Lucas in 1856 [1][2] and improved by Lehmer in the 1930s, now known as the Lucas-Lehmer test for Mersenne numbers. Specifically, it can be shown that (for n > 2) Mn = 2n − 1 is prime if and only if Mn divides Sn−2, where S0 = 4 and for k > 0, S_k=S_{k-1}^2-2.

Graph of number of digits in largest known Mersenne prime by year - electronic era. Note that the vertical scale is logarithmic.
Graph of number of digits in largest known Mersenne prime by year - electronic era. Note that the vertical scale is logarithmic.

The search for Mersenne primes was revolutionized by the introduction of the electronic digital computer. The first successful identification of a Mersenne prime, M521, by this means was achieved at 10:00 P.M. on January 30, 1952 using the U.S. National Bureau of Standards Western Automatic Computer (SWAC) at the Institute for Numerical Analysis at the University of California, Los Angeles, under the direction of Lehmer, with a computer search program written and run by Prof. R.M. Robinson. It was the first Mersenne prime to be identified in thirty-eight years; the next one, M607, was found by the computer a little less than two hours later. Three more — M1279, M2203, M2281 — were found by the same program in the next several months. M4253 is the first Mersenne prime that is titanic, M44497 is the first gigantic, and M6,972,593 was the first megaprime to be discovered, being a prime with at least 1,000,000 digits.[3] All three were the first known prime of any kind of that size.

[edit] Theorems about Mersenne numbers

c^n-d^n=(c-d)\sum_{k=0}^{n-1} c^kd^{n-1-k},

or

(2^a-1)\cdot \left(1+2^a+2^{2a}+2^{3a}+\dots+2^{(b-1)a}\right)=2^{ab}-1

by setting c = 2a, d = 1, and n = b

proof

(a-b)\sum_{k=0}^{n-1}a^kb^{n-1-k}
=\sum_{k=0}^{n-1}a^{k+1}b^{n-1-k}-\sum_{k=0}^{n-1}a^kb^{n-k}
=a^n+\sum_{k=1}^{n-1}a^kb^{n-k}-\sum_{k=1}^{n-1}a^kb^{n-k}-b^n
= an − bn
  • 2) If 2n − 1 is prime, then n is prime.

proof

By

(2^a-1)\cdot \left(1+2^a+2^{2a}+2^{3a}+\dots+2^{(b-1)a}\right)=2^{ab}-1

If n is not prime, or n = ab where 1 < a,b < n. Therefore, 2a − 1 would divide 2n − 1, or 2n − 1 is not prime.

  • 3) If p is an odd prime, then any prime q that divides 2p − 1 must be 1 plus a multiple of 2p. This holds even when 2p − 1 is prime. Example I: 25 − 1 = 31

is prime, and 31 is 1 plus a multiple of 2*5. Example II: 211 − 1=23*89, 23=1+2*11, and 89=1+8*11, and also 23*89=1+186*11.

proof

If q divides 2p − 1 then 2p is congruent to 1 mod q, so p divides the order of the multiplicative group mod q, by Lagrange's Theorem. This group has order q−1, so q−1=kp for some k, and q=1+kp. But q must be odd, and p is odd,(except for p=2) so k is even.

  • 4) If p is an odd prime, then any prime q that divides 2p − 1 must be congruent to \pm 1 \pmod 8. Proof: 2p + 1 = 2(mod q), so 2(p + 1) / 2 is a square root of 2 modulo q. By quadratic reciprocity, any prime modulo which two has a square root is congruent to \pm 1 \pmod 8.

[edit] History

Mersenne primes were considered already by Euclid, who found a connection with the perfect numbers. They are named after 17th century French scholar Marin Mersenne, who compiled a list of Mersenne primes with exponents up to 257. His list was only partially correct, as Mersenne mistakenly included M67 and M257 (which are composite), and omitted M61, M89, and M107 (which are prime). Mersenne gave no indication how he came up with his list, and its rigorous verification was completed more than two centuries later.

[edit] List of known Mersenne primes

The table below lists all known Mersenne primes (sequence A000668 in OEIS):

# n Mn Digits in Mn Date of discovery Discoverer
1 2 3 1 ancient ancient
2 3 7 1 ancient ancient
3 5 31 2 ancient ancient
4 7 127 3 ancient ancient
5 13 8191 4 1456 anonymous [4]
6 17 131071 6 1588 Cataldi
7 19 524287 6 1588 Cataldi
8 31 2147483647 10 1772 Euler
9 61 2305843009213693951 19 1883 Pervushin
10 89 618970019…449562111 27 1911 Powers
11 107 162259276…010288127 33 1914 Powers[5]
12 127 170141183…884105727 39 1876 Lucas
13 521 686479766…115057151 157 January 30, 1952 Robinson
14 607 531137992…031728127 183 January 30, 1952 Robinson
15 1,279 104079321…168729087 386 June 25, 1952 Robinson
16 2,203 147597991…697771007 664 October 7, 1952 Robinson
17 2,281 446087557…132836351 687 October 9, 1952 Robinson
18 3,217 259117086…909315071 969 September 8, 1957 Riesel
19 4,253 190797007…350484991 1,281 November 3, 1961 Hurwitz
20 4,423 285542542…608580607 1,332 November 3, 1961 Hurwitz
21 9,689 478220278…225754111 2,917 May 11, 1963 Gillies
22 9,941 346088282…789463551 2,993 May 16, 1963 Gillies
23 11,213 281411201…696392191 3,376 June 2, 1963 Gillies
24 19,937 431542479…968041471 6,002 March 4, 1971 Tuckerman
25 21,701 448679166…511882751 6,533 October 30, 1978 Noll & Nickel
26 23,209 402874115…779264511 6,987 February 9, 1979 Noll
27 44,497 854509824…011228671 13,395 April 8, 1979 Nelson & Slowinski
28 86,243 536927995…433438207 25,962 September 25, 1982 Slowinski
29 110,503 521928313…465515007 33,265 January 28, 1988 Colquitt & Welsh
30 132,049 512740276…730061311 39,751 September 19, 1983[6] Slowinski
31 216,091 746093103…815528447 65,050 September 1, 1985[7] Slowinski
32 756,839 174135906…544677887 227,832 February 19, 1992 Slowinski & Gage on Harwell Lab Cray-2 [8]
33 859,433 129498125…500142591 258,716 January 4, 1994 [9] Slowinski & Gage
34 1,257,787 412245773…089366527 378,632 September 3, 1996 Slowinski & Gage [10]
35 1,398,269 814717564…451315711 420,921 November 13, 1996 GIMPS / Joel Armengaud [11]
36 2,976,221 623340076…729201151 895,932 August 24, 1997 GIMPS / Gordon Spence [12]
37 3,021,377 127411683…024694271 909,526 January 27, 1998 GIMPS / Roland Clarkson [13]
38 6,972,593 437075744…924193791 2,098,960 June 1, 1999 GIMPS / Nayan Hajratwala [14]
39 13,466,917 924947738…256259071 4,053,946 November 14, 2001 GIMPS / Michael Cameron [15]
40* 20,996,011 125976895…855682047 6,320,430 November 17, 2003 GIMPS / Michael Shafer [16]
41* 24,036,583 299410429…733969407 7,235,733 May 15, 2004 GIMPS / Josh Findley [17]
42* 25,964,951 122164630…577077247 7,816,230 February 18, 2005 GIMPS / Martin Nowak [18]
43* 30,402,457 315416475…652943871 9,152,052 December 15, 2005 GIMPS / Curtis Cooper & Steven Boone [19]
44* 32,582,657 124575026…053967871 9,808,358 September 4, 2006 GIMPS / Curtis Cooper & Steven Boone [20]

*It is not known whether any undiscovered Mersenne primes exist between the 39th (M13,466,917) and the 44th (M32,582,657) on this chart; the ranking is therefore provisional.

To help visualize the size of the 44th known Mersenne prime, it would require 2,616 pages to display the number in base 10 with 75 digits per line and 50 lines per page.

[edit] Factorization of Mersenne numbers

Mersenne numbers are very good test cases for the special number field sieve algorithm, so often the largest number factorised has been a Mersenne number. As of March 2007, 21039 − 1 is the record-holder, after a calculation taking about a year on a couple of hundred computers, mostly at NTT in Japan and at EPFL in Switzerland. See integer factorization records for links to more information.

[edit] Perfect numbers

Mersenne primes are interesting to many for their connection to perfect numbers. In the 4th century BC, Euclid demonstrated that if Mp is a Mersenne prime then

2p−1×(2p−1) = Mp(Mp+1)/2

is an even perfect number. In the 18th century, Leonhard Euler proved that, conversely, all even perfect numbers have this form. It is unknown whether there are any odd perfect numbers.

[edit] Generalization

The binary representation of 2n − 1 is the digit 1 repeated n times, for example, 25 − 1 = 111112 in the binary notation. The Mersenne primes are therefore the base-2 repunit primes.

[edit] Notes

  1. ^ The largest prime has been a Mersenne prime since 1952, except between 1989 and 1992; see Caldwell, "The Largest Known Prime by Year: A Brief History", from The Prime Pages website, U of Tennessee at Martin.

[edit] See also

[edit] External links

[edit] Math World Links