Abstract 1 Introduction 2 Preliminaries 3 Chain Length Counting Bounds via Contractions 4 Leveraging Decomposition to Build Sparsifiers 5 Conclusion References Appendix A Proof of Corollary 21

Multiplicative Error Set System Sparsification:
A Simpler Proof via Chain Length Contraction

Joshua Brakensiek ORCID Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA, USA    Venkatesan Guruswami ORCID Simons Institute for the Theory of Computing and the University of California, Berkeley, CA, USA    Aaron Putterman ORCID School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
Abstract

The chain length of a set family 𝒮2[m] is the largest ascending sequence of sets in containment order in the union-closure of 𝒮. In this work, we provide a significantly simpler and more optimal characterization of the sparsifiability of set systems in terms of their chain length, improving on the work of Brakensiek and Guruswami [STOC 2025]. Our proof relies on a generalization of Karger’s [SODA 1993] famous contraction algorithm and its recent linear algebraic extensions [Khanna-Putterman-Sudan SODA 2024], and our resulting bounds show that, just as VC dimension characterizes the additive sparsifiability of a set system, chain length governs the multiplicative sparsifiability. As a corollary, we obtain improved bounds for weighted CSP sparsification.

Keywords and phrases:
constraint satisfaction problem, chain length, sparsification, VC dimension
Category:
Track A: Algorithms, Complexity and Games
Funding:
Joshua Brakensiek: Supported in part by the Simons Investigator award of Venkatesan Guruswami and NSF grants CCF-2211972 and DMS-2503280.
Venkatesan Guruswami: Supported in part by a Simons Investigator award and NSF award CCF-2211972.
Aaron Putterman: Supported in part by the Simons Investigator Awards of Madhu Sudan and Salil Vadhan and AFOSR award FA9550-25-1-0112.
Copyright and License:
[Uncaptioned image] © Joshua Brakensiek, Venkatesan Guruswami, and Aaron Putterman; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Constraint and logic programming
Acknowledgements:
We thank anonymous reviewers for many helpful comments improving the presentation of this paper. We also thank an anonymous reviewer of [7] whose comments partially inspired Theorem 22. A. P. thanks Sanjeev Khanna and Madhu Sudan for helpful conversations.
Related Version:
Full Version: https://arxiv.org/abs/2605.01508
Editors:
Sayan Bhattacharya, Danupon Nanongkai, Michael Benedikt, and Gabriele Puppis

1 Introduction

In theoretical computer science, sparsification describes a range of methods for reducing the storage size of a mathematical object while preserve some of its essential qualities. In this paper, we study the sparsification of set systems 𝒮2[m] where [m]:={1,2,,m} is a finite set of atoms. Here, we define sparsification analogously to that of Karger’s definition for graph cuts [19], where we find a reweighting w:[m]0 of the atoms such that for every set S𝒮, its weight before and after the reweighting match up to a multiplicative error of 1±ϵ (see Section 2 for a formal definition). In general, the goal is to find a map w such that the support supp(w):={i[m]:w(i)0} is as small as possible.

