Abstract 1 Introduction 2 Preliminaries 3 Non-explicit lower bounds for spiky rank 4 Connection to matrix rigidity 5 Connection to circuit complexity 6 Explicit lower bounds for spiky rank References

Spiky Rank and Its Applications to Rigidity and Circuits

Lianna Hambardzumyan The University of Copenhagen, Denmark    Konstantin Myasnikov EPFL, Lausanne, Switzerland    Artur Riazanov ORCID EPFL, Lausanne, Switzerland    Morgan Shirley Lund University, Sweden    Adi Shraibman The Academic College of Tel Aviv-Yaffo, Israel
Abstract

We introduce spiky rank, a new matrix parameter that enhances blocky rank by combining the combinatorial structure of the latter with linear-algebraic flexibility. A spiky matrix is block-structured with diagonal blocks that are arbitrary rank-one matrices, and the spiky rank of a matrix is the minimum number of such matrices required to express it as a sum. This measure extends blocky rank to real matrices and is more robust for problems with both combinatorial and algebraic character.

Our conceptual contribution is as follows: we propose spiky rank as a well-behaved candidate matrix complexity measure and demonstrate its potential through applications. We show that large spiky rank implies high matrix rigidity, and that spiky rank lower bounds yield lower bounds for depth-2 ReLU circuits, the basic building blocks of neural networks. On the technical side, we establish tight bounds for random matrices and develop a framework for explicit lower bounds, applying it to Hamming distance matrices and spectral expanders. Finally, we relate spiky rank to other matrix parameters, including blocky rank, sparsity, and the γ2-norm.

Keywords and phrases:
blocky rank, matrix rigidity, ReLU circuits, spiky rank
Category:
Track A: Algorithms, Complexity and Games
Funding:
Lianna Hambardzumyan: This work is funded by the European Research Council (ERC) under grant agreement no. 101125652 (ALBA). Most of the work was done while the author was a postdoctoral researcher at the University of Victoria, Canada, funded by NSERC and at Hebrew University of Jerusalem, Israel, funded by ISF grants 921/22 and 2635/19.
Konstantin Myasnikov: Supported by Swiss State Secretariat for Education, Research, and Innovation (SERI) under contract number MB22.00026.
Artur Riazanov: Supported by Swiss State Secretariat for Education, Research, and Innovation (SERI) under contract number MB22.00026.
Morgan Shirley: Supported by Knut and Alice Wallenberg grant KAW 2023.0116.
Copyright and License:
[Uncaptioned image] © Lianna Hambardzumyan, Konstantin Myasnikov, Artur Riazanov, Morgan Shirley,
and Adi Shraibman; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Circuit complexity
Related Version:
Full Version: https://arxiv.org/abs/2602.23503 [19]
Acknowledgements:
We thank Nathan Harms for helpful feedback that improved the presentation of the paper. We are grateful to Dmitry Sokolov for bringing ultra-lossless expanders to our attention, and to Amir Yehudayoff for pointing us to Warren’s theorem for the proof of Theorem 24. Finally, we thank Nati Linial, Shachar Lovett, Ryan Williams, and Amir Yehudayoff for many insightful discussions on spiky rank.
Editors:
Sayan Bhattacharya, Danupon Nanongkai, Michael Benedikt, and Gabriele Puppis

1 Introduction

Studying the structural parameters of matrices has proven to be a powerful way to capture the complexity of various computational models and to establish lower bounds. A recurring theme in complexity theory is identifying “well-behaved” matrix parameters: parameters strong enough to capture meaningful lower bounds, structured enough to admit combinatorial or algebraic analysis, and flexible enough to connect to multiple computational models. Rank, rigidity, γ2-norm, and blocky rank are all examples of such measures.

In this work we introduce and study spiky rank, a new matrix parameter that we argue meets the criteria for a well-behaved complexity measure. Spiky rank is a generalization of the blocky rank defined for communication complexity applications by Hambardzumyan, Hatami and Hatami [18] (Definitions 1 and 2) which later found a broad range of applications in complexity theory and beyond [7, 28, 36, 16, 13, 5].

Spiky rank adds a linear algebraic flavor to the combinatorial nature of blocky rank. This allows us to connect spiky rank to complexity measures that are both algebraic and combinatorial in their nature, such as matrix rigidity and circuits with “algebraic” gates.

Both blocky rank and spiky rank fit into a common paradigm for defining complexity measures: we begin by specifying a class of simple objects that we define to have complexity 1, and then define the complexity of a general object as the minimum number of such simple objects needed to construct it. To define blocky and spiky rank, we define our matrices of complexity 1 as follows:

Definition 1 (Blocky and Spiky Matrices).

A blocky matrix is a matrix that can be obtained from the identity matrix by permuting rows or columns, duplicating rows or columns, and adding all-zero rows or columns. Equivalently, a blocky matrix is a matrix whose rows and columns can be permuted so that it consists of diagonal blocks that are all-one submatrices (of possibly different sizes), with all off-diagonal blocks equal to zero. See Figure 1(a).

A spiky matrix is a blocky matrix in which each diagonal block is an arbitrary rank-one matrix (rather than an all-one matrix). Equivalently, a spiky matrix can be expressed as the entrywise product of a blocky matrix and a rank-one matrix. See Figure 1(c).

(a) A blocky matrix B.
(b) A rank-one matrix uv.
(c) A spiky matrix B(uv).
Figure 1: A spiky matrix is the entrywise product of a blocky matrix and a rank-one matrix.

We then define the complexity of general matrices as follows:

Definition 2 (Blocky and Spiky Rank).

A blocky matrix is said to have blocky rank one. The blocky rank of a matrix M, denoted 𝖻𝗋(M), is the minimum number r such that M can be written as a linear combination of r blocky rank one matrices.

A spiky matrix is said to have spiky rank one. The spiky rank of a matrix M, denoted 𝗌𝗉𝗋(M), is the minimum number r such that M can be written as a sum of r spiky rank one matrices.

1.1 Motivations for Blocky and Spiky Rank

Blocky rank.

The study of blocky matrices and related decompositions has a surprisingly long and diverse history. Although the terminology “blocky matrix” is recent, the same objects have been rediscovered many times under different names across theoretical computer science and mathematics. For instance, fat matchings [22, 29], equivalence graphs [12, 2], adjacency matrices of P4-free graphs [9, 5], and equality matrices [36]. Various versions of blocky decompositions have been studied; blocky rank, blocky covering (union of blocky matrices), blocky partitioning (sum of blocky matrices) and more. Each of these variants has turned out to capture useful phenomena in different research areas:

  • In communication complexity blocky decompositions have been used to study models of communication with access to the equality oracle [18, 28, 9, 16].

  • In operator theory, it was observed [18] that blocky matrices are the contractive idempotents of Schur multiplier algebras, making blocky rank a natural way to characterize the structure of idempotents in this setting (see also [13]).

  • In combinatorics, blocky coverings of different families of graphs have been studied in [12, 2, 9, 7].

  • In cryptography, Applebaum and Nir [5] studied blocky covering to characterize families of graphs that can be realized with linear secret sharing schemes of constant size.

  • In circuit complexity, blocky decompositions provide lower bounds: Jukna [22] used blocky partitioning to derive bounds against depth-3 circuits with parity gates, and Avraham and Yehudayoff [7] showed that blocky partitioning lower bounds depth-2 threshold circuits.

  • Branching program size lower bounds were obtained via blocky covering in [29].

