On vanishing of Kronecker coefficients

Christian Ikenmeyer, Ketan D. Mulmuley, Michael Walter

Introduction

One class of representation-theoretic obstructions in the context of the geometric complexity theory (GCT) approach to the permanent vs. determinant problem (Mulmuley & Sohoni, 2008; Mulmuley, 2011) is based on the existence of vanishing rectangular Kronecker coefficients (Mulmuley & Sohoni, 2008; Bürgisser et al., 2011a; Kumar, 2015). These are called occurrence-based obstructions, as opposed to the more general multiplicity-based obstructions. We refer to Bürgisser et al. (2011b); Bürgisser (2016) for introduction and background. It is now known that such occurrence-based obstructions based on vanishing of Kronecker coefficients cannot be used for proving superpolynomial lower bounds for the permanent (Ikenmeyer & Panova, 2016). However, they may still be useful for proving modest polynomial lower bounds. The partition triples associated with the rectangular Kronecker coefficients lie in the moment cone (Kirwan, 1984) associated with the Kronecker coefficients, called the Kronecker cone (Bürgisser et al., 2011a; Kumar, 2015). As pointed out in Kumar (2015), this makes the problem of showing the existence of such partition triples rather challenging, since the asymptotic techniques of algebraic geometry and representation theory, such as the ones based on the effective descriptions of the linear inequalities defining the Kronecker cone (Berenstein & Sjamaar, 2000; Klyachko, 2004; Ressayre, 2010; Vergne & Walter, 2017), cannot be used to prove this existence.

The main result in this article (1.5) establishes the existence of a superpolynomial number of partition triples with vanishing Kronecker coefficients, in the Kronecker cone for the given partition size, and satisfying a relaxed form of the additional shape restrictions that arise in GCT.

Its proof, based on the explicit proof strategy of GCT (Mulmuley, 2011, 2010a, 2010b), also yields results concerning the complexity of Kronecker coefficients that are of independent interest. The first such result (1.1) shows that the problem of deciding positivity of Kronecker coefficients is NP-hard. The second result (1.3) gives the first known instance of a positive (#P\#P) formula for a subclass of Kronecker coefficients whose positivity is NP-hard to decide.

We now state these results in more detail after the following preliminary section.

When we encode partitions as bit strings there are two fundamentally different ways of doing it: As a list of numbers in binary or as a list of numbers in unary. Note that in unary transposing a partition does not significantly change its encoding size, but in binary the 1-row partition (n)(n) can be encoded using O(log⁡n)O(\log n) bits, while (n)T=(1,1,…,1)(n)^{T}=(1,1,\ldots,1) requires O(n)O(n) bits. We will mostly encode partitions in unary.

It is natural to interpret partitions as vectors with integer entries, so that we have a well-defined addition and scalar multiplication with nonnegative integers. Moreover, dividing a partition by an integer results in a vector with rational entries.

Let G:=GL⁡rG:=\operatorname{GL}_{r} denote the general linear group, i.e., the group of invertible r×rr\times r matrices. Let VV be a finite dimensional vector space and let GL⁡(V)\operatorname{GL}(V) denote the set of linear isomorphisms of VV. A group homomorphism ϱ:G→GL⁡(V)\varrho:G\to\operatorname{GL}(V) is called a representation of GG. We say that VV is a representation if ϱ\varrho is clear from the context. We say that GG acts linearly on VV and use the short notation gv:=(ϱ(g))(v)gv:=(\varrho(g))(v) for g∈Gg\in G, v∈Vv\in V. If all the coordinate functions of ϱ\varrho are given by multivariate polynomials in the r2r^{2} coordinate variables of GL⁡r\operatorname{GL}_{r}, then we call ϱ\varrho a polynomial representation.

A linear subspace W⊆VW\subseteq V that satisfies ∀w∈W, g∈G:gw∈W\forall w\in W,\,g\in G:gw\in W, is called a subrepresentation. Subrepresentations of polynomial representations are always polynomial. For every representation VV, the zero vector space and VV itself are two subrepresentations. If VV has only these two subrepresentations, then VV is called irreducible. Given two representations (V,ρV)(V,\rho_{V}) and (W,ρW)(W,\rho_{W}), then a linear map φ:V→W\varphi:V\to W is called equivariant if gφ(v)=φ(gv)g\varphi(v)=\varphi(gv) for all g∈Gg\in G, v∈Vv\in V. If φ\varphi is an equivariant isomorphism of vector spaces, then φ\varphi is called a GG-isomorphism and the representations VV and WW are called isomorphic. The different types of isomorphic irreducible polynomial representations of GG have been classified completely: They are indexed by partitions of height at most rr. In a representation VV the sum of all subrepresentations of type λ\lambda is called the λ\lambda-isotypic component. Every representation decomposes into a direct sum of isotypic components.

The multiplicity of the type λ\lambda in a representation VV is the dimension of the vector space of highest weight vectors of type λ\lambda. If we decompose VV into a direct sum of irreducibles, then this multiplicity counts how often a copy of type λ\lambda appears in the decomposition.

If we have kk commuting actions of several copies of GG on VV, then we use the representation theory of the cartesian powers GkG^{k}, which is very similar to the representation theory of GG. We will mainly be concerned with k=2k=2 and k=3k=3. The types of irreducible representations of GkG^{k} are given by kk-tuples of partitions, weight vectors are defined by their scaling behavior under kk-tuples of diagonal matrices, and the Lie algebra action is defined via kk-tuples of matrices, where the raising operators are kk-tuples of matrices in which only one matrix is nonzero. Irreducible representations of GG are called Weyl modules, while irreducible representations of GkG^{k} are isomorphic to a kk-fold tensor product of Weyl modules.

The Kronecker coefficient arises as a multiplicity in several other representation theoretic decompositions, for example as the multiplicity of Vλ(GL⁡r)⊗Vμ(GL⁡r)⊗Vπ(GL⁡r)V_{\lambda}(\operatorname{GL}_{r})\otimes V_{\mu}(\operatorname{GL}_{r})\otimes V_{\pi}(\operatorname{GL}_{r}) in V(n)(GL⁡r3)V_{(n)}(\operatorname{GL}_{r^{3}}) via the group homomorphism GL⁡r3→GL⁡r3\operatorname{GL}_{r}^{3}\to\operatorname{GL}_{r^{3}}, (g,g′,g′′)↦g⊗g′⊗g′′(g,g^{\prime},g^{\prime\prime})\mapsto g\otimes g^{\prime}\otimes g^{\prime\prime}, where λ\lambda, μ\mu, and π\pi have exactly nn boxes.

For kμ,πλk^{\lambda}_{\mu,\pi} to be positive it is required that λ\lambda, μ\mu, and π\pi are partitions of the same number, i.e., ∣λ∣=∣μ∣=∣π∣=n|\lambda|=|\mu|=|\pi|=n for some nn. This implies that if kμ,πλ>0k^{\lambda}_{\mu,\pi}>0, then the rescaled partitions λ/n\lambda/n, μ/n\mu/n, and π/n\pi/n are three discrete probability distributions. Another necessary condition for kμ,πλ>0k^{\lambda}_{\mu,\pi}>0 is ht⁡(λ)≤ht⁡(μ)⋅ht⁡(π)\operatorname{ht}(\lambda)\leq\operatorname{ht}(\mu)\cdot\operatorname{ht}(\pi). The coefficient is invariant under permuting the three parameters, so

Moreover, transposing any two of the three parameters does not change the coefficient:

Given two partition triples (λ,μ,π)(\lambda,\mu,\pi) and (λ′,μ′,π′)(\lambda^{\prime},\mu^{\prime},\pi^{\prime}) such that kμ,πλ>0k^{\lambda}_{\mu,\pi}>0 and kμ′,π′λ′>0k^{\lambda^{\prime}}_{\mu^{\prime},\pi^{\prime}}>0, then kμ+μ′,π+π′λ+λ′>0k^{\lambda+\lambda^{\prime}}_{\mu+\mu^{\prime},\pi+\pi^{\prime}}>0. This is called the semigroup property. The convex cone defined by

Finding a combinatorial description of kμ,πλk^{\lambda}_{\mu,\pi} is an important outstanding problem (see 1.3 below). Only for some special cases is a combinatorial description known, for example for the Littlewood-Richardson coefficients. The Littlewood-Richardson coefficients are those kμ,πλk^{\lambda}_{\mu,\pi} for which λ\lambda, μ\mu, and π\pi have a sufficiently long first row such that ∣λ‾∣=∣μ‾∣+∣π‾∣|\overline{\lambda}|=|\overline{\mu}|+|\overline{\pi}|, where λ‾\overline{\lambda} is the partition λ\lambda with its longest row removed. The positivity of Littlewood-Richardson coefficients can be decided in strongly polynomial time (Knutson & Tao, 2001; Mulmuley et al., 2011).

Another important subclass of Kronecker coefficients are the rectangular Kronecker coefficients. For a given partition λ\lambda with size ∣λ∣|\lambda| divisible by rr, let δ(λ)\delta(\lambda) denote the rectangular partition (d,…,d)(d,\ldots,d) (rr times), where d=∣λ∣/rd=|\lambda|/r. We call the Kronecker coefficient kδ(λ),δ(λ)λk^{\lambda}_{\delta(\lambda),\delta(\lambda)} rectangular. Rectangular Kronecker coefficients play a special role in geometric complexity theory (see 1.4 below).

2 NP-hardness of deciding positivity of Kronecker coefficients

Let Kronecker be the problem of deciding positivity of kμ,πλk_{\mu,\pi}^{\lambda}, given as input the three partitions λ\lambda, μ\mu, and π\pi in unary. This problem is of fundamental interest in the context of the explicit proof strategy of GCT (Mulmuley, 2011, 2010a, 2010b). Our first result is the following:

It was conjectured in Mulmuley (2010b) that the problem of deciding positivity of Kronecker coefficients is in PP. 1.1 shows that this is not so, in general, assuming that P≠NPP\not=NP. This is in contrast to the special case of the Littlewood-Richardson coefficients, where positivity can be decided in strongly polynomial time, as explained above.

3 A #​P#𝑃\#P-formula for a subclass of partitions of type NP

To find a positive formula for Kronecker coefficients “akin to” the well known positive Littlewood-Richardson rule is an unsolved problem in classical representation theory. We refer to Stanley (2002) for the history and importance of this problem, where it is listed as one of the twenty-five “outstanding open problems”. In classical representation theory, the phrase “akin to” is used only informally. A formal complexity-theoretic version of this problem is to find a #P\#P-formula for Kronecker coefficients. By a #P\#P-formula for the Kronecker coefficient kμ,πλk^{\lambda}_{\mu,\pi}, we mean a formula of the form:

where, for a partition λ=(λ1,λ2,…,λl)\lambda=(\lambda_{1},\lambda_{2},\ldots,\lambda_{l}), ⟨λ⟩\langle\lambda\rangle denotes the total bitlength of the specification of λj\lambda_{j}’s in binary, p(⟨λ⟩,⟨μ⟩,⟨π⟩)p(\langle\lambda\rangle,\langle\mu\rangle,\langle\pi\rangle) is a polynomially-bounded function of the bit-lengths ⟨λ⟩,⟨μ⟩\langle\lambda\rangle,\langle\mu\rangle, and ⟨π⟩\langle\pi\rangle, and F(λ,μ,π,σ)F(\lambda,\mu,\pi,\sigma) is a polynomial-time-computable -11 function of λ,μ,π\lambda,\mu,\pi, and the bit-string σ\sigma. By a positive formula, we mean a #P\#P-formula henceforth.

Let Π\Pi be a class of partition triples. We say that Π\Pi is of type NP if the problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi}, with (λ,μ,π)∈Π(\lambda,\mu,\pi)\in\Pi, is NP-hard. (The problem mentioned here is a promise problem. That is, we are promised that the input triple is in the subclass Π\Pi.) Likewise, we say that Π\Pi is of type P if the problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi}, with (λ,μ,π)∈Π(\lambda,\mu,\pi)\in\Pi, is in P.