This framework precisely recover the graph cut sparsification framework of Karger, if we let [m] represent the set of edges of an undirected graph G=(V,E), and let 𝒮 be the family of all possible cuts of G (i.e., sets of the form (A×(VA))E, where AV. In this paper, we prove general bounds on the sparsifiability of all set systems, even when the initial atoms [m] are also weighted.

The general study of set system sparsifiers was initiated by Brakensiek and Guruswami [7] to help resolve questions in the emerging area of CSP sparsification, see Section 1.2 for further background. In particular, they showed that the optimal size of a (weighted) set system sparsifier is closely related to the chain length (coined by Bessiere, Carbonnel, and Katsirelos [5]) of the set system. For a set system 𝒮2[m], we define its chain length CL(𝒮) to be the length of the longest ascending chain of sets S1SCL(𝒮) contained in 𝒮, the closure of S with respect to set union–see Section 2 for a formal definition. We now state the main result of Brakensiek and Guruswami [7] on weighted set system sparsification.

Theorem 1 (Theorem 1.4 of [7]).

For a set system 𝒮2[m], weights w:[m]0 and ϵ>0, there exists a (1±ϵ) sparsifier of 𝒮 which retains only

O(CL(𝒮)log6mϵ2)

many atoms.

Of note, this bound is “near-optimal” in the sense that for any set system 𝒮2[m], there is a set of weights w such that any sparsifier needs at least CL(𝒮) atoms (see Lemma 8.9 of [7]111Lemma 8.9 appears in the full version of their paper at https://arxiv.org/abs/2411.03451.). Furthermore, by a data structure lower bound of Carlson, Kolla, Srivastava and Trevisan [9] for graph cuts, we know that the dependence on ϵ>0 is optimal in the worst case. The main inefficiency in the bound of Theorem 1 is thus the additional factor of log6(m).

The proof of Theorem 1 is rather complex, as the proof first seeks to find near-optimal sparsifiers in the unweighted setting by adapting techniques pioneered by Gilmer [14] in his recent breakthrough on the union-closed set conjecture (see [7] for a much more detailed history of this problem).

1.1 Main Result

Our main result is a much more direct proof of Theorem 1 using a “contraction”-style argument similar to those used in many papers in the sparsification literature [19, 4, 12, 24, 13, 23]. As a consequence, we also get a much sharper asymptotic analysis.

Theorem 2 (Main Result, see Theorem 22).

For a set system 𝒮2[m], weights w:[m]0 and ϵ>0, there exists a (1±ϵ) sparsifier of 𝒮 which retains only

O(CL(𝒮)log2(CL(𝒮)/ϵ)(loglog(CL(𝒮)/ϵ))2ϵ2)

many atoms.

Theorem 2 improves on Theorem 1 in two ways. First, the multiplicative overhead of log6m is reduced to a much smaller log2+o(1)(CL(𝒮)/ϵ). Second, the bound in Theorem 2 is independent222We observe that the bound in Theorem 1 could also be made independent of m by applying a recursive argument similar to that of Section 4.2. Even so, the asymptotics of Theorem 2 would still be superior. of m!

As an application of the above theorem, one can consider a graph G=(V,E), and the set system 𝒮2E which contains all the cuts in the graph G. I.e., for every set TV, one can consider the set of edges ETE which is cut by the set T and let 𝒮=TVET. Sparsifying the set system 𝒮 is thus equivalent to creating a cut-sparsifier of the graph G in the classical sense of [4]. Because the chain length of the set of cuts is bounded by |V|, Theorem 2 recovers (up to an extra O~(log(m)) factor) the sparsifier size bound of [4].

Note that Theorem 2 is however not efficiently implementable. This is in large part due to the fact that, given a set system 𝒮2[m], there is no known efficient algorithm for computing CL(𝒮).

1.1.1 Improvement to Code Sparsification

As an additional corollary to Theorem 2, we give the first bound for linear code sparsification (see [23]) which is independent of the underlying finite field. In more detail, given a finite field 𝔽q on q elements, we define a linear code to be an n-dimensional subspace C𝔽qm. A (1±ϵ) sparsifier of C is a reweighting w:[m]0 such that for every cC we have that

i=1mw(i)𝟏[ci0](1±ϵ)i=1m𝟏[ci0].

A main result of Khanna, Putterman, and Sudan [23] (see also [21]) is that C has a (1±ϵ) sparsifier of size333Size here is meant to be to the number of non-zero weights that are assigned; equivalently, the number of coordinates in [m] that are retained. |supp(w)|=nlogO(1)(n)logq/ϵ2. If we capture linear code sparsification by letting our set system 𝒮 be {suppc:cC}, then we obtain a linear code sparsifier of size CL(𝒮)log2+o(1)(CL(𝒮)/ϵ)/ϵ2. From existing results on chain length [5, 7], we know that CL(𝒮) is precisely the dimension n of the code C. This gives the first field-size independent bound of nlog2+o(1)(n/ϵ)/ϵ2 for code sparsification!

Lastly, recall that both [7], Theorem 1 and [23] used Theorem 1 and linear code sparsification respectively to design CSP sparsifiers. Because of our better parameters in Theorem 2, we immediately obtain tighter bounds for these applications.

1.2 Related Work

We now briefly discuss how our work connects to other parts of the sparsification and broader TCS literature.

1.2.1 Code, CSP, and Cayley graph sparsification

As previously mentioned, Khanna, Putterman and Sudan [23] pioneered the concept of code sparsification. In their original paper, [23] presented multiple applications of code sparsification including constraint satisfaction problem (CSP) sparsification and Cayley graph sparsification.

CSP sparsification was introduced by Kogan and Krauthgamer [24] to generalize graph cut sparsification to broader families of discrete structures, including hypergraph cut sparsification (e.g., [28, 10, 18, 17, 25, 16, 20]). A more systematic investigation of general CSP sparsification was started by Filtser–Krauthgamer [11] and Butti–Živný [8]. However, the case of linear equations over a finite field was not resolved until the method of code sparsification was developed [23]. In a follow-up Khanna, Putterman and Sudan [21] extended these methods to equations over Abelian groups among other CSPs. However, results concerning all CSPs were not established until Brakensiek and Guruswami [7] generalized linear code sparsification to non-linear codes (i.e., arbitrary set systems).

In a separate direction, the work of Khanna, Putterman and Sudan [23] introduced Cayley graph sparsification as a special case of graph cut sparsification where the input graph is a Cayley graph (i.e., the graph is defined by the generators of a suitable group). Here, the goal is not to construct any sparsifier, but rather a sparsifier which is itself a Cayley graph. Khanna, Putterman and Sudan [23] handled the case in which the Cayley graph corresponds to spanning vectors of a finite vector space. More recently, such methods were generalized to arbitrary groups [15], and even a more general theory of sparsifying sums of PSD matrices [2].

1.2.2 Chain length and query complexity

The notion of chain length was introduced by Bessiere, Carbonnel, and Katsirelos [5] in the context of analyzing the query complexity of constraint satisfaction problems in a model introduced by Bessiere et al. [6]. In this query model, one is not given explicit access to the constraints; rather, one queries a partial assignment to a subset of the variables, and the response is “YES” or “NO” depending on whether the partial assignment is consistent with all constraints induced by those variables. To connect this problem with the concept of chain length, consider a set system 𝒮2[m] where [m] represents the clauses of the CSP and S𝒮 if and only if there is a valid assignment to the CSP which satisfies the clauses indexed by S (and no others). Then, O~(CL(𝒮)) is an upper bound on the query complexity of this problem. However, it is an open question whether this chain length bound is tight. See [7] for further discussion.

1.2.3 Connections to Learning Theory

Recall that chain length of 𝒮2[m] is the length of the longest ascending chain in 𝒮. Another natural property of a set family is its VC-dimension. That is, the size of the largest set A[m] such that for all BA, there is some S with SA=B. Bessiere, Carbonnel, and Katsirelos [5] define the non-redundancy of 𝒮, denoted by NRD(𝒮) to be precisely the VC-dimension of 𝒮. Brakensiek and Guruswami [7] identify the non-redundancy of 𝒮 is closely related to the sparsifiability of 𝒮, assuming 𝒮 is unweighted.

More broadly, the VC-dimension of a set family 2[m] is well-known to characterize the sample complexity needed to obtain a small additive-error approximation of the set-sizes in the family. That is, with constant probability a random subset T[m] of size O(VC-dim()/ϵ2) [29, 1] satisfies |TS||S|±ϵm for every S. Conversely, no set of size smaller than Ωϵ(VC-dim()) can satisfy such a guarantee.

In this view, the chain length characterizes the sample complexity needed to obtain multiplicative-error approximation to set-sizes, in the weighted setting. Note that, unlike the VC-theorem which works with unweighted samples, reweighting of the chosen samples is necessary to obtain a multiplicative approximation. We also remark that an unweighted multiplicative sampling theorem (relative to VC-dimension) was proved by Li, Long, and Srinivasan [26]; however, their notion of multiplicative approximation significantly differs from ours.

Organization

In Section 2, we state some known results (Chernoff bounds, etc.). In Section 3, we prove a key counting bound by relating the technique of contractions to chain length. In Section 4, we use the results of Section 3 to construct a series of sparsifiers, culminating in Theorem 2. In Section 5, we give some concluding thoughts and open questions.

2 Preliminaries

We now present some background material on sparsification and related concepts.

2.1 Notation and Sparsification Definitions

For the purposes of our analysis, we represent set systems 𝒮2[m] as a set of vectors (also called a code) C{0,1}m, where S𝒮 if and only if cC where c{0,1}m is defined by

ci={1iS0iS.

Optionally, the coordinates [m] of the code may be given weights w:[m]0. Note that when no weights are provided for the code, we assume that the weights are all 1 (equivalently, an unweighted code). Occasionally, for a code C{0,1}m and a set S[m], we will use the notation C|S={c|S:cC} to be the coordinate restriction of all codewords in C to the set S.

Our goal, given a parameter ϵ(0,1) is to design a sparsifier:

Definition 3.

For a code C{0,1}m, weights w:[m]0, and an accuracy parameter ϵ>0, a (1±ϵ) sparsifier of C is a new set of weights w~:[m]0 such that, for every cC:

w~,c=i=1mw~(i)ci(1±ϵ)w,c=i=1mw(i)ci.

The goal is to design sparsifiers which reduce the support size of the starting code; i.e., which minimize |supp(w~)|.

 Remark 4.

Note that in some contexts, for an unweighted code C{0,1}m, we will use wC to refer to the code which assigns weight w to all coordinates in [m]. Similarly, for codes C1{0,1}A1, C2{0,1}A2 with A1A2=, we may use w1C1w2C2 to refer to the code which assigns weights w1 to coordinates in A1 and w2 to coordinates in A2.

For reference, we include the definitions of non-redundancy and chain length below. Note that non-redundancy essentially captures the largest diagonal matrix (after permuting the codewords) which one can find in a set of {0,1}-valued vectors.

Definition 5.

Let C{0,1}m be an arbitrary set of vectors. We say that the non-redundancy of C (denoted NRD(C)) is the largest size of a set S[m] such that for every jS, there is a codeword cC where cj=1, but for every other iS{j}, ci=0.

Likewise, we use the notion of chain length, which is essentially the largest upper triangular submatrix contained in a set of vectors:

Definition 6.

Let C{0,1}m be a code. A chain of length is a pair of injective maps a:[]m and c:[]C such that the following conditions hold:

  1. 1.

    i[]:c(i)a(i)=1

  2. 2.

    1i<j:c(i)a(j)=0.

The chain length of C, denoted by CL(C) is the length of the longest chain.

Note that NRD(C)CL(C) as any diagonal matrix is trivially an upper triangular matrix. With this, we can make use of the following bound (observed in [7]):

Claim 7.

|C|(m+1)NRD(C)(m+1)CL(C).

2.2 Concentration Bound

We will make use of the following concentration bound when sub-sampling to construct our sparsifiers:

Claim 8.

([12]) Let X1,X be random variables such that Xi takes on value 1/pi with probability pi, and is 0 otherwise. Also, suppose that minipip. Then, with probability at least 12e0.38ϵ2p,

iXi(1±ϵ).

3 Chain Length Counting Bounds via Contractions

In this section we show the following lemma via a simple contraction argument:

Lemma 9.

Let C{0,1}m be any code. Then, for any parameter d>0, there is a set T[m] of size |T|CL(C)d such that for any α+, the number of codewords in C|T¯ of weight αd is at most (CL(C)α)(m+1)α.

3.1 Contractions

To start, we have the following claim which governs how the chain length of a set of vectors behaves under contractions.

Claim 10.

Let C{0,1}m be an arbitrary set of vectors and let i[m] be any coordinate such that there exists a cC such that ci=1. Then, for C={cC:ci=0}, we have that CL(C)CL(C)1.

Proof.

For the code C, let the maps witnessing the chain length of C be denoted by a,c, and let the chain length be . Now, let us define a:[+1][m],c:[+1]C, such that for j[],c(j)=c(j),a(j)=a(j), and a(+1)=i and c(+1)=v, where vC is any codeword such that vi0.

Observe that we trivially have j[+1]:c(j)a(j)=1. Likewise, because every codeword c in C satisfies ci=0 and the first codewords mapped to by c are all in C, we also have that for any j[], c(j)a(+1)=c(j)i=0. Together with the conditions already guaranteed by c,a, this yields a chain of length +1 in C, and thus the claim.

3.2 Contraction Algorithm

Leveraging the above claim which governs chain length under contractions, we define the following iterative contraction procedure:

Algorithm 1 Contract(C,α).

Now, let us define a new quantity:

Definition 11.

The density of a code C{0,1}m is given by

Φ(C)=minCC|Supp(C)|CL(C),

where Supp(C)={i[m]:cC:ci0}.

We immediately get the following counting bound.

Claim 12.

Let C{0,1}m be a code. Then, for any positive integer α, the number of codewords of weight αΦ(C) is (CL(C)α)(m+1)α.

Proof.

Fix any codeword c of weight αΦ(C), and run Contract(C,α). For some intermediate code C achieved during the contraction procedure, observe that if CL(C)=k, then necessarily, |Supp(C)|Φ(C)k. Thus, in the next random choice of i, the probability that our codeword c satisfies ci=1 is at most

wt(c)Φ(C)kαΦ(C)Φ(C)kαk,

and thus c survives the contraction with probability 1αk. The probability that c survives all contractions until CL(C)α is thus

k=α+1CL(C)(1αk)=(CL(C)α)1.

Conditioned on surviving, the codeword c is then chosen uniformly at random among all surviving codewords, and thus returned by Contract(C,α) with probability (m+1)α, where we are using Claim 7 which bounds the number of codewords in a code with CL at most α. So, c is returned by Contract(C,α) with probability at least (m+1)α(CL(C)α)1. By taking the reciprocal, this yields the bound on the number of codewords of weight αΦ(C).

Importantly, we also have the following claim:

Claim 13.

Let C{0,1}m be a code, and let CC such that CL(C)=. Let T=Supp(C), and let C|T¯ denote the code C with the coordinates corresponding to the support of C removed. We have:

CL(C|T¯)CL(C).
Proof.

As before, let a1,c1 be the functions which witness the chain length of C|T¯, and let a2,c2 be the functions which witness the chain length of CL(C). Note that Im(a1)[m]T[m] and Im(a2)T[m].

We define functions a,c as follows: a:[+CL(C|T¯)][m], where a[j],c[j]=a2[j],c2[j] if j, and otherwise a[j],c[j]=a1[j],c1[j]. Importantly, because for j[], c(j)C, and a()T¯, we have that c(j)a()=0. The remaining conditions for chain length are then trivially satisfied by the fact that a1,c1,a2,c2 already satisfied the definition of chain length individually.

Pictorially, we have:

C=[CA0C|T¯],

thus we can trivially compose the upper triangular matrices corresponding to C,C|T¯.

Lemma 14.

Let C{0,1}m be any code. Then, for any parameter d>0, there is a set T[m] of size |T|CL(C)d such that for any α+, the number of codewords in C|T¯ of weight αd is at most (CL(C)α)(m+1)α.

Proof.

If Φ(C)>d, then we are immediately done by invoking Claim 12. Otherwise, there is some CC for which

|Supp(C)|CL(C)d. (1)

So, let T=Supp(C), and let us set C=C|T¯. Again, we can check whether Φ(C)>d, if so, we are done, and we simply return the set T as is. Otherwise, there is a new subcode C′′C for which |Supp(C′′)|CL(C′′)d. Again, we let TTSupp(C′′), let C=C|T¯ and continue on recursively.

Eventually, Φ(C)>d (or the entire support of the code is removed) at which point we terminate. All that remains is to bound the size of the set T. For this, every time we remove the support of a subcode C, we increase the size of the set T by at most dCL(C), as |Supp(C)|dCL(C) (by Equation 1). However, by Claim 13, every time we remove Supp(C) from C, CL(C) decreases by CL(C). Hence, the total number of coordinates removed in the set T can be at most dCL(C), before the chain length of the remaining code goes to 0. This yields the lemma.

4 Leveraging Decomposition to Build Sparsifiers

4.1 Basic Sparsifier Construction for Unweighted Codes

A consequence of Lemma 9 is that it allows us to build sparsifiers for unweighted codes which preserve only O~(CL(C)) many coordinates:

Theorem 15.

Let C{0,1}m be an unweighted code. Then, for any ϵ>0, there exists a (1±ϵ) sparsifier of C which preserves only

O(CL(C)(logm)2(loglogm)2/ϵ2)

many re-weighted coordinates.

The algorithm for achieving this sparsification is simple, and builds off of the intuition of [23, 21]. Namely, starting with a code C{0,1}m, we invoke the decomposition of Lemma 9 with parameter d=mCL(C). This yields two codes: CPeel which contains the CL(C)dmCL(C) many coordinates that Lemma 9 peels off, along with Cremaining, which contains all remaining coordinates, with the promise that in Cremaining, for any α+, the number of codewords in (Cremaining)|T¯ of weight αd is at most (CL(C)α)(m+1)α. We then sub-sample the coordinates of Cremaining, and then recurse. The formal algorithm is presented below.

Algorithm 2 Sparsify(C,ϵ,counter,m).

We now prove some basic claims regarding the above algorithm:

Claim 16.

When invoked on a code C{0,1}m, Line 7 in Algorithm 2 is always possible.

Proof.

After performing the decomposition of Lemma 9 with parameter d=mηCL(C), it must be the case that in Cremaining, for any α+, the number of codewords in (Cremaining) of weight αd is at most (CL(C)α)(m+1)α. Thus, if one randomly samples the coordinates of Cremaining at rate p=ηCL(C)m, and assigns weight 1/p to the sampled coordinates a simple Chernoff and union bound will show that every codeword’s weight is preserved to a (1±ϵ/20loglog(m)) factor with probability 11/poly(m). Indeed, for a fixed value of α, a codeword of weight [αd,2αd] will have its weight preserved to within a (1±ϵ/20loglog(m)) factor with probability (using Claim 8) at least

12e0.38(ϵ/20loglog(m))2αdηCL(C)m 12e0.38(ϵ/20loglog(m))2αη
12e0.38α1000log(m)
11m100α.

Then, we can take a union bound over all codewords of weight 2αd, at most (CL(C)2α)(m+1)2α(m+1)4α of them, and then a union bound over the at most m choices of α to conclude that every codeword has its weight preserved to a (1±ϵ/(20loglog(m))) factor with probability 11m95 (in particular such a sparsification exists).

At the same time, a Chernoff bound will show that when sampling at this rate, with probability Ω(1), at most 2CL(C)ηm many coordinates will survive the sampling, and thus there exists a set of coordinates |T~|2CL(C)ηm, and weight w such that wCT~ is a (1±ϵ/(20loglog(n)) code sparsifier of CRemaining.

Claim 17.

When invoked on a code C{0,1}m, wCT~CPeel as produced in Line 7 is a (1±ϵ/20loglog(m)) code sparsifier of C.

Proof.

This follows because CT¯CPeel is a perfect sparsifier of C. Then, when replacing CT¯ with a (1±ϵ/20loglog(m)) sparsifier of CT¯, it follows by composition that the resulting code is a (1±ϵ/20loglog(m)) code sparsifier of C.

Claim 18.

Let C{0,1}m be a code and let ϵ>0. Then Sparsify(C,ϵ,0) returns a code C~ which is a (1±ϵ) sparsifier of C.

Proof.

This follows by induction over the value of Counter. Indeed, when sparsifying the code C, we consider all the resulting codes for which Sparsify is recursively called. We let wi,1C1(i),wi,2iC2i(i) denote all 2i codes for which the Sparsify function is called with Counter=i. We inductively claim that wi,1C1(i)wi,2iC2i(i) is a (1±3iϵ/20loglog(m)) sparsifier of C.

The base case follows trivially, using Claim 17.

The inductive case follows by observing that for each code wi,jCj(i), by invoking Claim 17, we replace Cj(i) by two codes wi+1,2jC2j(i+1) and wi+1,2j+1C2j+1(i) such that

wi+1,2jC2j(i+1)wi+1,2j+1C2j+1(i)

is a (1±ϵ/20loglog(m)) sparsifier of wi,jCj(i). By composition over all the codes wi,1C1(i)wi,2iC2i(i), it follows then that

wi+1,1C1(i+1)wi+1,2i+1C2i+1(i+1)

will be a (1±ϵ/20loglog(m))(1±3ϵi/20loglog(m))(1±3ϵ(i+1)/20loglog(m))-sparsifier.

Taking i=loglog(m) yields our desired claim.

Claim 19.

Let C{0,1}m be a code and let ϵ>0. Then Sparsify(C,ϵ,0) returns a code C~ which retains only O(CL(C)log2(m)(loglog(m))2ϵ2) many coordinates.

Proof.

We claim that whenever Sparsify is invoked on a code C with Counter=i, then

|Supp(C)|4CL(C)(mCL(C))1/2iη.

The base case follows trivially: indeed, in the first invocation (i=1) of the algorithm, |Supp(CPeel)|CL(C)mηCL(C) which trivially satisfies the above bound, and |Supp(CT~)|2CL(C)mη2CL(C)mCL(C)η, which also trivially satisfies the above bound. Now, we assume the claim holds by induction.

Then, at level i of the algorithm

|Supp(C)|4CL(C)(mCL(C))1/2iη.

When we sparsify C, we obtain two codes, CPeel and CT~. We have that (by our choice of d in Algorithm 2)

|Supp(CPeel)|CL(C)4CL(C)(mCL(C))1/2iηηCL(C)4CL(C)η(mCL(C))1/2i+1,

as we desire. Likewise,

|Supp(CT~)|2CL(C)4CL(C)(mCL(C))1/2iηη4CL(C)η(mCL(C))1/2i+1
4CL(C)η(mCL(C))1/2i+1,

as we desire.

Now, once i=loglog(m), we have that

|Supp(C)|4CL(C)(mCL(C))1/2loglog(m)η8CL(C)η.

Thus, the returned code C~ is the union of 2loglog(m)=log(m) many codes, each of which retains 8CL(C)η many coordinates. Thus, the total sparsifier size is bounded by

log(m)8CL(C)η=O(CL(C)log2(m)(loglog(m))2ϵ2),

as we desire.

Together, these above claims give the following theorem:

Theorem 20.

For a code C{0,1}m and ϵ>0, there exists a (1±ϵ) code sparsifier of C which retains only O(CL(C)log2(m)(loglog(m))2ϵ2) many coordinates.

Proof.

We invoke Algorithm 2. The accuracy of the sparsifier follows from Claim 18 and the size follows from Claim 19.

Note that at this point, one can use Theorem 20 along with weighted to unweighted sparsifier reduction frameworks as in [23, 21, 7] to deduce sparsifiers for weighted codes which preserve O(CL(C)log2(m)(loglog(m))2ϵ2) many coordinates.

Corollary 21.

For a code C{0,1}m, weights w:[m]0 and ϵ>0, there exists a (1±ϵ) code sparsifier of C which retains only O(CL(C)log2(m)(loglog(m))2ϵ2) many coordinates.

For completeness, we include a proof of Corollary 21 in Appendix A.

4.2 Dimension-Free Bounds

In this section, we show that the dependence on log(m) can be removed, and instead replaced with dependence only on log(CL), thus constituting a dimension-free sparsification result. To do this, we merely repeatedly apply the sparsification result of Corollary 21:

Theorem 22.

For a code C{0,1}m, weights w:[m]0 and ϵ>0, there exists a (1±ϵ) code sparsifier of C which retains only

O(CL(C)log2(CL(C)/ϵ)loglog(CL(C)/ϵ)2ϵ2)

many coordinates.

Proof.

First, let us use K to denote the hidden constant in the O() notation above. In the below, we assume that m is at least a sufficiently large constant to begin with. Now, we consider two cases:

  1. 1.

    When m2CL(C)/(ϵ2), then after applying Corollary 21 with ϵ=ϵ/Qlogm (for a large constant Q), we obtain a (1±ϵ) sparsifier of C with support KCL(C)log4(m)loglog(m)2ϵ2log6(m).

  2. 2.

    When m<2CL(C)/(ϵ2), we sparsify C once with parameter ϵ=ϵ/Qlog(m), obtaining a code C which is a (1±ϵ) sparsifier of C with support KCL(C)log4(m)loglog(m)2ϵ2KCL(C)7/ϵ14. Then, we sparsify C with parameter ϵ/2, and obtain a code C′′ which is a (1±ϵ/2) sparsifier of C, and retains at most

    KCL(C)log2(CL(C)/ϵ)loglog(CL(C)/ϵ)2ϵ2

    many coordinates.

Thus, given an arbitrary codes C, we apply case (1) =O(log(m)) many times, yielding code C(1),C(). Then, we apply case (2) to C(), yielding a code C′′ which retains only

KCL(C)log2(CL(C)/ϵ)loglog(CL(C)/ϵ)2ϵ2

many coordinates.

All that remains is to show that C′′ is a (1±ϵ) sparsifier of C. This follows by composition of the sparsifier accuracy: we let m(i) denote the support size of the ith sparsifier we create in the above chain (with m(0) being the starting value of m). Then, C(1) is a (1±ϵ/Qlog(m(0))) sparsifier of C, C(2) is a (1±ϵ/log(m(1))) sparsifier of C(1), and so on. In particular, this implies that C() is a

i=01(1±ϵ/Qlog(m(i)))

sparsifier of C. Thus C() is a

(1±ϵO(i=011Qlog(m(i))))=(1±ϵ/10)

sparsifier of C, where we have chosen Q to be a large enough constant. In particular, here we are using the fact that, letting q=log(m(1)),

i=011log(m(i))1q+12q1/6+12(2q1/6)1/6+=O(1/q)=O(1).

This relation follows from the fact that m(i)log6(m(i1)) as established above.

Now C is a (1±ϵ/Q) sparsifier of C(), which implies that C is a (1±ϵ/10)(1±ϵ/Q)(1±ϵ/4) sparsifier of C.

Finally, C′′ is a (1±ϵ/2) sparsifier of C, and thus by composition, C′′ is a (1±ϵ) sparsifier of C, as we desire (note, here we are using that ϵ is a sufficiently small constant).

5 Conclusion

In this paper, we constructed sparsifiers for arbitrary set systems of near-optimal size. In particular, we simplify the proof strategy of Brakensiek and Guruswami [7] by presenting a contraction-based counting bound which results in an overall simpler proof in addition to better asymptotics. Furthermore, by recursively applying our sparsifier, our bound only depends on the accuracy ϵ and the chain length of the underlying set system. We conclude this paper with an exciting open direction and our conjecture regarding it.

Truly linear sparsifiers.

A celebrated result of Batson, Spielman and Srivastava [3] states that graphs on n vertices have (1±ϵ) cut sparsifiers with O(n/ϵ2)-edges. Furthermore, such asymptotics are known to be tight up to a constant factor [9]. Proving an analogue of [3] for arbitrary set systems seems quite difficult as there is no clear spectral analogue of set system sparsification (however, see [22]). Currently, the most general extensions of [3] are to the regime of PSD matrix sparsification in the work of [27] (which even includes arity 3 hypergraph cut sparsification). Nonetheless, we boldly conjecture that analogous constructions should exist in general.

Conjecture 23.

For a set system 𝒮2[m], weights w:[m]0 and ϵ>0, there exists a (1±ϵ) sparsifier of 𝒮 which retains only

O(CL(𝒮)ϵ2)

many coordinates.

Of note, even proving Conjecture 23 for linear code sparsification would be a significant result.

References

  • [1] Martin Anthony and Peter L. Bartlett. Neural Network Learning: Theoretical Foundations. Cambridge University Press, Cambridge, 1999.
  • [2] Arpon Basu, Pravesh K. Kothari, Yang P. Liu, and Raghu Meka. Sparsifying sums of positive semidefinite matrices. In Kasper Green Larsen and Barna Saha, editors, Proceedings of the 2026 Annual ACM-SIAM Symposium on Discrete Algorithms, SODA 2026, Vancouver, BC, Canada, January 11-14, 2026, pages 6042–6064. SIAM, 2026. doi:10.1137/1.9781611978971.216.
  • [3] Joshua D. Batson, Daniel A. Spielman, and Nikhil Srivastava. Twice-Ramanujan sparsifiers. In Michael Mitzenmacher, editor, Proceedings of the 41st Annual ACM Symposium on Theory of Computing, STOC 2009, Bethesda, MD, USA, May 31 - June 2, 2009, pages 255–262. ACM, 2009. doi:10.1145/1536414.1536451.
  • [4] András A. Benczúr and David R. Karger. Approximating s-t minimum cuts in Õ(n2) time. In Gary L. Miller, editor, Proceedings of the Twenty-Eighth Annual ACM Symposium on the Theory of Computing, Philadelphia, Pennsylvania, USA, May 22-24, 1996, pages 47–55. ACM, 1996. doi:10.1145/237814.237827.
  • [5] Christian Bessiere, Clément Carbonnel, and George Katsirelos. Chain Length and CSPs Learnable with Few Queries. Proceedings of the AAAI Conference on Artificial Intelligence, 34(02):1420–1427, April 2020. doi:10.1609/aaai.v34i02.5499.
  • [6] Christian Bessiere, Remi Coletta, Emmanuel Hebrard, George Katsirelos, Nadjib Lazaar, Nina Narodytska, Claude-Guy Quimper, and Toby Walsh. Constraint acquisition via partial queries. In Francesca Rossi, editor, IJCAI 2013, Proceedings of the 23rd International Joint Conference on Artificial Intelligence, Beijing, China, August 3-9, 2013, pages 475–481. IJCAI/AAAI, 2013. URL: http://www.aaai.org/ocs/index.php/IJCAI/IJCAI13/paper/view/6659.
  • [7] Joshua Brakensiek and Venkatesan Guruswami. Redundancy is all you need. In Michal Koucký and Nikhil Bansal, editors, Proceedings of the 57th Annual ACM Symposium on Theory of Computing, STOC 2025, Prague, Czechia, June 23-27, 2025, pages 1614–1625. ACM, 2025. doi:10.1145/3717823.3718212.
  • [8] Silvia Butti and Stanislav Živný. Sparsification of Binary CSPs. SIAM Journal on Discrete Mathematics, 34(1):825–842, January 2020. doi:10.1137/19M1242446.
  • [9] Charles Carlson, Alexandra Kolla, Nikhil Srivastava, and Luca Trevisan. Optimal lower bounds for sketching graph cuts. In Proceedings of the Thirtieth Annual ACM-SIAM Symposium on Discrete Algorithms, pages 2565–2569. SIAM, 2019. doi:10.1137/1.9781611975482.158.
  • [10] Yu Chen, Sanjeev Khanna, and Ansh Nagda. Near-linear size hypergraph cut sparsifiers. In Sandy Irani, editor, 61st IEEE Annual Symposium on Foundations of Computer Science, FOCS 2020, Durham, NC, USA, November 16-19, 2020, pages 61–72. IEEE, 2020. doi:10.1109/FOCS46700.2020.00015.
  • [11] Arnold Filtser and Robert Krauthgamer. Sparsification of two-variable valued constraint satisfaction problems. SIAM J. Discret. Math., 31(2):1263–1276, 2017. doi:10.1137/15M1046186.
  • [12] Wai Shing Fung, Ramesh Hariharan, Nicholas J.A. Harvey, and Debmalya Panigrahi. A general framework for graph sparsification. In Proceedings of the Forty-Third Annual ACM Symposium on Theory of Computing, STOC ’11, pages 71–80, New York, NY, USA, 2011. Association for Computing Machinery. doi:10.1145/1993636.1993647.
  • [13] Mohsen Ghaffari, David R. Karger, and Debmalya Panigrahi. Random contractions and sampling for hypergraph and hedge connectivity. In Philip N. Klein, editor, Proceedings of the Twenty-Eighth Annual ACM-SIAM Symposium on Discrete Algorithms, SODA 2017, Barcelona, Spain, Hotel Porta Fira, January 16-19, pages 1101–1114. SIAM, 2017. doi:10.1137/1.9781611974782.71.
  • [14] Justin Gilmer. A constant lower bound for the union-closed sets conjecture. arXiv preprint, 2022. arXiv:2211.09055.
  • [15] Jun-Ting Hsieh, Daniel Z. Lee, Sidhanth Mohanty, Aaron Putterman, and Rachel Yun Zhang. Sparsifying cayley graphs on every group. In Kasper Green Larsen and Barna Saha, editors, Proceedings of the 2026 Annual ACM-SIAM Symposium on Discrete Algorithms, SODA 2026, Vancouver, BC, Canada, January 11-14, 2026, pages 6029–6041. SIAM, 2026. doi:10.1137/1.9781611978971.215.
  • [16] Arun Jambulapati, Yang P Liu, and Aaron Sidford. Chaining, group leverage score overestimates, and fast spectral hypergraph sparsification. In Proceedings of the 55th Annual ACM Symposium on Theory of Computing, pages 196–206, 2023. doi:10.1145/3564246.3585136.
  • [17] Michael Kapralov, Robert Krauthgamer, Jakab Tardos, and Yuichi Yoshida. Spectral hypergraph sparsifiers of nearly linear size. In 62nd IEEE Annual Symposium on Foundations of Computer Science, FOCS 2021, Denver, CO, USA, February 7-10, 2022, pages 1159–1170. IEEE, 2021. doi:10.1109/FOCS52979.2021.00114.
  • [18] Michael Kapralov, Robert Krauthgamer, Jakab Tardos, and Yuichi Yoshida. Towards tight bounds for spectral sparsification of hypergraphs. In Proceedings of the 53rd Annual ACM SIGACT Symposium on Theory of Computing, pages 598–611, 2021. doi:10.1145/3406325.3451061.
  • [19] David R. Karger. Global min-cuts in rnc, and other ramifications of a simple min-cut algorithm. In Vijaya Ramachandran, editor, Proceedings of the Fourth Annual ACM/SIGACT-SIAM Symposium on Discrete Algorithms, 25-27 January 1993, Austin, Texas, USA, pages 21–30. ACM/SIAM, 1993. URL: http://dl.acm.org/citation.cfm?id=313559.313605.
  • [20] Sanjeev Khanna, Aaron Putterman, and Madhu Sudan. Near-optimal size linear sketches for hypergraph cut sparsifiers. In 2024 IEEE 65th Annual Symposium on Foundations of Computer Science (FOCS), pages 1669–1706. IEEE, 2024. doi:10.1109/FOCS61266.2024.00105.
  • [21] Sanjeev Khanna, Aaron Putterman, and Madhu Sudan. Efficient algorithms and new characterizations for CSP sparsification. In Michal Koucký and Nikhil Bansal, editors, Proceedings of the 57th Annual ACM Symposium on Theory of Computing, STOC 2025, Prague, Czechia, June 23-27, 2025, pages 407–416. ACM, 2025. doi:10.1145/3717823.3718205.
  • [22] Sanjeev Khanna, Aaron Putterman, and Madhu Sudan. A theory of spectral CSP sparsification. arXiv preprint, 2025. doi:10.48550/arXiv.2504.16206.
  • [23] Sanjeev Khanna, Aaron (Louie) Putterman, and Madhu Sudan. Code sparsification and its applications. In David P. Woodruff, editor, Proceedings of the 2024 ACM-SIAM Symposium on Discrete Algorithms, SODA 2024, Alexandria, VA, USA, January 7-10, 2024, pages 5145–5168. SIAM, 2024. doi:10.1137/1.9781611977912.185.
  • [24] Dmitry Kogan and Robert Krauthgamer. Sketching cuts in graphs and hypergraphs. In Tim Roughgarden, editor, Proceedings of the 2015 Conference on Innovations in Theoretical Computer Science, ITCS 2015, Rehovot, Israel, January 11-13, 2015, pages 367–376. ACM, 2015. doi:10.1145/2688073.2688093.
  • [25] James R. Lee. Spectral hypergraph sparsification via chaining. In Barna Saha and Rocco A. Servedio, editors, Proceedings of the 55th Annual ACM Symposium on Theory of Computing, STOC 2023, Orlando, FL, USA, June 20-23, 2023, pages 207–218. ACM, 2023. doi:10.1145/3564246.3585165.
  • [26] Yi Li, Philip M. Long, and Aravind Srinivasan. Improved bounds on the sample complexity of learning. Journal of Computer and System Sciences, 62(3):516–527, 2001. doi:10.1006/jcss.2000.1741.
  • [27] Marcel K De Carli Silva, Nicholas JA Harvey, and Cristiane M Sato. Sparse sums of positive semidefinite matrices. ACM Transactions on Algorithms (TALG), 12(1):1–17, 2015. doi:10.1145/2746241.
  • [28] Tasuku Soma and Yuichi Yoshida. Spectral sparsification of hypergraphs. In Proceedings of the Thirtieth Annual ACM-SIAM Symposium on Discrete Algorithms, pages 2570–2581. SIAM, 2019. doi:10.1137/1.9781611975482.159.
  • [29] V. N. Vapnik and A. Y. Chervonenkis. On the uniform convergence of relative frequencies of events to their probabilities. Theory of Probability and Its Applications, 16:264–280, 1971.

Appendix A Proof of Corollary 21

The specific method we use for proving Corollary 21 closely follows the methodology in Section 8 of [7] (cf. Theorem 8.4). To start, we prove a version of this corollary in the setting where all the weights are bounded:

Corollary 24.

For a code C{0,1}m, weights w:[m]0 such that w(i)m3 and ϵ>0, there exists a (1±ϵ) code sparsifier of C which retains only O(CL(C)log2(m)(loglog(m))2ϵ2) many coordinates.

Proof.

We adapt the proof of Lemma 8.15 in [7]. First, if ϵ1/m, then we simply return all the coordinates of the code C, as m1/ϵ2. So, for the rest of the proof, we assume that ϵ>1/m. We will also assume that miniw(i)=1, as otherwise we can simply rescale the weights by miniw(i) to ensure this is the case (if the min weight is 0, we can simply delete these coordinates).

The key intuition is now that the weights are bounded, we will simply duplicate the coordinates of the code a number of times proportional to their weight, thus transforming the weighted code into an unweighted code. With this unweighted code, we can then invoke Theorem 20.

To make this unweighted code for each i[m], we define b(i)=2w(i)ϵ2m4 to be the number of times we duplicate coordinate i. Thus, we create a new code C~{0,1}m~ where m~=i=1mb(i) where coordinate i of C is replaced with b(i) unweighted copies in C~. Now, we claim that for any codeword cC, its corresponding version c~ satisfies ϵ2i=1m~c~i(1±ϵ/2)i=1mw(i)ci. Indeed, this follows because

ϵ2i=1m~c~i=i=1mϵ22w(i)ϵcii=1mϵ2[2w(i)ϵ1,2w(i)ϵ]ci
i=1m[1ϵ/2,1]wici(1±ϵ/2)i=1mw(i)ci.

Thus, to conclude the proof, we must only construct this now unweighted code C~, sparsify it to accuracy (1±ϵ/3) using Theorem 20 to create a new, weighted code C^, and then return ϵ2C^. By the composition of sparsifier accuracy, this will be a (1±ϵ/2)(1±ϵ/3) sparsifier of C, which constitutes a (1±ϵ) sparsifier. The number of coordinates retained in the sparsifier is bounded by O(CL(C)log2(m)(loglog(m))2ϵ2), as we are sparsifying a code with m~=poly(m) many coordinates.

With this corollary established, we can now proceed to creating sparsifiers for codes with arbitrary weights.

Proof of Corollary 21.

We fix the code C{0,1}m, the weights w:[m]>0, and ϵ(0,1). We assume WLOG that CL(C)1 and that ϵ8/m (otherwise, the code C already satisfies the desired size bound).

We group the coordinates [m] according to their weights. Borrowing notation from the proof of Theorem 8.4 of [7], we define a function t:[m] such that

t(i):=logw(i)3logm,

and for each coordinate i[m], we assign its group to be t(i). For all t, we let It[n] denote the set of all coordinates in group t. For each codeword cC, we then say that the type of the codeword is the maximum group for any of its nonzero coordinates, i.e.,

type(c):=maxisupp(c)t(i).

For all t, we let Ct denote the set of all cC with type t.

The key observation, as used in [23, 21, 7] is that one must only consider sparsifying a codeword cCt with respect to the coordinates It; this is because the vast majority of its weight is contained in these coordinates. Furthermore, within the set of coordinates It, all weights are related by a factor of at most m3.

Define t to be proper if there exists z(CtCt+1)|It which has at least one coordinate equal to 1. For proper t, let w~t:[m]0 be an (ϵ/2)-sparsifier of (CtCt+1)|It{𝟏[supp(It)]}. Note that CL(CtCt+1)|It{𝟏[supp(It)]})1, so

CL((CtCt+1)|It{𝟏[supp(It)]}) CL(CtCt+1)|It{𝟏[supp(It)]})+1
2CL(CtCt+1)|It{𝟏[supp(It)]}).

Thus, by Corollary 24, we can ensure that

|supp(w~t)|O(CL((CtCt+1)|It)log2(m)loglog(m)2/ϵ2). (2)

Let w~:[m]0 be the sum of these sparsifiers {w~t:t,t proper}. This is well-defined as any coordinate appears in exactly one It and thus can be nonzero for only w~t. It remains only to show that w~ is indeed a (1±ϵ) sparsifier for C with weights w and that it satisfies our desired sparsity bound. Adapting a computation of [7], we have for all cC with t:=type(c) that

w~,c =w~t1,c|It1+w~t,c|It+tt2w~t,c|It
[1ϵ/2,1+ϵ/2]iIt1Itw(i)ci+tt2w~t,c|It
=[1ϵ/2,1+ϵ/2][w,cisupp(c)(It1It)w(i)]+tt2w~t,c|It.

First, we argue that tt2w~t,c|It is small. For this, suppose for the sake of contradiction that for some tt2, there is a choice of iIt such that w~(i)2m3t+4. Then, observe that for the codeword 𝟏[supp(It)], its weight reported by the sparsifier for It is at least 2m3t+4, which is 2iItw(i), as each w(i)m3t+3 (and |It|m). Thus, for this codeword 𝟏[supp(It)], w~ restricted to It (i.e,. w~t) would not have been a (1±ϵ) sparsifier. So, we see that w~(i)2m3t+4 for each iIt, and thus

tt2w~t,c|Itm2m3t+42m3t+52m3t1.

Since t=type(c), we have that w,cm3t. This immediately implies that tt2w~t,c|It2mw,cϵ4w,c since ϵ8/m.

Further, for all isupp(c)(It1It), we have that w(i)<m3t3. Therefore, the total contribution of w(i) for isupp(c)(It1It) is at most mm3t3ϵ4w,c. So,

w~,c [1ϵ/2,1+ϵ/2][1ϵ/4,1]w,c+[0,ϵ/4]w,c[1ϵ,1+ϵ]w,c,

as desired. It thus suffices to bound |supp(w~)|. By (Equation 2), we have that

|supp(w~)|t proper|supp(w~t)|tO(CL((CtCt+1)|It)log2(m)loglog(m)2/ϵ2). (3)

Analogous to Claim 8.16 of [7], we claim that

tCL((CtCt+1)|It)2CL(C).

To see this, let A1,AI1,I and B1,B(C1C2),(C2C3)(C,C+1) denote witnesses for the chain lengths of CL((CtCt+1)|It):t[]. The key observation is that Aodd=A1A3A5,Bodd=B1B3B5 and Aeven=A2A4A6,Beven=B2B4B6 are both valid witnesses of the chain length of C. To see why, we observe that because BiCiCi+1, it must be the case that every codeword cBi satisfies c|Ii+2=0. Thus, for any coordinate pAi+2, we see that cp=0. Importantly then, because Ai,Bi is a chain, and Ai+2,Bi+2 is defined on a disjoint support, their concatenation still remains a chain as per Claim 13.

With this, we then immediately have that

tCL((CtCt+1)|It)=t:t=0mod2CL((CtCt+1)|It)+t:t=1mod2CL((CtCt+1)|It)
CL(C)+CL(C)2CL(C),

and thus by (Equation 3) we can conclude that

|supp(w~)|O(CL(C)log2(m)loglog(m)2/ϵ2),

as desired.