How to refute a random CSP

Sarah R. Allen, Ryan O'Donnell, David Witmer

On refutation of random CSPs

Constraint satisfaction problems (CSPs) play a major role in computer science. There is a vast theory [BJK05] of how algebraic properties of a CSP predicate affect its worst-case satisfiability complexity; there is a similarly vast theory [Rag09] of worst-case approximability of CSPs. Finally, there is a rich range of research — from the fields of computer science, mathematics, and physics — on the average-case complexity of random CSPs; see [Ach09] for a survey just of random kk-SAT. This paper is concerned with random CSPs, and in particular the problem of efficiently refuting satisfiability for random instances. This is a well-studied algorithmic task with connections to, e.g., proof complexity [BB02], inapproximability [Fei02], SAT-solvers [SAT], cryptography [ABW10], learning theory [DLSS14], statistical physics [CLP02], and complexity theory [BKS13].

Historically, random CSPs are probably best studied in the case of kk-SAT, k≥3k\geq 3. The model here involves choosing a CNF formula I\mathcal{I} over nn variables by drawing mm clauses (ORs of kk literals) independently and uniformly at random. (The precise details of the random model are inessential; see Section 3.1 for more information.) This is one of the best known efficient ways of generating hard-seeming instances of NP\mathsf{NP}-complete and coNP\mathsf{coNP}-complete problems. The computational hardness depends crucially on the density, α=m/n\alpha=m/n. For each kk there is (conjecturally) a constant critical density αk\alpha_{k} such that I\mathcal{I} is satisfiable with high probability when α<αk\alpha<\alpha_{k}, and I\mathcal{I} is unsatisfiable with high probability when α>αk\alpha>\alpha_{k}. (Here and throughout, “with high probability (whp)” means with probability 1−o(1)1-o(1) as n→∞n\to\infty.) This phenomenon occurs for all nontrivial random CSPs; in the case of kk-SAT it’s been rigorously proven [DSS15] for sufficiently large kk.

There is a natural algorithmic task associated with each of the two regimes. When α<αk\alpha<\alpha_{k} one wants to find a satisfying assignment for I\mathcal{I}. When α>αk\alpha>\alpha_{k} one wants to refute I\mathcal{I}; i.e., find a certificate of unsatisfiability. Most heuristic SAT-solvers use DPLL-based algorithms; on unsatisfiable instances, they produce certificates that can be viewed as refutations within the Resolution proof system. More generally, a refutation algorithm for density α\alpha is any algorithm that: a) outputs “unsatisfiable” or “fail”; b) never incorrectly outputs “unsatisfiable”; c) outputs “fail” with low probability (i.e., probability o(1)o(1)). We caution the reader that in this paper we do not consider the related, but distinct, scenario of distinguishing planted random instances from truly random ones. Empirical work suggests that as α\alpha increases towards αk\alpha_{k}, finding satisfying assignments becomes more difficult; and conversely, as α\alpha increases beyond αk\alpha_{k}, finding certificates of unsatisfiability gradually becomes easier.

The special case of random 33-XOR has proved particularly important: it is related to 33-SAT refutation through Feige’s “3XOR Principle” (see [Fei02, FO05, FKO06]); it’s the basis for cryptographic schemes due to Alekhnovich [Ale03] (and is related to the “Learning Parities with Noise” problem); it’s used in the best known lower bounds for the SOS SDP hierarchy [Gri01, Sch08], which we discuss further in Section 6; and, Barak and Moitra [BM15] have shown it to be equivalent to a certain “tensor prediction problem” in learning theory.

Our results and techniques

Here we describe our main results and techniques at a high level. Precise theorem statements appear later in the work and the definitions of the terminology we use is given in Section 3. We also mention that in Section B we will generalize all of our results to the case of larger alphabets; but we’ll just discuss Boolean predicates P:{0,1}k→{0,1}P:\{0,1\}^{k}\to\{0,1\} for simplicity.

The proof of Theorem 2.1 follows ideas from [COGL07] and earlier works on “discrepancy” of random kk-SAT instances. The case of even kk is notably easier, and we present two “folklore” arguments for it. The case of odd kk is trickier. Roughly speaking we view the instance as a homogeneous degree-kk multilinear polynomial, which we want to certify takes on only small values on inputs in {−1,1}n\{-1,1\}^{n}. Considering separately the contributions based on the “last” of the kk variables in each constraint, and then using Cauchy–Schwarz, it suffices to bound the norm of a carefully designed quadratic form of dimension nk−1n^{k-1}, indexed by tuples of k−1k-1 variables. This is done using the trace method [Wig55, FK81]. Similar techniques, including the use of the trace method, date back to the 2001 Friedgman–Goerdt work [FG01] refuting random 33-SAT given m=n3/2+ϵm=n^{3/2+\epsilon} constraints.

With Theorem 2.1 in hand, the next step is certifying quasirandomness of random kk-ary CSP instances having m\geq\mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{k/2}) constraints. Roughly speaking we say that a CSP instance is quasirandom if, for every assignment x∈{0,1}nx\in\{0,1\}^{n}, the mm induced kk-tuples of literal values are close to being uniformly distributed over {0,1}k\{0,1\}^{k}. (Note that this is only a property of the instances’ constraint scopes/negations, and has nothing to do with PP.) Since the “Vazirani XOR Lemma” implies that a distribution on {−1,1}k\{-1,1\}^{k} is uniform if and only if all its 2k2^{k} XORs are have bias 12\frac{1}{2}, we are able to leverage Theorem 2.1 to prove:

If an instance is quasirandom, then no solution can be much better than a randomly chosen one. Thus by certifying quasirandomness we are able to strongly refute random instances of any CSP(P)(P):

In particular, this theorem improves upon [COCF10] by a factor of n\sqrt{n} whenever kk is odd; this savings is new even in the well-studied case of kk-SAT.

We remark that property of a predicate PP supporting a pairwise uniform distribution has played an important role in approximability theory for CSPs, ever since Austrin and Mossel [AM09] showed that such predicates are hereditarily approximation-resistant under the UGC. Also, note that the largest tt for which a predicate PP supports a tt-wise uniform distribution determines the minimum number of constraints required by our algorithm. This value is closely related to the notion of distribution complexity studied by Feldman, Perkins, and Vempala [FPV14, FPV15] in the context of planted random CSPs. Informally, the distribution complexity of a planted CSP is the largest tt for which the distribution over constraint inputs {−1,1}k\{-1,1\}^{k} induced by the planted assignment is tt-wise uniform. Despite this similarity, the algorithmic techniques used by Feldman, Perkins, and Vempala in the planted case [FPV14] do not seem to directly apply to refutation.

The idea behind the proof of Theorem 2.4 is that with \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{t/2}) constraints we can use the algorithm of Theorem 2.2 to certify quasirandomness (closeness to uniformity) for all subsets of tt out of kk coordinates. Thus for every assignment x∈{0,1}nx\in\{0,1\}^{n}, the induced distribution on constraint kk-tuples is (o(1)o(1)-close to) tt-wise uniform. Since PP does not support a tt-wise uniform distribution, this essentially shows that no xx can induce a fully-satisfying distribution on constraint inputs. To handle the o(1)o(1)-closeness caveat, we show that if PP does not support a tt-wise uniform distribution, then it is δ\delta-far from supporting such a distribution, for δ=Ωk(1)\delta=\Omega_{k}(1). The algorithm can then in fact (δ−o(1))(\delta-o(1))-refute random CSP(P)(P) instances.

To briefly illustrate the result, consider the Exactly-kk-out-of-2k2k-SAT CSP, studied previously in [BB02, GJ03]. The associated predicate supports a 11-wise uniform distribution, namely the uniform distribution over strings in {0,1}2k\{0,1\}^{2k} of Hamming weight kk. However, it is not hard to show that it does not support any pairwise uniform distribution. As a consequence, random instances of this CSP can be refuted with only \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n) constraints, independent of kk.

Recently, an exciting approach to proving hardness-of-learning results has been developed by Daniely, Linial, and Shalev-Shwartz [DLSS13, DLSS14, DSS14, Dan15]. The most general results appear in [DLSS14]. In this work, Daniely et al. prove computational hardness of several central learning theory problems, based on two assumptions concerning the hardness of random CSP refutation. While the assumptions made in [DSS14, Dan15] appear to be plausible, our work unfortunately shows that the more general assumptions made in [DLSS14] are false.

In [DLSS14] it is shown how to obtain three very notable hardness-of-learning results from these assumptions. However as stated, our work falsifies the SRCSP Assumptions. Indeed, the assumptions are false even in the three specific cases used by [DLSS14] to obtain hardness-of-learning results. We now describe these cases.

The Huang predicates (Hκ)(H_{\kappa}) are arity-Θ(κ3)\Theta(\kappa^{3}) predicates introduced in [Hua13]; they are hereditarily approximation resistant on satisfiable instances and have 00-variability Ω(κ)\Omega(\kappa). In [DLSS14] they are used in SRCSP Assumption 1 to deduce hardness of PAC-learning DNFs with ω(1)\omega(1) terms. However:

For all κ≥9\kappa\geq 9, the predicate HκH_{\kappa} does not support a 44-wise uniform distribution.

Thus by Theorem 2.4 we can efficiently refute random instances of CSP(HκH_{\kappa}) with just \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{2}) constraints, independent of κ\kappa. This contradicts SRCSP Assumption 1.

Finally, we also prove that SRCSP Assumption 1 is false for another family of predicates (Tk)(T_{k}) used by [DLSS14] to show hardness of PAC-learning intersections of 44 Boolean halfspaces.

Our results described in these three cases all use linear programming duality. Specifically, in Lemma 3.16 we show that PP is δ\delta-far from supporting a tt-wise uniform distribution if and only if there exists a kk-variable multilinear polynomial QQ that satisfies certain properties involving PP and δ\delta. We then explicitly construct these dual polynomials for the Huang, Majority, and TkT_{k} predicates.