All positive rules known so far for restricted classes of Kronecker coefficients have been for subclasses of partition triples that are either known or conjectured to be of type P. For example, the classical Littlewood-Richardson rule gives a positive rule for Littlewood-Richardson coefficients, which, as already mentioned, constitute a special class of Kronecker coefficients. The corresponding subclass of partition triples is of type P, since the problem of deciding positivity of Littlewood-Richardson coefficients is in PP (Knutson & Tao, 2001; Mulmuley et al., 2011). Blasiak et al. (2015) give a positive rule for Kronecker coefficients when two of the partitions have height at most two. The corresponding subclass of partition triples is of type P, since the Kronecker coefficient can be computed in this case (and more generally, for partitions of bounded height) in polynomial time (Christandl et al. (2012); Baldoni et al. (2017)). Blasiak (2017) gives a positive rule for Kronecker coefficients when one of the partitions is a hook. The corresponding subclass of partition triples is conjectured to be of type P, since the problem of deciding positivity of Kronecker coefficients, when one of the partitions is a hook, is believed to be in PP (in view of 6.10).

The following result gives the first known instance of a positive rule for Kronecker coefficients for a subclass of partition triples of type NP.

There exists a #P\#P-formula for Kronecker coefficients for a subclass of partition triples of type NP. Here the partition triples can be specified in unary or binary.

The proof of this result exhibits an explicit such subclass of partition triples of type NP (see 2 and 3).

1.3 provides good evidence in support of the conjecture in Mulmuley (2010b) that there exists a #P\#P-formula for Kronecker coefficients in general. This would in particular imply that Kronecker is in NP, which is not known so far.

4 Exceptional Kronecker coefficients

In order for a Kronecker coefficient kμ,πλk_{\mu,\pi}^{\lambda} to be useful for proving a polynomial lower bound for the permanent, the partition triple must have a number of exceptional properties (Mulmuley & Sohoni, 2008; Bürgisser et al., 2011b). This is captured by the following definition:

Fix any constant 0<ϵ≤10<\epsilon\leq 1, and a constant b>1b>1. We call a partition triple (λ,μ,π)(\lambda,\mu,\pi) with ∣λ∣=∣μ∣=∣π∣\lvert\lambda\rvert=\lvert\mu\rvert=\lvert\pi\rvert (ϵ,b)(\epsilon,b)-exceptional if:

μ=π=δ(λ)\mu=\pi=\delta(\lambda), with ∣λ∣=∣μ∣=∣π∣|\lambda|=|\mu|=|\pi| divisible by r:=ht⁡(μ)=ht⁡(π)r:=\operatorname{ht}(\mu)=\operatorname{ht}(\pi),

ht⁡(λ)≤rϵ\operatorname{ht}(\lambda)\leq r^{\epsilon},

(λ,μ,π)∈Kron⁡(r)(\lambda,\mu,\pi)\in\operatorname{Kron}(r),

λ0≥∣λ∣(1−rϵ/2−1)\lambda_{0}\geq|\lambda|(1-r^{\epsilon/2-1}).

We also call a partition tuple merely exceptional, without mentioning ϵ\epsilon and bb, if it is understood that ϵ\epsilon can be chosen to be arbitrarily small, with bb a large enough constant depending on ϵ\epsilon, and r→∞r\rightarrow\infty.

By Bürgisser et al. (2011a), 2 implies 4, assuming that the height of λ\lambda is ≤r2\leq r^{2}, which is so by 3.

The constraint 4 is significant. Proving existence of the partition triples as in 1.4 is delicate because of this constraint. Indeed, it may be possible to prove existence of superpolynomially many partition triples satisfying the constraints other than 2 and 4 using the known linear inequalities defining the Kronecker cone (Berenstein & Sjamaar, 2000; Klyachko, 2004; Ressayre, 2010; Vergne & Walter, 2017). But the constraint 4 implies that such asymptotic techniques based on the description of the Kronecker cone cannot be used to demonstrate existence of partition triples as in 1.4. This is the main significance of the results in Bürgisser et al. (2011a); Kumar (2015).

By the Saturation Theorem (Knutson & Tao, 1999; Derksen & Weyman, 2000), the Littlewood-Richardson coefficients cannot vanish for the partition triples that lie in the analogously defined Littlewood-Richardson cone. The constraint 4 also implies that in order to prove existence of the partition triples as in 1.4, one needs to understand the failure of the saturation property for the Kronecker coefficients in one way or another.

The constraint 7 is motivated by Kadish & Landsberg (2014). There, it is shown that this condition holds if Vλ(G)V_{\lambda}(G) is a representation-theoretic obstruction (Mulmuley & Sohoni, 2008).

It is a priori not at all clear that for any given constant 0<ϵ≤10<\epsilon\leq 1 and a large enough constant b>1b>1 depending on ϵ\epsilon, exceptional partition triples exist for arbitrary rr. The experimental evidence in Ikenmeyer (2012) for small values of rr (with suitable ϵ\epsilon and bb) suggests that they are very rare, though they do exist for these small values. In summary, although their density can be expected to be extremely small, it is a relevant and rather non-trivial problem in the context of GCT to show that exceptional partition triples exist and that their number is large enough.

5 Construction of superpolynomially many partition triples in the Kronecker cone with vanishing Kronecker coefficients

As the first step towards this goal, we relax the condition 2 to the weaker requirement that only μ=π\mu=\pi, the condition 6 to the weaker requirement weaker that λ\lambda is not a hook (since it can be shown that p(λ)=0p(\lambda)=0 if λ\lambda is a hook), and ignore the condition 7. A priori, it is not clear that partition triples with these properties exist even after this shape relaxation, since condition 4 is retained. The following result shows that the number of Kronecker coefficients with this relaxation of 1.4 is superpolynomial.

