Abstract 1 Introduction 2 Notations and Preliminaries 3 𝗿𝗼𝗔𝗕𝗣 Factor Non-Closure 4 𝗿𝗼𝗔𝗕𝗣 Complexity of Symmetric Polynomials 5 Non-closure corollaries for 𝗿𝗼𝗔𝗕𝗣 6 Conclusions and Open Problems References

On Closure Properties of Read-Once Oblivious Algebraic Branching Programs

Robert Andrews ORCID Cheriton School of Computer Science, University of Waterloo, Canada Jules Armand ORCID Université Savoie Mont Blanc, CNRS, LAMA, France Prateek Dwivedi ORCID IT University of Copenhagen, Denmark Magnus Rahbek Dalgaard
Hansen
ORCID IT University of Copenhagen, Denmark
Nutan Limaye ORCID IT University of Copenhagen, Denmark Srikanth Srinivasan ORCID Department of Computer Science, University of Copenhagen, Denmark Sébastien Tavenas ORCID Université Savoie Mont Blanc, CNRS, LAMA, France
Abstract

We investigate the closure properties of read-once oblivious Algebraic Branching Programs (roABPs) under various natural algebraic operations and prove the following.

  • Non-closure under factoring: There is a sequence of explicit polynomials (fn(x1,,xn))n that have poly(n)-sized roABPs such that some irreducible factor of fn requires roABPs of superpolynomial size in any order.

  • Non-closure under powering: There is a sequence of polynomials (fn(x1,,xn))n with poly(n)-sized roABPs such that any super-constant power of fn does not have roABPs of polynomial size in any order (and fnn requires exponential size in any order).

  • Non-closure under symmetric operations: There are symmetric polynomials (fn(e1,,en))n that have roABPs of polynomial size such that fn(x1,,xn) do not have roABPs of subexponential size. (Here, e1,,en denote the elementary symmetric polynomials in n variables.)

These results should be viewed in light of known results on models such as algebraic circuits, (general) algebraic branching programs, formulas and constant-depth circuits, all of which are known to be closed under these operations.

To prove non-closure under factoring, we construct hard polynomials based on expander graphs using gadgets that lift their hardness from sparse polynomials to roABPs. For symmetric compositions, we show that the circulant polynomial requires roABPs of exponential size in every variable order.

Keywords and phrases:
Factoring, Closure Properties, Sparsity Bounds, Symmetric Polynomials, roABP, Expander Graphs
Funding:
Srikanth Srinivasan: Supported by European Research Council (ERC) under grant agreement no. 101125652 (ALBA).
Sébastien Tavenas: Supported by ANR project VONBICA (ANR-22-CE48-0007).
Copyright and License:
[Uncaptioned image] © Robert Andrews, Jules Armand, Prateek Dwivedi, Magnus Rahbek Dalgaard Hansen, Nutan Limaye, Srikanth Srinivasan, and Sébastien Tavenas; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Complexity classes
; Theory of computation Problems, reductions and completeness ; Theory of computation Circuit complexity ; Theory of computation Algebraic complexity theory
Related Version:
Full Version: https://arxiv.org/abs/2509.10725
Acknowledgements:
This research in this paper was carried out during a visit of some of the authors to Université Savoie Mont Blanc. The authors thank the Laboratoire de Mathématiques (LAMA) at Université Savoie Mont Blanc, which is supported by the AAP GAFA project, for their hospitality and the conducive environment for research. NL would like to thank Université Savoie Mont Blanc for hosting her as a visiting researcher.
Funding:
PD, MH, and NL thank Independent Research Fund Denmark (grant agreement No. 10.46540/3103-00116B) and the support of Basic Algorithms Research Copenhagen (BARC), funded by VILLUM Foundation Grant 54451.
Editor:
Shubhangi Saraf

1 Introduction

Given any computational model, it is natural to study the closure properties of the model with respect to simple operations. In Boolean complexity, these simple operations typically take the form of Boolean operations such as union, intersection, complement etc. In the setting of algebraic complexity, the object of computation is a multivariate polynomial f𝔽[x1,,xn]. Here, it is intuitive to consider closure properties under algebraic operations.

In this paper, we study the closure properties of a very well studied model of algebraic computation, namely read-once oblivious Algebraic Branching programs (roABPs, see definition 6). The interest in this model stems from the fact that it is both expressive enough to capture many natural algorithmic paradigms while at the same time possible to analyze using standard “complexity measures” [44].

In particular, 𝗋𝗈𝖠𝖡𝖯s can efficiently compute several polynomials of interest, including elementary symmetric polynomials and iterated matrix multiplication111This does not imply any completeness result as roABP is not closed under projection., the latter being provably hard to compute for constant-depth circuits [40, 6, 43]. In addition, 𝗋𝗈𝖠𝖡𝖯s subsume well-studied models such as sparse polynomials, set-multilinear and diagonal depth-3 circuits [38], as well as polynomials with low partial derivative dimension [9]. On the other hand, this is also one of the few models where we have a perfect characterization of the complexity of any given polynomial (in the form of the rank of an associated matrix) and where we also have a perfect understanding of border complexity [25]. As a consequence, this model has played a central role in research on lower bounds, polynomial identity testing algorithms and “debordering” results [21, 22, 14].

We study the closure properties of this model under basic algebraic operations such as factorization, powering, and inversion under composition with an important algebraic map (the elementary symmetric polynomial map). Apart from being natural questions about any computational model, such investigations have played a vital role in understanding hardness-randomness tradeoffs [35, 23, 17, 10] and the complexity of basic algebraic problems such as the Resultant and GCD [5, 11] in other algebraic models.

1.1 Main Results

In contrast with what is known for other models, our results are mostly negative. Specifically, we show the following.

𝗿𝗼𝗔𝗕𝗣 factor non-closure

Our first main result shows that there are explicit polynomial sequences that have small roABPs but with an irreducible factor that has 𝗋𝗈𝖠𝖡𝖯 complexity super-polynomial in n. Specifically, we prove 𝗋𝗈𝖠𝖡𝖯 complexity lower bound for a root, which is an irreducible factor of the form xnf(x1,,xn1), even when the 𝗋𝗈𝖠𝖡𝖯 is allowed to scan the variables in any order. The formal statement is as follows.

Theorem 1 (𝗋𝗈𝖠𝖡𝖯 factor non-closure).

The following holds over any field. Let n be a parameter and dn. There exists an n-variate polynomial f of degree d computable by an 𝗋𝗈𝖠𝖡𝖯 of width w2O(n), such that one of its (irreducible) factors g requires an 𝗋𝗈𝖠𝖡𝖯 of width wΩ(logd) in every variable order.

Note that, an 𝗋𝗈𝖠𝖡𝖯 computing an n-variate polynomial by definition has only n layers. Hence, the size and the width of an 𝗋𝗈𝖠𝖡𝖯 are polynomially related. Secondly, the size and width parameters in the theorem above are not polynomial in the number of variables, but they can be easily made polynomial by padding with some additional “dummy” variables. In particular, one should think of n above as logarithmic in the number of “actual” variables and d as a growing parameter, up to a polynomial in the number of variables.

This is in contrast to other algebraic models such as algebraic circuits [36, 37], branching programs [52], formulas and constant-depth circuits [10], all of which satisfy the property that factors of a polynomial f have complexity comparable to that of f. An exception to this rule is the family of sparse polynomials [57], and our construction is based on “lifting” this example to the setting of 𝗋𝗈𝖠𝖡𝖯.

𝗿𝗼𝗔𝗕𝗣 complexity of Symmetric Composition

We study an analogue of the result of Bläser and Jindal [16] for 𝗋𝗈𝖠𝖡𝖯. More specifically, a classical result in the theory of symmetric functions says that any symmetric polynomial222A polynomial is symmetric if it is invariant under any permutation of its variables. fsym(x1,,xn) can be written as a unique polynomial combination f of the elementary symmetric polynomials 𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn, where 𝖤𝖲𝗒𝗆nd is the n-variate elementary symmetric polynomial of degree d. Looking for a computational analogue of this theorem, Lipton and Regan [41], asked: what is the complexity of fsym vis-à-vis that of f?

