Inequalities and tail bounds for elementary symmetric polynomial with applications

Parikshit Gopalan, Amir Yehudayoff

Introduction

The main message of this paper is that small total variance is sufficient to ensure that the product rule holds approximately even under kk-wise independence.

Let X1,…,XnX_{1},\ldots,X_{n} be random variables each distributed in the range $,withmean, with mean\mu_{i}andvarianceand variance\sigma^{2}_{i}repectively.Letrepectively. Let\sigma^{2}=\sum_{i}\sigma_{i}^{2}.Thereexistconstants. There exist constantsc_{1}>1andand1>c_{2}>0suchthatunderanysuch that under anyk−wiseindependentdistribution-wise independent distribution{\mathcal{D}}$,

Specifically, if σ<1/(2c1)\sigma<1/(2c_{1}) then k=O(log⁡(1/δ)/log⁡(1/σ))k=O(\log(1/\delta)/\log(1/\sigma))-wise independence suffices for Equation (1).

An important restriction that naturally arises is positivity, where each XiX_{i} lies in the interval $$. This setting of parameters (positive variables, small total variance) is important for the applications considered in this paper: pseudorandom generators for combinatorial rectangles [EGL+98, LLSZ97] and min-wise independent permutations [BCFM00]. The former is an important problem in the theory of unconditional pseudorandomness which has been studied intensively [EGL+98, LLSZ97, SSZZ99, ASWZ96, Lu02, GMR+12]. Min-wise independent hashing was introduced by Broder et al. [BCFM00] motivated by similarity estimation, and further studied by [Ind99, BCM98, SSZZ99]. [SSZZ99] showed that PRGs for rectangles give min-wise independent hash functions.

The results of [EGL+98, Ind99] tell us that under kk-wise independence, positivity and boundedness, the LHS of Equation (1) is bounded by exp⁡(−Ω(k))\exp(-\Omega(k)), hence k=O(log⁡(1/δ))k=O(\log(1/\delta)) suffices for error δ\delta. In contrast, we have seen that such a bound cannot hold in the case.However,oncethevarianceissmallerthansomeconstant,ourboundbeatsthisboundeveninthecase. However, once the variance is smaller than some constant, our bound beats this bound even in the setting. Concretely, when σ2<n−ε\sigma^{2}<n^{-\varepsilon} for some ε>0\varepsilon>0, our result says that O(1)O(1)-wise independence suffices for inverse polynomial error in Equation (1), as opposed to O(log⁡(n))O(\log(n))-wise independence. This improvement is crucial in analyzing PRGs and hash functions in the polynomially small error regime. A recent result of [GMR+12] achieves near-logarithmic seed-length for both these problems, even in the regime of inverse polynomial error. Their construction is simple, but its analysis is not. Using our results, we give a modular analysis of the pseudorandom generator construction for rectangles of [GMR+12], using the viewpoint of hash functions.

Our analysis is simpler and perhaps more intuitive. It also improves the seed-length of the construction, getting the dependence on the dimension nn down to O(log⁡log⁡(n))O(\log\log(n)) as opposed to O(log⁡(n))O(\log(n)), which (nearly) matches a lower bound due to [LLSZ97]. Given the basic nature of the question, we feel our results might find other applications. Very recently, [GKM15] constructed the first pseudorandom generators with near-logrithmic seed-length for several classes of functions including halfspaces, modular tests and combinatorial shapes. The key technical ingredient of their work is a generalization of Theorem 1 to the setting where each XiX_{i} takes values in the unit complex disc.

The main technical ingredient in our work is a new analytic inequality about symmetric polynomials in real variables which we believe is independently interesting. The kk’th symmetric polynomial in a=(a1,a2,…,an)a=(a_{1},a_{2},\ldots,a_{n}) is defined as

We give an overview of the new inequality, its use in the derivation of bounds under limited independence, and finally the application of these bounds to the construction of pseudorandom generators and hash functions.

The elementary polynomials appear as coefficients of a univariate polynomial with real roots, since ∏i∈[n](ξ+ai)=∑k=0nξkSn−k(a)\prod_{i\in[n]}(\xi+a_{i})=\sum_{k=0}^{n}\xi^{k}S_{n-k}(a). Symmetric polynomials have been well studied in mathematics, dating back to classical results of Newton and Maclaurin (see [Ste04] for a survey). This work focuses on their growth rates. Specifically, we study how local information on Sk(a)S_{k}(a) for two consecutive values of kk implies global information for all larger values of kk.

