Spectral Algorithms for Unique Games

Alexandra Kolla

Introduction

A Unique Game is defined in terms of a constraint graph G=(V,E)G=(V,E), an integer kk which is called the alphabet, a set of variables {xu}uV\{x_{u}\}_{u\in V}, one for each vertex uu, and a set of permutations (constraints) πuv:[k][k]\pi_{uv}:[k]\rightarrow[k], one for each edge (u,v)E(u,v)\in E. An assignment to the variables (labeling) is said to satisfy the constraint on the edge (u,v)E(u,v)\in E if πuv(xu)=xv\pi_{uv}(x_{u})=x_{v}. The edges are taken to be undirected and hence πuv=(πvu)1\pi_{uv}=(\pi_{vu})^{-1}. The goal is to assign a value from the set [k][k] to each variable xux_{u} so as to maximize the number of satisfied constraints.

Khot [Kho02] conjectured that it is NP-hard to distinguish between the cases when almost all the constraints of a Unique Game are satisfiable and when very few of the constraints are satisfiable:

(Unique Games Conjecture-UGC) For any constants ϵ,δ>0\epsilon,\delta>0, there is a k(ϵ,δ)k(\epsilon,\delta) such that for any k>k(ϵ,δ)k>k(\epsilon,\delta), it is NP-hard to distinguish between instances of Unique Games with alphabet size kk where at least 1ϵ1-\epsilon fraction of constraints are satisfiable and those where at most δ\delta fraction of constraints are satisfiable.

There are two reasons which make UGC particularly intriguing. First, it is a well-balanced question. Despite continuous efforts to prove or disprove it, there is still no consensus regarding its validity. This seems to indicate that UGC is more likely to be resolved in the near future in contrast to the P vs. NP problem for example, for which it is widely believed that P\neqNP but current techniques have been unable to prove it. Second, as seen by a series of works, the truth of UGC implies that the currently best known approximation algorithms for many important computational problems have optimal approximation ratios. Since its origin, UGC has been successfully used to prove often optimal hardness of approximation results for several important NP-hard problems such as Min-2Sat-Deletion [Kho02], Vertex Cover [KR03], Maximum Cut [KKMO04] and non-uniform Sparsest Cut [CKK+06, KV05]. In addition, in recent years, UGC has also proved to be intimately connected to the limitations of semidefinite programming (SDP). Making this connection precise, the authors in [Aus10] and [Rag08] show that if UGC is true, then for every constraint satisfaction problem (CSP) the best approximation ratio is given by a certain simple SDP.

Arguably, a seemingly strong reason for belief in UGC is the failure of several attempts to design efficient algorithms for Unique Games using current state-of-the-art techniques, such as linear and semidefinite programming (LP/SDP). Indeed, several works [KV05, RS09] show limitations of LPs and SDPs in solving Unique Games by exhibiting the existence of itegrality gap instances. Those are Unique Games instances that fool the SDP since, even though they have no good satisfying assignment, the respective SDP solution is high. The existence of such instances implies, in particular, that there is no hope for a “good” approximation algorithm based on those SDPs. Moreover, recently it was shown that solving Unique Games is at least as hard as another seemingly hard problem called the small set expansion problem [RS10]. Our work presents evidence that a different set of techniques might be more powerful when it comes to designing algorithms for Unique Games thus giving hope that such techniques can be used to potentially disprove UGC.

We present a purely spectral algorithm for Unique Games that finds highly satisfying assignments when they exist. The running time of our algorithm depends on spectral properties of the Label-Extended graph or the constraint graph associated with the instance of Unique Games. Our algorithm runs in polynomial or quasi-polynomial time for a large class of instances, including those where the standard SDP provably fails as well as instances where the constraint graph is an expander. At a high level, we show that given ϵ>0\epsilon>0, there is a δ=δ(ϵ)\delta=\delta(\epsilon) such that the algorithm is able to distinguish between the following two cases:

YES case: There exists an assignment that satisfies (1ϵ)(1-\epsilon) fraction of the constraints.

NO case: Every assignment satisfies less that δ\delta fraction of the constraints.

In particular, our algorithm runs in quasi-polynomial time on input the SDP integrality gap instance as it appears in [KV05] (hereafter referred to as KV\mathcal{K}\mathcal{V}). We note that the authors in [KV05] roughly showed that, when run on a certain highly unsatisfiable instance of Unique Games, the standard SDP relaxation has very high objective value and consequently fails to distinguish between the two cases above. As another special case, our algorithm runs in polynomial time when the constraint graph is an expander (therefore re-deriving results similar to [AKK+08] and [KT07]). Moreover, similarly to [AKK+08, MM09, KT07, AIMS09], the performance of the algorithm does not depend on the alphabet size kk.

(Main) Let U=(G,M,k)\mathcal{U}=(G,M,k) be a (1ϵ)(1-\epsilon) satisfiable instance of Unique Games on a dd-regular graph GG with alphabet size kk. Let MM be the adjacency matrix of the label-extended graph of GG. Let WW be the space spanned by eigenvectors of MM that have eigenvalue greater than (1γ)d(1-\gamma)d, for γ>8ϵ\gamma>8\epsilon. There is an algorithm that runs in time 2O(γϵdim(W))+poly(nk)2^{O(\frac{\gamma}{\epsilon}\text{dim}(W))}+\text{poly}(n\cdot k) and finds an assignment that satisfies at least (1O(ϵγ8ϵ+ϵ))(1-O(\frac{\epsilon}{\gamma-8\epsilon}+\epsilon)) fraction of the constraints.

We also show a similar theorem for the special case where all the constraints are of the form Γ\Gamma-Max-Lin and the graph satisfies some special properties. That theorem is then used to derive the above-mentioned algorithm for expanders.

Our key contribution is a technique of effectively using the full spectrum and eigenvectors of the relevant graph rather than just the first few eigenvalues. Interestingly, prior to this work, researchers had often attempted to take advantage of the full spectrum of a graph in the design of approximation algorithms, but no significant progress was made. We believe that our techniques are of independent interest and could contribute to developing better algorithms for a number of other problems. In a recent breakthrough paper, Arora, Barak and Steurer [ABS10] used our main technique as a crucial building block in their sub-exponential time algorithm for arbitrary instances of Unique Games. Consequently, we believe that our results give evidence that spectral techniques might be a more powerful tool than SDPs for attacking and, potentially even disproving, the Unique Games Conjecture.

