arxiv: v5 [math.nt] 6 Dec 2012

Transcription

1 Lerch Quotients, Lerch Primes, Fermat-Wilson Quotients, and the Wieferich-non-Wilson Primes 2, 3, arxiv: v5 [math.nt] 6 Dec 2012 Jonathan Sondow 209 West 97th Street, New York, NY jsondow@alumni.rinceton.edu Abstract The Fermat quotient q (a) :=(a 1 1)/, for rime a, and the Wilson quotient w :=(( 1)!+1)/ are integers. If w, then is a Wilson rime. For odd, Lerch roved that ( 1 a=1 q (a) w )/ is also an integer; we call it the Lerch quotient l. If l we say is a Lerch rime. A simle Bernoulli-number test for Lerch rimes is roven. There are four Lerch rimes 3,103,839,2237 u to ; we relate them to the known Wilson rimes 5,13,563. Generalizations are suggested. Next, if is a non-wilson rime, then q (w ) is an integer that we call the Fermat-Wilson quotient of. The GCD of all q (w ) is shown to be 24. If q (a), then is a Wieferich rime base a; we give a survey of them. Taking a=w, if q (w ) we say is a Wieferich-non-Wilson rime. There are three u to 10 7, namely, 2,3, Several oen roblems are discussed. 1 Introduction By Fermat s little theorem and Wilson s theorem, if is a rime and a is an integer not divisible by, then the Fermat quotient of base a, and the Wilson quotient of, q (a) := a 1 1, (1) w := ( 1)!+1, (2) are integers. (See [26,. 16 and 19] and [27, ].) For examle, the Fermat quotients of the rime = 5 base a = 1,2,3,4 are q 5 (a) = 0,3,16,51; the Fermat quotients of =3,5,7,11,13,17,19,23,29,31,... base a=2 are 1

2 2 Jonathan Sondow q (2)= = 1,3,9,93,315,3855,13797,182361, , ,... [29, sequence A007663]; and the Wilson quotients of =2,3,5,7,11,13,17,... are w = 1,1,5,103,329891, , ,... [29, sequence A007619]. A rime is called a Wilson rime [16, section A2], [26,. 277] if divides w, that is, if the suercongruence ( 1)!+1 0 (mod 2 ) holds. (A suercongruence is a congruence whose modulus is a rime ower.) For = 2,3,5,7,11,13, we find that w 1,1,0,5,1,0 (mod ) (see [29, sequence A002068]), and so the first two Wilson rimes are 5 and 13. The third and largest known one is 563, uncovered by Goldberg [14] in Crandall, Dilcher, and Pomerance [5] reorted in 1997 that there are no new Wilson rimes u to The bound was raised to by Costa, Gerbicz, and Harvey [4] in Vandiver in 1955 famously said (as quoted by MacHale [22,. 140]): It is not known if there are infinitely many Wilson rimes. This question seems to be of such a character that if I should come to life any time after my death and some mathematician were to tell me that it had definitely been settled, I think I would immediately dro dead again. As analogs of Fermat quotients, Wilson quotients, and Wilson rimes, we introduce Lerch quotients and Lerch rimes in Section 2, and Fermat-Wilson quotients and Wieferich-non-Wilson rimes in Section 3. We define them by combining Fermat and Wilson quotients in aarently new ways. 2 Lerch quotients and Lerch rimes In 1905 Lerch [21] roved a congruence relating the Fermat and Wilson quotients of an odd rime. Lerch s Formula. If a rime is odd, then 1 q (a) w a=1 (mod ), that is,