We conclude this section by emphasizing the importance of the Daniely–Linial–Shalev-Shwartz hardness-of-learning program, despite the above results. Indeed, subsequently to [DLSS14], Daniely and Shalev-Shwartz [DSS14] showed hardness of improperly learning DNF formulas with ω(log⁡n)\omega(\log n) terms under a much weaker assumption than SRCSP Assumption 11. Specifically, their work only assumes that for all dd there is a large enough kk such that refuting random kk-SAT instances is hard when there are m=ndm=n^{d} constraints. This assumption looks quite plausible to us, and may even be true with kk not much larger than 2d2d. Most recently, Daniely showed hardness of approximately agnostically learning halfspaces using the XOR predicate rather than majority [Dan15]. This result shows that there is no efficient algorithm that agnostically learns halfspaces to within a constant approximation ratio under the assumption that refuting random kk-XOR instances is hard when m=ncklog⁡km=n^{c\sqrt{k}\log k} for some c>0c>0. He also shows hardness of learning halfspaces to within an approximation factor of 2log⁡1−νn2^{\log^{1-\nu}n} for any ν>0\nu>0 assuming that there exists some constant c>0c>0 such that for all ss, refuting random kk-XOR instances with k=log⁡snk=\log^{s}n is hard when m=nckm=n^{ck}.

2 Evidence for the optimality of our results

It’s natural to ask whether the nt/2n^{t/2} dependence in our main Theorem 2.4 can be improved. As we don’t expect to prove unconditional hardness results, we instead merely seek good evidence that refuting (t−1)(t-1)-wise supporting predicates PP is hard when m≪nt/2m\ll n^{t/2}. A natural form of evidence would be showing that various strong classes of polynomial-time refutation algorithms fail when m≪nt/2m\ll n^{t/2}. To make sense of this we need to talk about the form of such algorithms; i.e., propositional proof systems.

Recently, there has been significant study of the “SOS” (Sum-Of-Squares) proof system, introduced by Grigoriev and Vorobjov [GV01]; see, e.g., [OZ13, BS14] for discussion. It has the following virtues: a) it is very powerful, being able to efficiently simulate other proof systems (e.g., Resolution, Lovász–Schrijver); b) it is automatizable [Las00, Par00], meaning that nn-variable “degree-dd SOS proofs” can be found in nO(d)n^{O(d)} time whenever they exist; c) we do know some examples of lower bounds for degree-dd SOS proofs. In Section 6 of this paper we show the following:

All of our refutation algorithms for kk-ary predicates can be extended to produce degree-2k2k SOS proofs.

We now return to the question of evidence for the optimality of constraint density used in our results. Dating back to Franco–Paull [FP83] and Chvátal–Szemerédi [CS88], there has been a long line of work in proof complexity showing lower bounds for refuting random 33-SAT instances, especially in the Resolution proof system. This culminated in the work of Ben-Sasson and Wigderson [BSW99], which showed that for random 33-SAT (and 33-XOR) with m=O(n3/2−ϵ)m=O(n^{3/2-\epsilon}), Resolution refutations require size 2nΩ(ϵ)2^{n^{\Omega(\epsilon)}} (whp). More recently, Schoenebeck [Sch08] showed (using the expansion analysis of [BSW99]) that random kk-XOR and kk-SAT instances with m≤nk/2−ϵm\leq n^{k/2-\epsilon} require SOS proofs of degree nΩ(ϵ)n^{\Omega(\epsilon)}, and hence take 2nΩ(ϵ)2^{n^{\Omega(\epsilon)}} time to refute by the “SOS Method”. See [Tul09, Cha13] for related larger-alphabet followups. These results show that the Barak–Moitra \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{k/2}) bound for refuting random kk-XOR (which also works in O(k)O(k)-degree SOS) and our bound for random kk-SAT are tight (up to subpolynomial factors) within the SOS framework. Given the power of the SOS framework, this arguably constitutes some reasonable evidence that no polynomial-time algorithm can refute random kk-SAT instances with m≪nk/2m\ll n^{k/2}.

We now discuss our main theorem’s nt/2n^{t/2} bound for predicates PP not supporting tt-wise uniform distributions. Suppose PP is a predicate that does support a (t−1)(t-1)-wise uniform distribution, where t>2t>2. In the context of inapproximability and SDP-hierarchy integrality gaps, this condition on PP has been significantly studied in the case of t=3t=3. For PP supporting pairwise uniformity, it is known [BGMT12, TW13] that the Sherali–Adams and Lovász–Schrijver+ SDP hierarchies require degree Ω(n)\Omega(n) to refute random CSP(P)(P) instances (whp) when m=O(n)m=O(n). This result was also recently proven for the stronger SOS proof system by Barak, Chan, and Kothari [BCK15], except that the CSP(P)(P) instances are not quite uniformly random; they are “slightly pruned” random instances. For any t>2t>2, the second and third authors recently essentially showed [OW14] that for the Sherali–Adams+ SDP hierarchy, degree nΩ(ϵ)n^{\Omega(\epsilon)} is (whp) necessary to refute random CSP(P)(P) instances when m≤nt/2−ϵm\leq n^{t/2-\epsilon}. As a caveat, again the instances are slightly pruned random instances, rather than being purely uniformly random. (The instances in [OW14] are also in the “Goldreich [Gol00] style”; i.e., there are no literals, but the “right-hand sides” are random. However it is not hard to show the proofs in [OW14] go through in the standard random model of this paper.) Future work [MWW15] is devoted to removing the pruning in these instances. Although the Sherali–Adams+ SDP hierarchy is certainly weaker than the SOS hierarchy, these works still constitute some evidence that our main theorem’s requirement of m\geq\mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{t/2}) for non-tt-wise supporting predicates may be essentially optimal.

3 Certifying independence number and chromatic number of random hypergraphs

Coja-Oghlan, Goerdt, and Lanka [COGL07] also use their CSP refutation techniques to certify that random 33- and 44-uniform hypergraphs have small independence number and large chromatic number. We extend these results to random kk-uniform hypergraphs.

For a random kk-uniform hypergraph HH, there is a polynomial time algorithm certifying that the independence number of HH is at most β\beta with high probability when HH has at least \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}\left(\frac{n^{5k/2}}{\beta^{2k}}\right) hyperedges.

For a random kk-uniform hypergraph HH, there is a polynomial time algorithm certifying that the chromatic number of HH is at least ξ\xi with high probability when HH has at least \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}\left(\xi^{2k}n^{k/2}\right) hyperedges.

The proofs of these theorems follow the outline of [COGL07]. We show Theorem 2.9 using a slightly more general form of Theorem 2.1. Theorem 2.10 follows almost directly from Theorem 2.9 using the fact that every color class of a valid coloring is an independent set. Details are given in Appendix C.

Preliminaries and notation

We next define a standard random model for CSPs. For P:{−1,1}k→{0,1}P:\{-1,1\}^{k}\to\{0,1\}, let FP(n,p)\mathcal{F}_{P}(n,p) be the distribution over CSP instances given by including each of the 2knk2^{k}n^{k} possible constraints independently with probability pp. Note that we may include constraints on different permutations of the same set of variables, constraints on the same tuple of variables with different negations cc, and constraints with the same variable occurring as more than one argument. It is reasonable to include such constraints in the case that the predicate PP is not a symmetric function. We use m‾\overline{m} to denote the expected number of constraints, namely 2knkp2^{k}n^{k}p. As noted in Fact 3.6 below, the number of constraints mm in a draw from FP(n,p)\mathcal{F}_{P}(n,p) is very tightly concentrated around m‾\overline{m}, and we often blur the distinction between these parameters. Appendix D explicitly describes a method for simulating an instance drawn from FP(n,p)\mathcal{F}_{P}(n,p) when the number of constraints is fixed.

We now introduce an important notion for this paper: that of a CSP instance being quasirandom. Versions of this notion originate in the works of Goerdt and Lanka [GL03] (under the name “discrepancy”), Khot [Kho06] (“quasi-randomness”), Austrin and Håstad [AH13] (“adaptive uselessness”), and Chan [Cha13] (“low correlation”), among other places. To define it, we first need to define the induced distribution of an instance and an assignment.

Here we use the notation UkU^{k} for the uniform distribution on {−1,1}k\{-1,1\}^{k} as well as the following:

An immediate consequence of an instance being quasirandom is that its optimum is close to P‾\overline{P}:

We conclude the discussion of CSPs by recording some facts that are proven easily with the Chernoff bound:

Let I∼FP(n,p)\mathcal{I}\sim\mathcal{F}_{P}(n,p). Then the following statements hold with high probability.

m=∣I∣∈m‾⋅(1±O(log⁡nm‾))m=\left\lvert\mathcal{I}\right\rvert\in\overline{m}\cdot\left(1\pm O\left(\sqrt{\frac{\log n}{\overline{m}}}\right)\right).

I\mathcal{I} is O(2k⋅nm‾)O\left(\sqrt{2^{k}\cdot\frac{n}{\overline{m}}}\right)-quasirandom.

2 Algorithms and refutations on random CSPs

where pp is defined by m‾=2knkp\overline{m}=2^{k}n^{k}p. Although A\mathcal{A} is often a deterministic algorithm, we do allow it to be randomized, in which case the above probability is also over the “internal random coins” of A\mathcal{A}.

3 tt-wise uniformity

An important notion for this paper is that of tt-wise uniformity. Recall:

Probability distribution D\mathcal{D} on {−1,1}k\{-1,1\}^{k} is said to be tt-wise uniform, 1≤t≤k1\leq t\leq k, if for all S⊆[k]S\subseteq[k] with ∣S∣=t|S|=t the random variable xSx_{S} is uniform on {−1,1}t\{-1,1\}^{t} when x∼Dx\sim\mathcal{D}. (We remark that this condition is sometimes inaccurately called “tt-wise independence” in the literature.)

We will also consider the more general notion of (ϵ,t)(\epsilon,t)-wise uniformity. This is typically defined using Fourier coefficients:

