arxiv: v1 [math.nt] 19 Oct 2010

Transcription

1 ASYMPTOTIC NORMALITY OF ADDITIVE FUNCTIONS ON POLYNOMIAL SEQUENCES IN CANONICAL NUMBER SYSTEMS MANFRED G. MADRITSCH AND ATTILA PETHŐ arxiv: v [math.nt] 9 Oct 200 Abstract. The objective of this paper is the study of functions which only act on the digits of an expansion. In particular, we are interested in the asymptotic distribution of the values of these functions. The presented result is an extension and generalization of a result of Bassily and Kátai to number systems defined in a quotient ring of the ring of polynomials over the integers.. Introduction In this paper we investigate the asymptotic behaviour of q-additive functions. But before we start we need an idea of additive functions and the number systems they are living in. Note that a function f is said to be q-additive if it acts only on the q-adic digits, i.e., f0 = 0 and fn = fa h nq h for n = a h nq h, where a h n N := 0,...,q } are the digits of the q-adic expansion of n. One of the first results dealing with the asymptotic behavior of such a q-additive function is the following, which is due to Bassily and Kátai [3]. Theorem. Let f be a q-additive function such that faq h = O as h and a N. Furthermore let m h,q := faq h, σh,q 2 := f 2 aq h m 2 q q h,q, and M q x := N N m h,q, Dq 2 x = with N = [log q x]. Assume that D q x/logx /3 as x and let P be a polynomial with integer coefficients, degree d and positive leading term. Then, as x, x # n < x fpn M q x d } D q x d < y y exp x 2 dx. 2π A first step towards a generalization of this concept is based on number systems living in an order in an algebraic number field. Definition.. Let R be an integral domain, b R, and N = n,...,n m } Z. Then we call the pair b,n a number system in R if every g R admits a unique and finite representation of the form g = a h gb h with a h g N and a h g 0 if h 0. We call b the base and N the set of digits. If N = N 0 = 0,,...,m} for m then we call the pair b,n a canonical number system. σ 2 h,q Date: October 20, Mathematics Subject Classification. K6 A63,60F05. Key words and phrases. additive functions, canonical number systems, exponential sums.

2 2 M. G. MADRITSCH AND A. PETHŐ When extending the number system to the complex plane one has to face effects such as amenability, i.e., there may exist two or more different expansions of one number. In fact, one can construct a graph the connection graph which characterizes all the amenable expansions. This has been done by Müller et al. [6] with a direct approach and by Scheicher and Thuswaldner [8] consideration of the odometer. A different view on digits in number systems is done by normal numbers. These are numbers in which expansion every possible block occurs asymptotically equally often. Constructions of such numbers have been considered by Dumont et al. [5] and the first author in [3, 4] In this paper we mainly concentrate on additive functions. Thus we define additive functions in these number systems as follows. Definition.2. Let b,n be a number system in the integral domain R. A function f is called b-additive if f0 = 0 and fg = h 0fa h gb h for g = a h gb h. The simplest version of an additive function is the sum-of-digits function s b defined by s b g := h 0a k g. The result by Bassily and Kátai was first generalized to number systems in the Gaussian integers by Gittenberger and Thuswaldner [7] who gained the following Theorem. Let b Z[i] and b,n be a canonical number system in Z[i]. Let f be a b-additive function such that fab h = O as h and a N. Furthermore let m h,b := fab h, σh,b 2 Nb := f 2 ab h m 2 h,b Nb, and M b x := L m h,b, Dbx 2 = with N the norm of an element over Q and L = [log Nb x]. Assume that D b x/logx /3 as x and let P be a polynomial of degree d with coefficients in Z[i]. Then, as N, #z Z[i] Nz < N} # Nz < N fpz M b N d D b N d L σ 2 h,b } < y y exp x 2 dx, 2π where N is the norm of an element over Q and z runs over the Gaussian integers. This build the base for further considerations of b-additive functions in algebraic number fields in general. Therefore let K = Qβ be an algebraic number field and denote by O K its ring of integers aka its maximal order. Furthermore let β O K then we set R = Z[β] to be an order in K. We now want to analyze additive functions for number systems in R. We need some more parameters in order to successfully generalize the theorem from above. Thus let K l l r be the real conjugates of K, while K m and K m+r2 r < m r +r 2 are the pairs of complex conjugates of K, where r +2r 2 = n. For γ K we denote by γ i i n the conjugates of γ. In order to extend the term of conjugation to the completion K of K we define for γ j K and x j R j n λ = j n x jγ j and λ i := j n x jγ i j. Next we have to guarantee that we choose the increasing set for our asymptotic distribution. In the integer case we had the logarithm of the value, since the length of expansion growth with the logarithm. Since R is of dimension n we need a way to enlarge the area under consideration such that the expansion growth also in a smooth way. This is motivated by the following