Bläser and Jindal [16] showed that the complexity of f and fsym are polynomially related in the algebraic circuit model. Recently, the work of Bhattacharjee, Kumar, Rai, Ramanathan, Saptharishi and Saraf [11] extended this result to formulas and constant-depth circuits to show that fundamental computations such as GCD, resultants and discriminants have efficient constant-depth circuits in any characteristic. This generalizes a similar result of Andrews and Wigderson [5] in characteristic 0.

We show in this paper that the 𝗋𝗈𝖠𝖡𝖯 complexity of a polynomial f and its symmetric counterpart f𝗌𝗒𝗆 can differ significantly. Taking fsym=d=0n𝖤𝖲𝗒𝗆nd(x1k,,xnk), we can show that fsym is easy but f is exponentially hard.

Theorem 2.

The following holds over fields of characteristic zero. Let n be a parameter. There exists an n-variate polynomial f such that the symmetric polynomial f𝗌𝗒𝗆f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn) is computable by an 𝗋𝗈𝖠𝖡𝖯 of constant width in every variable order, but any 𝗋𝗈𝖠𝖡𝖯 computing f in any variable order must have width 2Ω(n).

In the other direction, our next result shows that even if a polynomial f is easy to compute by 𝗋𝗈𝖠𝖡𝖯, its symmetric counterpart f𝗌𝗒𝗆 can still be hard for 𝗋𝗈𝖠𝖡𝖯 – once again in sharp contrast to the known results for circuits, formulas, and constant-depth circuits. Specifically, the lower bound for a power of the elementary symmetric polynomial yields an example where f is easy but fsym is exponentially hard.333This is an especially strong contrast to the other models where it is trivial to show that if f is easy, then so is fsym.

Theorem 3.

The following holds over fields of characteristic zero. Let n be a parameter. There exists an n-variate polynomial f computable by an 𝗋𝗈𝖠𝖡𝖯 of constant width such that its respective symmetric polynomial f𝗌𝗒𝗆=f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn) requires an 𝗋𝗈𝖠𝖡𝖯 of width 2Ω(n) in every variable order.

𝗿𝗼𝗔𝗕𝗣 non-closure corollaries

We also investigate the power of roABPs in relation to powering an efficiently computable polynomial. It is well-known that constant powers of such polynomials also have small roABPs (see e.g. [4, Lemma 2.5]). However, we show that for larger powers, a superpolynomial blow-up in width is unavoidable.

Corollary 4 (𝗋𝗈𝖠𝖡𝖯 powering non-closure).

The following holds over fields of characteristic zero. There exists an n-variate polynomial f computable by an 𝗋𝗈𝖠𝖡𝖯 of width O(n) such that for any d, any 𝗋𝗈𝖠𝖡𝖯 computing fd requires width at least (d+n/2n/2) in every variable order.

 Remark.

We give two example polynomials to prove the hardness of powering for 𝗋𝗈𝖠𝖡𝖯. The first is the elementary symmetric polynomial (this lower bound will also prove Theorem 3) and the second is a quadratic polynomial inspired by the proof of Theorem 1.

Another corollary of Theorem 3 is that computing the resultant and the discriminant is hard for 𝗋𝗈𝖠𝖡𝖯.

Corollary 5 (𝗋𝗈𝖠𝖡𝖯 discriminant non-closure).

The following holds over fields of characteristic zero. For all n, there exists an n-variate polynomial f(𝐱,y) computable by an 𝗋𝗈𝖠𝖡𝖯 of width O(n) such that any 𝗋𝗈𝖠𝖡𝖯 computing the discriminant Discy(f) requires width at least 2Ω(n) in every variable order.

 Remark.

As an immediate consequence of the corollary above, we get that 𝗋𝗈𝖠𝖡𝖯 is not closed under taking resultants.

Related Work

There have been many lines of investigation into roABPs from the point of view of lower bounds [44, 38, 4], PIT algorithms [46, 12, 2, 31, 27, 26, 4, 2, 32, 8, 49], border complexity [21, 22, 14, 13], algebraic meta-complexity [7, 9] and so on.

Our work is closely related to that of Kayal, Nair, and Saha [38], who proved separations between the power of 𝗋𝗈𝖠𝖡𝖯s and multilinear depth-3 circuits. Non-closure results of a similar flavour to ours have also been proved by Saha and Thankey [49, Appendix E.1]. They construct explicit families of polynomials that require 𝗋𝗈𝖠𝖡𝖯 of exponential size, but arise from applying invertible linear transformations to sparse polynomials f that have linear 𝗋𝗈𝖠𝖡𝖯 complexity. Some of their ideas, such as those involving the use of expander graphs, also appear in our work.

Similar separations between roABPs and other models (such as read-twice ABPs) were also addressed in the work of Anderson, Forbes, Saptharishi, Shpilka and Volk [4]. We re-prove a result from this work separating depth-2 algebraic circuits (products of linear polynomials) from roABPs in order to understand the roABP complexity of some explicit symmetric functions. Our lower bound is proved for the specific case of the determinant of a Circulant matrix, which is a naturally occurring mathematical object and hence may be independently interesting.

1.2 Proof Techniques

The main technique for understanding the roABP complexity of a polynomial is a characterization due to Nisan [44], who showed that the roABP complexity of a polynomial f (or more precisely the width of the smallest roABP computing f) in a given order is captured by the ranks of certain matrices related to f, also known as the evaluation dimension of f (formally defined in Definition 11). We also heavily rely on this notion in our work.

Factor non-closure

To construct our examples of polynomials that are efficiently computable by roABPs but hard to factor, we start with an analogous construction for a weaker setting, that of sparse polynomials. The following is a well-known construction due to [57, Example 5.1].

f(x1,,xn)=i=1n(xid1)=i=1n(xi1)i=1n(1+xi++xid1)g(x1,,xn).

Note that the polynomial f has 2n monomials while its factor g has dn monomials. This thus yields an example of a polynomial whose factors have many more monomials than the polynomial itself.

We would like to extend this to the setting of 𝗋𝗈𝖠𝖡𝖯. Unfortunately, the example above does not work as is, as the polynomial g is a product of univariate polynomials and hence has a small roABP. Our idea is to “lift” this sparsity lower bound to an roABP lower bound.

The basic idea of lifting, which has proven powerful in the area of Boolean complexity [20] and also Algebraic Proof complexity [28], is to start with a function f that is hard for a simpler computational model (in this case sparse polynomials) and convert it to a function g that is difficult for a much more powerful model by replacing the variables of f by functions (typically called “gadgets”) in a small number of new variables to obtain g. A version of this idea can be used to lift degree lower bounds on the multilinear representation for some functions to lower bounds for algebraic proof systems based on roABPs [28].

Inspired by [28], we replace the variables of the polynomial f by quadratic multilinear monomials in a new set of variables y1,,ym where m=Ω(n). We can associate this replacement with an undirected graph G on m vertices and n edges. We show that, as long as G is a sufficiently good constant-degree expander, the corresponding “lifted” polynomials fG and gG are easy and hard respectively for roABPs with similar parameters to the case of sparse polynomials.

The crucial property of expander graphs that allows us to prove a lower bound on gG is the Expander Mixing lemma. It can be used to show that given any balanced partition of the vertices of G, there is a large induced matching between the two sets in the parts. This allows us to find a large identity matrix as a submatrix of the evaluation matrix of gG, leading to strong bounds on its evaluation dimension.

The complexity of powering

We give two examples to demonstrate that roABPs are not closed under powering.

The first is a quadratic polynomial g whose monomials again correspond to a constant-degree expander graph as in the previous result. The Expander Mixing lemma can again be used to argue that large powers of g have large evaluation dimension.

The second example is just an elementary symmetric polynomial. Symmetric polynomials are particularly natural to study in the setting of roABPs, since the polynomials have the same complexity under any variable ordering. In particular, studying the complexity of a symmetric polynomial turns into understanding the ranks of combinatorially defined matrices. In the setting of a power of the elementary symmetric polynomial, we are able to show that this matrix has large rank.

No Bläser-Jindal type results for roABPs