For any 0<ϵ≤10<\epsilon\leq 1, there exists 0<a<10<a<1, such that, for all mm, there exist Ω(2ma)\Omega(2^{m^{a}}) partition triples (λ,μ,π)(\lambda,\mu,\pi) such that

ht⁡(μ)=m\operatorname{ht}(\mu)=m, and ht⁡(λ)≤mϵ\operatorname{ht}(\lambda)\leq m^{\epsilon},

(λ,μ,π)∈Kron⁡(m)(\lambda,\mu,\pi)\in\operatorname{Kron}(m),

Assuming coNP ≠\neq NP, the set of partition triples satisfying constraints 1–\short6 as well as

(λT,μT,μ)∈Kron⁡(m′)(\lambda^{T},\mu^{T},\mu)\in\operatorname{Kron}(m^{\prime}), where m′m^{\prime} is the maximum of the heights of λT\lambda^{T}, μT\mu^{T} and μ\mu, and

(λ,μT,μT)∈Kron⁡(m′′)(\lambda,\mu^{T},\mu^{T})\in\operatorname{Kron}(m^{\prime\prime}), where m′′m^{\prime\prime} is the maximum of the heights of λ\lambda and μT\mu^{T},

is superpolynomial in mm, as m→∞m\rightarrow\infty.

The constraints \short4, 7, and 8 together guarantee that the vanishing of kμ,μλk^{\lambda}_{\mu,\mu} cannot be directly shown using the defining inequalities of the Kronecker cone (Berenstein & Sjamaar, 2000; Klyachko, 2004; Ressayre, 2010; Vergne & Walter, 2017) in conjunction with the known symmetries 1 and 2.

Our proof of 1.5 shows that the partition triples (λ,μ,μ)(\lambda,\mu,\mu) satisfying the constraints therein can even be constructed explicitly. This means there is a one-to-one map from the set of Boolean strings of length ≤ma\leq m^{a} to the set of partitions triples (λ,μ,μ)(\lambda,\mu,\mu) with properties \short1–\short6 that can be computed in poly⁡(m)\operatorname{poly}(m) time.

While 1.5 shows existence of superpolynomially many partition triples satisfying the constraints therein, the density of such partition triples is exponentially small, since aa therein is much smaller than 11 (see 5.6). This may explain why vanishing Kronecker coefficients with the partition triples in the Kronecker cone occur so rarely in computer experiments, as observed in Ikenmeyer (2012).

6 Proof technique

1.1 is proved by extending the NP-completeness technique in Brunetti et al. (2001) in conjunction with the fundamental lower and upper bounds on Kronecker coefficients established in Manivel (1997); Bürgisser & Ikenmeyer (2013); Vallejo (2000). 1.3 is a byproduct of this proof.

A refined form of 1.1 lies at the heart of the proof of 1.5. Specifically, we show in 4.2 that the problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi} remains NP-hard under polynomial-time many-one reductions (Karp, 1972) even when the partitions (λ,μ,π)(\lambda,\mu,\pi) are required to satisfy the constraints 2–\short6. This is done by extending the proof technique of 1.5 using the result in Bürgisser et al. (2011a) that (λ,δ(λ),δ(λ))(\lambda,\delta(\lambda),\delta(\lambda)), for ∣λ∣|\lambda| divisible by rr, lies in the Kronecker cone whenever the height of λ\lambda is ≤r2\leq r^{2}. By Fortune (1979), if there exists a co-sparse NP-complete language under polynomial-time many-one-reductions, then P=NP. (Here we call a language sparse if the number of strings in it of bitlength ≤N\leq N is bounded by a fixed polynomial in NN. It is called co-sparse if its complement is sparse.) Hence 4.2 in conjunction with Fortune (1979) implies that the set of partition triples satisfying the constraints 1–\short6 is non-sparse, i.e., has size superpolynomial in mm, assuming that P≠\neqNP.

To prove 1.5, we have to discard of the assumption that P=NP and replace the superpolynomial bound by Ω(2ma)\Omega(2^{m^{a}}) bound for some a>0a>0. This is done in 5.3 by exhibiting a polynomial-time one-one reduction from the 3D Matching problem (Garey & Johnson, 1979) to the problem of deciding positivity of Kronecker coefficients, with the partition triples satisfying the constraints 2–\short6, where a polynomial-time one-one reduction means an injective polynomial-time many-one reduction. Hence, the Ω(2ma)\Omega(2^{m^{a}}) bound in 1.5 follows from a similar lower bound on the number of instances of the 3D Matching problem with a “NO” answer. The proof automatically shows that Ω(2ma)\Omega(2^{m^{a}}) partition triples satisfying the constraints in 1.5 can be constructed explicitly. This follows by fixing a suitable set of 2Nb2^{N^{b}} instances, for some constant b>0b>0, of the 3D Matching problem of bitlength ≤N\leq N with “NO” answer, and mapping them injectively, via a sequence of polynomial time one-one reductions, to Ω(2ma)\Omega(2^{m^{a}}) such partition triples.

1.6 is proved by extending the proof of 1.5 using an auxiliary result, 5.7, which extends the hardness vs. non-sparseness result in Fortune (1979), together with the result in Bürgisser et al. (2017) which asserts that the membership problem for the Kronecker cone is in \mboxNP∩\mboxcoNP\mbox{NP}\cap\mbox{coNP}.

7 Effectiveness of the explicit proof strategy

Perhaps the most novel aspect of this paper is the synthesis of the representation theory of Kronecker coefficients with the theory of NP-completeness to prove unconditionally existence of superpolynomially many partition triples in the Kronecker cone with vanishing Kronecker coefficients.

In principle, the existence of partition triples satisfying the constraints in 1.5 may be proved by a nonconstructive technique. Yet, the only way we can prove this existence at present is by constructing such partitions explicitly, using the theory of algorithms, as done in the proof of 1.5. Thus this proof illustrates effectiveness of the explicit proof strategy of GCT (Mulmuley, 2011, 2010a, 2010b) in a nontrivial setting.

8 Organization

The rest of this article is organized as follows. 2 describes the lower and upper bounds for the Kronecker coefficients that are needed for the proofs of 1.1 and 1.3. These proofs are given in 3. A refinement of 1.1, which is needed for the proof of 1.5, is proved in 4. 1.5 and 1.6 are proved in 5. 6 proves additional results in support of the conjecture in Mulmuley (2010b) that the problem of deciding positivity of rectangular Kronecker coefficients is in PP.

Lower and upper bounds for the Kronecker coefficient

In this section, we give representation-theoretic proofs of some known lower and upper bounds from Manivel (1997); Bürgisser & Ikenmeyer (2013); Vallejo (2000) for the Kronecker coefficients. These bounds as well as their representation-theoretic interpretation given here will play a crucial role in the proofs of 1.1 and 1.3 in 3. We begin with the following well-known result (whose proof we include for the sake of completeness):

Since [λT]=[λ]⊗[(1n)][\lambda^{T}]=[\lambda]\otimes[(1^{n})], [(1n)]⊗[(1n)]=[(n)][(1^{n})]\otimes[(1^{n})]=[(n)], and [(n)][(n)] is the trivial representation, 3 follows at once.

Given a (finite) point set P⊆{0,…,r−1}3P\subseteq\{0,\dots,r-1\}^{3}, let xP(i)x_{P}(i), 0≤i≤r−10\leq i\leq r-1, be the number of points in PP with the xx-coordinate ii. We call xP=(xP(0),…,xP(r−1))x_{P}=(x_{P}(0),\ldots,x_{P}(r-1)) the xx-marginal of PP. We similarly define the yy-marginal yPy_{P} and the zz-marginal zPz_{P}. The triple (xP,yP,zP)(x_{P},y_{P},z_{P}) is called the marginals of PP.

We define tμ,πλt^{\lambda}_{\mu,\pi} as the number of point sets P⊆{0,…,r−1}3P\subseteq\{0,\dots,r-1\}^{3} with marginals (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}).

Note that λiT\lambda^{T}_{i} is the number of boxes in the ii-th column of λ\lambda.

The coefficients tμ,πλt^{\lambda}_{\mu,\pi} have a pleasant representation-theoretical interpretation that is closely related to 2.1. To see this, observe that we can associate with any point set

Following Vallejo (2000), we call a subset P⊆{0,…,r−1}3P\subseteq\{0,\dots,r-1\}^{3} a pyramid if, for any (x,y,z)∈P(x,y,z)\in P and 0≤x′≤x0\leq x^{\prime}\leq x, 0≤y′≤y0\leq y^{\prime}\leq y, 0≤z′≤z0\leq z^{\prime}\leq z, we have that (x′,y′,z′)∈P(x^{\prime},y^{\prime},z^{\prime})\in P. (It would also be natural to call such PP a 3-partition; cf. Manivel (1997).)