Probability distribution D\mathcal{D} on {−1,1}k\{-1,1\}^{k} is said to be (ϵ,t)(\epsilon,t)-wise uniform if ∣D^(S)∣≤ϵ|\widehat{D}(S)|\leq\epsilon for all S⊆[k]S\subseteq[k] with 0<∣S∣≤t0<|S|\leq t, where D=2k⋅DD=2^{k}\cdot\mathcal{D} is the probability density associated with distribution D\mathcal{D}.

It is a simple fact (and it follows from Lemma 3.13 below) that (0,t)(0,t)-wise uniformity is equivalent to tt-wise uniformity.

Also important for us is a related but distinct notion, that of being ϵ\epsilon-close to a tt-wise uniform distribution. It’s easy to show that if D\mathcal{D} is ϵ\epsilon-close to a tt-wise uniform distribution then D\mathcal{D} is (2ϵ,t)(2\epsilon,t)-wise uniform. In the other direction, we have the following (see also [AAK+07] for some quantitative improvement):

(Alon–Goldreich-Mansour [AGM03, Theorem 2.1]). If D\mathcal{D} is an (ϵ,t)(\epsilon,t)-wise uniform distribution on {−1,1}k\{-1,1\}^{k}, then there exists a tt-wise uniform distribution D′\mathcal{D}^{\prime} on {−1,1}k\{-1,1\}^{k} with

In particular if t=kt=k we have the bound 2k⋅ϵ2^{k}\cdot\epsilon (and this can also be improved [Gol11] to 2k/2−1⋅ϵ2^{k/2-1}\cdot\epsilon).

A predicate P:{−1,1}k→{0,1}P:\{-1,1\}^{k}\to\{0,1\} is said to be tt-wise supporting if there is a tt-wise uniform distribution D\mathcal{D} whose support is contained in P−1(1)P^{-1}(1). We say PP is δ\delta-far from tt-wise supporting if every tt-wise uniform distribution D\mathcal{D} is δ\delta-far from being supported on PP; i.e., has probability mass at least δ\delta on P−1(0)P^{-1}(0).

4 A dual characterization of limited uniformity

It is known that the condition of PP supporting a tt-wise uniform distribution is equivalent to the feasibility of a certain linear program; hence one can show that PP is not tt-wise supporting by exhibiting a certain dual object, namely a polynomial. This appears, e.g., in work of Austrin and Håstad [AH09, Theorem 3.1]. Herein we will extend this fact by giving a dual characterization of being far from tt-wise supporting.

Q(z)≥δ−1∀z∈{−1,1}kQ(z)\geq\delta-1\quad\forall z\in\{-1,1\}^{k};

Q(z)≥δ∀z∈P−1(1)Q(z)\geq\delta\quad\forall z\in P^{-1}(1);

Q^(∅)=0\widehat{Q}(\emptyset)=0, i.e., QQ has no constant coefficient.

We now provide the quantitative version of the aforementioned [AH09, Theorem 3.1]:

Let P:{−1,1}k→{0,1}P:\{-1,1\}^{k}\rightarrow\{0,1\} and let 0<δ<10<\delta<1. Then PP is δ\delta-far from tt-wise supporting if and only if there is a δ\delta-separating polynomial for PP of degree at most tt.

The proof is an application of linear programming duality. Consider the following LP, which has variables D(z)\mathcal{D}(z) for each z∈{−1,1}kz\in\{-1,1\}^{k}.

minimize ∑z∈{−1,1}k(1−\displaystyle\sum_{\mathclap{z\in\{-1,1\}^{k}}}(1- P(z))D(z)\displaystyle P(z))\mathcal{D}(z) (1) s.t. ∑z∈{−1,1}kD(z)zS=2k⋅D^(S)\displaystyle\sum_{\mathclap{z\in\{-1,1\}^{k}}}\mathcal{D}(z)z^{S}=2^{k}\cdot\widehat{\mathcal{D}}(S) =0\displaystyle=0 ∀S⊆[k]0<∣S∣≤t\displaystyle\forall S\subseteq[k]\quad 0<|S|\leq t (2) ∑z∈{−1,1}kD(z)\displaystyle\sum_{{z\in\{-1,1\}^{k}}}\mathcal{D}(z) =1\displaystyle=1 (3) D(z)\displaystyle\mathcal{D}(z) ≥0\displaystyle\geq 0 ∀z∈{−1,1}k\displaystyle\forall z\in\{-1,1\}^{k}

Constraint (3) and the nonnegativity constraint ensure that D\mathcal{D} is a probability distribution on {−1,1}k\{-1,1\}^{k}. Constraint (2) expresses that D\mathcal{D} is tt-wise uniform (see Remark 3.12). The objective (1) is minimizing the probability mass that D\mathcal{D} puts on assignments in P−1(0)P^{-1}(0). Thus the optimal value of the LP is equal to the smallest γ\gamma such that PP is γ\gamma-close to tt-wise supporting; equivalently, the largest γ\gamma such that PP is γ\gamma-far from tt-wise supporting.

The following is the dual of the above LP. It has a variable c(S)c(S) for each 0<∣S∣≤t0<|S|\leq t as well as a variable ξ\xi corresponding to constraint (3).

maximize ξ\displaystyle\xi (4) s.t. ∑S⊆[k]0<∣S∣≤tc(S)zS\displaystyle\sum_{\begin{subarray}{c}S\subseteq[k]\\ 0<|S|\leq t\end{subarray}}c(S)z^{S} ≤1−P(z)−ξ\displaystyle\leq 1-P(z)-\xi ∀z∈{−1,1}k.\displaystyle\forall z\in\{-1,1\}^{k}. (5)

Observe that a feasible solution ({c(S)}S,ξ)(\{c(S)\}_{S},\xi) is precisely equivalent to a multilinear polynomial QQ of degree at most tt, namely Q(z)=−∑Sc(S)zSQ(z)=-\sum_{S}c(S)z^{S}, that ξ\xi-separates PP.

Thus PP is δ\delta-far from tt-wise supporting if and only if the LP’s objective (1) is at least δ\delta, if and only if the dual’s objective (4) is at least δ\delta, if and only if there is a δ\delta-separating polynomial for PP of degree at most tt. ∎

From this proof we can also derive that if PP fails to be tt-wise supporting then it must in fact be Ωk(1)\Omega_{k}(1)-far from being tt-wise supporting:

Suppose P:{−1,1}k→{0,1}P:\{-1,1\}^{k}\to\{0,1\} is not tt-wise supporting. Then it is in fact δ\delta-far from tt-wise supporting for \delta=2^{-\mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(k^{t})} (or \delta=2^{-\mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(2^{k})} when t=kt=k).

Let K=1+∑i=1t(kt)K=1+\sum_{i=1}^{t}\binom{k}{t} be the number of variables in the dual LP from Lemma 3.16, so K≤kt+1K\leq k^{t}+1 in general, with K≤2kK\leq 2^{k} when t=kt=k. By assumption, the objective (4) of the dual LP’s optimal solution is strictly positive. This optimum occurs at a vertex, which is the solution of a linear system given by a K×KK\times K matrix of ±1\pm 1 entries and a “right-hand side” vector with 0,10,1 entries. By Cramer’s rule, the solution’s entries are ratios of determinants of integer matrices with entries in {−1,0,1}\{-1,0,1\}. Thus any strictly positive entry is at least 1/N1/N, where NN is the maximum possible such determinant. By Hadamard’s inequality, N=KK/2N=K^{K/2} and the claimed result follows. ∎

Quasirandomness and its implications for refutation

For k≥2k\geq 2 and p≥n−k/2p\geq n^{-k/2}, let {w(T)}T∈[n]k\{w(T)\}_{T\in[n]^{k}} be independent random variables such that for each T∈[n]kT\in[n]^{k},

Then there is an efficient algorithm certifying that

In this form, the theorem is not really about CSP refutation at all. It says that the value of a polynomial with random coefficients is close to its expectation when its inputs are bounded.

We give the proof in Appendix A. It follows techniques from [COGL07] fairly closely and is essentially the same as the proof of [BM15]. We will use this theorem to prove our results in subsequent sections.

We obtain strong refutation of kk-XOR as a simple corollary.

so for a kk-XOR instance \mathcal{I}\sim\mathcal{F}_{\textup{k-XOR}}(n,p),

where w(T)=−2−k∑c∈{±1}k1{(T,c)∈I}∏i∈[k]ciw(T)=-2^{-k}\sum_{c\in\{\pm 1\}^{k}}1_{\{(T,c)\in\mathcal{I}\}}\prod_{i\in[k]}c_{i}. The w(T)w(T)’s are random variables depending on the choice of I\mathcal{I}; observe that E\/[w(T)]=0\mathop{\bf E\/}[w(T)]=0, Pr\/[w(T)≠0]≤2kp\mathop{\bf Pr\/}[w(T)\neq 0]\leq 2^{k}p, and ∣w(T)∣≤1\left\lvert w(T)\right\rvert\leq 1 for all T∈[n]kT\in[n]^{k}. By Theorem 4.1, there is an algorithm certifying that

with high probability when p≥n−k/2p\geq n^{-k/2}. Since m=(1+o(1))m‾m=(1+o(1))\overline{m} with high probability, choosing m‾≥2O(k)nk/2log⁡3nγ2\overline{m}\geq\frac{2^{O(k)}n^{k/2}\log^{3}n}{\gamma^{2}} gives the desired result. ∎

2 Quasirandomness and strong refutation of any kk-CSP

Next, we will use the algorithm of Theorem 4.1 to certify that an instance of CSP(P)(P) is quasirandom. This will immediately give us a strong refutation algorithm.

In order to certify quasirandomness, Lemma 3.13 implies that it suffices to certify each Fourier coefficient of DI,xD_{\mathcal{I},x} has small magnitude.

Let ∅≠S⊆[k]\emptyset\neq S\subseteq[k] with ∣S∣=s|S|=s. There is an algorithm that, with high probability, certifies that

for all x∈{−1,1}nx\in\{-1,1\}^{n}, assuming also that m‾≥max⁡{ns/2,n}\overline{m}\geq\max\{n^{s/2},n\}.