3 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS 3 Lemma. [2, Theorem]. Let lγ be the length of the expansion of γ to the base b. Then log γ i lγ max i n log b i C. Therefore we define RT,...,T r to be the set } RT,...,T r := λ R : λ i. Ti, i r. Now we use Lemma. to bound the area RT,...,T n such that we reach all elements of a certain length. Thus for a fixed T we set T i for i n such that logt i = logt log b i n.2 log Nb. Furthermore we will write for short RT := RT,...,T r with T i as in.2. Finally one can extend the definition of a number system also for negative powers of b. Then for γ K such that γ = a h b h with a h N h= we call h γ := h b a h and γ} := a h b h h the integer part and fractional part of γ, respectively. With all these tools we now can state the generalization of the theorem of Bassily and Kátai to arbitrary number fields. Theorem.2 [5]. Let b,n be a number system in R and f be a b-additive function such that fab h = O as h and a N. Furthermore let m h,b := fab h, σh,b 2 := f 2 ab h m 2 Nb Nb h,b, and M b x := L m h,q, Dbx 2 = with L = [log Nb x]. Assume that there exists an ε > 0 such that D b x/logx ε as x and let P K[X] be a polynomial of degree d. Then, as T let T i be as in.2, #RT # z RT f Pz M b T d D b T d } < y L 2. Definitions and result σ 2 h,q 2π y exp x 2 dx. The objective of this paper are generalizations of number systems to quotient rings of the ring of polynomials over the integers. Our aim is to extend Theorem.2 to such rings. To formulate our results we have to introduce the relevant notions. In particular we use the following definition in order to describe number systems in quotient rings of the ring of polynomials over the integers. Definition 2.. Let p Z[X] be monic of degree n and let N be a subset of Z. The pair p,n is called a number system if for every g Z[X] \ 0} there exist unique l N and a i N,h = 0,...,l;a l 0 such that 2. g a h gx h mod p. In this case a i are called the digits and l = la the length of the representation.

4 4 M. G. MADRITSCH AND A. PETHŐ Thisconceptwasintroducedin [7] andwasstudied amongothersin[, 2,, 2]. It wasproved in [2], that N must be a complete residue system modulo p0 including 0 and the zeroes of p are lying outside or on the unit circle. However, following the argument of the proof of Theorem 6. of [7], which dealt with the case p square free, one can prove that non of the zeroes of p are lying on the unit circle. If p is irreducible then we may replace X by one of the roots β of p. Then we are in the case of Z[X]/p = Z[β] being an integral domain in an algebraic number field cf. Section. Then we may also denote the number system by the pairβ,n instead of p,n. For example, let q 2 be a positive integer, then p,n with p = X q gives a number system in Z, which corresponds to the number systems q,n. Furthermore for n a positive integer and p = X 2 +2nX+n 2 + we get number systems in Z[i]. Now we want to come back to these more general number systems and consider additive functions within them. Definition 2.2. Let p,n be a number system. A function f is called additive if f0 = 0 and fg fa h gx h mod p for g a h gx h mod p. Since we have defined the analogues of number systems and additive functions to the definitions for number fields above, we now need to extend the length estimation of Lemma. in order to successfully state the result. But before we start we need a little linear algebra. We fix a number system p,n and factor p by p := t i= with p i Z[X] irreducible and degp i = n i. Furthermore we denote by β ik the roots of p i for i =,...,t and k =,...,n i. Then we define by R := Z[X]/p = t i= R i p mi i with R i = Z[X]/p mi i for i =,...,t the Z-module under consideration and in the same manner by t K := Q[X]/p = K i with K i = Q[X]/p mi i i= for i =,...,t the corresponding vector space. Finally we denote by K the completion of K according to the usual Euclidean distance. Obviously R is a free Z-module of rank n. Let λ : R R be a linear mapping and z,...,z n } be any basis of R. Then n λz j = a ij z i j =,...,n i= with a ij Z. The matrix Mλ = a ij is called the matrix of λ with respect to the basis z,...,z n }. For an element r R we define by λ r : R R the mapping of multiplication by r; that is λ r z = rz for every z R. Then we define the norm Nr and the trace Trr of an element r R as the determinant and the trace of Mλ r, respectively, i.e., Nr := detmλ r, Trr := TrMλ r. Note that these are unique despite of the used basis z,...,z n }. We can canonically extend these notions to K and K by everywhere replacing Z by Q and R, respectively. In the following we will need parameters which help us bounding the length of the expansion of an element g R. Therefore let g Z[X] be a polynomial, then we put B ijk g := dj g i =,...,t;j =,...,m i ;k =,...,n i. X=βik dx j