Several polynomial time approximation algorithms using linear and semidefinite programming have been developed for Unique Games on arbitrary graphs (see [Kho02], [Tre05], [GT06], [CMM06a], [CMM06b]). These algorithms start with an instance where the value of the SDP or LP relaxation is 1ϵ1-\epsilon and round it to a solution with value ν\nu. Here, value of the game refers to the maximum fraction of satisfiable constraints. For ν>δ\nu>\delta, this would give an algorithm to distinguish between the two cases. However, most of these algorithms give good approximations only when ϵ\epsilon is very small (ϵ=O(1/logn)\epsilon=O(1/\log n) or ϵ=O(1/logk)\epsilon=O(1/\log k)) and their approximation guarantee depends on the alphabet size kk It might be good to think of kk as O(logn)O(\log n) since this is the range of interest for most reductions.. For constant ϵ\epsilon, only the algorithm in [CMM06a] gives interesting parameters with νkϵ/(2ϵ)\nu\approx k^{-\epsilon/(2-\epsilon)}. We refer the reader to [CMM06a] for a comparison of parameters of various algorithms.

For some special cases of graphs, it has been shown that there are efficient algorithms that solve Unique Games. One such example is when the constraint graph is a spectral expander. In that case, the authors in [AKK+08] (and later [MM09] with improved parameters) showed that one can find a highly satisfying assignment (with, say, ν90%\nu\geq 90\%) in polynomial time. Notably, the approximation guarantee of their algorithm depends on the expansion parameters of the graph rather than the alphabet size. As another example, the authors in [AIMS09] presented efficient algorithms for constraint graphs with large local expansion.

Compared to the previous papers, our algorithm uses purely spectral methods as opposed to linear or semidefinite programming and runs in time that depends on spectral properties of the constraint graph or the label-extended graph. The approximation guarantee does not depend on the size of the graph or the size of the alphabet and thus our algorithm always finds a good assignment when one exists. We note that, in an unpublished manuscript, Kolla and Tulsiani [KT07] developed a spectral algorithm that runs in polynomial time and finds highly satisfying assignments on expander constraint graphs.

2 Overview of Our Techniques: Recovering Solutions by Spectral Methods

Our basic approach for the Unique Games algorithm is exhaustive search in a subspace spanned by several eigenvectors of a graph associated with the Unique Games instance. We identify a “good” subspace which contains a vector that “encodes” a highly satisfying assignment and then exhaustively search for this vector in (a discretization of) the subspace as we explain below.

Our algorithm simply looks at a set SW\mathcal{S}\subseteq W of appropriately many candidate vectors and “reads-off” an assignment.

Recover-SolutionS(U)\texttt{Recover-Solution}_{\mathcal{S}}(\mathcal{U})

For each xSx\in\mathcal{S}, construct an assignment LxL_{x} by assigning to each vertex uu, the index corresponding to the largest entry in the block (xu1,,xuk)(x_{u1},\ldots,x_{uk}) i.e. Lx(u)=argmaxixuiL_{x}(u)=\arg\max_{i}{x_{ui}}.

Out of all assignments LxL_{x} for xSx\in\mathcal{S}, choose the one satisfying the maximum number of constraints.

To choose S\mathcal{S},we identify a set of nice vectors NW\mathcal{N}\subseteq W such that the above algorithm works for any vector vv close to some vector in N\mathcal{N}. These are going to be the vectors that are close to the assignment eigenvectors. We then construct a set SW\mathcal{S}\subseteq W of test vectors such that at least one vector vSv\in\mathcal{S} is close to a vector in N\mathcal{N}. S\mathcal{S} is simply an epsilon-net for WW. Lastly, we go over all vectors in S\mathcal{S} until we find vv. The running time of the algorithm will depend on the size of S\mathcal{S} which, for an appropriately defined epsilon-net, it is exponential in the dimension of WW.

Certain “simple” graphs, including expander graphs and the Khot-Vishnoi instance, have only a few eigenvalues close to dd and thus the dimension of WW is small. Consequently, exhaustive enumeration in the subspace spanned by the corresponding eigenvalues will quickly find a good-enough assignment.

3 Organization of the Paper

The rest of the paper is organized as follows: in section 2 we give some preliminary notation and definitions that will be used in the rest of the paper. Section 3 contains the proofs of the main theorem above as well as of a theorem for the Γ\Gamma-Max-Lin case. More specifically, we prove the main theorem in subsection 3.1. The proof of the theorem for the Γ\Gamma-Max-Lin case appears in subsection 3.2. Subsection 3.3 contains the generalization of the main theorem for non-regular graphs. In section 4, we analyze our algorithm when run on the Khot-Vishnoi instance.

Preliminaries

We remind the reader that for a graph GG, the adjacency matrix A=AGA=A_{G} is defined as

For a graph GG, let DD be the diagonal matrix with diagonal entry D(u,u)=v:(u,v)EwuvD(u,u)=\sum_{v:(u,v)\in E}w_{uv} equal to the degree of node uu, namely the sum of the weights of edges adjacent to uu. The Laplacian of GG is defined as follows:

2 The Label-Extended Graph

For a given instance of Unique Games on a constraint graph G=(V,E)G=(V,E), let AA be the adjacency matrix of GG and let MM denote the nk×nknk\times nk matrix such that the k×kk\times k block MuvM_{uv} is equal to wuvΠuvw_{uv}\Pi_{uv}. We use Πuv\Pi_{uv} to denote the k×kk\times k matrix of the permutation πuv\pi_{uv}. Namely,

Note that πuv=(πvu)1\pi_{uv}=(\pi_{vu})^{-1} implies Πuv=ΠvuT\Pi_{uv}={\Pi_{vu}}^{T} therefore MM is symmetric. We can view MM as the adjacency matrix of the Label-Extended graph for that instance of Unique Games. Let kk denote the size of the alphabet. We will denote this instance by U=(G,M,k)\mathcal{U}=(G,M,k).

(Characteristic vector of a labeling) For any labeling of the vertices L={Lu}uVL=\{L_{u}\}_{u\in V} with Lu[k]L_{u}\in[k] we define the characteristic vector of that labeling as the nknk dimensional vector yLy^{L} with

We will often normalize such vectors to make them unit vectors. In this case, we have

With some abuse of notation we will refer to both vectors as characteristic vectors of labelings whenever it is clear from the context.

The following is an important observation:

Assume that a certain instance U\mathcal{U} is completely satisfiable. Namely, assume that there is a labeling L={Lu}uVL=\{L_{u}\}_{u\in V} with Lu[k]L_{u}\in[k] such that, for all (u,v)E(u,v)\in E we have πuv(Lu)=Lv\pi_{uv}(L_{u})=L_{v}. Then the characteristic vector yLy^{L} is a eigenvector of MM with eigenvalue dd, if GG is a dd-regular graph. It is also an eigenvector of the laplacian LML_{M} with eigenvalue .

3 ΓΓ\Gamma-Max-Lin Instances