Let pμ,πλp^{\lambda}_{\mu,\pi} denote the number of pyramids with marginals (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}).

From our representation-theoretic interpretation, we directly obtain the following fundamental bounds, which were proved previously using different methods in Manivel (1997); Bürgisser & Ikenmeyer (2013) (cf. Vallejo, 2000):

For all partitions λ,μ,π\lambda,\mu,\pi, we have pμ,πλ≤kμ,πλ≤tμ,πλp^{\lambda}_{\mu,\pi}\leq k^{\lambda}_{\mu,\pi}\leq t^{\lambda}_{\mu,\pi}.

is a pyramid with marginals (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}). We will show that ψP\psi_{P} is not only a weight vector, but in fact a highest weight vector. For this, we need to argue that ψ\psi is annihilated by all raising operators (cf. 1.1). Thus consider (Ex′,x,0,0)(E_{x^{\prime},x},0,0), where Ex′,xE_{x^{\prime},x} denotes the upper triangular matrix with a single 1 in the x′x^{\prime}-th row and xx-th column, and otherwise zero (here x′<xx^{\prime}<x). Its action on ψP\psi_{P} is given by

since each summand vanishes individually. Indeed, if x≠xjx\neq x_{j} then δx,xj=0\delta_{x,x_{j}}=0 and so the summand is zero. Otherwise, if x=xjx=x_{j} then (xj,yj,zj)∈P(x_{j},y_{j},z_{j})\in P and x′<xx^{\prime}<x imply that (x′,yj,zj)∈P(x^{\prime},y_{j},z_{j})\in P by the pyramid condition; therefore ex′⊗eyj⊗ezje_{x^{\prime}}\otimes e_{y_{j}}\otimes e_{z_{j}} appears twice in the wedge product and so the summand vanishes as well. The same argument applies to the other generators (0,Ey′,y,0)(0,E_{y^{\prime},y},0) and (0,0,Ez′,z)(0,0,E_{z^{\prime},z}) of n\mathfrak{n}. Thus we conclude that the pyramid condition ensures that ψP\psi_{P} is a highest weight vector.

Let λ\lambda, μ\mu, π\pi be partitions such that any point set with marginals (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}) is necessarily a pyramid. Then kμ,πλ=tμ,πλ=pμ,πλk^{\lambda}_{\mu,\pi}=t^{\lambda}_{\mu,\pi}=p^{\lambda}_{\mu,\pi}.

Kronecker coefficients with #P-formulae

For this, we first derive a sufficient condition on the marginals (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}) such that any compatible point set is necessarily a pyramid (and hence 2.8 is applicable). Adapting the approach of Brunetti et al. (2001), we consider a point set PP such that Pr⊆P⊊Pr+1P_{r}\subseteq P\subsetneq P_{r+1}, where

denotes the simplex of side length r≥1r\geq 1. Let nn denote the total number of points in PP. Then the projection of the barycenter bP:=∑p∈Ppb_{P}:=\sum_{p\in P}p of PP onto the diagonal (1,1,1)(1,1,1) can be computed as follows:

where brb_{r} denotes the barycenter of the simplex PrP_{r}. Note that this formula depends only nn, the number of points in the point set PP. We can thus define a function p(n)p(n) by 4, first for all nn such that ∣Pr∣≤n<∣Pr+1∣\lvert P_{r}\rvert\leq n<\lvert P_{r+1}\rvert, and then, by varying rr, for all nn. Explicitly,

where r(n)r(n) is the maximal rr such that ∣Pr∣=r(r+1)(r+2)/6≤n\lvert P_{r}\rvert=r(r+1)(r+2)/6\leq n.

Let us call (λ,μ,π)(\lambda,\mu,\pi), with ∣λ∣=∣μ∣=∣π∣=n≠0|\lambda|=|\mu|=|\pi|=n\not=0, simplex-like if there exists some rr such that the Young diagrams of λ,μ\lambda,\mu, and π\pi have at most r+1r+1 columns, and

Whether (λ,μ,π)(\lambda,\mu,\pi) is simplex-like can be checked in polynomial time (even assuming that λ,μ\lambda,\mu and π\pi are given in binary).

The following lemma justifies the term “simplex-like”.

Let (λ,μ,π)(\lambda,\mu,\pi) be simplex-like. Then any point set PP with marginals (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}) is necessarily of the form Pr⊆P⊊Pr+1P_{r}\subseteq P\subsetneq P_{r+1}, for some r≥1r\geq 1. In particular, PP is a pyramid.

The last step follows because (λ,μ,π)(\lambda,\mu,\pi) is simplex-like.

Let (λ,μ,π)(\lambda,\mu,\pi) be simplex-like. Then kμ,πλ=tμ,πλ=pμ,πλk^{\lambda}_{\mu,\pi}=t^{\lambda}_{\mu,\pi}=p^{\lambda}_{\mu,\pi}. In particular, this family of Kronecker coefficients has a #P\#P-formula. Here, λ,μ\lambda,\mu and π\pi can be given in unary or binary.

This follows from 3.2, 2.8, and the fact that tμ,πλt^{\lambda}_{\mu,\pi} has a #P\#P-formula. The last assertion follows because the bit-length of the unary specification of a simplex-like partition triple is polynomial in the bit-length of its binary specification.

One important class of simplex-like marginals is the following. Let (φT,φT,φT)(\varphi^{T},\varphi^{T},\varphi^{T}) denote the marginals of the simplex P2rP_{2r}, where r≥1r\geq 1. Define

Since ∣P2r∣≤n=∣P2r∣+(r+1)<∣P2r+1∣\lvert P_{2r}\rvert\leq n=\lvert P_{2r}\rvert+(r+1)<\lvert P_{2r+1}\rvert, this is indeed equal to p(n)p(n), where nn is the number of boxes of each of λ\lambda, μ\mu, and π\pi (compare with 5). These marginals arise when embedding permutation matrices on top of the simplex P2rP_{2r}, and in Brunetti et al. (2001) it was shown using this construction that:

The problem of deciding positivity of tμ,πλt^{\lambda}_{\mu,\pi}, given λ,μ,π\lambda,\mu,\pi in unary, is NP-hard with respect to polynomial-time many-one reductions, even when (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}) is restricted to be of the form 6.

Thus it follows at once from 3.4, in conjunction with this result, that:

The problem of deciding positivity of the Kronecker coefficient kμ,πλk^{\lambda}_{\mu,\pi}, given λ,μ,π\lambda,\mu,\pi in unary, is NP-hard with respect to polynomial-time many-one reductions, even when (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}) is restricted to be of the form 6.

This proves 1.1. 1.3 follows from this result and 3.4.

Refined NP-hardness result

Let Restricted Kronecker be the problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi}, when λ,μ,π\lambda,\mu,\pi satisfy the constraints 2–\short6. Specifically:

ht⁡(λ)≤mϵ\operatorname{ht}(\lambda)\leq m^{\epsilon},

(λ,μ,π)(\lambda,\mu,\pi) is in the Kronecker cone Kron⁡(m)\operatorname{Kron}(m),

For the proof of 1.5, we need the following refinement of 3.7:

Restricted Kronecker is NP-hard with respect to polynomial-time many-one reductions.

Let λ\lambda and μ\mu be partitions so that ht⁡(λ)≤h2\operatorname{ht}(\lambda)\leq h^{2}, where hh is the height of the smallest column of μ\mu. Then (λ,μ,μ)(\lambda,\mu,\mu) is in the Kronecker cone Kron⁡(l)\operatorname{Kron}(l), where l=max⁡{ht⁡(λ),ht⁡(μ)}l=\max\{\operatorname{ht}(\lambda),\operatorname{ht}(\mu)\}.

It is shown in Bürgisser et al. (2011a) that (λ,δ,δ)(\lambda,\delta,\delta) is in the Kronecker cone whenever δ\delta is a rectangle of height at least hh, and ht⁡(λ)≤h2\operatorname{ht}(\lambda)\leq h^{2}. As the Kronecker cone is a cone, this is also true if we rescale each of λ\lambda and δ\delta by an arbitrary positive number.

Let us write μ\mu as a sum of rectangles δ(1)+⋯+δ(k)\delta^{(1)}+\dots+\delta^{(k)}, where our assumption implies that each δ(j)\delta^{(j)} has height at least hh. It is easy to see that λ\lambda can be written as a sum λ(1)+⋯+λ(k)\lambda^{(1)}+\dots+\lambda^{(k)}, where each λ(j)\lambda^{(j)} is a rational partition with the same size as δ(j)\delta^{(j)} (i.e., ∣λ(j)∣=∣δ(j)∣|\lambda^{(j)}|=|\delta^{(j)}|), and with no more than h2h^{2} rows. By the preceding argument, each (λ(j),δ(j),δ(j))(\lambda^{(j)},\delta^{(j)},\delta^{(j)}) is in the Kronecker cone Kron⁡(l)\operatorname{Kron}(l). As cones are closed under addition, (λ,μ,μ)(\lambda,\mu,\mu) is likewise in the Kronecker cone.