3 Lerch Quotients and Primes, Fermat-Wilson quotients, and the WW Primes 2, 3, a 1 ( 1)! 0 (mod 2 ). (3) a=1 Proof. Relace a with ab in equation (1). Substituting a 1 = q (a)+1 and b 1 = q (b)+1, we deduce Eisenstein s logarithmic relation [11] q (ab) q (a)+q (b) (mod ) and Lerch s formula follows. For details, see [21] or [30]. Ribenboim [27,. 218] exlains the oint of Lerch s formula this way: Since the Fermat quotient is somehow hard to comute, it is more natural to relate their sums, over all the residue classes, to quantities defined by. Wilson quotients and Lerch s formula have been used (see [30]) to characterize solutions of the congruence 1 n + 2 n + + k n (k+1) n (mod k 2 ). 2.1 Lerch quotients Lerch s formula allows us to introduce the Lerch quotient of an odd rime, by analogy with the classical Fermat and Wilson quotients of any rime. Definition 1. The Lerch quotient of an odd rime is the integer For instance, l := 1 a=1 q (a) w = 1 a=1 a 1 ( 1)! 2. l 5 = = The Lerch quotients of =3,5,7,11,13,17,19,23,29,... are = 13. l = 0,13,1356, , , , , , ,..., and for rime the only Lerch quotientl that is itself a rime number is l 5 = 13 (see [29, Sequence A197630]). By contrast, the Wilson quotients w of the rimes = 5, 7, 11, 29, 773, 1321, 2621 are themselves rime [16, Section A2], [29, Sequence A050299].

4 4 Jonathan Sondow 2.2 Lerch Primes and Bernoulli Numbers We define Lerch rimes by analogy with Wilson rimes. Definition 2. An odd rime is a Lerch rime if dividesl, that is, if 1 a 1 ( 1)! 0 (mod 3 ). (4) a=1 For =3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83, 89,97,101,103,..., we find that l 0,3,5,5,6,12,13,3,7,19,2,21,34,33,52,31,51,38,32,25,25,25, 53,22,98,0,... (mod ) [29, Sequence A197631], and so the first two Lerch rimes are 3 and 103. We give a test for Lerch rimes involving Bernoulli numbers. Ubiquitous in number theory, analysis, and toology (see Dilcher [8]), they are rational numbers B n defined imlicitly for n 1 by the symbolic recurrence relation (B+1) n+1 B n+1 = 0. (Ribenboim [27,. 218] says, Treat B as an indeterminate and, after comuting the olynomial in the left-hand side, relace B k by B k. ) Thus for n=1, we have (B+1) 2 B 2 = 2B = 0, and so B 1 = 1/2. Now with n = 2, we see that (B+1) 3 B 3 = 3B 2 + 3B 1 + 1=0 leads to B 2 = 1/6. In this way, we get B 3 = 0,B 4 = 1 30,B 5 = 0,B 6 = 1 42,B 7 = 0,B 8 = 1 30,B 9 = 0,B 10 = 5 66,... In 1937 (before the era of high-seed comuters!) Emma Lehmer [20] showed that 5 and 13 are the only Wilson rimes 211. To do this, she used her husband D. H. Lehmer s table of Bernoulli numbers u to B 220, together with Glaisher s congruence [12] (see also [21]), which holds for any rime : Here recall the definition a b 0 w B (mod ). (5) (mod m) m a, GCD(a,b)=1. For examle, that 5 is a Wilson rime, but 7 is not, follows from the congruences

5 Lerch Quotients and Primes, Fermat-Wilson quotients, and the WW Primes 2, 3, w 5 B = 5 0 (mod 5), 6 w 7 B = 5 0 (mod 7). 6 Multilying Glaisher s congruence by and substituting w = ( 1)!+1 yields E. Lehmer s test: A rime is a Wilson rime if and only if B 1 1 (mod 2 ). We rovide an analogous test for Lerch rimes. Theorem 1 (Test For Lerch Primes). A rime >3 is a Lerch rime if and only if B 1 +( 1)! (mod 3 ). (6) Proof. We first establish the following Criterion: an odd rime is a Lerch rime if and only if (B+ ) 2 +! (mod 4 ). (7) To see this, recall the classical alication of Bernoulli numbers called Faulhaber s formula (also known as Bernoulli s formula Knuth [19] has insights on this): 1 n + 2 n + +(k 1) n = (B+k)n+1 B n+1. (8) n+1 (See Conway and Guy [2, ] for a lucid roof.) Now set k = and n = 1 in (8). It turns out that B = 0 (indeed, B 3 = B 5 = B 7 = B 9 = = 0; see [2,. 109], [17, section 7.9]), and it follows that the congruences (4) and (7) are equivalent. This roves the Criterion. To rove the Test, note that for any odd ositive integer, the vanishing of B 2k+1 for k 1 imlies (B+ ) = + 1 B 1 + ( 1)/2 k=1 ( ) 2k B 2k. (9) 2k The von Staudt-Clausen theorem [2,. 109], [17, section 7.9], [26,. 340] says in art that the denominator of B 2k is the roduct of all rimes q for which(q 1) 2k. (For instance, as (2 1) 2 and (3 1) 2, the denominator of B 2 is 2 3, agreeing with B 2 = 1/6.) Thus, if is an odd rime, then on the right-hand side of (9) only B 1 has denominator divisible by. From this we see, for 5, that 4 divides the numerator of each term excet 2 B 1. (For the k = ( 3)/2 term, this uses ( 3).) Therefore, the congruence (B+ ) 2 B 1 (mod 4 ) (10)

