| Note: If you do not see superscripts
and subscripts in this sentence,
then click here for a version of this page without super/subscripts. |
Enter French monk Marin Mersenne (1588-1648). Mersenne stated in the preface to his Cogitata Physica-Mathematica (1644) that the numbers 2n-1 were prime for
n = 2, 3, 5, 7, 13, 17, 19, 31, 67, 127 and 257and were composite for all other positive integers n < 257. Mersenne's (incorrect) conjecture fared only slightly better than Regius', but still got his name attached to these numbers.
Definition: When 2n-1 is prime it is said to be a Mersenne prime.It was obvious to Mersenne's peers that he could not have tested all of these numbers (in fact he admitted as much), but they could not test them either. It was not until over 100 years later, in 1750, that Euler verified the next number on Mersenne's and Regius' lists, 231-1, was prime. After another century, in 1876, Lucas verified 2127-1 was also prime. Seven years later Pervouchine showed 261-1 was prime, so Mersenne had missed this one. In the early 1900's Powers showed that Mersenne had also missed the primes 289-1 and 2107-1. Finally, by 1947 Mersenne's range, n < 258, had been completely checked and it was determined that the correct list is:
n = 2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107 and 127.See the table of known Mersenne primes below.
Definition: A positive integer n is called a perfect number if it is equal to the sum of all of its positive divisors, excluding n itself.For example, 6 is the first perfect number because 6=1+2+3. The next is 28=1+2+4+7+14. The next two are 496 and 8128. These four were all known before the time of Christ. Look at these numbers in the following partially factored form:
2.3, 4.7, 16.31, 64.127.Do you notice they all have the same form 2n-1(2n-1) (for n = 2, 3, 5, and 7 respectively)? And that in each case 2n-1 was a Mersenne prime? In fact it is easy to show the following theorems:
Theorem One: k is an even perfect number if and only if it has the form 2n-1(2n-1) and 2n-1 is prime. [Proof.]So the search for Mersennes is also the search for even perfect numbers!Theorem Two: If 2n-1 is prime, then so is n. [Proof.]
You may have also noticed that the perfect numbers listed above (6, 28, 496, 8128) all end with either the digit 6 or the digit 8--this is also very easy to prove (but no, they do not continue to alternate 6, 8, 6, 8,...). If you like that digit pattern, look at the first four perfect numbers in binary:
(The binary digit pattern is a consequence of Theorem One.) It is not known whether or not there is an odd perfect number, but if there is one it is big!110 11100 111110000 1111111000000
When checking to see if a Mersenne number is prime, we usually first look for any small divisors. The following theorem of Euler and Fermat is very useful in this regard.
Theorem Three: Let p and q be primes. If q divides Mp = 2p-1, thenFinally, we offer the following for your perusal:q = +/-1 (mod 8) and q = 2kp + 1for some integer k. [Proof.]
Theorem Four. Let p = 3 (mod 4) be prime. 2p+1 is also prime if and only if 2p+1 divides Mp. [Proof].
Theorem Five. If you sum the digits of any even perfect number (except 6), then sum the digits of the resulting number, and repeat this process until you get a single digit, that digit will be one. [Proof.]
| ## | p (exponent) | digits in Mp | digits in Pp |
year | discoverer | notes |
|---|---|---|---|---|---|---|
| 1 | 2 | 1 | 1 | ---- | ---- | |
| 2 | 3 | 1 | 2 | ---- | ---- | |
| 3 | 5 | 2 | 3 | ---- | ---- | |
| 4 | 7 | 3 | 4 | ---- | ---- | |
| 5 | 13 | 4 | 8 | 1456 | anonymous | |
| 6 | 17 | 6 | 10 | 1588 | Cataldi | |
| 7 | 19 | 6 | 12 | 1588 | Cataldi | |
| 8 | 31 | 10 | 19 | 1772 | Euler | |
| 9 | 61 | 19 | 37 | 1883 | Pervushin | |
| 10 | 89 | 27 | 54 | 1911 | Powers | |
| 11 | 107 | 33 | 65 | 1914 | Powers | note |
| 12 | 127 | 39 | 77 | 1876 | Lucas | |
| 13 | 521 | 157 | 314 | 1952 | Robinson | |
| 14 | 607 | 183 | 366 | 1952 | Robinson | |
| 15 | 1279 | 386 | 770 | 1952 | Robinson | |
| 16 | 2203 | 664 | 1327 | 1952 | Robinson | |
| 17 | 2281 | 687 | 1373 | 1952 | Robinson | |
| 18 | 3217 | 969 | 1937 | 1957 | Riesel | |
| 19 | 4253 | 1281 | 2561 | 1961 | Hurwitz | |
| 20 | 4423 | 1332 | 2663 | 1961 | Hurwitz | |
| 21 | 9689 | 2917 | 5834 | 1963 | Gillies | |
| 22 | 9941 | 2993 | 5985 | 1963 | Gillies | |
| 23 | 11213 | 3376 | 6751 | 1963 | Gillies | |
| 24 | 19937 | 6002 | 12003 | 1971 | Tuckerman | [Tuckerman71] |
| 25 | 21701 | 6533 | 13066 | 1978 | Noll & Nickel | [NN80] |
| 26 | 23209 | 6987 | 13973 | 1979 | Noll | " |
| 27 | 44497 | 13395 | 26790 | 1979 | Nelson & Slowinski | [Slowinski79] |
| 28 | 86243 | 25962 | 51924 | 1982 | Slowinski | [Ewing83] |
| 29 | 110503 | 33265 | 66530 | 1988 | Colquitt & Welsh | [CW91] |
| 30 | 132049 | 39751 | 79502 | 1983 | Slowinski | |
| 31 | 216091 | 65050 | 130100 | 1985 | Slowinski | |
| 32 | 756839 | 227832 | 455663 | 1992 | Slowinski & Gage | [Peterson92] |
| 33 | 859433 | 258716 | 517430 | 1994 | Slowinski & Gage | |
| 34 | 1257787 | 378632 | 757263 | 1996 | Slowinski & Gage | (web page) |
| 35 | 1398269 | 420921 | 841842 | 1996 | Armengaud, Woltman,
et. al. (GIMPS) | (web page) |
| ?? | 2976221 | 895932 | 1791864 | 1997 | Spence, Woltman, et. al. (GIMPS) | (web page) |
| ?? | 3021377 | 909526 | 1819050 | 1998 | Clarkson, Woltman, Kurowski et. al. (GIMPS, PrimeNet) | (web page) |
It is not known if these last prime prime below is the 37th Mersenne prime as the region between them and the previous primes has not been completely tested. (Notice that the 29th Mersenne was found while scanning the region between two primes found five years earlier.) GIMPS proved 2756839-1 is the 32nd (January 15, 1997), 2859433-1 is the 33rd (March 28, 1997), 21257787-1 is the 34th (August 28, 1997) and 21398269-1 is the 35th (October 11, 1997) by completing the tests for the remaining smaller exponents. All exponents less than 1,481,800 have now been tested at least once. See Mersenne Status for more current information (and how you may help).
Lucas-Lehmer Test: For p odd, the Mersenne number 2p-1 is prime if and only if 2p-1 divides S(p-1) where S(n+1) = S(n)2-2, and S(1) = 4. [Proof.](It is also possible to start with S(1)=10 and certain other values depending on p.) In pseudocode this test is:
The theory for this test was initiated by Lucas in the late 1870's and then made into this simple test about 1930 by Lehmer. The sequence S(n) is computed modulo 2p-1 to save time. This test is ideal for binary computers because the division by 2p-1 (in binary) can be done using rotation and addition only. (See the pages on proving primality for more information of proving numbers are prime.)Lucas_Lehmer_Test(p): s := 4; for i from 3 to p do s := s2-2 mod 2p-1; if s == 0 then 2p-1 is prime else 2p-1 is composite;
In 1811 Peter Barlow wrote in his text Theory of Numbers that 230(231-1) "is the greatest [perfect number] that will be discovered; for as they are merely curious, without being useful, it is not likely that any person will attempt to find one beyond it." I wonder what he would have made of the first attempts to climb Mount Everest, to run faster miles, or to jump a longer broad jump--other tasks that are curious but not useful. Obviously no one in the late 1800's had any idea of the power of modern computers. What might we know about the machines of 50 years from now? (See also "Why find big primes?")
After the 23rd Mersenne prime was found at the
University of Illinois, the mathematics department was so proud that they
had their postage meter changed to stamp "211213-1 is prime" on
each envelope.
The 25th and 26th Mersenne primes were found by high-school students Laura Nickel and Landon Curt Noll, who, though they had little understanding of the mathematics involved, used Lucas' simple test on the local university's mainframe (CSUH's CDC 174) to find the next two primes. Their discovery of the first prime made the national television news and the front page of the New York times. They went their separate ways after finding the first prime, but Noll kept the program running to find the second--so Noll claims complete ownership. Noll searched later, and though he never found another Mersenne prime, he is one of a team that holds the record for the largest non-Mersenne prime. He currently works for Silicon Graphics.
Slowinski, who works for Cray computers, has written a version of the Lucas test that he has convinced many Cray labs around the world to run in their spare time (time that would be lost otherwise). He had to delay announcing one of his prime records until he got permission to begin looking for it. Slowinski's search for record primes is "not so organized as you would suppose" (his words), as he does not search systematically. In fact, looking at the table of Mersennes you see he missed the 29th prime but found the 30th and 31st. Colquitt & Welsh worked to fill in the gaps and found the 29th.
Enter George Woltman, an excellent programmer and organizer. Starting in late 1995 he gathered up the disparate databases and combined them into one. Then he placed this database, and a free, highly optimized program for search for Mersennes onto the web. This began GIMPS (the Great Internet Mersenne Prime Search) which: has now found the three largest known Mersennes (via Clarkson, Armengaud and Spence), has scanned all of the regions left unexplored between the previous record primes, combines the efforts of dozens of experts and thousands of amateurs, and which offers free software for most computer platforms.
In late 1997 Scott Kurowski (and others) established PrimeNet to automate the selection of ranges and reporting of results for GIMPS, now almost anyone can join this search!
If k>1 and p=4k+3 is prime, then 2p+1 is prime if and only if 2p = 1 (mod 2p+1).So if p=4k+3 and 2p+1 are prime then the Mersenne number 2p-1 is composite (and it seems reasonable to conjecture that there are infinitely many primes pairs such p, 2p+1).
Let p be any odd natural number. If two of the following conditions hold, then so does the third:Notice how this conjecture is related to the theorem in the previous conjecture. This conjecture has been verified for all primes p<100,000.
- p = 2k+/-1 or p = 4k+/-3
- 2p-1 is a prime (obviously a Mersenne prime)
- (2p+1)/3 is a prime.
It seems very unlikely that C5 (or many of the larger terms) would be prime, so this is no doubt another example of Guy's strong law of small numbers [Guy94]. (Landon Curt Noll tells me he has used his program calc to verify that C5 has no prime divisors below 5.1050.)
C0 = 2 (prime) C1 = 3 (prime) C2 = 7 (prime) C3 = 127 (prime) C4 = 170141183460469231731687303715884105727 (prime) C5 > 1051217599719369681879879723386331576246 (is C5 prime ?)