Let Γ\Gamma be an abelian group. A Γ\Gamma-Max-Lin instance of Unique Games on graph G(V,E)G(V,E) has for each edge (u,v)E(u,v)\in E a constraint of the form xuxv=cuvx_{u}-x_{v}=c_{uv}, where xu,xvx_{u},x_{v} are variables taking values in Γ\Gamma and cuvΓc_{uv}\in\Gamma. The alphabet kk is the size of the group Γ\Gamma. We note that if a pair of labels (Lu,Lv)Γ×Γ(L_{u},L_{v})\in\Gamma\times\Gamma satisfies a constraint xuxv=cuvx_{u}-x_{v}=c_{uv} of the above form, then for all iΓi\in\Gamma, the pair (Lu+i,Lv+i)Γ×Γ(L_{u}+i,L_{v}+i)\in\Gamma\times\Gamma also satisfies the constraint. This implies that any labeling {Lu}uV\{L_{u}\}_{u\in V} of the vertices of GG is shift-invariant. Namely, for all iΓi\in\Gamma, the labelings {Lu}uV\{L_{u}\}_{u\in V} and {Lu+i}uV\{L_{u}+i\}_{u\in V} satisfy the same set of constraints.

4 The Spectrum of Cayley Graphs

For the purposes of this paper, we will use a generalized definition for Cayley graphs. For background on the standard definition and properties of Cayley graphs see, for example, [Kas02, Spi09].

We note that the standard definition of a Cayley graph, given a set SS of generators, is a special case of the above, where f(g1g2)=1f(g_{1}-g_{2})=1 if g1g2Sg_{1}-g_{2}\in S and otherwise.

Let C(Ω)C(\Omega) be Cayley graph of an abelian group Ω\Omega, as in definition 4, and AC(Ω)A_{C(\Omega)} be its adjacency matrix. Then for every gΩg\in\Omega, AC(Ω)χg=Ωf^(g)χgA_{C(\Omega)}\chi_{g}=|\Omega|\hat{f}(g)\chi_{g}. Namely, every character χg\chi_{g} is an eigenvector of AC(Ω)A_{C(\Omega)} and Ωf^(g)|\Omega|\hat{f}(g) is its corresponding eigenvalue. Here, with f^\hat{f} we denote fourier coefficients.

It is enough to prove the equality for each entry of AC(Ω)χgA_{C(\Omega)}\chi_{g}, that is for every vertex xΩx\in\Omega we will show {AC(Ω)χg}x=Ωf^(g)χg(x)\{A_{C(\Omega)}\chi_{g}\}_{x}=|\Omega|\hat{f}(g)\chi_{g}(x). We calculate:

5 The Khot-Vishnoi Graph

In [KV05], the authors considered the following family of graphs:

For parameters nn and ϵ\epsilon, let N=2nN=2^{n} and n=2kn=2^{k}. Denote by F\mathcal{F} the family of all boolean functions on {1,1}k\{-1,1\}^{k}. For f,gFf,g\in\mathcal{F} define the product fgfg as (fg)(x)=f(x)g(x)(fg)(x)=f(x)g(x). Let H={χSS[k]}H=\{\chi_{S}|S\subseteq[k]\} be the set of characters of the group F2n\mathbf{F}_{2}^{n}. Consider the equivalence relation \equiv on F\mathcal{F} defined as fgf\equiv g iff there is an S[k]S\subseteq[k] such that f=gχSf=g\chi_{S}. This relation partitions F\mathcal{F} into equivalence classes P1,,Pm\mathcal{P}_{1},\cdots,\mathcal{P}_{m}, where m=Nn=2nnm=\frac{N}{n}=\frac{2^{n}}{n}. For each equivalence class Pi\mathcal{P}_{i}, we could pick an arbitrary representative piPip_{i}\in\mathcal{P}_{i}, so that

For 0<ϵ10<\epsilon\leq 1 let fϵFf\in_{\epsilon}\mathcal{F} denote a random boolean function on {1,1}k\{-1,1\}^{k} where for every x{1,1}kx\in\{-1,1\}^{k}, independently, f(x)=1f(x)=1 with probability 1ϵ1-\epsilon and f(x)=1f(x)=-1 with probability ϵ\epsilon. For the given parameter ϵ\epsilon and boolean functions f,gFf,g\in\mathcal{F} let

The KVn,ϵ\mathcal{K}\mathcal{V}_{n,\epsilon} constraint graph with parameters nn and ϵ\epsilon can now be defined to have vertex set V={P1,,Pm}V=\{\mathcal{P}_{1},\cdots,\mathcal{P}_{m}\}. For every fPif\in\mathcal{P}_{i} and gPjg\in\mathcal{P}_{j}, there is an edge between the vertices Pi\mathcal{P}_{i} and Pj\mathcal{P}_{j} with weight wtϵ({f,g})\mathbf{wt}_{\epsilon}(\{f,g\})

We next describe the set of constraints. The set of labels for the instance will be 2[k]:={S:S[k]}2^{[k]}:=\{S:S\subseteq[k]\}, i.e. the set of labels [n][n] is identified with the set 2[k]2^{[k]} (and, thus, n=2kn=2^{k}). This identification will be used from now on. Let fPif\in\mathcal{P}_{i} with f=piχSf=p_{i}\chi_{S} and gPjg\in\mathcal{P}_{j} with g=piχTg=p_{i}\chi_{T} for some S,T[k]S,T\subseteq[k]. The constraint πe\pi_{e} for the edge e{f,g}e\{f,g\} corresponding to the pair of functions f,gf,g can now be defined as

Here, \vartriangle is the symmetric difference operator on sets.

We denote the label-extended graph of the instance described above with KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon}. The graph KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} has vertex set {1,1}n\{-1,1\}^{n}. Between two vertices f,g{1,1}nf,g\in\{-1,1\}^{n} there is an edge with weight nwtϵ({f,g})n\cdot\mathbf{wt}_{\epsilon}(\{f,g\}).

It will be useful, for the purposes of this paper, to consider an equivalent definition of KVn,ϵ\mathcal{K}\mathcal{V}_{n,\epsilon} and KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon}, by translating to the {0,1}\{0,1\} language. Let Hn={0,1}nH_{n}=\{0,1\}^{n} be the group F2n\mathbf{F}_{2}^{n} with addition (modulo 2) as group operation. Let n=2kn=2^{k}. Then the set HH is the Hadamard code on nn bits. Namely,

Note that H=n|H|=n and HH is a subgroup of HnH_{n}. The graph KVn,ϵ\mathcal{K}\mathcal{V}_{n,\epsilon} is just a Cayley graph of the quotient group of HnH_{n} by HH, namely of the group Q=HnHQ=H_{n}\diagup H. The vertex set V={P1,,Pm}V=\{\mathcal{P}_{1},\cdots,\mathcal{P}_{m}\} consists of the Nn\frac{N}{n} cosets of HH. For every ff in the coset Pi\mathcal{P}_{i} and for every gg in the coset Pj\mathcal{P}_{j}, there is an edge between the vertices Pi\mathcal{P}_{i} and Pj\mathcal{P}_{j} with weight wtϵ({f,g})\mathbf{wt}_{\epsilon}(\{f,g\}). For simplicity of notation, we will identify the cosets with their representatives, i.e. if x=Pix=\mathcal{P}_{i} is a vertex of KVn,ϵ\mathcal{K}\mathcal{V}_{n,\epsilon}, then we will use xx to refer both to the coset Pi\mathcal{P}_{i} as well as to its representative pip_{i}.