Next, we generalize 3.4 to a larger class of marginals. Let (λ,μ,π)(\lambda,\mu,\pi) be simplex-like, and let PP be a corresponding point set with marginals (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}), so that Pr⊆P⊊Pr+1P_{r}\subseteq P\subsetneq P_{r+1} for some rr (3.2). Let QQ denote the following point set obtained by adjoining PP to a rectangular box of size a×b×ca\times b\times c, where b,c≥r+1b,c\geq r+1:

Then we have the following generalization of 3.2.

We apply the pedestal construction to marginals of the form 6. Let us choose a rectangular box of size c×s×sc\times s\times s, where s=2r+1s=2r+1. That is, we set

The problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi}, with (λ,μ,π)(\lambda,\mu,\pi) restricted as above, is NP-hard, since by 4.8, these Kronecker coefficients agree with the ones in 3.7, and we can transform instances of the latter to instances of the former in polynomial time (as ϵ\epsilon is fixed).

We now verify that the five constraints in 4.1 are all satisfied. The first is clearly satisfied. For the second,

by our choice of c=c(s)c=c(s). The third follows from 4.3, as ht⁡(λ)=s2\operatorname{ht}(\lambda)=s^{2}, while every column in μ\mu is of height at least cs≥scs\geq s. The fourth follows, since

assuming that rr is large enough. Finally, it is clear that λ\lambda is not a hook.

Construction of vanishing Kronecker coefficients with partition triples in the Kronecker cone

By Fortune (1979), if there exists a co-sparse NP-complete language under polynomial-time many-one-reductions, then P=NP. 4.2, in conjunction with this result, implies that, assuming P≠\neqNP, the set of partitions triples satisfying the constraints 1–\short6 is non-sparse, i.e, its cardinality is superpolynomial in mm. (The result in Fortune (1979) applies to NP-complete sets, rather than NP-hard sets. But we can still apply this result to the NP-complete set of simplex-like partition triples (cf. 3.7) with positive Kronecker coefficients to get the desired conclusion.) To prove 1.5, we have to get rid of the assumption that P≠\neqNP and replace the superpolynomial bound by Ω(2ma)\Omega(2^{m^{a}}) bound, for some positive constant aa. This will be achieved by 5.1 and 5.3 below.

Recall from Garey & Johnson (1979) that the 3D Matching problem is to decide, given a set M⊆W×X×YM\subseteq W\times X\times Y, where WW, XX, YY are disjoint sets of size qq, whether MM contains a (perfect) matching, i.e., a subset M′⊆MM^{\prime}\subseteq M of size qq such that no two elements of M′M^{\prime} agree in any coordinate. Without loss of generality, we assume henceforth that each element in W∪X∪YW\cup X\cup Y appears in some triple of MM. We denote instances of 3D Matching by tuples (M,W,X,Y,q)(M,W,X,Y,q). It is known that the 3D Matching problem is NP-complete (Garey & Johnson, 1979).

The number of instances (M,W,X,Y,q)(M,W,X,Y,q) of the 3D Matching problem with total bit-length ≤N\leq N such that MM does not have a matching is Ω(2Nb)\Omega(2^{N^{b}}) for some positive constant b<1b<1.

Furthermore, such instances can be constructed explicitly. That is, for some positive constant b<1b<1, there is a polynomial-time-computable one-to-one function that maps any pair of the form (N,σ)(N,\sigma), where NN is a positive integer and σ\sigma is a binary string of length ≤Nb\leq N^{b}, to an instance of 3D Matching problem without matching of bitlength ≤N\leq N.

Consider any fixed instance (M0,W0,X0,Y0,q0)(M_{0},W_{0},X_{0},Y_{0},q_{0}) of 3D Matching, such that M0M_{0} does not have a matching. Its bitlength is thus a constant. Given any instance (M,W,X,Y,q)(M,W,X,Y,q) of 3D Matching, with W,XW,X and YY disjoint from W0,X0W_{0},X_{0} and Y0Y_{0}, consider the padded instance (M∪M0,W∪W0,X∪X0,Y∪Y0,q+q0)(M\cup M_{0},W\cup W_{0},X\cup X_{0},Y\cup Y_{0},q+q_{0}). Clearly, M∪M0M\cup M_{0} also does not have a matching. The number of instances of the form (M∪M0,W∪W0,X∪X0,Y∪Y0,q+q0)(M\cup M_{0},W\cup W_{0},X\cup X_{0},Y\cup Y_{0},q+q_{0}) with bitlength ≤N\leq N is clearly Ω(2Nb)\Omega(2^{N^{b}}) for some positive constant b<1b<1. This is because, for a given qq, the total number of instances of the form (M,W,X,Y,q)(M,W,X,Y,q) is 2q32^{q^{3}}, and the bit-length of the specification of any instance of this form is O(q3)O(q^{3}). (We assume that MM is specified by its q×q×qq\times q\times q adjacency matrix.) Furthermore, it is easy to show that the padded instances of the form (M∪M0,W∪W0,X∪X0,Y∪Y0,q+q0)(M\cup M_{0},W\cup W_{0},X\cup X_{0},Y\cup Y_{0},q+q_{0}) can be constructed explicitly.

Recall (cf. 1.6) that a polynomial-time one-one reduction means an injective polynomial-time many-one reduction.

There exists a polynomial-time one-one reduction ϕ\phi from the set of instances (M,W,X,Y,q)(M,W,X,Y,q) of 3D Matching of total bit-length nn to the set of partition triples (λ,μ,π)(\lambda,\mu,\pi) satisfying the conditions 2–\short6 (i.e., instances of Restricted Kronecker), with m=poly⁡(n)m=\operatorname{poly}(n), such that MM contains a matching iff the Kronecker coefficient associated with the partition triple ϕ(E)\phi(E) is positive.

Since 3D Matching is in NP, it follows from 4.2 that there exists a polynomial-time many-one reduction ϕ\phi from 3D Matching to the Restricted Kronecker problem of deciding positivity of the Kronecker coefficient kμ,πλk^{\lambda}_{\mu,\pi}, with (λ,μ,π)(\lambda,\mu,\pi) satisfying the constraints 2–\short6. We have to show that this reduction ϕ\phi can be chosen to be injective. We can obtain such an injective reduction ϕ\phi by composing the following sequence of polynomial-time computable one-one-reductions (1):

From 3D Matching to 4-Partition (cf. Theorem 4.3 in Garey & Johnson, 1979):

The 4-Partition problem is to decide, given a set AA of size 4m4m, a positive integer bound BB, a positive integer size s(a)s(a) for each a∈Aa\in A such that B/5<s(a)<B/3B/5<s(a)<B/3 and ∑a∈As(a)=mB\sum_{a\in A}s(a)=mB, whether AA can be partitioned into mm disjoint subsets A1,…,AmA_{1},\ldots,A_{m}, each of size four, such that, for each 1≤i≤m1\leq i\leq m, ∑a∈Ais(a)=B\sum_{a\in A_{i}}s(a)=B. We denote such an instance of 4-Partition by the tuple (A,m,B,s)(A,m,B,s).

The reduction in Garey & Johnson (1979) maps a given instance (M,W,X,Y,q)(M,W,X,Y,q) of 3D Matching to an instance (A,m,B,s)(A,m,B,s) of 4-partition, where:

The set AA has 4∣M∣=O(q3)4|M|=O(q^{3}) elements, one for each occurrence of a member of W∪X∪YW\cup X\cup Y in a triple in MM and one for each triple in MM.

Let W={w1,…,wq}W=\{w_{1},\ldots,w_{q}\}, X={x1,…,xq}X=\{x_{1},\ldots,x_{q}\}, and Y={y1,…,yq}Y=\{y_{1},\ldots,y_{q}\}. Given any z∈W∪X∪Yz\in W\cup X\cup Y, let N(z)N(z) denote the number of triples in MM that contain zz, and let z,z,…,z[N(z)]z,z,\ldots,z[N(z)] denote the elements in AA corresponding to zz. Let r=32qr=32q, and define

Let ulu_{l} denote the single element corresponding to a particular triple ml=(wi,xj,yk)∈Mm_{l}=(w_{i},x_{j},y_{k})\in M. For any such ulu_{l}, let s(ul)=10r4−kr3−jr2−ir+8s(u_{l})=10r^{4}-kr^{3}-jr^{2}-ir+8.

Note that max⁡{s(a)∣a∈A}≤216∣A∣4\max\{s(a)|a\in A\}\leq 2^{16}|A|^{4}. This means 4-Partition is NP-complete in the strong sense Garey & Johnson (1979). It can be checked that this reduction is injective.