6 6 Jonathan Sondow holds for all rimes > 3. Substituting (10) into Criterion (7) and dividing by, we arrive at Test (6). As a bonus, (10) affords a roof of Glaisher s congruence. Corollary 1. The congruence (5) holds. Equivalently, if is any rime, then B 1 +( 1)! (mod 2 ). (11) Proof. To see the equivalence, substitute (2) into (5) and multily by. To rove (11), first verify it for =2 and 3. If >3, use (3), (8), and the fact that B = 0 to get (B+ ) 2 +! (mod 3 ). Then (10) and division by yield (11). Notice that the congruences (6) and (11) are the same, excet that in (6) the modulus is 3, while in (11) it is 2. However, one cannot rove Corollary 1 trivially (by reducing (6) modulo 2 instead of 3 ), because (6) holds only for Lerch rimes, whereas (11) holds for all rimes. 2.3 Comuting Lerch rimes: a surrising crossover Let us comare two methods of comuting Lerch rimes: Definition (4) and Test (6). Both require, essentially, comutation modulo 3. The Test seems simler, but on the other hand it requires comuting B 1 modulo 2. To find out which is faster, we used the code If[Mod[Sum[PowerMod[a,-1,ˆ3], {a,1,-1}] - - (-1)!, ˆ3] == 0, Print[]] in a Mathematica (version 7.0.0) rogram for (4), and we used the code If[Mod[Numerator[*Mod[BernoulliB[-1],ˆ2] - - (-1)!], ˆ3] == 0, Print[]] in a rogram for (6). Here Mod[a,m] gives a mod m, PowerMod[a,b,m] gives a b mod m (and is faster than Mod[aˆb,m]), and BernoulliB[k] gives B k. Table 1 shows the CPU time (on a MacBook Air comuter with OS X 10.6 and 2.13GHz Intel rocessor) for each rogram to decide whether is a Lerch rime. Note the surrising crossover in the interval : before it, Test (6) is much faster than Definition (4), but after the interval the reverse is true. Notice also that for >10 4 the CPU times of (4) grow at about the same rate as, while those of (6) balloon at more than double that rate. The rograms for (4) and (6) searched u to 10 4 in about 47.3 and 0.6 seconds, resectively, and found the Lerch rimes 3, 103, 839, and 2237 (see [29, Sequence A197632]). There are no others u to 10 6, by the rogram for (4), which consumed about 160 hours. (To run the rogram for (6) that far u was not feasible.) Marek Wolf, using a modification of (4), has comuted that there are no Lerch rimes in the intervals and , as

7 Lerch Quotients and Primes, Fermat-Wilson quotients, and the WW Primes 2, 3, well as His comutation took six months of CPU time on a 64-bit AMD Oteron 2700 MHz rocessor at the cluster [18]. CPU time in seconds Definition vs. Test > > > > > < < < < < < Table 1: Time each of two rograms takes to comute whether is a Lerch rime. 2.4 Generalizations Euler and Gauss extended Fermat s little theorem and Wilson s theorem, resectively, to congruences with a comosite modulus n see [17, Theorems 71 and 129]. The corresonding generalizations of Fermat and Wilson quotients and Wilson rimes are called Euler quotients q n (a), generalized Wilson quotients w n, and Wilson numbers n w n (see [29, sequences A and A157250]). (The w n are not called Gauss quotients; that term aears in the theory of hyergeometric functions.) In 1998 Agoh, Dilcher, and Skula [1, Proosition 2.1] (see also Dobson [9] and Cosgrave and Dilcher [3]) extended Lerch s formula to a congruence between the q n (a) and w n. Armed with these facts, one can define generalized Lerch quotientsl n and Lerch numbers n l n. But that s another story for another time. 2.5 Oen Problems To conclude this section, we ose some oen roblems. 1. Isl 5 = 13 the only rime Lerch quotient? 2. Is there a fifth Lerch rime? Are there infinitely many? Of the rimes <10 6, only four are Lerch rimes. Thus the answer to the next question is clearly yes; the only thing lacking is a roof!