The total weight of all the edges between two cosets x=Pix=\mathcal{P}_{i} and y=Pjy=\mathcal{P}_{j} is

Note that the graph is nn-regular, namely the total weight of edges adjacent to any node xx is

The label-extended graph KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} is a Cayley graph of HnH_{n}, since for every two vertices u,v{0,1}nu,v\in\{0,1\}^{n} we have that the corresponding entry of the adjacency matrix of KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} depends only on the hamming distance of uu and vv, or, in other words, on uv=u+vu-v=u+v (modulo 2). Moreover, KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} can be described as an “ϵ\epsilon-perturbed” version of the standard hypercube graph: the weight between two vertices u,v{0,1}nu,v\in\{0,1\}^{n} is proportional to the probability that if one starts from the nn-bit string uu and flips each one of its bits independently with probability ϵ\epsilon, then one gets the nn-bit string vv. The following appears in [KV05]. We give the informal statement here, for simplicity. We refer the reader to [KV05] for the full statement and details.

(Integrality Gap for Unique Games, Informal Statement) Let nn be an integer and ϵ>0\epsilon>0 be a parameter. Then for the instance of Unique Games described above, with label-extended graph KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} the following hold:

Every labeling L:V[n]L:V\rightarrow[n] satisfies at most a 1nϵ\frac{1}{n^{\epsilon}} fraction of the total weight of the edges.

The standard SDP relaxation for Unique Games has objective value greater than 19ϵ1-9\epsilon.

6 Additional Notation

For simplicity of the presentation, we will assume that GG is a regular graph, namely vwuv=d\sum_{v}w_{uv}=d for all nodes uu. At the end of the next section, we show that the results can easily be generalized to non-regular graphs.

For the sake of convenience, in the rest of the paper we will often use the term eigenspace to refer to the space spanned by a set of eigenvectors that don’t necessarily have the same eigenvalue. The relevant set of eigenvectors that span this space will always be clear from the context. We use poly(n)\text{poly}(n) to refer to some polynomial function in nn.

We also use the following notation for the time needed to compute eigenvalues and eigenvectors:

Let U=(G,M,k)\mathcal{U}=(G,M,k) be an instance of Unique Games on a graph GG on nn nodes and with alphabet size kk. Let MM be the adjacency matrix of the label-extended graph of GG. Let WW be some space spanned by eigenvectors of MM. We denote by TW(M)T_{W}(M) the time needed to compute an eigenbasis and the corresponding eigenvalues of WW. Note that TW(M)T_{W}(M) is polynomial in the dimension of MM, namely TW(M)=poly(nk)T_{W}(M)=\text{poly}(n\cdot k).

We also note that in the rest of this paper we use the terms “assignment” and “labeling” interchangeably.

Recovering Solutions by Spectral Methods

In this section, we will show how, given a (1ϵ)(1-\epsilon) satisfiable instance of Unique Games U=(G,M,k)\mathcal{U}=(G,M,k) on graph GG and with alphabet size kk, the eigenvectors of the label-extended graph MM may be used to recover good assignments. Specifically, we show the following:

(Main) Let U=(G,M,k)\mathcal{U}=(G,M,k) be a (1ϵ)(1-\epsilon) satisfiable instance of Unique Games on a dd-regular graph GG on nn nodes with alphabet size kk. Let MM be the adjacency matrix of the label-extended graph of GG as above. Let WW be the space spanned by eigenvectors of MM that have eigenvalue greater than (1γ)d(1-\gamma)d, for γ>8ϵ\gamma>8\epsilon. There is an algorithm that runs in time 2O(γϵdim(W))+poly(nk)2^{O(\frac{\gamma}{\epsilon}\text{dim}(W))}+\text{poly}(n\cdot k) and finds an assignment that satisfies at least (1O(ϵγ8ϵ+ϵ))(1-O(\frac{\epsilon}{\gamma-8\epsilon}+\epsilon)) fraction of the constraints.

In particular, the theorem implies that for small enough ϵ\epsilon, for every 1ϵ1-\epsilon satisfiable instance of Unique Games that satisfies the assumptions of the theorem, one can find an assignment that satisfies more than 90 percent of the constraints in time that depends on the spectral profile of MM. We also remark the following:

Remark: If the algorithm of theorem 8 fails to find an assignment that satisfies at least (1O(ϵγ8ϵ+ϵ))(1-O(\frac{\epsilon}{\gamma-8\epsilon}+\epsilon)) fraction of the constraints, then for the input U\mathcal{U}, every assignment satisfies less than (1ϵ)(1-\epsilon) fraction of the constraints. Here, the constant in the O()O(\cdot) notation is the same as the one guaranteed by the theorem above.

We also consider the special case where the constraints are arbitrary Γ\Gamma-Max-Lin. For such constraints, we prove theorem 9 below.

Let U=(G,M,k)\mathcal{U}=(G,M,k) be a (1ϵ)(1-\epsilon) satisfiable Γ\Gamma-Max-Lin instance of Unique Games on a dd-regular graph GG on nn nodes with alphabet size kk. Let S(1γ)S_{(1-\gamma)} be the space spanned by eigenvectors of GG that have eigenvalue greater than (1γ)d(1-\gamma)d. Assume moreover that every unit-length vector ϕS(1γ)\phi\in S_{(1-\gamma)} has ϕCn\left\lVert\phi\right\rVert_{\infty}\leq\frac{C}{\sqrt{n}} for some constant CC. Then, for γ=Ω(ϵ)\gamma=\Omega(\sqrt{\epsilon}), there is an algorithm that runs in time 2O(kDS)+poly(nk)2^{O(k\cdot D_{S})}+\text{poly}(n\cdot k) and finds an assignment that satisfies at least (1ζ)(1-\zeta) fraction of the constraints for some ζ=O(ϵ)\zeta=O(\sqrt{\epsilon}). Here DSD_{S} denotes the dimension of S(1γ)S_{(1-\gamma)} and the constant in the Ω()\Omega(\cdot) notation depends on CC.

Assuming that TS(1γ)(M)<2O(kDS)T_{S_{(1-\gamma)}}(M)<2^{O(k\cdot D_{S})} the algorithm in theorem 9 has running time that solely depends on the spectral profile of GG, since roughly, the space of eigenvectors of MM with large eigenvalue has dimension equal to kk times the dimension of the corresponding eigenspace of GG.

