The determinant bound for discrepancy is almost tight

Jiri Matousek

Introduction

Let V=[n]:={1,2,…,n}V=[n]:=\{1,2,\ldots,n\} be a vertex set and F={F1,F2,…,Fm}\mathcal{F}=\{F_{1},F_{2},\ldots,F_{m}\} be a system of subsets of VV. The discrepancy of F\mathcal{F} is disc(F):=min⁡χdisc(F,χ){\mathop{\rm disc}\nolimits}(\mathcal{F}):=\min_{\chi}{\mathop{\rm disc}\nolimits}(\mathcal{F},\chi), where the minimum is over all colorings χ ⁣:V→{−1,+1}\chi\colon V\to\{-1,+1\}, and {\mathop{\rm disc}\nolimits}(\mathcal{F},\chi):=\max_{i=1,2,\ldots,m}\bigl{|}\sum_{j\in F_{i}}\chi(j)\bigr{|}.

The hereditary discrepancy of F\mathcal{F} is

Here F∣J\mathcal{F}|_{J} denotes the restriction of the set system F\mathcal{F} to the ground set JJ, i.e., {F∩J:F∈F}\{F\cap J:F\in\mathcal{F}\}.

Bounding the discrepancy or the hereditary discrepancy of a particular set system from below is usually challenging. One of the strongest known tools is a result known as the determinant lower bound. To formulate it we define, for a real matrix AA,

where the maximum is over all k×kk\times k submatrices BB of AA. For a set system F\mathcal{F}, we put detlb(F):=detlb(A){\mathop{\rm detlb}\nolimits}(\mathcal{F}):={\mathop{\rm detlb}\nolimits}(A), where AA is the incidence matrix of F\mathcal{F}.

For every (finite) set system F\mathcal{F} we haveThe bound in [LSV86] is stated without the 12\frac{1}{2} factor. This has two causes: first, their discrepancy is scaled by 12\frac{1}{2} compared to ours, and second, in their argument, at one step they seem to be multiplying by 22 where, in my opinion, one should divide by 2. herdisc(F)≥12detlb(F).{\mathop{\rm herdisc}\nolimits}(\mathcal{F})\geq{\textstyle\frac{1}{2}}{\mathop{\rm detlb}\nolimits}(\mathcal{F}).

Lovász et al. [LSV86] conjectured that the determinant lower bound is tight up to a constant factor, i.e., herdisc(F)=O(detlb(F)){\mathop{\rm herdisc}\nolimits}(\mathcal{F})=O({\mathop{\rm detlb}\nolimits}(\mathcal{F})) for all F\mathcal{F}. This was refuted by an example of Hoffmann,The vertex set in the example is the set of edges of the complete kk-ary tree TT of depth kk (so n≈kkn\approx k^{k}). The set system F\mathcal{F} is a union F1∪F2\mathcal{F}_{1}\cup\mathcal{F}_{2}, where F1\mathcal{F}_{1} consists of the edge sets of all root-to-leaf paths, and F2\mathcal{F}_{2} contains, for each non-leaf vertex vv, the set of the kk edges connecting vv to its successors. We have herdisc(F1)≤1{\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1})\leq 1, herdisc(F2)≤1{\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{2})\leq 1, detlb(F)=O(1){\mathop{\rm detlb}\nolimits}(\mathcal{F})=O(1), and disc(F)=k≈(log⁡n)/(log⁡log⁡n){\mathop{\rm disc}\nolimits}(\mathcal{F})=k\approx(\log n)/(\log\log n). See, e.g., [BS95] or [Mat10] for more details. which shows that herdisc(F)/detlb(F){\mathop{\rm herdisc}\nolimits}(\mathcal{F})/{\mathop{\rm detlb}\nolimits}(\mathcal{F}) can be of order (log⁡n)/(log⁡log⁡n)(\log n)/(\log\log n). A construction of Pálvölgyi [Pál10], also presented in Section 5 below, provides the slightly stronger lower bound of Ω(log⁡n)\Omega(\log n) for the same quantity.

Here we prove that herdisc(F)/detlb(F){\mathop{\rm herdisc}\nolimits}(\mathcal{F})/{\mathop{\rm detlb}\nolimits}(\mathcal{F}) cannot be much larger than in these examples, at least if ∣F∣|\mathcal{F}| is bounded by a polynomial function of nn.