8 8 Jonathan Sondow 3. Do infinitely many non-lerch rimes exist? As the known Lerch rimes 3,103,839,2237 are distinct from the known Wilson rimes 5,13,563, we may ask: 4. Is it ossible for a number to be a Lerch rime and a Wilson rime simultaneously? Denoting the nth rime by n, the known Wilson rimes are 3, 6, 103. The rimes among the indices 3,6,103, namely, 3 and 103, are Lerch rimes. This leads to the question: 5. If n is a Wilson rime and n is rime, must n be a Lerch rime? The answer to the converse question if n is a Lerch rime, must n be a Wilson rime? is no: 839 and 2237 lie strictly between 563 and , where according to [5] there are no Wilson rimes. In connection with Problem 5, comare Davis s Are there coincidences in mathematics? [6] and Guy s The strong law of small numbers [15]. 3 Fermat-Wilson quotients and the WW rimes 2, 3, Suose that a rime is not a Wilson rime, so that does not divide its Wilson quotient w. Then in the Fermat quotient q (a) of base a, we may take a=w. Definition 3. If is a non-wilson rime, then the Fermat-Wilson quotient of is the integer q (w )= w 1 1. For short we write g := q (w ). The first five non-wilson rimes are 2,3,7,11,17. Since w 2 = w 3 = 1, w 7 = 103, and w 11 = , the first four Fermat-Wilson quotients are g 2 = g 3 = 0, and g 7 = = , g 11 = = [29, Sequence A197633]. The fifth one, g 17, is a 193-digit number.

9 Lerch Quotients and Primes, Fermat-Wilson quotients, and the WW Primes 2, 3, The GCD of all Fermat-Wilson quotients We saw that at least one Lerch quotient and seven Wilson quotients are rime numbers. What about Fermat-Wilson quotients? Theorem 2. The greatest common divisor of all Fermat-Wilson quotients is 24. In articular, q (w ) is never rime. Proof. The rime factorizations of q (w )=g for =7 and 11 are and g 7 = g 11 = Since g 2 = g 3 = 0, we thus have GCD(g 2,g 3,g 7,g 11 )=2 3 3=24. To comlete the roof, we show that 24 divides g whenever >3. Since w =( 1)!+1, it is clear that if 5, then w, and hence w, is not divisible by 2 or 3. As even owers of such numbers are 1 (mod 8) and 1 (mod 3), and so 1 (mod 24), it follows that g (= w 1 1), and hence g, is divisible by Wieferich rimes base a Given an integer a, a rime is called a Wieferich rime base a if the suercongruence a 1 1 (mod 2 ) (12) holds. For instance, 11 is a Wieferich rime base 3, because =59048= Parahrasing Ribenboim [26,. 264], it should be noted that, contrary to the congruence a 1 1 (mod ) which is satisfied by every rime not dividing a, the Wieferich suercongruence (12) is very rarely satisfied. When it is, cannot divide a, and so the Fermat quotient q (a) is an integer. In fact, (1) shows that a rime is a Wieferich rime base a if and only if does not divide a but does divide q (a).