5 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS 5 In connection with these values we define the house function H as Hg := max t m i n i max max B ijkg. i= j= k= We want to investigate the elements with bounded maximum length of expansion. To this end we need a proposition which estimates the length of expansion in connection with properties of the number itself. The proof of this proposition will be presented in the following section. Proposition 2.. Assume that p,n is a number system. Let N = max a : a N} and we set } log Bijk g Mg := max : i =,...,t;j =,...,m i ;k =,...,n i. log β ik If g Z[X] is of degree at most n, then for any ε > 0 there exists L = Lε such that if lg > L then 2.2 lg Mg C. This provides us with an estimation for the length of the expansion and motivates us to look at subsets of R where the absolute values B ijk are bounded. For a vector T := T,...,T n = T,...,T n,t 2,...,T,mi,n,T 2,...,T t,mt,n t we denote by RT := g R : B ijk g T ijk } R i T := g R i : B ijk g T ijk }. We want to let the length of expansion to smoothly increase. Therefore we fix a T and set T ijk for i =,...,t, j =,...,m i, k =,...,n i such that log β ik n 2.5 logt ijk = logt log t i= k= β ik mi. Remark that T ijk is independent from j, which will be important in Lemma 4.2. In view of Proposition2.wegetthattheexpansionsoftheelementsinRTalmosthavethesamemaximum length. If not stated otherwise we denote by T the vector T,...,T t,mt,n t where the T ijk are as in 2.5. Since X is an invertible element in K we may extend the definition of a number system for negative powers of X. Then for γ K such that γ = a h X h with a h N we call γ := h= ni a h X h and γ} := h= the integer part and fractional part of γ, respectively. Now we have collected all the tools to state our main result. a h X h Theorem 2.2. Let p,n be a number system and f be an additive function such that fax h = O as h and a N. Furthermore let m h := fax h, σh 2 N := f 2 ax h m 2 h N, and Mx := L m h, D 2 x := L σh, 2 where L = log p0 x. Assume that there exists an ε > 0 such that Dx/logx ε as x and let P K[Y] be a polynomial of degree d. Then, as T let T ijk be as in 2.5, } #RT # z RT : f Pz MTd DT d < y y exp x 2 dx. 2π

6 6 M. G. MADRITSCH AND A. PETHŐ Our theorem shows that the distribution properties of patterns in the sequence of digits depend neither on the polynomial p nor on the quotient ring R. They are intimate properties of the backward division algorithm defined in [7]. We will show the main theorem in several steps. In the following section we will show properties of number systems which we need on the one hand to estimate the length of expansion and on the other hand to provide us with an Urysohn function, that helps us counting the occurrences of a fixed pattern of digits in the expansion. 2 Equipped with these tools we will estimate the exponential sums occurring in the proof in Section 4. Therefore we need to split the module R up into its components and consider each of them separately. We also show that we may neglect the nilpotent elements. 3 Now we take a closer look at the Urysohn function, which will count the occurences of our pattern in the expansion, and estimate the number of hits of the border of this function in Section 5. In particular, we count the number of hits of the area, where the function value lies between 0 and as this area corresponds to the error term. 4 In Section 6 we will show that any chosen patterns of digit and position occurs uniformly in the expansions. This will be our central tool in the proof of Theorem Finally we draw all the thinks together. The main idea here is to use the growth rate of the deviation together with the Fréchet-Shohat Theorem to cut of the head and the tail of the expansion. Then an application of the central Proposition 5. and juxtaposition of the moments will prove the result. 3. Number system properties In this section we want to show two properties we need in the sequel. The first deals with the above mentioned estimation of the length of an expansion Proposition 2.. We will need this result in order to justify our choice of T as in 2.5. Secondly we construct the Urysohn function for indicating the elements starting with a certain digit. The main idea is to embed the elements of R in R n and to use the properties of matrix number systems in this field. We start with the Proof of Proposition 2.. In the proof we combine ideas from [2] and [2]. We may assume g 0. As p,n} is a number system there are l = la and a h N for h = 0,...,l; a l 0 such that g a h X h mod p, i.e., g = a h X h +rp with a polynomial r Z[X]. For j this implies 3. d j g dx j = h=j h! h j +! a hx h j+ + j j s s=0 d s r dx s d j s p dx j s. Consider a zero β ik of p, which has multiplicity m i. As we noticed in the Introduction, the argument of the proof of Theorem 6.. of [7] allows to prove that β ik > for all i =,...,t;k =,...,n i. Inserting β ik into 3. we obtain B ijk g = dj g = X=βik dx j h=j h! h j +! a hβ h j+ ik