To prove this lemma, we need another lemma certifying that polynomials whose coefficients are sums of 00-mean random variables have small value.

The proof is straightforward and we defer it to Section 4.4.

Without loss of generality, assume 1∈S1\in S. Applying definitions, we see that

where we define wS(T,c′)=1{(T,(1,c′))∈I}(c′)S∖{1}−1{(T,(−1,c′))∈I}(c′)S∖{1}w_{S}(T,c^{\prime})=1_{\{(T,(1,c^{\prime}))\in\mathcal{I}\}}(c^{\prime})^{S\setminus\{1\}}-1_{\{(T,(-1,c^{\prime}))\in\mathcal{I}\}}(c^{\prime})^{S\setminus\{1\}} and recall that TST_{S} is the projection of TT onto the the coordinates in SS. It is clear that E\/[wS(T,c′)]=0\mathop{\bf E\/}[w_{S}(T,c^{\prime})]=0, Pr\/[wS(T,c′)≠0]≤p\mathop{\bf Pr\/}[w_{S}(T,c^{\prime})\neq 0]\leq p, and ∣wS(T,c′)∣≤1|w_{S}(T,c^{\prime})|\leq 1. There are τ=2k−1nk−s\tau=2^{k-1}n^{k-s} terms in each sum of wS(T,c′)w_{S}(T,c^{\prime})’s and we can apply Lemma 4.4. When s=2s=2, we plug in these values and see that we can certify that DI,x^(S)≤2O(s)ns/4log⁡5/2nm‾1/2\widehat{D_{\mathcal{I},x}}(S)\leq\frac{2^{O(s)}n^{s/4}\log^{5/2}n}{\overline{m}^{1/2}}. When s=1s=1, m‾≥n\overline{m}\geq n implies that τp≥12\tau p\geq\frac{1}{2} and we can certify that DI,x^(S)≤2O(s)nlog⁡nm‾1/2\widehat{D_{\mathcal{I},x}}(S)\leq\frac{2^{O(s)}\sqrt{n}\log n}{\overline{m}^{1/2}}. The lower bound can be proved in exactly the same way by considering the random variables −wS(T,c′)-w_{S}(T,c^{\prime}). ∎

The existence of an algorithm for certifying quasirandomness follows from Lemmas 3.13 and 4.3.

3 (ϵ,t)(\epsilon,t)-quasirandomness and Ω⁡(1)\Omega(1)-refutation of non-tt-wise-supporting CSPs

In the case that a predicate is not tt-wise supporting, a weaker notion of quasirandomness suffices to obtain Ω(1)\Omega(1)-refutation.

An instance I\mathcal{I} of CSP(P)\textrm{CSP}(P) is (ϵ,t)(\epsilon,t)-quasirandom if DI,x\mathcal{D}_{\mathcal{I},x} is (ϵ,t)(\epsilon,t)-wise uniform for every x∈{−1,1}nx\in\{-1,1\}^{n}.

Fact 3.6 shows that random instances with \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n) constraints are (o(1),t)(o(1),t)-quasirandom for all t≤kt\leq k with high probability. Lemma 4.3 directly implies that we can certify (ϵ,t)(\epsilon,t)-quasirandomness when m\geq\mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{t/2}).

There is an efficient algorithm that certifies that an instance I∼FP(n,p)\mathcal{I}\sim\mathcal{F}_{P}(n,p) of CSP(P)(P) is (γ,t)(\gamma,t)-quasirandom with high probability when m‾≥2O(t)nt/2log⁡5nγ2\overline{m}\geq\frac{2^{O(t)}n^{t/2}\log^{5}n}{\gamma^{2}} and t≥2t\geq 2.

We now reach the main result of this section, which states that if a predicate is δ\delta-far from tt-wise supporting, then we can almost δ\delta-refute instances of CSP(P)(P).

We give two proofs of this theorem. In Proof 1, the theorem follows directly from certification of (γ,t)(\gamma,t)-quasirandomness and Lemma 3.13.

Proof 2 gives a slightly weaker version of Theorem 4.9, requiring the stronger assumption that m‾≥2O(k)nt/2log⁡5nγ2\overline{m}\geq\frac{2^{O(k)}n^{t/2}\log^{5}n}{\gamma^{2}}. It is based on the dual polynomial characterization of being δ\delta-far from tt-wise supporting. While perhaps less intuitive than Proof 1, Proof 2 is more direct. It only uses the XOR refutation algorithm and bypasses [AGM03]’s connection between (ϵ,t)(\epsilon,t)-wise uniformity and ϵ\epsilon-closeness to a tt-wise uniform distribution. We were able to convert Proof 2 into an SOS proof (see Section 6.4), but we did not see how to translate Proof 1 into an SOS version. Proof 2 requires Plancherel’s Theorem, a fundamental result in Fourier analysis.

Since PP is δ\delta-far from tt-wise supporting, there exists a degree-tt polynomial QQ that δ\delta-separates PP. The definition of δ\delta-separating implies that P(z)−(1−δ)≤Q(z)P(z)-(1-\delta)\leq Q(z) for all z∈{−1,1}kz\in\{-1,1\}^{k}. Summing over all constraints, we get that for all x∈{−1,1}nx\in\{-1,1\}^{n},

It then remains to certify that E\/z∈DI,x[Q(z)]≤γ\mathop{\bf E\/}_{\boldsymbol{z}\in\mathcal{D}_{\mathcal{I},x}}[Q(\boldsymbol{z})]\leq\gamma. Observe that

where the second equality follows from Plancherel’s Theorem. Since Q≥−1Q\geq-1 and E\/[Q]=0\mathop{\bf E\/}[Q]=0, Q≤2kQ\leq 2^{k} and hence ∣Q^(S)∣≤2k|\widehat{Q}(S)|\leq 2^{k} for all SS. To finish the proof, we apply Lemma 4.3 to certify that ∣DI,x^(S)∣≤γ22k\left|\widehat{D_{\mathcal{I},x}}(S)\right|\leq\frac{\gamma}{2^{2k}} for all SS. ∎

4 Proof of Lemma 4.4

Let X1,…,XMX_{1},\ldots,X_{M} be independent 00-mean random variables such that ∣Xi∣≤B\left\lvert X_{i}\right\rvert\leq B. Then, for a>0a>0,

Observe that the vUv_{U}’s are independent and that each one is the sum of τ\tau mean-00, i.i.d. random variables with magnitude at most 11. Noting that ∑i=1τE\/[wU(i)2]≤τp\sum_{i=1}^{\tau}\mathop{\bf E\/}[w_{U}(i)^{2}]\leq\tau p, we can use Bernstein’s Inequality to show that the ∣vU∣|v_{U}|’s are not too big with high probability. If s≥2s\geq 2, Theorem 4.1 then implies that the desired algorithm exists. If s=1s=1, we are simply bounding a linear function over ±1\pm 1 variables. We consider two cases: Small pp and large pp.

Choosing a=2slog⁡na=2s\log n in Bernstein’s Inequality, we see that Pr\/[∣vU∣≥2slog⁡n]≤n−2s\mathop{\bf Pr\/}[|v_{U}|\geq 2s\log n]\leq n^{-2s}. A union bound over all UU then implies that \mathop{\bf Pr\/}[\text{any|v_{U}|>2s\log n}]\leq n^{-s}. If s≥2s\geq 2, we observe that Pr\/[vU≠0]≤τp\mathop{\bf Pr\/}[v_{U}\neq 0]\leq\tau p, scale the vUv_{U}’s down by 2slog⁡n2s\log n, and apply Theorem 4.1 to get the stated result. If s=1s=1, we obtain the second bound by observing that

We set a=4sτplog⁡na=4s\sqrt{\tau p}\log n and get that \mathop{\bf Pr\/}[\text{any\left\lvert v_{U}\right\rvert>4s\sqrt{\tau p}\log n}]\leq n^{-s} as above. If s≥2s\geq 2, we can then divide the vUv_{U}’s by 4sτplog⁡n4s\sqrt{\tau p}\log n and apply Theorem 4.1. If s=1s=1, we get a bound of 4τp⋅nlog⁡n4\sqrt{\tau p}\cdot n\log n in the same way as (10). ∎

Hardness of learning implications

Daniely et al. reduce the problem of distinguishing between random instances of CSP(P)\textrm{CSP}(P) and instances with value at least α\alpha as a PAC learning problem by transforming each constraint into a labeled example. To show hardness of improperly learning a certain hypothesis class in the PAC model, they define a predicate PP that is specific to the hypothesis class and assume hardness of distinguishing between random instances of CSP(P)\textrm{CSP}(P) and instances with ndn^{d} constraints and value at least α\alpha for all d>0d>0. They then demonstrate that the sample can be realized (or approximately realized) by some function in the hypothesis class if the CSP instance is satisfiable (or has value at least α\alpha). They also show that if the given CSP instance is random, the set of examples will have error at least 14\tfrac{1}{4} (in the agnostic case 15)\tfrac{1}{5}) for all hh in the hypothesis class with high probability. Using this approach, they obtain hardness results for the following problems: improperly learning DNF formulas, improperly learning intersections of 4 halfspaces, and improperly approximately agnostically learning halfspaces for any approximation factor.

The hardness assumptions made in [DLSS14] are the same as those presented in Section 2.1, except for a few minor differences. First, their model fixes the number of constraints rather than the probability with which each constraint is included in the instance. It is well-known that results in one model easily translate to the other. We include a proof in Appendix D for completeness. Additionally, SRCSP Assumptions 1 and 2 purport hardness of distinguishing random instances of CSP(P)\textrm{CSP}(P) from satisfiable instances, even when the algorithm is allowed to err with probability 14\tfrac{1}{4} over its internal coins. The algorithms presented in the preceding sections never err on satisfiable instances; further, they only fail to certify random instances with probability o(1)o(1). As a result, our refutation algorithms also falsify weaker versions of both SRCSP Assumptions, wherein the allowed probability of error is both lower and one-sided. For each predicate presented in [DLSS14], we falsify the appropriate SRCSP assumption using the following approach. For each predicate PP and corresponding δ>0\delta>0 , we define a degree-tt polynomial that δ\delta-separates PP. Using the refutation techniques presented in the preceding sections, we deduce that \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{t/2}) constraints are sufficient to distinguish random instances of CSP(P)\textrm{CSP}(P) from those that are satisfiable (or have value at least α\alpha). In order to simplify the presentation, we begin with simpler versions of the polynomials and then scale them to attain the appropriate values of δ\delta. The following lemma will be of use for this scaling.