Perhaps the most striking and far-reaching connection is due to Williams [36], who related blocky rank to the renowned Orthogonal Vectors Conjecture (OVC). In particular, Williams showed that to refute OVC – and consequently the Strong Exponential Time Hypothesis (SETH) – it would suffice to prove that the 2n×2n Disjointness matrix Dn satisfies 𝖻𝗋(Dn)2o(n). On the other hand, if one could show that 𝖻𝗋(Dn)2ω(n), this would immediately yield exponential lower bounds for depth-2 exact threshold circuits by [7, 36], resolving a notorious open problem in circuit complexity. This duality highlights the importance of studying blocky rank: either upper bounds or lower bounds on 𝖻𝗋(Dn) would have transformative implications in algorithms, fine-grained complexity, and circuit lower bounds.

Spiky rank.

As evident from this wide range of applications, blocky rank is a versatile measure. But it is fundamentally a combinatorial measure and, in particular, not a robust one. Small modifications to the matrix can cause its blocky rank to change drastically.

For example, in Figure 2, the identity matrix has 𝖻𝗋(In)=1, but if we place distinct weights on the diagonal, the blocky rank jumps to 𝖻𝗋(Dn)=n. Structurally these matrices are nearly the same – they are both just diagonal matrices – but blocky rank treats them in a completely different way. This lack of robustness makes blocky rank poorly suited for applications involving real-valued matrices or settings where algebraic properties matter. Spiky rank addresses this shortcoming by capturing the intuition that a diagonal matrix is expected to be simple regardless of its entries, in particular, 𝗌𝗉𝗋(Dn)=1.

In=[1000010000100001]Dn=[100002000030000n]
Figure 2: The identity matrix In (left) and the diagonal matrix Dn (right).

This extra flexibility makes spiky rank well-suited for problems that have both combinatorial and algebraic character. In this work, we demonstrate two such applications:

  • Spiky rank lower bounds matrix rigidity, a central problem in circuit and communication complexity (Section 1.2.1),

  • Spiky rank lower bounds the size of depth-2 ReLU circuits, which sit at the intersection of complexity theory and modern machine learning (Section 1.2.2).

Finally, note that spiky rank strictly strengthens blocky rank: every blocky decomposition is also spiky, so 𝗌𝗉𝗋(M)𝖻𝗋(M) for all M. Thus, lower bounds for spiky rank have the same implications as for blocky rank. On the other hand, upper bounds can be even more significant. As with blocky rank, upper bounds for spiky rank of the Disjointness matrix have far-reaching consequences: if the 2n×2n Disjointness matrix Disjn satisfies 𝗌𝗉𝗋(Disjn)2o(n), then Orthogonal Vector Conjecture (and thus SETH) would be false [37].

1.2 Results

Our main contribution is conceptual: we introduce spiky rank as a new matrix complexity measure, demonstrate its applicability to problems in complexity theory, and establish its basic properties. On the technical side, we prove lower bounds for both random and explicit matrix families and relate spiky rank to other classical parameters. Our main results are:

  1. 1.

    We demonstrate the natural applicability of spiky rank to two problems in complexity theory: matrix rigidity (Section 1.2.1) and ReLU circuits (Section 1.2.2).

  2. 2.

    We show that most Boolean matrices have spiky rank Ω(N/logN) (matching blocky rank), while most real matrices have spiky rank Ω(N) (Section 1.2.3).

  3. 3.

    We give the following bounds on the spiky rank of specific N×N matrices of interest:

    • Ω(logN) for the 1-Hamming Distance matrix (Section 6.2), simultaneously improving the previous loglogN bound for blocky rank as well [18].

    • Ω(logN) for adjacency matrices of certain expander graphs (Section 6.3).

    • Ω(logN/loglogN) for the Inner Product and Disjointness matrix (Section 6.4) (likely far from optimal but it represents a first step toward stronger bounds).

  4. 4.

    For Boolean matrices, we establish a dimension-free relation between blocky and spiky rank, implying in particular that constant spiky rank yields constant blocky rank (Section 1.2.5). In the full version, we further bound spiky rank in terms of matrix sparsity, separate it from the γ2-norm, and highlight gaps between the approximate and sign variants of blocky and spiky rank.

1.2.1 Application to Matrix Rigidity

Matrix rigidity, introduced by Valiant [32] and Grigoriev [17], asks how far a matrix is from having low rank: the rank-r rigidity M(r) is the minimum number of entry modifications required to reduce the rank of M to r.

Highly rigid matrices are central to long-standing open problems in complexity theory. Valiant showed that if MN×N satisfies M(ϵN)Ω(N1+δ) for some ϵ,δ>0, then any log-depth linear circuit computing the linear transformation xMx must have superlinear size. Additionally, Razborov [30] linked rigidity to communication complexity, showing that matrices in PHcc – the communication complexity analogue of the polynomial hierarchy – are upper-bounded in their rigidity. In particular, if a N×N Boolean matrix MPHcc, then its rigidity over 𝔽2 must satisfy

M𝔽2(2(loglogN)ω(1))N22(loglogN)O(1).

Hence, constructing sufficiently rigid matrices would yield lower bounds for PHcc and ultimately separate it from PSPACEcc – both of which remain major open problems in communication complexity.

Despite extensive study, progress has been limited on finding explicit rigid matrices: the best explicit constructions fall well short of the parameters needed for Valiant’s or Razborov’s applications.

A direct link between rigidity and spiky rank is rather straightforward to show: matrices with large spiky rank are also highly rigid.

Theorem 1.

Let M be a matrix and 0<r𝗌𝗉𝗋(M). Then,

M(r)(𝗌𝗉𝗋(M)r)24.

This has two key consequences. First, for Valiant’s target rank ϵN, an explicit matrix with 𝗌𝗉𝗋(M)=Ω(N) would achieve the strongest possible rigidity bound Ω(N2). Note that, since Boolean matrices have spiky rank at most N/logN, the only candidate matrices with spiky rank Ω(N) are the real matrices.

Second, for Razborov’s target rank 2(loglogN)ω(1), we would need to find an explicit matrix with spiky rank at least (see Corollary 33)

N2(loglogN)O(1).

This shows that spiky rank is not only a structural parameter, but also a potential tool for constructing explicit rigid matrices. With this goal in mind, we highlight the following open challenges to guide future research in this direction.

Major Open Problem 2.

Find an explicit real matrix MN×N with 𝗌𝗉𝗋(M)=Ω(N).

Major Open Problem 3.

Find an explicit Boolean matrix M{0,1}N×N with

𝗌𝗉𝗋(M)=N2(loglogN)O(1).

Although most Boolean matrices have spiky rank at least Ω(N/logN) and most real matrices have Ω(N) (see Section 1.2.3), constructing explicit examples that meet these bounds remains elusive. We expect that solving either problem above will be quite challenging. At present, the best lower bounds we can prove for explicit matrices are of order Ω(logN) (see Section 1.2.4).

Finally, we remark that Theorem 1 allows known non-rigidity results to translate into upper bounds on spiky rank. For instance, Alman and Williams [1] proved that for sufficiently small ε>0, the Inner Product matrix IPn of size 2n×2n satisfies

IPn(2n(1f(ε)))2n(1+ε),

where f(ε)=Θ(ε2/log(1/ε)). By Theorem 1, this yields an upper bound 𝗌𝗉𝗋(IPn)2(1δ)n for some constant δ>0, improving on the trivial bound of 2n/n (see Section 1.2.3).

1.2.2 Application to Circuit Complexity

Connections between blocky rank and circuit complexity were first observed by Jukna [22], and later by Avraham and Yehudayoff [7], who showed that the blocky rank of a matrix lower bounds the size of any ΣLTF circuit computing111In this paper, a circuit computes a matrix M if it takes as input a row label r and a column label c and outputs the value of the entry M[r,c]. it up to a certain loss factor (implicit in [20]). Here, ΣLTF denotes linear combinations of linear threshold functions (LTFs), where an LTF is a Boolean function of the form T(x)=𝟙[w,xα] for some wn and α.