Already the example of the elementary symmetric polynomial above shows that for the simple polynomial f(y1,,yn)=yid, the symmetric polynomial f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn) is hard to compute for roABPs.

To prove a converse result, we use the symmetric polynomial fsymd=0n𝖤𝖲𝗒𝗆nd(x1k,,xnk). In this case, we need to understand the complexity of the polynomial f (such as f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn)=fsym). It turns out that the polynomial f in this case is completely understood [1] and is closely related to the determinant of the Circulant matrix. To prove the lower bound, we prove an roABP lower bound on this determinant, which we believe is independently interesting.

Outline

We begin with preliminaries in Section 2. Section 3 contains the proof of Theorem 1. In Section 4, we prove Theorem 2 and Theorem 3. Finally, in Section 5, we prove all the corollaries.

2 Notations and Preliminaries

Throughout the paper, we will use a growing parameter n>0 to denote the number of variables in the polynomial. Let 𝐱=(x1,,xn) be the set of indeterminates. A monomial of the form x1e1xnen is denoted as 𝐱𝐞, where 𝐞=(e1,,en)n. The degree of a monomial 𝐱𝐞 is defined as deg(𝐱𝐞)e1++en. The degree of a polynomial f𝔽[x1,,xn] is defined as the maximum degree of its constituent monomials. We use coef𝐱𝐞(f) to denote the coefficient of the monomial 𝐱𝐞 in f.

2.1 Read-Once Oblivious Algebraic Branching Programs (𝗿𝗼𝗔𝗕𝗣)

Our model of interest arises as a natural restriction of Algebraic Branching Programs (ABPs), which we describe next. An Algebraic Branching Program (ABP) is a layered and directed graph with a source vertex s and a sink vertex t. All edges connect vertices from layer i to i+1. Further, the edges are labeled with affine polynomials over the underlying field 𝔽. For every path γ from s to t, wt(γ) is the product of labels on the edges of the path γ. The polynomial computed by the ABP is defined as

fpath γ:stwt(γ).

The depth of an 𝖠𝖡𝖯 is defined as the number of layers in the graph, and the width is the maximum number of nodes in a layer across the graph. The number of vertices used in the graph is the size of the 𝖠𝖡𝖯. The 𝗋𝗈𝖠𝖡𝖯 model is a restriction of 𝖠𝖡𝖯s, which we define below.

Definition 6 (𝗋𝗈𝖠𝖡𝖯).

Let n be arbitrary and fix a permutation π:[n][n]. An 𝗋𝗈𝖠𝖡𝖯 in the order π computing an n-variate polynomial f(𝐱) is an 𝖠𝖡𝖯 where in the i-th layer the edge labels are univariate polynomials over xπ(i).

The size of an 𝗋𝗈𝖠𝖡𝖯 is defined as the number of vertices it contains, and the width is defined as the maximum number of vertices in any layer.

In a foundational work, Nisan [44] characterized the complexity of an 𝖠𝖡𝖯 in the non-commutative setting with the rank of certain matrices. Remarkably, the characterization extends to 𝗋𝗈𝖠𝖡𝖯 as well. We define the relevant matrix to formally state this characterization.

Definition 7 (Nisan Matrix).

Consider an n-variate polynomial f(𝐱) and a variable partition YZ={x1,,xn}. The Nisan matrix of f with respect to Y,Z, denoted as MY,Z(f), is the matrix whose rows are indexed by monomials mY over Y and whose columns are indexed by monomials mZ over Z. Its entry at (mY,mZ) is defined as

MY,Z(f)[mY,mZ]=coefmYmZ(f).

Historically, the Nisan Matrix has also been referred to as the coefficient matrix or partial derivative matrix. The width of an 𝗋𝗈𝖠𝖡𝖯 computing a polynomial f can be exactly characterized by the rank of the Nisan matrix of f [29, Lemma 4.5.8].

Theorem 8 (𝗋𝗈𝖠𝖡𝖯 characterization).

Let f(𝐱) be an n-variate polynomial, and fix a permutation π on variables. For each i[n], consider the partition Yi{π(x1),,π(xi)} and Zi={π(xi+1),,π(xn)}, and let MYi,Zi(f) denote the corresponding Nisan matrix.
The width of the smallest 𝗋𝗈𝖠𝖡𝖯 computing f in the order π is exactly maxi[n]rank(MYi,Zi(f)). Moreover, the size of the smallest 𝗋𝗈𝖠𝖡𝖯 is exactly i[n]rank(MYi,Zi(f)).

We next prove a lemma to demonstrate the usefulness of the 𝗋𝗈𝖠𝖡𝖯 characterization, which will be used in our later proofs.

Observation 9.

Consider an n-variate polynomial as follows

fi[n](1+xi+xi2++xid1).

Let Y={y1,,yn} and Z={z1,,zn} be disjoint set of variables. Define a 2n-variate polynomial f~f(y1z1,,ynzn). The rank of the Nisan matrix MY,Z(f~) is dn.

 Remark 10.

Note that f itself can be computed by a constant width 𝗋𝗈𝖠𝖡𝖯 in any order.

Proof.

Observe that for every monomial mY over Y such that each variable has degree at most d1 in mY, there is a unique monomial mZ of the same form over Z such that coefmYmZ(f~) is not zero, and reciprocally, for any monomial mZ over Z there is a unique monomial mY. Consequently, the Nisan matrix MY,Z(f~) is a permutation matrix, and hence has rank dn.

Evaluation Dimension

An alternative perspective on the Nisan matrix was introduced by Saptharishi [27, Section 6]. As we will see in our proofs, this viewpoint often makes it easier to reason about 𝗋𝗈𝖠𝖡𝖯 complexity.

Definition 11 (Evaluation Dimension).

Let f(𝐱) be an n-variate polynomial on X={x1,,xn} over a field 𝔽, and a subset of variables Y and ZX\Y. The evaluation dimension of f with respect to the partition YZ is defined as

evalDimY,Z(f)rank({f(Y,𝐚)𝐚𝔽|Z|}).

Over large fields, the evaluation dimension is equivalent to the rank of the Nisan matrix. However, this equivalence does not hold when restricting the evaluation points, e.g. to the Boolean cube. Nevertheless, the evaluation dimension is always a lower bound of the rank of the Nisan matrix ([50, Lemma 11.9], and see also [29, Corollary 4.5.12]).

Theorem 12.

Let f(𝐱) be an n-variate polynomial, and fix a permutation π on variables. For a variable partition YZ with Y={xπ(1),,xπ(i)} and Z={xπ(i+1),,xπ(n)}, any 𝗋𝗈𝖠𝖡𝖯 that computed f in the order π has width at least evalDimY,Z(f).

Conversely, if the field 𝔽 is infinite, there is a 𝗋𝗈𝖠𝖡𝖯 computing f of width evalDimY,Z(f).

2.2 Elementary Symmetric Polynomials

Symmetric polynomials are those that are invariant under any permutation of the variables. A fundamental and well-studied family within symmetric polynomials is the elementary symmetric polynomials, which are defined as follows.

Definition 13 (Elementary Symmetric Polynomial).

The elementary symmetric polynomial of degree d, on variables x1,,xn, is defined as

𝖤𝖲𝗒𝗆nd(x1,,xn)1i1<<idnxi1xid.

Whenever clear from the context, we write ed𝖤𝖲𝗒𝗆nd to denote the degree-d elementary symmetric polynomial in n variables. A more convenient way to define these polynomials is via the following generating functions:

i[n](1+xit)=i=0n𝖤𝖲𝗒𝗆ni(𝐱)ti. (1)

These polynomials are called elementary because they form the fundamental building blocks for all symmetric polynomials. For any n-variate polynomial f(𝐱), we define the n-variate symmetric polynomial f𝗌𝗒𝗆f(e1,,en).

Theorem 14 (Fundamental Theorem of Symmetric Polynomials).

Let R be any commutative ring, and let gR[x1,,xn] be a symmetric polynomial. Then there exists a unique polynomial fR[y1,,yn] such that

g=f𝗌𝗒𝗆f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn).

We refer to [39, Theorem IV.6.1] for the proof of Fundamental Theorem of Symmetric Polynomials (see also [15]). We will also need the following variable partitioning lemma, which is a special case of [42, Theorem 1.1].