Next, we consider the situation where a set system F\mathcal{F} as above is a union of set systems F1,F2,…,Ft\mathcal{F}_{1},\mathcal{F}_{2},\ldots,\mathcal{F}_{t}, and we are interested in bounding herdisc(F){\mathop{\rm herdisc}\nolimits}(\mathcal{F}) in terms of the hereditary discrepancies of the Fi\mathcal{F}_{i}.

For t=2t=2, this problem was raised by Sós (it is cited, e.g., in [LSV86]). She asked whether herdisc(F1∪F2){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}\cup\mathcal{F}_{2}) can be estimated in terms of herdisc(F1){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}) and herdisc(F2){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{2}) for any two set systems F1\mathcal{F}_{1} and F2\mathcal{F}_{2} (on the same vertex set). Hoffman’s example mentioned above shows that herdisc(F1∪F2){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}\cup\mathcal{F}_{2}) cannot be bounded by a function of herdisc(F1){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}) and herdisc(F2){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{2}) alone.

The next theorem shows that a good bound is possible if we also allow for a moderate dependence on mm and nn. Namely, herdisc(F1∪F2){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}\cup\mathcal{F}_{2}) can exceed max⁡(herdisc(F1),herdisc(F2))\max({\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}),{\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{2})) at most by a factor polylogarithmic in nn and mm, not much more than in Hoffman’s example.A connection of Sós’s question to the tightness of the determinant lower bound was observed in [LSV86], although with a different proof, which would yield a quantitatively weaker result in our setting.

The only previous result in this direction, from [KMV05], shows that if F2\mathcal{F}_{2} consists of a single set, then herdisc(F1∪F2)=O(herdisc(F1)log⁡n){\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}\cup\mathcal{F}_{2})=O({\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1})\log n).

Let F\mathcal{F} be a system of mm sets on nn vertices, and let F=F1∪F2∪⋯∪Ft\mathcal{F}=\mathcal{F}_{1}\cup\mathcal{F}_{2}\cup\cdots\cup\mathcal{F}_{t}. Then

There exist systems F\mathcal{F} of nn sets on nn points with discrepancy of order n\sqrt{n} (e.g., systems derived from Hadamard matrices or random set systems; see, e.g., [Mat10]). If we let Fi\mathcal{F}_{i} be the set system consisting of the iith set of such an nn, i=1,2,…,ni=1,2,\ldots,n, then herdisc(Fi)=1{\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{i})=1, while herdisc(F)=Ω(n){\mathop{\rm herdisc}\nolimits}(\mathcal{F})=\Omega(\sqrt{n}). In this sense, the bound in Theorem 3 is tight up to a polylogarithmic factor, including the dependence on tt.

Theorem 3 is an immediate consequence of Theorems 1 and 2 and of the next linear-algebraic lemma.

Let A1,…,AtA_{1},\ldots,A_{t} be real matrices, each with nn columns, let D:=max⁡i=1,2,…,tdetlb(Ai)D:=\max_{i=1,2,\ldots,t}{\mathop{\rm detlb}\nolimits}(A_{i}), and let AA be a matrix in which each row is a row of some of the AiA_{i}. Then

In case where some of the detlb(Ai){\mathop{\rm detlb}\nolimits}(A_{i}) are much smaller than the others, for example, one can obtain a somewhat better upper bound by making finer estimates in the calculation in the proof. However, a general formulation of such a finer bound looks cumbersome, and it seems that in such cases, a similar improvement can usually be achieved by applying the lemma repeatedly in several stages with various values of DD.