It is easy to see that symmetric polynomials over the real numbers have the following property:

Over the real numbers, if S1(b)=S2(b)=0S_{1}(b)=S_{2}(b)=0 then b=0b=0.

This is equivalent to saying that if p(ξ)p(\xi) is a real univariate polynomial of degree nn with nn nonzero roots and p′(0)=p′′(0)=0p^{\prime}(0)=p^{\prime\prime}(0)=0 then p≡0p\equiv 0. This does not hold over all fields, for example, the polynomial p(ξ)=ξ3+1p(\xi)=\xi^{3}+1 has three nonzero complex roots and p′(0)=p′′(0)=0p^{\prime}(0)=p^{\prime\prime}(0)=0.

That is, if S1(a),S2(a)S_{1}(a),S_{2}(a) are small in absolute value, then so is everything that follows. We provide an essentially optimal bound.

The parameters promised by Theorem 2 are tight up to an exponential in kk which is often too small to matter (we do not attempt to optimise the constants). For example, if ai=(−1)ia_{i}=(-1)^{i} for all i∈[n]i\in[n] then ∣S1(a)∣≤1|S_{1}(a)|\leq 1 and ∣S2(a)∣≤n+1|S_{2}(a)|\leq n+1 but Sk(a)S_{k}(a) is roughly (n/k)k/2(n/k)^{k/2}.

A more general statement than Fact A actually holds (see Appendix A for a proof).

We prove a robust version of this fact as well: A twice-in-a-row bound on the increase of the symmetric functions implies a bound on what follows.

Theorem 3 is proved by reduction to Theorem 2. The proof of Theorem 2 is analytic and uses the method of Lagrange multipliers, and is different from that of [GMR+12] which relied on the Newton-Girrard identities. The argument is quite general, and similar bounds may be obtained for functions that are recursively defined. The proof can be found in Section 2.

Stronger bounds are known when the inputs are nonnegative. When ai≥0a_{i}\geq 0 for all i∈[n]i\in[n], the classical Maclaurin inequalities [Ste04] imply that Sk(a)≤(e/k)k(S1(a))kS_{k}(a)\leq(e/k)^{k}(S_{1}(a))^{k}. In contrast, when we do not assume non-negativity, one cannot hope for such bounds to hold under the assumption that ∣S1(a)∣|S_{1}(a)| or any single ∣Sk(a)∣|S_{k}(a)| is small (cf. the alternating signs example above).

2 Expectations of products under limited independence

We return to the question alluded to earlier about how much independence is required for the approximate product rule of expectation. This question arises in the context of min-wise hashing [Ind99], PRGs for combinatorial rectangles [EGL+98, GMR+12], read-once DNFs [GMR+12] and more.

We briefly outline our approach. We start from the results of [EGL+98, Ind99] who give an error bound of exp⁡(−k)\exp(-k). To prove this, they consider random variables Yi=1−XiY_{i}=1-X_{i}, so that

Our approach replaces inclusion-exclusion by a Taylor-series style expansion about the mean, as in [GMR+12]. Let us assume μi≠0\mu_{i}\neq 0 and let Xi=μi(1+Zi)X_{i}=\mu_{i}(1+Z_{i}). Thus,

Let D\cal{D} denote a distribution over Z=(Z1,…,Zn)Z=(Z_{1},\ldots,Z_{n}) as above where the ZiZ_{i}s are (2k+2)(2k+2)-wise independent. For t>0t>0 andA weaker but more technical assumption on t,σ,kt,\sigma,k suffices, see Equation (24). 16etσ≤116et\sigma\leq 1,

3 Applications to pseudorandom generators and hash functions

A hash function is a map h:[n]→[m]h:[n]\rightarrow[m]. Let U\mathcal{U} denote the family of all hash functions h:[n]→[m]h:[n]\rightarrow[m]. Let H⊆U\mathcal{H}\subseteq\mathcal{U} be a family of hash functions. For S⊆[n]S\subseteq[n], let min⁡h(S)=min⁡x∈Sh(x)\min h(S)=\min_{x\in S}h(x). The notion of min-wise independent hashing was introduced by Broder et al. [BCFM00] motivated by similarity estimation, and independently by Mulmuley [Mul96] motivated by computational geometry. The following generalization was introduced by Broder et al. [BCM98]:

