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 -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 -SAT, . The model here involves choosing a CNF formula over variables by drawing clauses (ORs of 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 -complete and -complete problems. The computational hardness depends crucially on the density, . For each there is (conjecturally) a constant critical density such that is satisfiable with high probability when , and is unsatisfiable with high probability when . (Here and throughout, “with high probability (whp)” means with probability as .) This phenomenon occurs for all nontrivial random CSPs; in the case of -SAT it’s been rigorously proven [DSS15] for sufficiently large .
There is a natural algorithmic task associated with each of the two regimes. When one wants to find a satisfying assignment for . When one wants to refute ; 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 is any algorithm that: a) outputs “unsatisfiable” or “fail”; b) never incorrectly outputs “unsatisfiable”; c) outputs “fail” with low probability (i.e., probability ). 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 increases towards , finding satisfying assignments becomes more difficult; and conversely, as increases beyond , finding certificates of unsatisfiability gradually becomes easier.
The special case of random -XOR has proved particularly important: it is related to -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 for simplicity.
The proof of Theorem 2.1 follows ideas from [COGL07] and earlier works on “discrepancy” of random -SAT instances. The case of even is notably easier, and we present two “folklore” arguments for it. The case of odd is trickier. Roughly speaking we view the instance as a homogeneous degree- multilinear polynomial, which we want to certify takes on only small values on inputs in . Considering separately the contributions based on the “last” of the variables in each constraint, and then using Cauchy–Schwarz, it suffices to bound the norm of a carefully designed quadratic form of dimension , indexed by tuples of 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 -SAT given constraints.
With Theorem 2.1 in hand, the next step is certifying quasirandomness of random -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 , the induced -tuples of literal values are close to being uniformly distributed over . (Note that this is only a property of the instances’ constraint scopes/negations, and has nothing to do with .) Since the “Vazirani XOR Lemma” implies that a distribution on is uniform if and only if all its XORs are have bias , 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:
In particular, this theorem improves upon [COCF10] by a factor of whenever is odd; this savings is new even in the well-studied case of -SAT.
We remark that property of a predicate 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 for which a predicate supports a -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 for which the distribution over constraint inputs induced by the planted assignment is -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 out of coordinates. Thus for every assignment , the induced distribution on constraint -tuples is (-close to) -wise uniform. Since does not support a -wise uniform distribution, this essentially shows that no can induce a fully-satisfying distribution on constraint inputs. To handle the -closeness caveat, we show that if does not support a -wise uniform distribution, then it is -far from supporting such a distribution, for . The algorithm can then in fact -refute random CSP instances.
To briefly illustrate the result, consider the Exactly--out-of--SAT CSP, studied previously in [BB02, GJ03]. The associated predicate supports a -wise uniform distribution, namely the uniform distribution over strings in of Hamming weight . 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 .
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 are arity- predicates introduced in [Hua13]; they are hereditarily approximation resistant on satisfiable instances and have -variability . In [DLSS14] they are used in SRCSP Assumption 1 to deduce hardness of PAC-learning DNFs with terms. However:
For all , the predicate does not support a -wise uniform distribution.
Thus by Theorem 2.4 we can efficiently refute random instances of CSP() with just \mathchoice{\hbox{\displaystyle\widetilde{O}}}{\hbox{\textstyle\widetilde{O}}}{\hbox{\scriptstyle\widetilde{O}}}{\hbox{\scriptscriptstyle\widetilde{O}}}(n^{2}) constraints, independent of . This contradicts SRCSP Assumption 1.
Finally, we also prove that SRCSP Assumption 1 is false for another family of predicates used by [DLSS14] to show hardness of PAC-learning intersections of Boolean halfspaces.
Our results described in these three cases all use linear programming duality. Specifically, in Lemma 3.16 we show that is -far from supporting a -wise uniform distribution if and only if there exists a -variable multilinear polynomial that satisfies certain properties involving and . We then explicitly construct these dual polynomials for the Huang, Majority, and 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 terms under a much weaker assumption than SRCSP Assumption . Specifically, their work only assumes that for all there is a large enough such that refuting random -SAT instances is hard when there are constraints. This assumption looks quite plausible to us, and may even be true with not much larger than . 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 -XOR instances is hard when for some . He also shows hardness of learning halfspaces to within an approximation factor of for any assuming that there exists some constant such that for all , refuting random -XOR instances with is hard when .
2 Evidence for the optimality of our results
It’s natural to ask whether the 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 -wise supporting predicates is hard when . A natural form of evidence would be showing that various strong classes of polynomial-time refutation algorithms fail when . 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 -variable “degree- SOS proofs” can be found in time whenever they exist; c) we do know some examples of lower bounds for degree- SOS proofs. In Section 6 of this paper we show the following:
All of our refutation algorithms for -ary predicates can be extended to produce degree- 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 -SAT instances, especially in the Resolution proof system. This culminated in the work of Ben-Sasson and Wigderson [BSW99], which showed that for random -SAT (and -XOR) with , Resolution refutations require size (whp). More recently, Schoenebeck [Sch08] showed (using the expansion analysis of [BSW99]) that random -XOR and -SAT instances with require SOS proofs of degree , and hence take 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 -XOR (which also works in -degree SOS) and our bound for random -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 -SAT instances with .
We now discuss our main theorem’s bound for predicates not supporting -wise uniform distributions. Suppose is a predicate that does support a -wise uniform distribution, where . In the context of inapproximability and SDP-hierarchy integrality gaps, this condition on has been significantly studied in the case of . For supporting pairwise uniformity, it is known [BGMT12, TW13] that the Sherali–Adams and Lovász–Schrijver+ SDP hierarchies require degree to refute random CSP instances (whp) when . This result was also recently proven for the stronger SOS proof system by Barak, Chan, and Kothari [BCK15], except that the CSP instances are not quite uniformly random; they are “slightly pruned” random instances. For any , the second and third authors recently essentially showed [OW14] that for the Sherali–Adams+ SDP hierarchy, degree is (whp) necessary to refute random CSP instances when . 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--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 - and -uniform hypergraphs have small independence number and large chromatic number. We extend these results to random -uniform hypergraphs.
For a random -uniform hypergraph , there is a polynomial time algorithm certifying that the independence number of is at most with high probability when 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 -uniform hypergraph , there is a polynomial time algorithm certifying that the chromatic number of is at least with high probability when 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 , let be the distribution over CSP instances given by including each of the possible constraints independently with probability . 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 , 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 is not a symmetric function. We use to denote the expected number of constraints, namely . As noted in Fact 3.6 below, the number of constraints in a draw from is very tightly concentrated around , and we often blur the distinction between these parameters. Appendix D explicitly describes a method for simulating an instance drawn from 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 for the uniform distribution on as well as the following:
An immediate consequence of an instance being quasirandom is that its optimum is close to :
We conclude the discussion of CSPs by recording some facts that are proven easily with the Chernoff bound:
Let . Then the following statements hold with high probability.
.
is -quasirandom.
2 Algorithms and refutations on random CSPs
where is defined by . Although 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 .
3 tt-wise uniformity
An important notion for this paper is that of -wise uniformity. Recall:
Probability distribution on is said to be -wise uniform, , if for all with the random variable is uniform on when . (We remark that this condition is sometimes inaccurately called “-wise independence” in the literature.)
We will also consider the more general notion of -wise uniformity. This is typically defined using Fourier coefficients:
Probability distribution on is said to be -wise uniform if for all with , where is the probability density associated with distribution .
It is a simple fact (and it follows from Lemma 3.13 below) that -wise uniformity is equivalent to -wise uniformity.
Also important for us is a related but distinct notion, that of being -close to a -wise uniform distribution. It’s easy to show that if is -close to a -wise uniform distribution then is -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 is an -wise uniform distribution on , then there exists a -wise uniform distribution on with
In particular if we have the bound (and this can also be improved [Gol11] to ).
A predicate is said to be -wise supporting if there is a -wise uniform distribution whose support is contained in . We say is -far from -wise supporting if every -wise uniform distribution is -far from being supported on ; i.e., has probability mass at least on .
4 A dual characterization of limited uniformity
It is known that the condition of supporting a -wise uniform distribution is equivalent to the feasibility of a certain linear program; hence one can show that is not -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 -wise supporting.
;
;
, i.e., has no constant coefficient.
We now provide the quantitative version of the aforementioned [AH09, Theorem 3.1]:
Let and let . Then is -far from -wise supporting if and only if there is a -separating polynomial for of degree at most .
The proof is an application of linear programming duality. Consider the following LP, which has variables for each .
minimize (1) s.t. (2) (3)
Constraint (3) and the nonnegativity constraint ensure that is a probability distribution on . Constraint (2) expresses that is -wise uniform (see Remark 3.12). The objective (1) is minimizing the probability mass that puts on assignments in . Thus the optimal value of the LP is equal to the smallest such that is -close to -wise supporting; equivalently, the largest such that is -far from -wise supporting.
The following is the dual of the above LP. It has a variable for each as well as a variable corresponding to constraint (3).
maximize (4) s.t. (5)
Observe that a feasible solution is precisely equivalent to a multilinear polynomial of degree at most , namely , that -separates .
Thus is -far from -wise supporting if and only if the LP’s objective (1) is at least , if and only if the dual’s objective (4) is at least , if and only if there is a -separating polynomial for of degree at most . ∎
From this proof we can also derive that if fails to be -wise supporting then it must in fact be -far from being -wise supporting:
Suppose is not -wise supporting. Then it is in fact -far from -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 ).
Let be the number of variables in the dual LP from Lemma 3.16, so in general, with when . 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 matrix of entries and a “right-hand side” vector with entries. By Cramer’s rule, the solution’s entries are ratios of determinants of integer matrices with entries in . Thus any strictly positive entry is at least , where is the maximum possible such determinant. By Hadamard’s inequality, and the claimed result follows. ∎
Quasirandomness and its implications for refutation
For and , let be independent random variables such that for each ,
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 -XOR as a simple corollary.
so for a -XOR instance \mathcal{I}\sim\mathcal{F}_{\textup{k-XOR}}(n,p),
where . The ’s are random variables depending on the choice of ; observe that , , and for all . By Theorem 4.1, there is an algorithm certifying that
with high probability when . Since with high probability, choosing 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 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 has small magnitude.
Let with . There is an algorithm that, with high probability, certifies that
for all , assuming also that .
To prove this lemma, we need another lemma certifying that polynomials whose coefficients are sums of -mean random variables have small value.
The proof is straightforward and we defer it to Section 4.4.
Without loss of generality, assume . Applying definitions, we see that
where we define and recall that is the projection of onto the the coordinates in . It is clear that , , and . There are terms in each sum of ’s and we can apply Lemma 4.4. When , we plug in these values and see that we can certify that . When , implies that and we can certify that . The lower bound can be proved in exactly the same way by considering the random variables . ∎
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 -wise supporting, a weaker notion of quasirandomness suffices to obtain -refutation.
An instance of is -quasirandom if is -wise uniform for every .
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 -quasirandom for all with high probability. Lemma 4.3 directly implies that we can certify -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 of CSP is -quasirandom with high probability when and .
We now reach the main result of this section, which states that if a predicate is -far from -wise supporting, then we can almost -refute instances of CSP.
We give two proofs of this theorem. In Proof 1, the theorem follows directly from certification of -quasirandomness and Lemma 3.13.
Proof 2 gives a slightly weaker version of Theorem 4.9, requiring the stronger assumption that . It is based on the dual polynomial characterization of being -far from -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 -wise uniformity and -closeness to a -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 is -far from -wise supporting, there exists a degree- polynomial that -separates . The definition of -separating implies that for all . Summing over all constraints, we get that for all ,
It then remains to certify that . Observe that
where the second equality follows from Plancherel’s Theorem. Since and , and hence for all . To finish the proof, we apply Lemma 4.3 to certify that for all . ∎
4 Proof of Lemma 4.4
Let be independent -mean random variables such that . Then, for ,
Observe that the ’s are independent and that each one is the sum of mean-, i.i.d. random variables with magnitude at most . Noting that , we can use Bernstein’s Inequality to show that the ’s are not too big with high probability. If , Theorem 4.1 then implies that the desired algorithm exists. If , we are simply bounding a linear function over variables. We consider two cases: Small and large .
Choosing in Bernstein’s Inequality, we see that . A union bound over all then implies that \mathop{\bf Pr\/}[\text{any|v_{U}|>2s\log n}]\leq n^{-s}. If , we observe that , scale the ’s down by , and apply Theorem 4.1 to get the stated result. If , we obtain the second bound by observing that
We set 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 , we can then divide the ’s by and apply Theorem 4.1. If , we get a bound of in the same way as (10). ∎
Hardness of learning implications
Daniely et al. reduce the problem of distinguishing between random instances of and instances with value at least 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 that is specific to the hypothesis class and assume hardness of distinguishing between random instances of and instances with constraints and value at least for all . 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 ). They also show that if the given CSP instance is random, the set of examples will have error at least (in the agnostic case for all 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 from satisfiable instances, even when the algorithm is allowed to err with probability 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 . 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 and corresponding , we define a degree- polynomial that -separates . 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 from those that are satisfiable (or have value at least ). In order to simplify the presentation, we begin with simpler versions of the polynomials and then scale them to attain the appropriate values of . The following lemma will be of use for this scaling.
Define . Clearly is also unbiased and has degree . Then for all , . Similarly, for all , . ∎
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 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 -variability is [DLSS14].
Daniely et al. reduce the problem of distinguishing between random instances of with constraints and satisfiable instances to the problem of improperly PAC learning the class of DNF formulas with terms on a sample of training examples with error with probability at least . Here we show that there exists a polynomial time algorithm that refutes random instances of by demonstrating that does not support a 4-wise uniform distribution and applying Theorem 4.9.
As a notational shorthand, write for . Define as follows:
Observe that for each monomial of , for every , . Further, for each with , appears exactly once in . Let be the set of all ordered 6-tuples of distinct elements of . For an ordered tuple , we use to denote membership in .
Observe that does not depend on any of . By construction, contains no constant term, so . Clearly for all because (11) is always at least .
Now we lower bound the value of on all that satisfy . We first show that for any that strongly satisfies the Huang predicate, , then bound for any with Hamming distance at most from . By definition, for each , we have that . So for each monomial of ,
where the last line follows from the fact that Because there are monomials in , their sum is .
Now we consider the case where does not strongly satisfy the Huang Predicate, but . Any singleton index on which and differ will not change the value of . Let . We lower bound by counting the number of monomials in which each appears and
For fixed , the number of monomials containing the variables of is
because there are exactly ways to permute the three indices of in I and the remaining indices are permuted in the remaining 3 positions of . So
If we instead choose to scale by a factor of rather than substituting into (12), we can achieve a better separation of . For , this expression is strictly increasing and it approaches as 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 we present a multivariate polynomial that -separates for some . 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 -separates its predicate for the appropriate value of . We give the construction below.
For where , define and call the Hamming weight of .
Note that this is analogous to the notion of a Hamming weight of a vector in , but differs in that it is not simply the count of the number of ’s. We define a general predicate that is satisfied when is at least some fixed threshold value .
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 for all .
where denotes the Krawtchouk polynomial of degree [Kra29, KL96]. Substituting , 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 and , we can work with a univariate polynomial while ensuring that its multivariate analogue is unbiased and has degree at most 4 (degree 3 when ).
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 , . It therefore suffices to lower bound both when and for all .
First we show that is monotonically increasing in .
which is evidently positive for .
Because is monotonically increasing in , for all .
which is clearly negative for . Now it just remains to lower-bound for . Again, since is monotonically increasing in , we use the value :
Because is evidently and for all , we have the following Corollary.
For odd and sufficiently large , there exists an efficient algorithm that distinguishes between random instances of 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 , . To simplify , let . Then we can rewrite (18) as follows:
where the last inequality follows from the fact that and the first two terms are always nonnegative.
Next we lower-bound for .
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 -wise supporting, one may attempt to find another threshold-based predicate. Here we show that a symmetric threshold predicate that is -wise supporting must be satisfied by all strings with Hamming weight at least . Furthermore, there exists a symmetric threshold predicate that is -wise supporting with a threshold of and we sketch its construction.
Again, for simplicity we set and obtain the following expression:
Observe that for , . We now lower-bound the value of for , or equivalently, :
The first three terms are always nonnegative, so .
Applying Lemma 5.1, is -far from supporting a -wise uniform distribution. ∎
Assume for some integer . Then there exists a -wise uniform distribution supported only on such that .
Let be a binary BCH code of length with designed distance and let be its dual. Then the uniform distribution on the codewords of is -wise uniform [ABI86, MS77]; see also [AS04, Ch 16.2].
Let be a codeword of where each . The Carlitz-Uchiyama bound [MS77, page 280] states that for all ,
Observe that the quantity simply maps from to so that we can write the bound to match the presentation in [MS77]. Therefore,
Setting , we can obtain -wise uniformity on this distribution and each string in the support of the distribution has Hamming weight at least ∎
In order to construct a -wise uniform distribution for any value of , one could simply express 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 , after which point we use the uniform distribution, and obtain a -wise uniform distribution. The total Hamming weight of a vector supported by this distribution would then be at least .
SOS refutation proofs
with , for all , and for all . If it also holds that , we will write .
It is well-known that a degree- SOS proof can be found using an SDP of size if it exists [Sho87, Par00, Las00, Las01].
In this section, we will take the set to be , enforcing that variables are -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 and , let be independent random variables such that for each ,
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 , the lemma follows by using Theorem 6.2 instead of Theorem 4.1. If , it suffices to show that
for any since summing over all as in (10) finishes the proof. If , observe that
If , we use instead of . ∎
The lemma implies an SOS version of Lemma 4.3. To make this precise, we define a specific polynomial representation of :
where . Note that this is a polynomial in the ’s.
We can show these polynomials are not too large.
Let with . Then
with high probability, assuming also that .
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().
with high probability when .
Note that this is just Plancherel’s Theorem in SOS. The theorem then follows from Lemma 6.4 and the observation that . ∎
4 Ω(1)\Omega(1)-refutation of non-tt-wise supporting CSPs
with high probability when and .
The proof is an SOS version of Proof 2 of Theorem 4.9 above. Claim 6.7 implies that for of degree at most that -separates ,
Rearranging terms as in the proof of Theorem 6.5, we see that the right hand side is equal to
Since has mean , and . 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 .
Directions for future work
Additionally, it would be good to investigate whether our refutation algorithms can be extended from the purely random CSP setting to the “smoothed”/“semi-random” setting of Feige [Fei07], in which the constraints scopes are worst-case and only the negation pattern for literals is random. Feige showed how to efficiently refute random -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 -SAT when (as well as a subexponential-time deterministic algorithm). This raises the question of whether there exist polynomial-size refutations for random CSP 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 -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 and , let be independent random variables such that for each ,
Then there is an efficient algorithm certifying that
When is even, we can think of as a quadratic form:
where . We give two methods to certify that the value of this quadratic form is at most . The first method uses an SDP-based approximation algorithm and works only for . The second method uses ideas from random matrix theory and works for any with .
If , we can apply an approximation algorithm of Charikar and Wirth [CW04] for quadratic programming. They prove the following theorem:
[CW04, Theorem 1] Let be any matrix with all diagonal elements . There exists an efficient randomized algorithm that finds such that
By Markov’s Inequality, this statement holds with probability at least . We can run the algorithm 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 . For the first term, we will need the following claim.
This follows from applying Bernstein’s Inequality (Theorem 4.12) for fixed and then taking a union bound over all . 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 .
We will use the next claim to bound the second term of (29).
Since the and , the claim follows from the Chernoff Bound. The second term of (29) is upper bounded by and we can compute this quantity in polynomial time to certify that its value is at most .
Observe that (28) is for a matrix indexed by so that . Then . To certify that is small, we compute . We need to show that is small with high probability. First, note that is equal to the norm of the symmetric matrix
[Tao12, Proposition 2.3.13] Let be a random symmetric matrix whose upper triangular entries with are independent random variables with mean , variance at most , and magnitude at most . Then, with high probability,
A.2 The odd arity case
Fix an assignment . For , the monomials containing can contribute at most to the objective if is set optimally. By Cauchy-Schwarz,
so it suffices to bound . We will write this as a quadratic polynomial and then bound it using spectral methods:
Define the matrix indexed by :
The first term is at most since the variables are bounded. We can compute to certify this. With high probability, is not too big.
Let and . Let be indepedent random variables satisfying conditions (6), (7), and (8) above. Let be defined as in (32). With high probability,
We can therefore certify that the first term is . We will prove the lemma in Appendix A.4.
The second term of (33) is at most . We can easily compute this and the Chernoff Bound implies that its value is at most with high probability.
So far, with high probability we can certify that . Plugging this bound into (30) concludes the proof.
It would have been more natural to have written for such that . However, could be too large because of the contribution of the second term in (33). We use the additional assumption that to get around this issue.
A.3 An SOS version
In this section, we will prove the SOS version of Theorem 4.1.
For and , let be independent random variables such that for each ,
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 . When , there exists a matrix such that . 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 and the statement of the lemma.
Let and . Let be indepedent random variables satisfying conditions (6), (7), and (8) above. Let be the indexed by that is defined as follows:
Recall that we index by elements of divided into two blocks of coordinates each. First, note that
Expanding this out using the definition of and setting , we get that
To bound this sum, we will start by bounding . We will need two claims.
The two claims also imply two other facts we will need below.
If , then .
If , then .
This can be proved in exactly the same manner.
Next, observe that the number of choices of with is at most . The number of choices of with is at most . All together, we can write
Recall that we assumed . Since and , the claim follows. ∎
If we did not have conditions (35), (36), and (37), we would only have been able to show that . This would have led to a weaker bound of .
Appendix B Extension to larger alphabets
Let . Then the following statements hold with high probability.
.
is -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 is equal to the degree of .
by the assumption that . The degree of is therefore . ∎
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 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 of CSP is -quasirandom with high probability when .
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 -refutation for predicates that are -far from -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 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 are small, we again need to define a specific polynomial representation of .
with high probability, assuming also that .
Given an instance of CSP,
with high probability when .
First, use the Fourier expansion of to write
Summing over all , , and , we obtain the following.
with high probability when and .
To prove this theorem, we need a version of Claim 6.7 for larger alphabets.
The first term is equal to . For the second term, note that each of the products must contain factors with since . We have the axiom , so the second term is . Then 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- SOS as
Since and , we know that and therefore . We can then apply Lemma B.12 for each to complete the proof. ∎
Appendix C Certifying that random hypergraphs have small independence number and large chromatic number
We define to be the distribution over -vertex, -uniform (unordered) hypergraphs in which each of the possible hyperedges is included independently with probability . Let be the expected number of hyperedges .
Coja-Oghlan, Goerdt, and Lanka used CSP refutation techniques to show the following results [COGL07]:
(Coja-Oghlan–Goerdt–Lanka [COGL07, Theorem 3]). For , there is a polynomial time algorithm certifying that with high probability for any constant when and .
(Coja-Oghlan–Goerdt–Lanka [COGL07, implicit in Section 4]). For , there is a polynomial time algorithm certifying that with high probability for any constant when .
(Coja-Oghlan–Goerdt–Lanka [COGL07, Theorem 4]). For , there is a polynomial time algorithm certifying that with high probability for constant when .
We generalize these results to -uniform hypergraphs:
For , there is a polynomial time algorithm certifying that with high probability when , assuming that .
For , there is a polynomial time algorithm certifying that with high probability when , assuming that .
The proofs are simple extensions of the and cases from [COGL07]. We will first prove Theorem C.4 using Theorem 4.1 and this will almost immediately imply Theorem C.5.
For , we define the random variable as follows:
Let be the indicator vector of an independent set so that if and otherwise. First, observe that
where the second term is because is an independent set. The ’s satisfy conditions (6), (7), and (8) and , 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 from the statement of the theorem completes the proof. ∎
For a coloring of a hypergraph , each color class is an independent set of . If , then there exists a color class of size at least and therefore . We can then certify that 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 . 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 is generated as follows. For each and each , constraint is included with probability , so the expected number of constraints is .
In the model where the number of constraints is fixed, the instance is guaranteed to have distinct constraints for some value of . The instance is chosen uniformly from all subsets of with size exactly
On a random instance with constraints, we can generate an instance that simulates this behavior by choosing an appropriate value for , drawing and then discarding of the constraints. For brevity, let Algorithm 1 describes the behavior of .
Furthermore, the probability of failing to refute an instance with value at most due to exiting at step 2 is . We treat as a sum of independent Bernoulli variables with probability and denote by . Applying a Chernoff bound yields the following.