arxiv: v1 [math.nt] 20 Dec 2018

Transcription

1 arxiv: v1 [math.nt] 20 Dec 2018 DISTRIBUTION OF SHORT SUBSEQUENCES OF INVERSIVE CONGRUENTIAL PSEUDORANDOM NUMBER GENERATOR LÁSZLÓ MÉRAI AND IGOR E. SHPARLINSKI Abstract. In this paper we study the distribution of very short sequences of inversive congruential pseudorandom numbers modulo 2 t. We derive a new bound on exponential sums with such sequences and use it to give estimate their discrepancy. The technique we use, based the method of N. M. Korobov (1972) of estimating double Weyl sums and a fully explicit form of the Vinogradov mean value theorem due to K. Ford (2002), has never been used in this area and is very likely to find further applications. 1. Introduction 1.1. Background on the Möbius map. Let t ě 3 be an integer and write U t R t for the group of units of the residue ring R t Z{2 t Z modulo 2 t. Then #U t 2 t 1. It will often be convenient to identify elements of R t with the corresponding elements of the least residue system modulo 2 t. We fix a matrix ˆm11 m M 12 m 21 m 22 with entries from R t, which is nonsingular modulo 2. We then consider sequences generated by iterations of the Möbius map (1.1) M : x ÞÑ m 11x `m 12 m 21 x `m 22 over U t. That is for u 0 P R t we consider the trajectory (1.2) u n M pu n 1 q M n pu 0 q, n 1,2,..., generated by iterations of the Möbius map (1.1) associated with M Mathematics Subject Classification. 11K38, 11K45, 11L07. Key words and phrases. Inversive congruential pseudorandom numbers, prime powers, exponential sums, Vinogradov mean value theorem. 1

2 2 L. MÉRAI AND I. E. SHPARLINSKI Assume that the characteristic polynomial of M has two distinct eigenvalues ϑ 1 and ϑ 2 from the algebraic closure Q 2 of the field of 2- adic fractions Q 2. It is not difficult to prove by induction on n that there is an explicit representation of the form (1.3) u n γ 11ϑ n 1 `γ 12ϑ n 2 γ 21 ϑ n 1 `γ 22 ϑ n 2 with some coefficients γ ij P Q 2, i,j 1,2. Here we consider the split case when the eigenvalues ϑ 1,ϑ 2 P Z 2 are 2-adic integers, in which case, interpolating, we also have γ i,j P Z 2, i,j 1,2. It is easy to see that in this case when γ 2,1 and γ 2,2 are odd, defining g P U t by the equation g ϑ 1 {ϑ 2 in R t (recall that M is invertible in R 2 ) for the sequence generated by (1.2), the representation (1.3) we have (1.4) u n a g n b `c with some coefficients a,b,c, P R t. Furthermore, it is also easy to see that b 0 mod Motivation. The sequences (1.2) are interesting in their own rights but it has also been used as a source of pseudorandom number generation where the map is known under the name of inversive generator, for example, see [4] for the period length and [10] for the distributional properties. Moreprecisely, let τ bethemultiplicative order ofg modulo2 t. Then pu n q is a periodic sequence with period length τ, provided, that a is odd. Niederreiter and Winterhof [10], extending the results of [9] from odd prime powers to powers of 2, obtained nontrivial result for rather long segments of these sequences, namely for sequences of length N satisfying (1.5) τ ě N ě 2 p1{2`ηqt for any fixed η ą 0 and sufficiently large t. Here using very different techniques we will significantly reduce the range (1.5) and obtain results which are nontrivial provided (1.6) τ ě N ě 2 ct3{2

3 DISTRIBUTION OF INVERSIVE GENERATOR 3 for some absolute constant c ą 0. We also consider this as an opportunity to introduce new techniques into the area of pseudorandom number generation which we believe may have more applications and lead to new advances Our results. Here we establish upper bounds for the exponential sums S h pl,nq L`N 1 n L e `hu n {2 t, 1 ď N ď τ, where, as usual, we denote epzq expp2πizq and, as before, τ be the multiplicative order of g modulo 2 t, Using the method of Korobov [8] together with the use of the Vinogradov mean value theorem in the explicit form given by Ford [6], we can estimate S h pl,nq for the values N in the range (1.6). Throughout the paper we always introduce the parameter (1.7) ρ logn t which control the size of N relatively to the modulus 2 t on the logarithmic scale. Theorem 1.1. Let gcdpg,2q gcdpa,2q 1 and write Then for 2 8β ă N ď τ we have g 2 1 `w 0 2 β, gcdpw 0,2q 1. S h pl,nq ď cn 1 ηρ2 where ρis given by (1.7), for some absolute constants c,η ą 0 uniformly over all integer h with gcdph,2q 1. From a sequence pu n q defined by (1.4) we derive the inversive congruential pseudorandom numbers with modulus 2 t : u 0 {2 t,u 1 {2 t,...,u N 1 {2 t P r0,1q. The discrepancy D N of these numbers is defined by D N sup ApJ, Nq ˇ J ˇ N ˇ, JĂr0,1q where the supremum is taken over all subintervals J of r0,1q, ApN,Jq is the number of point u n {2 t in J for 0 ď n ă N, and J is the length of J. The Erdős Turán inequality (see [5, Theorem 1.21]) allows us to give upper bound on the discrepancy D N in terms of S h pnq.

4 4 L. MÉRAI AND I. E. SHPARLINSKI Theorem 1.2. Let pu n q as in Theorem 1.1 and assume that 2 32β ă N ď τ. Then we have D N ď c 0 N η 0ρ 2 where ρ is given by (1.7), for some constants c 0,η 0 ą Preparation 2.1. Notation. We recall that the notations U! V, and V " U are equivalent to the statement that the inequality U ď cv holds with some absolute constant c ą Multiplicative order of integers. The following assertion describes the order of elements modulo a power of two modulus. Lemma 2.1. Let g 1 be an odd integer and write g 2 1 `w 0 2 β, gcdpw 0,2q 1. Then for s ě β the multiplicative order τ s of g modulo 2 s is τ s 2 s β`1 and (2.1) g τs 1 `w s 2 s, gcdpw s,2q 1. Proof. Fist we note that β ě 2. If (2.1) holds for some s ě β, then by squaring it we get g 2τs 1 `w s 2 s`1 `w 2 s 22s`2 1 `w s`1 2 s`1, with w s`1 1 `w s 2 s 1 1 mod Explicit form of the Vinogradov mean value theorem. Let N k,n pmq be the number of integral solutions of the system of equations x j 1 `... `xj k yj 1 `... `yj k, 1 ď x i,y i ď M, i 1,...,k. j 1,...,n, Our application of Lemma 2.3 below rests on a version of the Vinogradov mean value theorem which gives a precise bound on N k,n pmq. We use its fully explicit version given by Ford [6, Theorem 3], which we present here in the following weakened and simplified form. Lemma 2.2. For every integer n ě 129 there exists an integer k P r2n 2,4n 2 s such that for any integer M ě 1 we have N k,n pmq ď n 3n3 M 2k 0.499n2.

5 DISTRIBUTION OF INVERSIVE GENERATOR 5 We note that the recent striking advances in the Vinogradov mean value theorem due to Bourgain, Demeter and Guth [3] and Wooley [11] are not suitable for our purposes here as they contain implicit constants that depend on k and n, while in our approach k and n grow together with M Double exponential sums with polynomials. Our main tool to bound the exponential sum S h pnq is the following result of Korobov [8, Lemma 3]. Lemma 2.3. Assume that ˇ al ˇαl q l ˇ ď 1 and gcdpa ql 2 l,q l q 1, for some real α l and integer a l,q l, l 1,...,n. Then for the sum M S e pα 1 xy `... `α n x n y n q we have where x,y 1 S 2k2 ď `64k 2 logp3qq n{2 M 4k2 2k N k,n pmq nź min "M l,? * q l `? Ml, ql l 1 Q maxtq l : 1 ď l ď nu. We also need the following simple result which allows to reduce single sums to double sums. Lemma 2.4. Let f : R Ñ R be an arbitrary function. Then for any integers M,N and a we have N 1 epf pxqq ˇ ˇ ď 1 N M epfpx `ayzqq ˇ ˇ `2aM2. x 0 M 2 x 1 y,z 1 Proof. Examining the non-overlapping parts of the sums below, we see that for any positive integers y and z N 1 N 1 epfpxqq epfpx `ayzqq ˇ ˇ ď 2ayz ď 2aM2. x 0 N 1 Henceˇˇˇ ˇ ˇM2 epfpxqq x 0 x 0 M N 1 y,z 1 x 0 epfpx `ayzqq ˇ ď 2ayz ď 2aM2.

6 6 L. MÉRAI AND I. E. SHPARLINSKI Changing the order of summation and using the triangle inequality, and the result follows Sums of binomial coefficients. We need results of certain sums of binomial coefficients. The first ones are immediate and we leave the proof for the reader. Lemma 2.5. For any integer n ě 1 (1) for any integers k ď n we have n ˆi ˆn `1 ; k k `1 i k (2) for any polynomial P pxq P ZrXs of degree degp ă n we have n p 1q jˆn P pjq 0. j j 0 Lemma 2.6. For any n,k with k ď n we have Proof. As l 1`...`l k n l 1,...,l k ě1 k n! l 1!...l k! l 1`...`l k n i 0 p 1q k iˆk n! l 1!...l k! kn i n. i the result follows directly from the inclusion exclusion principle. 3. Proofs of the main results 3.1. Proof of Theorem 1.1. As a ag L u n`l g n`l b `c `c, gn bg L we can assume, that L 0 and we put S h p0,nq S h pnq. We can also assume, that a 1 and c 0. Finally we assume, that N ě 2 6t1{2 since otherwise the result is trivial. Define r tlog2 logn ρ 1 log2

7 DISTRIBUTION OF INVERSIVE GENERATOR 7 where ρ is given by (1.7). First assume, that and put Then s Z ^ t 4r if N is large enough. Indeed, r ě 129 and κ R V t 1. s s ą β, 2 s ď N 1{4, r ď κ ă s, s ě t logn 1 4r 4log2 1 ě 2β 1 ą β and 2s ď 2 t 4r N 1{4. Moreover, κ ě t 1 ě 4r 1 ě r s and κ ď t s ď plognq2 36plog2q 2 s t2 36r 2 s ď s. Let τ s be the order of g modulo 2 s. As s ą β, g τs 1 `w 2 s with gcdpw,2q 1 by Lemma 2.1. Clearly, for all even x, we have 1 1 x 1 `x`... `x2t 1 mod 2 t, thus Define u n τs 1 2t 1 b g 1 n τs 1 p1 b`g nτs q F κ pnq Then 2 t 1 p1 b`p1 `w 2 s q n q l κ l 1 l 0 2 t 1 l 0 2 b` 2t 1 pw 2 s q l n i 1 j j 0 ν 1 l 0 ˆn pw i l 2 s q mod 2 t. i ˆj p2 bq j ν ν u n τs F κ pnq mod 2 t. i 1`...`i ν l i 1,...,i νě1 p1 b`g n τs q l ˆn i 1... ˆn i ν.

8 8 L. MÉRAI AND I. E. SHPARLINSKI The expression κ!f κ pnq is a polynomial of 2 s n of degree at most κ. Thus we can define the integers a 0,...,a κ by κ κ!f κ pnq a l 2 ls n l. l 1 Then the coefficients satisfy a l κ! 2 t 1 l! wl j j 1 ν 1 ˆj p2 bq j ν ν i 1`...`i ν l i 1,...,i νě1 l! i 1!...i ν! mod 2s. We have v 2 pa l q v 2 pκ!{l!q. Indeed, as w is odd and b is even, by Lemmas 2.6 and 2.5 we get 2 t 1 j j 1 ν 1 ˆj p2 bq j ν ν l j 1 i 1`...`i j l i 1,...,i j ě1 l l p 1q i i l i 0 l`1 i 1 j i i 1`...`i ν l i 1,...,i νě1 l! l i 1!...i j! ˆj i p 1q iˆl `1 i l! i 1!...i ν! j j 1i 0 p 1q j iˆj l p 1q i i lˆl `1 i `1 i 0 ˆl `1 pi 1q l 0 i i l p 1q l 1 mod 2 (we note that several last congruences are actually equations). Write ω l v 2 pa l q. Then Y κ ] Y κ ] ω l ď v 2 pκ!q ď ` `... ă κ for l ă κ 2 4 and ω κ 0.

9 DISTRIBUTION OF INVERSIVE GENERATOR 9 To conclude the proof observe, that by Lemma 2.4 we have S h pnq ď 1 N 1 2 s ˆh e 2 2s ˇ 2 tu n`τ sxy ˇˇˇˇ ˇ `2τ s2 2s n 0 x,y 1 ď 1 N 1 2 s ˆh e ˇˇˇˇ g n 2 2s ˇ 2 t gτsxy bg n ˇ `23s n 0 x,y 1 ď 1 N 1 2 s ˆhg n pa 1 2 s xy `... `a κ 2 κs pxyq κ q e ˇˇˇˇ 2 2s ˇ κ!2 t ˇ `N3{4, n 0 x,y 1 where the coefficients a l a pbg n l q for l 1,...,κ, are determined as above with bg n instead of b. Write hg n a l 2 ls r l gcdpr κ!2 t l,q l q 1, l 1,...,κ q l with (3.1) 2 t ls ω l ď q l ď κ!2 t ls ω l l 1,...,κ. Then (3.2) S h pnq ď 1 2 2s where Put N 1 n 0 ˇ 2 s x,y 1 e pf n px,yqqˇ f n py,yq r 1 q 1 xy `... ` rκ q κ pzyq κ. σ n 2 s x,y 1 e pf n px,yqq. ˇ `N3{4 For κ, there exists a k P r2κ 2,4κ 2 s such that for N k,κ we have the bound of Lemma 2.2 (with κ instead of n). Then by Lemma 2.3 we have (3.3) where σ n 2k2 ď `64k 2 logp3qq κ{2 2 p4k2 2kqs N k,κ p2 s q κź min "2 ls,? * q l `? 2ls, ql l 1 (3.4) logq ď logpκ!2 t q ď κtlogp2κq. by (3.1).

10 10 L. MÉRAI AND I. E. SHPARLINSKI By the choice of κ we have sκ ă t ď spκ `1q. As ω l ď κ ď s, under κ `1 ď l ă κ 2 we have by (3.1) thus Whence κź (3.5) l 1 q l ď κ!2 spκ`1 lq ď κ!2 ls and q l ą 2 spκ 1 lq? 1 ql 1 `κ!? ` ď? ď κ κ 2 s 2 pκ 1 lq. ql 2ls ql min "2 ls,? * q l `? 2ls ql By Lemma 2.2 we have (3.6) 2 sκpκ`1q{2 κź l 1 ď 2 sκpκ`1q{2 ź " 1 min 1,? ` ql κ κ 2 spκ 1 lq{2 κ 2 ălăκ ď κ κ2 2 sκpκ`1q{2 spκ 2qpκ 4q{16. N k,κ p2 s q ď κ 3κ3 2 2ks 0.499κ2 s Combining (3.3), (3.4), (3.5) and (3.6) we have? ql σ n 2k2 ď `65κ 3 logpκqt κ{2 κ 4κ 3 2 4k2 s`sκpκ`1q{2 spκ 2qpκ 4q{ κ 2 s then σ n! t 1{p16κ3q 2 2s s{p32770κ2q. Since tκ 2 ă p t s q3 s ă p6rq 3 s, then Moreover 2 s{κ2 N rs{ptκ2q ą N 1{p216r2q. t 1{κ3 ď N logt{p129r2 lognq ď N loglogn{p387r2 lognq. Whence σ n! 2 2s N ηρ2, for some η ą 0 if N is large enough. Thus by (3.2) we have S h pnq ď 1 N 1 σ 2 2s n `N 3{4! N 1 ηρ2 `N 3{4! N 1 η{r2 n 0 which gives the result for r ě 129. If r ă 129, define N 0 X 2 t{129\ ρ 0 logn 0 t log2 129 `Op1{tq. 2 ls *

11 As N ď τ ă 2 t, we have (3.7) Then DISTRIBUTION OF INVERSIVE GENERATOR 11 S h pnq ď logn 0 logn ą ďkăN{N 0 ˇ pk`1qn 0 1 n kn 0 ephu n {2 t qˇ ˇ. Applying the previous argument to the inner sums, we get S h pnq! N N 0 N 1 ηρ2 0 0! N ηρ 2 0 by (3.7). Thus replacing η to η{129 3, we conclude the proof Proof of Theorem 1.2. By the Erdős-Turán inequality, see [5] for any integer H ě 1 we have (3.8) D N! 1 H ` 2 H 1 N h S hp0,nq. Define H h 1 Z ^ τt?. N For a given 1 ď h ď H, write h 2 d j with odd j and d ď log 2 H. Then consider the sequence pu n q modulo 2 t d. Then clearly S h p0,nq S d,j p0,nq. where S d,j p0,nq is defined as S j p0,nq, however with respect to the modulus 2 t d. By the above choice of parameters, we have (3.9) t d ě t log 2 H ě 1 2 log 2N `β ą 17β by Lemma 2.1, thus (3.10) τ t d 2 t d β. Using (3.8), we have (3.11) D N! 1 H ` 1 N! 1 H ` 1 N H 1 h S hp0,nq d j S d,jp0,nq. 1ďjďH{2 d j odd h 1 0ďdďlog 2 H

12 12 L. MÉRAI AND I. E. SHPARLINSKI For fixed d and j put R N N d τ t d V and K d N N d τ t d. Then (3.12) S d,j p0,nq ď As N d 2 i 0 If K d ă 2 8β, we use the trivial estimate by (3.9), we get S d,j piτ t d,τ t d q ` S d,j ppn d 1qτ t d,k d q S d,j pn d τ t d,k d q ď K d ă 2 8β. 8β ă 1 pt d βq 2 (3.13) S d,j pn d τ t d,k d q ď τ 1 ηpt dq 2 plogτ t d q 2 t d. If K d ě 2 8β, then as K d ď τ t d we also have (3.13) by Theorem 1.1. Thus by (3.12) we have S d,j p0,nq! N d τ 1 ηpt dq 2 plogτ t d q 2 t d! N 1 ηpt dq 2 plogτ t d q 2 {logn. By (3.9) and (3.10) we have plogτ t d q 3 pt d βq3 pt d βq3 ě ě 1 lognpt dq2 lognpt dq2 lognt 2 8 whence Then by (3.11), D N! 1 H ` 0ďdďlog 2 H S d,j p0,nq! N 1 ηρ2 { d j N ηρ2 {8 1ďjďH{2 d j odd ˆlogN t 2 ρ 2 {8,! 2 pt βq{2 `N ηρ2 {8 logh! 1 t `N ηρ2 {8 logh! N ηρ2 {16 if N is large enough.

13 DISTRIBUTION OF INVERSIVE GENERATOR Comments Wenotethatanextensionofourresultstothecaseofsequences (1.4) modulo prime powers p t with a prime p ě 3 is immediate and can be achieved at the cost of merely typographical changes. It is certainly natural to study the multidimensional distribution of the sequence generated by (1.2), that is, the s-dimensional vectors pu n,...,u n`s 1 q, n 1,...,N. Our method is capable to address this problem, however investigating the 2-divisibility of the corresponding polynomial coefficients which is an important part of our argument in Section 3.1 is more difficult and may require new arguments. We also use this as an opportunity to pose a question about studying short segments of the inversive generator modulo a large prime p. While results of Bourgain [1,2] give a non-trivial bound on exponential sums for very short segments of sequence ag n mod p, n 1,...,N, see also [7, Corollary 4.2], their analogues for even simplest rational expressions like 1{pg n bq mod p are not known. Obtaining such results beyond the easy range N ě p 1{2`ε (with any fixed ε ą 0) is apparently a difficult questions requiring new ideas. Acknowledgement During the preparation of this wok L. M. was partially supported by the Austrian Science Fund FWF Projects P30405 and I. S. by the Australian Research Council Grants DP and DP References [1] J. Bourgain, Multilinear exponential sums in prime fields under optimal entropy condition on the sources, Geom. and Funct. Anal., 18(2009), [2] J. Bourgain, On exponential sums in finite fields, Bolyai Soc. Math. Stud., 21, János Bolyai Math. Soc., Budapest, 2010, [3] J. Bourgain, C. Demeter and L. Guth, Proof of the main conjecture in Vinogradov s mean value theorem for degrees higher than three, Ann. Math., 184 (2016), [4] W.-S. Chou, The period lengths of inversive congruential recursions, Acta Arith., 73 (1995), [5] M. Drmota and R. F. Tichy, Sequences, discrepancies and applications, Springer-Verlag, Berlin, , 11 [6] K. Ford, Vinogradov s integral and bounds for the Riemann zeta function, Proc. London Math. Soc., 85 (2002), , 4

14 14 L. MÉRAI AND I. E. SHPARLINSKI [7] M. Z. Garaev, Sums and products of sets and estimates of rational trigonometric sums in fields of prime order, Uspekhi Mat. Nauk, 65 (2010), no.4, 5 66 (in Russian), translated in Russian Math. Surveys, 65 (2010), [8] N. M. Korobov, The distribution of digits in periodic fractions, Matem. Sbornik, 89 (1972), (in Russian), translated in Math. USSR-Sb., 18 (1974), , 5 [9] H. Niederreiter and I. E. Shparlinski, Exponential sums and the distribution of inversive congruential pseudorandom numbers with prime-power modulus, Acta Arith., 92 (2000), [10] H. Niederreiter and A. Winterhof, Exponential sums and the distribution of inversive congruential pseudorandom numbers with power of two modulus, Int. J. Number Theory, 1 (2005), [11] T. D. Wooley, The cubic case of the main conjecture in Vinogradov s mean value theorem, Adv. in Math., 294 (2016), L.M.: Johann Radon Institute for Computational and Applied Mathematics, Austrian Academy of Sciences, Altenberger Straße 69, A Linz, Austria address: laszlo.merai@oeaw.ac.at I.E.S.: School of Mathematics and Statistics, University of New South Wales. Sydney, NSW 2052, Australia address: igor.shparlinski@unsw.edu.au