arxiv:math/ v1 [math.nt] 3 Dec 2004

Transcription

1 arxiv:math/ v1 [math.nt] 3 Dec 2004 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD Abstract. We consider rimitive divisors of terms of integer sequences defined by quadratic olynomials. Aart from some small counterexamles, when a term has a rimitive divisor, that rimitive divisor is unique. It seems likely that the number of terms with a rimitive divisor has a natural density. We discuss two heuristic arguments to suggest a value for that density, one using recent advances made about the distribution of roots of olynomial congruences. 1. Mersenne Numbers and Primitive Prime Divisors A notorious roblem from elementary number theory is the Mersenne Prime Conjecture. This says that the Mersenne sequence M = (M n ) n 1, defined by M n = 2 n 1, will contain infinitely many rime terms known as Mersenne rimes. The sequence roduces some rimes initially, for examle M 2 = 3, M 3 = 7, M 5 = 31, M 7 = 127,... However, the aearance of Mersenne rimes quickly thins out; only 41 are known, the largest of which, M , has over 7 million decimal digits. This was discovered by Josh Findley on May 15, 2004 as art of the GIMPS roject [21], which harnesses idle time on thousands of comuters all over the world to run a distributed version of the Lucas- Lehmer test for rimality of M n. Only 41 rimes might seem rather a small return for such an effort. Anybody looking for gold or gems who made such a small return for so much effort would surely abandon the search. It seems fair to ask why should we believe there are infinitely many. In the absence of a rigorous roof, our exectations may be informed by heuristic arguments. A 1991 Mathematics Subject Classification. 11A41, 11B37. Key words and hrases. rime, rimitive divisor, quadratic olynomial, recurrence sequence, recurrence relation. 1

2 2 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD later section will discuss heuristic arguments for this roblem and other more or less tractable roblems in Number Theory. In 1892, Zsigmondy [22] discovered a beautiful argument which shows that the sequence M does yield infinitely many rime numbers but in a less restrictive sense. Given any sequence S = (S n ) n 1, define a rimitive divisor of the term S n 0 to be a rime number which divides S n but which does not divide any nonzero term S m with m < n. Factorizing the first few terms of the Mersenne sequence reveals several rimitive divisors, which are underlined below: M 1 = 1, M 2 = 3, M 3 = 7, M 4 = 3 5, M 5 = 31, M 6 = 3 2 7, M 7 = 127,.... Notice that the term M 6 has no rimitive divisor. Zsigmondy [22] roved that all the terms M n for n > 6 have a rimitive divisor. More generally, Zsigmondy roved a similar result for sequences of the form U n = a n b n, whenever a > b are ositive corime integers: he showed that a term U n will always have a rimitive divisor unless a = 2, b = 1 and n = 6 or a + b is a ower of 2 and n = 2. Any sequence of this form is of articular interest because it satisfies a linear recurrence relation. A sequence U = (U n ) is called a binary linear recurrence sequence if it satisfies the linear recurrence relation U n+2 = cu n+1 + du n for all n 1 for constants c and d. Zsigmondy s sequences U n = a n b n satisfy this relation with c = a + b and d = ab. Although Zsigmondy s result is much weaker than the Mersenne Prime Conjecture, it initiated a great deal of interest in the arithmetic of general linear recurrence sequences (see [10, Cha. 6]). Schinzel [15] extended Zsigmondy s result, giving further insight into the finer arithmetic of sequences like M. He established a large class of indices for which M n has a comosite rimitive divisor (that is, for which M n has more than one rimitive divisor). As an examle of his results, he roved that M 4k has a comosite rimitive divisor for all odd k > 5. Other researchers have removed the constraint that a and b be integers in Zsigmondy s theorem. Let u and v denote conjugate quadratic integers; in other words, zeros of a monic irreducible olynomial with integer coefficients. The integer sequences C and D defined by C n (u, v) = (u n v n )/(u v) and D n (u, v) = u n + v n are both binary linear recurrence sequences. Bilu, Hanrot and Voutier [3] used owerful methods from Diohantine analysis to rove that both C n and D n have a rimitive divisor for any n > 30. The two striking asects of this result are the uniform nature of the bound and its small