Lemma 15 (𝖤𝖲𝗒𝗆 Variable Partition).

Let YZ be a partition of the variables. Then,

𝖤𝖲𝗒𝗆|YZ|d(Y,Z)=i=0d(𝖤𝖲𝗒𝗆|Y|i(Y)𝖤𝖲𝗒𝗆|Z|di(Z)).
Proof.

Let Y={y1,,ym} and Z={z1,,zn}. Then using Equation 1 we can write,

i=0m𝖤𝖲𝗒𝗆mi(Y)ti =i=1m(1+yit)
i=0n𝖤𝖲𝗒𝗆ni(Z)ti =i=1n(1+zit)

Taking the product of the two polynomials above, and comparing the coefficients of td on both sides proves the lemma.

As a direct consequence of extracting the coefficient of td from Equation 1, Shpilka and Wigderson [51] (crediting Ben-Or) presented the following identity for elementary symmetric polynomials, which yields a near-optimal 𝗋𝗈𝖠𝖡𝖯 of width O(n) computing 𝖤𝖲𝗒𝗆nd in any variable order:

Proposition 16 ([51, Theorem 5.1]).

For any n and dn, let ω be a primitive n-th root of unity. There exist β0,,βn1 such that

𝖤𝖲𝗒𝗆nd(x1,,xn)=0j<nβj(1+ωjx1)(1+ωjx2)(1+ωjxn).
 Remark 17.

𝖤𝖲𝗒𝗆nd can also be computed by a provably tight 𝗋𝗈𝖠𝖡𝖯 of width min(d+1,nd+1) in any variable order using only coefficients 0 and 1 (see [45, Construction 1.2]).

2.3 Resultant and Discriminant

We recall the definitions and properties of resultant and discriminant from factorization literature. We encourage readers to refer to [56, Chapter 6] for a more detailed textbook treatment of these concepts.

Definition 18 (Resultant).

Consider two n-variate polynomials f,g𝔽[𝐱][y] as follows:

fi=0d1fi(𝐱)yiandgi=0d2gi(𝐱)yi.

Define the Sylvester matrix of f and g as the following (d1+d2)×(d1+d2) matrix:

𝐒y(f,g)=(fd1gd2fd11fd1gd21gd2fd11gd21fd1gd2f0fd11g0gd21f0g0f0g0)

Then the resultant of the two polynomials with respect to y is defined as the determinant of the Sylvester matrix as:

Resy(f,g)𝖣𝖾𝗍(𝐒y(f,g)).

The resultant of two polynomials is non-zero if and only if their gcd is 1. A well-known case of resultant relevant for factoring algorithms is the discriminant.

Definition 19 (Discriminant).

Consider a n-variate polynomial f. The discriminant with respect to y of f is defined as the resultant, with respect to y, of f and its y-derivative, i.e.,

Discy(f)Resy(f,yf).

The following well-known observation will be useful in the analyses of complexity of the resultant and the discriminant.

Observation 20 (see [19, Chapter 3]).

Let f=i[n](yαi) and g=i[m](yβi) be two univariate polynomials. Then the resultant of f and g with respect to y is given by

Resy(f,g)=i[n]g(αi).

3 𝗿𝗼𝗔𝗕𝗣 Factor Non-Closure

To prove Theorem 1, we need a polynomial of low 𝗋𝗈𝖠𝖡𝖯 complexity, that has a factor of high 𝗋𝗈𝖠𝖡𝖯 complexity. We will use explicit expander graphs for this purpose. The only property we require from the expander graph is that, for any sufficiently large partition of its vertex set into two parts, it contains a large induced matching between the two parts.

Lemma 21 (Induced Matching Lemma).

For every n there exists a constant degree graph Gn=(V,E) on n vertices such that the following holds: for any partition (S,T) of V with |S|=εn and |T|=(1ε)n where ε[13,23], the graph contains Ω(n) edges between S and T that form an induced matching.

Proof Sketch.

There exists an absolute constant δ(0,1) such that, for any k with k1, we can construct explicit k-regular expander graphs Gn=(V,E) whose second-largest eigenvalue is at most kδ; see [47].

When k is chosen to be sufficiently large such that the second-largest eigenvalue of Gn is strictly smaller than k/3, then the lemma follows as an easy consequence of the Expander Mixing Lemma [3] (see also [33, Lemma 2.5]). See, for example, [34, Claim 4].

Define an n-variate polynomial PG associated with constant degree graph Gn=(V,E) guaranteed by Lemma 21 as follows:

PG(i,j)E((xixj)d1). (2)

Since the degree of the graph is constant, the sparsity of PG is 2|E|=2O(n)=:w. Therefore, PG can be computed by an 𝗋𝗈𝖠𝖡𝖯 of width w in every variable order. To prove the hardness of its factor, consider the following polynomial QG:

QG(i,j)E(1+(xixj)+(xixj)2++(xixj)d1). (3)

It is well known that (1+x+x2++xd1)(x1)=xd1. Using the identity, we immediately obtain

PG=QG(i,j)E(xixj1).
Lemma 22.

The polynomial QG defined in Equation 3 requires an 𝗋𝗈𝖠𝖡𝖯 of width dΩ(n) in every variable order.

Proof.

Let π be any variable order on the variables X={x1,,xn}, and consider the partition Y={xπ(1),,xπ(n/2)} and Z=XY.

Let Y and Z also denote the partition of vertices of Gn. Then from Lemma 21, we know there exists an induced matching between Y and Z of size Ω(n). Define

f~ (i,j)(1+(xixj)+(xixj)2++(xixj)d1)
=i[t](1+(yizi)+(yizi)2++(yizi)d1),

where for every i[t], yi is a variable in Y and zi is a variable in Z, and t=Ω(n). Here we have used the fact that is an induced matching. In particular, f~ is obtained from QG by setting to zero the variables which are not in the matching . Hence, the rank of the Nisan matrix can only decrease. Finally, by Observation 9,

rank(MY,Z(QG))rank(MY,Z(f~))dΩ(n).

We obtain the claimed lower bound for width of 𝗋𝗈𝖠𝖡𝖯 computing QG by Theorem 8.

We will now use the discussion so far to give the complete proof of the factor non-closure result.

Theorem 1 (𝗋𝗈𝖠𝖡𝖯 factor non-closure). [Restated, see original statement.]

The following holds over any field. Let n be a parameter and dn. There exists an n-variate polynomial f of degree d computable by an 𝗋𝗈𝖠𝖡𝖯 of width w2O(n), such that one of its (irreducible) factors g requires an 𝗋𝗈𝖠𝖡𝖯 of width wΩ(logd) in every variable order.

Proof.

Consider an n-variate polynomial gQGn1+z, where z is an auxiliary variable and QGn1 is defined as in Equation 3 using a constant-degree graph Gn1. We then define

fg(i,j)E(xixj1)=PG+z(i,j)E(xixj1).

As argued after Equation 2, both PG and (i,j)E(xixj1) have sparsity 2O(n) and hence we can compute them by an 𝗋𝗈𝖠𝖡𝖯 of width w=2O(n) in every variable order. Therefore, f itself admits an 𝗋𝗈𝖠𝖡𝖯 of width w in every variable order.

Observe that g is an irreducible polynomial because it is linear in the auxiliary variable z.444See [57, Example 5.1] where the hardness is lifted to irreducible factor by considering g=QG+n. Further, by Lemma 22, any 𝗋𝗈𝖠𝖡𝖯 computing QG+z must have width at least dΩ(n) in every variable order. Since dn, the claimed width lower bound for 𝗋𝗈𝖠𝖡𝖯 computing g follows.

4 𝗿𝗼𝗔𝗕𝗣 Complexity of Symmetric Polynomials

In the following two sections we prove Theorem 2 and Theorem 3 along with their corollaries.

4.1 𝒇𝘀𝘆𝗺 is easy, but 𝒇 is hard

In this section, we work over fields of characteristic zero. For the proof of Theorem 2, we consider f𝗌𝗒𝗆=d=0n𝖤𝖲𝗒𝗆nd(x1k,,xnk) for a suitable choice of k[n] to be fixed later.

