Minimal dispersion of large volume boxes in the cube
Abstract
In this note we present a construction which improves the best known bound on the minimal dispersion of large volume boxes in the unit cube. Let . The dispersion of is defined as the supremum of the volume taken over all axis parallel boxes in the cube which do not intersect . The minimal dispersion of points in the cube is defined as the infimum of the dispersion taken over all such that . Define the “large volume” regime as the class of all volumes . The inverse of the minimal dispersion is denoted as . When the volume is large, the best known upper bound on is of the order The construction presented in this note yields an upper bound given by
Some of our intermediate estimates are sharp given the condition that , where is a positive constant which depends only on the volume .
1 Introduction
The dispersion of is defined as the supremum of the volume taken over all axis parallel boxes in the cube which do not intersect . The minimal dispersion of points in the cube is defined as the infimum of the dispersion taken over all such that . The inverse of the minimal dispersion, denoted as , is the size of the smallest set of points which intersects any axis parallel box with volume exceeding . The problem of estimating the minimal dispersion (which was originally defined in [19] as modification of a concept in [6]) has been given considerable attention in recent years. Some contemporary works such as [1], [3], [5], [7], [8] [10], [11], [12], [13], [14], [15], [17], [21], [22], are focused on this problem, or modifications (general -dispersion, spherical dispersion, etc.) thereof. We will refer to these works when we discuss the known results, and the best known bounds on the minimal dispersion. The dispersion of some particular sets has been studied in [9], [18], [20]. Dispersion is of interest in studying topics in subjects such as Discrete Geometry and Approximation Theory and in studying particular random point configurations as in [15]. Define the “large volume” regime as the class of all volumes . When the volume is large, the best known upper bound given in [17], is of the order In this note we present a construction which improves the best known bound in this setting. First, consider the situation when . The minimal dispersion is attained with one point at the center of the cube. Thus, it is clear that in the large volume regime we are interested in estimating the minimal dispersion when . Theorem 1.1 improves the best known bound previously given by Sosnovec in [17].
Theorem 1.1.
Let , and let . Then
(1) |
Let , and let . We construct a set of points on the diagonal in the following way. Consider the fractional sequence
Denote the subsequent elements in the sequence as . We show that there exists a smallest number such that
Consider the configuration given by
We can see some examples of the configurations when (respectively). They are given below.
Define the function
(2) |
where denotes the cardinality.
It may not be obvious yet, but we can see what the function looks like (a monotone decreasing step function) on some inclusion-ordered intervals in the examples given below.
We obtain Theorem 1.1 by showing that this function is an upper bound for the minimal dispersion.
Finally, we show that some of our estimates are sharp given that , where is a positive constant dependent on .
2 Previous results
Recall that the inverse of the minimal dispersion, denoted as is defined as the smallest number of points that are needed to intersect any axis parallel box whose volume exceeds .
In their work [1], Aistleitner, Hinrichs, and Rudolf proved a lower bound for given by
(3) |
This result gives a non-trivial estimate showing that the dispersion asymptotically increases with dimension when the volume is not large. The upper bound given by
(4) |
was proved by Larcher, and is also presented in [1]. This is an improvement on the bound given by Rote and Tichy in [15]. The upper bound given by
(5) |
is a consequence of a more general result given in [2]. The authors present an argument which uses the dimension of (which is ) instead of the ambient dimension . In the paper [16], Rudolf presented a probabilistic argument which yields the bound in (5). The bound in (5) is an improvement on (4) under the assumption that where is an absolute constant. Sosnovec in [17] obtained another upper bound which is better when grows to infinity, namely
(6) |
The constant obtained in [17] grows extremely fast with . This constant was improved by Ullrich and Vybìral [21] who showed that
Litvak in [11] gives an improvement on the known bounds when , and showed that
Litvak also established that when that
which is an improvement on the bound given by Ullrich and Vybìral. Recently, Litvak and Livshyts in [12], proved an upper bound for , and given by
(7) |
They also showed that a random choice of points with respect to the uniform distribution in the cube gives the result with high probability. We mention here that Hinrichs, Krieg, Kunsch, and Rudolph in [7] showed that using a random choice of points uniformly distributed within the cube, one cannot expect to get anything better than
Thus we can see that in the probabilistic setting, (7) is close to the best possible. Now we turn our attention to the large volume regime . Sosnovec in [17] gave a dimension independent upper bound for . In particular, he proved that for , and ,
The goal of this paper is to improve this bound. This is established in Theorem 1.1.
2.1 Definitions and Notation
Let . Denote unit cube as . Let denote the cardinality, let denote the domain, and let . The set of axis parallel boxes is defined as
The dispersion of is defined as
The minimal dispersion is defined as
Let . The inverse of the minimal dispersion is defined as
Let . Inductively define a sequence of functions in the following way. Let . Define as the identity . Given functions , and numbers , define
(8) |
and define by
(9) |
Proposition 3.2 shows that the infimum in (8) is attained for all . Let . It is clear that , hence . Let (we include the the right endpoint here since the intermediate estimates are still valid here). The following definition is a complete description of the step function introduced in (2). Define
(10) |
and define
(11) |
3 Auxiliary tools
Proposition 3.1.
Let . The function is strictly increasing on its domain.
Proof.
Employ induction on . Let . By definition . For all , . The base case is seen to be trivial. Assume that the proposition holds for . Let be such that . Then by domain inclusion . By the induction hypothesis . It follows from the definition of as in (9) that
This proves the proposition. ∎
Proposition 3.2.
Let . For each there exists a unique number
with the properties
(12) |
(13) |
(14) |
Proof.
Note that . Employ induction on . For each we produce a number with the properties (12), (13), (14).
Let . Recall the definition of given in (9), by
We will show that there exists such that
Equivalently, we find the solution to the equation It is clear that , and that . This implies that Hence, has properties (12), (13), (14).
Let . Recall the definition of given in (9), by
We show that there exists such that
Equivalently, we show that there exists a solution to the equation
The function is strictly increasing by Proposition 3.1. Note that , and . Apply the Intermediate Value Theorem to . This yields a unique solution such that
It follows that , and that Hence has properties (12), (13), (14).
Let . Assume that there exist numbers with the properties (12), (13), (14). Under the assumption that where , and . Recall the definition of given in (9), by
We show that there exists , such that
Equivalently, we show that there exists a solution to the equation
The function is strictly increasing by Proposition 3.1. Note that and by the induction hypothesis, Apply the Intermediate Value Theorem to . This yields a unique solution , such that
It follows that , and that This yields a number with the properties (12), (13), (14). This proves the proposition. ∎
Remark 3.3.
Proposition 3.4.
Let . Let be such that for all . Then for all ,
Proof.
Fix . Employ induction on . Let be such that for all . Let . Recall the definitions of . Since , it follows that
Let . Assume as the induction hypothesis that for all ,
By assumption , then by definition of it follows that
This proves the proposition. ∎
Corollary 3.5.
Let . Let be such that . Then for all ,
Proof.
Remark 3.6.
Corollary 3.5 gives the property that for all , . Property (14) in Proposition 3.2 gives that for all , . Let , and recall the definition in (10) given by
Then we have that for all ,
This gives the integral values of the step function in (2) evaluated at the left endpoints of the interval partition given by
(16) |
The following proposition shows that the function in (10) is constant over any interval in the partition (16).
Proposition 3.7.
Let . Let be such that . Then
Proof.
First, let , and let be such that . Since , it follows that . Hence, .
Let , and let be such that . Since , Corollary 3.5 implies that for all , . Apply Proposition 3.1 on to get . By Proposition 3.2,
It follows that
This means that for all ,
It follows that . Since , apply Proposition 3.1 on to get . Recall that by Proposition 3.2,
It follows that
Thus, . It follows that Therefore, This proves the proposition. ∎
Remark 3.8.
A brief discussion on Geometric Rational Sequences follows. We use the results herein to obtain explicit values for the numbers as in (15). The paper [4] provides results which can be applied to sequences of the form defined below.
Definition 3.9.
Let . A Geometric Rational Sequence is defined by setting an initial condition , and recursively defining
If , then define , , so that .
Definition 3.10.
Let . The reduced form of a Geometric Rational Sequence is defined by setting an initial condition
and recursively defining
If , then define , , so that .
Remark 3.11.
Note that if , then . Hence, and . This occurs if and only if the sequence is cyclical.
Remark 3.12.
Let , and let . Then,
Continuing in this way, we see that for all and ,
Proposition 3.13.
Let . Then the reduced sequence is cyclical with cycle length .
Proof.
From Remark 3.12, we have that for all ,
Apply Proposition 3.4 to yield
for all . This guarantees no repetition in the first terms of the reduced sequence . By Proposition 3.2, . It follows that
Recall the reduced sequence given in Definition 3.10. Then by definition
This shows that . It follows that the sequence is cyclical with cycle length . ∎
The following theorem from [4] will be used. Note that there is a typographical error in the condition .
Theorem 3.14.
Let , with , and . A sequence satisfying , , has a finite or infinite number of cluster points depending on whether or not is rational. Moreover, when is irreducible, the sequence takes on distinct values which are thereafter repeated in this order.
Proposition 3.15.
Let . Let
Then the reduced sequence has cycle length .
Proof.
Remark 3.16.
Proposition 3.17.
Recall the sequence defined in (15). Then for all ,
Proof.
Apply induction on . Let . It is easy to check that
Let . Recall that is decreasing. Hence, . By Proposition 3.2 it follows that . By Proposition 3.15, the sequence has cycle length . From Remark 3.11, it follows that . Then
This implies that . Thus, .
4 Main results
In this section we give an upper bound for the minimal dispersion.
4.1 Upper bound
Let , and let be as in (11). We define the following configurations on the diagonal
(17) |
It is clear that .
We define two fundamental box types (Type 1 and Type 2). The definitions are given below, but first some examples in the plane should be elucidating. One can always think about the dimensional projection of a higher dimensional box onto one of its faces.
(Example of Type 1 boxes, )
![[Uncaptioned image]](https://cdn.awesomepapers.org/papers/971f618e-9ad7-471e-8264-29b581f639f7/x9.png)
![[Uncaptioned image]](https://cdn.awesomepapers.org/papers/971f618e-9ad7-471e-8264-29b581f639f7/x10.png)
(Example of Type 2 boxes, )
![[Uncaptioned image]](https://cdn.awesomepapers.org/papers/971f618e-9ad7-471e-8264-29b581f639f7/x11.png)
Definition 4.1.
Let .
The box is Type , if one of the following conditions holds.
There exists , such that , or such that .
There exists , and , such that .
The box is Type , if
there exist , and , such that , and .
Lemma 4.2.
Let . Let be a box of Type or Type . Then Vol.
Proof.
Let . First assume that is Type . Assume that . Then the lemma trivially holds. Assume that there exists , such that . Since is the smallest integer such that , it follows that
Assume that there exists , and such that . Then
Recall that is the smallest integer such that . Since , it follows that . Then
Hence, if is Type , then Vol.
Let be a Type box. By definition there exist and such that and . Recall the definition of as in (9) given by
It follows that
This proves the lemma. ∎
Lemma 4.3.
Let . Let be such that . Then is either Type or Type .
Proof.
Let
(18) |
Let . Notice that for all . Then it is clear that there exists such that . Define
If , then is Type . From here we list the remaining possible cases.
Case 1: In this case . Then . It is clear that . Since , there exists such that or . If , then by definition is Type . If , then since , by definition is Type .
Case 2: Denote . In the second case let
and define
The following algorithm shows that is Type or Type . Let . Then it is clear that . Recall that . Then there exists , such that or . If , then is Type . If , then since , is Type . If , then denote . Let
and define
If , then appeal to Case 1. Since there exists , such that or . If , then is Type . If , then since is Type . If , then appeal to Case . If denote , and continue the algorithm. At step of the algorithm let
and define
If , then is Type . Assume . If , then appeal to Case . If there exists an interval , such that , then since by definition is Type . The algorithm terminates after at most steps.
Case 3: In the last case . Since there exists or . If , then is Type . Finally, if , then apply the results in Case and Case . This proves the lemma. ∎
Corollary 4.4.
Let . Then
Proof.
Corollary 4.4 gives an upper bound for the minimal dispersion.
4.2 Formula derivation
We derive a simple formula for the upper bound given in Corollary 4.4.
Corollary 4.5.
Let . Then
Proof.
Theorem 4.6.
Let . Then
5 Sharp estimates
We prove some results about our configurations on the diagonal and diagonal analogues in the cube.
5.1 The diagonal
Proposition 5.1.
Let , and let be on the diagonal such that
Then
Proof.
Let . Let be on the diagonal such that
Let be as in (11), and let be the configuration as defined in (17). Let be as defined in (18). Let
such that
Assume toward a contradiction that . Partition into intervals
(20) |
Since , there exist at least two intervals which do not intersect . Assume . Then . Construct the box
Then . However, Vol which contradicts the assumption that disp.
Now assume
Then . Remove from the partition in (20). Then two intervals in the partition
must not intersect . If there exists such that , then construct the box
It follows that , however,
This contradicts the assumption that disp. Let be the smallest integer such that
Assume is the smallest integer such that . Construct the box
Since is on the diagonal , however,
This contradicts the assumption that disp.
The proposition follows. ∎
Proposition 5.2.
Proof.
Let . Let be as defined in (15). Employ an inductive argument on . Let . Then by Proposition 3.2 it follows that
Let . Then by definition of the functions, and by Proposition 3.2
By Proposition 3.2, . It follows that
Then
in particular
Fix , and assume the induction hypothesis, for all . Namely that
We show that
By the induction hypothesis
By construction of the functions, we have that
Therefore,
It follows that
The proposition follows. ∎
5.2 Extended diagonal
Definition 5.3.
Let , and . Let be given by (18), and denoted as . Define the Extended Diagonal as
Definition 5.4.
Let . Let . Define
Proposition 5.5.
Let , and . Let be such that
Then
Proof.
Assume that the hypothesis holds. Let be as in (11), and let be the configuration as in (17). Define
For define
Assume toward a contradiction that . The points contained in must lie in the disjoint sets in composing . There exist at least two integers , such that .
First assume that . Since , it follows that for each . Let
Construct the box . The magnitude of the components of each is bounded below by . It follows that , however, Vol. This contradicts the assumption that .
Now assume that . Let be the smallest number such that Assume that for all , . Then set Construct a box
Since , it follows that , however,
This contradicts the assumption that . Assume that for some smallest . Then set
Construct a box
Since , it follows that , however,
This contradicts the assumption that . It follows that
∎
Remark 5.6.
Proof.
Assume the hypothesis. Note . For all define to be as in the proof of Proposition 5.5. Assume that for all , . Then it follows that for all , . Hence, . We state here that for all , . This follows directly from the proof of Proposition 5.5: thus we shall omit the details to avoid repetition. Then for all , . Now apply induction on to show that for all , .
Let . Assume toward a contradiction that
Let . There exists a maximum component of which is less than . Without loss of generality, assume that this is the first component. Denote the magnitude of this component as . Construct the box
Since , it follows that . However,
This contradicts (21). Therefore, .
Fix . As an induction hypothesis assume that for all , . Assume toward a contradiction that . Let . There exists a maximum component of which is less than . Without loss of generality, assume that this is the first component. Denote the magnitude of this component as . Construct the box
Since , it follows that . Then
This contradicts (21). Therefore, it follows that . Hence, for all , . It follows that . ∎
Now we show that the bound in Corollary 4.5 is sharp, given that is large enough. Recall that for each , there exists such that , and
Theorem 5.8.
Let . Then
given that , where .
Proof.
Let , then for all ,
Let be such that . Let . Define , , and denote . Assume toward a contradiction that there exists such that
Then either or . Assume . Construct the box . Then , and
This contradicts the assumption that Assume . Construct the box . Then , however,
This contradicts the assumption that Then and by Corollary 4.4, It follows that
Let be such that . Since , . Define , , . Select two points . Assume toward a contradiction that . The components of and are contained in the intervals . Denote the components of as , . Let denote the largest number of components of contained in a single interval, denoted as . Denote the corresponding indices as Let denote the largest number of components of contained in a single interval, denoted as . Denote the corresponding indices as . It follows that
Project onto the components,
and
The points are projected onto a -dimensional face of , given by
Note that
and
The components are contained in . By Proposition 5.7 there exists such that
where . However, . Let
It is clear that . However, . This contradicts the assumption that . Then , and by Corollary 4.4, It follows that
Fix , and let . Since ,
Let be arbitrary points in the cube. Assume toward a contradiction that Define the partition
For all , denote the components of as . Denote . Let
denote the largest number of components of which are contained in a single interval, denoted as . Denote the corresponding indices as
Let
denote the largest number of components of which are contained in a single interval, denoted as . Denote the corresponding indices as .
For all let
denote the largest number of components of which are contained in . Denote the corresponding indices as . This guarantees that . Define a projection on , such that
This embeds into . Then by Proposition 5.7, there exists a box such that with , such that
This can be extended to a box
It is clear that and that . This contradicts the assumption that . Then , and by Corollary 4.4, It follows that
∎
6 Concluding remarks
When , and is small the following configurations are better than the best known bound which asymptotically is . We present a configuration which is easy to describe and visualize.
Proposition 6.1.
Let . Then .
Proof.
Let let
Every box inside is contained in one of the following
(22) |
This gives the result in the dimensional case. Let . Let
Let . Then each box in such that is contained in a product of intervals with one of the boxes in (22). Each of the boxes have volume . Therefore,
and ∎
6.1 Acknowledgments
The author would like to thank his supervisor A.E. Litvak for introducing him to the dispersion problem, for his valuable expertise, and helpful comments while the author was working on this note. The author would also like to thank the reviewers for their comments and suggestions.
References
- [1] C. Aistleitner, A. Hinrichs, D. Rudolf, On the size of the largest empty box amidst a point set, Discrete Appl. Math. 230 (2017), 146-150.
- [2] A. Blumer, A. Ehrenfeucht, D. Haussler, M. Warmuth, Learnability and the Vapnik–Chervonenkis dimension, J. Assoc. Comput. Mach. 36 (1989), 929-965.
- [3] S. Breneis, A. Hinrichs Fibonacci lattices have minimal dispersion on the two-dimensional torus preprint, (2019), arXiv:1905.03856
- [4] L. Brand, A Sequence Defined by a Difference Equation, The American Mathematical Monthly. 62 (7) (1955), pp. 489-492.
- [5] A. Dumitrescu, M. Jiang, On the largest empty axis-parallel box amidst n points, Algorithmica. 66 (2013), 225-248.
- [6] E. Hlawka, Abschatzung von trigonometrischen Summen mittels diophantischer Approximationen, Osterreich. Akad. Wiss. Math.- Naturwiss. Kl. S.-B. II, 185 (1976), 43–50.
- [7] A. Hinrichs, D. Krieg, R.J. Kunsch, D. Rudolf, Expected dispersion of uniformly distributed points, J. Complexity, to appear.
- [8] A. Hinrichs, J. Prochno, M. Ullrich, J. Vybìral, The minimal k-dispersion of point sets in high dimensions, J. Complexity 53 (2019), 68-78.
- [9] D. Krieg, On the dispersion of sparse grids, J. Complexity 45 (2018), 115–119.
- [10] D. Krieg, D. Rudolf, Recovery algorithms for high-dimensional rank one tensors, Journal of Approximation Theory 237 (2019): 17-29.
- [11] A.E. Litvak, A remark on the minimal dispersion, Communications in Contemporary Mathematics, to appear.
- [12] A.E. Litvak, G. V. Livshyts, New bounds on the minimal dispersion, J. Complexity, to appear.
- [13] E. Novak, D. Rudolf, Tractability of the approximation of high- dimensional rank one tensors Constructive Approximation 43.1 (2016): 1-13.
- [14] J. Prochno, D. Rudolph, The minimal spherical dispersion preprint, (2021), arXiv:2103.11701
- [15] G. Rote, R.F. Tichy, Quasi-Monte Carlo methods and the dispersion of point sequences, Math. Comput. Modelling. 23 (1996), 9–23.
- [16] D. Rudolf, An upper bound of the minimal dispersion via delta covers, Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan, Springer-Verlag. (2018), 1099-1108.
- [17] J. Sosnovec, A note on the minimal dispersion of point sets in the unit cube, European J. of Comb. 69 (2018), 255–259.
- [18] V.N. Temlyakov, Dispersion of the Fibonacci and the Frolov point sets, preprint, (2017), arXiv:1709.08158.
- [19] M. Ullrich, A lower bound for the dispersion on the torus, Mathematics and Computers in Simulation, 143 (2018), 186–190.
- [20] M. Ullrich, A note on the dispersion of admissible lattices, Discrete Appl.Math, 257 (2019), 385–387.
- [21] M. Ullrich, J. Vybìral, An upper bound on the minimal dispersion, Jour- nal of Complexity. 45 (2018), 120–126.
- [22] M. Ullrich, J. Vybìral, Deterministic constructions of high-dimensional sets with small dispersion, Preprint, (2019), arXiv:1901.06702
Kurt S. MacKay
Dept. of Math. and Stat. Sciences,
University of Alberta,
Edmonton, AB, Canada, T6G 2G1.
e-mail: [email protected]