3 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES 3 numerical value. In articular, for any given sequence it is easy to check the first 30 terms for rimitive divisors, arriving at a comlete icture. For instance, the Fibonacci sequence A 1 = 1, A 2 = 1, A 3 = 2,... has A n = C n ( 1 2 (1 + 5), 1 2 (1 5) ), so a calculation reveals that A n does not have a rimitive divisor if and only if n = 1, 2, 6 or 12. The theory of linear recurrence sequences has a bilinear analogue. For examle, the Somos-4 sequence V = (V n ) satisfies the bilinear recurrence relation This sequence begins V n+4 V n = V n+3 V n+1 + Vn+2 2, for all n 1. 1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529,8209, , , ,... Remarkably, all the terms are integers even though calculating V n+4 a riori involves dividing by V n. This sequence was discovered by Michael Somos [17] and it is known to be associated to the arithmetic of ellitic curves. See [10, Sect and 11.1] for a summary of this, and further references, including a remarkable observation due to Pro et al. that the terms of the sequence must be integers because they count matchings in a sequence of grahs. Silverman [16] showed a qualitative analogue of Zsigmondy s result for ellitic curves which alies, in articular, to this sequence. An exlicit form of this result roved by Everest, McLaren and Ward [8] guarantees that from V 5 onwards, all terms have a rimitive divisor. 2. Polynomial Sequences Consider now the sequence P defined by P n = n 2 + β. Any integer sequence defined by a monic quadratic olynomial also satisfies a linear recurrence relation; P satisfies P n+3 = 3P n+2 3P n+1 + P n for all n 1. It has long been susected that for any fixed β with β not a square, the sequence P will contain infinitely many rime terms. No roof is known even for one value of β. It seems reasonable to ask the aarently simler question about the existence of rimitive divisors of terms. Clearly any rime term will itself be a rimitive divisor but do the comosite terms have rimitive divisors? Comutations suggest the following.

4 4 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD Conjecture 2.1. Suose f denotes a quadratic olynomial with integer coefficients, and let π f (N) denote the number of terms f(n) with n < N having a rimitive divisor in the sequence (f(n)). If the zeros of f are irrational then π f (N) c f N for some constant c f > 0 deending on f. In the next section we find a necessary and sufficient condition for a term to have a rimitive divisor. It is not clear why the terms of the sequence should fail, infinitely often, to roduce rimitive divisors. This will be roved in a subsequent section. Finally, an attemt will be made to quantify this failure, using heuristic arguments. To simlify the analysis, we will restrict our attention to quadratic olynomials of the form f(x) = x 2 + β. From now on P will denote the sequence defined by P n = n 2 + β. 3. Losided numbers Since we are mainly interested in asymtotic behaviour, assume from now on that n > β. The terms with n β are not guaranteed to fit the behaviour described below. Consider first a rime < n dividing P n. Then P n 0 (mod ), so P m 0 (mod ) for some m < simly by choosing m to be the residue of n (mod ). By assumtion, n so m < n. In others words, is not a rimitive divisor of P n. This means that to find rimitive divisors of P n we have to look for rime divisors which are greater than n. We can say more: We can guarantee a solution of P m 0 (mod ) for some m /2. Thus, to find rimitive divisors we have only to look amongst the rime divisors which are bigger than 2n. This turns out to characterize those values of P which do have a rimitive divisor. Definition 3.1. An integer k is losided if it has a rime factor q with q > 2 k. Thus any rime greater than 3 is losided; 22, 26, 34 are losided, while 24 and 28 are not. Proosition 3.2. For all n > β, the term P n has a rimitive divisor if and only if P n is losided. For all n > β, if P n has a rimitive divisor, then that rimitive divisor is unique. The requirement that n > β is necessary: if β is rime then P β has rimitive divisor β but is not losided. Also, terms with small n may have more than one rimitive divisor: for examle, the sequence of values of the olynomial n begins 7, 10,... so the second term