7 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS 7 for i =,...,t and j =,...,m i. This implies by taking absolute value h! B ijk g N h j +! β ik h j+ h=j l! N l j +! βl j+ N lj β ik l β ik, k h=j hh...h j + ll...l j + β ik h l which verifies the lower bound for l, because β ik >. Now we turn to prove the upper bound. Denote by V = V p the following mapping: for g Z[X] of degree at most n choose an a N such that g0 a mod p0. Such an a exists by Theorem 6. of [7]. Putting q = g0 a p0, let Vg = g q p a X. Obviously Vg Z[X] and has degree at most n, thus V can be iterated. Moreover we have u 3.2 g a h X h +X u+ V u+ g mod p with a h N for h = 0,...,u. Choose u the largest integer satisfying B ijk g uj β ik u j+ β ik for all i =,...,t, j =,...,m i and k =,...,n i. Then u +ε/2ma. Proceeding like in the previous case we get B ijk g = dj g u h! dx j = X=βik h j +! a hβ h j ik + h=j j j u+! s u+ s! βu+ s ik s=0 d j s V u+ g dx j s. X=βik By its definition V u+ g has integer coefficients. Dividing the last equation by u+! u+ j! βu+ j ik and consider the obtained equations for i =,...,t and k =,...,n i and for fixed i and k for j =,...,m i successively, then using the choice of u and that β ik > we conclude that d j s V u+ g dx j s < c, X=βik wherecisaconstantdependingonlyonn aswellasthesizeandthemultiplicitiesofthezeroesofp. These can be considered as n inequalities for the n unknown coefficients of V u+ g. Furthermore the determinant of the coefficient matrix is not zero c.f. [2]. Thus the solutions are bounded. As they are integers there are only finitely many possibilities for V u+ g. As p,n} is a number system, V u+ g has a representation, which length is bounded by a constant, say c, which depends only on N as well as the size and the multiplicities of the zeroes of p. Thus lg u+c +εmg and the proposition is proved. Nowweturnourattentionbacktothecountingofthenumbersandinparticulartotheconstruction of the Urysohn function. In order to properly count the elements we need the fundamental domain, which is defined as the set of all numbers whose integer part is zero. Since this is not so easy to define in this context we want to consider its embedding in R n. The main idea is to use the corresponding matrix of the polynomial p and to use properties of matrix number systems. This idea essentially goes back to Gröchenig and Haas [8]. The following definitions are standard in that area and we mainly follow Gittenberger and Thuswaldner [7] and Madritsch [5]. Wenotethatifp,NisanumbersystemthenX isaintegralpowerbaseofk, i.e.,,x,...,x n } is an R-basis for K. Thus we define the embedding φ by φ : K R n, a +a 2 X + +a n X n a,...,a n.

8 8 M. G. MADRITSCH AND A. PETHŐ Now let p = b n X n + +b X +b 0. Then we define the corresponding matrix B by 0 0 b B = b n One easily checks that φx A = B φa. Since B is invertible we can extend the definition of φ by setting for an integer h 3.4 φx h A := B h φa. By this we define the embedded fundamental domain by F := z Rn z = B h a h,a h φn. h Following Gröchenig and Haas [8] we get that λf +g F +g 2 = 0 for every g,g 2 Z n with g g 2, where λ denotes the n dimensional Lebesgue measure. Thus B,φN is a matrix number system and a so called just touching covering system. Therefore we are allowed to apply the results of the paper by Müller et al. [6]. We now follow the lines of Madritsch [5] where the ideas of Gittenberger and Thuswaldner [7] were combined with the results of Kátai and Környei [0] and Müller et al. [6]. Our main interest is the fundamental domain consisting of all numbers whose first digit equals a N, i.e., F a = B F +φa. Imitating the proof of Lemma 3. of [7] we get the following. Lemma 3.. For all a N and all v N there exist a µ < detb and an axe-parallel tube P v,a with the following properties: F a P v,a for all v N, the Lebesgue measure of P v,a is a O µ v detb, v P v,a consists of Oµ v axe-parallel rectangles, each of which has Lebesgue measure O detb v, where λ denotes the Lebesgue measure. As in the proof of Lemma 3. of [7] we can construct for each pair v,a an axe-parallel polygon Π v,a and the corresponding tube P v,a := z R n z Πv,a 2c p detb v}, where c p is an arbitrary constant. Furthermore we denote by I v,a the set of all points inside Π v,a. Now we define our Urysohn function u a by 3.5 where 3.6 u a x,...,x n = κ n κ 2 κ 2 κ 2 κ 2 ψ a x +y,...,x n +y n dy dy n, κ := 2c u detb v