The result for expander graphs as it appears in [KT07] can be derived as a corollary of theorem 9.

(Unique Games are Easy on Expanders) Let U=(G,M,k)\mathcal{U}=(G,M,k) be a (1ϵ)(1-\epsilon) satisfiable Γ\Gamma-Max-Lin instance of Unique Games. Assume moreover that GG is a dd-regular spectral expander on nn nodes. Namely, the second eigenvalue of the adjacency matrix of GG is λ(1γ)d\lambda\leq(1-\gamma)d, for γ=Ω(ϵ)\gamma=\Omega(\sqrt{\epsilon}). Then, there is a polynomial time algorithm that finds an assignment that satisfies at least (1ζ)(1-\zeta) fraction of the constraints for some ζ=O(ϵ)\zeta=O(\sqrt{\epsilon}).

We remark that the Γ\Gamma-Max-Lin requirement is necessary for corollary 10. The proof fails to produce similar guarantee for the general case. This is due to the fact that the spectral properties of label-extended graphs that correspond to arbitrary Unique Games are poorly understood.

Proof Overview. Assume we are given a (1ϵ)(1-\epsilon) satisfiable instance U=(G,M,k)\mathcal{U}=(G,M,k) that satisfies the assumptions of theorem 8. We define a completely satisfiable game U~=(G,M~,k)\mathcal{\widetilde{U}}=(G,\widetilde{M},k) (with value 1) that “corresponds” to U\mathcal{U} as follows:

For the sake of the proof, we will assume that an almost satisfiable instance U=(G,M,k)\mathcal{U}=(G,M,k) is constructed as follows:

Let U~=(G,M~,k)\mathcal{\widetilde{U}}=(G,\widetilde{M},k) be a completion of U\mathcal{U}.

Let an adversary pick the ϵ\epsilon fraction of edges that were unsatisfied in U\mathcal{U} and change their constraints back to the original ones. We can now think of MM as a “perturbation” of M~\widetilde{M} and U\mathcal{U} as a “perturbed” game of U~\mathcal{\widetilde{U}}.

Let WW be the span of the eigenvectors of MM with eigenvalue at least (1γ)d(1-\gamma)d, for some γ>8ϵ\gamma>8\epsilon. The algorithm simply looks at a set SW\mathcal{S}\subseteq W of appropriately many candidate vectors and “reads-off” an assignment. We describe later on how to choose the set S\mathcal{S}.

Recover-SolutionS(U)\texttt{Recover-Solution}_{\mathcal{S}}(\mathcal{U})

For each xSx\in\mathcal{S}, construct an assignment LxL_{x} by assigning to each vertex uu, the index corresponding to the largest entry in the block (xu1,,xuk)(x_{u1},\ldots,x_{uk}) i.e. Lx(u)=argmaxixuiL_{x}(u)=\arg\max_{i}{x_{ui}}.

Out of all assignments LxL_{x} for xSx\in\mathcal{S}, choose the one satisfying the maximum number of constraints.

To choose S\mathcal{S}, we look at the highest eigenvectors for the matrix M~\widetilde{M}. Those are the assignment eigenvectors, namely the characteristic vectors of the (perfectly) satisfying assignments of U~\mathcal{\widetilde{U}} (since they all have eigenvalue equal to dd). We will first observe that every such eigenvector yy is close to some vector in WW, and the length of the projection of yy onto WW depends on ϵ,γ\epsilon,\gamma. We then identify a set of nice vectors NW\mathcal{N}\subseteq W such that the above algorithm works for any vector vv close to some vector in N\mathcal{N}. These are going to be the vectors that are close to the assignment eigenvectors. We then construct a set SW\mathcal{S}\subseteq W of test vectors such that at least one vector vSv\in\mathcal{S} is close to a vector in N\mathcal{N}. S\mathcal{S} is simply an epsilon-net for WW. Lastly, we go over all vectors in S\mathcal{S} until we find vv. The running time of the algorithm will depend on the size of S\mathcal{S} which, for an appropriately defined epsilon-net, it is exponential in the dimension of WW.

Our first step is to identify an appropriate eigenspace WW of MM such that characteristic vectors of good labelings have large projection onto WW.

For the matrix M~\widetilde{M}, let L={L(u)}uV(G)\mathcal{L}=\{L(u)\}_{u\in V(G)} be a completely satisfying assignment. Define vector y(L)y^{(\mathcal{L})} as

It is easy to see that y(L)y^{(\mathcal{L})} is an eigenvector of M~\widetilde{M} with eigenvalue dd. Since y(L)y^{(\mathcal{L})} corresponds to the satisfying assignment L\mathcal{L}, it can be seen as the characteristic vector of the assignment. We refer to such vectors as the “assignment” eigenvectors. Our next goal is to show that, for some appropriate choice of γ\gamma, the eigenspace WW contains vectors close to an assignment vector.

(Closeness) For every completely satisfying assignment L\mathcal{L} there is a unit vector vLWv_{\mathcal{L}}\in W such that vL=αy(L)+βy(L)v_{\mathcal{L}}=\alpha y^{(\mathcal{L})}+\beta{y^{(\mathcal{L})}}_{\perp}, with α>0\alpha>0 and β2ϵγ|\beta|\leq\sqrt{\frac{2\epsilon}{\gamma}}. Here both y(L){y^{(\mathcal{L})}} and y(L){y^{(\mathcal{L})}}_{\perp} are unit vectors and y(L)y(L){y^{(\mathcal{L})}}_{\perp}\perp y^{(\mathcal{L})}. By taking, for example, γ200ϵ\gamma\geq 200\epsilon, we have β110|\beta|\leq\frac{1}{10}.

We can easily see that (y(L))TMy(L)d(12ϵ)(y^{(\mathcal{L})})^{T}My^{(\mathcal{L})}\geq d(1-2\epsilon). We can now write y(L)=avL+b(vL)y^{(\mathcal{L})}=av_{\mathcal{L}}+b(v_{\mathcal{L}})_{\perp}, with a>0a>0, vLW{v_{\mathcal{L}}}\in W and (vL)W({v_{\mathcal{L}}})_{\perp}\in W^{\perp}. We calculate:

from which we get that b2ϵγ|b|\leq\sqrt{\frac{2\epsilon}{\gamma}}.

Now, we can in turn express vL=αy(L)+βy(L)v_{\mathcal{L}}=\alpha y^{(\mathcal{L})}+\beta{y^{(\mathcal{L})}}_{\perp}, where α=vL,y(L)=a=1b212ϵγ\alpha=\langle v_{\mathcal{L}},y^{(\mathcal{L})}\rangle=a=\sqrt{1-b^{2}}\geq\sqrt{1-\frac{2\epsilon}{\gamma}} or, equivalently, β=1α22ϵγ|\beta|=\sqrt{1-\alpha^{2}}\leq\sqrt{\frac{2\epsilon}{\gamma}}. ∎