5 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES 5 has two rimitive divisors. The kind of results discussed here are asymtotic results, so this restriction is unimortant. Proof of Proosition 3.2. As shown above, a rime dividing P n is a rimitive divisor only if > 2n. Conversely, suose that is a rime dividing P n which is not a rimitive divisor. Then n 2 + β 0 (mod ) and there is an integer m < n with m 2 + β 0 (mod ) so (by subtracting the two congruences) m 2 n 2 0 (mod ). It follows that m ± n 0 (mod ), and in articular m + n < 2n. Thus a rime is a rimitive divisor of P n if and only if divides P n and > 2n. To comlete the roof of the first statement we claim that, for n greater than β, P n has such a rime divisor if and only if P n is losided: If is a rime dividing P n and > 2n then 2n + 1 > 2 n 2 + n > 2 n 2 + β. Conversely, if > 2 n 2 + β then 2 n 2 n + 1 > 2n 1 so > 2n (since 2n cannot be rime). The uniqueness of the rimitive divisor follows at once. If is a rime dividing P n and > 2 P n, then no other rime divisor can be as large, so cannot be rimitive. Proosition 3.2 gives some insight into the arithmetic of the seqence P = (P n ) but does not say anything about how often rimitive divisors should occur. 4. Terms With and Terms Without Primitive Divisors Recalling that P n = n 2 + β, write N Q N = P n. n=1 Theorem 4.1. Denote by ω(q N ) the number of rime divisors of Q N. Given any constant C > 0, ω(q N ) < N CN for all sufficiently large N. There is a constant c > 0 such that N 2 cn < ω(q N) for all sufficiently large N.

6 6 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD Corollary 4.2. There are infinitely many terms of the sequence P which do not have a rimitive divisor. Asymtotically, at least one half of all terms of the sequence have a rimitive divisor. The roof of the corollary follows by the earlier remark that with finitely many excetions, a rimitive divisor is unique. It will become clear during the roof of Theorem 4.1 that a slightly stronger result is rovable with these methods, namely ω(q N ) < N N log for all sufficiently large N. (1) This refinement does not hel us directly, so we will not rove it here. Proof of Theorem 4.1. The roof uses little aart from well-known estimates for sums over rimes, which can be found in Aostol [1]. Aart from a finite number of rimes, any rime that divides n 2 +β has the roerty that β is a quadratic residue modulo. So let R denote the set of odd rimes for which β is a quadratic residue; notice that R comrises the intersection of a finite union of arithmetic rogressions with the set of rimes. Given integers K > 1 and N > K, fix three subsets of R as follows: S T U = { R < N}; = { R Q N, N < < KN}; = { R Q N, KN < }. Note that all S divide Q N. Write s = S, t = T and u = U. We wish to bound above and below the exression s + t + u. Let A(x) = { R < x}. Dirichlet s Theorem on rimes in arithmetic rogression imlies that A(x) x 2 log x, so asymtotically. s = N ( ) N 2 + o (2)

7 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES 7 Notice that log Q N = so by Stirling s Formula N log n 2 + β n=1 On the other hand, we may write N ( = 2 log n + O ( ) ) 1 n 2 n=1 ( ) N = 2 log n + O(1), n=1 log Q N = 2N 2N + O(1). (3) log Q N = Q N e log, (4) corresonding to the rime decomosition Q N = Q N e for ositive integers e. The sum on the right of (4) decomoses according to the definitions of S, T and U to give log Q N = S e log + T log + U log + O(N), (5) noting that e = 1 whenever > N. The term O(N) takes account of the finite number of rimes we excluded. Now for every S, e. Hence there is a constant c > 0 such that the right-hand side of (5) is bounded below by 2N log + t + u log(kn) cn. S By Aostol [1, Th. 7.3], S 2N 2N Hence there is a constant a > 0 such that S 2N log = N + O(N). (6) log > N an. Putting this together and using (3) gives 2N 2N + O(1) > N (a + c)n + t + u log(kn),