Define Q(z)=Q(z)θ1−θ0\mathcal{Q}(z)=\frac{Q(z)}{\theta_{1}-\theta_{0}}. Clearly Q\mathcal{Q} is also unbiased and has degree tt. Then for all z∈P1z\in P_{1}, Q(z)θ1−θ0≥θ1θ1−θ0\frac{Q(z)}{\theta_{1}-\theta_{0}}\geq\frac{\theta_{1}}{\theta_{1}-\theta_{0}}. Similarly, for all zz, Q(z)θ1−θ0≥θ0θ1−θ0=−θ1−θ0θ1−θ0+θ1θ1−θ0=−1+θ1θ1−θ0\frac{Q(z)}{\theta_{1}-\theta_{0}}\geq\frac{\theta_{0}}{\theta_{1}-\theta_{0}}=-\frac{\theta_{1}-\theta_{0}}{\theta_{1}-\theta_{0}}+\frac{\theta_{1}}{\theta_{1}-\theta_{0}}=-1+\frac{\theta_{1}}{\theta_{1}-\theta_{0}}. ∎

We now demonstrate that the above can be applied to the predicates suggested in [DLSS14] by defining separating polynomials and applying Theorem 4.9

2 Huang’s predicate and hardness of learning DNF formulas

In order to obtain hardness of improperly learning DNF formulas with ω(1)\omega(1) terms, Daniely et al. use the following predicate, introduced by Huang [Hua13]. Huang showed that it is hereditarily approximation resistant; Daniely et al. also observed that its 00-variability is Ω(k1/3)\Omega(k^{1/3}) [DLSS14].

Daniely et al. reduce the problem of distinguishing between random instances of CSP(Hκ)\textrm{CSP}(H_{\kappa}) with 2nd2n^{d} constraints and satisfiable instances to the problem of improperly PAC learning the class of DNF formulas with ω(1)\omega(1) terms on a sample of O(nd)O(n^{d}) training examples with error ϵ=1/5\epsilon=1/5 with probability at least 34\tfrac{3}{4}. Here we show that there exists a polynomial time algorithm that refutes random instances of CSP(Hκ)\textrm{CSP}(H_{\kappa}) by demonstrating that HkH_{k} does not support a 4-wise uniform distribution and applying Theorem 4.9.

As a notational shorthand, write zabcz_{abc} for z{ia,ib,ic}z_{\{i_{a},i_{b},i_{c}\}}. Define ζ:[κ]6×{−1,1}k→\zeta:[\kappa]^{6}\times\{-1,1\}^{k}\rightarrow as follows:

Observe that for each monomial zT1zT2zT3zT4z_{T_{1}}z_{T_{2}}z_{T_{3}}z_{T_{4}} of ζ\zeta, for every j∈j\in, ∑i=141{Ti∋j}=2\sum_{i=1}^{4}1_{\{T_{i}\ni j\}}=2. Further, for each T⊆T\subseteq with ∣T∣=3|T|=3, zTz_{T} appears exactly once in ζ\zeta. Let Z6\mathcal{Z}_{6} be the set of all ordered 6-tuples of distinct elements of [κ][\kappa]. For an ordered tuple II, we use ∈()\in_{()} to denote membership in II.

Observe that QQ does not depend on any of z{1},…z{κ}z_{\{1\}},\ldots z_{\{\kappa\}}. By construction, QQ contains no constant term, so Q^(∅)=0\widehat{Q}(\emptyset)=0. Clearly Q(z)≥−5Q(z)\geq-5 for all zz because (11) is always at least −5-5.

Now we lower bound the value of QQ on all zz that satisfy HκH_{\kappa}. We first show that for any z′z^{\prime} that strongly satisfies the Huang predicate, Q(z′)=5Q(z^{\prime})=5, then bound Q(z′)−Q(z)Q(z^{\prime})-Q(z) for any zz with Hamming distance at most κ\kappa from z′z^{\prime}. By definition, for each zTi′z^{\prime}_{T_{i}}, we have that zTi′∏j∈Tizj′=1z^{\prime}_{T_{i}}\prod_{j\in T_{i}}z^{\prime}_{j}=1. So for each monomial of QQ,

where the last line follows from the fact that ∑1=141{Ti∋j}=2.\sum_{1=1}^{4}1_{\{T_{i}\ni j\}}=2. Because there are 5⋅∣Z6∣5\cdot|\mathcal{Z}_{6}| monomials in QQ, their sum is 55.

Now we consider the case where zz does not strongly satisfy the Huang Predicate, but Hκ(z)=1H_{\kappa}(z)=1. Any singleton index on which zz and z′z^{\prime} differ will not change the value of QQ. Let N={T:zT≠zT′}N=\{T:z_{T}\neq z^{\prime}_{T}\}. We lower bound QQ by counting the number of monomials in which each zTz_{T} appears and

For fixed TT, the number of monomials containing the variables of zTz_{T} is

because there are exactly 120120 ways to permute the three indices of TT in I and the remaining κ−3\kappa-3 indices are permuted in the remaining 3 positions of II. So

If we instead choose to scale QQ by a factor of 15⋅κ2−3κ+22κ2−6κ−44\tfrac{1}{5}\cdot\frac{\kappa^{2}-3\kappa+2}{2\kappa^{2}-6\kappa-44} rather than substituting κ=9\kappa=9 into (12), we can achieve a better separation of δ=κ2−3κ−462κ2−6κ−44\delta=\tfrac{\kappa^{2}-3\kappa-46}{2\kappa^{2}-6\kappa-44}. For κ≥9\kappa\geq 9, this expression is strictly increasing and it approaches 12\tfrac{1}{2} as κ\kappa grows.

3 Hamming weight predicates

The remaining predicates we would like to examine are symmetric, meaning they are functions only of their Hamming weights. Again for each predicate PP we present a multivariate polynomial that δ\delta-separates PP for some 0≤δ≤10\leq\delta\leq 1. Each of these polynomials can also be written as a univariate polynomial on the Hamming weight of its input, which we will use to show that each of the following polynomials δ\delta-separates its predicate for the appropriate value of δ\delta. We give the construction below.

For z∈{−1,1}kz\in\{-1,1\}^{k} where z=z1,…,zkz=z_{1},\ldots,z_{k}, define Sz=∑i=1kziS_{z}=\sum_{i=1}^{k}z_{i} and call SzS_{z} the Hamming weight of zz.

Note that this is analogous to the notion of a Hamming weight of a vector in {0,1}k\{0,1\}^{k}, but differs in that it is not simply the count of the number of 11’s. We define a general predicate that is satisfied when SzS_{z} is at least some fixed threshold value θ\theta.

Because the multilinear separating polynomials we will use are symmetric, we present a transformation to an equivalent univariate polynomial on the Hamming weight of the original input.

Then Q(z)=Qu(Sz)Q(z)=Q_{u}(S_{z}) for all z∈{−1,1}kz\in\{-1,1\}^{k}.

where Ki(ν;k)=∑j=0i(−1)j(νi)(k−νi−j)\mathscr{K}_{i}(\nu;k)=\sum_{j=0}^{i}(-1)^{j}{\nu\choose i}{k-\nu\choose i-j} denotes the Krawtchouk polynomial of degree ii [Kra29, KL96]. Substituting ν=k−Sz2\nu=\frac{k-S_{z}}{2}, yields the following expressions. In [KL96] the first three expressions are given explicitly and the fourth can be easily obtained by applying their recursive formula.

Finally, substituting these expressions into (14) and by some algebra,

As a consequence, by choosing values of a,b,c,a,b,c, and dd, we can work with a univariate polynomial while ensuring that its multivariate analogue is unbiased and has degree at most 4 (degree 3 when d=0d=0).

Daniely et al. define the following predicate in order to show hardness of improperly learning intersections of four halfspaces.

Then by Lemma 5.8, for all z∈{−1,1}kz\in\{-1,1\}^{k}, Q(z)=Qu(Sz)Q(z)=Q_{u}(S_{z}). It therefore suffices to lower bound Qu(s)Q_{u}(s) both when s≥−1s\geq-1 and for all s∈[−k,k]s\in[-k,k].

First we show that QuQ_{u} is monotonically increasing in ss.

which is evidently positive for k≥5k\geq 5.

Because QQ is monotonically increasing in ss, Qu(s)≥Qu(−k)Q_{u}(s)\geq Q_{u}(-k) for all s∈[−k,k]s\in[-k,k].

which is clearly negative for k≥5k\geq 5. Now it just remains to lower-bound Qu(s)Q_{u}(s) for s≥−1s\geq-1. Again, since QuQ_{u} is monotonically increasing in ss, we use the value Qu(−1)Q_{u}(-1):

Because VAR0(I8k)\textrm{VAR}_{0}(I_{8k}) is evidently Ω(k)\Omega(k) and I8k‾<17\overline{I_{8k}}<\tfrac{1}{7} for all k≥5k\geq 5, we have the following Corollary.

For odd k≥5k\geq 5 and sufficiently large nn, there exists an efficient algorithm that distinguishes between random instances of CSP(I8k)\textrm{CSP}(I_{8k}) with \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{3/2}) constraints and satisfiable instances with high probability.

4 Majority and hardness of approximately agnostically learning halfspaces

Then by Lemma 5.8, for all z∈{−1,1}kz\in\{-1,1\}^{k}, Q(z)=Qu(Sz)Q(z)=Q_{u}(S_{z}). To simplify QQ, let σ=sk−1/2\sigma=sk^{-1/2}. Then we can rewrite (18) as follows:

