Deterministic constructions of high-dimensional sets with small dispersion

TL;DR

This paper presents explicit deterministic algorithms for constructing low-dispersion point sets in high dimensions, leveraging coding theory, with applications in numerical integration and discrepancy theory.

Contribution

It introduces a novel deterministic construction method for low-dispersion sets based on coding theory, improving over randomized approaches.

Findings

Constructs low-dispersion point sets deterministically

Algorithms run in polynomial time in dimension d

Time complexity is super-exponential in 1/ε

Abstract

The dispersion of a point set $P\subset[0,1]^d$ is the volume of the largest box with sides parallel to the coordinate axes, which does not intersect $P$. Here, we show a construction of low-dispersion point sets, which can be deduced from solutions of certain $k$-restriction problems, which are well-known in coding theory. It was observed only recently that, for any $\varepsilon>0$, certain randomized constructions provide point sets with dispersion smaller than $\varepsilon$ and number of elements growing only logarithmically in $d$. Based on deep results from coding theory, we present explicit, deterministic algorithms to construct such point sets in time that is only polynomial in $d$. Note that, however, the running-time will be super-exponential in $\varepsilon^{-1}$.

Equations

$$\mathrm{disp}(P):=\underset{B:B\cap P=\varnothing }{\sup }|B|,$$