9 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS 9 with c u a constant and if x,...,x n I v,a ψ a x,...,x n = 2 if x,...,x n Π v,a 0 otherwise. Thusu a isthe desiredurysohnfunction whichequalsforz I v,a \P v,a, 0forz R n \I v,a P v,a, and linear interpolation in between. We now do a Fourier transform of u a and estimate the coefficients in the same way as in Lemma 3.2 of [7]. Lemma 3.2. Let u a x,...,x n = m,...,m n Z n c m,...,m n em x + +m n x n be the Fourier series of u a. Then the Fourier coefficients c m,...,m n can be estimated by with c 0,...,0 = detb, c m,...,m n µ v rm i = κm i m i 0, m i = 0. n i= rm i 4. Estimation of the Weyl Sum Before we continue with the estimation of the number of points inside the fundamental domain and those hitting the boarder, we want to estimate the exponential sums, which will occur in the following sections. In particular we want to prove the following. Proposition 4.. Let T 0 and T ijk as in 2.5. Let L be the maximum length of the b-adic expansion of z RT and let C and C 2 be sufficiently large constants. Furthermore let l,...,l h be positions and h,...,h h be corresponding n-dimensional vectors. If 4. C logl l < l 2 < < l h dl C 2 logl and 4.2 z RT h r logt σ for r h, then we have h e hr,b lr φpz T n logt tσ0 r= where σ 0 depends on σ, C and C 2. Our main idea consists in several steps. First we will split the ring R up into the R i and consider each of them separately. Then we distinguish two cases according to whether m i > or not. The latter reduces to an estimation of the sum in an algebraic number field. Whereas for the case of m i > we have to deal with nilpotent elements. Therefore we divide R i into the radical and the nilpotent elements. Thus we define R i as R i := Z[X]/p i and Ri T := g R } 4.3 i : B ijk g T ijk and the set N i to be the nilpotent elments, i.e., 4.4 N i := g R i : g 0 mod p i } and N i T := g R i T : g 0 mod p i }. But before we start with the proof we need to show that the estimation is good compared with the trivial one. Thus we will show the following.

10 0 M. G. MADRITSCH AND A. PETHŐ Lemma 4.2. Let T ijk for i =,...,t, j =,...,m i, k =,...,n i be positive reals. Then t #RT = #R i T, i= ni #R i T = c i k= ni #N i T = c i k= T ik mi +O T mini 0, mi T ik +O where T 0 = max t max,t i T i,mi,n i i= and the constants c i will be defined in Lemma 4.3. T mi ni 0 m i n i Proof. The first assertion follows immediately from the definition of RT. Since the R i are independent we fix an i and focus on R i T. Obviously we have that 4.5 R i = Z[X]/p mi i = Z[X]/p i mi. Thus we concentrate on Z[X]/p i Z[X] which is an order in a number field K i of degree n i over Q. For γ K i let γ l l r be the real conjugates and γ m and γ m+r2 r + m r +r 2 be the pairs of complex conjugates of γ. Note that r + 2r 2 = n i. We will apply the following lemma. Lemma 4.3 [5, Lemma 3.3]. Let T k k r +r 2 be positive integers and set T r+r 2+k = T r+k for k r 2. Then } # a Z[X]/p i : a k Tk = c i T T ni +O T ni 0, where T 0 = max,t T ni /ni and ci is a constant depending on Z[X]/p i. and Furthermore since 4.5 holds, we get that there exists a Z linear mapping M i such that m i #R i T = M i T i,...,t i,mi,n i = T i,..., T i,mi,n i j= # a Z[X]/p i :, a k Tijk k n i }. As the value of T ijk is independent from j, an application of Lemma 4.3 yields ni #R i T = c i k= T ik mi where c i depends on c i and M i and T i0 = max,t i T i,mi,n i For the estimate involving N i T we note that +O T mini i0, m i n i N i = g R i : g 0 mod p i } = Z[X]/p i mi and the result follows in the same way as for R i T. With help of all these tools we can show Proposition 4.. Proof of Proposition 4.. The first step consists in splitting the sum over RT up into those over R i T. Therefore let π i : R R i be the canonical projections. Then π := π,...,π t is an isomorphism by the Chinese Remainder Theorem. Furthermore let φ i be the embedding defined by K φ i : i R nimi, a +a 2 X + +a min i X nimi a,...,a min i.

11 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS for i =,...,t. Finally we define the matrix M to be such that M φz := φ π z,...,φ t π t z. Furthermore we note that for P i := π i P and l Z we have Thus z RT M φpzx l = φ P z X l,...,φ t Pz t X l. h e hr,φpzx lr = r= t i= z i R it h e hri,φ i P i z i X lr where h r,...,h rt := h r M. Now we will consider each sum over R i T separately. Therefore we fix until the end of the proof an i t and distinguish two cases according to whether m i = or not. 4.6 Case, m i = : In this case we set β = β i and observe that K = K i = Qβ i and R i = Z[β]. Furthermore let OK be the maximum order aka the ring of integers of K, then clearly Z[β i ] O K. We denote by β ik = β k the conjugates of β. Now we will proceed as in the proof of Proposition 6. of Madritsch [5]. Therefore we need some parameters of the field K and for short we set n = n i during this case. Then we order the conjugates by denoting with β k for k r the real conjugates, whereas β k and β k+r2 denote the pairs of complex conjugates, where n = r +2r 2. Let Tr be the trace of an element of K over Q, then we define r= τz := Trz,Trβz,...,Trβ n z = Ξφ i z, where Ξ = VV T and V is the Vandermonde matrix β β 2 β n n V =.... β n β 2 n β n n We set h r,..., h rn := h ri Ξ and note that hri,φ i P i z i X lr n = h T ri Ξ τp i z i β lr = Tr h rk β k lr 2 P i z i. Thus we may rewrite the sum under consideration as follows h e hri,φ i P i z i X lr = h n e Tr h rk β k lr 2 P i z i. z R it r= z R it k= r= k= Now we need an approximation lemma which essentially goes back to Siegel [9]. Therefore let δ be the different of K over Q and be the absolute value of the discriminant of K. Then we have the following. Lemma4.4. Let N,...,N r+r 2 be real numbers andlet N = n N N r N r+ N r+r 2 2 be their geometric mean. Suppose that N > n, then, corresponding to any ξ K, there exist q O K and a δ such that q k ξ k a k < N q k, 0 < k N k, k r +r 2, q max N k k ξ k a k, q k 2, k r +r 2, and Nq,aδ 2.