8 8 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD and hence an uer bound for u, u < N log(kn) t (a + c 2)N + log(kn) log(kn) which is bounded above by N + (a + c 2 log K)N + o(1), ( t 1 log K ) ( ) N + o. Using this and invoking (2) gives ( ) log K (a + c 3/2 log K)N s + t + u < N + t + Now so t A(KN) A(N) < s + t + u < N + ( ) KN N 2 log(kn) + o, (a + c 3/2 log K)N ( ) N + o. ( ) N + o. Given C, choose K so that a + c 3/2 log K < C, then for all sufficiently large N ω(q N ) = s + t + u + O(1) < N CN. Next we obtain a lower bound for s+t+u. For each S and k N, denote by ord (k) the index of the greatest ower of dividing k and ut B (N) = {n < N ord (n 2 + β) > 1}. Then e log = 2N log + ord (n 2 + β) log. S S S n B (N) For each S, β has two -adic square roots, and ord (n 2 +β) = r+1 if and only if the -adic exansion of n agrees with one of these square roots u to the term in r and no further. Hence ( ) ord (n 2 N + β) r 2 log S n B (N) log = S < 2N S r 1 log ( 1) 2, r+1

9 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES 9 which is O(N) since the last sum converges. Putting this together with (6) shows that there is a b > 0 such that e log < N + bn. S Together with (3) and (5), this gives a d > 0 such that N (2 + b + d)n < (t + u) log(n 2 + β), using N 2 + β, and therefore N 2 (2 + b + d)n 2 Invoking (2) again gives the lower bound ( N (1 + b + d)n N + o 2 2 ( ) N + o < t + u. ) < s + t + u. Finally, choosing any c with c > (1 + b + d)/2 gives the result. As an exercise, we suggest roving the stronger statement (1) using a similar argument. 5. Heuristic Arguments There are a number of ways that mathematicians form a view on which statements are likely to be true. These views inform research directions and hel to concentrate effort on the most fruitful areas of enquiry. The only certainty in mathematics comes from rigorous roofs that adhere to the rules of logic; the discourse of logos. When such a roof is not available, other kinds of arguments can make mathematicians exect that statements will be true, even though these arguments fall well short of a roof. These are called heuristic arguments the word comes from the Greek root ÙÖ (Eureka) meaning, I have found it. It usually means the rinciles used to make decisions in the absence of comlete information or the ability to examine all ossibilities. In informal ways, mathematicians use heuristic arguments all the time when they discuss mathematics, and these are art of the mythos discourse in mathematics. One consequence of the Prime Number Theorem is the following statement: the robability a large integer N is rime is roughly 1/. What this means is that if an element of the set {1,...,N} is chosen at random using a fair N-sided die for each N, then the robability ρ N that it is rime satisfies ρ N 1 as N. This crude estimate

10 10 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD has been used several times to argue heuristically in favour of the lausibility of difficult roblems. Some examles follow in each case the argument resented falls well short of a roof, yet it still seems to have some redictive ower and has suggested lines of attack Fermat rimes. Hardy and Wright [11, Sect. 2.5] argued along these lines that there ought to be only finitely many Fermat rimes. A Fermat rime is a rime number of the form F n = 2 2n + 1. Fermat demonstrated that F 1,...,F 4 are all rimes; Euler showed that F 5 is comosite using congruence arguments. Since then many Fermat numbers have been shown to be comosite and quite a few have been comletely factorized. The number of Fermat rimes F n with n < N, if they are no more or less likely to be rime than a random number of comarable size, should be roughly 1 1 log F n 2 n log 2 < 1 for all N 1. log 2 n<n n<n Statements like this cannot be taken too literally the numbers F n have many secial roerties, not all of which are understood. However, this kind of argument tends to suort the belief that there are only finitely many Fermat rimes, and would incline many mathematicians to attemt to rove that statement rather than its negation. Massive advances in comuting ower suggest that we know indeed that Fermat knew all the Fermat rimes Mersenne rimes. The Prime Number Theorem can also be used to argue in suort of the Mersenne Prime Conjecture as follows. The exected number of Mersenne rimes M n with n < N is roughly 1 log(2 n 1) n<n if M n is no more or less likely to be a rime than a random number of similar size. This sum is roughly 1 1 log 2 n. n<n Since the Harmonic Series diverges, this suggests that there are infinitely many Mersenne rimes. More can be said: This is asymtotically / log 2 so we could ask if this suggests there are roughly that number of Mersenne rimes M n with n < N. This seems to be a dream too far. Nonetheless, Wagstaff [20], then Pomerance and Lenstra [14]

