arxiv: v2 [math.nt] 13 Dec 2013

Transcription

1 THE LEAST COMMON MULTIPLE OF RANDOM SETS OF POSITIVE INTEGERS JAVIER CILLERUELO, JUANJO RUÉ, PAULIUS ŠARKA, AND ANA ZUMALACÁRREGUI arxiv: v2 [math.nt] 13 Dec 2013 Abstract. We study the typical behavior of the least common multiple of the elements of a random subset A {1,...,n}. For example we prove that lcm{a : a A} = 2 n(1+o(1 for almost all subsets A {1,...,n}. 1. Introduction The function ψ(n = loglcm{m : 1 m n} was introduced by Chebyshev in his study on the distribution of the prime numbers. It is a well nown fact that the asymptotic relation ψ(n n is equivalent to the Prime Number Theorem, which was proved independently by J. Hadamard and C.J. de la Vallée Poussin. In the present paper, instead of considering the whole set {1,...,n}, we study the typical behavior of the quantity ψ(a := loglcm{a : a A} for a random set A in {1,...,n} when n. We define ψ( = 0. We consider two natural models. In the first one, denoted by B(n,δ, each element in A is chosen independently at random in {1,...,n} with probability δ = δ(n, typically a function of n. Theorem 1.1. If δ = δ(n < 1 and δn then ψ(a n δlog(δ 1 asymptotically almost surely in B(n,δ when n. The case δ = 1 corresponds to the classical Chebyshev function and its asymptotic estimate appears as the limiting case, as δ tends to 1, in Theorem 1.1, since lim δ 1 δlog(δ 1 = 1. When δ = 1/2 all the subsets A {1,...,n} are chosen with the same probability and Theorem 1.1 gives the following result. Corollary 1.1. For almost all sets A {1,...,n} we have that lcm{a : a A} = 2 n(1+o(1. For a given positive integer = (n, again typically a function of n, we consider the second model, where each subset of elements is chosen uniformly at random among all sets of size in {1,...,n}. We denote this model by S(n,. When δ = /n the heuristic suggests that both models are quite similar. Indeed, this is the strategy we follow to prove Theorem

2 2 J. CILLERUELO, J. RUÉ, P. ŠARKA, AND A. ZUMALACÁRREGUI Theorem 1.2. For = (n < n and we have ψ(a = log(n/ 1 /n ( 1+O(e C log almost surely in S(n, when n for some positive constant C. The case = n, which corresponds to Chebyshev s function, is also obtained as a limiting case in Theorem 1.2 in the sense that lim /n 1 log(n/ 1 /n = 1. This wor has been motivated by a result of the first author about the asymptotic behavior of ψ(a when A = A q,n := {q(m : 1 q(m n} for a quadratic polynomial q(x Z[x]. We wondered if that behavior was typical among the sets A {1,...,n} of similar size. We analyze this issue in the last section. 2. Chebyshev s function for random sets in B(n, δ. Proof of Theorem 1.1 The following lemma provides us with an explicit expression for ψ(a in terms of the Mangoldt function logp if m = p for some 1 Λ(m = 0, otherwise. Lemma 2.1. For any set of positive integers A we have ψ(a = m Λ(mI A(m, where Λ denotes the classical Von Mangoldt function and { 1 if A {m,2m,3m,...}, I A (m = 0 otherwise. Proof. We observe that for any positive integer l, the number logl can be written as logl = p l logp, where the sum is taen over all the powers of primes. Thus, using that p lcm{a : a A} if and only if A {p,2p,3p,...}, we get loglcm(a : a A = logp = (logpi A (p = Λ(mI A (m. p m p lcm(a: a A Note that if A = {1,...,n} then ψ(a = m nλ(m is the classical Chebychev function ψ(n Expectation. First of all we give an explicit expression for the expected value of the random variable X = ψ(a where A is a random set in B(n,δ. Proposition 2.1. For the random variable X = ψ(a in B(n,δ we have E(X = n δlog(δ 1 +δ ( r 1R r 1, where R(x = ψ(x x denotes the error term in the Prime Number Theorem.

3 THE LEAST COMMON MULTIPLE OF RANDOM SETS OF POSITIVE INTEGERS 3 Proof. The ambiguous case δ = 1 must be understood as the limit as δ 1, which recovers the equality ψ(n = n + R(n. In the following we assume that δ < 1. By linearity of the expectation, Lemma 2.1 clearly implies E(X = m nλ(me(i A (m. Since E(I A (m = P(A {m,2m,...} = 1 r n/m P(rm A = 1 (1 δ n/m, we obtain (1 E(X = m nλ(m ( 1 ( n/m. n We observe that n/m = r whenever r+1 < m n r, so we split the sum into intervals J r = r+1, n r ], obtaining E(X = (1 ( r Λ(m m J r r 1 = ( r 1(1 ( r ψ = δ ( r 1ψ r 1 = δn r 1 = n δlog(δ 1 ( r 1 Corollary 2.1. If δ = δ(n < 1 and δn then for some constant C > 0. r E(X = n δlog(δ 1 ψ r +1 +δ ( r 1R r 1. +δ ( r 1R r 1. ( 1+O ( e C log(δn. Proof. We estimate the absolute value of sum appearing in Proposition 2.1. For any positive integer T and using that R(y < 2y for all y > 0 we have R(n/ ( r 1 = R(n/ ( r 1 + R(n/ ( r 1 r 1 1 r T n 1 r T ( n max x n/t n log(δ 1 ( r T+1 R(n/ ( r 1 +2n ( r 1 (n/ r r r T+1 R(x ( r 1 +2n ( r 1 x r r 1 r T r T+1 ( R(x max + 2n ( T x n/t x T +1 δ Taing into account that ( T < e δt and the nown estimate max x>y R(x x e C1 logy

4 4 J. CILLERUELO, J. RUÉ, P. ŠARKA, AND A. ZUMALACÁRREGUI for the error term in the PNT, we have r 1 R(n/ ( r 1 n log(δ 1 ( e C1 log(n/t +n e δt δt. Thus we have proved that for any positive integer T we have ( 1+O (e C1 log(n/t +O E(X = n δlog(δ 1 ( e δt log(δ 1. δt We tae T δ 1 log(δn to minimize the error term. To estimate the first error term we observe that log(n/t log(δn/ log(δn log(δn, so e C1 log(n/t e C log(δn for some constant C. To bound the second error term we simply observe that δt > 1 and that log(δ 1 1 and we get a similar upper bound Variance. Proposition 2.2. For the random variable X = ψ(a in B(n,δ we have V(X δnlog 2 n. Proof. By linearity of expectation we have that V(X = E ( X 2 E 2 (X = Λ(mΛ(l(E(I A (mi A (l E(I A (me(i A (l. m,l n We observe that if Λ(mΛ(l 0 then l m, m l or (m,l = 1. Let us now study the term E(I A (mi A (l in these cases. (i If l m then (ii If (l,m = 1 then E(I A (mi A (l = 1 ( n/m. E(I A (mi A (l = 1 ( n/m ( n/l +( n/m + n/l n/ml. Both of these relations are subsumed in E(I A (mi A (l = 1 ( n/m ( n/l +( n/m + n/l n(m,l/ml. Therefore, it follows from (1 that for each term in the sum we have Λ(mΛ(l(E(I A (mi A (l E(I A (me(i A (l = Λ(mΛ(l( n/m + n/l n(m,l/ml ( 1 ( n(m,l/ml. Finally, by using the inequality 1 (1 x r rx we have and therefore: Λ(mΛ(l(E(I A (mi A (l E(I A (me(i A (l δn Λ(l Λ(m l m (m,l, V(X 2δn 1 l m n Λ(l Λ(m l m (m,l.

5 THE LEAST COMMON MULTIPLE OF RANDOM SETS OF POSITIVE INTEGERS 5 We now split the sum according to l m or (l,m = 1 and estimate each one separately. 1 l m n l m 1 l m n (l,m=1 Λ(l Λ(m l m (m,l = p n 1 i Λ(l Λ(m l m (m,l 1 l n logp logp p i p i p p n 1 i Λ(l l 1 m n ilog 2 p p i log 2 n, Λ(m log 2 n, m as we wanted to prove. We finish the proof of Theorem 1.1 by observing that V(X = o(e(x 2 when δn, so X E(X asymptotically almost surely. 3. Chebyshev s function for random sets in S(n,. Proof of Theorem 1.2 Let us consider again the random variable X = ψ(a, but in the model S(n,. From now on E (X and V (X will denote the expected value and the variance of X in this probability space. Clearly, for s = 1,2 we have E (X s = V (X = Lemma 3.1. For s = 1,2 and 1 < we have that 1 ψ s (A A = 1 (ψ(a E (X 2 A = E (X s E (X s E (X s +( s s log s n. Proof. In order to prove the lower bound it is enough to consider the case = 1. Observe that the function ψ is monotone with respect to inclusion, i.e. ψ(a {a} ψ(a for any A,{a} [n]. Using this we get A = 1 ψ s (A 1 n +1 a [n]\a Inequality then follows from ( n 1 = n (n +1(. ψ s (A {a} = (n +1 A = ψ s (A. and any partition into two sets For the second inequality we observe that for any set A ( [n] A = A A with A =, A = wehavethatψ(a ψ(a +ψ(a ψ(a +( logn. Similarly, ψ 2 (A (ψ(a +( logn 2 = ψ 2 (A +2ψ(A ( logn+( 2 log 2 n ψ 2 (A +2( log 2 n+( 2 log 2 n = ψ 2 (A +( 2 2 log 2 n.

6 6 J. CILLERUELO, J. RUÉ, P. ŠARKA, AND A. ZUMALACÁRREGUI Thus, for s = 1,2 we have Then, A = ψ s (A ψ s (A = = = and the second inequality holds. ( 1 (ψ s (A +( s s log s n A A A = 1( ( ( 1 A A A = A = A A A = ψ s (A +( s s log s n. ψ s (A + ( 1 ψ s (A 1+ A = A A A = ( 1 ( n ( n ψ s (A + A = A = ψ s (A + ( n ( s s log s n ( n ( s s log s n ( s s log s n, ( n ( s s log s n Proposition 3.1. For s = 1,2 we have that E (X s = E(X s +O( s 1/2 log s n where E(X s denotes the expectation of X s in B(n,/n and E (X s the expectation in S(n,. Proof. Observe that for s = 1,2 we have E(X s E (X s = E (X s + = E (X s + = n =0 for s = 1,2. Using Lemma 3.1 we get ( n (2 E (X s E(X s log s n =0 n =0 n =0 ( 1 n n =0 ( ( 1 n ψ s (A n n A = ( ( 1 n n n E (X s n ( n (E (X s E (X s, ( ( 1 n ( n s s. n n The sum in (2 for s = 1 is E( Y E(Y, where Y Bin(n,/n is the binomial distribution of parameters n and /n. Chauchy-Schwarz inequality for the expectation implies that this quantity is bounded by the standard deviation of the binomial distribution. n ( ( (3 1 n ( n n(/n(1 /n, n n which proves Proposition 3.1 for s = 1.

7 THE LEAST COMMON MULTIPLE OF RANDOM SETS OF POSITIVE INTEGERS 7 To estimate the sum in (2 for s = 2, we split the expression in two terms: the sum indexed by 2 and the one with > 2. We use (3 to get ( ( 1 n ( n n ( ( n ( n n n n n 2 On the other hand, >2 =0 3 3/2. ( ( 1 n ( n 2 2 n n (l l 2 l< (l+1 (l P(Y > l l 2 ( ( 1 n ( n n n where, once again, Y Bin(n,/n. Chernoff s Theorem implies that for any ǫ > 0 we have P(Y > (1+ǫ e ǫ2 /3. Applying this inequality to P(Y > l we get ( ( 1 n ( n 2 2 n n >2 (l e (l 12/3 2 e /3 3/2. l 2 The next corollary proves the first part of Theorem 1.2. Corollary 3.1. If = (n < n and then E (X = log(n/ 1 /n ( 1+O (e C log Proof. Proposition 3.1 for s = 1 and Corollary 2.1 imply that E (X = log(n/ ( ( 1+O e C log +O 1 /n ( and clearly 1/2 = o e C log when. ( 1/2 To conclude the proof of Theorem 1.2 we combine Proposition 2.2 and Proposition 3.1 to estimate the variance V (X in S(n,: V (X = E (X 2 E 2 (X = V(X+ ( E (X 2 E(X 2 +(E(X E (X(E(X+E (X ( log 2 n+ 1/2 logn (logn 3/2 log 2 n. The second assertion of Theorem 1.2 is a consequence of the estimate V (X = o ( E 2 (X when.

8 8 J. CILLERUELO, J. RUÉ, P. ŠARKA, AND A. ZUMALACÁRREGUI 3.1. The case when is constant. The case when is constant and n is not relevant for our original motivation but we give a brief analysis for the sae of completeness. In this case Fernández and Fernández [3] have proved that E (ψ(a = logn+c +o(1 where C = + =2( ( 1 ζ ( ζ(. Actually, they consider the probabilistic model with independent choices in {1,...,n}, but when is fixed it does not mae big differences because the probability of a repetition between the choices is tiny. It is easy to prove that with probability 1 o(1 we have that ψ(a logn. To see this we observe that a 1 a (a i,a 1 lcm(a 1,...,a a 1 a n, i< so i=1 loga i i< log(a i,a ψ(a logn. Now, let us note that P(a i n/logn for some i = 1,..., /logn and that P((a i,a logn d>logn P(d a i, d a d>logn 1 d < 1 2 logn. These observations imply that with probability at least 1 +( 2 logn we have that logn(1 O(loglogn/logn ψ(a logn. The analysis in the model B(n,δ when δn c can be done using again Proposition 2.1. E(ψ(A = n δlog(δ 1 +δ r<n/logn R ( r 1 +δ n/logn r n R ( r 1 We use the estimate R(x x/logx in the first sum and the estimate R(x x in the second one. We have E(ψ(A = clogn+o(1+o c ( r 1 ( +O c r 1 loglogn r r r< n log n ( clogδ = clogn+o +O(cloglogn loglogn = clogn(1+o(1. n log n r n Of course in this model we cannot expect concentration around the expectation because for example the probability that A is the empty set tends to a positive constant, P(A = e c, and then P(ψ(A = 0 e c. 4. The least common multiple of the values of a polynomial Chebyshev s function could be also generalized to ψ q (n = loglcm{q( : 1, 1 q( n} for a given polynomial q(x Z[x] and it is natural to try to obtain the asymptotic behavior for ψ q (n. Some progress has been made in this direction. While the Prime Number Theorem is equivalent to the asymptotic ψ q (n n for q(x = x, Paul Bateman noticed that the Prime Number Theorem for arithmetic progressions could be exploited to obtain the asymptotic estimate when q(x = a 1 x + a 0 is a linear polynomial and proposed it as a problem [1] in the

9 THE LEAST COMMON MULTIPLE OF RANDOM SETS OF POSITIVE INTEGERS 9 American Mathematical Monthly: ψ q (n n a 1 m φ(m 1 l m (l,m=1 where m = a 1 /(a 1,a 0. The first author [2] has extended this result to quadratic polynomials. For a given irreducible quadratic polynomial q(x = a 2 x 2 +a 1 x+a 0 with a 2 > 0 the following asymptotic estimate holds: (4 ψ q (n = 1 2 (n/a 2 1/2 log(n/a 2 +B q (n/a 2 1/2 +o(n 1/2, where the constant B q depends only on q. In the particular case of q(x = x 2 + 1, he got ψ q (n = 1 2 n1/2 logn+b q n 1/2 +o(n 1/2 with B q = γ 1 log2 2 1 l, ( 1 p 1 2 logp, p 1 p 2 where γ is the Euler constant and the sum is considered over( all odd prime numbers. It has been proved [4] that the error term in (4 for q(x = x 2 +1 is O n 1/2 (logn 4/9+ǫ for each ǫ > 0. When q(x is a reducible polynomial the behavior is, however, different. In this case it is nown (see Theorem 3 in [2] that: ψ q (n cn 1/2 where c is an explicit constant depending only on q. For example for q(x = x 2 1 the constant is c = 1. The asymptotic behavior of ψ q (n remains unnown for irreducible polynomials of higher degree. Conecture 1 (Cilleruelo [2]. Let q(x be an irreducible polynomial of degree d 3. Then (5 ψ q (n (1 1/d(n/a d 1/d log(n/a d, where a d > 0 is the coefficient of x d in q(x. For example, this conecture would imply ψ q (n 2 3 n1/3 logn for q(x = x We observe that ψ q (n = ψ(a q,n where A q,n := {q( : 1, 1 q( n}. It is natural to wonder whether for a given polynomial q(x the asymptotic E (X ψ q (n holds, when n, where = A q,n and X = ψ(a for a random set A of elements in {1,...,n}. However, consider for example the polynomials q(x = x 2 1 and q(x = x In both cases A q,n n but the asymptotic behaviors of ψ q (n are distinct: n when q(x = x 2 1 ψ q (n 1 2 nlogn when q(x = x So, what is the typical behavior of ψ(a when A n? Is it lie in the reducible case or lie in the irreducible one? Maybe neither of them represent the typical behavior of a random set. This question was the original motivation of this wor. Theorem 1.2 with = A q,n = n/a2 +O(1 gives E (X = log(n/ ( ( 1+O e C log = 1 ( 1 /n 2 (n/a 2 1/2 log(n/a 2 +o n 1/2.

10 10 J. CILLERUELO, J. RUÉ, P. ŠARKA, AND A. ZUMALACÁRREGUI This shows that, when q(x is an irreducible quadratic polynomial, the asymptotic behavior of ψ q (n coincides with ψ(a, for almost all sets of size A q,n. Theorem 1.2 also supports Conecture 1 for any q(x = a d x d + +a 0 irreducible polynomial of degree d 3. Nevertheless, there are some differences in the second term. For example, if q(x = x 2 + 1, we have ψ q (n = 1 2 n1/2 logn+b q n 1/2 +o(n 1/2, for B = On the other hand, Theorem 1.2 implies that in corresponding model S(n, with = A q,n = n 1 we have that ψ(a = 1 2 n1/2 logn+o(n 1/2 almost surely. In other words, when q(x is an irreducible quadratic polynomial, the asymptotic behavior of ψ q (n is the same that ψ(a in the corresponding model S(n,. But, the second term is not typical unless B q = 0. Probably B q 0 for any irreducible quadratic polynomial q(x but we have not found a proof. Acnowledgments: This wor was supported by grants MTM of MICINN, ICMAT Severo Ochoa proect SEV and Research Council of Lithuania Grant No. MIP- 068/2013/LSS Part of this wor was made during third author s visit at Universidad Autónoma de Madrid, he would lie to than the people of the Mathematics Department and especially Javier Cilleruelo for their warm hospitality. The authors are also grateful to Surya Ramana and Javier Cárcamo for useful comments, and especially to Ben Green whose suggestions improved the clarity and exposition of the results. References [1] P. Bateman, J. Kalb and A. Stenger, A limit involving Least Common Multiples, American Mathematical Monthly 10797, 109 (2002, no. 4, [2] J. Cilleruelo, The least common multiple of a quadratic sequence, Compositio Mathematica 147 (2011 no. 4, [3] J.L. Fernández and P. Fernández, On the probability distribution of the gcd and lcm of r-tuples of integers. Preprint. [4] J. Rué, P. Šara, A. Zumalacárregui, On the error term of the logarithm of the lcm of a quadratic sequence, Journal de Théorie des Nombres de Bordeaux 25(2013, no. 2, p

11 THE LEAST COMMON MULTIPLE OF RANDOM SETS OF POSITIVE INTEGERS 11 J. Cilleruelo: Instituto de Ciencias Matemáticas (CSIC-UAM-UC3M-UCM and Departamento de Matemáticas, Universidad Autónoma de Madrid, Madrid, Spain address: J. Rué: Institut für Mathemati, Freie Universität Berlin, Arnimallee 3-5, D Berlin, Germany address: P. Šara: Institute of Mathematics and Informatics, Aademios 4, Vilnius LT-08663, Lithuania and Department of Mathematics and Informatics, Vilnius University, Naugarduo 24, Vilnius LT-03225, Lithuania address: A. Zumalacárregui: Instituto de Ciencias Matemáticas (CSIC-UAM-UC3M-UCM and Departamento de Matemáticas, Universidad Autónoma de Madrid, Madrid, Spain address: