The Rudin-Shapiro sequence and similar sequences are normal along squares
Clemens M\"ullner

TL;DR
This paper proves that certain digital sequences, including the Rudin-Shapiro sequence, are normal along squares, providing explicit and efficiently generatable normal numbers for any base.
Contribution
It establishes the normality of digital sequences along squares, encompassing sequences like sum of digits and Rudin-Shapiro, which was previously unknown.
Findings
Sequences like sum of digits modulo m are normal along squares.
Rudin-Shapiro sequence is normal along squares.
Normal numbers can be explicitly constructed and efficiently generated.
Abstract
We prove that digital sequences modulo along squares are normal, which covers some prominent sequences like the sum of digits in base modulo , the Rudin-Shapiro sequence and some generalizations. This gives, for any base, a class of explicit normal numbers that can be efficiently generated.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
The Rudin-Shapiro sequence and similar sequences are normal along squares
Clemens Müllner
Institut für Diskrete Mathematik und Geometrie TU Wien
Wiedner Hauptstr. 8-10, 1040 Wien, Austria
Abstract.
We prove that digital sequences modulo along squares are normal, which covers some prominent sequences like the sum of digits in base modulo , the Rudin-Shapiro sequence and some generalizations. This gives, for any base, a class of explicit normal numbers that can be efficiently generated.
2010 Mathematics Subject Classification:
Primary: 11A63, 11L03, 11B85; Secondary: 11N60, 60F05
This research is supported by the project F55-02 of the Austrian Science Fund FWF which is part of the Special Research Program “Quasi-Monte Carlo Methods: Theory and Applications”, by Project F5002-N15 (FWF), which is a part of the Special Research Program “Algorithmic and Enumerative Combinatorics” and by Project I1751 (FWF), called MUDERA (Multiplicativity, Determinism, and Randomness).
1. Introduction
This paper deals with digital sequences modulo . Such sequences are “simple” in the sense that they are deterministic and uniformly recurrent sequences. We show that the situation changes completely when we consider the subsequence along squares, i.e., we show that this subsequence is normal.
Thus, we describe a new class of normal numbers that can be efficiently generated, i.e., the first digits of the normal number can be generated by using elementary operations.
In this paper we let denote the set of positive integers and we let denote the set of prime numbers. We let denote the set of complex numbers of modulus and we use the abbreviation exp for any real number .
For two functions, and that take only strictly positive real values we write or if is bounded.
We let denote the floor function and denote the fractional part of . Furthermore, we let denote the indicator function for in .
Moreover we let denote the number of divisors of , denote the number of distinct prime factors of and denote the number of positive integers smaller than that are co-prime to .
Furthermore, let denote the -th digit in the base expansion of a non-negative integer , i.e., , where . We usually omit the superscript, as we work with arbitrary but fixed base .
1.1. Digital Sequences
The main topic of this paper are digital sequences modulo . We use a slightly different definition of digital functions than the one found in [1].
Definition 1.1**.**
We call a function a strongly block-additive -ary function or digital function if there exist and such that and
[TABLE]
where we define for all . The difference to the usual definition is the range of the sum ( or ) which does not matter for all appearing examples.
Remark 1.2**.**
The name strongly block-additive -ary function was inspired by (strongly) -additive functions. Bellman and Shapiro [3] and Gelfond [9] denoted a function to be -additive if
[TABLE]
holds for all , and . Mendès France [14] denoted a function to be strongly -additive if
[TABLE]
holds for all , and . Thus, we can write for a strongly -additive function ,
[TABLE]
A quite prominent example of a strongly block-additive function is the sum of digits function in base . This is a strongly block-additive function with and . In particular, gives the well-known Thue–Morse sequence.
Another prominent example is the Rudin-Shapiro sequence which is given by the parity of the blocks of the form “” in the digital expansion in base . Let be the digital sequence corresponding to and , then we find . This can be generalized to functions that are given by the parity of blocks of the form “” for fixed length of the block; these functions have for example been mentioned and studied in [13].
Digital sequences are regular sequences (see for example [5]). Consequently we find that digital sequences modulo are automatic sequences (see [1, Corollary 16.1.6]) which implies some interesting properties. For a detailed treatment of automatic sequences see [1].
We define the subword complexity of a sequence , that takes only finitely many different values, as
[TABLE]
It is well known that the subword complexity of automatic sequences is sub-linear (see [1, Corollary 10.3.2]), i.e. for every automatic sequence we have
[TABLE]
For a random sequence one finds that with probability one. Thus, automatic sequences are far from being random.
1.2. Main Result
It is well known that these properties are preserved when considering arithmetic subsequences of automatic sequences and, therefore, digital sequences modulo . However, the situation changes completely when one considers the subsequence along squares.
Definition 1.3**.**
A sequence is normal if, for any and any , we have
[TABLE]
Drmota, Mauduit and Rivat showed a first example for that phenomenon [6]. They considered the classical Thue–Morse sequence and showed that not only , but were able to show that is normal. The fact that had already been proven by Moshe [15], who was able to give exponentially growing lower bounds for extractions of the Thue–Morse sequence along polynomials of degree at least . In this paper we go one step further than Drmota, Mauduit and Rivat and show a similar result for general digital sequences.
Theorem 1.4**.**
*Let be a digital function and with and
. Then is normal.*
There are only few known explicit constructions of normal numbers in a given base (see for example [4, Chapters 4 and 5]). This result provides us with a whole class of normal sequences for any given base that can be generated efficiently, i.e. it takes elementary operations to produce the first elements.
The easiest construction for normal sequences is the Champernowne construction that is given by concatenating the base expansion of successive integers. This gives for example for base : . Using the first integers takes elementary operations and gives a sequence of length .
Scheerer [17] analyzed the runtime of some algorithms that produce absolutely normal numbers, i.e. real numbers in whose expansion in base are normal for every base . Algorithms by Sierpinski [19] and Turing [20] use double exponentially many operations and algorithms by Levin [11] and Schmidt [18] use exponentially many operations. Moreover, Becher, Heiber and Slaman [2] gave an algorithm that takes just above operations to produce the first digits.
Digital sequences modulo have interesting (dynamical) properties. Firstly, they are primitive and, therefore, uniformly recurrent ([1, Theorem 10.9.5]) , i.e., every block that occurs in the sequence at least once, occurs infinitely often with bounded gaps.
There is a natural way to associate a dynamical system - the symbolic dynamical system - to a sequence that takes only finitely many values.
Definition 1.5**.**
The symbolic dynamical system associated to a sequence is the system , where is the shift on and the closure of the orbit of under the action of for the product topology of .
Some of the mentioned properties of automatic sequences also imply important properties for the associated symbolic dynamical system.
The fact that every digital sequence modulo , denoted by , is uniformly recurrent implies that the associated symbolic dynamical system is minimal; i.e., the only closed invariant sets in are and - see for example [8] or [16].
Furthermore, the entropy of symbolic dynamical systems to a sequence , that takes only finitely many values, is equal to
[TABLE]
(see for example [10] or [7]). Consequently, we know that the entropy of the symbolic dynamical system associated to a digital sequence modulo equals [math], and, therefore, the dynamical system is deterministic.
1.3. Outline of the proof
In order to prove our main result, we will work with exponential sums. We present here the main theorem on exponential sums and further show its connection to Theorem 1.4.
Theorem 1.6**.**
For any integer and such that , there exists such that
[TABLE]
Lemma 1.7**.**
Theorem 1.6 implies Theorem 1.4.
Proof.
Let be an arbitrary sequence of length . We count the number of occurrences of this sequence in . Assuming that (1.1) holds, we obtain by using the well known identity for and [math] otherwise
[TABLE]
[TABLE]
with the same as in Theorem 1.6.
To obtain the last equality we separate the term with . ∎
The structure of the rest of the paper is presented below.
In Section 2 we discuss some properties of digital sequences. These properties will be very important for the estimates of the Fourier terms.
In Section 3, we derive the main ingredients of the proof of Theorem 1.6 which are upper bounds on the Fourier terms
[TABLE]
where with some special properties defined in Section 3.2 and is a truncated version of which is properly defined in Definition 2.1.
The main results of Section 3 are Propositions 3.7 and 3.8. Proposition 3.7 yields a bound on averages of Fourier transforms and Proposition 3.8 yields a uniform bound on Fourier transforms.
In Section 4, we discuss how Proposition 3.7 and Proposition 3.8 are used to prove Theorem 1.6. The approach is very similar to [6] and we will mainly describe how it has to be adapted. We use Van-der-Corput-like inequalities in order to reduce our problem to sums depending only on few digits of . By detecting these few digits, we are able to remove the quadratic terms, which allows a proper Fourier analytic treatment. After the Fourier analysis, the remaining sum is split into two sums. The first sum involves quadratic exponential sums which are dealt with using the results from Section 5.2.
The Fourier terms appear in the second sum and Propositions 3.7 and 3.8 provide the necessary bounds.
We have to distinguish the cases and . Sections 4.1 and 4.2 tackle one of these cases each. In Section 4.1, we prove that – if – we deduce Theorem 1.6 from Proposition 3.7. For , Section 4.2 shows that we can deduce Theorem 1.6 from Proposition 3.8.
In Section 5, we present some auxiliary results also used in [6].
2. Digital Functions
In this section we discuss some important properties of digital functions. We start with some basic definitions.
Definition 2.1**.**
We define for the truncated function and the two-fold restricted function by
[TABLE]
We see directly that is a periodic function and we extend it to a ( periodic) function which we also denote by .
For any , we define .
Since , we can rewrite and for as follows
[TABLE]
We show that for any block-additive function, we can choose without loss of generality such that it fulfills a nice property.
Lemma 2.2**.**
Let be a strongly block-additive function corresponding to . Then, there exists another function such that also corresponds to and
[TABLE]
holds for all .
Proof.
We start by defining a new function
[TABLE]
This already allows us to define the function :
[TABLE]
We find directly that . It remains to show that corresponds to and that (2.1) holds, which are simple computations,
[TABLE]
Furthermore, we find
[TABLE]
∎
We assume from now on that for any strongly block-additive function (2.1) holds. This allows us to find an easier expression for :
Corollary 2.3**.**
Let be a digital function fulfilling (2.1). Then
[TABLE]
and
[TABLE]
holds for all .
We easily find the following recursion.
Lemma 2.4**.**
Let and . Then
[TABLE]
holds for all and
[TABLE]
Proof.
We compute
[TABLE]
The second case can be treated analogously. ∎
As we are dealing with the distribution of digital functions along a special subsequence, we will start discussing some distributional result for digital functions.
Lemma 2.5**.**
Let be a strongly block-additive function and . Then the following three statements are equivalent.
** 2.
** 3.
**
Proof.
Obviously . Next we show :
Let be the smallest natural number such that . By Lemma 2.4 holds
[TABLE]
By the definition of , holds and, therefore, .
It remains to prove :
Let be the smallest natural number such that . By , we have . We compute ,
[TABLE]
as for implies that . ∎
Remark 2.6**.**
The following example shows that we can not replace by in Lemma 2.5:
Let and . We see that for all and also for all . However, and .
Next, we show a technical result concerning block-additive functions, which will be useful later on.
Lemma 2.7**.**
Let be a strongly block-additive function in base and such that and . Then there exist integers such that
[TABLE]
holds.
Proof.
Without loss of generality we can restrict ourselves to the case where . Let us assume on the contrary that there exists such that
[TABLE]
holds for all . Under this assumption, we find a new expression for , where :
[TABLE]
The last equality holds since is a periodic function in . This gives
[TABLE]
By comparing this expression for and - note that - we find
[TABLE]
as .
Together with (2.5), this implies that for all . This is a contradiction to by Lemma 2.5. ∎
We will use this result in a different form.
Corollary 2.8**.**
Let be a strongly block-additive function in base and such that and . For every exist and such that and
[TABLE]
Proof.
Let where and . We apply Lemma 2.7 for and find such that
[TABLE]
where
[TABLE]
This implies
[TABLE]
∎
3. Bounds on Fourier Transforms
The goal of this section is to prove Propositions 3.7 and 3.8. To find the necessary bounds we first need to recall one important result on the norm of matrix products which was first presented by Drmota, Mauduit and Rivat [6].
Afterwards, we deal with Fourier estimates and formulate Proposition 3.7 and Proposition 3.8. The following Sections 3.3 and 3.4 give proofs of Proposition 3.7 and Proposition 3.8, respectively.
3.1. Auxiliary Results for the Bounds of the Fourier Transforms
In this section we state necessary conditions under which the product of matrices decreases exponentially with respect to the matrix row-sum norm.
Lemma 3.1**.**
Let , , be -matrices with complex entries , for , and absolute row sums
[TABLE]
Furthermore, we assume that there exist integers and and constants and such that
- (1)
every product of consecutive matrices has the property that,
[TABLE] 2. (2)
every product of consecutive matrices has the property
[TABLE]
Then there exist constants and such that
[TABLE]
uniformly for all and (where denotes the matrix row-sum norm).
Proof.
See [6]. ∎
Lemma 3.2**.**
Let . Then
[TABLE]
Proof.
The proof is a straight-forward computation and can be found for example at the end of the proof of [13, Lemma 12]. ∎
3.2. Fourier estimates
In this section, we discuss some general properties of the occurring Fourier terms.
For any , we denote by the set of integer vectors with and for .
Furthermore, we denote by the set of integer vectors with and .
This set obviously consists of elements. For any , and , we define
[TABLE]
for fixed coefficients . This sum can then be seen as the discrete Fourier transform of the function
[TABLE]
which is periodic.
Furthermore, we define the important parameter
[TABLE]
We would like to find a simple recursion for in terms of . Instead we relate it to a different function for which the recursion is much simpler:
[TABLE]
This sum can then be seen as the discrete Fourier transform of the function
[TABLE]
which is periodic. We show now how and are related.
Lemma 3.3**.**
*Let and . It holds *
[TABLE]
where
[TABLE]
Proof.
One checks easily that . We evaluate :
[TABLE]
∎
Next we define a transformation on and a weight function .
Definition 3.4**.**
Let and . Then, we define for
[TABLE]
We see immediately that for all possible values of and . Furthermore, we extend the definition of for arbitrary by
[TABLE]
The next Lemma shows some basic properties of these functions.
Lemma 3.5**.**
Let , and . Then, the following facts hold.
- •
- •
- •
.
Proof.
The first two facts are direct consequences of basic properties of the floor function and the last fact is just a reformulation of the definition of in terms of . ∎
Now we can find a nice recursion for the Fourier transform .
Lemma 3.6**.**
Let and . We have
[TABLE]
Proof.
We evaluate and use (2.2):
[TABLE]
∎
The following propositions are crucial for our proof of the main Theorem 1.6.
Proposition 3.7**.**
If and , then there exists such that for any
[TABLE]
holds uniformly for all integers .
Proposition 3.8**.**
If , then there exists such that for any
[TABLE]
holds uniformly for all non-negative integers and .
Proofs for Proposition 3.7 and 3.8 are given in the following sections.
3.3. Proof of Proposition 3.7
This section is dedicated to prove Proposition 3.7. We start by reducing the problem from to for which we have found a nice recursion.
Proposition 3.9**.**
For and , we find such that for any
[TABLE]
holds uniformly for all integers .
Lemma 3.10**.**
Proposition 3.9 implies Proposition 3.7.
Proof.
We see by (3.5) that
[TABLE]
Thus we find
[TABLE]
∎
Using Lemma 3.6, it is easy to establish a recursion for
[TABLE]
where , and . For and it yields for the following expression
[TABLE]
To find this recursion, one has to split up the sum over into the equivalence classes modulo .
This identity gives rise to a vector recursion for . We use the recursion for :
[TABLE]
where the -matrix is independent of and . By construction, all absolute row sums of are bounded by .
It is useful to interpret these matrices as weighted directed graphs. The vertices are the pairs and, starting from each vertex, there are directed edges to the vertices
- where - with corresponding weights
[TABLE]
Products of such matrices correspond to oriented paths of length in these graphs, which are weighted with the corresponding products. The entries at position of such product matrices correspond to the sum of weights along paths from to . Lemma 3.6 allows us to describe this product of matrices directly.
Lemma 3.11**.**
The entry of equals
[TABLE]
Proof.
Follows directly by Lemma 3.6. ∎
This product of matrices corresponds to oriented paths of length . They can be encoded by the triple and they correspond to a path from to with unimodular weight .
To simplify further computations we define
[TABLE]
and find directly that
[TABLE]
and the absolute value of the entry of is bounded by .
In order to prove Proposition 3.7, we will use Lemma 3.1 uniformly for with . Therefore, we need to check Conditions (3.1) and (3.2).
Note that, since , we have
[TABLE]
Lemma 3.12**.**
The matrices defined above fulfill Condition (3.1) of Lemma 3.1.
Proof.
We need to show that there exists an integer such that every product
[TABLE]
of consecutive matrices verifies Condition (3.1) of Lemma 3.1.
We define . It follows directly from the definition, that for all . In the graph interpretation this means that for every vertex there is a path of length from to . Fix a row indexed by in the matrix . We already showed that the entry is the sum of at least one term of absolute value , i.e., .
There are two possible cases. If the absolute row sum is at most
[TABLE]
with then we are done.
In case the absolute row sum is strictly greater than , we show that : The inequality implies that is the sum of at least two terms of absolute value , i.e. . Thus, we can use the triangle inequality to bound the absolute row sum by
[TABLE]
Since
[TABLE]
we find
[TABLE]
This contradicts the assumption that the absolute row sum is strictly greater than
[TABLE]
Consequently, we find
[TABLE]
∎
Lemma 3.13**.**
The matrices fulfill Condition (3.2) of Lemma 3.1.
Proof.
We need to show that there exists an integer such that for every product
[TABLE]
of consecutive matrices the absolute row-sum of the first row is bounded by . We concentrate on the entry , i.e. we consider all possible paths from to of length in the corresponding graph and show that a positive saving for the absolute row sum is just due to the structure of this entry.
Since , we have at least two paths from to and it follows that the entry is certainly a sum of terms of absolute value (for every ). This means that there are paths from to of length in the corresponding graph, or in other words
Our goal is to construct two paths from to such that
[TABLE]
holds for all .
We construct a path from to with exactly times (where ). We set and find the following lemma.
Lemma 3.14**.**
Let and be as above. Then
[TABLE]
Proof.
This follows directly by the definitions and simple computations. ∎
By applying Lemma 3.14 we find a transformation from to . This gives a path from to by applying this transformation component-wise. We concatenate this path with another path of length where . The weight of the concatenation of these two paths equals
[TABLE]
We denote by the -th coordinate of and see that
[TABLE]
Thus, we have found for each a path from to .
We can use the special structure of to make the weight of this path more explicit: At first, we note that
[TABLE]
by the definition of . Furthermore, we use the condition to find
[TABLE]
We find by the definition of that for each ,
[TABLE]
We find by Corollary 2.8 that there exist such that
[TABLE]
and .
We now compare the following two paths from to of length :
- •
: We split up this path into the path of length from to and the path of length from to : The first path can be described by the triple and its weight is obviously .
The second path - i.e. the path from to - can be described by the triple and its weight equals
[TABLE]
Thus, the overall weight of the path from to has weight
[TABLE]
- •
: we compute directly the weight of this path:
[TABLE]
We recall quickly that for all and, therefore, also . We finally see that
[TABLE]
Thus we have
[TABLE]
Therefore condition (3.2) of Lemma 3.1 is verified with and . ∎
At the end of this section, we want to recall the important steps of the proof of Proposition 3.7. At first we observe that
[TABLE]
Thus Proposition 3.7 is equivalent to . Next we considered the vector and find the recursion
[TABLE]
Then we defined and showed that we can apply Lemma 3.1. Therefore we know that – since
[TABLE]
with and obtained by Lemma 3.1. Thus we know that with uniformly for all . This concludes the proof of Proposition 3.7.
3.4. Proof of Proposition 3.8
We start again by reducing the problem from to , for possibly different values of and .
Proposition 3.15**.**
For there exists such that for any
[TABLE]
holds uniformly for all non-negative integers and .
Lemma 3.16**.**
Proposition 3.15 implies Proposition 3.8.
Proof.
Follows directly by (3.5). ∎
We assume from now on that holds.
We formulate Lemma 3.6 as a matrix vector multiplication:
[TABLE]
where for any and we have
[TABLE]
To prove Proposition 3.15 we aim to show that
[TABLE]
Indeed, we find that this is already sufficient to show Proposition 3.15.
Lemma 3.17**.**
(3.7) implies Proposition 3.15.
Proof.
We first note that
[TABLE]
holds for all , and by definition.
Next we split the digital expansion of - read from left to right - into parts of length and possible one part of length . We denote the first parts by and the last part by , i.e.,
[TABLE]
Thus we find
[TABLE]
where . ∎
Throughout the rest of this section, we aim to prove (3.7).
Therefore, we try to find for each and a pair and such that for all holds
[TABLE]
Let us assume for now that (3.8) holds. Indeed we find
[TABLE]
However, we find for each some fulfilling (3.8). This gives
[TABLE]
Thus, we find in total
[TABLE]
It just remains to find fulfilling (3.8) and this turns out to be a rather tricky task.
We fix now some arbitrary and . We start by defining for and
[TABLE]
and show some basic properties of .
Lemma 3.18**.**
For every exists such that
[TABLE]
Proof.
One finds easily that
[TABLE]
which means that is a partition of for each . Thus, we find for every
[TABLE]
and the proof follows easily. ∎
Lemma 3.19**.**
Let and .
Then, there exists such that for each exists such that
[TABLE]
Remark 3.20**.**
This is equivalent to the statement that
[TABLE]
implies
[TABLE]
Proof.
We have already seen that is a partition of . Furthermore, we find for and that
[TABLE]
This implies that is a refinement of and we find
[TABLE]
It is well known that the maximal length of a chain in the set of partitions of is . This means that there exists such that . ∎
Furthermore, we define
[TABLE]
We can now choose , where is given by Lemma 3.19. We consider and provided by Lemma 3.18 and Lemma 3.19 and know that . Therefore we apply Corollary 2.8 and find such that
[TABLE]
and .
We are now able to define
[TABLE]
It just remains to check (3.8) which we split up into the following two lemmata.
Lemma 3.21**.**
Let be defined as above. Then
[TABLE]
holds.
Proof.
We need to show that
[TABLE]
holds for all and . We know that belongs to for some . Thus, we find for
[TABLE]
Therefore, (3.9) does hold, unless
[TABLE]
We find
[TABLE]
We first consider the case :
[TABLE]
For :
[TABLE]
However
[TABLE]
as . Thus, (3.9) holds. ∎
Lemma 3.22**.**
There exists only depending on such that for and defined as above holds
[TABLE]
for all .
Proof.
We start by computing the weights . For arbitrary , we find:
[TABLE]
where
[TABLE]
Note that only depends on .
We can describe this product by using the weights defined above.
[TABLE]
Furthermore, we can rewrite every for which as some where . This gives then
[TABLE]
Thus we find for that:
[TABLE]
For , we find
[TABLE]
as .
Consequently, we find
[TABLE]
where
[TABLE]
and
[TABLE]
Also, we find
[TABLE]
where
[TABLE]
This implies
[TABLE]
It remains to apply Lemma 3.2 to find that (3.10) holds with . ∎
At the end of this section, we recall the important steps of the proof of Proposition 3.15.
We started to rewrite our recursion for into a matrix vector multiplication
[TABLE]
We then split up this matrix into a product of many matrices , where . Thereafter, we showed that , where . This implies then Proposition 3.15.
To show that , we found two different such that
[TABLE]
holds for all .
4. Proof of the Main Theorem
In this section, we complete the proof of Theorem 1.6 following the ideas and structure of [6]. As the proof is very similar, we only outline it briefly and comment on the important changes.
The structure of the proof is similar for both cases: At first we want to substitute the function by . This can be done by applying Lemma 5.5 and Lemma 5.7 in the case . For the case we have to use Lemma 5.7 first.
Thereafter, we apply Lemma 5.6 to detect the digits between and . Next, we use characteristic functions to detect suitable values for . Lemma 5.9 allows us to replace the characteristic functions by exponential sums. We split the remaining exponential sum into a quadratic and a linear part and find that the quadratic part is negligibly small. For the remaining sum, we apply Proposition 3.7 or 3.8 – depending on whether . The case needs more effort to deal with.
4.1. The case
In this section, we show that, if , Proposition 3.7 provides an upper bound for the sum
[TABLE]
Let be the unique integer such that and we choose all appearing exponents - i.e. , etc. - as in [6].
By using Lemma 5.5, and the same arguments as in [6], we find
[TABLE]
where
[TABLE]
Now we use Lemma 5.7 - with and - to relate to a sum in terms of :
[TABLE]
where
[TABLE]
and
[TABLE]
where is an interval included in (which we do not specify).
Next we use Lemma 5.6 to detect the digits of and between and - with a negligible error term. Therefore, we have to take the digits between and into account, where will be chosen later.
We set the integers , , , , and to satisfy the conditions of Lemma 5.6 and detect them by characteristic functions. Thus, we find
[TABLE]
where
[TABLE]
where is defined by (5.10) and . Lemma 5.9 allows us to replace the characteristic functions by trigonometric polynomials. More precisely, using (5.16) with and for some suitable (which is a fraction of chosen later), we have
[TABLE]
where and are the error terms specified in (5.16) and
[TABLE]
where we use the last sum to detect the correct value of .
The error terms , , can easily be estimated with the help of Lemma 5.4, just as in [6].
By using the representations of and , we obtain
[TABLE]
We now distinguish the cases and . For , we can estimate the exponential sum by using Lemma 5.4 and the following estimate
[TABLE]
Thus, we find
[TABLE]
This gives then
[TABLE]
where denotes the part of with .
We set and (where ). Furthermore, we define . As is contained in , we have - by the same arguments as in [6] -
[TABLE]
Using the estimate and the Cauchy-Schwarz inequality, yields
[TABLE]
We now replace by , by and apply Proposition 3.7.
[TABLE]
Next we average over and , as in [6], by applying Lemma 5.2. Thus we have a factor compared to . Combining all the estimates as in [6] gives then
[TABLE]
– provided that the following conditions hold
[TABLE]
For example, the choice
[TABLE]
ensures that the above conditions are satisfied.
Summing up we proved that for - where is given by Proposition 3.7 - holds
[TABLE]
which is precisely the statement of Theorem 1.6.
4.2. The case
In this section, we show that, for , Proposition 3.8 provides an upper bound for the sum
[TABLE]
Let , , and be integers satisfying
[TABLE]
to be chosen later - just as in [6]. Since we can not use Lemma 5.5 directly. Therefore, we apply Lemma 5.7 with and . Summing trivially for yields
[TABLE]
where
[TABLE]
and is an interval included in . By Lemma 5.5 we conclude that for all but values of . Therefore, we see that
[TABLE]
with
[TABLE]
This leads to
[TABLE]
and, by using the Cauchy-Schwarz inequality to
[TABLE]
For we can use Lemma 5.7 again: Let to be chosen later such that . After applying Lemma 5.7 with and
[TABLE]
we observe that for any we have
[TABLE]
and thus
[TABLE]
with
[TABLE]
where is an interval included in .
We now make a Fourier analysis similar to the case - as in [6]. We set and . We apply Lemma 5.6 and detect the correct values of by characteristic functions. This gives
[TABLE]
Furthermore, we use Lemma 5.9 to replace the characteristic functions by trigonometric polynomials. Using (5.16) with , and , and integers , verifying
[TABLE]
we obtain
[TABLE]
for the error terms obtained by (5.16) and obtained by replacing the characteristic function by trigonometric polynomials. We now reformulate by expanding the trigonometric polynomials, detecting the correct value of and restructuring the sums:
[TABLE]
One can estimate the error terms just as in [6] and finds that they are bounded by either or . In conclusion we deduce that
[TABLE]
We now split the sum into two parts:
[TABLE]
where denotes the contribution of the terms for which while denotes the contribution of the terms for which . We can estimate just as in [6] and find
[TABLE]
and it remains to consider . Setting , and , (where ) we can replace the two-fold restricted block-additive function by a truncated block-additive function
[TABLE]
Using the periodicity of modulo , we replace the variable by such that . Furthermore we introduce a new variable such that
[TABLE]
We then follow the arguments of [6] and find
[TABLE]
with
[TABLE]
The next few steps are again very similar to the corresponding ones in [6] and we skip the details. We find
[TABLE]
where
[TABLE]
Here we introduce the integers and such that
[TABLE]
This leads to
[TABLE]
where , and denote the contribution of the terms , and respectively.
Estimate of
By (5.2) we have
[TABLE]
and, therefore,
[TABLE]
By Proposition 3.8 (replacing by and by ), we find some such that
[TABLE]
By Parseval’s equality and recalling that , it follows that
[TABLE]
We obtain
[TABLE]
uniformly in , , , , and .
The remaining proof is analogue to the corresponding proof in [6]. The only difference is again that by using Lemma 5.2 we obtain a factor instead of . This gives
[TABLE]
which concludes this part.
Estimate of and
By following the arguments of [6] and applying the same changes as in the estimate of we find
[TABLE]
and
[TABLE]
Combining the estimates for
It follows from (4.16), (4.17) and (4.18) that
[TABLE]
Choosing
[TABLE]
we obtain
[TABLE]
Since , we obtain using (4.12) and (4.11), that
[TABLE]
We recall by (4.8) that and by (4.7) that , and insert the estimation from above in (4.9):
[TABLE]
For and , we obtain
[TABLE]
for all . Therefore we have seen that Proposition 3.8 implies the case of Theorem 1.6.
5. Auxiliary Results
In this last section, we present some auxiliary results which are used in Section 4, to prove the main theorem. For this proof, it is crucial to approximate characteristic functions of the intervals where by trigonometric polynomials. This is done by using Vaaler’s method - see Section 5.5. As we deal with exponential sums we also use a generalization of Van-der-Corput’s inequality which we have already seen in Section 5.4. In Section 5.1, we acquire some results dealing with sums of geometric series which we use to bound linear exponential sums. Section 5.2 is dedicated to one classic result on Gauss sums and allows us to find appropriate bounds on the occurring quadratic exponential sums in Section 4. The last part of this section deals with carry propagation. We find a quantitative statement that carry propagation along several digits is rare, i.e. exponentially decreasing.
We would like to note that all these auxiliary results have already been presented in [6].
5.1. Sums of geometric series
We will often make use of the following upper bound for geometric series with ratio and , :
[TABLE]
which is obtained from the formula for finite geometric series.
The following results allow us to find useful estimates for special double and triple sums involving geometric series.
Lemma 5.1**.**
Let with , and . For any real number , we have
[TABLE]
Proof.
This is [6, Lemma 6]. ∎
Lemma 5.2**.**
Let and be integers and . For any real number , we have
[TABLE]
and, if , we have an even sharper bound
[TABLE]
where denotes the number of divisors of .
Proof.
See [6]. ∎
5.2. Gauss sums
In the proof of the main theorem, we will meet quadratic exponential sums. We first consider Gauss sums which are defined by:
[TABLE]
In this section, we recall one classic result on Gauss sums, namely Theorem 5.3.
Theorem 5.3**.**
For all with ,
[TABLE]
holds.
Proof.
This form was for example obained in [12, Proposition 2]. ∎
Consequently we obtain the following result for incomplete quadratic Gauss sums.
Lemma 5.4**.**
For all with and , we have
[TABLE]
Proof.
This is Lemma 9 of [6]. ∎
5.3. Carry Lemmas
As mentioned before, we want to find a quantitative statement on how rare carry propagation along several digits is.
Lemma 5.5**.**
Let such that . For any integer with , the number of integers for which there exists an integer with is . Hence, we find for any block-additive function , that the number of integers with
[TABLE]
is also .
Proof.
A proof for the Thue-Morse sequence can be found in [6] and it is easy to adapt it for this more general case. ∎
The next lemma helps to replace quadratic exponential sums depending only on few digits.
Lemma 5.6**.**
Let such that , and and set . For integers , and we set
[TABLE]
where the integers , , , , , and satisfy the above conditions. Then for any integer the number of integers for which one of the following conditions
[TABLE]
is satisfied is .
Proof.
A proof for the sum of digits function in base can be found in [6] and it is straight forward to adapt it to fit this more general case. ∎
5.4. Van-der-Corput’s inequality
The following lemma is a generalization of Van-der-Corput’s inequality.
Lemma 5.7** ([12]).**
For all complex numbers and all integers and , we have
[TABLE]
where denotes the real part of .
5.5. Vaaler’s method
The following theorem is a classical method to detect real numbers in an interval modulo by means of exponential sums developed by Vaaler [21]. For with , we denote by the characteristic function of the interval modulo :
[TABLE]
The following theorem is a consequence of the mentioned paper by Vaaler. The presented form was first published by Mauduit and Rivat [13].
Theorem 5.8**.**
For all with and all integer , there exist real-valued trigonometric polynomials and such that for all
[TABLE]
The trigonometric polynomials are defined by
[TABLE]
with coefficients and satisfying
[TABLE]
Using this method we can detect points in a -dimensional box (modulo ):
Lemma 5.9**.**
For and with ,…, , we have for all
[TABLE]
where and are the real valued trigonometric polynomials defined by (5.12).
Proof.
See again [13]. ∎
Let with ,…, and define ,…,. For and we have
[TABLE]
Let with , and such that . If , we can express the sum
[TABLE]
as
[TABLE]
We now define with ,…, ,
[TABLE]
Lemma 5.10**.**
With the notations from above, we have
[TABLE]
Proof.
See again [13]. ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] J.-P. Allouche and J. Shallit. Automatic sequences. Theory, applications, generalizations. Cambridge: Cambridge University Press, 2003.
- 2[2] V. Becher, P. A. Heiber, and T. A. Slaman. A polynomial-time algorithm for computing absolutely normal numbers. Inform. and Comput. , 232:1–9, 2013.
- 3[3] R. Bellman and H. N. Shapiro. On a problem in additive number theory. Ann. of Math. (2) , 49:333–340, 1948.
- 4[4] Y. Bugeaud. Distribution modulo one and Diophantine approximation , volume 193 of Cambridge Tracts in Mathematics . Cambridge University Press, Cambridge, 2012.
- 5[5] E. Cateland. Suites digitales et suites k-régulières. Ph D Thesis, Université Bordeaux I, 1992.
- 6[6] M. Drmota, C. Mauduit, and J. Rivat. The Thue-Morse sequence along squares is normal. http://www.dmg.tuwien.ac.at/drmota/alongsquares.pdf , 2013. [Online; first accessed 26.2.2014].
- 7[7] S. Ferenczi. Complexity of sequences and dynamical systems. Discrete Math. , 206(1-3):145–154, 1999. Combinatorics and number theory (Tiruchirappalli, 1996).
- 8[8] N. P. Fogg. Substitutions in dynamics, arithmetics and combinatorics , volume 1794 of Lecture Notes in Mathematics . Springer-Verlag, Berlin, 2002. Edited by V. Berthé, S. Ferenczi, C. Mauduit and A. Siegel.