where the last inequality follows from the fact that k≥24k\geq 24 and the first two terms are always nonnegative.

Next we lower-bound Qu(s)Q_{u}(s) for s>0s>0.

5 Predicates satisfied by strings with Hamming weight at least −Θ⁡(k)-\Theta(\sqrt{k}).

In light of the fact that the threshold based predicates above are not 44-wise supporting, one may attempt to find another threshold-based predicate. Here we show that a symmetric threshold predicate that is 44-wise supporting must be satisfied by all strings with Hamming weight at least −k2-\tfrac{\sqrt{k}}{2}. Furthermore, there exists a symmetric threshold predicate that is 44-wise supporting with a threshold of −Θ(k)-\Theta(\sqrt{k}) and we sketch its construction.

Again, for simplicity we set σ=sk−1/2\sigma=sk^{-1/2} and obtain the following expression:

Observe that for k≥99k\geq 99, 1k(83σ2+23σ−2)=23k((2σ−14)2−4916)>−148\tfrac{1}{k}\left(\tfrac{8}{3}\sigma^{2}+\tfrac{2}{3}\sigma-2\right)=\tfrac{2}{3k}\left(\left(2\sigma-\tfrac{1}{4}\right)^{2}-\tfrac{49}{16}\right)>-\tfrac{1}{48}. We now lower-bound the value of QuQ_{u} for s≥−12k1/2s\geq-\tfrac{1}{2}k^{1/2}, or equivalently, σ≥−12\sigma\geq-\tfrac{1}{2}:

The first three terms are always nonnegative, so Qu(s)≥−143Q_{u}(s)\geq-\tfrac{14}{3}.

Applying Lemma 5.1, Tk−12kT_{k}^{-\tfrac{1}{2}\sqrt{k}} is 1255\tfrac{1}{255}-far from supporting a 44-wise uniform distribution. ∎

Assume k=2m−1k=2^{m}-1 for some integer m≥3m\geq 3. Then there exists a 44-wise uniform distribution supported only on z∈{−1,1}kz\in\{-1,1\}^{k} such that Sz≥1−2k+1S_{z}\geq 1-2\sqrt{k+1}.

Let C\mathcal{C} be a binary BCH code of length kk with designed distance 2ι+12\iota+1 and let C⊥\mathcal{C}^{\bot} be its dual. Then the uniform distribution on the codewords of C\mathcal{C} is 2ι2\iota-wise uniform [ABI86, MS77]; see also [AS04, Ch 16.2].

Let c=c1…ckc=c_{1}\ldots c_{k} be a codeword of C⊥,\mathcal{C}^{\bot}, where each ci∈{−1,1}c_{i}\in\{-1,1\}. The Carlitz-Uchiyama bound [MS77, page 280] states that for all c∈C⊥c\in\mathcal{C}^{\bot},

Observe that the quantity 12(1−ci)\tfrac{1}{2}(1-c_{i}) simply maps cic_{i} from {−1,1}\{-1,1\} to {0,1}\{0,1\} so that we can write the bound to match the presentation in [MS77]. Therefore,

Setting ι=2\iota=2, we can obtain 44-wise uniformity on this distribution and each string in the support of the distribution has Hamming weight at least −1−2k+1.-1-2\sqrt{k+1}. ∎

In order to construct a 44-wise uniform distribution for any value of kk, one could simply express kk as a sum of powers of 2, construct separate distributions on disjoint variables as described above for each power of 2 (down to the minimum length for which we can achieve distance at least 55, after which point we use the uniform distribution, and obtain a 44-wise uniform distribution. The total Hamming weight of a vector supported by this distribution would then be at least −O(k)-O(\sqrt{k}).

SOS refutation proofs

with deg⁡(u0)≤d\deg(u_{0})\leq d, deg⁡(uiqi)≤d\deg(u_{i}q_{i})\leq d for all i∈[m]i\in[m], and deg⁡(viri)≤d\deg(v_{i}r_{i})\leq d for all i∈[m′]i\in[m^{\prime}]. If it also holds that u0,u1,…,um=0u_{0},u_{1},\ldots,u_{m}=0, we will write A⊢ds=0A\vdash_{d}s=0.

It is well-known that a degree-dd SOS proof can be found using an SDP of size nO(d)n^{O(d)} if it exists [Sho87, Par00, Las00, Las01].

In this section, we will take the set AA to be {xi2=1}i∈[n]\{x_{i}^{2}=1\}_{i\in[n]}, enforcing that variables are ±1\pm 1-valued. We show that with high probability there exists a low-degree SOS proof that a polynomial representing the value of a CSP instance is close to its expectation.

For more information on the SOS proof system and its applications to approximation algorithms, see, e.g., [OZ13, Lau09].

2 SOS certification of quasirandomness

All of our SOS results rely on the following theorem, which is the SOS version of Theorem 4.1.

For k≥2k\geq 2 and p≥n−k/2p\geq n^{-k/2}, let {w(T)}T∈[n]k\{w(T)\}_{T\in[n]^{k}} be independent random variables such that for each T∈[n]kT\in[n]^{k},

This theorem was essentially proven by Barak and Moitra [BM15]. We give a proof in Appendix A.3. We first use this theorem to show that an SOS version of Lemma 4.4 holds.

We sketch the differences from the proof of Lemma 4.4 given in Section 4.4. For s≥2s\geq 2, the lemma follows by using Theorem 6.2 instead of Theorem 4.1. If s=1s=1, it suffices to show that

for any vv since summing over all ii as in (10) finishes the proof. If vi≥0v_{i}\geq 0, observe that

If v(i)<0v(i)<0, we use (xi+1)2(x_{i}+1)^{2} instead of (xi−1)2(x_{i}-1)^{2}. ∎

The lemma implies an SOS version of Lemma 4.3. To make this precise, we define a specific polynomial representation of DI,x^(S)\widehat{D_{\mathcal{I},x}}(S):

where xTS=∏i∈SxTix_{T}^{S}=\prod_{i\in S}x_{T_{i}}. Note that this is a polynomial in the xix_{i}’s.

We can show these polynomials are not too large.

Let ∅≠S⊆[k]\emptyset\neq S\subseteq[k] with ∣S∣=s|S|=s. Then

with high probability, assuming also that m‾≥max⁡{ns/2,n}\overline{m}\geq\max\{n^{s/2},n\}.

Based on Lemma 3.13, we will think of Lemma 6.4 as giving an SOS proof of quasirandomness. Below, we use it to prove SOS versions of Theorems 4.6 and 4.9.

3 Strong refutation of any kk-CSP

We can then give an SOS proof strongly refuting CSP(PP).

with high probability when m‾≥2O(k)nk/2log⁡5nγ2\overline{m}\geq\frac{2^{O(k)}n^{k/2}\log^{5}n}{\gamma^{2}}.

Note that this is just Plancherel’s Theorem in SOS. The theorem then follows from Lemma 6.4 and the observation that ∑S⊆[k]∣P^(S)∣≤2O(k)\sum_{S\subseteq[k]}|\widehat{P}(S)|\leq 2^{O(k)}. ∎

4 Ω⁡(1)\Omega(1)-refutation of non-tt-wise supporting CSPs

with high probability when m‾≥2O(k)nt/2log⁡5nγ2\overline{m}\geq\frac{2^{O(k)}n^{t/2}\log^{5}n}{\gamma^{2}} and t≥2t\geq 2.

The proof is an SOS version of Proof 2 of Theorem 4.9 above. Claim 6.7 implies that for QQ of degree at most tt that δ\delta-separates PP,

Rearranging terms as in the proof of Theorem 6.5, we see that the right hand side is equal to

Since QQ has mean 00, ∣Q∣≤2k|Q|\leq 2^{k} and ∑S⊆[k]∣P^(S)∣≤2O(k)\sum_{S\subseteq[k]}|\widehat{P}(S)|\leq 2^{O(k)}. The theorem then follows from Lemma 6.4. ∎

with high probability when \overline{m}\geq 2^{\mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(k^{t})}n^{t/2}\log^{5}n and t≥2t\geq 2.

Directions for future work

Additionally, it would be good to investigate whether our refutation algorithms can be extended from the purely random CSP(P)(P) setting to the “smoothed”/“semi-random” setting of Feige [Fei07], in which the mm constraints scopes are worst-case and only the negation pattern for literals is random. Feige showed how to efficiently refute random 33-SAT instances with m\geq\mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{3/2}) constraints even in this model.

Followup work on the very interesting paper [FKO06] of Feige, Kim, and Ofek also seems warranted. Recall that it gives a nondeterministic refutation algorithm for random 33-SAT when m≥O(n1.4)m\geq O(n^{1.4}) (as well as a subexponential-time deterministic algorithm). This raises the question of whether there exist polynomial-size refutations for random CSP(P)(P) instances that are nevertheless hard to find efficiently.

The authors would like to thank Amin Coja–Oghlan for help with the literature, and Boaz Barak and Ankur Moitra for permission to reprint the proof of the strong kk-XOR refutation result. The last author would like to thank Anupam Gupta for several helpful discussions.

References

Appendix A Proof of Theorem 4.1

For k≥2k\geq 2 and p≥n−k/2p\geq n^{-k/2}, let {w(T)}T∈[n]k\{w(T)\}_{T\in[n]^{k}} be independent random variables such that for each T∈[n]kT\in[n]^{k},

Then there is an efficient algorithm certifying that

When kk is even, we can think of ∑T∈[n]kw(T)xT\sum_{T\in[n]^{k}}w(T)x^{T} as a quadratic form:

where yU=xUy_{U}=x^{U}. We give two methods to certify that the value of this quadratic form is at most Ok(pn3k/4log⁡n)O_{k}(\sqrt{p}n^{3k/4}\log{n}). The first method uses an SDP-based approximation algorithm and works only for x∈{−1,1}nx\in\{-1,1\}^{n}. The second method uses ideas from random matrix theory and works for any xx with ∥x∥∞≤1\left\lVert x\right\rVert_{\infty}\leq 1.