The bound in Lemma 4 is generally tight up to a constant factor; this can be seen from an example similar to the one mentioned below Theorem 3. Namely, let AA be an n×nn\times n Hadamard matrix (i.e., a matrix with pairwise orthogonal rows and with ±1\pm 1 entries, which is well known to exist for infinitely many values of nn), and let AiA_{i} be the single-row matrix made of the iith row of AA, i=1,2,…,ni=1,2,\ldots,n. Then, obviously, detlb(Ai)=1{\mathop{\rm detlb}\nolimits}(A_{i})=1, while detlb(A)≥det⁡(A)1/n=n{\mathop{\rm detlb}\nolimits}(A)\geq\det(A)^{1/n}=\sqrt{n}. Moreover, if we partition the rows of this AA into tt blocks A1,…,AtA_{1},\ldots,A_{t} by n/tn/t rows each, then the Hadamard bound, stating that the determinant of a matrix is at most the product of the Euclidean norms of the rows, implies detlb(Ai)≤n/t{\mathop{\rm detlb}\nolimits}(A_{i})\leq\sqrt{n/t}. This shows the tightness of Lemma 4 for all t≤nt\leq n.

Vector discrepancy and Bansal’s algorithm. The proof of Theorem 2 is based on a recent breakthrough—an algorithm of Bansal [Ban10]. The algorithm produces a low-discrepancy coloring of a given set system, using semidefinite programming and a clever randomized rounding strategy. To state the consequence of Bansal’s work that we will use, we first introduce another notion of discrepancy.

where ∥⋅∥\|\cdot\| is the Euclidean norm. So, for vector discrepancy, one colors by unit vectors instead of ±1\pm 1’s. Since a ±1\pm 1 coloring can also be regarded as a vector coloring by the vectors e1=(1,0,…,0){\bf e}_{1}=(1,0,\ldots,0) and −e1-{\bf e}_{1}, we have vecdisc(F)≤disc(F){\mathop{\rm vecdisc}\nolimits}(\mathcal{F})\leq{\mathop{\rm disc}\nolimits}(\mathcal{F}).

The hereditary vector discrepancy hervecdisc(F){\mathop{\rm hervecdisc}\nolimits}(\mathcal{F}) is the maximum vector discrepancy of a restriction of F\mathcal{F} to a subset J⊆VJ\subseteq V.

We conjecture that the claim of Theorem 5 actually holds with log⁡(mn)\sqrt{\log(mn)} instead of log⁡(mn)\log(mn). If true, this would yield a similar improvement in Theorem 2 and get close to an asymptotically optimal bound, at least assuming mm bounded by a polynomial in nn.

A dual formulation of the vector discrepancy

Proof. We will use the duality of semidefinite programming. Dualizing a semidefinite program is a routine procedure, but unfortunately, I am not aware of an explicit recipe for the general case in the literature. Rather than converting the relevant semidefinite program to a standard form, it seems more convenient to use a duality theorem for conic programing from Duffin [Duf56], which we now introduce.

Let V,WV,W be real vector spaces (for our purposes we may assume that they are finite-dimensional), each with a scalar product, which we denote by ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle in both cases. Let K⊆VK\subseteq V and L⊆WL\subseteq W be closed convex cones,A convex cone is a convex set KK such that x∈K{\bf x}\in K implies λx∈K\lambda{\bf x}\in K for all λ≥0\lambda\geq 0. The dual cone of KK is K∗={y:⟨x,y⟩≥0\mboxforallx∈K}K^{*}=\{{\bf y}:\langle{\bf x},{\bf y}\rangle\geq 0\mbox{ for all }{\bf x}\in K\}. A simple property we will often use is (K⊕L)∗=K∗⊕L∗(K\oplus L)^{*}=K^{*}\oplus L^{*}, where ⊕\oplus denotes direct sum; we assume K⊆VK\subseteq V, L⊆WL\subseteq W, where V,WV,W are disjoint vector spaces, and K⊕L⊆V⊕WK\oplus L\subseteq V\oplus W. let b∈W{\bf b}\in W and c∈V{\bf c}\in V be vectors, and let F ⁣:V→WF\colon V\to W be a linear map. We consider the primal cone program (P)

Here FT ⁣:W→VF^{T}\colon W\to V denotes the adjoint of FF; if we fix orthonormal bases in VV and WW, then the matrix of FTF^{T} is the transpose of the matrix of FF. The duality theorem asserts that if the maximum in (P) is a finite number γ\gamma and if the set of feasible solutions of (P) has an interior point, then (D) is feasible as well and its minimum equals γ\gamma. (Since (P) and (D) are dual to one another, one can also interchange their role in the theorem.)