Naturally this connection extends to spiky rank and to ΣReLU circuits, where ReLU functions generalize threshold functions by outputting max{0,w,xα}. The following proposition admits a simple proof, yet it is conceptually important.

Proposition 4.

For M:{0,1}n×{0,1}n, let s be the size of a ΣReLU circuit that computes M. Then

s𝗌𝗉𝗋(M)3(n+1).

Observe that ΣLTF circuits which compute Boolean functions can be simulated by ΣReLU circuits with only constant factor blowup in size. This follows from the fact that LTF gates for Boolean functions can have wn and α without loss of generality (see e.g. [26]) and the relation

𝟙[w,xα]=max{0,w,xα+1}max{0,w,xα}

which is true when wn and α.

Linear threshold and ReLU circuits are central models in both complexity theory and machine learning. ΣLTF circuits capture linear combinations of halfspaces and have been studied since the early days of threshold logic. On the other hand, ReLU gates have become the dominant activation function in modern neural networks, thus ΣReLU circuits are the basic building blocks of modern neural networks. However, our understanding of their expressive power from a Boolean complexity perspective remains limited. Understanding the power and limitations of these models is therefore both an important task in circuit complexity and in machine learning.

Despite decades of work, the best known lower bounds for these circuits are very weak: even for ΣLTF, essentially no superlinear lower bounds are known for explicit functions [31]. Williams [35] gave the strongest results to date, proving that for every k there exists a function in NP that requires circuit size nk in both the ΣLTF and ΣReLU models. For explicit functions, the strongest known result is due to Kane and Williams [23], who showed that even the stronger class LTFLTF requires size Ω(n3/2) to compute the n-variate Andreev function. Mukherjee and Basu [27] extend this to prove a lower bound of Ω(n1δ) for LTFReLU circuits computing the Andreev function.

Our theorem shows that spiky rank serves as a lower bound measure for depth-two neural networks with ReLU activation, just as blocky rank does for ΣLTF circuits. In particular, proving spiky rank lower bounds of ω((logN)5/2) for explicit N×N matrices would immediately yield new circuit lower bounds. At present, our lower bounds are of order logN; the first milestone toward this goal would be to establish any super-logarithmic spiky rank lower bounds for explicit matrices.

Open Problem 5.

Find an explicit real matrix MN×N with 𝗌𝗉𝗋(M)=ω((logN)5/2).

1.2.3 Bounds for Random Matrices

How large can the spiky rank of a matrix be? Trivially, both the spiky rank and blocky rank of an N×N Boolean matrix are at most N. For blocky rank, this bound can be improved to NlogN [29, 7], and therefore the same upper bound applies to spiky rank as well. A natural question is whether spiky rank admits a substantially stronger upper bound. Avraham and Yehudayoff [7] showed that most matrices are maximally hard for blocky rank. We prove that the same phenomenon holds for spiky rank.

Theorem 6 (Informal).

A random Boolean matrix of size N×N has spiky rank of order Ω(NlogN) with high probability.

For real matrices, the situation is more extreme: the trivial upper bound of N turns out to be essentially tight.

Theorem 7 (Informal).

A random real matrix of size N×N has spiky rank at least N2 with high probability.

This mirrors what is known for many other natural measures in complexity theory, such as circuit complexity and communication complexity, where random functions or matrices almost always exhibit the largest possible complexity.

1.2.4 Bounds for Specific Matrices

The lower bounds discussed above, while strong, are for random (non-explicit) matrices. For applications in complexity theory, however, it is essential to establish lower bounds for explicit families of matrices.

For blocky rank, the best known lower bounds for explicit families of matrices are of order Ω(logN) (Greater-Than [7], Inner-Product, Disjointness [18]).

For spiky rank, we develop a general lower bound framework (Theorem 37) and apply it to two examples: adjacency matrices of the Hamming cube and expander graphs. The common properties of these graphs that we use are:

Thinness:

every small subgraph has smaller edge-to-node ratio compared to the entire graph;

Induced matchings:

there are large induced matchings in all large subgraphs.

This framework can only yield lower bounds up to Θ(logN), the bottleneck is in the use of a variation of a lemma from [6, Lemma 12].

Hamming distance. Consider the matrix HD1n of size N×N with N=2n, whose rows and columns are indexed by all n-bit strings, and where