11 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES 11 have sought to refine this heuristic argument to obtain an asymtotic estimate that does more closely match the available evidence. On the basis of these refined heuristics, they conjecture that the number of Mersenne rimes M n with n < N is asymtotically e γ, log 2 where γ is the Euler-Mascheroni constant. Caldwell s Prime Page [4] is an excellent online resource which gives more details about these arguments and the hunt for new Mersenne rimes Bilinear Recurrence Sequences. Consider now the Somos-4 sequence defined earlier. Amongst the early terms, several are rime. It seems natural to ask if there are infinitely many rime terms. Somos- 4 can be generalized to integer sequences S satisfying a relation S n+4 S n = es n+3 S n+1 + fs 2 n+2, where e and f are integral constants not both zero. Such sequences are often called Somos Sequences and Swart showed [18] how they are related to the arithmetic of ellitic curves. Note that some care is needed because, for examle, a binary linear recurrence sequence will always satisfy a bilinear recurrence relation of this kind: we will refer to a Somos sequence as non-degenerate if it does not satisfy any linear recurrence relation. These are natural generalizations of linear recurrence sequences, so erhas we should exect them to contain infinitely many rime terms. Comutational evidence in [5] tended to suort that belief because of the relatively large rimes discovered. The kind of argument above using the Prime Number Theorem was adated in [7] and it suggests that a non-degenerate Somos sequence should contain only finitely many rime terms. To see why this is so, note that general results about heights on ellitic curves show that the growth rate of S n is quadratic-exonential. In other words n<n log S n cn 2, where c is a ositive constant. Thus the exected number of rime terms with n < N should be aroximately log S n c n π2 for all N. 2 6c n<n This resembles the argument of Hardy and Wright, and suggests that only finitely many rime terms should be exected. Proofs of the finiteness in many secial cases have subsequently been found [9] and

12 12 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD the search for these roofs was motivated in art by the heuristic arguments. Interestingly, it is known that the constant c is uniformly bounded below. Thus the style of heuristic argument suggests that erhas the total number of rime terms is uniformly bounded above across all Somos sequences. Much calculation has failed to yield a Somos sequence with more than a dozen rime terms Quadratic olynomials. Suose β is an integer which is not the negative of a square. Again, the Prime Number Theorem redicts roughly 1 log(n 2 + β) n<n rime terms in the sequence (n 2 + β) with n < N if n 2 + β is neither more nor less likely to be rime than a random number of that size. The sum is asymtotically N/2 which suorts the belief that there are infinitely many rime terms in the sequence (n 2 +β). Comutation suggests that for fixed β there will be cn/ rime terms with n < N, where c = c(β) is a constant which deends uon β. To see that c should deend uon β is easy. For examle, one half of the terms in the sequence will be even, so the robability that such a term is rime is reduced. Bateman and Horn [2] gave a heuristic argument and numerical evidence to suggest that the number of rime values f(n) of an irreducible olynomial f with n < N that takes on infinitely many rime values should be CN/, where C = 1 deg(f) ( 1 1 ) 1 ( 1 w() where the roduct is taken over all rimes and w() denotes the number of solutions x (mod ) to the congruence f(x) 0 (mod ). We aim to suort Conjecture 2.1 by showing firstly that an asymtotic formula can be obtained for counting losided numbers. Theorem 5.1. Let π ı (N) denote the number of losided numbers less than or equal to N. Then π ı (N) N log 2 as N. Theorem 5.1 can be alied to give the following heuristic argument in suort of Conjecture 2.1. The robability that a large integer is losided is roughly log 2. Hence the exected number of losided values of n 2 + β with n < N is asymtotically N log 2. ),