In our case, we start with a vector program defining vecdisc(F){\mathop{\rm vecdisc}\nolimits}(\mathcal{F}), namely,

For the constraint ∥∑j∈Fiuj∥2≤t\|\sum_{j\in F_{i}}{\bf u}_{j}\|^{2}\leq t, we expand the left-hand side to ∑j,k∈FiujTuk\sum_{j,k\in F_{i}}{\bf u}_{j}^{T}{\bf u}_{k}, and then it translates to

where ai{\bf a}_{i} is the iith row of the incidence matrix AA of F\mathcal{F}, regarded as an n×1n\times 1 matrix (so aiaiT{\bf a}_{i}{\bf a}_{i}^{T} is an n×nn\times n matrix), and ∙\bullet denotes the scalar product of matrices (given by X∙Y=∑i,jxijyijX\bullet Y=\sum_{i,j}x_{ij}y_{ij}). The constraint ∥uj∥=1\|{\bf u}_{j}\|=1 then reads

where EjE_{j} is the matrix with 11 at position (j,j)(j,j) and s elsewhere.

It remains to verify that the duality theorem can be applied to these (P) and (D). We will check that (D) is bounded and has a feasible interior point. To verify boundedness, which means that ∑j=1nzj\sum_{j=1}^{n}z_{j} cannot be arbitrarily large, we use (1) with x=1n{\bf x}={\bf 1}_{n}: then the left-hand side is bounded, and so ∑j=1nzj\sum_{j=1}^{n}z_{j} is bounded as well.

For an interior feasible point, we can take, e.g., w1=⋯=wm=12mw_{1}=\cdots=w_{m}=\frac{1}{2m}, z1=⋯=zn=−1z_{1}=\cdots=z_{n}=-1. Then (1) obviously holds for all x{\bf x}.

Thus, the duality theorem applies and shows that if vecdisc(F)≥D{\mathop{\rm vecdisc}\nolimits}(\mathcal{F})\geq D, then the maximum in (4) is at least D2D^{2}. This yields the existence of the desired wiw_{i} and zjz_{j}, and the lemma is proved. □\Box

Proof of Theorem 2

We begin with a simple and probably standard lemma.

where y[K]{\bf y}[K] denotes the ∣K∣|K|-component vector (yj:j∈K)(y_{j}:j\in K).

Proof. Let ymax:=max⁡j∣yj∣y_{\rm max}:=\max_{j}|y_{j}|, and for i=0,1,2,…i=0,1,2,\ldots, let Ki:={j:∣yj∣∈(2−i−1ymax,2−iymax]}K_{i}:=\{j:|y_{j}|\in(2^{-i-1}y_{\rm max},2^{-i}y_{\rm max}]\}. The contribution to ∥y∥\|{\bf y}\| of the components of y{\bf y} with indices in KiK_{i} for i≥2log⁡ni\geq 2\log n, say, is negligible, and so there exists some i0i_{0} for which ∑j∈Ki0yj2=Ω(∥y∥2/log⁡n)\sum_{j\in K_{i_{0}}}y_{j}^{2}=\Omega(\|{\bf y}\|^{2}/\log n). Then K:=Ki0K:=K_{i_{0}} will do. □\Box

Theorem 2 will follow from Bansal’s result (Theorem 5) and the next lemma.

Let F={F1,…,Fm}\mathcal{F}=\{F_{1},\ldots,F_{m}\} be a set system on [n][n] with vecdisc(F)=D{\mathop{\rm vecdisc}\nolimits}(\mathcal{F})=D. Then detlb(F)=Ω(D/log⁡n ){\mathop{\rm detlb}\nolimits}(\mathcal{F})=\Omega(D/\sqrt{\log n}\,).

Proof. We begin with the dual formulation of vector discrepancy from Lemma 6. For more convenient notation, we will write the nonnegative weight wiw_{i} as βi2\beta_{i}^{2}. Moreover, we let J⊆[n]J\subseteq[n] consist of the indices jj with zj>0z_{j}>0, and we will use the inequality (1) in Lemma 6 only for vectors x{\bf x} that are zero outside JJ. Writing zj=γj2z_{j}=\gamma_{j}^{2} for j∈Jj\in J, we arrive at the inequality