$$\mathrm{disp}(n,d):=\underset{\begin{array}{c}P\subset [0,1{]}^d:\\ \#P=n\end{array}}{\inf }\mathrm{disp}(P)$$

$$N(\varepsilon,d)\;:=\;\min\Bigl{\{}n\colon\operatorname{disp}(n,d)\leq\varepsilon\Bigr{\}}.$$

$$\frac{{\log }_2(d)}{8\epsilon }\le N(\epsilon ,d)\le \frac{C^d}{\epsilon }$$

$$\frac{\log_{2}(d)}{8\varepsilon}\;\leq\;N(\varepsilon,d)\;\leq\;2^{7}\,\log_{2}(d)\,\frac{\bigl{(}1+\log_{2}(\varepsilon^{-1})\bigr{)}^{2}}{\varepsilon^{2}},$$

$$\#P\le C{\left(\frac{1+{\log }_2({\epsilon }^{-1})}{\epsilon }\right)}^6\log (d^{\ast }),$$

$$c=c(\mathcal{C}):=\underset{1\le j\le M}{\min }\#C_j.$$

$$\mathbb{P}(v\text{ satisfies }{C}_j\text{ at }S)=\frac{\#C_j}{b^k}\ge \frac{c}{b^k}.$$

$${\mathbb{P}}(\text{no }v_{1},\dots,v_{N}\text{ satisfies }C_{j}\text{ at }S)\,=\,\Bigl{(}1-\frac{\#C_{j}}{b^{k}}\Bigr{)}^{N}\leq\Bigl{(}1-\frac{c}{b^{k}}\Bigr{)}^{N}.$$

$$\binom{n}{k}\cdot M\cdot\Bigl{(}1-\frac{c}{b^{k}}\Bigr{)}^{N}.$$

$$N\,>\,\frac{\log(n^{k}M)}{\log\bigl{(}\frac{b^{k}}{b^{k}-c}\bigr{)}}.$$

$$N\ge \lceil \frac{b^k}{c}\log (n^{k}M)\rceil ,$$

$$\lceil \frac{b^k}{c}\log (n^{k}M)\rceil ,$$

$$\mathcal{O}\left(\frac{b^k}{c}{\left(\frac{e{n}^2}{k}\right)}^{k}M\mathcal{T}\right),$$

$$R\ge {\gamma }_0\left(1-H_{k^2}\left(1-\frac{1,5}{k^2}\right)\right)\left(1-\frac{1-2/k^2}{1-1,5/k^2}\right),$$

$$\mathcal{O}\bigl{(}k^{4}\log(n)\bigr{)}$$

$$\#P\le c_{\epsilon }\log (d).$$

$${\mathbb{M}}_{m}=\Bigl{\{}\frac{1}{2^{m}},\dots,\frac{2^{m}-1}{2^{m}}\Bigr{\}}.$$

$$\Omega_{m}:=\Big{\{}\mathbb{B}=I_{1}\times\dots\times I_{d}\subset[0,1]^{d}\,:\,\mathrm{vol}(\mathbb{B})>\frac{1}{2^{m}}\Big{\}}$$

$$\mathcal{A}(\mathbb{B})=\{j\in \{1,\dots ,d\}:{\mathbb{M}}_m\not\subset I_j\},$$

$$\left(\genfrac{}{}{0ex}{}{d}{A_m}\right)(2^m-1{)}^{A_m}[1-(2^m-1{)}^{-A_m}{]}^N.$$

$$\#P\le c_{\epsilon }{\log }_2(d).$$

$$N=2^{m{A}_m+O({\log }^2(m{A}_m))}\log (d)\le 2^{m^2{2}^m+O(m^2)}{\log }_2(d).$$

$$\#P\le 2^7{\log }_2(d)\frac{(1+{\log }_2({\epsilon }^{-1}){)}^2}{{\epsilon }^2}.$$

$${\Omega }_m(s,p)$$

$$\displaystyle\qquad\qquad\text{and}\quad\inf I_{\ell}\in\Big{[}p_{\ell}-\frac{1}{2^{m}},p_{\ell}\Big{)}\bigg{\}}.$$

$$\#{\mathbb{I}}_{m}\leq\exp\Bigl{(}m2^{m}\log(2^{m+3}d)\Bigr{)},$$

$$\mathbb{B}_{m}(s,p)\,:=\,\bigcap_{\mathbb{B}\in\Omega_{m}(s,p)}\mathbb{B}\,=\,\prod_{\ell=1}^{d}\Bigl{[}p_{\ell},p_{\ell}+\frac{s_{\ell}-1}{2^{m}}\Bigr{]}.$$

$${\mathbb{P}}\big{(}z\in\mathbb{B}_{m}(s,p)\big{)}\geq\frac{1}{2^{m+4}}\,.$$

$${\mathbb{P}}\bigl{(}\forall\ell\in\{1,\dots,N\}\colon x^{\ell}\not\in{\mathbb{B}}_{m}(s,p)\bigr{)}\,\leq\,(1-2^{-m-4})^{N}.$$

$${\mathbb{P}}\bigl{(}\exists(s,p)\in{\mathbb{I}}_{m}\colon\forall\ell\in\{1,\dots,N\}\colon x^{\ell}\not\in{\mathbb{B}}_{m}(s,p)\bigr{)}\,\leq\,\#{\mathbb{I}}_{m}\cdot(1-2^{-m-4})^{N}.$$

Full text

It was observed only recently that, for any $\varepsilon>0$, certain randomized constructions provide point sets with dispersion smaller than $\varepsilon$ and number of elements growing only logarithmically in $d$. Based on deep results from coding theory, we present explicit, deterministic algorithms to construct such point sets in time that is only polynomial in $d$. Note that, however, the running-time will be super-exponential in $\varepsilon^{-1}$.

We would like to express our gratitude to the Erwin Schrödinger International Institute for Mathematics and Physics for its hospitality during the programme on “Tractability of High Dimensional Problems and Discrepancy”, where some part of this research was carried out. We also gratefully acknowledge the support of the Oberwolfach Research Institute for Mathematics, where initial discussion were held during the workshop “Perspectives in High-Dimensional Probability and Convexity”. We thank also Michael Gnewuch, Daniel Král’ and Hemant Tyagi for fruitful discussions. JV was supported by the grant P201/18/00580S of the Grant Agency of the Czech Republic and by the European Regional Development Fund-Project “Center for Advanced Applied Science” (No. CZ.02.1.01/0.0/0.0/16_019/0000778)

1. Introduction and Results

For $d\in{\mathbb{N}}$ and a point set $P\subset[0,1]^{d}$ we define the dispersion of $P$ by

[TABLE]

where the supremum is over all axis-parallel boxes $B=I_{1}\times\dots\times I_{d}$ with intervals $I_{\ell}\subset[0,1]$, and $|B|$ denotes the (Lebesgue) volume of $B$. As we are interested in point sets which make the above quantity as small as possible, we additionally define, for $n,d\in{\mathbb{N}}$, the $n$-th-minimal dispersion

[TABLE]

and, for $\varepsilon\in(0,1)$, its inverse function

[TABLE]

Hence, $N(\varepsilon,d)$ is the minimal cardinality of a point set $P\subset[0,1]^{d}$ that has dispersion smaller than $\varepsilon$.

Besides the fact that the above geometric quantities are interesting in its own right, they also attracted attention in the numerical analysis community, especially when it comes to very high dimensional applications. The reason is that bounds on the (minimal) dispersion lead to bounds on worst-case errors, and hence the complexity, for some numerical problems, including optimization and approximation in various settings, see [5, 27, 30, 32, 36, 37, 40]. This is a similar situation as for the much more studied discrepancy, which corresponds to certain numerical integration problems, see e.g. [8, 9, 12, 28, 29, 31].

Moreover, bounding the dispersion is clearly also related to the problem of finding the largest empty box. In dimension two, this is the Maximum Empty Rectangle Problem, which is one of the oldest problems in computational geometry. For the state of the art and further references we refer to [13, 14, 15, 16, 24].

Regarding bounds on the inverse of the minimal dispersion, we have

[TABLE]

for some $C<\infty$ and all $\varepsilon<1/8$, see [1], where the upper bound is attained by certain digital nets. See also [14, 32] for a related upper bound. Although these bounds show the correct dependence on $\varepsilon$, they are rather bad with respect to the dimension $d$. This gap was narrowed in the past years by several authors, see  [21, 33, 35, 39], including the important work of Sosnovec who proved that the logarithmic dependence on $d$ is optimal. With this respect, the best bound at present is

[TABLE]

see [39] and Theorem 4.1 below. Note that the logarithmic dependence is special for the cube, as it is known that, for the same problem on the torus, we have a lower bound linear in $d$, see [38].

The main drawback of these results is that they only show the existence of point sets with small dispersion. The only explicit constructions we are aware of are the above mentioned digital nets, which lead to a bad $d$-dependence, and sparse grids, which satisfy the upper bound $N(\varepsilon,d)\leq(2d)^{\log_{2}(1/\varepsilon)}$, see [21]. It is already clear from the number of points, that both do not lead to constructions of point sets with small dispersion that can be carried out in time that is polynomial in $d$. Moreover, as the existence proofs are based on random points on a finite grid, one could use the naive algorithm: try each of the possible configurations, calculate its dispersion (if possible) and output a set that satisfies the requested bound. The running-time of this algorithm is (in the worst case) clearly at least exponential in $d$.

Remark 1.1**.**

Note that the decision problem, if a given point set has discrepancy smaller than $\varepsilon$, is known to be NP-hard, see [18] or [11, Section 3.3], and the same is true for the dispersion, if the dimension $d$ is part of the input, cf. [6].

For upper bounds on the cost of constructing points with small discrepancy and further literature, see [10, 11, 17]. However, all known algorithms for this problem so far have running-time at least exponential in $d$. We hope that the results of this paper will lead to some progress also for this problem.

Reconsidering the bound (1.1), one may hope that points with small dispersion may be constructable also in very high dimensions. Ideally, we would like to have algorithms for the construction of point sets of size $N=\mathcal{O}(N(\varepsilon,d))$ with dispersion at most $\varepsilon>0$, whose computational cost is polynomial in $\varepsilon^{-1}$ and $d$. However, this seems to be out of reach. (Note that the output already costs $d\cdot N$.)

Here, we focus on the dependence on $d$ and show, using deep results from the theory of error-correcting codes, that point sets with small dispersion of size $\sim\log(d)$ can be constructed in time that is polynomial in $d$. Unfortunately, we do not have a good control on the dependence on $\varepsilon^{-1}$. It remains an open problem to find fully-polynomial constructions for the dispersion.

Our two main results are derandomized versions of the results from [35] and [39]. Both use different approaches and lead to somewhat different results. The first one, as discussed in Section 3, leads to point sets of size $\mathcal{O}_{\varepsilon}(\log d)$ that can be constructed in time linear in the size of the output, i.e., in time $\mathcal{O}_{\varepsilon}(d\log d)$, see Algorithm 1 and Theorem 3.3. Here, and in the following, $\mathcal{O}_{\varepsilon}(\log d)$ means that the implied constant depends on $\varepsilon$ in an unspecified way. A second, and much more involved, construction will be given by Algorithm 2 in Section 4. The corresponding result reads as follows.

Theorem 4.4 Let $\varepsilon\in(0,\frac{1}{4}]$ and $d\geq 2$. Then, there is an absolute constant $C<\infty$, such that Algorithm 2 constructs a set $P\subset[0,1]^{d}$ with $\operatorname{disp}(P)\leq\varepsilon$ and

[TABLE]

where $d^{*}:=\max(d,2/\varepsilon)$.

*The running-time of Algorithm 2 is $\mathcal{O}_{\varepsilon}(d^{c})$ for some $c<\infty$. *

Note that, in contrast to Algorithm 1, the point set that is constructed by Algorithm 2 has size that is polynomial in $\varepsilon^{-1}$. However, as the proof shows, its computational cost is much larger.

Our general approach is as following. We start with a detailed inspection of the random constructions of point sets with small dispersion. This will allow us to clearly separate the setting of the construction and its randomized part. It will turn out then, which properties are crucial for each of the approaches and what exactly is the role of randomness. Afterwards, we replace the randomized part by a deterministic one.

2. Basics from Coding Theory

Our deterministic constructions of point sets with small dispersion will be essentially obtained by certain “derandomization” of recent proofs from [35, 39]. As we will rely on rather deep tools from coding theory, we summarize in this section the necessary definitions and results for later use.

2.1. Universal sets

First, we introduce the concept of $(n,k)$-universal sets, which is also known in coding theory under the name of $t$-independent set problem. It has its roots in testing of logical circuits [34].

Definition 2.1** ($(n,k)$-universal sets).**

Let $1\leq k\leq n$ be positive integers. We say that $T\subset\{0,1\}^{n}$ is an $(n,k)$-universal set, if for every index set $S\subset\{1,\dots,n\}$ with $\#S=k$, the projection of $T$ on $S$ contains all possible $2^{k}$ configurations.

Naturally one is interested in (randomized and deterministic) constructions of small $(n,k)$-universal sets. The straightforward randomized construction provides the existence of an $(n,k)$-universal set of size $\lceil k2^{k}\log(n)\rceil$. On the other hand, [20] gives a lower bound on the size of an $(n,k)$-universal set of the order $\Omega(2^{k}\log(n))$. There exist several deterministic constructions of $(n,k)$-universal sets in the literature (cf. [2, 3, 25]) and we shall rely on the results given in [26].

Theorem 2.2** ([26, Theorem 6]).**

There is a deterministic construction of an $(n,k)$-universal set of size $2^{k+O(\log^{2}(k))}\log(n)$, which can be listed in linear time of the length of the output.

Although the notion of an $(n,k)$-universal set is not very flexible and comes from a different area of mathematics, we will see in Section 3, that there is indeed a link to sets with small dispersion. Reusing the known results from coding theory, it will already allow us to obtain our first deterministic construction of a point set with cardinality of order $\log(d)$. However, in this approach we have only very limited control of the dependence on $\varepsilon$.

We also need the following natural generalization of $(n,k)$-universal sets. We could not find this concept in the literature, but we assume that this and the proceeding lemma are known.

Definition 2.3** ($(n,k,b)$-universal sets).**

Let $1\leq k\leq n$ and $b\geq 2$ be positive integers. We say that $T\subset\{0,1,\dots,b-1\}^{n}$ is an $(n,k,b)$-universal set, if for every index set $S\subset\{1,\dots,n\}$ with $\#S=k$, the projection of $T$ on $S$ contains all possible $b^{k}$ configurations.

If $b=2$, Definitions 2.1 and 2.3 coincide, i.e. $(n,k,2)$-universal sets are just the usual $(n,k)$-universal sets. We use the following two observations to transfer the known results about $(n,k)$-universal sets to our setting.

Lemma 2.4**.**

Let $1\leq k\leq n$ and $b\geq 2$ be positive integers.

1. $(i)$

Let $T\subset\{0,1,\dots,b\}^{n}$ be an $(n,k,b+1)$-universal set. Then there is an $(n,k,b)$-universal set $T^{\prime}\subset\{0,1,\dots,b-1\}^{n}$ of at most the same size. 2. $(ii)$

Let $m\in{\mathbb{N}}$ and $T\subset\{0,1\}^{mn}$ be an $(mn,mk)$-universal set. Then there is an $(n,k,2^{m})$-universal set of the same size.

Proof.

The proof is quite straightforward. To show $(i)$, just replace all occurrences of $b$ among the coordinates of $T$ by zero. For the proof of the second part it is enough to interpret each $x\in T\subset\{0,1\}^{mn}$ as a digital representation of $\tilde{x}\in\{0,1,\dots,2^{m}-1\}^{n}$. ∎

The direct random construction yields the existence of an $(n,k,b)$-universal set of the size $\lceil kb^{k}\log(ebn/k)\rceil$. Using Theorem 2.2, we can easily obtain a deterministic construction of an $(n,k,b)$-universal set of only a slightly larger size.

Theorem 2.5**.**

There is a deterministic construction of an $(n,k,2^{m}-1)$-universal set of size $2^{mk+O(\log^{2}(mk))}\log(n)$, which can be listed in linear time of the length of the output.

Proof.

By Theorem 2.2, there is a construction of an $(mn,mk)$-universal set with at most $2^{mk+O(\log^{2}(mk))}\log(mn)$ elements. The result then follows from Lemma 2.4. ∎

2.2. $k$-restriction problems

For the derandomization of the analysis of [39] we need a more flexible notion of the so-called $k$-restriction problems, see [26, Section 2.2]. Solutions to these problems will be one of the building blocks of our deterministic construction of sets with small dispersion whose size is polynomial in $1/\varepsilon$ and, still, logarithmic in $d$.

Definition 2.6** ($k$-restriction problems).**

Let $b,k,n,M$ be positive integers and let ${\mathcal{C}}=\{C_{1},\dots,C_{M}:C_{i}\subset\{0,1,\dots,b-1\}^{k}\}$ be invariant under the permutations of the index set $\{1,\dots,k\}$. We say that $T=\{x^{1},\dots,x^{N}\}\subset\{0,1,\dots,b-1\}^{n}$ satisfies the $k$-restriction problem with respect to ${\mathcal{C}}$, if for every $S\subset\{1,\dots,n\}$ with $\#S=k$ and for every $j\in\{1,\dots,M\}$, there exists $x^{\ell}\in T$ with $x^{\ell}|_{S}\in C_{j}.$

Definitions 2.1 and 2.3 are indeed special cases of Definition 2.6. To show this, let us choose $M=b^{k}$ and let $C_{1},\dots,C_{M}$ be all the different singleton subsets of $\{0,\dots,b-1\}^{k}$. Then, $T$ satisfies the $k$-restriction problem with respect to ${\mathcal{C}}$ if for every index set $S\subset\{1,\dots,n\}$ and every possible $v\in\{0,\dots,b-1\}^{k}$ there exists an $x\in T$ with $x|_{S}=v$, i.e. if the restriction of $T$ to every index set $S$ with $k$ elements attains all $b^{k}$ possible values.

An important parameter of the $k$-restriction problems is the minimal size of each of the restriction sets $C_{j}$, i.e.

[TABLE]

Random constructions of sets satisfying the $k$-restriction problem with parameters $(b,k,n,M)$ and ${\mathcal{C}}=\{C_{1},\dots,C_{M}\}$ are based on a simple union bound. Indeed, let $1\leq j\leq M$ and $S\subset\{1,\dots,n\}$ with $\#S=k$ be fixed. The probability, that a randomly chosen vector $v\in\{0,1,\dots,b-1\}^{n}$ satisfies $C_{j}$ on $S$ is at least

[TABLE]

If we choose $N$ random vectors $v_{1},\dots,v_{N}$ independently, the probability that none of them satisfies $C_{j}$ at $S$ is at most

[TABLE]

Finally, the probability that there is a set $S\subset\{1,\dots,n\}$ with $\#S=k$ and $1\leq j\leq M$, such that no $v_{1},\dots,v_{N}$ satisfies $C_{j}$ at $S$ is, by the union bound, at most

[TABLE]

This expression is smaller than one if

[TABLE]

This means that there exist solutions to a $k$-restriction problem with parameters $(b,k,n,M)$ of size $N$ whenever

[TABLE]

where $c$ is from (2.1). Theorem 1 of [26] states that there is a deterministic algorithm that outputs such a solution of size equaling the union bound. The main idea of its proof is that the random sampling can be replaced by an extensive search through a $k$-wise independent probability space with $n$ random variables with values in $\{1,\dots,b\}.$

Theorem 2.7** ([26, Theorem 1]).**

For any $k$-restriction problem with parameters $(b,k,n,M)$ with $b\leq n$, there is a deterministic algorithm that outputs a collection obeying the $k$-restrictions, with the size of the collection equaling

[TABLE]

where $c$ is from (2.1). The time taken to output the collection is

[TABLE]

where $\mathcal{T}$ is the time complexity of the membership oracle.

Here, the membership oracle is a procedure which, for given $v\in\{1,\dots,b\}^{n}$, $S\subset\{1,\dots,n\}$ with $\#S=k$ and $j\in\{1,\dots,M\}$, outputs if the restriction of $v$ on $S$ belongs to $C_{j}.$ In what follows it can be executed in ${\mathcal{T}}=O(k)$ time.

2.3. Splitters

The last ingredient of our derandomization procedure are splitters. They played a central role in [26] as the basic building blocks of all deterministic constructions given there. Essentially, they allow to split a large problem into smaller problems which can then be treated by the extensive search of Theorem 2.7.

Definition 2.8** ($(n,k,l)$-splitter).**

Let $n,k,l$ be positive integers. An $(n,k,l)$-splitter $H$ is a family of functions from $\{1,\dots,n\}$ to $\{1,\dots,l\}$, such that for every $S\subset\{1,\dots,n\}$ with $\#S=k$ there is $h\in H$, which splits $S$ perfectly. It means that the sets $h^{-1}(\{j\})\cap S$ are of the same size for all $j\in\{1,\dots,l\}$ (or as similar as possible if $l\nmid k$).

Similarly to [26] and [4], we will rely on $(n,k,k^{2})$-splitters. By Definition 2.8, $A(n,k)$ is an $(n,k,k^{2})$-splitter, if it is a collection of mappings $a:\{1,\dots,n\}\to\{1,\dots,k^{2}\}$ such that for every $S\subset\{1,\dots,n\}$ with $\#S=k$, there is an $a\in A(n,k)$, which is injective on $S$.

Small $(n,k,k^{2})$-splitters can be obtained from asymptotically good error correcting codes, see [4, Lemma 3]. Indeed, let $c_{1},\dots,c_{n}\in\{1,\dots,k^{2}\}^{L}$ denote the codewords of an error correcting code of length $n$ over the alphabet $\{1,\dots,k^{2}\}$ and normalized Hamming distance at least $1-\frac{2}{k^{2}}$. This is, $c_{i}$ and $c_{j}$ ($i\neq j$) can be equal on at most $2L/k^{2}$ coordinates. If now $S\subset\{1,\dots,m\}$ with $\#S=k$, then there are $k(k-1)/2$ pairs $(i,j)$ with $i\neq j\in S$. As $\frac{2L}{k^{2}}\cdot\frac{k(k-1)}{2}<L$, there must be a coordinate, where all the codewords differ from each other. Finally, if we consider the mappings $h_{j}:i\to c_{i}(j)$, $i=1,\dots,n$ we observe that $H=\{h_{j}:j=1,\dots,L\}$ is an $(n,k,k^{2})$-splitter of size $L$.

Such explicit codes exist by [3] with $L=\mathcal{O}(k^{4}\log n)$. To see this, note that the rate $R$ of a code as above is defined by $R:=\log_{k^{2}}(n)/L=\log(n)/(L\log(k^{2}))$. By [3, eq. (5)], see also [4, Lemma 3], we obtain that an explicit code with normalized Hamming distance at least $\delta=1-\frac{2}{k^{2}}$ exist with

[TABLE]

where $\gamma_{0}>0$ is an absolute constant and $H_{q}(x):=-x\log_{q}(x)-(1-x)\log_{q}(1-x)+x\log_{q}(q-1)$. Simple computations show that this implies $R\geq c/(k^{4}\log(k^{2}))$ for some $c>0$. This, in turn, implies that we can choose $L=\mathcal{O}(k^{4}\log n)$.

The explicit construction of [3] yields a linear code that is based on a two-fold concatenation code that combines the Wozencraft ensemble, Justesen codes and expander codes, which in turn rely on famous deterministic constructions of expander graphs [23]. This construction is ’uniformly constructive’ (see [3]), i.e., the construction can be done in time growing only polynomially in $n$. Furthermore, the code satisfies [3, eq. (5)] for all $\delta<1-1/q$, where $q=k^{2}$ in our case. Note also that the running-time depends only polynomially also on $k$, cf. [26, Thm. 3(iv)]. For the details we refer to Sections 3 and 4 of [3].

The following lemma summarizes the discussion above and shows that $(n,k,k^{2})$-splitters of relatively small size can be constructed explicitly in polynomial time.

Lemma 2.9** (cf. [4, Lemma 3]).**

There is an explicit $(n,k,k^{2})$-splitter of size

[TABLE]

that can be constructed in polynomial time in $n$ and $k$.

The $(n,k,k^{2})$-splitters can be used whenever $n$ is “very large” compared to $k$. Roughly speaking, and in the context of the present paper, we will transform a $k$-restriction problem of size $n$ to a $k$-restriction problem of size $k^{2}$, which can then be solved using the results of Section 2.2. This will lead to construction algorithms with an apparently optimal dependence of their running time on the original problem size $n$. This approach was already used to prove [26, Theorem 6], see Theorem 2.5.

3. Derandomization of Sosnovec’s proof

First, we consider the construction of Sosnovec [35], which gives logarithmic dependence of $N(\varepsilon,d)$ on $d$ but involves no special control of its dependence on $\varepsilon^{-1}$. His main theorem was the following.

Theorem 3.1**.**

([35, Theorem 2]) For every $\varepsilon\in(0,\frac{1}{4}]$, there exists a constant $c_{\varepsilon}>0$, such that for every $d\geq 2$ there is a point set $P\subset[0,1]^{d}$ with $\operatorname{disp}(P)\leq\varepsilon$ and

[TABLE]

We will see that this result can be essentially derandomized using results from coding theory, while loosing only a negligible factor. The drawback of this approach is the extremely bad dependence of $c_{\varepsilon}$ on $\varepsilon$. We sketch the main ideas of Sosnovec’s proof.

3.1. Sosnovec’s proof - the setting

For an integer $m\geq 2$ with $2^{-m}\leq\varepsilon<2^{-m+1}$, we define

[TABLE]

The point set constructed will be a subset of ${\mathbb{M}}_{m}^{d}$. Furthermore, we define

[TABLE]

to be the set of all boxes with sides parallel to the coordinate axes and volume larger than $2^{-m}.$

Let ${\mathbb{B}}=I_{1}\times\dots\times I_{d}\in\Omega_{m}$. The key observation of [35] is that the number of indices $j\in\{1,\dots,d\}$ with ${\mathbb{M}}_{m}\not\subset I_{j}$ is bounded from above by $m2^{m}$, a quantity independent on $d$. To be more specific, if we denote

[TABLE]

then $\#{\mathcal{A}}({\mathbb{B}})\leq A_{m}:=\min(m2^{m},d)$ for every $\mathbb{B}\in\Omega_{m}$. We will refer to ${\mathcal{A}}({\mathbb{B}})$ as the set of “active indices” of ${\mathbb{B}}.$ If ${\mathcal{A}}({\mathbb{B}})$ is not of the full possible size, we enlarge it by adding any of the other indices to obtain a set with cardinality equal to $A_{m}$. Therefore, we can associate to each ${\mathbb{B}}\in\Omega_{m}$ (possibly in a non-unique way) a set ${\mathcal{A}}$ with $\#{\mathcal{A}}=A_{m}$ and a vector $z\in{\mathbb{M}}_{m}^{A_{m}}$ such that any $x\in{\mathbb{M}}_{m}^{d}$ with $x|_{\mathcal{A}}=z$ lies in ${\mathbb{B}}.$

Vice versa, if we have a point set $P=\{x^{1},\dots,x^{N}\}\subset{\mathbb{M}}_{m}^{d}$, such that for every ${\mathcal{A}}\subset\{1,\dots,d\}$ with $\#{\mathcal{A}}=A_{m}$ and to every $z\in{\mathbb{M}}_{m}^{A_{m}}$, there is some $x^{j}\in P$ with $x^{j}|_{\mathcal{A}}=z$, then, by what we just said, $P$ intersects every ${\mathbb{B}}\in\Omega_{m}$. Therefore, the dispersion of $P$ can not be larger than $2^{-m}$, i.e., $\operatorname{disp}(P)\leq 2^{-m}$ and hence $N(2^{-m},d)\leq N$.

To simplify the combinatorial part later on, we multiply all coordinates by $2^{m}$, which results to vectors with integer components. This motivates the following definition.

Definition 3.2**.**

Let $m\geq 2$. We say that $T=\{x^{1},\dots,x^{N}\}\subset\{1,2,\dots,2^{m}-1\}^{d}$ satisfies the condition (S) of the order $m$ if for every ${\mathcal{A}}\subset\{1,\dots,d\}$ with $\#{\mathcal{A}}=A_{m}:=\min(m2^{m},d)$, the set of restrictions $x^{1}|_{\mathcal{A}},\dots,x^{N}|_{\mathcal{A}}$ contains all $(2^{m}-1)^{A_{m}}$ possible values.

By what we said above, anytime $m\geq 2$ and $T=\{x^{1},\dots,x^{N}\}$ satisfies the condition (S) of the order $m$, then $P\subset{\mathbb{M}}_{m}^{d}$ with $P=2^{-m}\cdot T$ satisfies $\operatorname{disp}(P)\leq 2^{-m}$. The proof of Theorem 3.1 is therefore finished, once we find a set $T$ with $\#T\leq c_{m}\log(d)$ satisfying the condition (S) of order $m$.

3.2. Sosnovec’s proof - randomized construction

The rest of the proof in [35] can now be understood as a randomized construction of a small set satisfying the condition (S). Indeed, there are $\binom{d}{A_{m}}$ subsets $S$ of $\{1,\dots,d\}$ with $A_{m}$ elements. We now fix one such set $S$ and one vector $z\in\{1,\dots,2^{m}-1\}^{A_{m}}$. The probability that a point $x\in\{1,\dots,2^{m}-1\}^{d}$ chosen at random (from the uniform distribution) fulfills $x|_{S}=z$ is $(2^{m}-1)^{-A_{m}}$. Therefore, the probability that none of the $N$ randomly chosen points $x^{1},\dots,x^{N}$ fulfills this restriction is $[1-(2^{m}-1)^{-A_{m}}]^{N}$. Finally, the probability that there is a set $S\subset\{1,\dots,d\}$ with $\#S=A_{m}$ and a vector $z\in\{1,\dots,2^{m}-1\}^{A_{m}}$ such that no $\{x^{1},\dots,x^{N}\}$ satisfies $x^{j}|_{S}=z$ is, by the union bound, at most

[TABLE]

By simple calculus, if $N$ is of the order $c_{m}\log_{2}(d)$ for $c_{m}$ large enough, this expression is smaller than one. Hence, with positive probability, the randomly chosen point set $\{x^{1},\dots,x^{N}\}$ satisfies the condition (S).

3.3. Derandomization using universal sets

Definition 3.2 resembles very much the concept of $(n,k)$-universal sets, see Section 2.1. In particular, it is easy to see that every $(d,A_{m},2^{m}-1)$-universal set -after adding 1 to each coordinate- satisfies condition (S) of the order $m$.

Therefore, we can use Theorem 2.5 to replace the random arguments of the last section by a deterministic algorithm. We glue all the components together in the form of an algorithm.

Algorithm 1

(1)

For $\varepsilon\in(0,\frac{1}{4}]$ and $d\geq 2$, choose a positive integer $m$

with $2^{-m}\leq\varepsilon<2^{-m+1}$ and set $A_{m}:=\min(m2^{m},d)$;

(2)

Generate an $(md,mA_{m})$-universal set as in [26, Theorem 6];

(3)

Interpret these vectors as digital decompositions to obtain

an $(d,A_{m},2^{m})$-universal set;

(4)

Replace $2^{m}-1$ by [math] in all coordinates to obtain

an $(d,A_{m},2^{m}-1)$-universal set;

(5)

Increase all the coordinates by one; (This set satisfies (S).)

(6)

Finally, divide all the coordinates by $2^{m}$ and output the point set.

It remains to consider the running-time of this algorithm.

Theorem 3.3**.**

Let $\varepsilon\in(0,\frac{1}{4}]$ and $d\geq 2$. Then there is a positive constant $c_{\varepsilon}>0$, such that Algorithm 1 constructs a set $P\subset[0,1]^{d}$ with $\operatorname{disp}(P)\leq\varepsilon$ and

[TABLE]

*The running time of Algorithm 1 is linear in the length of the output. *

Proof.

The first part of the theorem is proven by what we said above. Concerning the running-time, we obtain from Theorem 2.5, that the $(d,A_{m},2^{m}-1)$-universal set, which we generate in steps (2)–(4) of Algorithm 1, can be constructed in linear time and with size

[TABLE]

The remaining operations can be done in a linear time, without enlarging the point set.

∎

As expected, the dependence of the size of $P$ on $2^{m}\approx\varepsilon^{-1}$ is rather bad (as it was in [35]), but there is indeed only a logarithmic dependence on $d$.

4. Improving the dependence in $\varepsilon^{-1}$

The main aim of [39] was to refine the analysis of [35] and to achieve a better dependence of $N(\varepsilon,d)$ on $\varepsilon$, without sacrificing the logarithmic dependence on $d$. The main theorem of [39] was the following.

Theorem 4.1**.**

Let $d\geq 2$ be a natural number and let $\varepsilon\in(0,1/2)$. Then there exists a point set $P\subset[0,1]^{d}$ with $\operatorname{disp}(P)\leq\varepsilon$ and

[TABLE]

Also this result can be derandomized using results from coding theory. By doing this, we will lose some power of $1/\varepsilon$ in the size of the point set. However, it will still be of order $\log(d)$.

4.1. Enhanced analysis of the random construction

The main novelty of [39] was a more careful splitting of $\Omega_{m}$ (see (3.1)) into subgroups. To be more specific (and using the notation of [19]), for $s=(s_{1},\dots,s_{d})\in\{1,\dots,2^{m}-1\}^{d}$ and $p=(p_{1},\dots,p_{d})\in\mathbb{M}_{m}^{d}$, we denoted $\Omega_{m}(s,p)$ to be those cubes from $\Omega_{m}$, which have $I_{\ell}$ approximately of the length $\frac{s_{\ell}}{2^{m}}$ and its left point around $p_{\ell}$ for all $\ell=1,\dots,d$, i.e.

[TABLE]

We denote by ${\mathbb{I}}_{m}$ the pairs $(s,p)$, for which $\Omega_{m}(s,p)$ is non-empty. It is easy to see that their number is bounded from above by

[TABLE]

see [19, eq. (2.1)]. If $(s,p)\in{\mathbb{I}}_{m}$, we further set

[TABLE]

The advantage of dividing of $\Omega_{m}$ into the groups $\Omega_{m}(s,p)$ is the surprisingly good control of the probability that a randomly chosen point $z\in{\mathbb{M}}_{m}^{d}$ lies in the intersection of all the cubes from $\Omega_{m}(s,p).$ It is actually, up to a constant, of the same order as the volume of each of the cubes in $\Omega_{m}$, i.e. of $2^{-m}$.

Lemma 4.2**.**

([39, Lemma 3] and [19, Lemma 2.1]) Let $m\in{\mathbb{N}}$, $s\in\{0,1,\dots,2^{m}-1\}^{d}$ and $p\in\{1/2^{m},\dots,1-1/2^{m}\}^{d}$ be such that $(s,p)\in{\mathbb{I}}_{m}$. Let $z$ be uniformly distributed in $\mathbb{M}_{m}^{d}$. Then

[TABLE]

The aim of [39] was to combine (4.1) with Lemma 4.2 and the union bound. Indeed, the probability that a randomly chosen point $z\in{\mathbb{M}}_{m}^{d}$ avoids ${\mathbb{B}}_{m}(s,p)$ is at most $1-2^{-m-4}.$ Therefore, the probability that a set $P=\{x^{1},\dots,x^{N}\}\subset{\mathbb{M}}_{m}^{d}$ of $N$ randomly and independently generated points does not intersect ${\mathbb{B}}_{m}(s,p)$ is at most $(1-2^{-m-4})^{N}$, i.e.

[TABLE]

By the union bound over all $(s,p)\in{\mathbb{I}}_{m}$, we get further

[TABLE]

As ${\mathbb{B}}_{m}(s,p)$ was defined in (4.2) as the intersection of all cubes from $\Omega_{m}(s,p)$, finding a point $x^{\ell}\in{\mathbb{B}}_{m}(s,p)$ means that the same point may be found in all cubes in $\Omega_{m}(s,p).$ We conclude that if $N$ is large enough to ensure that

[TABLE]

i.e., $N>m2^{2m+4}\log(2^{m+3}d)$, then the randomly generated $P=\{x^{1},\dots,x^{N}\}\subset{\mathbb{M}}_{m}^{d}$ intersects every ${\mathbb{B}}\in\Omega_{m}$ with positive probability. Hence, there exists $P$ with $\#P\leq N$ such that $\operatorname{disp}(P)\leq 2^{-m}$. This is essentially the result of [39].

4.2. Connection to $k$-restriction problems

By what we said above, if a point set $P\subset[0,1]^{d}$ intersects ${\mathbb{B}}_{m}(s,p)$ for all $(s,p)\in{\mathbb{I}}_{m}$, then $\operatorname{disp}(P)\leq 2^{-m}.$ The randomized construction in [39], which we now want to replace by a deterministic one, was restricted in its choice of points to ${\mathbb{M}}_{m}^{d}$. Therefore we define

[TABLE]

Definition 4.3**.**

Let $m\geq 2$. We say that $T\subset\{1,2,\dots,2^{m}-1\}^{d}$ satisfies the condition (S’) of the order $m$ if, for every $(s,p)\in{\mathbb{I}}_{m}$ it intersects ${\mathbb{C}}_{m}(s,p).$

The rest of [39] then provides a randomized construction of a small set $T$, which satisfies the condition $(S^{\prime})$ of order $m$. This task has two things in common with $k$-restriction problems. First, the system $\{{\mathbb{C}}_{m}(s,p):(s,p)\in{\mathbb{I}}_{m}\}$ is invariant under permutations of $\{1,\dots,d\}$ and second, the number of active coordinates is, for every $(s,p)\in{\mathbb{I}}_{m}$, bounded from above by $A_{m}:=\min(m2^{m},d).$

To build the connection between the condition $(S^{\prime})$ and the $k$-restriction problems, we choose the quadruplet of parameters $(b,k,n,m)$, see Definition 2.6, as $(2^{m}-1,A_{m},d,M)$. The system ${\mathcal{C}}$ collects the sets ${\mathbb{C}}_{m}(s,p)$ for those $(s,p)\in{\mathbb{I}}_{m}$ which have the corresponding active coordinates in $\{1,\dots,A_{m}\}$. Finally, $M$ is the cardinality of ${\mathcal{C}}.$

More formally, let $(s,p)\in{\mathbb{I}}_{m}$ with $s_{j}=2^{m}-1$ for $j>A_{m}$. Then, we set

[TABLE]

and define

[TABLE]

We observe that a set $T$ satisfies the condition $(S^{\prime})$ of order $m$ if, and only if, the set $T-1$ satisfies the $k$-restriction problem with respect to ${\mathcal{C}}.$

The parameter $M$, which is just the cardinality of ${\mathcal{C}}$, can be estimated from above in a way similar to (4.1), but note that we do not have to choose the subset of active indices anymore. Each $C_{m}{(s,p)}\in{\mathcal{C}}$ is characterized by $s\in\{0,1,\dots,2^{m}-1\}^{A_{m}}$ and $p\in\{1/2^{m},\dots,(2^{m}-1)/2^{m}\}^{A_{m}}$. Therefore,

[TABLE]

The second important parameter of a $k$-restriction problem is the minimal size $c:=c(\mathcal{C})$ of each of the restriction sets $C_{m}(s,p)$, see (2.1). A lower bound on $c$ follows directly from Lemma 4.2 and we obtain

[TABLE]

With the choice of parameters as given above, we have $c/b^{k}\geq 2^{-m-4}$.

4.3. A first attempt for derandomization

Using the arguments of the last section, one could use the construction from Theorem 2.7 directly to solve the corresponding $k$-restriction problem with parameters $(2^{m}-1,A_{m},d,2^{2mA_{m}})$, whenever $d>2^{m}$. This leads to a point set $P\subset\mathbb{M}_{m}^{d}$ with $\operatorname{disp}(P)\leq 2^{-m}$ and

[TABLE]

Note that this bound matches the union bound from Section 4.1. However, the running-time of the algorithm, as given by Theorem 2.7, is

[TABLE]

where $\mathcal{T}$ is the time complexity of the membership oracle, which can be assumed of the order $\mathcal{O}(A_{m})=\mathcal{O}(m2^{m})$ in this case.

4.4. Derandomization using splitters

We now describe how we can improve the construction of an explicit solution to the desired $k$-restriction problem with parameters $(b,k,n,m)$ equal to $(2^{m}-1,A_{m},d,2^{2mA_{m}})$ and the set system ${\mathcal{C}}$ defined by (4.3). We use the approach of [26] to obtain solutions of the $k$-restriction problem which are ’small’ in size and running-time of the corresponding algorithm. ’Small’ means here, that the dependence on the original problem dimension $d$ is as small as possible.

In the heart of the constructions are splitters, see Section 2.3. As already indicated in Section 2.3, we use a $(d,A_{m},A_{m}^{2})$-splitter, say $A(m,d)$, to map the original $d$-dimensional problem to a $k$-restriction problem in dimension $A_{m}^{2}$, which can then be solved with cost independent of $d$.

Recall that, by Definition 2.8, $A(m,d)$ is an $(d,A_{m},A_{m}^{2})$-splitter, if it is a collection of mappings $a:\{1,\dots,d\}\to\{1,\dots,A_{m}^{2}\}$ such that for every $S\subset\{1,\dots,d\}$ with $\#S=A_{m}$, there is an $a\in A(m,d)$, which is injective on $S$, i.e., $a(S)\subset\{1,\dots,A_{m}^{2}\}$ has $A_{m}$ elements.

Further, let $T(m)\subset\{1,2,\dots,2^{m}-1\}^{A_{m}^{2}}$ be the solution of the $k$-restriction problem with parameters $(2^{m}-1,A_{m},A_{m}^{2},2^{2mA_{m}})$ with respect to the original system of restrictions ${\mathcal{C}}$, see (4.3). This means, that $T(m)=\{\tau_{1},\dots,\tau_{K}\}\subset\{0,1,\dots,2^{m}-2\}^{A_{m}^{2}}$ such that for every $S^{\prime}\subset\{1,\dots,A_{m}^{2}\}$ with $\#S^{\prime}=A_{m}$ and any $C\in\mathcal{C}$ there is $\tau\in T(m)$ with $\tau|_{S^{\prime}}\in C$.

Now we are in the position to define the solution to the restriction problem with parameters $(2^{m}-1,A_{m},d,2^{2mA_{m}})$ and the system ${\mathcal{C}}$. Indeed, we define

[TABLE]

to be the set of concatenations of any splitter with any element of the solution to the restriction problem. Here, we switch between the notion of vectors of length $d$ (resp. $A_{m}^{2}$) and mappings from $\{1,\dots,d\}$ (resp. $\{1,\dots,A_{m}^{2}\}$) to $\{0,1,\dots,2^{m}-2\}$. This should not lead to any confusion.

To show that $T^{*}$ is indeed a solution to our $k$-restriction problem, let $S\subset\{1,\dots,d\}$ with $\#S=A_{m}$. Then, there exists $a\in A(m,d)$, such that $S^{\prime}=a(S)\subset\{1,2,\dots,A_{m}^{2}\}$ has $A_{m}$ mutually different elements, i.e., $\#S^{\prime}=A_{m}$. Now, for every $C\in\mathcal{C}$, there is some $\tau\in T(m)$, such that $C\ni\tau|_{S^{\prime}}=(\tau\circ a)|_{S}$. Hence, $T^{*}$, which satisfies $\#T^{*}=\#T(m)\cdot\#A(m,d)$, is a solution to the restriction problem with parameters $(2^{m}-1,A_{m},d,2^{2mA_{m}})$.

We merge all the components together in a form of an algorithm.

Algorithm 2

(1)

For $\varepsilon\in(0,\frac{1}{4}]$ and $d\geq 2$, choose a positive integer $m$

with $2^{-m}\leq\varepsilon<2^{-m+1}$ and set $A_{m}:=\min(m2^{n},d)$;

(2)

Generate a $(d,A_{m},A_{m}^{2})$-splitter $A(m,d)$ as in Lemma 2.9;

(3)

Generate a solution $T(m)$ to the $k$-restriction problem with parameters $(2^{m}-1,A_{m},A_{m}^{2},2^{2mA_{m}})$ and restrictions $\mathcal{C}$ from (4.3) as in Theorem 2.7;

(4)

Set $T^{*}=T(m)\circ A(m,d)\subset\{0,\dots,2^{m}-2\}^{d}$;

(5)

Increase all the coordinates by one, then divide them by $2^{m}$;

(6)

Output the resulting point set $P$.

Theorem 4.4**.**

Let $\varepsilon\in(0,\frac{1}{4}]$ and $d\geq 2$. Then, there is an absolute constant $C<\infty$, such that Algorithm 2 constructs a set $P\subset[0,1]^{d}$ with $\operatorname{disp}(P)\leq\varepsilon$ and

[TABLE]

where $d^{*}:=\max(d,2/\varepsilon)$.

*The running-time of Algorithm 2 is $\mathcal{O}_{\varepsilon}(d^{c})$ for some $c<\infty$. *

Proof.

In case that $d<\lfloor 2/\varepsilon\rfloor$, replace $d$ by $\lfloor 2/\varepsilon\rfloor$ in all what follows and delete at the end the last $\lfloor 2/\varepsilon\rfloor-d$ coordinates of the constructed point set.

By Lemma 2.9, there is an explicit, deterministic construction of a $(d,A_{m},A_{m}^{2})$-splitter $A(m,d)$ of size $\mathcal{O}(A_{m}^{4}\log(d))=\mathcal{O}(m^{4}2^{4m}\log(d))$. Furthermore, by Theorem 2.7 and since $d\geq\lfloor 2/\varepsilon\rfloor>2^{m}-1$, there is a deterministic solution $T(m)$ to the $k$-restriction problem with parameters $(2^{m}-1,A_{m},A_{m}^{2},2^{2mA_{m}})$ and restrictions $\mathcal{C}$ from (4.3) of size

[TABLE]

Here, we also used (4.5). Altogether, we get

[TABLE]

which implies the result since $m\leq 1+\log_{2}(1/\varepsilon)$.

For the running-time, we need $\mathcal{O}_{m}(d^{c})$ for the construction of $A(m,d)$, see Lemma 2.9, and

[TABLE]

for the construction of $T(m)$, see Theorem 2.7, where, again, the cost $\mathcal{T}$ of the membership oracle can be assumed to be $\mathcal{O}(A_{m})$. Finally, we need $\mathcal{O}(\#T^{*})$ time to build up $T^{*}$. The remaining steps are less expensive, which implies the result.

∎