10 10 Jonathan Sondow In 1909, while still a graduate student at the University of Münster in Germany, Wieferich created a sensation with a result related to Fermat s Last Theorem: If x + y = z, where is an odd rime not dividing any of the integers x,y, or z, then is a Wieferich rime base 2. One year later, Mirimanoff roved that is also a Wieferich rime base 3. (See [7, ], [27, Chater 8], and [31,. 163].) The only known Wieferich rimes base 2 (also simly called Wieferich rimes) are 1093 and 3511, discovered by Meissner in 1913 and Beeger in 1922, resectively. In 2011 Dorais and Klyve [10] comuted that there are no others u to It is unknown whether infinitely many exist. (Neither is it known whether there are infinitely many non-wieferich rimes base 2. However, Silverman has roved it assuming the abc-conjecture see his leasantly-written aer [28].) Likewise, only two Wieferich rimes base 3 (also known as Mirimanoff rimes) have been found, namely, 11 and The second one was uncovered by Kloss in An unanswered question is whether it is ossible for a number to be a Wieferich rime base 2 and base 3 simultaneously. (See [16, section A3] and [26, , ].) For tables of all Wieferich rimes base a with 2< <2 32 and 2 a 99, see Montgomery [23]. 3.3 The Wieferich-non-Wilson rimes 2, 3, Let us consider Wieferich rimes base a where a is the Wilson quotient of. Definition 4. Let be a non-wilson rime, so that its Fermat-Wilson quotient q (w ) is an integer. If divides q (w ) equivalently, if the suercongruence w 1 1 (mod 2 ) (13) holds then is a Wieferich rime base w, by definition (12). In that case, we call a Wieferich-non-Wilson rime, or WW rime for short. For the non-wilson rimes =2,3,7,11,17,19,23,29,31,37,41,43,47,53,59, 61,67,71,73,79,83,..., the Fermat-Wilson quotients q (w ) = g are congruent modulo to g 0,0,6,7,9,7,1,6,18,17,30,11,25,30,24,46,64,16,18,4,29,... (mod ) [29, Sequence A197634]. In articular, 2 and 3 are WW rimes. But they are trivially so, because g 2 and g 3 are equal to zero. Is there a non-trivial WW rime? Perhas surrisingly, the answer is yes but the smallest one is In the next subsection, we give some details on using a comuter to show that is a WW rime. It is non-trivial because g In fact, taking logarithms, one finds that

11 Lerch Quotients and Primes, Fermat-Wilson quotients, and the WW Primes 2, 3, g = ( 14770!+1 ) > , so that the number g has more than 800 million decimal digits. 3.4 Comuter search To search for WW rimes, one can use a comuter to calculate whether or not a given rime satisfies condition (13). Exlicitly, if the number ( ) ( 1)!+1 1 mod 2 (14) is equal to 1, then is a WW rime. Mathematica s function Mod[a,m] can comute (14) when is small. But if is large, an Overflow message results. However, it is easy to see that in (14) one may relace( 1)! with( 1)! mod 3, a much smaller number. For examle, it takes just a few seconds for a rogram using the code If[PowerMod[(Mod[(-1)!, ˆ3] + 1)/, -1, ˆ2] == 1, Print[]] to test the first 2000 rimes and rint the WW rimes 2,3,14771 (see [29, Sequence A197635]). Michael Mossinghoff, emloying the GMP library [13], has comuted that there are no other WW rimes u to More oen roblems We conclude with three more oen roblems. 6. Can one rove that is a WW rime (i.e., that divides g ) without using a comuter? Such a roof would be analogous to those given by Landau and Beeger that 1093 and 3511, resectively, are Wieferich rimes base 2. (See Theorem 91 and the notes on Chater VI in [17], and History and search status in [32].) However, roofs for Wieferich rimes are comaratively easy, because (high) owers are easy to calculate in modular arithmetic, whereas factorials are unlikely to be calculable in logarithmic time. 7. Is there a fourth WW rime? Are there infinitely many? Comments similar to those receding Problem 3 also aly to the next question. 8. Do infinitely many non-ww rimes exist?