Let C:=A[∗,K]C:=A[*,K] be the m×km\times k incidence matrix of the system F∣K\mathcal{F}|_{K} (consisting of the columns of AA whose indices belong to KK), and let Cˇ\check{C} be the m×km\times k matrix obtained from CC by multiplying the iith row by βi\beta_{i}. Then (6) can be rewritten as

where the summation is over all kk-element subsets I⊆[m]I\subseteq[m] and Cˇ[I,∗]\check{C}[I,*] consists of the rows of Cˇ\check{C} whose indices lie in II.

We have det⁡(Cˇ[I,∗])=det⁡(C[I,∗])∏i∈Iβi\det(\check{C}[I,*])=\det(C[I,*])\prod_{i\in I}\beta_{i}. Setting M:=max⁡I∣det⁡(C[I,∗])∣M:=\max_{I}|\det(C[I,*])|, we can rewrite the right-hand side of (7) and estimate it as follows:

where the penultimate inequality follows because every term ∏i∈Iβi2\prod_{i\in I}\beta_{i}^{2} occurs k!k! times in the multinomial expansion of (β12+⋯+βm2)k(\beta_{1}^{2}+\cdots+\beta_{m}^{2})^{k}.

Proof of Theorem 2. By Theorem 5, there is a subset J⊆[n]J\subseteq[n] with vecdisc(F∣J)=Ω(herdisc(F)/log⁡(mn)){\mathop{\rm vecdisc}\nolimits}(\mathcal{F}|_{J})=\Omega({\mathop{\rm herdisc}\nolimits}(\mathcal{F})/\log(mn)). Theorem 2 follows by applying Lemma 8 to F∣J\mathcal{F}|_{J}. □\Box

Proof of Lemma 4

On the other hand, by the Binet–Cauchy formula, we have

Putting (8), (9), and (10) together, we arrive at

Then we estimate, using the concavity of the function x↦xln⁡(1/x)x\mapsto x\ln(1/x) and Jensen’s inequality,

Thus det⁡(B)1/k≤Det\det(B)^{1/k}\leq D\sqrt{et} as claimed. □\Box

Remark. The case t=2t=2 has a somewhat simpler proof using the Laplace expansion of det⁡B\det B, which asserts that

where the sum is over all ∣I∣|I|-element subsets J⊆[k]J\subseteq[k], I‾=[k]∖I\overline{I}=[k]\setminus I, and sgn(I,J)∈{±1}\mathop{\rm sgn}\nolimits(I,J)\in\{\pm 1\} is a sign depending on II and JJ in a way that is of no concern for us (see, e.g., [BJN83, Theorem 4.3]).

A bound on detlb(A){\mathop{\rm detlb}\nolimits}(A) for larger tt can also be obtained by iterating this argument, but this method apparently leads only to detlb(A)=O(tD){\mathop{\rm detlb}\nolimits}(A)=O(tD).

Pálvölgyi’s example

As was pointed out by Pálvölgyi (private communication, 2011), his geometric construction in [Pál10] actually provides a slight quantitative improvement over Hoffman’s example. Translated to the setting of set systems, the construction yields, for every k≥1k\geq 1, two systems F1,F2\mathcal{F}_{1},\mathcal{F}_{2} of kk-element subsets of [n][n], with n=(2kk)−1<4kn={2k\choose k}-1<4^{k}, such that

herdisc(F1),herdisc(F2)≤1{\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{1}),{\mathop{\rm herdisc}\nolimits}(\mathcal{F}_{2})\leq 1, and

under every two-coloring of [n][n], F1∪F2\mathcal{F}_{1}\cup\mathcal{F}_{2} contains a monochromatic set (and consequently, disc(F1∪F2)=k{\mathop{\rm disc}\nolimits}(\mathcal{F}_{1}\cup\mathcal{F}_{2})=k).

Since the construction in [Pál10] is presented geometrically, and property (i) is not entirely obvious, we provide a short self-contained exposition.

I would like to thank Nikhil Bansal for enlightening e-mail discussions concerning his algorithm and for help with dualizing the semidefinite program, and Dömötör Pálvölgyi for explaining me how his construction in [Pál10] improves on Hoffmann’s example.

References