Let us consider the polynomial

g(y1,,yn,t,z0,,zk1)j=0k1(1+i[n]yi(tzj)i).

The polynomial g is symmetric in the variables z0,,zk1. So by Theorem 14, there exists a polynomial g~[y1,,yn,t,z0,,zk1] such that g(𝐲,t,𝐳)=g~(𝐲,t,e1(𝐳),,ek(𝐳)). Notice that if ω is a k-th primitive root of the unity, Equation 1 implies

{ei(ω0,,ωk1)=0for 1i<k,ek(ω0,,ωk1)=1.

Let us define

f(y1,,yn)g~(𝐲,1,0,,0,1)[𝐲]. (4)

The previous paragraph ensures that for any k-th primitive root of the unity ω, we have

f(y1,,yn)=j=0k1(1+i[n]yiωji). (5)

The following lemma shows that f is indeed the unique polynomial inducing the symmetric polynomial d=0n𝖤𝖲𝗒𝗆nd(x1k,,xnk). The following is an argument in [1], which we reproduce here for completeness.

Lemma 23 (Circulant Polynomial).

For any n and odd positive integer kn,

f𝗌𝗒𝗆(𝐱)f(e1(𝐱),,en(𝐱))=d=0n𝖤𝖲𝗒𝗆nd(x1k,,xnk).
Proof.

Let ω be a k-th primitive root of the unity. Using the factorization identity (1tk)=j(1ωjt) together with Equation 1 we obtain the following:

i[n](1xik(t)k) =i=1nj=0k1(1+ωjxit)
=j=0k1(1+i=1nei(𝐱)(ωjt)i).

By Equation 5 and instantiating t by 1, we obtain

i[n](1+xik)=f(e1(𝐱),,en(𝐱)).

We call the polynomial f a circulant polynomial because it is closely related to the determinant of a Circulant matrix.555Specifically, in the case k=n, the homogeneous component of degree k of the polynomial f is exactly the determinant of the circulant matrix of first row (x1,,xn).

In the following lemma, we show that the polynomial f is hard for 𝗋𝗈𝖠𝖡𝖯 in every variable order over any field 𝔽 of characteristic 0.

Lemma 24.

For any prime k with 2kn, the n-variate polynomial f defined in Equation 4 requires 𝗋𝗈𝖠𝖡𝖯 width at least 2(k1)/2 in any variable order.

Proof.

By instantiating a variable, the roABP width can only decrease. So it is sufficient to consider f(y1,,yk)f(y1,,yk,0,,0).

By the standard evaluation-dimension lower bound for 𝗋𝗈𝖠𝖡𝖯, it suffices to show the following. For any variable order π and the variable partition (U,V) where U={yπ(1),,yπ((k1)/2)} and V={yπ((k1)/2+1),,yπ(k)} we have

evalDimU,V(f)=rank({f(𝐮,𝐚)𝐚𝔽|V|}) 2Ω(k).

Re-writing Equation 5 in terms of variable partition (U,V), we have

f(𝐮,𝐚)=j=0k1(j(𝐮)+j(𝐚)+1), (6)

where j(𝐮) and j(𝐯) are linear polynomials in U and V, respectively.

Arranging the coefficients in the linear forms {j(𝐮)+j(𝐯)}j as the rows of a k×k matrix yields a matrix M whose (j,i)-th entry is ωjπ(i) for j{0,,k1} and i[k]. Note that M can be obtained from the standard k×k DFT matrix (ωji)i,j[k] by permuting columns.

Figure 1: The matrix M with (j,i)-th entry ωjπ(i), corresponding to the natural variable order. The highlighted submatrices MA and MB are used in the analysis of evaluation dimension.

When k is prime, Chebotarev’s theorem on roots of unity [53] states that every square submatrix of the DFT matrix (and hence also M) is nonsingular; see also [54, Lemma 1.3] and [30]. Since |U|=(k1)/2, we can fix a subset A{0,,k1} of size (k1)/2 such that the set {j(𝐮)}jA corresponds to a square submatrix MA of M. Such a submatrix MA is non-singular due to Chebotarev’s theorem. Consequently, the set {j(𝐮)}jA is linearly independent. We can assume that i(𝐮)=ui for iA, since an invertible linear transformation on the variables in U does not change the evaluation dimension evalDimU,V(f).

Consider any BA. By Chebotarev’s theorem, the submatrix MB with rows indexed by B and columns v2,,v|B|+1 is invertible (the choice |A|=(k1)/2 ensures that V contains enough variables). So, for any b𝔽, there is a unique point βB,b𝔽|B| such that for any j in B, j(b,βB,b,0,,0)+1=0. Similarly, for any other row index ȷ~{0,,k1}B, there exists a unique point γB,ȷ~𝔽|B|+1 such that for any j in B{ȷ~}, we have j(γB,ȷ~,0,,0)+1=0. It follows that γB,ȷ~ is of the form (bȷ~,βB,bȷ~) for a particular bȷ~𝔽.

Since 𝔽 is infinite, we can choose b in 𝔽 outside of {bȷ~ȷ~{0,,k1}B}, and define αB(b,βB,b,0,,0). For any j in {0,,k1}:

j(αB)+1=0jB.

Consequently f(𝐮,αB)=(jBuj)(j[k]B(j(U)+cB,j)) where the (cB,j)jB are non-zero constants. Since for each B, f(𝐮,αB) has a distinct lowest degree monomial (jBuj), the set {f(𝐮,αB)BA} is linearly independent. Therefore,

evalDimU,V(f)) dim{f(𝐮,αB)BA}=2(k1)/2.

By the evaluation-dimension lower bound of Theorem 12, any 𝗋𝗈𝖠𝖡𝖯 computing f (in any order) must have width at least 2(k1)/2.

 Remark 25.

A close look at the above proof reveals that the lower bound also applies to the circulant polynomial j=0k1(i=1nyiωij) which is exactly the determinant of the circulant matrix.

The lower bound established in Lemma 24 serves as the key technical ingredient needed to prove Theorem 2.

Theorem 2. [Restated, see original statement.]

The following holds over fields of characteristic zero. Let n be a parameter. There exists an n-variate polynomial f such that the symmetric polynomial f𝗌𝗒𝗆f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn) is computable by an 𝗋𝗈𝖠𝖡𝖯 of constant width in every variable order, but any 𝗋𝗈𝖠𝖡𝖯 computing f in any variable order must have width 2Ω(n).

Proof.

Fix k to be a prime number between n/2 and n. We consider the symmetric polynomial f𝗌𝗒𝗆d=0n𝖤𝖲𝗒𝗆nd(x1k,,xnk). By applying Equation 1 with each xi replaced by xik, we obtain that f𝗌𝗒𝗆 admits an 𝗋𝗈𝖠𝖡𝖯 of constant width in every variable order. Moreover, by Lemma 24, the width of any 𝗋𝗈𝖠𝖡𝖯 computing f is at least 2(k1)/2=2Ω(n).

4.2 𝒇 is easy, but 𝒇𝘀𝘆𝗺 is hard

To prove Theorem 3, it suffices to show the following technical lemma, which shows that taking powers of elementary symmetric polynomials is hard for 𝗋𝗈𝖠𝖡𝖯s. This lemma also implies Corollary 4 from the introduction. In Subsection 5.1, we will present an alternative proof of Theorem 3 using a quadratic polynomial based on graph-based polynomial from Section 3.

Lemma 26 (Powers of 𝖤𝖲𝗒𝗆).

Let kn/2. Any 𝗋𝗈𝖠𝖡𝖯 computing (𝖤𝖲𝗒𝗆nk)d in any variable order requires width at least (k+dk).

Theorem 3. [Restated, see original statement.]

The following holds over fields of characteristic zero. Let n be a parameter. There exists an n-variate polynomial f computable by an 𝗋𝗈𝖠𝖡𝖯 of constant width such that its respective symmetric polynomial f𝗌𝗒𝗆=f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn) requires an 𝗋𝗈𝖠𝖡𝖯 of width 2Ω(n) in every variable order.

Proof.