From 4-Partition to 3-Partition (cf. Theorem 4.4 in Garey & Johnson (1979)):

The 3-Partition problem is to decide, given a set AA of size 3m3m, a positive integer bound BB, a positive integer size s(a)s(a) for each a∈Aa\in A such that B/4<s(a)<B/2B/4<s(a)<B/2 and ∑a∈As(a)=mB\sum_{a\in A}s(a)=mB, whether AA can be partitioned into mm disjoint subsets A1,…,AmA_{1},\ldots,A_{m}, each of size three, such that, for each 1≤i≤m1\leq i\leq m, ∑a∈Ais(a)=B\sum_{a\in A_{i}}s(a)=B. We denote an instance of 3-Partition by the tuple (A,m,B,s)(A,m,B,s).

The reduction in Garey & Johnson (1979) maps an instance (A,m,B,s)(A,m,B,s) of 4-Partition, with ∣A∣=4m|A|=4m and max⁡{s(a)∣a∈A}≤216∣A∣4\max\{s(a)|a\in A\}\leq 2^{16}|A|^{4}, to the instance (A′,m′,B′,s′)(A^{\prime},m^{\prime},B^{\prime},s^{\prime}) of 3-Partition, where A′A^{\prime} has m′=O(m2)m^{\prime}=O(m^{2}) elements: one element wiw_{i} for each element aia_{i} of AA, two elements ui,ju_{i,j} and uˉi,j\bar{u}_{i,j} for each pair (ai,aj)(a_{i},a_{j}) of elements from AA, and 8m2−3m8m^{2}-3m filler elements uk∗u^{*}_{k}, 1≤k≤8m2−3m1\leq k\leq 8m^{2}-3m. Their sizes are:

We let B′=64B+4B^{\prime}=64B+4. It can be checked that this reduction is injective.

From 3-Partition to Machine Flow, the decision version of the two-machine flow scheduling problem with unit processing times defined in Chapter 3 in Yu (1996) (where it is called F2UD’):

The Machine Flow problem is to decide, given two machines M1 and M2, each of which can process at most one job at a time, and nn jobs jj, 1≤j≤n1\leq j\leq n, where each job takes unit processing time and the job jj is assigned a delay ljl_{j} that describes the minimum amount of time between the completion of the job jj on M1 and its start on M2, and a threshold yy, whether there exists a feasible schedule of the jobs so that the last job is completed before time yy.

The reduction in Yu (1996) from 3-Partition to Machine Flow goes as follows. Without loss of generality, we consider a modified version of 3-Partition by multiplying the partition elements by 4m4m. Thus we are given a set of positive integers A={a1,…,a3m}A=\{a_{1},\ldots,a_{3m}\} and a positive integer BB such that (1) B<ai<2BB<a_{i}<2B for all ii, (2) ∑jaj=4mB\sum_{j}a_{j}=4mB, (3) ai=0a_{i}=0 (mod mm) for all ii, and (4) 4B=04B=0 (mod mm). The problem is to decide if AA can be partitioned into mm disjoint 3-element subsets A1,…,AmA_{1},\ldots,A_{m} such that ∑aj∈Aiaj=4B\sum_{a_{j}\in A_{i}}a_{j}=4B, for all ii. An instance of this modified version of 3-Partition is mapped to an instance of Machine Flow with delays (1) lj=ajl_{j}=a_{j} for 1≤j≤3m1\leq j\leq 3m, (2) lj=0l_{j}=0 for 3m+1≤j≤4mB3m+1\leq j\leq 4mB, (3) lj=u+1l_{j}=u+1 for 4mB+1≤j≤mu4mB+1\leq j\leq mu, where u=4(m+1)Bu=4(m+1)B, and (4) the threshold y=n+4mB+2y=n+4mB+2, where n=mun=mu is the total number of jobs. It can be checked that this reduction is injective.

From Machine Flow to RN3DM (Restricted Numerical 3-Dimensional Matching, cf. page 31 in Yu, 1996):

The RN3DM problem is to decide, given a positive integer set U={u1,…,un}U=\{u_{1},\ldots,u_{n}\} and a positive integer ee such that ∑j=1nuj+n(n+1)=ne\sum_{j=1}^{n}u_{j}+n(n+1)=ne, whether there exist two nn-permutations λ\lambda and μ\mu such that j+λ(j)+uμ(j)=ej+\lambda(j)+u_{\mu(j)}=e for 1≤j≤n1\leq j\leq n. (It can be assumed that each ui<e−1u_{i}<e-1).

The reduction (Corollary 3 on page 32 in Yu, 1996) maps an instance of Machine Flow to that of RN3DM given by uj=lju_{j}=l_{j} for 1≤j≤n1\leq j\leq n, and e=ye=y. We assume that the instance of Machine Flow here arises in the reduction from 3-Partition to Machine Flow given in (III) above. This will ensure that ∑juj+n(n+1)=ne\sum_{j}u_{j}+n(n+1)=ne and each ui<e−1u_{i}<e-1, which we require for (V) below to be injective. It can be checked that this reduction is injective.

From RN3DM to RNMTS (Restricted Numerical Matching with Target Sums, cf. Brunetti et al., 2008):

The RNMTS problem is to decide, given positive integers y1,…,yny_{1},\ldots,y_{n} such that 2≤y1≤y2≤⋯≤yn≤2n2\leq y_{1}\leq y_{2}\leq\cdots\leq y_{n}\leq 2n and ∑iyi=n(n+1)\sum_{i}y_{i}=n(n+1), if there exist nn-permutations σ\sigma and π\pi such that σ(k)+π(k)=yk\sigma(k)+\pi(k)=y_{k} for 1≤k≤n1\leq k\leq n.

RN3DM is mapped to RNMTS by letting yj=e−ujy_{j}=e-u_{j} and then reordering the yjy_{j} as per their values. It can be checked that this reduction is injective.

The Permutation problem (page 69 in Brunetti et al., 2008, where it is called Permutation (S3S_{3})) is to decide, given non-negative integers z2,…,z2n∈{0,…,n}z_{2},\ldots,z_{2n}\in\{0,\ldots,n\}, whether there exists an n×nn\times n permutation matrix PP such that ∑i,j:i+j=lPi,j=zl\sum_{i,j:i+j=l}P_{i,j}=z_{l} for 2≤l≤2n2\leq l\leq 2n.

The reduction in Brunetti et al. (2008) maps an instance y=(y1,…,yn)y=(y_{1},\ldots,y_{n}) of RNMTS to an instance z=(z2,…,z2n)z=(z_{2},\ldots,z_{2n}) of Permutation by setting zl=∣{k≤n ∣ yk=l}∣z_{l}=|\{k\leq n\ |\ y_{k}=l\}|. It can be checked that this reduction is injective.

Special Consistency is the problem addressed in 3.6, namely, the problem of deciding positivity of tμ,πλt^{\lambda}_{\mu,\pi}, given λ,μ\lambda,\mu and π\pi in unary, when (λT,μT,πT)(\lambda^{T},\mu^{T},\pi^{T}) is restricted to be of the form 6.

The reduction in Brunetti et al. (2001) maps an instance z=(z2,…,z2n)z=(z_{2},\ldots,z_{2n}) of Permutation to an instance (λ,μ,π)(\lambda,\mu,\pi) of Special Consistency satisfying 6, with r=n−1r=n-1 and di=zi+2d_{i}=z_{i+2}, 0≤i≤2r0\leq i\leq 2r (it can be shown that ∑kdk=r+1\sum_{k}d_{k}=r+1, and ∑kkdk=r(r+1)\sum_{k}kd_{k}=r(r+1)). This reduction is injective.

From Special Consistency to Restricted Kronecker:

The reduction given in the proof of 4.2 is also injective.

1.5 follows from 5.1 and 5.3. This proof also shows that the superpolynomially many partition triples in 1.5 can be constructed explicitly (as defined in 1.5).

In the preceding proof, we can use, in place of 3D Matching, any problem in NP which has a polynomial-time-computable padding function (Berman & Hartmanis, 1977), and which can be reduced by a polynomial-time one-one reduction to Restricted Kronecker. For example, Sat also has a polynomial-time-computable padding function, and it can also be reduced by a polynomial-time one-one reduction to Restricted Kronecker. This reduction is obtained by composing the injective reduction from Sat to 3D Matching given in Garey & Johnson (1979) with the injective reduction from 3D Matching to Restricted Kronecker given in the proof of 5.3.

Though the reduction ϕ\phi in 5.3 is polynomial-time computable, the blow-up in size can be substantial. For example, let us start with a trivial instance of the 3D Matching problem, wherein q=2q=2, W={w1,w2}W=\{w_{1},w_{2}\}, X={x1,x2}X=\{x_{1},x_{2}\}, Y={y1,y2}Y=\{y_{1},y_{2}\}, and M={(w1,x1,y1),(w2,x1,y2),(w1,x2,y2)}M=\{(w_{1},x_{1},y_{1}),(w_{2},x_{1},y_{2}),(w_{1},x_{2},y_{2})\}. Clearly, MM does not contain a matching.