13 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES 13 Proof of Theorem 5.1. Write a losided number as m where is the largest rime factor, so that the losided condition is > 4m. To comute the number of losided numbers below N, note that the counting can be achieved by dividing the set into two arts. Let denote a variable rime. For < 2 N there are /4 losided integers m < N. For 2 N, each m < N is losided so there are N losided integers m. Hence the total number is <2 N N <N N. The first sum is O(N/ ) so it can be ignored asymtotically. The second sum differs from N 1 2 N <N by an amount which is O(N/ logn) by the Prime Number Theorem. To estimate this sum use Mertens Formula [1, Th. 4.12], 1 = log log x + A + o(1). <x Alying this gives [ N log + A log(log ] N + log 2) A + o(1) which is asymtotically N log 2. Comutational evidence suggests that the number of losided terms in n 2 + β is asymtotically cn for some constant c. Comutations with β < 10 suggest the constant c looks reasonably close to log 2 in each case. In the next section we will give another kind of argument, still heuristic, yet closer to known arguments in analytic number theory. Towards the end, this argument resonates with the critical cancellation of the two log terms which occurred in the revious argument. 6. Another Heuristic Argument The style of argument in the roof of Theorem 4.1 is suggestive. In this section, we will adat it, taking great liberties with error terms, to give a more refined heuristic argument in suort of our main conjecture. For every odd R, there is an associated smallest ositive solution n < of the congruence 2 n2 β (mod ). We are going to use a statement about the distribution of the ratios n as varies in R. The study of the distribution of sequences of this tye goes back

14 14 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD a long way. Hooley [12] originally studied and roved a uniform distribution result for a general modulus (that is, where the rime is relaced by a natural number m). In [13] he gave more refined results, again for a general modulus. The first uniform distribution result for a rime modulus was given as recently as 1995 by Duke, Friedlander and Iwaniec [6], but only for β > 0. The general result was roved by Tóth [19] in 2000, and the next result is taken from his aer. Theorem 6.1. For every δ [0, 1/2], let A δ (x) = { R n/ < δ and < x}, where n is the smallest ositive solution n < of 2 n2 β (mod ). Then A δ (x) δx log x. Let B(N) = {n < N n 2 + β is losided} and let k > 1 be a real number. Recall that A(x) = { R < x}. For < 2N, there are A(2N) rimes in R and each has n < < N so each value of n is included. 2 For 2N < < 2kN, there are A 1/2k (2kN) A 1/2k (2N) rimes in R with n < 1. For each of these rimes, the corresonding n 2k has n < < N so it is included. We may nevertheless be missing some 2k rimes here, when 1 < n < 1 ; however, if n > 1 then n > N, so we 2 2k 2k k will only be missing rimes from terms with large n. Now we reeat this for 2kN < < 2k 2 N. We count A 1/2k 2(2k 2 N) A 1/2k 2(2kN) rimes in R for which n < < N. For the missing rimes, note 2k 2 again that if n > 1 then n > N. 2k 2 k Assuming that the rimes not counted do not occur very often (they can only come from terms with n > N ), we obtain the following estimate for k B(N): B(N) A(2N) + A 1/2k (2kN) A 1/2k (2N) +A 1/2k 2(2k 2 N) A 1/2k 2(2kN) +. (7) The right-hand side of (7) simlifies to an asymtotic estimate for ( ) ( ) k 1 1 N k log(2n) + 1 log(2kn) +.