We have now managed to show that there exists a set of “nice” vectors N={vL}W\mathcal{N}=\{v_{\mathcal{L}}\}\subseteq W. The next claim shows that, if we knew the vLv_{\mathcal{L}}’s then we could set S=N\mathcal{S}=\mathcal{N} and the algorithm Recover-SolutionN(U)\texttt{Recover-Solution}_{\mathcal{N}}(\mathcal{U}) would return a highly satisfying assignment for U\mathcal{U}.

If xx is a vector such that x=αy(L)+βy(L)x=\alpha y^{(\mathcal{L})}+\beta{y^{(\mathcal{L})}}_{\perp} for some y(L)y^{(\mathcal{L})} with α>0\alpha>0 and y(L)y(L){y^{(\mathcal{L})}}_{\perp}\perp y^{(\mathcal{L})}, then the coordinate xuL(u){x_{uL(u)}} is maximum in absolute value in at least (12β2α2)n(1-\frac{2\beta^{2}}{\alpha^{2}})n blocks.

Within each block uu, in order for coordinate L(u)L(u) to be no longer the maximum one, it must happen that for some jj

Since y(L)=1\left\lVert{y^{(\mathcal{L})}}_{\perp}\right\rVert=1, this can only happen for at most 2β2α2n\frac{2\beta^{2}}{\alpha^{2}}n blocks. ∎

1.2 Eigenspace Enumeration: Finding the Set 𝒮𝒮\mathcal{S} of Test Vectors

Our next step is to show that the search for a good assignment can be done in time exponential to the number of large eigenvalues of MM, namely the dimension of WW.

Since we don’t know the y(L)y^{(\mathcal{L})} (if we did, we would be done), and the space WW contains infinitely many unit vectors, we cannot identify N\mathcal{N}. To get around this, we discretize WW with an appropriate epsilon-net. If we let w(0),,w(dim(W)1)w^{(0)},\ldots,w^{(\text{dim}(W)-1)} be an eigenbasis for WW. We define the set S\mathcal{S} as

It can be calculated (see, for instance [FO05]) that the number of points in the set S\mathcal{S} is at most 2O(γϵdim(W))2^{O(\frac{\gamma}{\epsilon}\text{dim}(W))}.

To conclude the proof of theorem 8, it remains to show that S\mathcal{S} has a vector close to a nice vector vLNv_{\mathcal{L}}\in\mathcal{N}. By construction, S\mathcal{S} contains at least one vector close to every vector in N\mathcal{N} and thus it also contains at least one vector vv such that v=αvL+β(vL)v=\alpha v_{\mathcal{L}}+\beta({v_{\mathcal{L}}})_{\perp} for some L\mathcal{L} and β2ϵγ|\beta|\leq\sqrt{\frac{2\epsilon}{\gamma}}. Together with claim 12, this implies that we can also write v=ay(L)+by(L)v=ay^{(\mathcal{L})}+by^{(\mathcal{L})}_{\perp}, with b2ϵγ+2ϵγ=22ϵγ|b|\leq\sqrt{\frac{2\epsilon}{\gamma}}+\sqrt{\frac{2\epsilon}{\gamma}}=2\sqrt{\frac{2\epsilon}{\gamma}}. Thus, by claim 13, for this vector vv, Recover-SolutionS(U)\texttt{Recover-Solution}_{\mathcal{S}}(\mathcal{U}) recovers an assignment which agrees with y(L)y^{(\mathcal{L})} in (1O(ϵγ8ϵ))(1-O(\frac{\epsilon}{\gamma-8\epsilon})) fraction of the blocks. We will consider that a constrain is violated if one of the two following events happen: either the assignment we recovered does not agree with the initial perfectly satisfying assignment or the constraint is one of the (at most) ϵ\epsilon fraction of constraints that were changed. It follows that the assignment which we recovered, violates constraints on edges that have total weight at most O(ϵγ8ϵ)nd+ϵndO(\frac{\epsilon}{\gamma-8\epsilon})nd+\epsilon nd. Since the total weight of constraints is nd/2nd/2, theorem 8 follows. The algorithm runs in 2O(γϵdim(W))+TW(M)=2O(γϵdim(W))+poly(nk)2^{O(\frac{\gamma}{\epsilon}\text{dim}(W))}+T_{W}(M)=2^{O(\frac{\gamma}{\epsilon}\text{dim}(W))}+\text{poly}(n\cdot k) time.

2 Proof of Theorem 9

Let θ>0\theta>0 such that γθΩ(ϵγ)\gamma\geq\theta\geq\Omega(\epsilon\gamma). Let WW be the span of eigenvectors of MM with eigenvalue at least (1θ)d(1-\theta)d. Let YY be the span of eigenvectors of M~\widetilde{M} with eigenvalue at least (1γ)d(1-\gamma)d. In particular, YY contains all the assignment eigenvectors, namely the characteristic vectors of the (perfectly) satisfying assignments of U~\mathcal{\widetilde{U}} (since they all have eigenvalue equal to dd). Let YY_{\perp} be the orthogonal complement of YY. For the proof, we will use theorem 8 together with a bound on the dimension of W. Roughly we show that the dimension of WW is at most as large as the dimension of YY which, in turn, is at most kk times the dimension of S(1γ)S_{(1-\gamma)}. To conclude, we apply theorem 8 for the eigenspace WW.

We will assume w.l.o.g that the graphs we are dealing with are connected. Otherwise, we apply the results to each connected component individually.

Assume γθΩ(ϵγ)\gamma\geq\theta\geq\Omega(\epsilon\gamma) for some appropriately large constant. Then any unit-length vector wWw\in W can be expressed as αy+βy\alpha y+\beta y_{\perp} where yYy\in Y and yYy_{\perp}\in Y_{\perp} are both unit length vectors, and βO(θγ3)|\beta|\leq O(\sqrt{\frac{\theta}{\gamma^{3}}}). By taking, for example, γΩ(θ)\gamma\geq\Omega(\sqrt{\theta}) and, consequently, γΩ(ϵ)\gamma\geq\Omega(\sqrt{\epsilon}), we get β18|\beta|\leq\frac{1}{8}.

Combining claims 15 and 16, we can proceed as follows: From claim 16 we obtain that WW has dimension dim(W)dim(Y)\text{dim}(W)\leq\text{dim}(Y). Otherwise, we would find a vector orthogonal to all the vectors in YY which cannot be close to their span. From claim 15 we obtain dim(Y)=kDS\text{dim}(Y)=k\cdot D_{S}. To conclude the proof of theorem 9, we apply theorem 8 with WW being the (at most) kDSk\cdot D_{S} dimensional eigenspace of MM with eigenvalues (1θ)d\geq(1-\theta)d.