It can be checked that ϕ(M,W,X,Y,q)\phi(M,W,X,Y,q), with ϵ=1\epsilon=1 in condition 3, is a partition triple whose height is >1016>10^{16} and the total size is >1046>10^{46}. By 5.3, the Kronecker coefficient associated with this partition triple is zero. One cannot verify this fact directly using a computer, since computation of Kronecker coefficients for partition triples of this height and size is far beyond the reach of computer algebra systems. Thus 5.3 maps instances of 3D Matching which do not contain matching for trivial reasons to partition triples whose associated Kronecker coefficients vanish for highly nontrivial reasons. Thus the image of something trivial is highly nontrivial. This happens because of the nontriviality of the sequence of reductions that produce the image.

2 Proof of 1.6

For the proof of 1.6, we need the following lemma, which proves a variant of the result in Fortune (1979) that coNP-complete languages cannot be sparse unless P=NP.

Let L\mathcal{L} be a coNP-hard language given as a disjoint union

where L′\mathcal{L}^{\prime} is sparse (i.e., there are only poly⁡(n)\operatorname{poly}(n) words of length nn in L′\mathcal{L}^{\prime}) and L′′∈\mboxNP∩\mboxcoNP\mathcal{L}^{\prime\prime}\in\mbox{NP}\cap\mbox{coNP}. Then \mboxcoNP=\mboxNP\mbox{coNP}=\mbox{NP}.

We will show that the assumptions imply that \mboxSATc\mbox{SAT}^{c} (the complement of SAT) is in NP – this would imply that \mboxcoNP⊆\mboxNP\mbox{coNP}\subseteq\mbox{NP}, and hence, \mboxcoNP=\mboxNP\mbox{coNP}=\mbox{NP}. For this, we adapt the proof in Mahaney (1982); Fortune (1979).

Since L\mathcal{L} is coNP-hard, there exists a polynomial-time many-one reduction RR such that R(\mboxSAT)⊆LcR(\mbox{SAT})\subseteq{\mathcal{L}}^{c} and R(\mboxSATc)⊆LR(\mbox{SAT}^{c})\subseteq\mathcal{L}. Since L′′\mathcal{L}^{\prime\prime} is in \mboxNP∩\mboxcoNP\mbox{NP}\cap\mbox{coNP}, there exist non-deterministic Turing machines M1M_{1} and M2M_{2} such that, given input xx, M1M_{1} halts (in polynomial time) if and only if x∈L′′x\in\mathcal{L}^{\prime\prime}, while M2M_{2} halts (in polynomial time) if and only if x∉L′′x\not\in\mathcal{L}^{\prime\prime}.

Let FF be a formula for which we have to decide unsatisfiability. We perform depth-first search on the binary tree obtained by self-reducing FF (the root of this tree is FF, and the children of a node GG are G0G_{0} and G1G_{1}, the formulas of smaller size obtained by specializing the first variable in GG to true or false, and applying trivial simplifications), starting at the root node. We maintain a table U\mathcal{U} of labels (RR-values) of unsatisfiable formulae, starting with U:={R(\mboxfalse)}\mathcal{U}:=\{R(\mbox{false})\}. At each node GG, we first compute R(G)R(G) and then do one of the following:

If R(G)∈UR(G)\in\mathcal{U}, prune the subtree and return to the parent node.

Otherwise, if G=\mboxtrueG=\mbox{true}, enter an infinite loop.

Otherwise, run both non-deterministic Turing machines M1M_{1} and M2M_{2} in parallel on the input R(G)R(G) until one of the two halts (which will always happen, for some sequence of non-deterministic choices, in polynomial time):

If M1M_{1} halts (in which case R(G)∈L′′⊆LR(G)\in\mathcal{L}^{\prime\prime}\subseteq\mathcal{L}, and hence, GG is unsatisfiable), add R(G)R(G) to U\mathcal{U}, prune the subtree and return to the parent node.

If M2M_{2} halts, visit both children G0G_{0} and G1G_{1}. Upon return (if this happens), it will always be true that G0G_{0} and G1G_{1} are unsatisfiable, and hence R(G0),R(G1)∈UR(G_{0}),R(G_{1})\in\mathcal{U} and GG is unsatisfiable. Thus add R(G)R(G) to U\mathcal{U} and return to the parent node.

It is clear that this algorithm can be understood as a non-deterministic Turing machine that halts if and only if FF is unsatisfiable.

It suffices to show that, if FF is unsatisfiable, this algorithm halts in polynomial time. For this, it suffices to show that the number of interior nodes that are visited by the algorithm is polynomial in the size ∣F∣|F| of the formula FF (since the tree is binary, the number of visited leaves is at most twice the number of visited interior nodes). Now observe that interior nodes only arise in the case where M2M_{2} halts on input R(G)R(G), in which case R(G)∈L′R(G)\in\mathcal{L}^{\prime}. Thus any interior node is necessarily labeled by an element of the sparse set L′\mathcal{L}^{\prime}. We can thus conclude the argument precisely as in Lemma 2.2 of Mahaney (1982): If GG and G′G^{\prime} are two interior nodes that have the same label, R(G)=R(G′)∈L′R(G)=R(G^{\prime})\in\mathcal{L}^{\prime}, then they necessarily ought to appear in the same branch of the search tree (because we proceed by depth-first search). As the depth of the tree is no more than mm – the number of variables in FF – we find that each label can occur at most mm times. Therefore, the number of visited interior nodes can be upper bounded by m⋅p(q(∣F∣))m\cdot p(q(\lvert F\rvert)), where qq is a polynomial that bounds the increase in length induced by the reduction RR and p=p(n)p=p(n) is a polynomial that bounds the number of strings of length ≤n\leq n in the sparse set L′\mathcal{L}^{\prime}. We conclude that \mboxSATc∈\mboxNP\mbox{SAT}^{c}\in\mbox{NP}.

Another ingredient needed for the proof of 1.6 is the following result.

The problem of deciding if (λ,μ,π)∈Kron⁡(m)(\lambda,\mu,\pi)\in\operatorname{Kron}(m) is in \mboxNP∩\mboxcoNP\mbox{NP}\cap\mbox{coNP}. Here, mm denotes the maximum height of λ\lambda, μ\mu, or π\pi, and the partition triple (λ,μ,π)(\lambda,\mu,\pi) is given in unary.

For given mm, let L\mathcal{L} be the set of partition triples (λ,μ,μ)(\lambda,\mu,\mu) satisfying the constraints 1–\short6. Let L′\mathcal{L}^{\prime} be the set of partition triples (λ,μ,μ)(\lambda,\mu,\mu) satisfying both the constraints 1–\short6 and 7–\short8. Let L′′\mathcal{L}^{\prime\prime} be the set of partition triples satisfying the constraints 1–\short6 such that either \short7 or \short8 in 1.6 are violated. Then, clearly, L=L′∪L′′\mathcal{L}=\mathcal{L}^{\prime}\cup\mathcal{L}^{\prime\prime}.

In the definition of L′′\mathcal{L}^{\prime\prime}, we can drop the constraint \short1, since it is automatically satisfied if \short7 or \short8 are violated (as kμ,μλ=kμt,μλt=kμt,μtλk^{\lambda}_{\mu,\mu}=k^{\lambda^{t}}_{\mu^{t},\mu}=k^{\lambda}_{\mu^{t},\mu^{t}}). By 5.9, the problem of deciding whether a partition triple belongs to the Kronecker cone is in \mboxNP∩\mboxcoNP\mbox{NP}\cap\mbox{coNP}. It follows that L′′∈\mboxNP∩\mboxcoNP\mathcal{L}^{\prime\prime}\in\mbox{NP}\cap\mbox{coNP}. By 4.2, L\mathcal{L} is coNP-hard. It now follows from 5.7 that L′\mathcal{L}^{\prime} is not sparse, assuming \mboxcoNP≠\mboxNP\mbox{coNP}\not=\mbox{NP}. This proves 1.6.

There seems to be a surprising correlation between the complexities of kμ,πλk^{\lambda}_{\mu,\pi} and tμ,πλt^{\lambda}_{\mu,\pi}. On the one hand, positivity of kμ,πλk^{\lambda}_{\mu,\pi} is, in general, NP-hard to decide (3.7), just as it is for tμ,πλt^{\lambda}_{\mu,\pi} (3.6). On the other hand, suppose Π\Pi is a subclass of partition triples such that the problem of deciding positivity of tμ,πλt^{\lambda}_{\mu,\pi}, for (λ,μ,π)∈Π(\lambda,\mu,\pi)\in\Pi, is in PP. While the corresponding problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi}, for (λ,μ,π)∈Π(\lambda,\mu,\pi)\in\Pi, may not always be in PP, the results in this section suggest that it may indeed be so for many “natural” subclasses Π\Pi. In particular, 6.13 proved in this section suggests that the problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi}, when μ\mu and π\pi are rectangular (=δ(λ)=\delta(\lambda)), is in PP, as conjectured in Mulmuley (2010b).