Combinatorial rectangles are a well-studied class of tests in pseudorandomness [EGL+98, LLSZ97, SSZZ99, ASWZ96, Lu02, GMR+12]. In addition to being a natural class of statistical tests, constructing generators for them with optimal seeds (up to constant factors) will improve on Nisan’s generator for logspace [ASWZ96], a long-standing open problem in derandomization.

A combinatorial rectangle is a function f:[m]n→{0,1}f:[m]^{n}\rightarrow\{0,1\} which is specified by nn co-ordinate functions fi:[m]→{0,1}f_{i}:[m]\rightarrow\{0,1\} as f(x1,…,xn)=∏i∈mfi(xi)f(x_{1},\ldots,x_{n})=\prod_{i\in m}f_{i}(x_{i}). A map G:{0,1}r→[m]n\mathcal{G}:\{0,1\}^{r}\rightarrow[m]^{n} is a PRG{\mathsf{PRG}} for combinatorial rectangles with error ε\varepsilon if for every combinatorial rectangle f:[m]n→{0,1}f:[m]^{n}\rightarrow\{0,1\},

We take the view of GMR\mathcal{G_{MR}} as a collection of hash functions g:[n]→[m]g:[n]\to[m], based on iterative applications of an alphabet squaring step. We describe the generator formally in Section 5. We start by observing that fooling rectangles is easy when mm is small; O(log⁡(1/δ))O(\log(1/\delta))-wise independnce suffices, and this requires O(log⁡(1/δ)log⁡(m))=O(log⁡(1/δ))O(\log(1/\delta)\log(m))=O(\log(1/\delta)) random bits for m=O(1)m=O(1).

The key insight in [GMR+12] is that gradually increasing the alphabet is also easy (in that it requires only logarithmic randomness). Assume that we have a hash function g0:[n]→[m]g_{0}:[n]\to[m] and from it, we define g1:[n]→[m2]g_{1}:[n]\to[m^{2}]. To do this, we pick a function g1′:[n]×[m]→[m2]g_{1}^{\prime}:[n]\times[m]\to[m^{2}] and set g1(i)=g1′(i,g0(i))g_{1}(i)=g_{1}^{\prime}(i,g_{0}(i)). The key observation is that it suffices to pick g1′g_{1}^{\prime} using only O(log⁡(1/δ)/log⁡(m))O(\log(1/\delta)/\log(m))-wise independence (rather than the O(log⁡(1/δ))O(\log(1/\delta))-wise independence needed for one shot).

Let GMR\mathcal{G_{MR}} be the family of hash functions from [n][n] to [m][m] defined in Section 5.1 with error parameter δ>0\delta>0. The seed length is at most O((log⁡log⁡(n)+log⁡(m/δ))log⁡log⁡(m/δ))O((\log\log(n)+\log(m/\delta))\log\log(m/\delta)). Then, for every S1,…,Sn⊆[m]S_{1},\ldots,S_{n}\subseteq[m],

This improves the [GMR+12] bound in the dependence on nn and δ\delta (their bound was O(log⁡(mn/δ)log⁡log⁡(m)+log⁡(1/δ)log⁡log⁡(1/δ)log⁡log⁡log⁡(1/δ))O(\log(mn/\delta)\log\log(m)+\log(1/\delta)\log\log(1/\delta)\log\log\log(1/\delta))). In particular, the dependence on nn reduces from log⁡(n)\log(n) to log⁡log⁡(n)\log\log(n)The reason log⁡log⁡(n)\log\log(n) seedlength is possible is because every rectangle can be ε\varepsilon-approximated by one that depends only on O(mlog⁡(1/ε))O(m\log(1/\varepsilon)) co-ordinates. Hence the number of functions to fool grows polynomially in nn, rather than exponentially.. [LLSZ97] showed a lower bound of Ω(log⁡(m)+log⁡(1/ε)+log⁡log⁡(n))\Omega(\log(m)+\log(1/\varepsilon)+\log\log(n)) even for hitting sets, so our bound is tight upto the log⁡log⁡(m/δ)\log\log(m/\delta) factor. While [LLSZ97] constructed hitting-set generators for rectangles with near-optimal seedlength, we are unaware of previous constructions of pseudorandom generators for rectangles where the dependence of the seedlength on nn is o(log⁡(n))o(\log(n)).