Finally, it remains to argue about the running time of the algorithm. Since dim(W)kDS\text{dim}(W)\leq k\cdot D_{S}, the number of points is exponential in kDSkD_{S}. Hence, the algorithm runs in 2O(kDS)+TW(M)=2O(kDS)+poly(nk)2^{O(k\cdot D_{S})}+T_{W}(M)=2^{O(k\cdot D_{S})}+\text{poly}(n\cdot k) time.

We next prove claim 16. Towards this end, we will apply some results from matrix perturbation theory. To find appropriate γ\gamma and θ\theta such that the eigenspaces WW and YY are close, we use the following claim which essentially appears in [DK70] as the sin θ\theta theorem and was used in [KT07].

Let ww be a unit length eigenvector of MM with eigenvalue λ(1θ)d\lambda\geq(1-\theta)d, and let λs\lambda_{s} denote the largest eigenvalue of M~\widetilde{M} which is smaller than (1γ)d(1-\gamma)d. Then, ww can be written as αy+βy\alpha y+\beta y_{\perp} with β(M~M)w(λλs)|\beta|\leq\frac{\left\lVert(\widetilde{M}-M)w\right\rVert}{(\lambda-\lambda_{s})}. Here yYy\in Y and yYy_{\perp}\in Y^{\perp} are unit length vectors.

As seen by the previous lemma, in order to prove that the space YY does not change by much due to the perturbation, we simply need to bound (M~M)w\left\lVert(\widetilde{M}-M)w\right\rVert. We will need the fact that ww is somewhat “uniform” across blocks. To formalize this, let wˉ\bar{w} be the nn-dimensional vector such that wˉu=wu\bar{w}_{u}=\left\lVert w_{u}\right\rVert where wuw_{u} is the kk-dimensional vector (wu1,,wuk)T(w_{u1},\ldots,w_{uk})^{T}. We then show that wˉ\bar{w} is very close to a vector in S(1γ)S_{(1-\gamma)}.

If ww is a unit-length eigenvector of MM with eigenvalue more than (1θ)d(1-\theta)d and wˉ\bar{w} as above, then we can write wˉ\bar{w} as wˉ=aϕ+bϕ\bar{w}=a\phi+b\phi_{\perp} with bθγ|b|\leq\sqrt{\frac{\theta}{\gamma}} and ϕS(1γ)\phi\in S_{(1-\gamma)}, ϕS(1γ)\phi_{\perp}\in S_{(1-\gamma)}^{\perp}, both having unit length.

Since, ww corresponds to a large eigenvalue, we have that

Here the second inequality follows from the fact that the operator norm of Πuv\Pi_{uv} is at most 1. Writing wˉ\bar{w} as aϕ+bϕa{\phi}+b{\phi}_{\perp}, we get

Using the above, and the fact that the matrix M~\widetilde{M} is only perturbed in ϵ\epsilon fraction of the edges, we can now bound (M~M)w\left\lVert(\widetilde{M}-M)w\right\rVert as follows.

(M~M)wO(θγ)d\left\lVert(\widetilde{M}-M)w\right\rVert\leq O\left(\sqrt{\frac{\theta}{\gamma}}\right)d

Define the n×nn\times n matrix RR as Ruv=wuv1R_{uv}=w_{uv}\leq 1 when the block (M~M)uv(\widetilde{M}-M)_{uv} has any non-zero entry, and Ruv=0R_{uv}=0 otherwise. Note that if (M~M)uv(\widetilde{M}-M)_{uv} is non-zero, then it must be the (scaled) difference of two permutation matrices. Thus, for all vv we have (M~M)uvwv2Ruvwv\left\lVert(\widetilde{M}-M)_{uv}w_{v}\right\rVert\leq 2R_{uv}\left\lVert w_{v}\right\rVert. We calculate:

To estimate Rwˉ\left\lVert R\bar{w}\right\rVert, we break it up as

Since each row of RR has total sum of entries at most dd, bRϕθγd|b|\left\lVert R\phi_{\perp}\right\rVert\leq\sqrt{\frac{\theta}{\gamma}}d. Also,

Since RR has a total sum of entries at most ϵnd\epsilon nd, this expression is maximized when it has dd 1s in ϵn\epsilon n rows. This gives RϕCϵd\left\lVert R\phi\right\rVert\leq C\sqrt{\epsilon}d. And, putting everything together we obtain

We note that the O()O(\cdot) in the above expression depends on CC. We obtained the last inequality by using the assumption of the lemma θΩ(ϵγ)\theta\geq\Omega(\epsilon\gamma). ∎

Combining the above bound with claim 17, we get that any unit-length vector wWw\in W can be expressed as αy+βy\alpha y+\beta y_{\perp} where yYy\in Y and βO(θγd1(1θ)dλs)|\beta|\leq O(\sqrt{\frac{\theta}{\gamma}}d\cdot\frac{1}{(1-\theta)d-\lambda_{s}}). Recall that λs\lambda_{s} was smaller than (1γ)d(1-\gamma)d, which implies

To conclude the proof of claim 16 we take θ<<γ\theta<<\gamma. We calculate:

It is sufficient to consider γΩ(θ)\gamma\geq\Omega(\sqrt{\theta}) and observe that this also implies, by the assumption θΩ(ϵγ)\theta\geq\Omega(\epsilon\gamma) that γΩ(ϵ)\gamma\geq\Omega(\sqrt{\epsilon}).

3 Generalizing to Non-Regular Graphs

It remains to show that the above results hold for non-regular graphs. We will consider eigenvectors and eigenvalues of the Laplacian matrix of the label-extended graph. Let GG and MM as before. Let DD be the nk×nknk\times nk diagonal matrix with block Duu=deg(u)IkD_{uu}=deg(u)\cdot I_{k}, where IkI_{k} is the k×kk\times k identity matrix. The Laplacian of the label-extended graph is the matrix LM=DML_{M}=D-M. Let dd be the average degree of GG, and therefore of MM, namely d=uVdund=\frac{\sum_{u\in V}d_{u}}{n}. We re-state the main theorem 8 to capture the non-regular case.

Let U=(G,M,k)\mathcal{U}=(G,M,k) be a (1ϵ)(1-\epsilon) satisfiable instance of Unique Games and WW the eigenspace of LML_{M} with eigenvalues less than γd\gamma d, for γ8ϵ\gamma\geq 8\epsilon. There is an algorithm that runs in time 2O(γϵdim(W))+poly(nk)2^{O(\frac{\gamma}{\epsilon}\text{dim}(W))}+\text{poly}(n\cdot k) and finds an assignment that satisfies at least (1O(ϵγ8ϵ+ϵ))(1-O(\frac{\epsilon}{\gamma-8\epsilon}+\epsilon)) fraction of the constraints.