12 2 M. G. MADRITSCH AND A. PETHŐ 4.7 For r h we set ξ r to be the leading coefficient of n k= h rk P i z. Then we apply Lemma 4.4 with N k = Ti,,k d for k r logt σ2 +r 2 in order to get that there exist a δ and q O K such that h ξ k r β k lr+ qk a k < logtσ2 and 0 < q k Ti,,k d < for k n. logt σ2 r= Ti,,k d Lemma 4.5 [5, Proposition 3.2]. Suppose that QX = α d X d + +α X is a polynomial of degree d with coefficients in K. If for the leading coefficient α d there exist a δ and q O K as in Lemma 4.4 with N k = Tik d and logt σ2 logt σ2 q k T d ik logt σ2 k r +r 2, then x R it with σ 2 2 d σ 0 +r2 2d. etrqx T ni logt σ0 Now we distinguish several cases according to the quality of approximation by Lemma 4.4, which is represented by the size of Hq: Case., Hq logt σ2 : We apply Lemma 4.5 and get h n h rk P i z i E T n logt σ0. z i R it r= k= β lr+ Case.2, 2 Hq < logt σ2 : In the last two cases we need Minkowski s lattice theory cf. [9]. Let λ be the first successive minimum of the Z-lattice δ. Then we get h ξ k r β k lr+ r= a k q k q k 2 λ q k q k 2 which implies h ξ k logt σ2 r β k lr+ r= β k l+. h r=ξ r k Since β k > and 4.2 we have β k l + ξ k logt σ2 nlogt σ2+σ, which yields l σ 2 +σ log β k loglogt λ 2 q k logt σ2, contradicting the lower bound of l for sufficiently large C in 4.. Case.3, 0 < Hq < 2: In this case we will again use Minkowski s lattice theory cf. [9]. Let λ be the first successive minimum of the Z-lattice δ, then we have to consider two subcases Case.3., H then h r= ξ r β lr+ q λ 2 h ξ r k β k lr+qk r= : Let k n be such that ξ r k q k β k l+, λ 2 h r= l + loglogt

13 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS 3 again contradicts the lower bound of l for sufficiently large C in 4.. h ξ Case.3.2, H r r= q < λ β lr+ 2 : By Minkowski s theorem cf. [9] we get that a = 0. Thus for k n h ξ r k β k lr+qk = h ξ k β k lr+ r β k l h l r q k logtσ2 Tik d r= r= which implies taking the norm of the left side l h + ndlog Nβ T cloglog Nβ T contradicting the upper bound for sufficiently large C 2. Case 2, m i > : Now we have to go one step further and to take a closer look at R i. In particular we divide every element z i R i into its radical and its nilpotent part. We fix an element z R and set z i := π i z. On the one hand, since R i = R i N i we have for z i R i the unique representation 4.8 z i = z i +z i2 z i R it = with z i R i and z i2 N i. This motivates the definition of the linear map π ij such that π ij z := z ij for i =,...,t and j =,2. On the other hand, since R i = Z[X]/pi mi we have for every z i R i the unique representation m i z i = a ij p j i = j= m i n i j= k= a ijk X k p j i with a ij Z[X] and a ijk Z, respectively. We clearly have n i m i n i π i z = a ik X k and π i2 z = a ijk X k p j i. k= j=2 k= Thus we define for z i R i the embeddings ψ i and ψ i2 by ψ i π i z i = a i,...,a ini and ψ i2 π i2 z i = a i2,...,a i,mi,n i. Then there exists an invertible matrix M i such that M i φ i π i z = ψ i π i z,ψ i2 π i2 z. Now we can divide the sum up as follows. h e hri,φ i P i z i X lr r= z i R it z i2 N it = z i R it h e hri,φ i P i z i +z i2 X lr r= h e hri,ψ i P i z i X lr r= z i2 N it h e hri2,ψ i2 P i2 z i2 X lr, where we have set P ij = π ij P for j =,2. Since for the first sum we have that m i = we may follow Case above and use Lemma 4.3 for trivially estimating the second one to prove the proposition for this case. r=

14 4 M. G. MADRITSCH AND A. PETHŐ 5. Treatment of the border In Section 3 above we have constructed the Urysohn function we need in order to properly count the number of elements within the fundamental domain. In this construction we also used an axe-parallel tube in order to cover the border of the fundamental domain. The number of hits of this tube gives rise to the error term which we will consider in this section. We fix a positive integer v, which will be chosen later, and a real vector T. Furthermore for l 0 we define F l to be the number of hits of the border of the Urysohn function which is F l := # z RT B l φpz } 5. P v,a mod B Z n. As indicated above we are interested in an estimation of F l. Proposition 5.. Let µ < detb be as in Section 3 and C and C 2 be sufficiently large positive reals. Suppose that l is a positive integer such that 5.2 Then for any positive σ 3 we have C loglogt l dlog detb T C 2 loglogt. F l µ v T n detb v +logt tσ3. In order to estimate F l we need the Erdős-Turán-Koksma Inequality. Lemma 5.2 [4, Theorem.2]. Let x,...,x S be points in the n-dimensional real vector space R n and H an arbitrary positive integer. Then the discrepancy D S x,...,x S fulfills the inequality D S x,...,x S 2 H + + S eh,x s rh S, 0< h H where h Z n and rh = n i= max, h i. Proof of Proposition 5.. We want to proceed in three steps. First we subdivide the tube P v,a into rectangles in order to apply the Erdős-Turán-Koksma Inequality in the second step. Finally we put them together to gain the desired result. Recall that the tube P v,a defined in Lemma 3. consists of a family of rectangles. Let R a be one of them, then we want to estimate F l R a := # z RT B l φpz } R a mod B Z n. Using the definition of the discrepancy we get that F l R a T n B λr a +D l S φpz } 5.3, z RT where λ is the n-dimensional Lebesgue measure and S is the number of elements in RT. By Lemma 4.2 we have that t ni mi S = c i T ik +O T0 n Now we apply Lemma 5.2 to get D S B l φpz } z RT i= k= 2 H + + 0< h H rh S s= z RT The next step consists in an application of Proposition 4. which yields e h,b l φpz 5.6 Tn logt tσ0 z RT e h,b l φpz.

15 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS 5 Putting 5.4, 5.5, and 5.6 togetter in 5.3 gives F l R a T n λr a + T n λr a + T n logt σ +Tn logt tσ0 0< h H T n logt σ +Tn logt tσ0 loglogt n. Setting σ := tσ 0 /2 and summing over all rectangles R a yields F l µ v T n Nb v +logt tσ0/2. rh Finally we set σ 3 = tσ 0 /2 which proves the proposition. 6. The main proposition The main idea is to understand the additive function as putting weights on the digits. Thus if we can show that the digits areuniformly distributed the same is true for the values of the additive functions. Therefore we look at patterns in the expansion of Pz. In particular, we count the number of occurrences of certain digits at certain positions in the expansions. Proposition 6.. Let f be an additive function. Let T 0 and T ijk as in 2.5. Let L be the maximum length of the b-adic expansion of z RT and let C and C 2 be sufficiently large. Then for 6. we have, C logl l < l 2 < < l h dl C 2 logl Θ := #z RT a lr fz = b r,r =,...,h} = c c t detb htn +O T n logt tσ0. uniformly for T, where l r,b r N N are given pairs of position and digit and σ 0 is an arbitrary positive constant. Proof. We recall our Urysohn function u a defined in 3.5 and set for ν R n tν = u b B l ν u bh B l h ν, where B is the matrix defined in 3.3. Now we want to apply the Fourier transformation, which we developed in Lemma 3.2. Therefore we set M := M = h,...,h h h r Z n, for r =,...,h}. An application of Lemma 3.2 yields 6.2 tν = M M h T M e h r B lr ν, where T M = h r= c m r,...,m rn. Combining this with the definition of F l in 5. we get Θ 6.3 tφpz F l + +F lh. z RT r= Plugging 6.2 into 6.3 together with an application of Lemma 3.2 for the coefficients yields Θ = c c t detb htn + h h T M e h r,b lr φpz +O F lr. 0 M M r= r=

16 6 M. G. MADRITSCH AND A. PETHŐ Now an application of Proposition 4. to treat the exponential sums, of Proposition 5. for the border F l with v loglogt and the observation that T M κ 2h detb 2hv logt tσ0/2, M M where we used the definition of κ in 3.6, proves the proposition. 7. Proof of Theorem 2.2 For this proof we mainly follow the proof of the Theorem of Bassily and Kátai [3]. In the same manner we cut of the head and tail of the expansion and show the theorem for a truncated version of the additive function. In particular we set C := maxc,c 2, A := [ClogL] and B := L A, where L, C and C 2 are defined in the statement of Proposition 6.. Furthermore we define the truncated function f to be B f Pz = fa j Pzb j. j=a By the definition of A and fab j with a N we get that f Pz = fpz+ologl. In the same manner we define the truncated mean and standard deviation M T := B m j and D 2 T := j=a B σj 2. At this point we need that the deviation D tends sufficiently fast do infinity. In particular, we could refine the statement, if we shrink the part, which is cut of. Since MT M T = OlogL and D 2 T D 2 T = OlogL we get that it suffices to show that #RT # z RT f Pz M T d D T d j=a } < y Φy. By the Fréchet-ShohatTheoremcf. [6, Lemma.43] this holds true ifand only if the moments ξ k T := #RT z RT f Pz M T d D T d converge to the moments of the normal law for T. We will show the last statement by comparing the moments ξ k with η k T := #RT d z NT d where T d = T d,...,td n = Td,...,Td t,n t,m t. An application of Proposition 6. gives that k f z M T d k D T d, ξ k T η k T 0 for T. Furthermore we get by Proposition 2. that these sums consist of independently identically distributed random variables with possible 2C exceptions. By the central limit theorem we get that their distribution converges to the normal law. Thus the η k T converge to the moments of the normal law. This yields lim ξ kt = lim η kt = x k dφ. T T We apply the Fréchet-Shohat Theorem again to prove the theorem.

17 ASYMPTOTIC NORMALITY IN CANONICAL NUMBER SYSTEMS 7 Acknowledgment This paper was written while M. Madritsch was a visitor at the Faculty of Informatics of the University of Debrecen. He thanks the centre for its hospitality. During his stay he was supported by the project HU 04/200 founded by the ÖAD. The second author was supported by the Hungarian National Foundation for Scientific Research Grant No.T67580 and by the TÁMOP 4.2./B-09//KONV project. The second project is implemented through the New Hungary Development Plan co-financed by the European Social Fund, and the European Regional Development Fund. References [] S. Akiyama, T. Borbély, H. Brunotte, A. Pethő, and J. M. Thuswaldner, Generalized radix representations and dynamical systems. I, Acta Math. Hungar , no. 3, [2] S. Akiyama and H. Rao, New criteria for canonical number systems, Acta Arith. 2004, no., [3] N. L. Bassily and I. Kátai, Distribution of the values of q-additive functions on polynomial sequences, Acta Math. Hungar , no. 4, [4] M. Drmota and R. F. Tichy, Sequences, discrepancies and applications, Lecture Notes in Mathematics, vol. 65, Springer-Verlag, Berlin, 997. [5] J. M. Dumont, P. J. Grabner, and A. Thomas, Distribution of the digits in the expansions of rational integers in algebraic bases, Acta Sci. Math. Szeged , no. 3-4, [6] P. D. T. A. Elliott, Probabilistic number theory. I, Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Science], vol. 239, Springer-Verlag, New York, 979, Mean-value theorems. [7] B. Gittenberger and J. M. Thuswaldner, Asymptotic normality of b-additive functions on polynomial sequences in the Gaussian number field, J. Number Theory , no. 2, [8] K. Gröchenig and A. Haas, Self-similar lattice tilings, J. Fourier Anal. Appl. 994, no. 2, [9] L.-K. Hua, Introduction to number theory, Springer-Verlag, Berlin, 982, Translated from the Chinese by Peter Shiu. [0] I. Kátai and I. Környei, On number systems in algebraic number fields, Publ. Math. Debrecen 4 992, no. 3-4, [] B. Kovács and A. Pethő, Number systems in integral domains, especially in orders of algebraic number fields, Acta Sci. Math. Szeged 55 99, no. 3-4, [2], On a representation of algebraic integers, Studia Sci. Math. Hungar , no. -2, [3] M. G. Madritsch, A note on normal numbers in matrix number systems, Math. Pannon , no. 2, [4], Generating normal numbers over Gaussian integers, Acta Arith , no., [5], Asymptotic normality of b-additive functions on polynomial sequences in number systems, Ramanujan J , no. 2, [6] W. Müller, J. M. Thuswaldner, and R. F. Tichy, Fractal properties of number systems, Period. Math. Hungar , no. -2, [7] A. Pethő, On a polynomial transformation and its application to the construction of a public key cryptosystem, Computational number theory Debrecen, 989, de Gruyter, Berlin, 99, pp [8] K. Scheicher and J. M. Thuswaldner, Canonical number systems, counting automata and fractals, Math. Proc. Cambridge Philos. Soc , no., [9] C. L. Siegel, Generalization of Waring s problem to algebraic number fields, Amer. J. Math , M. G. Madritsch Department of Analysis and Computational Number Theory Graz University of Technology A-800 Graz, Austria address: madritsch@math.tugraz.at A. Pethő Department of Computer Science, University of Debrecen, Number Theory Research Group, Hungarian Academy of Sciences and University of Debrecen P.O. Box 2, H-400 Debrecen, Hungary address: petho.attila@inf.unideb.hu