Combining this with Theorem 19, we get the following corollary.

4 Subsequent work

Very recently, Gopalan, Kane and Meka [GKM15] constructed the first pseudorandom generators with seed-length O((log⁡(n/δ)log⁡log⁡(n/δ)2)O((\log(n/\delta)\log\log(n/\delta)^{2}) for several classes of functions including halfspaces, modular tests and combinatorial shapes. The key technical ingredient of their work is a generalization of Theorem 1 to the setting where the XiX_{i}s are complex valued random variables lying in the unit disc. Their proof however is very different from ours, and in particular it does not imply the inequalities and tail bounds for symmetric polynomials that are proved here.

We present the proofs of our inequalities for symmetric polynomials in Section 2 and tail bounds for symmetric polynomials in Section 3. We use these bounds to prove Theorem 1 on products of low-variance variables in Section 4 and to analyze the [GMR+12] generator in Section 5.

Inequalities for symmetric polynomials

under the constraint that S1(a)S_{1}(a) is fixed. Since ϕk\phi_{k} is projectively defined, its supremum is attained in the (compact) unit sphere, and is therefore a maximum. Choose a≠0a\neq 0 to be a point that achieves the maximum of ϕk\phi_{k}. We assume, without loss of generality, that S1(a)S_{1}(a) is non-negative (if S1(a)<0S_{1}(a)<0, consider −a-a instead of aa). There are two cases to consider:

The first case is that for all i∈[n]i\in[n],

In this case we do not need the induction hypothesis and can in fact replace each aia_{i} by its absolute value. Let P⊆[n]P\subseteq[n] be the set of i∈[n]i\in[n] so that ai≥0a_{i}\geq 0. Then by Equation (9),

The second case is that there exists i0∈[n]i_{0}\in[n] so that

Hence, for all δ\delta close enough to zero so that ∑iδi=0\sum_{i}\delta_{i}=0,

For the above inequality to hold for all such δ\delta, it must be that there is λ\lambda so that for all i∈[n]i\in[n],

To see why this is true, set λi=aiSk(a)k−(S12(a)+E2(a))Sk−1(−i)\lambda_{i}=a_{i}S_{k}(a)k-(S^{2}_{1}(a)+E_{2}(a))S_{k-1}(-i) . We now have λ1,…,λn\lambda_{1},\ldots,\lambda_{n} so that

for every δ1,…,δn\delta_{1},\ldots,\delta_{n} of sufficiently small norm where ∑iδi=0\sum_{i}\delta_{i}=0. We claim that this implies that in fact λi=λ\lambda_{i}=\lambda for every ii. To see this, assume for contradiction that λ1≠λ2\lambda_{1}\neq\lambda_{2} and ∣λ1∣>∣λ2∣|\lambda_{1}|>|\lambda_{2}|. Set

for μ>0\mu>0 sufficiently small. It follows that ∑iδi=0\sum_{i}\delta_{i}=0 and ∑iλiδi=μ(λ1λ2−λ12)<0\sum_{i}\lambda_{i}\delta_{i}=\mu(\lambda_{1}\lambda_{2}-\lambda_{1}^{2})<0 so Equation (12) is violated.

This specifically holds for i0i_{0}, so using (10) we have

To apply induction we need to bound S12(−i0)+E2(−i0)S_{1}^{2}(-i_{0})+E_{2}(-i_{0}) from above. Since

The proof is by reduction to Theorem 2. Assume a1,…,ama_{1},\ldots,a_{m} are nonzero and am+1,…,ana_{m+1},\ldots,a_{n} are zero. Denote a′=(a1,…,am)a^{\prime}=(a_{1},\ldots,a_{m}) and notice that for allFor k>mk>m we have Sk(a)=0S_{k}(a)=0 so there is nothing to prove. k∈[n]k\in[n],

Tail bounds under limited independence

The goal is proving a tail bound on the behaviour of the symmetric functions under limited independence.

We start by obtaining tail estimates, under full independence. Let U\cal{U} denote the distribution over X=(X1,…,Xn)X=(X_{1},\ldots,X_{n}) where X1,…,XnX_{1},\ldots,X_{n} are independent.

Since the expectation of XiX_{i} is zero for all i∈[n]i\in[n],

If 2e1/2tσ≤k1/22e^{1/2}t\sigma\leq k^{1/2} then by the union bound

In the following the underlying probability distribution over XX is D\cal{D}. By Lemma 9, for i∈{k,k+1}i\in\{k,k+1\},

which occurs with probability at least 1−2t−2k1-2t^{-2k}. Fix x=(x1,…,xn)x=(x_{1},\ldots,x_{n}) such that Equation (17) holds.

We claim that there must exist k0∈{0,…,k−1}k_{0}\in\{0,\ldots,k-1\} for which the following bounds hold:

To see this, mark point j∈{0,…,k+1}j\in\{0,\ldots,k+1\} as high if

A point is marked both high and low if equality holds. Observe that is marked high (and low) since S0(x)=1S_{0}(x)=1 and kk and k+1k+1 are marked low by Equation (17). This implies the existence of a triple k0,k0+1,k0+2k_{0},k_{0}+1,k_{0}+2 where the first point is high and the next two are low.

Let γ>0\gamma>0 be the smallest number so that the following inequalities hold:

By definition, one of Equations (21) and (22) holds with equality so

Observe further that γ≤tσ\gamma\leq t\sigma by Equations (18), (19) and (20). Combining this with the bounds in Equations (19) and (20)

Equations (21) and (22) let us apply Theorem 3 with C=γk0+1C=\gamma\sqrt{k_{0}+1} and h≥3h\geq 3 to get

Bounding ∣Sk0∣|S_{k_{0}}| by Equation (23), we get

As in Lemma 11, fix x=(x1,…,xn)x=(x_{1},\ldots,x_{n}) such that Equation (17) holds (the random vector XX has this property with D\cal{D}-probability at least 1−2t−2k1-2t^{-2k}). By the proof of lemma, since by assumption 6etσ<1/26et\sigma<1/2,

For this section, let X1,…,XnX_{1},\ldots,X_{n} be so that each XiX_{i} is uniform over {−1,1}\{-1,1\}. Thus σ2=∑iVar[Xi]=n\sigma^{2}=\sum_{i}\mathsf{Var}[X_{i}]=n. By Lemma 9, we have

implies that for any (2k+2)(2k+2)-wise independent distribution,

Limited independence fools products of bounded variables

In this section we work with the following setup. We have nn random variables X1,…,XnX_{1},\ldots,X_{n} each distributed in the interval $.Let. Let\mu_{i}andand\sigma_{i}^{2}denotethemeanandvarianceofdenote the mean and variance ofX_{i},andlet, and let\sigma^{2}=\sum_{i=1}^{n}\sigma_{i}^{2}.Wewilltypicallyuse. We will typically use\mathcal{U}todenotethedistributionwheretheto denote the distribution where theX_{i}sarefullyindependent,ands are fully independent, and{\mathcal{D}}$ to denote distributions with limited independence.

There exist constants c,c′>0c,c^{\prime}>0 such that under any ckck-wise independent distribution D{\mathcal{D}},

Define H⊂[n]H\subset[n] to be the set of indices such that ∣μi∣≤σ|\mu_{i}|\leq\sqrt{\sigma}. Note that if H≥2kH\geq 2k, then we are done since if c≥2c\geq 2, then

Further, since the variables are bounded in $$, we have

The same bound also holds under U\mathcal{U}, hence

So now assume that ∣H∣≤2k|H|\leq 2k. Let T=H∖[n]T=H\setminus[n]. Even after conditioning on the outcome of variables in HH, the resulting distribution on TT is (c−2)k=c′′k(c-2)k=c^{\prime\prime}k-wise independent. Since the product of variables in HH has absolute value at most 11, it suffices to show that for a c′′kc^{\prime\prime}k-wise independent distribution D{\mathcal{D}},

For ease of notation, we shall assume that T=[m]T=[m] for some m≤nm\leq n. We may assume that m>c′′km>c^{\prime\prime}k else there is nothing to prove.

Let us write Xi=μi(1+Zi)X_{i}=\mu_{i}(1+Z_{i}), so that ZiZ_{i} has mean and variance σi2/μi2\sigma_{i}^{2}/\mu_{i}^{2}. We write

For a c′′kc^{\prime\prime}k-wise independent distribution D{\mathcal{D}},

We first show how to finish the proof of Theorem 13 with this claim. We have

The first two are bounded by (c′σk)/2(c^{\prime}\sigma^{k})/2 by the claim, and the last is since c′′kc^{\prime\prime}k-wise independence fools degree 4k4k polynomials for c′′>4c^{\prime\prime}>4.

Recall that the XiX_{i}s for i∈[m]i\in[m] have expectation μi\mu_{i} where ∣μi∣≥σ|\mu_{i}|\geq\sqrt{\sigma}. We let Xi=μi(1+Zi)X_{i}=\mu_{i}(1+Z_{i}), where ZiZ_{i} has mean and variance σˉi2\bar{\sigma}_{i}^{2} where

Hence the total variance of the ZiZ_{i}s can be bounded by

Let GG denote the event that ∣P(Z)−P′(Z)∣≤2(6eσˉ)4k|P(Z)-P^{\prime}(Z)|\leq 2(6e\sqrt{\bar{\sigma}})^{4k}. Letting t=1/σˉt=1/\sqrt{\bar{\sigma}} and applying Theorem 4, for c′′>8k+2c^{\prime\prime}>8k+2

Analyzing the [GMR+12] generator

Gopalan et al. [GMR+12] proposed and analyzed a PRG{\mathsf{PRG}} for combinatorial rectangles, which we denote by GMR\mathcal{G_{MR}}. In this section, we provide a different analysis of their construction, which is based on our results concerning the symmetric polynomials. Our analysis is simpler and follows the intuition that products of low variance events are easy to fool using limited independence. It also improves one their seedlength in the dependence on n,δn,\delta (see the discussion following Theorem 7).

Let U\mathcal{U} denote the uniform distribution on [m]n[m]^{n}, and let D{\mathcal{D}} be a distribution on [m]n[m]^{n}. For x∈[m]nx\in[m]^{n} and K⊆[n]K\subseteq[n], let xK=(xi:i∈K)x_{K}=(x_{i}:i\in K). We sometimes abuse notation and write xKx_{K} instead of the probability distribution of xKx_{K}. We denote by dTVd_{TV} the total variation distance.

A distribution D{\mathcal{D}} on [m]n[m]^{n} is (k,ε)(k,\varepsilon)-wise independent if for every K⊆[n]K\subseteq[n] of size kk, and x∈D,y∈Ux\in{\mathcal{D}},y\in\mathcal{U}, we have dTV(xK,yK)≤εd_{TV}(x_{K},y_{K})\leq\varepsilon.

Such distributions can be generated using seed length O(log⁡log⁡(n)+klog⁡(m)+log⁡(1/ε))O(\log\log(n)+k\log(m)+\log(1/\varepsilon)) when mm is a power of 22 using standard constructions [NN93]. We can also assume that every co-ordinate is uniformly random in [m][m]. See the appendix for details.

(by adding the string (a,a,…,a)(a,a,\ldots,a) modulo mm, where a∈[m]a\in[m] is uniformly random).

Being (k,ε)(k,\varepsilon)-wise independent is equivalent to saying that for every K⊆[n]K\subseteq[n] of size kk and every f:[m]k→{0,1}f:[m]^{k}\rightarrow\{0,1\},

The following more general property holds. Let PP be a real linear combination of combinatorial rectangles,

We use an alternate view of GMR\mathcal{G_{MR}} as a collection of hash functions g:[n]→[m]g:[n]\rightarrow[m]. The generator GMR\mathcal{G_{MR}} is based on iterative applications of an alphabet increasing step. The first alphabet m0m_{0} is chosen to be large enough, and at each step t>1t>1 the size of the alphabet mtm_{t} is squared mt=mt−12m_{t}=m_{t-1}^{2}. There is a constant C>0C>0 so that the following holds. Denote by δ\delta the error parameter of the generator. Let T≤Clog⁡log⁡(m)T\leq C\log\log(m) be the first integer so that mT≥mm_{T}\geq m. Let δ′=δ/T\delta^{\prime}=\delta/T.

Base Case: Let m0≥Clog⁡(1/δ)m_{0}\geq C\log(1/\delta) be a power of 22. Sample g0:[n]→[m0]g_{0}:[n]\rightarrow[m_{0}] using a (k0,ε0)(k_{0},\varepsilon_{0})-wise independent distribution on [m0]n[m_{0}]^{n} with

This requires seed length O(log⁡log⁡(n)+log⁡(log⁡log⁡(m)/δ)log⁡log⁡(log⁡log⁡(m)/δ))O(\log\log(n)+\log(\log\log(m)/\delta)\log\log(\log\log(m)/\delta)).

Squaring the alphabet: Pick gt′:[mt−1]×[n]→[mt]g^{\prime}_{t}:[m_{t-1}]\times[n]\rightarrow[m_{t}] using a (kt,εt)(k_{t},\varepsilon_{t})-wise independent distribution over [mt]mt−1×n[m_{t}]^{m_{t-1}\times n} with

Define a hash function gt:[n]→[mt]g_{t}:[n]\rightarrow[m_{t}] as

This requires seed length O(log⁡log⁡(n)+log⁡(mt)+log⁡(log⁡log⁡(m)/δ))O(\log\log(n)+\log(m_{t})+\log(\log\log(m)/\delta)).

2 Analyzing the generator

We first analyze the base case using the inclusion-exclusion approach of [EGL+98]. We need to extend their analysis to the setting where the co-ordinates are only approximately kk-wise independent.

Let D{\mathcal{D}} be a (k,ε)(k,\varepsilon)-wise independent distribution on [m]n[m]^{n} with kk odd. Then,

Let pi=∣Si∣/mp_{i}=|S_{i}|/m, and qi=1−piq_{i}=1-p_{i}. Observe that

We consider two cases based on ∑iqi\sum_{i}q_{i}.

Case 1: When ∑iqi≤k/(2e)\sum_{i}q_{i}\leq k/(2e). Since every non-zero qiq_{i} is at least 1/m1/m, there can be at most mk/(2e)mk/(2e) indices ii so that qi>0q_{i}>0. For ii so that qi=0q_{i}=0, we have Si=[m]S_{i}=[m], so we can drop such indices and assume n≤mk/(2e)n\leq mk/(2e). By Bonferroni inequality, since kk is odd,

A similar bound holds for hh. The (k,ε)(k,\varepsilon)-wise independence thus implies

The second term is twice Sk(q1,…,qn)S_{k}(q_{1},\ldots,q_{n}), which we can bound by Maclaurin’s identity as

Case 2: When ∑iqi>k/2e\sum_{i}q_{i}>k/2e. Once again, we drop indices ii so that qi=0q_{i}=0. Consider the largest n′n^{\prime} such that

Repeating the argument from Case 1 for this n′n^{\prime},

To analyze the iterative steps, we use the following lemma:

There is C>0C>0 so that the following holds for δ>0\delta>0 small enough. Assume

If ∣H∣≥k|H|\geq k, we show that the probabilities are small which means that they are close. Indeed, let H′H^{\prime} be the first kk indices in HH. First,

Fooling the Tail:

We may assume that qi≥1/mq_{i}\geq 1/m and pi>0p_{i}>0 for all i∈Ti\in T, since otherwise SiS_{i} is trivial and we can drop such an index. As in the proof of Lemma 16, by restricting to a subset if necessary, we can also assume that

For simplicity of notation, we denote ∣T∣|T| by nn. Therefore, n≤Cmlog⁡(1/δ)n\leq Cm\log(1/\delta).

Since g′(a,i)g^{\prime}(a,i) is uniform over SiS_{i},

We will show that Pk(A)P_{k}(A) is a good approximation to Pn(A)P_{n}(A) under (O(k),εO(1))(O(k),\varepsilon^{O(1)})-wise independence, hence under both D{\mathcal{D}} and U\mathcal{U}.

Plugging in the bounds from Equations (34):

We argue for D{\mathcal{D}}, the same argument holds for U\mathcal{U}. Write

If A1,…,AnA_{1},\ldots,A_{n} are (O(k),0)(O(k),0)-wise independent, then, by Lemma 9,

Hence, under (O(k),ε)(O(k),\varepsilon)-wise independence,

Denote by ¬G\neg G the complement of GG. Write

It remains to bound the second term. Bound

We are ready to prove the main theorem of this section.

The proof uses an hybrid argument. The GMR\mathcal{G_{MR}} generator chooses g0:[n]→[m0]g_{0}:[n]\rightarrow[m_{0}], and then g1′,…,gT′g^{\prime}_{1},\ldots,g^{\prime}_{T} where gt′=[mt−1]×[n]→[mt]g^{\prime}_{t}=[m_{t-1}]\times[n]\rightarrow[m_{t}] has error δ′=δ/T\delta^{\prime}=\delta/T and defines

Let h0,h1′,…,ht′h_{0},h^{\prime}_{1},\ldots,h^{\prime}_{t} be truly random hash functions with similar domains and ranges. For 0≤t,l≤T0\leq t,l\leq T, define the hybrid family Gtl={ftl:[m]→[n]}\mathcal{G}^{l}_{t}=\{f^{l}_{t}:[m]\rightarrow[n]\} as follows: for t=0t=0 and every ll,

For every ll, let Gl=GTl\mathcal{G}^{l}=\mathcal{G}^{l}_{T}. Thus, G0=GMR\mathcal{G}^{0}=\mathcal{G_{MR}} and GT=U\mathcal{G}^{T}=\mathcal{U}. We will show by induction on l≥1l\geq 1 that

The desired bound then follows by the triangle inequality.

In the base case when l=1l=1, couple G0\mathcal{G}^{0} and G1\mathcal{G}^{1} by picking the same g1′,…,gT′g^{\prime}_{1},\ldots,g^{\prime}_{T}, and use them to define the function f′:[m1]×[n]→[m]f^{\prime}:[m_{1}]\times[n]\rightarrow[m] so that

by applying Lemma 16 with k=O(log⁡(1/δ′))k=O(\log(1/\delta^{\prime})) and ε=δ′⋅m0−O(k)\varepsilon=\delta^{\prime}\cdot m_{0}^{-O(k)}.

For the inductive case l>1l>1, couple Gl\mathcal{G}^{l} and Gl−1\mathcal{G}^{l-1} by picking the same gl+1′,…,gT′g^{\prime}_{l+1},\ldots,g^{\prime}_{T}, and pick x∈[ml−1]nx\in[m_{l-1}]^{n} uniformly at random. There is a function f′:[ml]×[n]→[m]f^{\prime}:[m_{l}]\times[n]\rightarrow[m] so that

Acknowledgements

We thank Nati Linial, Raghu Meka, Yuval Peres, Dan Spielman, Avi Wigderson and David Zuckerman for helpful discussions. We thank an anonymous referee for pointing out an error in the statement of Theorem 4 in a previous version of the paper.

References

Appendix A Missing Proofs

Consider p(n−k−1)(ξ)p^{(n-k-1)}(\xi) which is the (n−k−1)th(n-k-1)^{th} derivative of p(ξ)p(\xi). Since Sk(b)=Sk+1(b)=0S_{k}(b)=S_{k+1}(b)=0 for k>0k>0, it follows that ξ2\xi^{2} divides p(n−k−1)(ξ)p^{(n-k-1)}(\xi) and hence mult(p(n−k−1),0)≥2\mathsf{mult}(p^{(n-k-1)},0)\geq 2. Applying the above fact n−k−1n-k-1 times, we get mult(p,0)≥n−k+1\mathsf{mult}(p,0)\geq n-k+1 so Sn(b)=…=Sk(b)=0S_{n}(b)=\ldots=S_{k}(b)=0. ∎

Finally we discuss how to generate the (k,ε)(k,\varepsilon)-wise independent distributions on [m]n[m]^{n} with seed length O(log⁡log⁡(n)+klog⁡(m)+log⁡(1/ε))O(\log\log(n)+k\log(m)+\log(1/\varepsilon)). We claim that it suffices to take a k′=klog⁡(m)k^{\prime}=k\log(m)-wise ε\varepsilon-independent string of length n′=nlog⁡(m)n^{\prime}=n\log(m). Naor and Naor [NN93] showed that such distributions can be generated using seed-length O(log⁡log⁡(n)+klog⁡(m)+log⁡(1/ε))O(\log\log(n)+k\log(m)+\log(1/\varepsilon)). We can also assume that every co-ordinate is uniformly random in [m][m] by adding the string (a,…,a)(a,\ldots,a) where a∈[m]a\in[m] is chosen randomly.