15 PRIMITIVE DIVISORS OF QUADRATIC POLYNOMIAL SEQUENCES 15 The sum of the first R terms of this gives ( ) k 1 R 1 N k log(2n) + r log k. r=0 This sum can be estimated by Euler summation to give ( ) k 1 R ( ) dt k 1 [ N k log(2n) + t log k =N ] R log(log(2n+ t log k) k log k 0 which is asymtotically ( k 1 N k log k 0 ) ( log + log ( 1 + R log k ) ) log. The two log log terms cancel just as before! Suose we were allowed to sum out to R with R ( / ) log k as N. This would give an asymtotic estimate of N log 2. k 1 k log k Now, letting k 1, so that we do not miss out any rimes in our counting, we get N log 2 as exected. However, there is a reason why we may not build such an assumtion into the argument. The estimate for A δ (x) is for a fixed δ. We have no justification for alying it when δ deends on x without a uniformity result which is not available. On the other hand, we may argue heuristically as above in the hoe that the error terms do actually behave well in ractice. References [1] T. M. Aostol, Introduction to analytic number theory, Sringer-Verlag, New York, 1976, Undergraduate Texts in Mathematics. MR 55 #7892 [2] P. T. Bateman and R. A. Horn, A heuristic asymtotic formula concerning the distribution of rime numbers, Math. Com. 16 (1962), MR 26 #6139 [3] Y. Bilu, G. Hanrot, and P. M. Voutier, Existence of rimitive divisors of Lucas and Lehmer numbers, J. Reine Angew. Math. 539 (2001), , With an aendix by M. Mignotte. MR 2002j:11027 [4] C. Caldwell, The Prime Pages, [5] D. V. Chudnovsky and G. V. Chudnovsky, Sequences of numbers generated by addition in formal grous and new rimality and factorization tests, Adv. in Al. Math. 7 (1986), no. 4, MR 88h:11094 [6] W. Duke, J. B. Friedlander, and H. Iwaniec, Equidistribution of roots of a quadratic congruence to rime moduli, Ann. of Math. (2) 141 (1995), no. 2, MR 95k:11124 [7] M. Einsiedler, G. R. Everest, and T. Ward, Primes in ellitic divisibility sequences, LMS J. Comut. Math. 4 (2001), MR 2002e:11181 [8] G. Everest, G. McLaren, and T. Ward, Primitive divisors of ellitic divisibility sequences, arxiv:math.nt/ , 2004.

16 16 G. EVEREST, S. STEVENS, D. TAMSETT AND T. WARD [9] G. Everest, V. Miller, and N. Stehens, Primes generated by ellitic curves, Proc. Amer. Math. Soc. 132 (2004), no. 4, MR 2005a:11076 [10] G. R. Everest, A. J. van der Poorten, I. Sharlinski, and T. Ward, Recurrence sequences, Mathematical Surveys and Monograhs, vol. 104, American Mathematical Society, Providence, RI, MR 2004c:11015 [11] G. H. Hardy and E. M. Wright, An introduction to the theory of numbers, fifth ed., The Clarendon Press Oxford University Press, New York, MR 81i:10002 [12] C. Hooley, On the distribution of the roots of olynomial congruences, Mathematika 11 (1964), MR 29 #1173 [13], On the greatest rime factor of a quadratic olynomial, Acta Math. 117 (1967), MR 34 #4225 [14] H. W. Lenstra, Jr., Primality testing, Mathematics and comuter science (Amsterdam, 1983), CWI Monogr., vol. 1, North-Holland, Amsterdam, 1986, MR 88b:11087 [15] A. Schinzel, On rimitive rime factors of a n b n, Proc. Cambridge Philos. Soc. 58 (1962), MR 26 #1280 [16] J. H. Silverman, Wieferich s criterion and the abc-conjecture, J. Number Theory 30 (1988), no. 2, MR 89m:11027 [17] M. Somos, Problem 1470, Crux Mathematicorum 15 (1989), 208. [18] C. Swart, Ellitic Curves and Related Sequences, PhD Thesis, University of London (2003). [19] Á. Tóth, Roots of quadratic congruences, Internat. Math. Res. Notices (2000), no. 14, MR 2001m:11159 [20] S. S. Wagstaff, Jr., Divisors of Mersenne numbers, Math. Com. 40 (1983), no. 161, MR 84j:10052 [21] G. W. Woltman, GIMPS, [22] K. Zsigmondy, Zur Theorie der Potenzreste, Monatsh. Math. 3 (1892), (GE) School of Mathematics, University of East Anglia, Norwich NR4 7TJ, UK address: g.everest@uea.ac.uk