12 12 Jonathan Sondow Is it ossible to solve Problem 3 or Problem 8 assuming the abc-conjecture? (See the remark in Section 3.2 about Silverman s roof.) Acknowledgments I am grateful to Wadim Zudilin for suggestions on the Test, for a simlification in comuting WW rimes, and for verifying that there are no new ones u to 30000, using PARI/GP [24]. I thank Marek Wolf for comuting Lerch rimes, and Michael Mossinghoff for comuting WW rimes. References 1. Agoh, T., Dilcher, K., Skula, L.: Wilson quotients for comosite moduli. Math. Com. 67, (1998) 2. Conway, J.H., Guy, R.K.: The Book of Numbers. Sringer-Verlag, New York (1996) 3. Cosgrave, J.B., Dilcher, K.: Extensions of the Gauss-Wilson theorem. Integers 8, article A39 (2008) 4. Costa, E., Gerbicz, R., Harvey, D.: A search for Wilson rimes (2012, rerint); available at htt://arxiv.org/abs/ Crandall, R., Dilcher, K., Pomerance, C.: A search for Wieferich and Wilson rimes. Math. Com. 66, (1997) 6. Davis, P.J.: Are there coincidences in mathematics? Am. Math. Mon. 88, (1981) 7. Dickson, L.E.: History of the Theory of Numbers, vol. 1. Carnegie Institution of Washington, Washington, D. C. (1919); rerinted by Dover, Mineola, NY (2005) 8. Dilcher, K.: A bibliograhy of Bernoulli numbers (2011); available at htt:// dilcher/bernoulli.html 9. Dobson, J.B.: On Lerch s formula for the Fermat quotient (2012, rerint); available at htt://arxiv.org/abs/ Dorais, F.G., Klyve, D.W.: A Wieferich rime search u to J. Integer Seq. 14, Article (2011); available at htt:// 11. Eisenstein, F.: Eine neue Gattung zahlentheoretischer Funktionen, welche von zwei Elementen abhaengen und durch gewisse lineare Funktional-Gleichungen definiert werden, Verhandlungen der Koenigl. Preuss. Akademie der Wiss. zu Berlin (1850), 36 42; rerinted in Mathematische Werke, vol. 2, Chelsea, New York (1975) 12. Glaisher, J.W.L.: A congruence theorem relating to Eulerian numbers and other coefficients. Proc. Lond. Math. Soc. 32, (1900) 13. GMP: The GNU Multile Precision Arithmetic Library (2011); available at htt://gmlib.org/ 14. Goldberg, K.: A table of Wilson quotients and the third Wilson rime. J. Lond. Math. Soc. 28, (1953) 15. Guy, R.K.: The strong law of small numbers. Am. Math. Mon. 95, (1988) 16. Guy, R.K.: Unsolved Problems in Number Theory, 3rd ed. Sringer, New York (2004) 17. Hardy, G.H., Wright, E.M.: An Introduction to the Theory of Numbers, 5th ed. Oxford University Press, Oxford (1989) 18. Klaster Instytutu Fizyki Teoretycznej UWr (2012); htt://zero.ift.uni.wroc.l/ 19. Knuth, D.: Johann Faulhaber and sums of owers. Math. Com. 61, (1993) 20. Lehmer, E.: A note on Wilson s quotient. Am. Math. Mon. 44, (1937)

13 Lerch Quotients and Primes, Fermat-Wilson quotients, and the WW Primes 2, 3, Lerch, M.: Zur Theorie des Fermatschen Quotienten a 1 1 = q(a). Math. Ann. 60, (1905) 22. MacHale, D.: Comic Sections: The Book of Mathematical Jokes, Humour, Wit and Wisdom. Boole Press, Dublin (1993) 23. Montgomery, P.L.: New solutions of a 1 1 (mod 2 ). Math. Com. 61, (1993) 24. PARI/GP (2011); available at htt://ari.math.u-bordeaux.fr/ 25. Ribenboim, P.: Math. Intelligencer 5, (1983) 26. Ribenboim, P.: The Book of Prime Number Records, 2nd. ed. Sringer-Verlag, New York (1989) 27. Ribenboim, P.: My Numbers, My Friends: Poular Lectures on Number Theory. Sringer- Verlag, New York (2000) 28. Silverman, J.H.: Wieferich s criterion and the abc-conjecture. J. Number Theory 30 (2), (1988) 29. Sloane, N.J.A.: The On-Line Encycloedia of Integer Sequences (2011); ublished at htt://oeis.org/ 30. Sondow, J., MacMillan, K.: Reducing the Erdős-Moser equation 1 n + 2 n + + k n =(k+1) n modulo k and k 2. Integers 11, article A34 (2011); exanded version available at htt://arxiv.org/abs/ v1 31. Wells, D.: The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, London (1986) 32. Wikiedia: Wieferich rime (2012); htt://en.wikiedia.org/wiki/wieferich rime