Let k=n/2. Consider the polynomial f(x1,,xn)=xkk. It is easy to see that f can be computed by an 𝗋𝗈𝖠𝖡𝖯 of constant width. However, by Lemma 26, any 𝗋𝗈𝖠𝖡𝖯 computing the symmetrisation

f𝗌𝗒𝗆=f(𝖤𝖲𝗒𝗆n1,,𝖤𝖲𝗒𝗆nn)=(𝖤𝖲𝗒𝗆nk)k

must have width at least

(2n/2n/2)=Ω(2n/n).

 Remark 27.

We recall that 𝖤𝖲𝗒𝗆nn/2 can be expressed as a sum of n many products of univariate polynomials (see Proposition 16). Consequently, using the multinomial theorem, it follows that (𝖤𝖲𝗒𝗆nn/2)n/2 can be expressed as a sum of at most O(21.5n) many products of univariate polynomials. Hence, the bound we obtain in Theorem 3 is almost optimal.

Proof of Lemma 26.

Assume that (ek)d is computed by an 𝗋𝗈𝖠𝖡𝖯 of width w and variable order π. Let Y={xπ(1),,xπ(k)} and Z=XY be a partition of the variables X. By Theorem 12, we know that

wevalDimY,Z((ek)d).

Using Lemma 15, and the multinomial theorem, we can write the powers of the elementary symmetric polynomial ek as follows:

(ek(X))d =(t=0ket(Y)ekt(Z))d
=t0++tk=dti0(dt0,,tk)(e0t0(Y)ektk(Y))(ek0t0(Z)ekktk(Z)). (7)

To argue about the evaluation dimension of (ek(X))d, we will need the following elementary fact from linear algebra.

Proposition 28.

If a matrix M=i=1ruiviT where {u1,,ur} and {v1,,vr} are linearly independent sets of vectors, then M has rank exactly r.

To use the above fact, we note that the algebraic independence of the elementary symmetric polynomials (a consequence of Theorem 14) implies that the sets

E={e0t0(Y)ektk(Y):t0++tk=d} and E~={ek0t0(Z)ekktk(Z):t0++tk=d}

are both linearly independent sets of polynomials (E~ can be obtained from E by just changing the underlying variable). Further, each term on the right-hand side of Equation 7 (corresponding to a tuple (t0,,tk) summing to d) has an evaluation matrix that is the outer product of the coefficient vectors of the corresponding polynomials in E and E~, scaled by a suitable multinomial coefficient (which is non-zero because we have assumed that the characteristic of the underlying field is 0). This implies that the evaluation matrix of (ek(X))d has rank exactly the number of terms which is (k+dk).

5 Non-closure corollaries for 𝗿𝗼𝗔𝗕𝗣

We will now give the proofs of corollaries stated in Subsection 1.1. We use observations from the earlier sections to show that operations such as powering, computing resultant, and discriminant can be hard for 𝗋𝗈𝖠𝖡𝖯.

5.1 Hardness of Powering: a second example

In this section, we give a second proof of (a slightly weaker form of) Corollary 4. Note that we already proved this in the form of Lemma 26. By Proposition 16 and the following remark, we know that 𝖤𝖲𝗒𝗆2nn admits an 𝗋𝗈𝖠𝖡𝖯 of width O(n) in any variable order. On the other hand, any 𝗋𝗈𝖠𝖡𝖯 that computes (𝖤𝖲𝗒𝗆2nn)d must have width at least (n+dn).

Inspired by the graph-based polynomial which was used to prove factor non-closure in Section 3, we can even define a quadratic polynomial Q and prove that powering Q is hard for this polynomial 𝗋𝗈𝖠𝖡𝖯. The lower bound we obtain is slightly weaker, but the example is even simpler since Q is just a quadratic polynomial, as opposed to the high-degree and high-sparsity elementary symmetric polynomial.

Corollary 29 (Variant of Corollary 4).

The following holds over fields of characteristic zero. There exists an n-variate quadratic polynomial Q computable by an 𝗋𝗈𝖠𝖡𝖯 of width O(n) such that for any d, any 𝗋𝗈𝖠𝖡𝖯 computing Qd requires width at least (d+mm) in every variable order where m=Ω(n).

Proof.

Let G=(V,E) be a constant degree graph on n vertices such that Lemma 21 holds. Define the quadratic polynomial:

QG=(i,j)Exixj (8)

where variables xi correspond to the vertices of G. It is easy to observe that QG can be computed by an 𝗋𝗈𝖠𝖡𝖯 of width |E|=O(n) in any variable order. We will prove that any 𝗋𝗈𝖠𝖡𝖯 computing QGd must have large width.

Let π be any variable order on the variables X={x1,,xn}, and consider the partition Y={xπ(1),,xπ(n/2)} and Z=XY. Let Y and Z also denote the partition of vertices on G. By Lemma 21, there exists an induced matching between Y and Z of size Ω(n). By renaming the variables if necessary we assume that the matching is between the vertices corresponding to yi and zi where i[t] and t=Ω(n).

Define the polynomial:

Q~d=((i,j)xixj)d=(i[t]yizi)d.

In particular, Q~d is obtained from QGd by setting to zero the variables which are not in the matching . Hence, the rank of the Nisan matrix corresponding to Q~ is a lower bound on the evaluation dimension of Q w.r.t. the partition (Y,Z).

By construction, for every monomial mY of degree exactly d over Y, there exists a unique monomial mZ over Z such that the coefficient of mYmZ in Q~d is nonzero (cf. Observation 9). Therefore,

rank(MY,Z(QGd))rank(MY,Z(Q~d))(d+t1t1).

Applying Theorem 8, we obtain the desired lower bound on the width of any 𝗋𝗈𝖠𝖡𝖯 computing QGd.

5.2 Hardness of computing resultant and discriminant

We design a polynomial that is simple for 𝗋𝗈𝖠𝖡𝖯, but which turns out to be difficult for 𝗋𝗈𝖠𝖡𝖯 when one computes its discriminant. This, in turn, immediately implies that computing resultant is also hard for 𝗋𝗈𝖠𝖡𝖯.

Corollary 5 (𝗋𝗈𝖠𝖡𝖯 discriminant non-closure). [Restated, see original statement.]

The following holds over fields of characteristic zero. For all n, there exists an n-variate polynomial f(𝐱,y) computable by an 𝗋𝗈𝖠𝖡𝖯 of width O(n) such that any 𝗋𝗈𝖠𝖡𝖯 computing the discriminant Discy(f) requires width at least 2Ω(n) in every variable order.

Proof.

Let g be an (n1)-variate polynomial to which the lower bound in Corollary 4 is applicable. Fix any d=Ω(n). Define

fydg(𝐱)y.

Then we have yf=dyd1g. It is easy to see that the roots of f are α0=0 and αi=ωig1/(d1) for 1id1, where ω is a primitive (d1)-th root of unity. Here we work over a suitable field extension of the base field 𝔽 to ensure that we have an (d1)-th root of unity.

The discriminant of f is defined as the resultant of f and yf with respect to y, i.e., Discy(f)=Resy(f,yf). Then using Observation 20 we can compute:

Discy(f) =i=0d1yf(αi)=gi=1d1(d1)g
=(d1)d1gd.

Thus, computing the discriminant of f amounts to powering the polynomial g, for which we have the required lower bound by Corollary 4. The upper bound on the 𝗋𝗈𝖠𝖡𝖯 complexity of f follows from the one for g.

6 Conclusions and Open Problems

In this work, we proved that a width-w 𝗋𝗈𝖠𝖡𝖯 computes a polynomial whose irreducible factor requires 𝗋𝗈𝖠𝖡𝖯s of width at least wlogd, yielding a quasipolynomial separation. This showed that 𝗋𝗈𝖠𝖡𝖯s are not closed under factoring (see Section 3). A natural next step is to search for polynomials that exhibit an exponential separation between the 𝗋𝗈𝖠𝖡𝖯 complexity of a polynomial and that of its factor.

Our non-factor closure proof relied on the idea that polynomials that are hard for the simpler sparse model but easy for 𝗋𝗈𝖠𝖡𝖯s can be transformed, using simple gadgets, into polynomials that are hard even for 𝗋𝗈𝖠𝖡𝖯s. This raises an intriguing question about the scope of such hardness lifting. Specifically, given a polynomial f(𝐱) of sparsity s, can we always lift using a gadget ϕ such that the composed polynomial f(ϕ𝐱) requires an 𝗋𝗈𝖠𝖡𝖯 of size Ω(s)?