We begin with a lemma that is needed for proving these results.

Let λ,μ,π\lambda,\mu,\pi be partitions of dd. An obstruction predesign is defined to have type (λ,μ,π)(\lambda,\mu,\pi) if the number of columns in λ\lambda of length kk equals the number of hyperedges in layer 1 with kk vertices, and the number of columns in μ\mu of length kk equals the number of hyperedges in layer 2 with kk vertices, and the number of columns in π\pi of length kk equals the number of hyperedges in layer 3 with kk vertices.

To each vertex we can assign its triple of hyperedges. An obstruction design is an obstruction predesign such that no two vertices have the same triple of hyperedges.

2 Littlewood-Richardson coefficients

3 Partitions of constant height

It is known that positivity of kμ,πλk^{\lambda}_{\mu,\pi} can be decided in polynomial time when λ,μ\lambda,\mu, and π\pi have constant heights. This is consistent with the following result.

The algorithm is a hybrid algorithm based on the number of boxes ∣λ∣|\lambda|. The values for tt in the case ∣λ∣<(c+2)c|\lambda|<(c+2)c are stored in a database of constant size. The case ∣λ∣≥(c+2)c|\lambda|\geq(c+2)c is trivial, as the following lemma shows.

If ht⁡(λ)≤c\operatorname{ht}(\lambda)\leq c, ht⁡(μ)≤c\operatorname{ht}(\mu)\leq c, ht⁡(π)≤c\operatorname{ht}(\pi)\leq c, and ∣λ∣=∣μ∣=∣π∣≥(c+2)c|\lambda|=|\mu|=|\pi|\geq(c+2)c, then tμ,πλ>0t^{\lambda}_{\mu,\pi}>0.

If ∣λ∣|\lambda| is divisible by cc, then we arrange the vertices in a rectangular array whose columns contain cc vertices each. Otherwise we add an extra column containing less than cc vertices. The crucial property is that since ∣λ∣≥(c+2)c|\lambda|\geq(c+2)c each row contains at least c+2c+2 vertices.

Proceeding column-wise from top to bottom and from left to right, we greedily assign vertices to hyperedges of the first layer according to the column lengths of μ\mu. Note that each hyperedge constructed thus lies either in a single column or in two adjacent columns. Likewise, proceeding row-wise from left to right and from top to bottom, we greedily assign vertices to hyperedges of the second layer according to the column lengths of π\pi. Since each row contains at least c+2c+2 vertices, a layer 2 hyperedge cannot contain two vertices from the same or adjacent columns.

Note that tμ,πλ>0t_{\mu,\pi}^{\lambda}>0 and tμ′,π′λ′>0t_{\mu^{\prime},\pi^{\prime}}^{\lambda^{\prime}}>0 implies tμ+μ′,π+π′λ+λ′>0t_{\mu+\mu^{\prime},\pi+\pi^{\prime}}^{\lambda+\lambda^{\prime}}>0, and that Lemma 6.5 shows that for constant height the semigroup of triples with positive tμ,πλt_{\mu,\pi}^{\lambda} is finitely generated, as it is known for kμ,πλk_{\mu,\pi}^{\lambda} as well.

4 When one partition is a hook

Blasiak (2017) has a given a #P\#P-formula for kμ,πλk^{\lambda}_{\mu,\pi} when λ\lambda is a hook. The problem of deciding positivity of kμ,πλk^{\lambda}_{\mu,\pi} in this case may be conjectured to be in PP in view of the following result.

Positivity of tμ,πλt_{\mu,\pi}^{\lambda}, given λ,μ\lambda,\mu, and π\pi in unary, can be decided in polynomial time if λ\lambda is a hook.

Fix α\alpha and β\beta. If kk is larger than the number of (α,β)(\alpha,\beta)-equivalence classes, then by the pigeonhole principle the λ\lambda-hyperedge of size kk must contain two (α,β)(\alpha,\beta)-equivalent vertices. Therefore this construction does not yield an obstruction design. If kk is smaller than the number of (α,β)(\alpha,\beta)-equivalence classes, the λ\lambda-hyperedge of size kk can be chosen to contain pairwise (α,β)(\alpha,\beta)-nonequivalent vertices. The other vertices are singletons in the λ\lambda layer, so this construction yields an obstruction design.

A solution with flow at least kk exists iff there exist α\alpha and β\beta such that the number of (α,β)(\alpha,\beta)-equivalence classes is at least kk.

Given α\alpha and β\beta with at least kk equivalence classes we construct a solution to the flow problem by sending one flow unit for each equivalence class: For the equivalence class corresponding to the iith μ\mu-column and jjth π\pi-column we send a unit from the source vertex to the iith μ\mu-vertex, from there to the jjth π\pi-vertex and then to the sink vertex. This satisfies the capacity constraints and is a solution to the flow problem that sends at least kk flow units.

From a solution of the max flow problem we readily generate a solution to a relaxed max flow problem where we remove the capacities on the edges from the μ\mu-vertices to the π\pi-vertices. We send flow units on additional arbitrary paths from the source to the sink. Once all capacities are saturated we are guaranteed to send exactly DD flow units. From this new solution we construct set partitions α\alpha and β\beta by defining that the size of the (α,β)(\alpha,\beta)-equivalence class corresponding to the iith μ\mu-column and the jjth π\pi-column is the amount of flow from the iith μ\mu-vertex to the jjth π\pi-vertex. So if the original solution had at least kk flow units, then there are at least kk (α,β)(\alpha,\beta)-equivalence classes in our construction.

5 Rectangular Kronecker coefficients

It is conjectured in Mulmuley (2010b) that the problem deciding positivity of the rectangular Kronecker coefficient kδ(λ),δ(λ)λk^{\lambda}_{\delta(\lambda),\delta(\lambda)} is in PP. This is supported by the following result.

Let λ\lambda be any partition with drdr boxes and at most min⁡(d2,r2)\min(d^{2},r^{2}) rows, and let δ=δ(λ)=(d,…,d)\delta=\delta(\lambda)=(d,\ldots,d) (rr times). Then tδ,δλ>0t_{\delta,\delta}^{\lambda}>0. In particular, the problem of deciding positivity of tδ,δλt_{\delta,\delta}^{\lambda} is trivial.

Since kδ,δλ=0k_{\delta,\delta}^{\lambda}=0 if ht⁡(λ)>min(d2,r2)\operatorname{ht}(\lambda)>min(d^{2},r^{2}), the constraint on λ\lambda here is very natural.

The case d≥rd\geq r is easier, so we handle this case first. We have to construct an obstruction design with drdr vertices and go about it as follows. Let iremdi\mathop{\text{rem}}d denote the remainder when dividing ii by dd. The vertex set VV is a subset of the d×dd\times d grid {(i,j)∣0≤i,j<d}\{(i,j)\mid 0\leq i,j<d\}. We have (i,j)∈V(i,j)\in V iff (i+j)remd∈{0,1,…,r−1}(i+j)\mathop{\text{rem}}d\in\{0,1,\ldots,r-1\}. For example, for r=4r=4 and d=6d=6, the vertex set is arranged as follows (row is at the top, column is at the left):

Note that every row and every column has exactly rr boxes. The rows correspond to the hyperedges of the first layer, where the columns correspond to the hyperedges of the second layer. Note that no matter how the hyperedges of the third layer are placed, no two vertices can share all three hyperedges, because no two vertices even share their two hyperedges in layer one and two. Therefore an arbitrary placement of the third layer shows tδ,δλ>0t_{\delta,\delta}^{\lambda}>0.

For d<rd<r an analogous construction can be made, but several vertices share a location, see the example r=6r=6 and d=4d=4 below.

Note that in this construction, if r/d>2r/d>2, then three or more vertices lie at the same position. As in the case d≥rd\geq r, the rows correspond to the hyperedges of the first layer, where the columns correspond to the hyperedges of the second layer. But now the third layer cannot be placed arbitrarily, but care has to be taken. The hyperedges can be placed in any order, but not at arbitrary positions. When a hyperedge is placed, it first uses those places where several vertices are grouped together (and of course only uses one from each such place). If there are places with more than two vertices, the hyperedge first takes vertices from those places with the most vertices. This greedy method ensures that no hyperedge contains a pair of vertices from the same place, because by the length restriction on λ\lambda a hyperedge cannot use more than d2d^{2} vertices.

This work was supported by NSF grant CCF-1017760. MW acknowledges support by the Simons Foundation, FQXI, and AFOSR (grant no. FA9550-16-1-0082). A part of this work was done at the Simons Institute for the Theory of Computing, Berkeley.

References