Let WW be the span of the eigenvectors of LML_{M} with eigenvalue less than γd\gamma d, for some γ8ϵ\gamma\geq 8\epsilon. We again recover highly satisfying assignments using our algorithm Recover-SolutionS(U)\texttt{Recover-Solution}_{\mathcal{S}}(\mathcal{U}). Similarly to the previous section, the set SW\mathcal{S}\subseteq W will be a discretization of our new WW.

All we need to show is the following analog of claim 12. We use the same notation as in section 3.1.

For every completely satisfying assignment L\mathcal{L} there is a unit vector vLWv_{\mathcal{L}}\in W such that vL=αy(L)+βy(L)v_{\mathcal{L}}=\alpha y^{(\mathcal{L})}+\beta{y^{(\mathcal{L})}}_{\perp}, with β2ϵγ|\beta|\leq\sqrt{\frac{2\epsilon}{\gamma}}. Here y(L)y^{(\mathcal{L})} and y(L){y^{(\mathcal{L})}}_{\perp} are both unit length vectors and y(L)y(L){y^{(\mathcal{L})}}_{\perp}\perp y^{(\mathcal{L})}. By taking, for example, γ200ϵ\gamma\geq 200\epsilon, we have β110|\beta|\leq\frac{1}{10}.

We first observe that the assignment eigenvectors yLy^{\mathcal{L}} are eigenvectors of LM~L_{\widetilde{M}} with eigenvalue .

We can now write y(L)=avL+b(vL)y^{(\mathcal{L})}=av_{\mathcal{L}}+b(v_{\mathcal{L}})_{\perp}, with vLW{v_{\mathcal{L}}}\in W and (vL)W({v_{\mathcal{L}}})_{\perp}\in W^{\perp}. We calculate:

from which we get that b2ϵγ|b|\leq\sqrt{\frac{2\epsilon}{\gamma}}.

Now, we can in turn express vL=αy(L)+βy(L)v_{\mathcal{L}}=\alpha y^{(\mathcal{L})}+\beta{y^{(\mathcal{L})}}_{\perp}, where α=vL,y(L)=a=1b212ϵγ\alpha=\langle v_{\mathcal{L}},y^{(\mathcal{L})}\rangle=|a|=\sqrt{1-b^{2}}\geq\sqrt{1-\frac{2\epsilon}{\gamma}} or, equivalently, β=1α22ϵγ|\beta|=\sqrt{1-\alpha^{2}}\leq\sqrt{\frac{2\epsilon}{\gamma}}. ∎

To conclude the proof of the theorem, we observe that the analog of claim 13 follows immediately. We define S\mathcal{S} as in section 3.1.2, namely

where w(i)w^{(i)} are now eigenvectors of our new WW. The theorem follows.

Solving Unique Games on the Khot-Vishnoi Instance

In this section, we show that, when run on the Khot-Vishnoi integrality gap instance, our main algorithm runs in quasi-polynomial time and correctly decides the (un)satisfiability of the instance. Namely, we show the following:

Our main algorithm (as in theorem 8) when given input the integrality gap instance of Khot and Vishnoi as described in the preliminaries section, with label-extended graph KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon}, runs in time npoly(logn)n^{\text{poly}(\log n)} and correctly decides that the instance is highly unsatisfiable. Here N=2nN=2^{n} the number of nodes of HnH_{n} and nn is the alphabet size.

The proof goes by showing that the eigenspace WW as per theorem 8 has relatively low dimension. The theorem guarantees that exhaustive search in WW would find a good assignment if such assignment existed and consequently, failure to find such an assignment, will correctly classify the Khot-Vishnoi instance as unsatisfiable.

The rest of this section is devoted to the proof of claim 22. We use the equivalent definition of KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} as a Cayley graph of HnH_{n} which is an “ϵ\epsilon-perturbed” version of the hypercube graph. We will need the following claim for the spectrum of KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon}. We denote the adjacency matrix of KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} by MKVM_{\mathcal{K}\mathcal{V}}.

The following is true for the spectrum of KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon}.

The eigenvectors are the characters of the group HnH_{n}.

The eigenvalue that corresponds to χω\chi_{\omega} is (12ϵ)rn(1-2\epsilon)^{r}n, where r=ωr=|\omega| is the hamming weight of ω\omega, and appears with multiplicity Cr=(nr)C_{r}={n\choose r}.

The first item above is immediate from the definition of the graph as it appears in the preliminaries section.

For the second item, we just need to calculate the eigenvalues corresponding to χω\chi_{\omega} or, equivalently, the quadratic form

To conclude the proof, we need to argue about the multiplicity of each eigenvalue or, equivalently, the number of characters χω\chi_{\omega} of a given hamming weight rr. It is easily seen that there are Cr=(nr)C_{r}={n\choose r} characters of HnH_{n} that correspond to hamming weight rr. ∎

We will directly apply theorem 8. For this purpose we need to calculate the dimension of the eigenspace WW of KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} for some appropriate γΩ(ϵ)\gamma\geq\Omega(\sqrt{\epsilon}).

Let γ=Cϵ\gamma=C\sqrt{\epsilon} for some sufficiently large constant CC. The dimension of the eigenspace WW of the graph KV~n,ϵ\widetilde{\mathcal{K}\mathcal{V}}_{n,\epsilon} with eigenvalues (1γ)\geq(1-\gamma) is at most DW1ϵnO(1ϵ)=poly(n)=poly(logN)D_{W}\leq\frac{1}{\epsilon}n^{O(\frac{1}{\epsilon})}=\text{poly}(n)=\text{poly}(\log N).

We first need to identify an upperbound on rr such that (12ϵ)r(1γ)(1-2\epsilon)^{r}\leq(1-\gamma), for γ=Cϵ\gamma=C\sqrt{\epsilon}. It is enough to take r=log1γlog12ϵr=\frac{\log 1-\gamma}{\log 1-2\epsilon} . In order to approximate rr, we use the Taylor series expansion:

And for ϵ,γ\epsilon,\gamma small enough we get

We conclude by arguing, according to theorem 8, that our algorithm will run in time bounded by 2DW=21ϵnO(1ϵ)=N1ϵ(logN)O(1ϵ)=N(poly(logN))2^{D_{W}}=2^{\frac{1}{\epsilon}\cdot n^{O(\frac{1}{\epsilon})}}=N^{\frac{1}{\epsilon}\cdot(\log N)^{O(\frac{1}{\epsilon})}}=N^{(\text{poly}(\log N))}. It will fail to find a highly satisfying assignment and thus decide that the above instance is highly unsatisfiable. Here we used the fact that for the graph in question, TSW(M)2DWT_{S_{W}}(M)\leq 2^{D_{W}}.

Acknowledgements

The author would like to thank Madhur Tulsiani for his contribution to the ideas and techniques for developing spectral algorithms for Unique Games. The author would also like to thank the anonymous reviewers for their comments and suggestions.

References