If x∈{−1,1}nx\in\{-1,1\}^{n}, we can apply an approximation algorithm of Charikar and Wirth [CW04] for quadratic programming. They prove the following theorem:

[CW04, Theorem 1] Let MM be any n×nn\times n matrix with all diagonal elements 00. There exists an efficient randomized algorithm that finds y∈{−1,1}ny\in\{-1,1\}^{n} such that

By Markov’s Inequality, this statement holds with probability at least 1/21/2. We can run the algorithm O(log⁡n)O(\log n) times to get a high probability result. To apply Theorem A.1, we separate out the diagonal terms of (28), rewriting it as

We can certify that each of the two terms in this expression is at most O(pn3k/4log⁡n)O(\sqrt{p}n^{3k/4}\log n). For the first term, we will need the following claim.

This follows from applying Bernstein’s Inequality (Theorem 4.12) for fixed yy and then taking a union bound over all y∈{−1,1}nk/2y\in\{-1,1\}^{n^{k/2}}. Using the claim, we see that Theorem A.1 allows us to certify that the value of the first term in (29) is at most O(pn3k/4log⁡n)O(\sqrt{p}n^{3k/4}\log n).

We will use the next claim to bound the second term of (29).

Since the ∣wT∣≤1|w_{T}|\leq 1 and Pr\/[wT≠0]≤p\mathop{\bf Pr\/}[w_{T}\neq 0]\leq p, the claim follows from the Chernoff Bound. The second term of (29) is upper bounded by ∑U∈[n]k/2∣wU,U∣\sum_{U\in[n]^{k/2}}|w_{U,U}| and we can compute this quantity in polynomial time to certify that its value is at most O(pn3k/4)O(\sqrt{p}n^{3k/4}).

Observe that (28) is y⊤Byy^{\top}By for a matrix BB indexed by U∈[n]k/2U\in[n]^{k/2} so that BU1,U2=w(U1,U2)B_{U_{1},U_{2}}=w(U_{1},U_{2}). Then y⊤By≤∥B∥∥y∥2y^{\top}By\leq\left\lVert B\right\rVert\left\lVert y\right\rVert^{2}. To certify that y⊤Byy^{\top}By is small, we compute ∥B∥\left\lVert B\right\rVert. We need to show that ∥B∥\left\lVert B\right\rVert is small with high probability. First, note that ∥B∥\left\lVert B\right\rVert is equal to the norm of the 2nk/2×2nk/22n^{k/2}\times 2n^{k/2} symmetric matrix

[Tao12, Proposition 2.3.13] Let MM be a random symmetric matrix n×nn\times n whose upper triangular entries MijM_{ij} with i≥ji\geq j are independent random variables with mean 00, variance at most 11, and magnitude at most KK. Then, with high probability,

A.2 The odd arity case

Fix an assignment x∈nx\in^{n}. For i∈[n]i\in[n], the monomials containing xix_{i} can contribute at most Wi:=∣∑T∈[n]k−1w(T,i)xT∣W_{i}:=\left\lvert\sum_{T\in[n]^{k-1}}w(T,i)x^{T}\right\rvert to the objective if xix_{i} is set optimally. By Cauchy-Schwarz,

so it suffices to bound ∑i∈[n]Wi2\sum_{i\in[n]}W_{i}^{2}. We will write this as a quadratic polynomial and then bound it using spectral methods:

Define the nk−1×nk−1n^{k-1}\times n^{k-1} matrix AA indexed by [n]k−1[n]^{k-1}:

The first term is at most ∥A∥nk−1\left\lVert A\right\rVert n^{k-1} since the variables are bounded. We can compute ∥A∥\left\lVert A\right\rVert to certify this. With high probability, ∥A∥\left\lVert A\right\rVert is not too big.

Let k≥3k\geq 3 and p≥n−k/2p\geq n^{-k/2}. Let {w(T)}T∈[n]k\{w(T)\}_{T\in[n]^{k}} be indepedent random variables satisfying conditions (6), (7), and (8) above. Let AA be defined as in (32). With high probability,

We can therefore certify that the first term is 2O(k)pn3k/2−1log⁡3n2^{O(k)}pn^{3k/2-1}\log^{3}n. We will prove the lemma in Appendix A.4.

The second term of (33) is at most ∑T∈[n]kw(T)2\sum_{T\in[n]^{k}}w(T)^{2}. We can easily compute this and the Chernoff Bound implies that its value is at most pn3k/2−1pn^{3k/2-1} with high probability.

So far, with high probability we can certify that ∑i∈[n]Wi2=2O(k)pn3k/2−1log⁡3n\sum_{i\in[n]}W_{i}^{2}=2^{O(k)}pn^{3k/2-1}\log^{3}n. Plugging this bound into (30) concludes the proof.

It would have been more natural to have written ∑i∈[n]Wi2=(x⊗k−1)⊤A′x⊗k−1\sum_{i\in[n]}W_{i}^{2}=(x^{\otimes k-1})^{\top}A^{\prime}x^{\otimes k-1} for A′A^{\prime} such that AT,U′=∑i∈[n]w(T,i)w(U,i)A^{\prime}_{T,U}=\sum_{i\in[n]}w(T,i)w(U,i). However, ∥A′∥\left\lVert A^{\prime}\right\rVert could be too large because of the contribution of the second term in (33). We use the additional assumption that ∥x∥∞≤1\left\lVert x\right\rVert_{\infty}\leq 1 to get around this issue.

A.3 An SOS version

In this section, we will prove the SOS version of Theorem 4.1.

For k≥2k\geq 2 and p≥n−k/2p\geq n^{-k/2}, let {w(T)}T∈[n]k\{w(T)\}_{T\in[n]^{k}} be independent random variables such that for each T∈[n]kT\in[n]^{k},

Rather than writing out the full proof, we will indicate the small changes required to convert the above proof of Theorem 4.1 into SOS form.

The random matrix proof for the even case can easily be converted into an SOS proof with degree kk. When O(kpnk/4log⁡n)I−B⪰0O(k\sqrt{p}n^{k/4}\log{n})I-B\succeq\mathbf{0}, there exists a matrix MM such that M⊤M=O(kpnk/4log⁡n)I−BM^{\top}M=O(k\sqrt{p}n^{k/4}\log{n})I-B. Then

A couple of additional issues arise in the odd case. First of all, the square root in (30) is not easily expressed in SOS, so we instead prove the squared version

By a simple extension of [OZ13, Fact 3.3], (34) implies (27) in SOS :

Secondly, we do not know how to prove the Cauchy-Schwarz inequality (30) in SOS. However, O’Donnell and Zhou show that a very similar inequality can be proved in SOS [OZ13, Fact 3.8]:

Using this fact instead of Cauchy-Schwarz to prove the squared version of (30), we can follow the argument above to show that

The norm bound can be proven in SOS exactly as in the even case.

A.4 Proof of Lemma A.5

We restate the definition of the matrix AA and the statement of the lemma.

Let k≥3k\geq 3 and p≥n−k/2p\geq n^{-k/2}. Let {w(T)}T∈[n]k\{w(T)\}_{T\in[n]^{k}} be indepedent random variables satisfying conditions (6), (7), and (8) above. Let AA be the [n]k−1×[n]k−1[n]^{k-1}\times[n]^{k-1} indexed by [n]k−1[n]^{k-1} that is defined as follows:

Recall that we index AA by elements of nk−1n^{k-1} divided into two blocks of k−12\frac{k-1}{2} coordinates each. First, note that

Expanding this out using the definition of AA and setting wT=w(T)w_{T}=w(T), we get that

To bound this sum, we will start by bounding E\/[PJ,L]\mathop{\bf E\/}[P_{J,L}]. We will need two claims.

The two claims also imply two other facts we will need below.

If ∣L∣>r|L|>r, then E\/[PJ,L]=0\mathop{\bf E\/}[P_{J,L}]=0.

If ∣J∣>2r+2|J|>2r+2, then E\/[PJ,L]=0\mathop{\bf E\/}[P_{J,L}]=0.

This can be proved in exactly the same manner.

Next, observe that the number of choices of JJ with ∣J∣=a|J|=a is at most na(k−1)2a4r≤na(k−1)2(4r)4rn^{\frac{a(k-1)}{2}}a^{4r}\leq n^{\frac{a(k-1)}{2}}(4r)^{4r}. The number of choices of LL with ∣L∣=b|L|=b is at most nbb2r≤nb(2r)2rn^{b}b^{2r}\leq n^{b}(2r)^{2r}. All together, we can write

Recall that we assumed nkp2≥1n^{k}p^{2}\geq 1. Since a≤2r+2a\leq 2r+2 and b≤rb\leq r, the claim follows. ∎

If we did not have conditions (35), (36), and (37), we would only have been able to show that ∣L∣≤2r|L|\leq 2r. This would have led to a weaker bound of O(n)O(\sqrt{n}).

Appendix B Extension to larger alphabets

Let I∼Fq,P(n,p)\mathcal{I}\sim\mathcal{F}_{q,P}(n,p). Then the following statements hold with high probability.

m=∣I∣∈m‾⋅(1±O(log⁡nm‾))m=\left\lvert\mathcal{I}\right\rvert\in\overline{m}\cdot\left(1\pm O\left(\sqrt{\frac{\log n}{\overline{m}}}\right)\right).

I\mathcal{I} is O(qklog⁡q⋅nm‾)O\left(\sqrt{q^{k}\log q\cdot\frac{n}{\overline{m}}}\right)-quasirandom.

Orthonormality once again gives us Plancherel’s Theorem in this setting:

See [O’D14, Aus08] for more background on Fourier analysis over larger domains.

B.2 Conversion to Boolean functions

The degree of fbf^{b} is equal to the degree of ff.

by the assumption that v∈Ωkv\in\Omega_{k}. The degree of fbf^{b} is therefore ∣σ∣|\sigma|. ∎

B.3 Quasirandomness and strong refutation