One consequence of our study of graph-based polynomials and symmetric compositions is the proof that powering is hard for 𝗋𝗈𝖠𝖡𝖯s (see Subsection 5.1). This naturally raises the question in the other direction: does there exist a polynomial fgd that is easy to compute by an 𝗋𝗈𝖠𝖡𝖯, while g is hard for 𝗋𝗈𝖠𝖡𝖯? An affirmative answer would, once again, stand in sharp contrast to other models such as circuits, algebraic branching programs, and formulas, where low complexity of f leads to low complexity of g. Interestingly, the analogous question for sparse polynomials was answered affirmatively for d=2 in classical works by Rényi [48] and Erdős [24], and was subsequently extended to arbitrary d in later works [55, 18].

References

  • [1] achille hui. Evaluating elementary symmetric polynomial at nth powers. Mathematics Stack Exchange, 2021. URL: https://math.stackexchange.com/a/4259092.
  • [2] Manindra Agrawal, Rohit Gurjar, Arpita Korwar, and Nitin Saxena. Hitting-sets for ROABP and sum of set-multilinear circuits. SIAM J. Comput., 44(3):669–697, 2015. doi:10.1137/140975103.
  • [3] Noga Alon and Fan R. K. Chung. Explicit construction of linear sized tolerant networks. Discret. Math., 306(10-11):1068–1071, 2006. doi:10.1016/J.DISC.2006.03.025.
  • [4] Matthew Anderson, Michael A. Forbes, Ramprasad Saptharishi, Amir Shpilka, and Ben Lee Volk. Identity testing and lower bounds for read-k oblivious algebraic branching programs. ACM Trans. Comput. Theory, 10(1):3:1–3:30, 2018. doi:10.1145/3170709.
  • [5] Robert Andrews and Avi Wigderson. Constant-depth arithmetic circuits for linear algebra problems. In 2024 IEEE 65th Annual Symposium on Foundations of Computer Science—FOCS 2024, pages 2367–2386. IEEE Computer Society, 2024. doi:10.1109/FOCS61266.2024.00138.
  • [6] C. S. Bhargav, Sagnik Dutta, and Nitin Saxena. Improved lower bound, and proof barrier, for constant depth algebraic circuits. In 47th International Symposium on Mathematical Foundations of Computer Science, MFCS 2022, August 22-26, 2022, Vienna, Austria, volume 241 of LIPIcs, pages 18:1–18:16. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2022. doi:10.4230/LIPIcs.MFCS.2022.18.
  • [7] Vishwas Bhargava, Pranjal Dutta, Sumanta Ghosh, and Anamay Tengse. The complexity of order-finding for roabps. CoRR, abs/2411.18981, 2024. doi:10.48550/arXiv.2411.18981.
  • [8] Vishwas Bhargava and Sumanta Ghosh. Improved hitting set for orbit of roabps. Comput. Complex., 31(2):15, 2022. doi:10.1007/S00037-022-00230-9.
  • [9] Vishwas Bhargava and Anamay Tengse. Explicit commutative roabps from partial derivatives. In 44th IARCS Annual Conference on Foundations of Software Technology and Theoretical Computer Science, FSTTCS 2024, December 16-18, 2024, Gandhinagar, Gujarat, India, volume 323 of LIPIcs, pages 10:1–10:15. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2024. doi:10.4230/LIPIcs.FSTTCS.2024.10.
  • [10] Somnath Bhattacharjee, Mrinal Kumar, Shanthanu S. Rai, Varun Ramanathan, Ramprasad Saptharishi, and Shubhangi Saraf. Closure under factorization from a result of Furstenberg. CoRR, abs/2506.23214, 2025. doi:10.48550/arXiv.2506.23214.
  • [11] Somnath Bhattacharjee, Mrinal Kumar, Shanthanu S. Rai, Varun Ramanathan, Ramprasad Saptharishi, and Shubhangi Saraf. Constant-depth circuits for polynomial GCD over any characteristic. CoRR, abs/2506.23220, 2025. doi:10.48550/arXiv.2506.23220.
  • [12] Pranav Bisht and Nitin Saxena. Blackbox identity testing for sum of special roabps and its border class. Comput. Complex., 30(1):8, 2021. doi:10.1007/S00037-021-00209-Y.
  • [13] Markus Bläser, Julian Dörfler, and Christian Ikenmeyer. On the complexity of evaluating highest weight vectors. In 36th Computational Complexity Conference, CCC 2021, July 20-23, 2021, Toronto, Ontario, Canada (Virtual Conference), volume 200 of LIPIcs, pages 29:1–29:36. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2021. doi:10.4230/LIPIcs.CCC.2021.29.
  • [14] Markus Bläser, Christian Ikenmeyer, Meena Mahajan, Anurag Pandey, and Nitin Saurabh. Algebraic branching programs, border complexity, and tangent spaces. In 35th Computational Complexity Conference, CCC 2020, July 28-31, 2020, Saarbrücken, Germany (Virtual Conference), volume 169 of LIPIcs, pages 21:1–21:24. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2020. doi:10.4230/LIPIcs.CCC.2020.21.
  • [15] Ben Blum-Smith and Samuel Coskey. The Fundamental Theorem on Symmetric Polynomials: History’s First Whiff of Galois Theory. College Mathematics Journal, 48(1):18–29, 2017. doi:10.4169/college.math.j.48.1.18.
  • [16] Markus Bläser and Gorav Jindal. On the Complexity of Symmetric Polynomials. In 10th Innovations in Theoretical Computer Science Conference (ITCS 2019), volume 124 of Leibniz International Proceedings in Informatics (LIPIcs), pages 47:1–47:14. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2019. doi:10.4230/LIPIcs.ITCS.2019.47.
  • [17] Chi-Ning Chou, Mrinal Kumar, and Noam Solomon. Closure results for polynomial factorization. Theory Comput., 15:Paper No. 13, 34, 2019. doi:10.4086/toc.2019.v015a013.
  • [18] Don Coppersmith and James Davenport. Polynomials whose powers are sparse. Acta Arith., 58:79–87, 1991. URL: https://eudml.org/doc/206337.
  • [19] David A. Cox, John Little, and Donal O’Shea. Using Algebraic Geometry, volume 185 of Graduate Texts in Mathematics. Springer, New York, NY, 2 edition, 2005. doi:10.1007/b138611.
  • [20] Susanna F. de Rezende, Mika Göös, and Robert Robere. Guest column: Proofs, circuits, and communication. SIGACT News, 53(1):59–82, 2022. doi:10.1145/3532737.3532746.
  • [21] Pranjal Dutta, Prateek Dwivedi, and Nitin Saxena. Demystifying the border of depth-3 algebraic circuits. In 2021 IEEE 62nd Annual Symposium on Foundations of Computer Science (FOCS), pages 92–103, 2022. doi:10.1109/FOCS52979.2021.00018.
  • [22] Pranjal Dutta and Nitin Saxena. Separated borders: Exponential-gap fanin-hierarchy theorem for approximative depth-3 circuits. In 63rd IEEE Annual Symposium on Foundations of Computer Science, FOCS 2022, Denver, CO, USA, October 31 - November 3, 2022, pages 200–211. IEEE, 2022. doi:10.1109/FOCS54457.2022.00026.
  • [23] Zeev Dvir, Amir Shpilka, and Amir Yehudayoff. Hardness-randomness tradeoffs for bounded depth arithmetic circuits. SIAM J. Comput., 39(4):1279–1293, 2009. doi:10.1137/080735850.
  • [24] Paul Erdős. On the number of terms of the square of a polynomial. Nieuw Arch. Wiskunde (2), 23:63–65, 1949. URL: https://users.renyi.hu/˜p_erdos/1949-08.pdf.
  • [25] Michael Forbes. Some concrete questions on the border complexity of polynomials. Video lecture, Workshop on Algebraic Complexity Theory (WACT), Tel Aviv, 2016. URL: https://www.youtube.com/watch?v=1HMogQIHT6Q.
  • [26] Michael A. Forbes, Ramprasad Saptharishi, and Amir Shpilka. Hitting sets for multilinear read-once algebraic branching programs, in any order. In Symposium on Theory of Computing, STOC 2014, New York, NY, USA, May 31 - June 03, 2014, pages 867–875. ACM, 2014. doi:10.1145/2591796.2591816.
  • [27] Michael A. Forbes and Amir Shpilka. Quasipolynomial-time identity testing of non-commutative and read-once oblivious algebraic branching programs. In 2013 IEEE 54th Annual Symposium on Foundations of Computer Science—FOCS 2013, pages 243–252. IEEE Computer Soc., Los Alamitos, CA, 2013. doi:10.1109/FOCS.2013.34.
  • [28] Michael A. Forbes, Amir Shpilka, Iddo Tzameret, and Avi Wigderson. Proof complexity lower bounds from algebraic circuit complexity. Theory Comput., 17:Paper No. 10, 88, 2021. doi:10.4086/toc.2021.v017a010.
  • [29] Michael Andrew Forbes. Polynomial Identity Testing of Read-Once Oblivious Algebraic Branching Programs. PhD thesis, Massachusetts Institute of Technology, 2014. URL: https://dspace.mit.edu/handle/1721.1/89843.
  • [30] P. E. Frenkel. Simple proof of chebotarev’s theorem on roots of unity, 2004. doi:10.48550/arXiv.math/0312398.
  • [31] Zeyu Guo and Rohit Gurjar. Improved explicit hitting-sets for roabps. In Approximation, Randomization, and Combinatorial Optimization. Algorithms and Techniques, APPROX/RANDOM 2020, August 17-19, 2020, Virtual Conference, volume 176 of LIPIcs, pages 4:1–4:16. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2020. doi:10.4230/LIPIcs.APPROX/RANDOM.2020.4.
  • [32] Rohit Gurjar, Arpita Korwar, Nitin Saxena, and Thomas Thierauf. Deterministic identity testing for sum of read-once oblivious arithmetic branching programs. Comput. Complex., 26(4):835–880, 2017. doi:10.1007/S00037-016-0141-Z.
  • [33] Shlomo Hoory, Nathan Linial, and Avi Wigderson. Expander graphs and their applications. Bulletin of the American Mathematical Society, 43(4):439–561, 2006. doi:10.1090/S0273-0979-06-01126-8.
  • [34] Stasys Jukna. Expanders and time-restricted branching programs. Theoretical computer science, 409(3):471–476, 2008. doi:10.1016/j.tcs.2008.09.012.
  • [35] Valentine Kabanets and Russell Impagliazzo. Derandomizing polynomial identity tests means proving circuit lower bounds. Comput. Complexity, 13(1-2):1–46, 2004. doi:10.1007/s00037-004-0182-6.
  • [36] Erich Kaltofen. Factorization of polynomials given by straight-line programs. Adv. Comput. Res., 5:375–412, 1989. URL: https://api.semanticscholar.org/CorpusID:14414372.
  • [37] Erich Kaltofen and Barry M. Trager. Computing with polynomials given by black boxes for their evaluations: Greatest common divisors, factorization, separation of numerators and denominators. J. Symbolic Comput., 9(3):301–320, 1990. doi:10.1016/S0747-7171(08)80015-6.
  • [38] Neeraj Kayal, Vineet Nair, and Chandan Saha. Separation between read-once oblivious algebraic branching programs (roabps) and multilinear depth-three circuits. ACM Trans. Comput. Theory, 12(1):2:1–2:27, 2020. doi:10.1145/3369928.
  • [39] Serge Lang. Algebra, volume 211 of Graduate Texts in Mathematics. Springer, New York, NY, revised 3rd edition, 2002. doi:10.1007/978-1-4613-0041-0.
  • [40] Nutan Limaye, Srikanth Srinivasan, and Sébastien Tavenas. Superpolynomial lower bounds against low-depth algebraic circuits. Journal of the ACM, 72(4):1–35, 2025. doi:10.1145/3734215.
  • [41] Richard J. Lipton and Ken W. Regan. Arithmetic complexity and symmetry. Blog post, Gödel’s Lost Letter and P=NP, 2009. URL: https://rjlipton.com/2009/07/10/arithmetic-complexity-and-symmetry/.
  • [42] Mircea Merca. A convolution for complete and elementary symmetric functions. Aequationes Mathematicae, 86(3):217–229, 2013. doi:10.1007/s00010-012-0170-x.
  • [43] N. Nisan and A. Wigderson. Lower bounds on arithmetic circuits via partial derivatives. In Proceedings of IEEE 36th Annual Foundations of Computer Science, pages 16–25, 1995. doi:10.1109/SFCS.1995.492458.
  • [44] Noam Nisan. Lower bounds for non-commutative computation (extended abstract). In Proceedings of the 23rd Annual ACM Symposium on Theory of Computing, May 5-8, 1991, New Orleans, Louisiana, USA, pages 410–418. ACM, 1991. doi:10.1145/103418.103462.
  • [45] C. Ramya and Anamay Tengse. On finer separations between subclasses of read-once oblivious abps. In 39th International Symposium on Theoretical Aspects of Computer Science (STACS 2022), volume 219 of Leibniz International Proceedings in Informatics (LIPIcs), pages 53:1–53:23, Dagstuhl, Germany, 2022. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.STACS.2022.53.
  • [46] Ran Raz and Amir Shpilka. Deterministic polynomial identity testing in non-commutative models. Comput. Complex., 14(1):1–19, 2005. doi:10.1007/S00037-005-0188-8.
  • [47] Omer Reingold, Salil P. Vadhan, and Avi Wigderson. Entropy waves, the zig-zag graph product, and new constant-degree expanders and extractors. In 41st Annual Symposium on Foundations of Computer Science, FOCS 2000, 12-14 November 2000, Redondo Beach, California, USA, pages 3–13. IEEE Computer Society, 2000. doi:10.1109/SFCS.2000.892006.
  • [48] Alfréd Rényi. On the minimal number of terms of the square of a polynomial. Hungarica Acta Math, 1:30–34, 1947. URL: https://static.renyi.hu/renyi_cikkek/1947_On_the_minimal_number_of_terms_of_the_square_of_a_polynomial.pdf.
  • [49] Chandan Saha and Bhargav Thankey. Hitting sets for orbits of circuit classes and polynomial families. ACM Trans. Comput. Theory, 16(3):14:1–14:50, 2024. doi:10.1145/3665800.
  • [50] Ramprasad Saptharishi. A survey of lower bounds in arithmetic circuit complexity. GitHub repository, version 9.0.3, 2021. A community curated survey. URL: https://github.com/dasarpmar/lowerbounds-survey.
  • [51] Amir Shpilka and Avi Wigderson. Depth-3 arithmetic circuits over fields of characteristic zero. Comput. Complex., 10(1):1–27, 2001. doi:10.1007/PL00001609.
  • [52] Amit Sinhababu and Thomas Thierauf. Factorization of polynomials given by arithmetic branching programs. Comput. Complexity, 30(2):Paper No. 15, 47, 2021. doi:10.1007/s00037-021-00215-0.
  • [53] P. Stevenhagen and H. W. Lenstra. Chebotarëv and his density theorem. The Mathematical Intelligencer, 18(2):26–37, 1996. doi:10.1007/BF03027290.
  • [54] Terence Tao. An uncertainty principle for cyclic groups of prime order. 2004. doi:10.48550/arXiv.math/0308286.
  • [55] W. Verdenius. On the number of terms of the square and the cube of polynomials. Indag. Math., 11:459–465, 1949.
  • [56] Joachim von zur Gathen and Jürgen Gerhard. Modern Computer Algebra. Cambridge University Press, Cambridge, third edition, 2013. doi:10.1017/CBO9781139856065.
  • [57] Joachim von zur Gathen and Erich L. Kaltofen. Factoring sparse multivariate polynomials. J. Comput. Syst. Sci., 31(2):265–287, 1985. doi:10.1016/0022-0000(85)90044-3.