HD1n(x,y)={1if x and y differ in exactly 1 coordinate,0otherwise.

Applying our framework to this matrix yields the following bound:

Theorem 8.

𝗌𝗉𝗋(HD1n)Ω(logN).

This bound is not far from the optimal; upper bound on the spiky rank is 𝗌𝗉𝗋(HD1n)𝖻𝗋(HD1n)logN, where the last upper bound follows from the following blocky matrix decomposition: Let Mi be the matrix such that Mi(x,y)=1 iff xiyi and all the other entries of x and y coincide. Then Mi is a permutation matrix, so 𝖻𝗋(Mi)=1 and HD1n=M1++Mn.

This lower bound is quite remarkable considering that HD1n matrix is very simple with respect to many other matrix complexity measures, as we discuss in Section 1.2.5.

Open Problem 9.

Is 𝗌𝗉𝗋(HD1n)Ω(logN)?

Expander graphs. As a second application, let G be an (N,d,λ)-spectral expander: a d-regular graph with the second largest eigenvalue of the adjacency matrix not exceeding λ, denote its adjacency matrix with MG. We prove:

Theorem 10.

𝗌𝗉𝗋(MG)Ω(min(dλ,logNd2)).

This implies (see Corollary 40) an explicit Ω(logN) lower bound for the spiky rank of adjacency matrices of (N,Θ(log2N),Θ(logN))-spectral expanders.

Inner Product. The Boolean inner product matrix is defined by IPn(x,y)=i[n]xiyimod2 for x,y{0,1}n. For IPn we obtain bounds that are weak but improve over the trivial ones:

nlogn𝗌𝗉𝗋(IPn)2(1δ)n,

for some constant δ>0. The lower bound is established in Section 6.4, while the upper bound follows from non-rigidity results (Section 1.2.1).

Disjointness. The Disjointness matrix Disjn is the 2n×2n Boolean matrix whose rows and columns are indexed by {0,1}n, with Disjn(x,y)=1 iff x and y have no common 1’s. For the Disjointness matrix we observe the same lower bound as for IPn (Section 6.4), so we have

nlogn𝗌𝗉𝗋(Disjn)5n/5

where the upper bound follows from the 𝖻𝗋(Disjn) upper bound by Williams [36].

We believe the true lower bounds for Disjn and IPn are much closer to the upper bounds, and in particular ask:

Open Problem 11.

Is it true that 𝗌𝗉𝗋(IPn)=ω(n) or 𝗌𝗉𝗋(Disjn)=ω(n)?

Figure 3: The picture illustrates the upper and lower bounds for spiky rank: above the scale there are known lower bounds and upper bounds. Below the scale are consequences of the corresponding lower bounds. The dashed part of the scale represents the state of the art of explicit lower bounds. Here N is always the size of the matrix, so for HD1n, Disjn, and IPn, N=2n.

1.2.5 Relations to other matrix parameters

When introducing a new matrix parameter, a natural and useful question to understand is how it compares to existing ones: is it stronger, weaker, or qualitatively different? We defer the proofs of these relations to the full version.

Blocky rank

Since spiky rank can be seen as an algebraic strengthening of blocky rank, the first step is to ask whether this additional flexibility actually increases its power. Is blocky rank equivalent to spiky rank?

Question 12.

Does every Boolean matrix M satisfy 𝖻𝗋(M)poly(𝗌𝗉𝗋(M))?

We answer this question qualitatively: there is a dimension-free relation between spiky and blocky rank:

Theorem 13.

For any Boolean matrix M, 𝖻𝗋(M)𝗌𝗉𝗋(M)O(𝗌𝗉𝗋(M)).

This situation is reminiscent of tantalizing log-rank conjecture. There, too, a combinatorial parameter – the partition number – is vastly strengthened by adding algebraic property via matrix rank. Similar to Theorem 13, an exponential upper bound is known, and the challenge is whether this relation can be improved to a quasi-polynomial one.

𝜸𝟐-norm

Since our main goal is to establish lower bounds on spiky rank, it is natural to compare it with matrix parameters that are more amenable to analysis – particularly matrix norms. A well-studied example, with many applications in complexity theory, is the γ2-norm (see [8] and references within). For a matrix M, it is defined as

γ2(M)=minU,V:M=UVUrowVcol,

where Urow denotes the maximum 2-norm of the rows of U, and Vcol denotes the maximum 2-norm of the columns of V.

Strong lower bounds are known for γ2-norm for functions of interest; for instance, γ2(IPn)NΩ(1), where IPn(a,b)=i[n]aibimod2 (see e.g. [25]). This motivates comparing γ2 with rank-based parameters such as blocky rank and spiky rank. For blocky rank, it is known that for every Boolean matrix M, 𝖻𝗋(M)logγ2(M) [18, Proposition 3.1]. For spiky rank, while an analogous inequality currently seems out of reach, we observe a simple non-explicit separation:

Theorem 14.

There exists a Boolean matrix M such that 𝗌𝗉𝗋(M)Ω(γ22(M)).

The reverse direction – whether blocky rank is upper-bounded by a function of γ2-norm – is a conjecture asked in [18], and is equivalent to a conjecture in operator theory. Very recently, Goh and Hatami [13] proved that for N×N Boolean matrix M, 𝖻𝗋(M)2γ2(M)(logN)2, and the same bound also applies to spiky rank.

Sparsity

One of the most basic and easily computable parameters of a matrix is its sparsity, the number of non-zero entries. We show that sparse matrices necessarily have low spiky rank, a fact that becomes useful when relating spiky rank to matrix rigidity (Lemma 29). In particular, for any real matrix M with sparsity m, we prove

𝗌𝗉𝗋(M)m.

For Boolean matrices, this bound can be sharpened, and even expressed directly in terms of blocky rank. Specifically, for any Boolean matrix M with sparsity m, it holds

𝖻𝗋(M)O(mloglogmlogm).
Sign and approximate measures

Some computational models are naturally characterized not by exact matrix parameters but by their approximate or sign variants. For instance, randomized communication complexity is lower-bounded by approximate rank, while unbounded-error communication complexity is captured by sign-rank. Formally, for a matrix measure μ:{0,1}N×N, we define its approximate version by

μϵ(M)inf{μ(M)MMϵ , MN×N}.

Similarly, the sign version of μ is μ±(M)inf{μ(M)sign(M)=sign(M)}.

Approximate measures tend to be stronger than their exact counterparts, and often defy intuition developed for them. For example, the approximate version of the log-rank conjecture is false [10]. In case of spiky rank, it is sharply separated from all related approximate and sign measures of a matrix: 𝗌𝗉𝗋(HD1n)logN (Theorem 8), while rk±(HD1n)=O(1) [15], γ2,ε(HD1n)=O(1) [11], 𝖻𝗋ϵ(HD1n)=poly(loglogN).

This highlights an interesting contrast between the spiky and blocky rank. For blocky rank, there is a dimension-free relation with its sign version: we show that 𝖻𝗋(M)2𝖻𝗋±(M). By contrast, for spiky rank we have 𝗌𝗉𝗋(HD1n)Ω(logN) but 𝗌𝗉𝗋±(HD1n)rk±(HD1n)=O(1), ruling out any dimension-free relation between 𝗌𝗉𝗋 and 𝗌𝗉𝗋±.

Thus, proving lower bounds for sign spiky rank is an even harder challenge. At the same time, this opens a promising direction to seek explicit versions of the known non-explicit bounds below the logarithmic barrier of our framework. For example, we know that the adjacency matrix of a random degree-d graph has sign spiky rank Ω(d). A natural candidate for derandomization is a spectral expander, leading to the following question:

Question 15.

What is 𝗌𝗉𝗋±(MG) for a spectral expander G?

Organization.

Section 2 introduces preliminaries. In Section 3, we prove the existence of Boolean and real matrices with high spiky rank. The applications to matrix rigidity and circuit complexity are presented in Section 4 and Section 5, respectively. Section 6 establishes spiky rank bounds for specific explicit matrices. Additional material may be found in the full version.

2 Preliminaries

Our convention in this paper is to use capital letters for matrix dimensions, e.g. AM×N. We use the notation for the entrywise product. The notation [k] represents the set of integers {1,,k}.

For a field 𝔽, we denote 0𝔽 for the additive identity and 1𝔽 for the multiplicative identity. Addition, multiplication, and the entrywise product of matrices over 𝔽 are denoted +𝔽, 𝔽, and 𝔽 respectively.

For a matrix M𝔽N×N its sparsity 𝗌𝗉𝖺𝗋(M) is defined as the number of its non-0𝔽 elements.

We write 𝟙S for S[N]×[N] for the matrix in {0,1}N×N, where 𝟙S[i,j]=1 iff (i,j)S.

Blocky and spiky measures.

As some of our results generalize to arbitrary fields, we rewrite Definitions 1 and 2 here with respect to an arbitrary field 𝔽.

Definition 16.

A matrix M𝔽M×N is blocky if there is an integer k and two collections of disjoint sets S1,,Sk[M] and T1,,Tk[N] such that M[i,j]=1𝔽 if (i,j)[k]Si×Ti and M[i,j]=0𝔽 otherwise.

A matrix M𝔽M×N is spiky if there exists a blocky matrix B and a rank-one matrix C (where rank is taken over 𝔽) such that M=B𝔽C.

Definition 17.

The blocky rank of a matrix M𝔽M×N, denoted 𝖻𝗋𝔽(M), is the minimum r such that M=i[r]αiBi where each αi𝔽 and each Bi is blocky.

The spiky rank of a matrix M𝔽M×N, denoted 𝗌𝗉𝗋𝔽(M), is the minimum r such that M=i[r]Si where each Si is spiky.

When the underlying field is , we write simply 𝖻𝗋(M) and 𝗌𝗉𝗋(M).

Here are some easy facts about spiky rank:

Fact 18 (Subadditivity).

For any field 𝔽 and matrices A,B𝔽M×N, 𝗌𝗉𝗋𝔽(A+B)𝗌𝗉𝗋𝔽(A)+𝗌𝗉𝗋𝔽(B).

Fact 19 (Upper bounds).

For any field 𝔽 and matrix A𝔽M×N, 𝗌𝗉𝗋𝔽(A)rank𝔽(A) (where rank𝔽 is the rank over 𝔽) and 𝗌𝗉𝗋𝔽(A)𝖻𝗋𝔽(A).

Sign blocky rank and sign spiky rank.

We study the sign versions blocky rank and spiky rank over the reals. In the following definition, we assume that the matrices have no 0-entries. For a matrix A, the value sign(A) is the sign function applied element-wise:

sign(A)[i,j]={1if A[i,j]>01if A[i,j]<0.
Definition 20.

Let A be a matrix in M×N.

  • The sign blocky rank of A, denoted 𝖻𝗋±(A), is defined as

    𝖻𝗋±(A)=minBM×N{𝖻𝗋(B):sign(A)=sign(B)}.
  • The sign spiky rank of A, denoted 𝗌𝗉𝗋±(A), is defined as

    𝗌𝗉𝗋±(A)=minBM×N{𝗌𝗉𝗋(B):sign(A)=sign(B)}.
Fact 21 (Upper bounds for sign measures.).

For any matrix AM×N, 𝖻𝗋±(A)𝖻𝗋(A) and 𝗌𝗉𝗋±(A)𝗌𝗉𝗋(A).

3 Non-explicit lower bounds for spiky rank

3.1 Lower bound for random Boolean matrices

We first count the number of N×N Boolean matrices that have spiky rank r for a fixed number rN. It is sufficient to count the number of ±1-matrices of with spiky rank r. Indeed, if M{0,1}N×N, then its signed version M{±1}N×N defined as M=J2M with J being the N×N all-ones matrix, has 𝗌𝗉𝗋(M)r+1.

In fact, we will actually count the number of ±1-matrices with sign spiky rank r (Definition 20). As this measure is upper-bounded by spiky rank, this gives us an estimate for the number of matrices with spiky rank r as well.

The idea is to use Warren’s theorem [34] to count the number of possible sign patterns. The same idea and argument were used to count the number of matrices with sign rank r in [4]. Before stating the theorem, let us set it up first.

Let P1,P2,,Pm be real polynomials in variables each. For x, such that no polynomial among {Pi} vanishes at x, its sign pattern is the vector

(sign(P1(x)),sign(P2(x)),,sign(Pm(x))){1,1}m.
Theorem 22 ([34]).

Let P1,P2,,Pm be real polynomials, each in variables and of degree at most k. If m, then the number of different sign pattern for x ranging over all is at most (4ekm/).

Lemma 23.

The number of N×N sign matrices with sign spiky rank rN/2 is at most N6rN.

Proof.

Let M be a N×N sign matrix with sign spiky rank r. Let AA1++ArN×N be the sign representation of M: M=sign(A1++Ar). Let us represent each Ai=BiCi where Bi is a blocky matrix and Ci is a rank-1 matrix. There are at most N4rN collections B1,,Br, so let us fix one and count the number of potential matrices M that appear as we change C1,,Cr.

With fixed B1,,Br each entry of A can be written as a degree-2 polynomial in 2Nr. Let for each i[r], Ci=xiTyi for some vectors xi,yiN. Then for j,k[N] we have Ajk=i[r](Bi)jkxijyik. Then by Theorem 22 we get that the number of different sign matrices M of form sign(B1C1++BrCr) with fixed B1,,Br is at most (4e2N22rN)2rNN2rN, so in total there are at most N4rNn2rN matrices of sign spiky rank r.

Combining the lemma with the observation that 𝗌𝗉𝗋±(M)𝗌𝗉𝗋(M), we get that a uniformly random Boolean matrix has high spiky rank with high probability. More precisely,

Theorem 24.

If M is a N×N uniformly random Boolean matrix, then

Pr[𝗌𝗉𝗋(M)N12logN]Pr[𝗌𝗉𝗋±(M)N12logN]12N22.
Proof.

The number of matrices of sign spiky rank less than N/(12logN) is at most N6N2/(12logN)2N2/2 by Lemma 23 and the total number of N×N sign matrices is 2N2.

3.2 Lower bound for real matrices

The spiky rank of any N×N real matrix is trivially bounded above by N, since the real rank of a matrix is at most N, and every rank-one matrix is also a spiky matrix. For Boolean matrices, even blocky rank can be bounded above by O(N/logN). However, for general real matrices, we show that the trivial upper bound is essentially tight: there exist matrices whose spiky rank is Ω(N). The key point in the proof as in Lemma 23 is that all spiky matrices with a fixed blocky pattern can be seen as evaluations of a polynomial.

We will need the following standard probabilistic notions: a distribution μ over m is absolutely continuous if there exists an integrable density function f:m+ such that for every Lebesgue-measurable set Am we have μ(A)=Af𝑑λ where λ is the Lebesgue measure. Examples of absolutely continuous distributions include the uniform distribution over [0,1]m and the Gaussian distribution.

Theorem 25.

Let μ be any absolutely continuous probability distribution over N×N. Then

PrMμ[𝗌𝗉𝗋(M)N/2]=1.

In particular this holds for M[0,1]N×N.

Proof.

It suffices to show that the Lebesgue measure of the set of all matrices with spiky rank at most N/21 is zero. This immediately implies the existence of a matrix with spiky rank at least N/2.

Let MN×N be a matrix with 𝗌𝗉𝗋(M)=r. Then M=i[r]Si, where each Si is a spiky matrix. Each spiky matrix Si can be viewed as an entrywise product of a blocky matrix Bi and a rank-1 matrix aiTbi for some vectors ai,bin. Fix a collection of blocky matrices B1,,Br{0,1}N×N, and define the corresponding set of spiky rank-r matrices with the blocky pattern B1,,Br:

𝒮B1,,Br={i[r](aiTbi)Bia1,,ar,b1,,brN},

where is the Hadamard (entrywise) product. Then, the set of all matrices in N×N of spiky rank at most r is equal to the union over all such choices of B1,,Br: B1,Br𝒮B1,,Br. Since the number of blocky matrices is finite (trivially it does not exceed 2N2), to show that the union set has measure zero it is sufficient to show that each set 𝒮B1,,Br has measure zero whenever rN/21.

To show that 𝒮B1,,Br has Lebesgue measure zero, define a smooth function f:2rnn×n as follows: for a1,ar,b1,,brN and i,j[N], let

f(a1,b1,,ar,br)[i,j]=[r]1B[i,j]=1a[i]b[j].

The image of f is exactly 𝒮B1,,Br.

Fact 26.

Suppose f1,,fn:m are polynomials and Am is a set with Lebesgue measure zero. Then f(A){(f1(x),,fn(x))xA} also has Lebesgue measure zero.

Proof.

Since polynomials are C1-differentiable (in fact, infinitely differentiable) we apply the more general standard result (see e.g. [14, Proposition 1.3]) that states that any C1-differentiable transformation maps a measure-zero set to a measure-zero set.

Since 2rN has measure zero when 2rN<N2 and 𝒮B1,,Br=f(2rN), by Fact 26 𝒮B1,,Br has measure zero for all r<N/2, and so does their finite union.

4 Connection to matrix rigidity

Definition 27 (Rigidity).

For a matrix M𝔽N×N over a field 𝔽, the rank-r rigidity M𝔽(r) is the minimum number of entries that must be modified to reduce the rank of M to r. Formally,

M𝔽(r)=min{sM=A+𝔽C,rank𝔽(A)r,𝗌𝗉𝖺𝗋(C)s}.

We show that matrices with large spiky rank must also be highly rigid. Our main result is the following generalization of Theorem 1 to arbitrary fields.

Theorem 28.

Let M be a matrix over a field 𝔽 and 0<r𝗌𝗉𝗋𝔽(M). Then,

M𝔽(r)(𝗌𝗉𝗋𝔽(M)r)24.

The key ingredient is that sparse matrices have small spiky rank.

Lemma 29.

Let M𝔽N×N be a matrix over an arbitrary field 𝔽 with at most m non-zero entries. Then,

𝗌𝗉𝗋𝔽(M)2m.

We defer the proof to the full version. This immediately yields the following relationship between spiky rank and matrix rigidity.

Proof of Theorem 28.

Assume M𝔽(r)=s for some rN. Then M can be written as M=A+𝔽C, where A𝔽N×N has rank𝔽(A)=r and C𝔽N×N has at most s non-zero entries. By the subadditivity of spiky rank over 𝔽, we have

𝗌𝗉𝗋𝔽(M)=𝗌𝗉𝗋𝔽(A+𝔽C)𝗌𝗉𝗋𝔽(A)+𝗌𝗉𝗋𝔽(C)r+2s,

where the last inequality follows from Lemma 29 and the trivial observation that 𝗌𝗉𝗋𝔽(M)rank𝔽(M) for all M. Rearranging gives s(𝗌𝗉𝗋(M)r)2/4, which completes the proof.

 Remark 30.

The condition 0<r𝗌𝗉𝗋𝔽(M) in Theorem 28 limits its applicability for proving intermediate rigidity bounds. For instance, for Valiant’s target rank ϵN, the theorem applies only to matrices already known to have spiky rank Ω(N).

 Remark 31.

Lemma 29, thus also Theorem 28, is tight for real matrices: a random real matrix has spiky rank Ω(N) by Theorem 25, and thus sparsity Ω(N2) and rigidity Ω(N2), the maximum possible. For Boolean matrices, we can slightly improve the sparsity and spiky rank tradeoff (see the full version).

To discuss the implications of Theorem 28 for communication complexity, we use the following result of Razborov [30] and Wunderlich [38].

Theorem 32 ([38, Theorem 5.15]).

Let 𝔽 be a finite field and let M be a N×N matrix. If

M𝔽(2(loglogN)ω(1))Ω(N22(loglogN)O(1)).

Then MPHcc.

As an immediate consequence we obtain:

Corollary 33.

Suppose a N×N Boolean matrix M satisfies

𝗌𝗉𝗋𝔽(M)N2(loglogN)O(1),

where 𝔽 is the finite field. Then MPHcc.

Proof.

Applying Theorem 28 with r=𝗌𝗉𝗋𝔽(M)/2, we obtain

M𝔽(N2(loglogN)O(1))Ω(N22(loglogN)O(1)).

Since

M𝔽(2(loglogN)ω(1))M𝔽(N2(loglogN)O(1)),

the claim follows directly from Theorem 32.

5 Connection to circuit complexity

Definition 34.

A ReLU gate (rectified linear unit) is a function :{0,1}n+ such that there exist a vector wn and a scalar a such that for all x{0,1}n

(x)=max{0,w,x+a}.

Intuitively, a ReLU function is a weighted linear threshold function. For our purposes it will be convenient to view a ReLU gate as a matrix: a ReLU gate is a matrix L:{0,1}n×{0,1}n+ such that there exist vectors w1,w2n and a scalar α for which

L(x,y)=max{0,w1,x+w2,y+α}.

We now show that ΣReLU circuit complexity – the number of ReLU gates – of any function is upper bounded by its spiky rank up to an O(n) factor. As a first step, we bound the spiky rank of a single ReLU gate.

Claim 35.

Let L:{0,1}n×{0,1}n+ be a ReLU gate. Then

𝗌𝗉𝗋(L)3n+3.
Proof.

Let w1,w2n and α be such that L(x,y)=max{0,w1,x+w2,y+α}. Denote the underlying function L(x,y)=w1,x+w2,y+α. Observe that L is a sum of three rank-1 matrices w1,x, w2,y and α (each depending only on x, y or being constant), so rank(L)3.

Define TL as the threshold function corresponding to L as follows:

TL(x,y)={1, if L(x,y)>00, otherwise.

Avraham and Yehudayoff [7] show that any threshold function on {0,1}n×{0,1}n can be written as a disjoint union of at most n+1 blocky matrices:

TL=B1++Bn+1,

where B1,,Bn+1 are disjoint blocky matrices, i.e., their supports are disjoint rectangles in the matrix. We now define weighted matrices R1,,Rn+1 corresponding to Bi’s by placing the original values of L on the support of each Bi:

Ri(x,y)={L(x,y), if Bi(x,y)=10, otherwise.

Now since Bi’s are disjoint, we have L=R1++Rn+1. The only problem now is that Ri might not be a spiky matrix as it might have blocks of rank greater than 1. However, the blocks of Ri are submatrices of L, so each block of Ri has rank at most 3. If there is a block Ri of rank greater than 1, then we can decompose Ri into at most three matrices that have rank-1 on that block. If multiple blocks in the same Ri have rank greater than 1, then we can still partition Ri into at most 3 matrices that have rank-1 on all blocks, because the blocks are on disjoint rows and columns. This shows that each Ri can be decomposed into at most three spiky matrices. Hence, 𝗌𝗉𝗋(M)3(n+1).

Proposition 36.

For M:{0,1}n×{0,1}n, let s be the size of a ΣReLU circuit that computes M. Then

𝗌𝗉𝗋(M)3(n+1)s.
Proof.

For the lower bound, let M=i=1sLi, where Li are ReLU functions. Then,

𝗌𝗉𝗋(M)i=1s𝗌𝗉𝗋(Li)s(3n+3).

6 Explicit lower bounds for spiky rank

6.1 Lower bound framework

In this section we give a framework for proving spiky ranks lower bounds that applies to thin matrices such as HD1n and adjacency matrices of expander graphs.

Theorem 37.

Let MN×N, s,k,D be positive integers and γ(0,1). If

  1. (P1)

    For every S[N] and T[N] with |S|s and |T|s we have 𝗌𝗉𝖺𝗋(M|S×T)k(|S|+|T|).

  2. (P2)

    For every MM (entrywise) such that 𝗌𝗉𝖺𝗋(MM)γ𝗌𝗉𝖺𝗋(M), there exist sets S[N] and T[N] such that |S|=|T|=D and M|S×T is a permutation matrix.

then

𝗌𝗉𝗋(M)Ω(min(γ𝗌𝗉𝖺𝗋(M)kN,log(Ds2N))).

We need the following lemma for the proof. The proof of it is identical to the proof of [7, Lemma 12], but we include it here for completeness.

Lemma 38.

Let IN be identity matrix of size N×N. If there is a blocky (or spiky) matrix B, such that BIN=0 (i.e. there is no entry where both IN and B are nonzero). Then there exists V[N], s.t. |V|N/4 and B|V×V=0.

Proof.

Let blocks of B be A1×B1,A2×B2,,Ak×Bk. Define two random subsets of [N] as follows. Let ξ1,ξ2,,ξk be i.i.d. uniformly distributed on {0,1}.

Let S be the complement of a:ξa=0Sa, and T be the complement of a:ξa=1Ta. Let I={i[N](i,i)S×T}. Note that B|S×T is zero matrix. For each i[N], the probability that (i,i)S×T is 1/4, because Bii=0. Thus there is a choice for S×T so that |I|N/4. For such S and T, we define V=ST.

Proof of Theorem 37.
Figure 4: The picture illustrates the proof for the matrix HD15. The parts Ti are colored red, Hi are colored green. By 1 not too many 1-cells (black ones) are covered by the small blocks, so by 2 we find a permutation submatrix S×T with very few covered 1-cells, and then further shrink it using Lemma 38 to A×B that does not contain any red cells at all. Then the rank of HD15|A×B is maximal, but the total rank of large blocks is small, which is a contradiction.

Let M=i[t]Si. We call a block small, if the sum of its dimensions is at most s, and big otherwise. Assume the contrary and let tγ𝗌𝗉𝖺𝗋(M)2kN and t14log2(Ds2N).

We represent Si=Ti+Hi where Ti is comprised of the small blocks of the matrix Ti and Hi of the big ones (i.e. the sets of non-zero columns of Ti and Hi are disjoint, same for the non-zero rows).

By 1, each Ti covers (blocky matrix covers matrix entry, if this entry is non-zero in this matrix) at most 2kN 1-entries of M, because the sum of block sizes in a spiky matrix does not exceed 2N. This implies that all small matrices cover at most 2Ntkγ𝗌𝗉𝖺𝗋(M) ones.

By 2 we get a permutation submatrix S×T and applying Lemma 38 to it t times, we obtain a permutation submatrix A×B of M with |A|=|B|=D/22t, s.t. for each i the intersection of Ti and A×B is empty. Thus, we obtain

rk(Mi[t]Ti)D22t (1)

On the other hand, each Hi contains at most 2N/s blocks, that means

rk(i[t]Hi)i[t]rk(Hi)t2Ns. (2)

Combining (1) and (2), we get D22t2Nt/s, which implies Ds2N23t, which contradicts our assumption.

6.2 Hamming-Distance-1

HD1n is the function which determines whether the Hamming distance between the two input strings is 1. Its matrix is the adjacency matrix of the n-dimensional Hamming cube.

HD1n(x,y)={1if xy1=10otherwise

See 8 We will make use of the following lemma.

Lemma 39 (Lemma 3.10 in [24]).

If G=(V,E) is a subgraph of a Hamming cube, then

|V|2|E|/|V|.
Proof of Theorem 8.

We use Theorem 37 with N=2n, M=HD1n, s=2n, k=2N, γ=0.5, D=N/4. We now prove two conditions of the theorem. 1 directly from Lemma 39 applied to ST.

We proceed to prove 2. The 1-entries of HD1n are partitioned into the following 2n subsets:

K2ij={(x,y)HD1n(x,y)=1,xi=j,yi=1j} for i[n]; j{0,1}.

Thus, each 1-entry in HD1n belongs to exactly one of these subsets, and each subset has size 2n1. The subsets have the following property: if (x,y) and (v,w) are distinct elements of Kt, then xw1>1 and vy1>1, because we have xi=viyi=wi for some i. Thus, we have at least half of the 1-entries present in M, which means one of the groups has at least 122n1 elements in the intersection with M. Without loss of generality, we assume that it is the subset K1. Then we set

S{1}×{0,1}n1,T{0}×{0,1}n1,

such that S×T covers exactly 2n2 1-entries of K1, which are present in M, and they form a permutation submatrix. So, the desired Ω(n) bound is immediately implied by the theorem.

6.3 Expanders

For a graph G, we denote by MG{0,1}N×N its adjacency matrix, i.e.,

(MG)u,v={1if (u,v)E,0otherwise.

We recall that an (N,d,λ)-spectral expander is a d-regular graph G=(V,E) on N vertices such that the second largest eigenvalue of its adjacency matrix MG is at most λ in absolute value (see, for example, [21]).

Theorem 10. [Restated, see original statement.]

𝗌𝗉𝗋(MG)Ω(min(dλ,logNd2)).

Proof.

We use Theorem 37 with M=MG, s=Nλ10d, k=λ, D=N4d, γ=1/2. We check the theorem conditions. 1 holds by the expander mixing lemma

e(S,T)d|S||T|N+λ|S||T|dsN|S||T|+λ|S||T|2λ|S||T|λ(|S|+|T|).

Let us prove 2. We pick edges greedily. At each step, when edge e={u,v} is picked, we throw away all vertices that are connected either to u or v. At each step, at most 2d1 vertices are used, so there will be at least γN2dN4d steps and we get at least D×D permutation submatrix.

Applying the theorem, we get the desired bound. As an implication, we get an Ω(logN) lower bound for (N,Θ(log2N),Θ(logN))-expanders:

Corollary 40.

There exists an explicit graph G with N nodes and degree Θ(log2N) such that

𝗌𝗉𝗋(MG)=Ω(logN).
Proof.

It is shown in [3, Section 2.3] that for every prime p with pmod4=1 there exist two primes q1,q2=pΘ(1) such that for every 1 there exists a degree p+1 Ramanujan graph G with Θ((q1q2)3) nodes, i.e. G is a (Θ((q1q2)3),p+1,2p)-spectral expander. Then the number of nodes in G is N=pΘ(), so taking =p/logp we get that N=2Θ(p), so p=Θ(log2N). Finally, we apply Theorem 10 to get the desired result.

We get a stronger corollary for ultra-lossless vertex expanders. These are known to exist, but there are currently no explicit constructions.

Theorem 41.

Let G be a d-regular graph with N vertices and vertex expansion dO(1) for at most Nδ vertex subsets, where δ(0,1) (i.e. for any SV(G), it has at least (dO(1))|S| neighbors). Let MG be its adjacency matrix. Then 𝗌𝗉𝗋(MG)Ω(min(d,logN)).

Proof.

The proof is similar to the proof of Theorem 10, except for the proof of 1. It is proved in the following way.

Let SV with |S|Nδ. By the assumption of the claim, it has at least (dO(1))|S| neighbors. On the other hand the sum of degrees of vertices in |S| is d|S|. Thus, for any TV(G) with |T|=|S| we have e(S,T)d|S|(dO(1))|S|+|T|=O(|S|).

6.4 Inner Product and Disjointness

A simple consequence of the lower bound for random Boolean matrices (Theorem 24) is the following. Suppose an explicit N×N matrix M contains a k×2k submatrix with all columns distinct. Then 𝗌𝗉𝗋±(M)Ω(klogk). Indeed, such a submatrix necessarily contains every possible k×k Boolean matrix (up to column permutations, which do not affect 𝗌𝗉𝗋±), and hence in particular one with maximum sign spiky rank Ω(klogk).

The largest k for which this holds is exactly the Vapnik–Chervonenkis dimension VC(M) of the matrix M [33]. Thus we obtain the general bound

𝗌𝗉𝗋±(M)Ω(VC(M)logVC(M)).

This observation immediately implies sign spiky rank lower bounds for Inner Product and Disjointness matrices. Recall that IPn(x,y)=i[n]xiyimod2, where x,y{0,1}n, and Disjn{0,1}N×N is defined with Disjn(x,y)i[n]¬xi¬yi where x and y are identified with their binary representations.

Proposition 42.

𝗌𝗉𝗋±(IPn)Ω(nlogn) and 𝗌𝗉𝗋±(Disjn)Ω(nlogn).

Proof.

The VC-dimension of IPn and Disjn is at least n. Consider the n×2n submatrix of Disjn (or IPn) consisting of the n rows indexed by the standard basis vectors

ei(0i1,1,0ni),i[n].

For a column indexed by x{0,1}n, its entries in this submatrix are as follows:

  • In the inner product matrix IPn, the column is exactly the vector x.

  • In the Disjointness matrix Disjn, the column is 1nx, where 1n is the all-ones vector.

Thus the submatrix captures all distinct column patterns, corresponding either to x (for IPn) or its complement (for Disjn).

References

  • [1] Josh Alman and Ryan Williams. Probabilistic rank and matrix rigidity. In Proceedings of the 49th Annual ACM SIGACT Symposium on Theory of Computing, STOC 2017, pages 641–652, New York, NY, USA, 2017. Association for Computing Machinery. doi:10.1145/3055399.3055484.
  • [2] Noga Alon. Covering graphs by the minimum number of equivalence relations. Combinatorica, 6(3):201–206, 1986. doi:10.1007/BF02579381.
  • [3] Noga Alon. Explicit expanders of every degree and size. Combinatorica, 41(4):447–463, 2021. doi:10.1007/s00493-020-4429-x.
  • [4] Noga Alon, Shay Moran, and Amir Yehudayoff. Sign rank versus VC dimension. In Conference on Learning Theory, pages 47–80. PMLR, 2016. URL: http://proceedings.mlr.press/v49/alon16.html.
  • [5] Benny Applebaum and Oded Nir. The meta-complexity of secret sharing. In Proceedings of the 57th Annual ACM Symposium on Theory of Computing, pages 965–976, 2025. doi:10.1145/3717823.3718161.
  • [6] Daniel Avraham and Amir Yehudayoff. On blocky ranks of matrices. Technical Report TR22-137, Electronic Colloquium on Computational Complexity (ECCC), September 2022. URL: https://eccc.weizmann.ac.il/report/2022/137/.
  • [7] Daniel Avraham and Amir Yehudayoff. On blocky ranks of matrices. computational complexity, 33, March 2024. doi:10.1007/s00037-024-00248-1.
  • [8] Igor Balla, Lianna Hambardzumyan, and István Tomon. Factorization norms and an inverse theorem for MaxCut. Mathematische Annalen, 394(3):52, February 2026. doi:10.1007/s00208-026-03355-2.
  • [9] Alexander R. Block, Simina Brânzei, Hemanta K. Maji, Himanshi Mehta, Tamalika Mukherjee, and Hai H. Nguyen. P4-free partition and cover numbers & applications. In Stefano Tessaro, editor, 2nd Conference on Information-Theoretic Cryptography (ITC 2021), volume 199 of Leibniz International Proceedings in Informatics (LIPIcs), pages 16:1–16:25, Dagstuhl, Germany, 2021. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.ITC.2021.16.
  • [10] Arkadev Chattopadhyay, Nikhil S. Mande, and Suhail Sherif. The log-approximate-rank conjecture is false. Journal of the ACM, 67, 2020. doi:10.1145/3396695.
  • [11] Tsun Ming Cheung, Hamed Hatami, Rosie Zhao, and Itai Zilberstein. Boolean functions with small approximate spectral norm. Discrete Analysis, September 18 2024. doi:10.19086/da.122971.
  • [12] Péter Frankl. Covering graphs by equivalence relations. In North-Holland Mathematics Studies, volume 60, pages 125–127. Elsevier, 1982.
  • [13] Marcel K Goh and Hamed Hatami. Block complexity and idempotent Schur multipliers. International Mathematics Research Notices, 2025(24):356, 2025.
  • [14] Martin Golubitsky and Victor Guillemin. Stable mappings and their singularities, volume 14. Springer Science & Business Media, 2012.
  • [15] Mika Göös, Nathaniel Harms, Valentin Imbach, and Dmitry Sokolov. Sign-rank of k-Hamming Distance is constant. arXiv preprint arXiv:2506.12022, 2025. doi:10.48550/arXiv.2506.12022.
  • [16] Mika Göös, Nathaniel Harms, and Artur Riazanov. Equality is far weaker than constant-cost communication. arXiv preprint arXiv:2507.11162, 2025. doi:10.48550/arXiv.2507.11162.
  • [17] Dima Grigoriev. An application of separability and independence notions for proving lower bounds of circuit complexity. J.Soviet Math., 14:1450–1456, January 1980.
  • [18] Lianna Hambardzumyan, Hamed Hatami, and Pooya Hatami. Dimension-free bounds and structural results in communication complexity. Israel Journal of Mathematics, 253(2):555–616, 2023.
  • [19] Lianna Hambardzumyan, Konstantin Myasnikov, Artur Riazanov, Morgan Shirley, and Adi Shraibman. Spiky rank and its applications to rigidity and circuits, 2026. doi:10.48550/arXiv.2602.23503.
  • [20] Kristoffer Arnsfelt Hansen and Vladimir V Podolskii. Exact threshold circuits. In 2010 IEEE 25th Annual Conference on Computational Complexity, pages 270–279. IEEE, 2010. doi:10.1109/CCC.2010.33.
  • [21] Shlomo Hoory, Nathan Linial, and Avi Wigderson. Expander graphs and their applications. Bull. Amer. Math. Soc., 43(04):439–562, August 2006. doi:10.1090/s0273-0979-06-01126-8.
  • [22] Stasys Jukna. On graph complexity. Combinatorics, Probability and Computing, 15(6):855–876, 2006. doi:10.1017/S0963548306007620.
  • [23] Daniel M Kane and Ryan Williams. Super-linear gate and super-quadratic wire lower bounds for depth-two and depth-three threshold circuits. In Proceedings of the forty-eighth annual ACM symposium on Theory of Computing, pages 633–643, 2016. doi:10.1145/2897518.2897636.
  • [24] Rohan Karthikeyan, Siddharth Sinha, and Vallabh Patil. On the resolution of the sensitivity conjecture. Bulletin of the American Mathematical Society, 57:1, March 2020. doi:10.1090/bull/1697.
  • [25] Troy Lee and Adi Shraibman. Lower bounds in communication complexity. Foundations and Trends® in Theoretical Computer Science, 3(4):263–399, 2009. doi:10.1561/0400000040.
  • [26] Meena Mahajan. Depth-2 threshold circuits. Resonance, 24(3):371–380, March 2019. doi:10.1007/s12045-019-0786-4.
  • [27] Anirbit Mukherjee and Amitabh Basu. Lower bounds over Boolean inputs for deep neural networks with relu gates. arXiv preprint arXiv:1711.03073, 2017. arXiv:1711.03073.
  • [28] Toniann Pitassi, Morgan Shirley, and Adi Shraibman. The strength of equality oracles in communication. In Yael Tauman Kalai, editor, 14th Innovations in Theoretical Computer Science Conference (ITCS 2023), volume 251 of Leibniz International Proceedings in Informatics (LIPIcs), pages 89:1–89:19, Dagstuhl, Germany, 2023. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.ITCS.2023.89.
  • [29] P. Pudlák and V. Rödl. Some combinatorial-algebraic problems from complexity theory. Discrete Mathematics, 136(1):253–279, 1994. doi:10.1016/0012-365X(94)00115-Y.
  • [30] Alexander Razborov. On rigid matrices (in russian), 1989. URL: https://eccc.weizmann.ac.il/report/2022/137/.
  • [31] Vwani P Roychowdhury, Alon Orlitsky, and Kai-Yeung Siu. Lower bounds on threshold and related circuits via communication complexity. IEEE Transactions on Information Theory, 40(2):467–474, 1994. doi:10.1109/18.312169.
  • [32] Leslie G Valiant. Graph-theoretic arguments in low-level complexity. In Mathematical Foundations of Computer Science 1977: Proceedings, 6th Symposium, Tatranská Lomnica September 5–9, 1977 6, pages 162–176. Springer, 1977. doi:10.1007/3-540-08353-7_135.
  • [33] V. N. Vapnik and A. Ya. Chervonenkis. On the uniform convergence of relative frequencies of events to their probabilities. Theory of Probability & Its Applications, 16(2):264–280, 1971.
  • [34] Hugh E Warren. Lower bounds for approximation by nonlinear manifolds. Transactions of the American Mathematical Society, 133(1):167–178, 1968.
  • [35] Ryan Williams. Limits on representing Boolean functions by linear combinations of simple functions: Thresholds, ReLUs, and low-degree polynomials. arXiv preprint arXiv:1802.09121, 2018. arXiv:1802.09121.
  • [36] Ryan Williams. The Orthogonal Vectors Conjecture and non-uniform circuit lower bounds. In 2024 IEEE 65th Annual Symposium on Foundations of Computer Science (FOCS), pages 1372–1387, October 2024. doi:10.1109/FOCS61266.2024.00088.
  • [37] Ryan Williams. Personal communication, 2024.
  • [38] Henning Wunderlich. On a theorem of Razborov. computational complexity, 21:431–477, 2012. doi:10.1007/S00037-011-0021-5.