To prove quasiandomness and strong refutation results for CSPs over larger alphabets, we proceed exactly as in the binary case. We used the t=kt=k case of Lemma 3.13 (the Vazirani XOR Lemma [Vaz86, Gol11]) to certify quasirandomness for binary CSPs. A generalization of this case holds for Abelian groups [Rao07, Lemma 4.2].

The proof is essentially identical to the proof of Lemma 4.3. We highlight the differences. First of all, we can write

These two lemmas then imply the larger alphabet versions of the quasirandomness certification and strong refutation results above.

There is an efficient algorithm that certifies that an instance I∼Fq,P(n,p)\mathcal{I}\sim\mathcal{F}_{q,P}(n,p) of CSP(P)(P) is γ\gamma-quasirandom with high probability when m‾≥qO(k)nk/2log⁡5nγ2\overline{m}\geq\frac{q^{O(k)}n^{k/2}\log^{5}n}{\gamma^{2}}.

B.4 Refutation of non-tt-wise supporting CSPs

The proof uses the following dual linear programs exactly as in the proof of Lemma 3.16.

The rest of the proof is exactly as in the binary case. ∎

We can again use these separating polynomials to obtain almost δ\delta-refutation for predicates that are δ\delta-far from tt-wise supporting.

The proof is essentially identical to Proof 2 of Theorem 4.9.

Corollary 4.11 also extends to larger alphabets.

This follows directly from Theorem B.9 and the following extension of Corollary 3.17 to larger alphabets.

The proof is essentially identical to the proof of Corollary 3.17: Observe that the LP (39) has at most qtktq^{t}k^{t} variables and proceed exactly as before.

B.5 SOS proofs

Here we give SOS versions of our refutation results for larger alphabets.

To give an SOS proof that Fourier coefficients of DI,ybD^{b}_{\mathcal{I},y} are small, we again need to define a specific polynomial representation of DI,yb^(σ)\widehat{D^{b}_{\mathcal{I},y}}(\sigma).

with high probability, assuming also that m‾≥max⁡{ns/2,n}\overline{m}\geq\max\{n^{s/2},n\}.

Given an instance I∼Fq,P(n,p)\mathcal{I}\sim\mathcal{F}_{q,P}(n,p) of CSP(P)(P),

with high probability when m‾≥qO(k)nk/2log⁡5nγ2\overline{m}\geq\frac{q^{O(k)}n^{k/2}\log^{5}n}{\gamma^{2}}.

First, use the Fourier expansion of PP to write

Summing over all TT, cc, and σ\sigma, we obtain the following.

with high probability when m‾≥qO(k)nt/2log⁡5nγ2\overline{m}\geq\frac{q^{O(k)}n^{t/2}\log^{5}n}{\gamma^{2}} and t≥2t\geq 2.

To prove this theorem, we need a version of Claim 6.7 for larger alphabets.

The first term is equal to fb(σ)f^{b}(\sigma). For the second term, note that each of the products ∏i∈[k]v(i,αi′)v(i,αi′′)\prod_{i\in[k]}v(i,\alpha^{\prime}_{i})v(i,\alpha^{\prime\prime}_{i}) must contain factors v(i,a)v(i,b)v(i,a)v(i,b) with a≠ba\neq b since α′≠α′′\alpha^{\prime}\neq\alpha^{\prime\prime}. We have the axiom v(i,a)v(i,b)=0v(i,a)v(i,b)=0, so the second term is 00. Then fb=(gb)2f^{b}=(g^{b})^{2} and the claim follows. ∎

With this claim, the proof of the theorem exactly follows that of Theorem 4.9.

Summing over all constraints, we get that

Just as in the proof of Theorem B.14, we can rewrite this in degree-kk SOS as

Since E\/[Q]=0\mathop{\bf E\/}[Q]=0 and Q≥−1Q\geq-1, we know that ∣Q∣≤qO(k)|Q|\leq q^{O(k)} and therefore ∣Q^(σ)∣≤qO(k)|\widehat{Q}(\sigma)|\leq q^{O(k)}. We can then apply Lemma B.12 for each σ\sigma to complete the proof. ∎

Appendix C Certifying that random hypergraphs have small independence number and large chromatic number

We define H(n,p,k)\mathcal{H}(n,p,k) to be the distribution over nn-vertex, kk-uniform (unordered) hypergraphs in which each of the (nk)\binom{n}{k} possible hyperedges is included independently with probability pp. Let m‾\overline{m} be the expected number of hyperedges p(nk)p\binom{n}{k}.

Coja-Oghlan, Goerdt, and Lanka used CSP refutation techniques to show the following results [COGL07]:

(Coja-Oghlan–Goerdt–Lanka [COGL07, Theorem 3]). For H∼H(n,p,3)H\sim\mathcal{H}(n,p,3), there is a polynomial time algorithm certifying that α(H)<ϵn\alpha(H)<\epsilon n with high probability for any constant ϵ>0\epsilon>0 when m‾>n3/2ln⁡6n\overline{m}>n^{3/2}\ln^{6}n and m‾=o(n2)\overline{m}=o(n^{2}).

(Coja-Oghlan–Goerdt–Lanka [COGL07, implicit in Section 4]). For H∼H(n,p,4)H\sim\mathcal{H}(n,p,4), there is a polynomial time algorithm certifying that α(H)<ϵn\alpha(H)<\epsilon n with high probability for any constant ϵ>0\epsilon>0 when m‾≥O(n2ϵ4)\overline{m}\geq O\left(\frac{n^{2}}{\epsilon^{4}}\right).

(Coja-Oghlan–Goerdt–Lanka [COGL07, Theorem 4]). For H∼H(n,p,4)H\sim\mathcal{H}(n,p,4), there is a polynomial time algorithm certifying that χ(H)>ξ\chi(H)>\xi with high probability for constant ξ\xi when m‾≥O(ξ4n2)\overline{m}\geq O(\xi^{4}n^{2}).

We generalize these results to kk-uniform hypergraphs:

For H∼H(n,p,k)H\sim\mathcal{H}(n,p,k), there is a polynomial time algorithm certifying that α(H)<β\alpha(H)<\beta with high probability when m‾≥Ok(n5k/2log⁡3nβ2k)\overline{m}\geq O_{k}\left(\frac{n^{5k/2}\log^{3}n}{\beta^{2k}}\right), assuming that β≥n3/4log⁡n\beta\geq n^{3/4}\log n.

For H∼H(n,p,k)H\sim\mathcal{H}(n,p,k), there is a polynomial time algorithm certifying that χ(H)>ξ\chi(H)>\xi with high probability when m‾≥Ok(ξ2knk/2log⁡3n)\overline{m}\geq O_{k}\left(\xi^{2k}n^{k/2}\log^{3}n\right), assuming that ξ≤n1/4log⁡n\xi\leq\frac{n^{1/4}}{\log n}.

The proofs are simple extensions of the k=3k=3 and k=4k=4 cases from [COGL07]. We will first prove Theorem C.4 using Theorem 4.1 and this will almost immediately imply Theorem C.5.

For T∈[n]kT\in[n]^{k}, we define the random variable w(T)w(T) as follows:

Let x∈{0,1}nx\in\{0,1\}^{n} be the indicator vector of an independent set II so that xT=1x^{T}=1 if T⊆IT\subseteq I and xT=0x^{T}=0 otherwise. First, observe that

where the second term is 00 because II is an independent set. The w(T)w(T)’s satisfy conditions (6), (7), and (8) and ∥x∥∞≤1\|x\|_{\infty}\leq 1, so Theorem 4.1 implies we can certify that

with high probability. Simplifying, we see that we can certify

and plugging in the value of m‾\overline{m} from the statement of the theorem completes the proof. ∎

For a coloring of a hypergraph HH, each color class is an independent set of HH. If χ(H)≤ξ\chi(H)\leq\xi, then there exists a color class of size at least nξ\frac{n}{\xi} and therefore α(H)≥nξ\alpha(H)\geq\frac{n}{\xi}. We can then certify that α(H)<nξ\alpha(H)<\frac{n}{\xi} using Theorem C.4. ∎

Appendix D Simulating ℱP​(n,p)\mathcal{F}_{P}(n,p) with a fixed number of constraints

The setting of [DLSS14] fixes the number of constraints in a CSP instance, whereas the model described in Section 3 includes each possible constraint in the instance with some probability pp. Here we show that results from our setting easily extend to that of [DLSS14] by giving an algorithm that simulates the behavior of our model when the number of constraints is fixed.

Recall that an instance I∼FP(n,p)\mathcal{I}\sim\mathcal{F}_{P}(n,p) is generated as follows. For each S∈[n]kS\in[n]^{k} and each c∈{−1,1}kc\in\{-1,1\}^{k}, constraint (c,S)(c,S) is included with probability pp, so the expected number of constraints is p⋅(2n)kp\cdot(2n)^{k}.

In the model where the number of constraints is fixed, the instance is guaranteed to have mm distinct constraints for some value of mm. The instance J\mathcal{J} is chosen uniformly from all subsets of {−1,1}k×[n]k\{-1,1\}^{k}\times[n]^{k} with size exactly mm

On a random instance with μ\mu constraints, we can generate an instance I\mathcal{I} that simulates this behavior by choosing an appropriate value for pp, drawing m∼Binomial(p,(2n)k)m\sim\textrm{Binomial}\left(p,(2n)^{k}\right) and then discarding μ−m\mu-m of the constraints. For brevity, let d=(μ−1ln⁡μ)1/2.d=\left(\mu^{-1}\ln\mu\right)^{1/2}. Algorithm 1 describes the behavior of A\mathcal{A}.

Furthermore, the probability of failing to refute an instance with value at most 1−η+2d1-\eta+2d due to exiting at step 2 is ok,t(1)o_{k,t}(1). We treat mm as a sum of (2n)k(2n)^{k} independent Bernoulli variables with probability pp and denote E\/[m]\mathop{\bf E\/}[m] by m‾\overline{m}. Applying a Chernoff bound yields the following.