This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Computational complexity of counting coincidences

Swee Hong Chan Department of Mathematics, Rutgers University, Piscatway, NJ 08854. [email protected]  and  Igor Pak Department of Mathematics, UCLA, Los Angeles, CA 90095. [email protected]
Abstract.

Can you decide if there is a coincidence in the numbers counting two different combinatorial objects? For example, can you decide if two regions in 3\mathbb{R}^{3} have the same number of domino tilings? There are two versions of the problem, with  2×1×12\times 1\times 1  and  2×2×12\times 2\times 1  boxes. We prove that in both cases the coincidence problem is not in the polynomial hierarchy unless the polynomial hierarchy collapses to a finite level. While the conclusions are the same, the proofs are notably different and generalize in different directions.

We proceed to explore the coincidence problem for counting independent sets and matchings in graphs, matroid bases, order ideals and linear extensions in posets, permutation patterns, and the Kronecker coefficients. We also make a number of conjectures for counting other combinatorial objects such as plane triangulations, contingency tables, standard Young tableaux, reduced factorizations and the Littlewood–Richardson coefficients.

Key words and phrases:
Domino tiling, permanent, #\#P-completeness, graph matching, linear extensions of posets, order ideals of posets, matroid bases, Kronecker coefficient, standard Young tableau

1. Introduction

1.1. Tilings

In this paper we consider coincidences of combinatorial counting functions. Consider two bounded regions in the plane. Do they have the same number of domino tilings? Here we are assuming that the regions are finite subsets of unit squares on a square grid, we write  Γ2\Gamma\subset\mathbb{Z}^{2}, and the dominos are the usual  2×12\times 1  rectangles. For example, a  2×32\times 3  rectangle has three domino tilings, and both regions on the right have four domino tilings:

[Uncaptioned image]

Algorithmically, the coincidence problem is easy, since we can simply compute the number of domino tilings of each region in polynomial time, and then compare the numbers. Indeed, the Kasteleyn formula  gives the number  τ(Γ)\tau(\Gamma)  of domino tilings as an  n×nn\times n determinant which can be computed in time polynomial in nn, where  n=|Γ|n=|\Gamma|  is the area of region Γ\Gamma, see e.g. [Ken04, LP09].

Now consider what happens to domino tilings in 3\mathbb{R}^{3}. Should we expect that the coincidence problem remains easy? It turns out, this is a much harder problem. In fact, there are two versions of 33-dimensional dominoes: a  2×2×12\times 2\times 1  box which we call a slab, and a  2×1×12\times 1\times 1  box which we call a brick. For example, there are three tilings of a  2×2×22\times 2\times 2  box with a slab and nine with a brick, shown below up to rotations:

[Uncaptioned image]

For both versions, there is no analogue of the Kasteleyn formula for the number of tilings. Indeed, for tilings with slabs or with bricks, Jed Yang and the second author proved that computing the number of tilings is #P-complete [PY13a]. But this is where the similarities end: the coincidence problems have very different nature in these two cases. The reason for the difference is the type of reduction used in the proofs (see Section 4 and §\S6.2).

Let  τs(Γ)\tau_{s}(\Gamma)  be the number of tilings of a region  Γ3\Gamma\subset\mathbb{Z}^{3}  with slabs. Denote by  CST{}_{\text{ST}}  the slab tiling coincidence problem:

(1.1) CST:={τs(Γ)=?τs(Γ),whereΓ,Γ3}.\text{\sc C${}_{\text{ST}}$}\ :=\ \big{\{}\hskip 1.70709pt\tau_{s}(\Gamma)\hskip 0.85355pt=^{?}\hskip 0.85355pt\tau_{s}(\Gamma^{\prime}),\ \ \text{where}\ \ \Gamma,\Gamma^{\prime}\subset\mathbb{Z}^{3}\big{\}}.
Theorem 1.1.

Problem  CST{}_{\text{ST}}  is  coNP-hard. Furthermore,  CST{}_{\text{ST}}  is not in the polynomial hierarchy unless the polynomial hierarchy collapses to a finite level:   CST{}_{\text{ST}}PHPH=Σmp\in\textup{{{PH}}}\ \Rightarrow\ \textup{{{PH}}}=\Sigma^{{\textup{p}}}_{m}  for some mm.

In particular, the theorem implies that  CST{}_{\text{ST}}  does not have a polynomial time algorithm (unless P=NP\textup{{{P}}}=\textup{{{NP}}}). Nor does it have a probabilistic polynomial time algorithm (unless PH collapses), since  BPPPH\textup{{{BPP}}}\subseteq\textup{{{PH}}}. Nor does there exist a polynomial size witness for the coincidence (ditto).

The first part of the theorem follows immediately from NP-completeness of the slab tileability [PY13a, Thm 1.1], which is a special cases of the  CST{}_{\text{ST}}. The second part follows from the parsimonious reduction in the proof of #P-completeness of the  CST{}_{\text{ST}}, combined with some known computational complexity (see §\S4.1).

Let  τb(Γ)\tau_{b}(\Gamma)  be the number of tilings of a region  Γ3\Gamma\subset\mathbb{Z}^{3}  with bricks. We similarly denote by  CBT{}_{\text{BT}}  the brick tiling coincidence problem:

(1.2) CBT:={τb(Γ)=?τb(Γ),whereΓ,Γ3}.\text{\sc C${}_{\text{BT}}$}\ :=\ \big{\{}\hskip 1.70709pt\tau_{b}(\Gamma)\hskip 0.85355pt=^{?}\hskip 0.85355pt\tau_{b}(\Gamma^{\prime}),\ \ \text{where}\ \ \Gamma,\Gamma^{\prime}\subset\mathbb{Z}^{3}\big{\}}.
Theorem 1.2.

Problem  CBT{}_{\text{BT}}  is not in the polynomial hierarchy unless the polynomial hierarchy collapses to a finite level:   CBT{}_{\text{BT}}PHPH=Σmp\in\textup{{{PH}}}\ \Rightarrow\ \textup{{{PH}}}=\Sigma^{{\textup{p}}}_{m}  for some mm.

Note that we make no claims that  CBT{}_{\text{BT}}  is  coNP-hard. Proof of that would require a major advance in computational complexity. However, the theorem does imply that  CBT{}_{\text{BT}}  is not in  P,  BPP,  NP,  coNP, etc., unless  PH=Σmp\textup{{{PH}}}=\Sigma^{{\textup{p}}}_{m}. For the proof we again need some standard results in computational complexity, combined with a curious combinatorial result of independent interest (Theorem 3.1).

1.2. Beyond tilings

The paper starts with proofs of two theorems above which allows us to develop tools to prove similar results for the coincidence problem of many other combinatorial counting functions. We group the results into two, by analogy with the two theorems above.

Let  f#Pf\in\textup{{{\#P}}}  be a counting function. The coincidence problemCf  is defined as follows:

Cf:={f(x)=?f(y)}.\text{\sc C}_{f}\ :=\ \big{\{}\hskip 0.85355ptf(x)\hskip 1.70709pt=^{?}\hskip 0.85355ptf(y)\hskip 0.85355pt\big{\}}.
Theorem 1.3.

The coincidence problem  Cf  is not in the polynomial hierarchy unless the polynomial hierarchy collapses to a finite level, where  ff  is either one of the following:

(0)(0)  the number of satisfying assignments of a  3SAT  formula,

(1)(1)  the number of proper 33-colorings of a planar graph,

(2)(2)  the number of Hamiltonian cycles in a graph,

(3)(3)  the number  PCπ(σ)\mathrm{PC}_{\pi}(\sigma)  of patterns π\hskip 0.85355pt\pi\hskip 0.85355pt in a permutation σ\hskip 0.85355pt\sigma,

(4)(4)  the Kronecker coefficient  g(λ,μ,ν)g(\lambda,\mu,\nu)  for partitions  λ,μ,νn\lambda,\mu,\nu\vdash n  given in unary.

Furthermore, the coincidence problem  Cf  is  coNP-hard in all these cases.

As we explain, the theorem follows from known complexity results. This is in sharp contrast with our next theorem where each part requires additional work. We need just one definition.

A rational matroid  is a matroid that is realizable over \mathbb{Q}. We assume that this matroid is given by a set of nn vectors in d\mathbb{Q}^{d}. In this presentation, computing the number  b(M)b(M)  of bases of a rational matroid MM is known to be #P-complete [Sno12] (cf. §\S6.3).

Theorem 1.4 (main theorem).

The coincidence problem  Cf  is not in the polynomial hierarchy unless the polynomial hierarchy collapses to a finite level, where  ff  is either of the following:

(0)(0)  the number of perfect matchings  PM(P)\mathrm{PM}(P)  in a simple bipartite graph,

(1)(1)  the number of satisfying assignments of a  MONOTONE 2SAT  formula,

(2)(2)  the number of independent sets  λ(G)\lambda(G)  in a planar bipartite graph,

(3)(3)  the number of order ideals  μ(P)\mu(P)  of a poset,

(4)(4)  the number of linear extensions  e(P)e(P)  of a 22-dimensional poset,

(5)(5)  the number of matchings  ma(G){\rm ma}(G)  in a planar bipartite graph,

(6)(6)  the number of bases  b(M)b(M)  of a rational matroid.

Note that all counting functions in Theorems 1.3 and 1.4 are  #P-complete. Unfortunately, that by itself does not imply that the corresponding coincidence problems are not in PH (see §\S6.18). Example 5.15 (see also §\S6.4) has an especially notable variation on part (4)(4) in the theorem for posets of height two. While the number of linear extensions is known to be  #P-complete in this case, it is open whether the corresponding coincidence problem is in PH. Similarly, a variation on part (6)(6) for bicircular matroids gives another example of this type (see §\S6.3). The number of bases is known to be  #P-complete in this case, but corresponding coincidence problem remains unexplored.

Paper structure

We start with basic definitions and notation in a short Section 2. In Section 3, we give some preliminary results on domino tilings. In Section 4, we prove Theorems 1.1, 1.2 and 1.3. In a lengthy Section 5, we discuss further examples, and prove Main Theorem 1.4. We conclude with many final remarks and open problems in Section 6.

2. Definitions and notation

General notation

Let  [n]={1,,n}[n]=\{1,\ldots,n\}  and  ={0,1,2,}\mathbb{N}=\{0,1,2,\ldots\}. For a sequence  a=(a1,,am)\textbf{{a}}=(a_{1},\ldots,a_{m}), denote  |a|:=a1++am|a|:=a_{1}+\ldots+a_{m} . Similarly, for the integer partitions  μλ\mu\subset\lambda, let the size|λ/μ||\lambda/\mu|  be the number of squares in the skew Young diagram λ/μ\lambda/\mu. For  |λ|=n|\lambda|=n  we also write  λn\lambda\vdash n.

Combinatorics

We think of  dd\mathbb{Z}^{d}\subset\mathbb{R}^{d}  both as a lattice and a collection of the corresponding unit dd-cube. A regionΓd\Gamma\subset\mathbb{Z}^{d}  is a subset of dd-cubes. Denote by  |Γ||\Gamma|  the size of Γ\Gamma, which can also be viewed as the volume of the union of the corresponding unit cubes. Region  Γ\Gamma  is called simply connected if the union of the corresponding (closed) dd-cubes is simply connected. A tile  in  d\mathbb{R}^{d}  is a finite simply connected region. For a set of tiles  T={t1,,tm}T=\{t_{1},\ldots,t_{m}\}, a tiling  is a disjoint union of copies of tiles tit_{i} (unless stated otherwise, parallel translations, rotations and reflections are allowed).

We assume that the reader is familiar with basic notions in algebraic combinatorics, such as  standard Young tableaux, Kostka numbers, Littlewood–Richardson  and Kronecker coefficients. Defining them, explaining their importance, combinatorial interpretations and properties would take too much space and be a distraction. We refer the reader to [Mac95, Man01, Sta12] and further references sprinkled throughout the paper.

Complexity

We assume that the reader is familiar with basic notions and results in computational complexity and only recall a few definitions. We use standard complexity classes  P,  NP, coNP,  #P,  Σmp\Sigma^{{\textup{p}}}_{m}  and  PH. The notation  {a=?b}\{a=^{?}b\}  is used to denote the decision problem whether  a=ba=b. We use the oracle notationKL{}^{\text{\sf L}}  for two complexity classes  K, L PSPACE\subseteq\textup{{{PSPACE}}}, and the polynomial closure  \langleA\rangle for a problem  A PSPACE\in\textup{{{PSPACE}}}. We will also use less common classes

GapP:={fgf,g#P}andC=P:={f(x)=?g(y)f,g#P}.\textup{{{GapP}}}:=\{f-g\mid f,g\in\textup{{{\#P}}}\}\quad\text{and}\quad{\textup{{{C${}_{=}$P}}}}:=\{f(x)=^{?}g(y)\mid f,g\in\textup{{{\#P}}}\}.

The distinction between binary  and unary  presentation will also be important. We refer to [GJ78] and [GJ79, §\S4.2] for the corresponding notions of NP-completeness and strongNP-completeness.

We also that assume the reader is familiar with standard decision and counting problems, such as  2SAT,  MONOTONE 2SAT,  3SAT,  1-in-3 SAT,  HAMILTON CYCLE,  #2SAT,  #3SAT,  PERMANENT, etc. Occasionally, we conflate counting functions ff and the problems of computing ff. We hope this does not lead to a confusion.

We refer to [AB09, Gol08, MM11, Pap94] for definitions and standard results in computational complexity. See [GJ79] for the classical introduction and a long list of NP-complete problems. See also [Pak22, §\S13] for a recent overview of #P-complete problems in combinatorics. For surveys on counting complexity, see [For97, Sch90].

3. Counting tilings

3.1. Domino tilings

Denote by  𝒯(n)\mathcal{T}(n)  the set of numbers of domino tilings over all regions of size 2n2n:

𝒯(n):={τ(Γ),whereΓ2,|Γ|=2n}.\mathcal{T}(n)\ :=\ \big{\{}\hskip 1.70709pt\tau(\Gamma)\hskip 1.70709pt,\ \ \text{where}\ \ \Gamma\subset\mathbb{Z}^{2},\ |\Gamma|\hskip 0.85355pt=\hskip 0.85355pt2n\hskip 1.70709pt\big{\}}.

Clearly,  𝒯(n){0,1,,4n}\mathcal{T}(n)\subseteq\{0,1,\ldots,4^{n}\}  since each domino tilings is determined by the 44 choices for a domino at every even square. The following result proves a converse:

Theorem 3.1.

There is a constant  c>1c>1, such that  𝒯(n){0,1,,cn}\mathcal{T}(n)\supseteq\{0,1,\ldots,c^{n}\}, for all  n1n\geq 1. Moreover, for all  kcnk\leq c^{n}, a region  Γ2\Gamma\subset\mathbb{Z}^{2}  with  τ(Γ)=k\tau(\Gamma)=k  and  |Γ|=2n|\Gamma|=2n, can be constructed in time polynomial in nn.

Remark 3.2.

It was known before that  n𝒯(n)=\cup_{n}\mathcal{T}(n)=\mathbb{N}, i.e. that every nonnegative integer is the number of domino tilings of some  region. This was shown by Philippe Nadeau with an elegant explicit construction.111See mathoverflow.net/a/178100 . Unfortunately, this construction has  τ(Γ)=Θ(n)\tau(\Gamma)=\Theta(n)  where  n=|Γ|n=|\Gamma|, and thus does not give our theorem.

In a different direction, Brualdi and Newman in [BN65] gave (in the language of permanents), an explicit construction of a simple bipartite graph on  nn  vertices with exactly kk perfect matchings, for all  0k2n10\leq k\leq 2^{n-1}. These graphs have unbounded degree and thus very far from being grid graphs (or even planar graphs). There was a recent constant factor improvement in [GT18], but relatively little attention otherwise (see [OEIS, A089477]), compared to the corresponding determinant problem (see §\S6.10).

Proof of Theorem 3.1.

For an integer  k0k\geq 0, we give an explicit construction of a region  Γ2\Gamma\subset\mathbb{Z}^{2}  with  τ(G)=k\tau(G)=k  and  |Γ|=O(logk)|\Gamma|=O(\log k). This implies the result.

A square  (i,j)2(i,j)\in\mathbb{Z}^{2}  is called even (odd ) if  i+ji+j  is even (odd). Denote by  Γe=Γe2\Gamma_{e}=\Gamma\cap\mathbb{Z}_{e}^{2}  and  Γo:=Γo2\Gamma_{o}:=\Gamma\cap\mathbb{Z}_{o}^{2}  the sets of even and odd squares in Γ\Gamma, respectively. Clearly, if  |Γe||Γo||\Gamma_{e}|\neq|\Gamma_{o}|  then  τ(Γ)=0\tau(\Gamma)=0. We use  Γxy\Gamma-x-y  to denote Γ\Gamma with squares x,yx,y removed. Denote by  𝒟(a,b)\mathcal{D}(a,b)  the set of regions  Γ\Gamma such that  τ(Γ)=a\tau(\Gamma)=a  and  τ(Γxy)=b\tau(\Gamma-x-y)=b  for some  xΓex\in\Gamma_{e}  and  yΓoy\in\Gamma_{o}.

Below we give constructions which prove the following implications:

(3.1) 𝒟(a,1)𝒟(2a,1),\displaystyle\mathcal{D}(a,1)\neq\varnothing\ \Longrightarrow\ \mathcal{D}(2a,1)\neq\varnothing\hskip 0.85355pt,
𝒟(a,1)𝒟(2a+1,1).\displaystyle\mathcal{D}(a,1)\neq\varnothing\ \Longrightarrow\ \mathcal{D}(2a+1,1)\neq\varnothing\hskip 0.85355pt.

Starting with  𝒟(1,1)\mathcal{D}(1,1)\neq\varnothing  and iterating these in  O(logk)O(\log k)  times gives the desired  Γ𝒟(k,1)\Gamma\in\mathcal{D}(k,1).

For a region  Γ\Gamma  and squares  x,yΓx,y\in\Gamma, we say that  (Γ,x,y)(\Gamma,x,y)  is a CC-triple  if

\circ   x=(i,j)Γex=(i,j)\in\Gamma_{e}  and  y=(i,j+7)Γoy=(i,j+7)\in\Gamma_{o} ,

\circ   (u,v)Γ(u,v)\notin\Gamma  for all  u>iu>i,

\circ   (u,v)Γ(u,v)\notin\Gamma  for all  i1uii-1\leq u\leq i  and  j+1vj+6j+1\leq v\leq j+6.

Refer to caption
Figure 3.1. Region Γ𝒟(1,1)\Gamma\in\mathcal{D}(1,1), and two transformations  Γ𝒟(a,1)Γ𝒟(2a,1)\Gamma\in\mathcal{D}(a,1)\hskip 0.85355pt\Rightarrow\hskip 0.85355pt\Gamma^{\prime}\in\mathcal{D}(2a,1), and  Γ𝒟(a,1)Γ′′𝒟(2a+1,1)\Gamma\in\mathcal{D}(a,1)\hskip 0.85355pt\Rightarrow\hskip 0.85355pt\Gamma^{\prime\prime}\in\mathcal{D}(2a+1,1).

We start with a CC-triple  (Γ,x,y)𝒟(1,1)(\Gamma,x,y)\in\mathcal{D}(1,1)  as in Figure 3.1. We then define two transformations  (Γ,x,y)(Γ,x,y)(\Gamma,x,y)\to(\Gamma^{\prime},x^{\prime},y^{\prime})  and  (Γ,x,y)(Γ′′,x′′,y′′)(\Gamma,x,y)\to(\Gamma^{\prime\prime},x^{\prime\prime},y^{\prime\prime})  as in the figure, which prove the implications (3.1). Note that in both cases we obtain two CC-triples by adding at most 30 squares. This completes the construction. ∎

Corollary 3.3.

In notation of the proof above, we have  𝒟(a,b)\mathcal{D}(a,b)\neq\varnothing  for all  a,b0a,b\geq 0.

Proof.

By modifying our two transformations, one can show that

(3.2) 𝒟(a,b)𝒟(b,a),\displaystyle\mathcal{D}(a,b)\neq\varnothing\ \Longrightarrow\ \mathcal{D}(b,a)\neq\varnothing\hskip 0.85355pt,
𝒟(a,b),𝒟(a,b)𝒟(aa+bb,bb).\displaystyle\mathcal{D}(a,b)\neq\varnothing\hskip 0.85355pt,\ \mathcal{D}(a^{\prime},b^{\prime})\neq\varnothing\ \Longrightarrow\ \mathcal{D}(aa^{\prime}+bb^{\prime},bb^{\prime})\neq\varnothing\hskip 0.85355pt.

The first of these is given by  (Γ,x,y)(Γ,x,y)(\Gamma,x,y)\to(\Gamma^{\prime},x^{\prime},y^{\prime})  as in Figure 3.2. Similarly, the second is given by (Γ,x,y),(Γ,x,y)(Γ′′,x,y)(\Gamma,x,y),(\Gamma,x^{\prime},y^{\prime})\to(\Gamma^{\prime\prime},x,y^{\prime})  as in the figure. Note that  (Γ′′,x,y)(\Gamma^{\prime\prime},x,y^{\prime})  is no longer a CC-triple, so this transformation can be used only once.

In the proof above, we showed that  𝒟(n,1)\mathcal{D}(n,1)\neq\varnothing  for all  n0n\geq 0. By the first transformation in (3.2), this implies that  𝒟(1,m)\mathcal{D}(1,m)\neq\varnothing  for all  m0m\geq 0. Therefore, by the second transformation in (3.2), we have  𝒟(m+n,m)\mathcal{D}(m+n,m)\neq\varnothing. We then have  𝒟(m,m+n)\mathcal{D}(m,m+n)\neq\varnothing  by the first (after a modification where xx^{\prime} and yy^{\prime} are above and placed below xx^{\prime} and yy^{\prime}, respectively). Note that the second transformation is used only here, and thus the construction is well defined.

The remaining cases in  𝒟(n,0)\mathcal{D}(n,0)  and  𝒟(0,n)\mathcal{D}(0,n), follow from the two transformations above applied to  𝒟(1,0)\mathcal{D}(1,0). Alternatively, they follow the last two transformations in Figure 3.2; the details are straightforward. This finishes the proof. ∎

Refer to caption
Figure 3.2. Four transformations:  Γ𝒟(a,b)Γ𝒟(b,a)\Gamma\in\mathcal{D}(a,b)\hskip 0.85355pt\Rightarrow\hskip 0.85355pt\Gamma^{\prime}\in\mathcal{D}(b,a),   Γ𝒟(a,b)\Gamma\in\mathcal{D}(a,b),  Γ𝒟(a,b)Γ′′𝒟(aa+bb,bb)\Gamma^{\prime}\in\mathcal{D}(a^{\prime},b^{\prime})\hskip 0.85355pt\Rightarrow\hskip 0.85355pt\Gamma^{\prime\prime}\in\mathcal{D}(aa^{\prime}+bb^{\prime},bb^{\prime}),  Γ𝒟(a,b)Γ𝒟(0,a)\Gamma\in\mathcal{D}(a,b)\hskip 0.85355pt\Rightarrow\hskip 0.85355pt\Gamma^{\prime}\in\mathcal{D}(0,a)   and   Γ𝒟(a,b)Γ′′𝒟(a,0)\Gamma\in\mathcal{D}(a,b)\hskip 0.85355pt\Rightarrow\hskip 0.85355pt\Gamma^{\prime\prime}\in\mathcal{D}(a,0).
Corollary 3.4.

Let  𝒟(a,b)\mathcal{D}^{\prime}(a,b)  be the set of regions  Γ\Gamma such that  τ(Γ)=a\tau(\Gamma)=a  and  τ(Γxy)=b\tau(\Gamma-x-y)=b  for some domino  (x,y)(x,y), where  x,yΓx,y\in\Gamma. Then  𝒟(a,b)\mathcal{D}^{\prime}(a,b)\neq\varnothing,  for all  0ba0\leq b\leq a.

Proof.

Note that if  (x,y)(x,y)  is a domino in Γ\Gamma, then  τ(Γxy)τ(Γ)\tau(\Gamma-x-y)\leq\tau(\Gamma). Thus the assumption  bab\leq a  in the claim. Now, for  Γ𝒟(ab,1)\Gamma\in\mathcal{D}(a-b,1)  where  b1b\geq 1, and  Γ𝒟(1,b)\Gamma^{\prime}\in\mathcal{D}(1,b), the second transformation in (3.2) gives a region  Γ′′𝒟(a,b)\Gamma^{\prime\prime}\in\mathcal{D}(a,b)  as in Figure 3.2. Removing one white domino and keeping the other, gives the desired region in  𝒟(a,b)\mathcal{D}^{\prime}(a,b).

Similarly, taking a region  Γ2\Gamma\subset\mathbb{Z}^{2}  with  τ(Γ)=a\tau(\Gamma)=a, attaching a domino (x,y)(x,y) to a top right square  zΓz\in\Gamma  gives region  Γ𝒟(a,0)\Gamma^{\prime}\in\mathcal{D}^{\prime}(a,0)  since  τ(Γ+x+y)=τ(Γ)\tau(\Gamma+x+y)=\tau(\Gamma)  and  τ(Γxz)=0\tau(\Gamma-x-z)=0. The details are straightforward. ∎

3.2. Slab tilings

Denote by  𝒯s(n)\mathcal{T}_{s}(n)  the set of numbers of tilings with slabs:

𝒯s(n):={τs(Γ),whereΓ3,|Γ|=4n}.\mathcal{T}_{s}(n)\ :=\ \big{\{}\hskip 1.70709pt\tau_{s}(\Gamma)\hskip 1.70709pt,\ \ \text{where}\ \ \Gamma\subset\mathbb{Z}^{3},\ |\Gamma|\hskip 0.85355pt=\hskip 0.85355pt4n\hskip 1.70709pt\big{\}}.
Theorem 3.5.

There is a constant  c>1c>1, such that  𝒯s(n){0,1,,cn}\mathcal{T}_{s}(n)\supseteq\{0,1,\ldots,c^{n}\}, for all  n1n\geq 1. Moreover, for all  kcnk\leq c^{n}, a region  Γ3\Gamma\subset\mathbb{Z}^{3}  with  τs(Γ)=k\tau_{s}(\Gamma)=k  and  |Γ|=2n|\Gamma|=2n, can be constructed in time polynomial in nn.

Note that the corresponding result for the set  𝒯b(n)\mathcal{T}_{b}(n)  of numbers of tiling with bricks, follows trivially from Theorem 3.1 since  𝒯(n)𝒯b(n)\mathcal{T}(n)\subseteq\mathcal{T}_{b}(n).

Proof of Theorem 3.5.

The result follows from the proof of Theorem 3.1. Indeed, for every region  Γ2\Gamma\subset\mathbb{Z}^{2}  we can take a 2-layered region  Γ2:=Γ×{0,1}3\Gamma_{2}:=\Gamma\times\{0,1\}\subset\mathbb{Z}^{3}. Assuming  Γ\Gamma  does not have a  2×22\times 2  square inside, we have  τs(Γ2)=τ(Γ)\tau_{s}(\Gamma_{2})=\tau(\Gamma). The result now follows from reductions in Figure 3.1, where the  2×22\times 2  square in the middle is replaced by a  3×33\times 3  square without a center square (see an example below). Note that the notion of CC-triples also needs to be adjusted accordingly. The details are straightforward. ∎

[Uncaptioned image]

4. Complexity of coincidences

4.1. Parsimonious reductions

Let  f#Pf\in\textup{{{\#P}}}  be a counting function. As in the introduction, denote by

Cf:={f(x)=?f(y)}\text{\sc C}_{f}\ :=\ \big{\{}\hskip 0.85355ptf(x)\hskip 1.70709pt=^{?}\hskip 0.85355ptf(y)\hskip 0.85355pt\big{\}}

the coincidence problem  for ff  (cf. 6.1). This problem naturally belongs to the complexity class

(4.1) C=P:={f(x)=?g(y)wheref,g#P }.{\textup{{{C${}_{=}$P}}}}:=\big{\{}\hskip 0.85355ptf(x)\hskip 1.70709pt=^{?}\hskip 0.85355ptg(y)\ \text{where}\ f,\hskip 1.70709ptg\hskip 0.85355pt\in\hskip 0.85355pt\textup{{{\#P}}}\hskip 0.85355pt\big{\}}.

Note that  coNPC=P\textup{{{coNP}}}\subseteq{\textup{{{C${}_{=}$P}}}}, by definition. The  #3SAT  coincidence problem  C#3SAT{}_{\text{\sc\#3SAT}}  is a complete problem in this class.

Function  ff  is said to have a parsimonious reduction  to gg, if there is an injection  ρ:xy\rho:x\to y  from the instances of ff to the instances of gg, which maps values of the functions:  f(x)=g(y)f(x)=g(y), and such that ρ\rho can be computed in polynomial time. If  #3SAT  has a parsimonious reduction to ff, then the coincidence problem  Cf  is  coNP-hard because the decision problem  3SAT  is  NP-complete. Moreover,  Cf  is  C=P-complete by the definition of  C=P. In general, there are other ways for a problem to be  C=P-complete, but having a parsimonious reduction is the most straightforward.

Proposition 4.1 (cf. Theorem 1.3, part (0)(0)).

Let  f#Pf\in\textup{{{\#P}}}  be a function with a parsimonious reduction from  #3SAT. Then the coincidence problem  Cf  is not in  PH, unless  PH=Σmp\textup{{{PH}}}=\Sigma^{{\textup{p}}}_{m}  for some mm.

Proof.

Suppose that  CfPH{}_{f}\in\textup{{{PH}}}. Then  CfΣnp{}_{f}\in\hskip 0.85355pt\Sigma^{{\textup{p}}}_{n}  for some nn. We then have:

(4.2) PH NPC=PNPC#3SATNPCfNPΣnpΣn+1p,\textup{{{PH}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{{\textup{{{C${}_{=}$P}}}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc C}_{\text{\sc\#3SAT}}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{{\sc C}}_{f}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\Sigma^{{\textup{p}}}_{n}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\Sigma^{{\textup{p}}}_{n+1}\,,

where the first inclusion is by Tarui [Tar91] (see also [Gre93]), the second inclusion follows since  C#3SAT{}_{\text{\sc\#3SAT}}  is a complete problem in C=P, the third follows from the parsimonious reduction of  #3SAT  to ff, and the last inclusion is by definition. This proves the desired collapse of the polynomial hierarchy. ∎

Proof of Theorem 1.1.

The NP-completeness of  SlabTileability  is proved in [PY13a] by a bijective reduction from  1-in-3 SAT. Since Schaefer’s reduction  of  1-in-3 SAT  from  3SAT  is also by a bijection [Sch78] (see also [GJ79, Problem LO4]), we conclude that the number of slab tilings has a parsimonious reduction from  #3SAT. The result now follows from Proposition 4.1. ∎

Proof of Theorem 1.3, parts (1)(1), (2)(2).

These follow immediately from Proposition 4.1 and parsimonious reductions from #3SAT  given in [GJ79]. ∎

4.2. Pattern containment

Let  πSk\pi\in S_{k}  and  σSn\sigma\in S_{n}. We say that  σ\sigmacontainsπ\pi  if there is  A={a1<<ak}[n]A=\{a_{1}<\ldots<a_{k}\}\subset[n], such that  σ(a1),,σ(ak)\sigma(a_{1}),\ldots,\sigma(a_{k})  has the same relative order as π\hskip 0.85355pt\pi. Denote by  PCπ(σ)\mathrm{PC}_{\pi}(\sigma)  the number of subsets AA as above. This counting functions is well studied in various settings, see e.g. [Kit11, Vat15].

Define the pattern containment coincidence problem:

CPC={PCπ(σ)=?PCπ(σ)}.\text{\sc C}_{\text{\sc PC}}\,=\,\big{\{}\hskip 0.85355pt\mathrm{PC}_{\pi}(\sigma)\hskip 1.70709pt=^{?}\hskip 0.85355pt\mathrm{PC}_{\pi^{\prime}}(\sigma^{\prime})\hskip 0.85355pt\big{\}}.
Proof of Theorem 1.3, part (3)(3).

It was shown in [BBL98, Cor. 4], that computing the pattern containment function  PCπ(σ)\mathrm{PC}_{\pi}(\sigma)  is #P-complete, and the proof uses a parsimonious reduction from  #3SAT. Now the theorem follows from Proposition 4.1. ∎

4.3. Kronecker coefficients

Consider the following application of the tools above to algebraic combinatorics. Let  g(λ,μ,ν)g(\lambda,\mu,\nu), where λ,μ,νn\lambda,\mu,\nu\vdash n, denote the Kronecker coefficients:

g(λ,μ,ν):=χλχμ,χν=1n!σSnχλ(σ)χμ(σ)χν(σ),g(\lambda,\mu,\nu)\,:=\,\langle\chi^{\lambda}\chi^{\mu},\chi^{\nu}\rangle\,=\,\frac{1}{n!}\hskip 1.70709pt\sum_{\sigma\in S_{n}}\hskip 1.70709pt\chi^{\lambda}(\sigma)\hskip 1.70709pt\chi^{\mu}(\sigma)\hskip 1.70709pt\chi^{\nu}(\sigma)\hskip 1.70709pt,

where  χλ\chi^{\lambda}  denote the irreducible SnS_{n} character corresponding to λn\lambda\vdash n. By definition,  g(λ,μ,ν)g(\lambda,\mu,\nu)\in\mathbb{N}  and is symmetric with respect to permutations of  λ,μ,ν\lambda,\mu,\nu. It is known that  gg  is in  GapP:=#P#P\textup{{{GapP}}}:=\textup{{{\#P}}}-\textup{{{\#P}}}, i.e. can be written as the difference of two functions in #P. It is a major open problem whether  gg  is in #P. Here and below we are assuming that partitions are given in unary. We refer to recent surveys [Pak22, Pan23] for further background.

Define the Kronecker coefficients coincidence problem:

CKRON={g(λ,μ,ν)=?g(α,β,γ)}.\text{\sc C}_{\text{\sc KRON}}\,=\,\big{\{}\hskip 0.85355ptg(\lambda,\mu,\nu)\hskip 1.70709pt=^{?}\hskip 0.85355ptg(\alpha,\beta,\gamma)\hskip 0.85355pt\big{\}}.
Proof of Theorem 1.3, part (4)(4).

It was shown in [IMW17] that the Kronecker vanishing problem{g(λ,μ,ν)=?0}\{g(\lambda,\mu,\nu)=^{?}0\}  is  coNP-hard. This proves the first part. For the second part, recall that computing  gg  is  #P-hard (in the unary) [IMW17], and that the original proof gives an explicit (although rather involved) parsimonious reduction from  #3SAT  to g\hskip 0.85355ptg. Now the theorem follows from Proposition 4.1. ∎

4.4. The permanent

We start with a more general problem which will be used to demonstrate our approach. Recall that the  0/1 Permanent  is the benchmark #P-complete problem, which corresponds to counting perfect matchings in simple bipartite graphs.

Let  CPER{}_{\text{\sc PER}}  denote the  0/1 Permanent Coincidence problem:

CPER:={per(M)=?per(M)M,Mn},\text{\sc C${}_{\text{\sc PER}}$}\ :=\ \big{\{}\mathrm{per}(M)\,=^{?}\,\mathrm{per}(M^{\prime})\hskip 1.70709pt\mid\hskip 1.70709ptM,\hskip 0.85355ptM^{\prime}\in\mathcal{M}_{n}\big{\}},

where  n\mathcal{M}_{n}  is the set of n×nn\times n  matrices with entries in {0,1}\{0,1\}.

Theorem 4.2 (= Theorem 1.4, part (0)(0)).

CPER{}_{\text{\sc PER}} PH\notin\textup{{{PH}}}  unless  PH=Σmp\textup{{{PH}}}=\Sigma^{{\textup{p}}}_{m}  for some m\hskip 0.85355ptm.

Proof.

Let  VPER{}_{\text{\sc PER}}  denote the  0/1 Permanent Verification problem:

VPER:={per(M)=?k},\text{\sc V${}_{\text{\sc PER}}$}\ :=\ \big{\{}\mathrm{per}(M)\,=^{?}\,k\big{\}},

where  MM  is a 0/1 matrix and  kk\in\mathbb{N}  is given in binary. We have:

(4.3) PH P#PNPVPERNPCPER.\textup{{{PH}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{P}}}^{\textup{{{\#P}}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc V${}_{\text{\sc PER}}$}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc C${}_{\text{\sc PER}}$}\rangle}.

The first inclusion is Toda’s theorem [Toda91]. The second inclusion follows because  0/1 PERMANENT  is  #P-complete [Val79a, Val79c]. Indeed, transform every query to a #P oracle to a 0/1 permanent instance, guess the value kk of that permanent, then call  VPER{}_{\text{\sc PER}}  to check that this guess is correct.222Here we are using Ryan Williams’s approach in  cstheory.stackexchange.com/a/53024/

For the third inclusion in (4.3), use Theorem 3.1 to construct a region  Γ2\Gamma\subset\mathbb{Z}^{2}  of size  |Γ|=O(logk)c|\Gamma|=O(\log k)^{c}  and with exactly kk domino tilings. The dual bipartite graph GG then has exactly kk perfect matchings.333We can also use the Brualdi–Newman construction here (see Remark 3.2), instead of Theorem 3.1. This corresponds to an instance of the 0/1 permanent equal to kk. Now call  CPER{}_{\text{\sc PER}}  to simulate  VPER{}_{\text{\sc PER}}.

Finally, suppose  CPER{}_{\text{\sc PER}} PH\in\textup{{{PH}}}. Then  CPER{}_{\text{\sc PER}} Σnp\in\hskip 0.85355pt\Sigma^{{\textup{p}}}_{n}  for some nn. By (4.3), this implies:

(4.4) PH NPCPERNPΣnpΣn+1p,\textup{{{PH}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc C${}_{\text{\sc PER}}$}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\Sigma^{{\textup{p}}}_{n}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\Sigma^{{\textup{p}}}_{n+1}\,,

as desired. ∎

Proof of Theorem 1.2.

Recall that the problem of counting  2×1×12\times 1\times 1  brick tilings is #P-complete [PY13a, Thm 1.3], see also [Val79b]. Now, note that the  CBT{}_{\text{BT}}  restricted to the case of regions in 2\mathbb{Z}^{2} still satisfies conclusions of Theorem 3.1. From this point, the proof follows verbatim the proof of Theorem 4.2. ∎

Remark 4.3.

The proof in [PY13a] is via a parsimonious reduction from the number of perfect matchings in 33-regular bipartite graphs, which in turn is proved  #P-complete via a non-parsimonious reduction from the  0/1 PERMANENT[DL92, Thm 6.2]. This is why we cannot proceed by analogy with the proof of Theorem 1.1.

5. Variations on the theme

In this section we use the tools that we developed to solve other coincidence problems.

5.1. Complete functions

Let  𝒳=𝒳n\mathcal{X}=\sqcup\hskip 0.85355pt\mathcal{X}_{n}, where  𝒳n{0,1}γ(n)\mathcal{X}_{n}\subseteq\{0,1\}^{\gamma(n)}  be a set of combinatorial objects, which means that the input size  γ(n)\gamma(n)  is polynomial in nn. A counting function  f#Pf\in\textup{{{\#P}}}  can be viewed as a function  f:𝒳f:\mathcal{X}\to\mathbb{N}. For example, let  𝒳\mathcal{X}  be the set of simple graphs  G=(V,E)G=(V,E)  where  n=|V|n=|V|, and let  f(G)f(G)  be the number of Hamiltonian cycles in GG.

Denote

𝒯f(n):={f(x)x𝒳n}and𝒯f:={f(x)x𝒳}\mathcal{T}_{f}(n)\,:=\,\{f(x)\hskip 1.70709pt\mid\hskip 1.70709ptx\in\mathcal{X}_{n}\}\quad\text{and}\quad\mathcal{T}_{f}\,:=\,\{f(x)\hskip 1.70709pt\mid\hskip 1.70709ptx\in\mathcal{X}\}

the set of all values of the function ff. We say that  ff  is complete  if  𝒯f=\mathcal{T}_{f}=\mathbb{N}. We say that  ff  is almost complete  if  𝒯f\mathcal{T}_{f}  contain all but a finite number of integers  kk\in\mathbb{N}.

For example, by Theorem 3.1, the number of domino tilings in the plane is a complete function. Trivial examples of complete functions include the number of 33-cycles in a simple graphs, or the number  inv(σ)\operatorname{{\rm inv}}(\sigma)  of inversions in a permutation σ\sigma. Similarly, the number of spanning trees in a simple graph is in  {0,1,3,4,}\{0,1,3,4,\ldots\}, and thus almost complete.

Example 5.1 (Linear extensions).

Let  P=(X,)P=(X,\prec)  be a finite poset with  n=|X|n=|X|  elements. A linear extension of PP is a bijection  ρ:X[n]\rho:X\to[n], such that  ρ(x)<ρ(y)\rho(x)<\rho(y)  for all  xyx\prec y. Denote by  e(P)e(P)  the number of linear extensions of PP. We use  #LE  to denote the problem of computing  e(P)e(P). It is known that  #LE  is  #P-complete [BW91].

Clearly,  e(P)1e(P)\geq 1  for all PP. Note that the set of numbers of linear extensions  𝒯e={1,2,3,}\mathcal{T}_{e}=\{1,2,3,\ldots\}, since  e(Ck+C1)=ke(C_{k}+C_{1})=k, where we take a parallel sum of an element and a chain of length kk. In particular, the function ee is almost complete.

Example 5.2 (Symmetric Kronecker coefficients).

For a slightly nontrivial example, consider the symmetric Kronecker coefficientgs(λ):=g(λ,λ,λ)g_{s}(\lambda):=g(\lambda,\lambda,\lambda), see [PP22]. Note that function  gs:{λ}g_{s}:\{\lambda\}\to\mathbb{N}  is complete:

gs(12)=0,gs(1)=1andgs(4k,2k)=k+1for allk1,g_{s}(1^{2})\hskip 0.85355pt=\hskip 0.85355pt0\hskip 0.85355pt,\quad g_{s}(1)\hskip 0.85355pt=\hskip 0.85355pt1\quad\text{and}\quad g_{s}(4k,2k)\hskip 1.70709pt=\hskip 1.70709ptk+1\ \ \text{for all}\ \ k\geq 1,

where  (4k,2k)6k=n(4k,2k)\vdash 6k=n, see e.g. [Ste14].

Example 5.3 (Littlewood–Richardson coefficients).

The  Littlewood–Richardson ((LR)) coefficientscμνλc^{\lambda}_{\mu\nu}  can be defined as structure constants for the ring of Schur functions:

sμsν=λcμνλsλ,s_{\mu}\hskip 0.85355pt\cdot\hskip 0.85355pts_{\nu}\ =\ \sum_{\lambda}\,c^{\lambda}_{\mu\nu}\hskip 1.70709pts_{\lambda}\hskip 1.70709pt,

see e.g. [Sta12, Ch. 7]. It remains open whether computing  cμνλc^{\lambda}_{\mu\nu}  is #P-complete (in unary), see an extensive discussion in [Pak22] and [Pan23]. We note that this is a complete function:

c1,123=0,c1,12=1andck(21),k(21)k(321)=k+1for allk1,c^{3}_{1,1^{2}}\hskip 0.85355pt=\hskip 0.85355pt0,\quad c^{2}_{1,1}\hskip 0.85355pt=\hskip 0.85355pt1\quad\text{and}\quad c^{k(321)}_{k(21),\hskip 0.85355ptk(21)}\hskip 1.70709pt=\hskip 1.70709ptk+1\ \ \text{for all}\ \ k\geq 1,

where the last equality follows form  c21,21321=2c^{321}_{21,21}=2  combined with [Ras04, Rem. 5.2].

Example 5.4 (Contingency tables).

Let  a=(a1,,ar)r\textbf{{a}}=(a_{1},\ldots,a_{r})\in\mathbb{N}^{r}  and  b=(b1,,bs)s\textbf{{b}}=(b_{1},\ldots,b_{s})\in\mathbb{N}^{s}. Denote by  CT(a,b)\text{CT}(\textbf{{a}},\textbf{{b}})  the number of contingency tablesM=(xij)rsM=(x_{ij})\in\mathbb{N}^{rs}, defined by

j=1sxij=aifor all i ,i=1rxij=bjfor all j ,andxij0for all i,j .\sum_{j=1}^{s}\hskip 0.85355ptx_{ij}\hskip 1.70709pt=\hskip 1.70709pta_{i}\ \ \text{for all \ $i$}\hskip 0.85355pt,\quad\sum_{i=1}^{r}\hskip 0.85355ptx_{ij}\hskip 1.70709pt=\hskip 1.70709ptb_{j}\ \ \text{for all \ $j$}\hskip 0.85355pt,\quad\text{and}\quad x_{ij}\hskip 0.85355pt\geq\hskip 0.85355pt0\ \ \text{for all \ $i,j$}\hskip 0.85355pt.

Note that  CT(a,b)1\text{CT}(\textbf{{a}},\textbf{{b}})\geq 1  for all  |a|=|b||\textbf{{a}}|=|\textbf{{b}}| .

Computing the number  CT(a,b)\text{CT}(\textbf{{a}},\textbf{{b}})  of contingency tables (with the input in unary) is conjectured to be  #P-complete [Pak22, §\S13.4.1]. On the other hand,  CT()\text{CT}(\cdot)  is clearly a complete function, since e.g.  CT(a,a)=k+1\text{CT}(\textbf{{a}},\textbf{{a}})=k+1, for  m=n=2m=n=2  and  a=(k,k)\textbf{{a}}=(k,k). Using standard reductions this also implies that the Kostka number  is also a complete function, see e.g. [PV10].

Example 5.5 (Pattern containment).

Let  σSn\sigma\in S_{n}. Clearly, the pattern containment function  PC21(σ)=inv(σ)\mathrm{PC}_{21}(\sigma)=\operatorname{{\rm inv}}(\sigma)  is complete. The following result is a generalization:

Proposition 5.6.

Fix  πSk\pi\in S_{k} , where  k2k\geq 2. Then function  PCπ(σ)\mathrm{PC}_{\pi}(\sigma)  is complete.

Proof.

Without loss of generality, assume that  π\pi  starts with an ascent. Let  σ\sigma  be a decreasing sequence  m,m1,,π(1)=am,m-1,\ldots,\pi(1)=a. Append it with  π(2),,π(k)\pi(2),\ldots,\pi(k)  which are shifted above  mm  if they are strictly greater than aa. Suppose exactly  rr  elements are shifted. Let  m=nrm=n-r  be so that the resulting σ\hskip 0.85355pt\sigma  is in  SnS_{n} , and observe that  PCπ(σ)=nra+1\mathrm{PC}_{\pi}(\sigma)=n-r-a+1. This implies the result.

For example, let  π=(2,5,1,3,6,4)\pi=(2,5,1,3,6,4). The construction above gives:

PC251364(n4,n5,,3,2,n1,1,n3,n,n2)=n5,\mathrm{PC}_{251364}\hskip 0.85355pt(n-4,n-5,\ldots,\hskip 0.85355pt3,\hskip 0.85355pt2,n-1,\hskip 0.85355pt1,n-3,n,n-2)\hskip 1.70709pt=\hskip 1.70709ptn-5,

where  k=6k=6,  a=2a=2,  r=4r=4, and  m=n4m=n-4. ∎

5.2. Concise functions

Note that being complete is neither necessary nor sufficient for our approach to the complexity of the coincidence problems. The following purely combinatorial definition gets us closer to the goal.

Definition 5.7.

Function  f:𝒳f:\mathcal{X}\to\mathbb{N}  is called concise  if there exist some fixed constants  C,c>0C,c>0, such that for all  k𝒯fk\in\mathcal{T}_{f}  there is an element  x𝒳nx\in\mathcal{X}_{n}  with  f(x)=kf(x)=k  and  n<C(logk)cn<C\hskip 0.85355pt(\log k)^{c}.

Our Theorem 3.1 shows that number  τ(Γ)\tau(\Gamma)  of domino tilings of regions  Γ2\Gamma\subset\mathbb{Z}^{2}  is concise. On the other hand, the numbers of patterns  PCπ(σ)(nk)\mathrm{PC}_{\pi}(\sigma)\leq\binom{n}{k}  for every  πSk\pi\in S_{k}  and  σSn\sigma\in S_{n} (see Example 5.5), so the function  PCπ\mathrm{PC}_{\pi}  is not concise for a fixed π\pi. We now present several less obvious examples of concise functions.

Example 5.8 (Independent sets).

Let  G=(V,E)G=(V,E)  be a finite simple graph, and let  λ(G)\lambda(G)  be the number of independent sets in GG, i.e. subsets  XVX\subseteq V  such that XX contains no two adjacent vertices. Recall that computing  λ\lambda  is  #P-complete, see [PB83]. Moreover, this holds even for planar bipartite graphs [Vad01].

We now show that  λ\lambda  is concise. Denote by GG^{\prime} the graph obtained from GG by adding a new vertex ww and adding all edges from ww to VV. Similarly, denote by G′′G^{\prime\prime} the graph obtained from GG by adding a new vertex ww disconnected from VV. Observe that  λ(G)=λ(G)+1\lambda(G^{\prime})=\lambda(G)+1  and  λ(G′′)=2λ(G)\lambda(G^{\prime\prime})=2\lambda(G). Iterating these two operations we obtain the desired graph GG on nn vertices, with  λ(G)=k\lambda(G)=k  and  n=O(logk)n=O(\log k).

Example 5.9 (Order ideals).

Let  P=(X,)P=(X,\prec)  be a finite poset with  n=|X|n=|X|  elements. A subset  YXY\subseteq X  is a lower order ideal  if for all  xyx\prec y  where  xXx\in X  and  yYy\in Y, we also have  xYx\in Y. Denote by  μ(P)\mu(P)  the number of lower order ideals in PP. Recall that computing  μ\mu  is  #P-complete, see [PB83].

We can similarly show that  μ\mu  is concise by constructing posets  P,P′′P^{\prime},\hskip 0.85355ptP^{\prime\prime}  with an extra element xx that is either smaller than all elements in XX, or incomparable to XX. We then have  μ(P)=μ(P)+1\mu(P^{\prime})=\mu(P)+1,  μ(P′′)=2μ(P)\mu(P^{\prime\prime})=2\mu(P). Iterating these two operations proves the claim.

Remark 5.10.

For closely related problems on the smallest topology  with a given number of open sets, and the shortest addition chain  of a given integer, see [RT10], [Knu98, §\S4.6.3], and sequences [OEIS, A137814], [OEIS, A003064].

Example 5.11 (Satisfiability).

Recall that satisfiability decision problem  2SAT  is in P, while the corresponding counting problems  #MONOTONE 2SAT  is #P-complete [Val79a]. Here monotone  refers to boolean formulas which do not have negative variables.

Observe that the number of independent sets  λ(G)\lambda(G)  has an obvious parsimonious reduction to  #MONOTONE 2SAT:

G=(V,E)ΦG:=(v,w)E(xvxw),G=(V,E)\ \ \longrightarrow\ \ \Phi_{G}\,:=\,\bigwedge_{(v,w)\in E}\hskip 1.70709pt(x_{v}\vee x_{w}),

since complements to independent sets  VXV\smallsetminus X  are in natural bijection with satisfying assignments of ΦG\Phi_{G}. Since  λ(G)\lambda(G)  is concise, then so is  #MONOTONE 2SAT. Similarly,  #3SAT  is also concise, via the standard reduction:

G=(V,E)ΦG:=(zzz)(v,w)E(xvxwz¯).G=(V,E)\ \ \longrightarrow\ \ \ \Phi_{G}^{\prime}\,:=\,(z\vee z\vee z)\bigwedge_{(v,w)\in E}(x_{v}\vee x_{w}\vee\overline{z}).

The approach in this example can be distilled in the following basic observation:

Proposition 5.12.

Suppose a concise counting function  gg  has a parsimonious reduction to f\hskip 0.85355ptf. Then  ff  is also concise.

5.3. Further examples

We start with the following general notion of exponential growth which will prove useful in several examples.444This definition is somewhat non-standard as we make no distinction between weakly exponential, exponential, factorial and superexponential growths:  ene^{\sqrt{n}},  ene^{n},  enlogne^{n\log n}  and  en2e^{n^{2}}. We refer, e.g., to [Gri14] for a more refined treatment of growth functions.

Definition 5.13.

We say that a counting function  f:𝒳f:\mathcal{X}\to\mathbb{N}  has exponential growth  if there exist  A,α>0A,\alpha>0, such that  f(x)>Aexp(nα)f(x)>A\exp(n^{\alpha})  for all  x𝒳nx\in\mathcal{X}_{n} .

For example, the number of independent sets of a bipartite graph G\hskip 0.85355ptG\hskip 0.85355pt on n\hskip 0.85355ptn\hskip 0.85355pt vertices has exponential growth:  λ(G)2n/2\lambda(G)\geq 2^{n/2}. Note also that every  f#Pf\in\textup{{{\#P}}}  satisfies the opposite inequality:  f(x)<Bexp(nβ)f(x)<B\exp(n^{\beta})  for all  x𝒳nx\in\mathcal{X}_{n}  and some  B,β>0B,\beta>0.

Proposition 5.14.

Let  𝒳=𝒳n\mathcal{X}=\sqcup\hskip 0.85355pt\mathcal{X}_{n}  be a set of combinatorial objects. Suppose function  f:𝒳f:\mathcal{X}\to\mathbb{N}  has exponential growth. Suppose also that  ff  is almost complete. Then  ff  is concise.

Proof.

Since  ff  is almost complete, there exists an integer  NN, s.t. for every  k>Nk>N, we have  f(x)=kf(x)=k  for some  x𝒳nx\in\mathcal{X}_{n}. Since  ff  has exponential growth, we must have  k>Aexp(nα)k>A\exp(n^{\alpha}). Let

C:=max{n:f(x)Nfor somex𝒳n}.C\,:=\,\max\big{\{}\hskip 0.85355ptn\hskip 1.70709pt:\hskip 1.70709ptf(x)\leq N\ \ \text{for some}\ \ x\in\mathcal{X}_{n}\hskip 0.85355pt\big{\}}.

Then,  n<C+(logkA)1/αn<C+\big{(}\log\frac{k}{A}\big{)}^{1/\alpha}  for all  kk, so  ff  is concise. ∎

Example 5.15 (Linear extensions of restricted posets).

The height  and width  of a  PP  is the size of the longest chain and antichain, respectively. For posets of bounded width,  #LE is in P (via dynamic programming). On the other hand, for posets of height two,  #LE  is  #P-complete [DP18].

For a permutation  σSn\sigma\in S_{n}  consider a partial order  Pσ=([n],)P_{\sigma}=([n],\prec), where  iji\prec j  if and only if  1i<jn1\leq i<j\leq n  and  σ(i)<σ(j)\sigma(i)<\sigma(j). Such posets  PσP_{\sigma}  are called two-dimensional, and  #LE  is also  #P-complete for this family [DP18]. Note that all posets of width two also have dimension two, see e.g. [Tro95].

In [KS21], Kravitz and Sah show that function ee is concise on posets of width two (see also [CP24]). Additionally, they prove that

(5.1) 𝒯e(n){1,,cn/(logn)}for somec>1.\mathcal{T}_{e}(n)\,\supseteq\,\big{\{}1,\ldots,c^{n/(\log n)}\big{\}}\quad\text{for some}\ \ c>1.

In particular, they prove a  O(logkloglogk)O(\log k\log\log k)  bound on the minimal size of a poset with kk linear extensions, cf.  [OEIS, A160371]  and  [OEIS, A281723]. They also conjecture the  O(logk)O(\log k)  bound for posets of width two, and thus for general posets [KS21, Conj. 7.37.3 and 7.47.4].

On the other hand, for posets of height two, we have  n/2!n/2!e(P)n!\lfloor n/2\rfloor!\cdot\lceil n/2\rceil!\hskip 0.85355pt\leq\hskip 0.85355pte(P)\hskip 0.85355pt\leq\hskip 0.85355ptn! . Thus, we have  k=exp(nlogn+O(n))k=\exp\big{(}n\log n+O(n)\big{)}, i.e. function  ee  restricted to height two posets has exponential growth. Now Proposition 5.14 gives a somewhat unexpected result:

Proposition 5.16.

If the function  ee  restricted to height two posets is almost complete, then it is also concise.

This suggests the following conjecture:

Conjecture 5.17.

The function  ee  restricted to height two posets is almost complete.

By Proposition 5.16 this gives a strong improvement over (5.1) for general posets:

Proposition 5.18.

Conjecture 5.17 implies that

(5.2) 𝒯e(n){1,,enlogncn}for somec>0.\mathcal{T}_{e}(n)\,\supseteq\,\big{\{}1,\ldots,e^{n\log n-c\hskip 0.85355ptn}\big{\}}\quad\text{for some}\ \ c>0.

In particular, Conjecture 5.17 would imply [KS21, Conj. 7.37.3] mentioned above. See §\S6.4 for more on this conjecture, and [CP24] for connections to other conjectures in number theory.

Example 5.19 (Matchings).

Let  G=(V,E)G=(V,E)  be a simple graph on  n=|V|n=|V|  vertices. A matching  is a subset  XEX\subseteq E  of pairwise nonadjacent edges. Denote by  ma(G){\rm ma}(G)  number of matchings. In [Jer87], Jerrum showed that computing  ma{\rm ma}  is #P-complete, even when restricted to planar graphs. Vadhan extended this to planar bipartite graphs [Vad01]. The following result shows that  ma{\rm ma}  is concise even when restricted to forests.

Proposition 5.20.

The function   ma{\rm ma}  restricted to forests is concise.

Proof.

By analogy with the proof of Theorem 3.1, let  𝒟(a,b)\mathcal{D}(a,b)  be the set of trees TT such that  ma(T)=a{\rm ma}(T)=a  and  ma(Tx)=b{\rm ma}(T-x)=b  for some vertex xx in TT. Let TT^{\prime} be a tree obtained by adding a vertex yy and an edge (xy)(xy). Clearly,  ma(T)=a+b{\rm ma}(T^{\prime})=a+b  and  ma(Tx)=b{\rm ma}(T^{\prime}-x)=b. We conclude:

𝒟(a,b)𝒟(a+b,a)and𝒟(a+b,b),\mathcal{D}(a,b)\neq\varnothing\ \ \Longrightarrow\ \ \mathcal{D}(a+b,a)\neq\varnothing\quad\text{and}\quad\mathcal{D}(a+b,b)\neq\varnothing,

where in the first implication we have  (T,x)(T,y)(T,x)\to(T^{\prime},y), and the second implication we have  (T,x)(T,x)(T,x)\to(T^{\prime},x). Iterating this procedure starting with a single edge, we obtain a tree TabT_{a\hskip 0.85355ptb} for every relatively prime (a,b)(a,b), ab1a\geq b\geq 1. Following [KS21], we can think of pairs  (a,b)(a,b)  as vertices of the Calkin–Wilf tree, and conclude that for every prime aa there is an integer bb, such that some tree in  𝒟(a,b)\mathcal{D}(a,b)  has  O(logalogloga)O(\log a\hskip 0.85355pt\log\log a)  vertices.

For a given integer kk, take a prime factorization  k=a1ak=a_{1}\cdots a_{\ell}  and let  FF  be the union of the corresponding trees  Ti𝒟(ai,bi)T_{i}\in\mathcal{D}(a_{i},b_{i}). Continuing the analysis in [KS21], we conclude that the forest  FF  has  O(logkloglogk)O(\log k\hskip 0.85355pt\log\log k)  vertices and  ma(F)=k{\rm ma}(F)=k. ∎

Example 5.21 (Spanning trees).

Let  G=(V,E)G=(V,E)  be a simple graph and let  st(G){\rm st}(G)  be the number of spanning trees in GG. Only recently, Stong proved in [Sto22] that function  st{\rm st}  is concise using a technical number theoretic argument. This resolved an open problem which goes back to [Sed70]. More precisely, Stong proved that for all  kk  sufficiently large, there is a simple planar  graph  GG  on nn vertices with exactly  kk  spanning trees, and  n<C(logk)3/2/(loglogk)n<C\hskip 0.85355pt(\log k)^{3/2}/(\log\log k). Note also that the function  stFP{\rm st}\in\textup{{{FP}}}  since it can be computed by the matrix-tree theorem.

Example 5.22 (Pattern containment).

As we discussed in §\S4.2, the pattern containment function  PC(σ,π)\mathrm{PC}(\sigma,\pi)  has a parsimonious reduction from  #3SAT. From above, this immediately implies that function PC\hskip 0.85355pt\mathrm{PC}\hskip 0.85355pt is concise. It would be interesting to see a more direct construction of this result.

Example 5.23 (Young tableaux).

Denote by  SYT(λ/μ)\operatorname{{\rm SYT}}(\lambda/\mu)  the number of standard Young tableaux  of a skew shape λ/μ\lambda/\mu. Recall that computing  SYT(λ/μ)\operatorname{{\rm SYT}}(\lambda/\mu)  is in FP by the Aitken–Feit determinant formula, see e.g. [Sta12, Eq. (7.71)(7.71)]. Note also  SYT\operatorname{{\rm SYT}}  is almost complete, since  SYT(n1,1)=n1\operatorname{{\rm SYT}}(n-1,1)=n-1.

Question 5.24.

Is  SYT\operatorname{{\rm SYT}}  concise? In other words, does for all  k>0k>0  there exist partitions  μλ\mu\subset\lambda  which satisfy  SYT(λ/μ)=k\operatorname{{\rm SYT}}(\lambda/\mu)=k  and   |λ/μ|C(logk)c|\lambda/\mu|\leq C(\log k)^{c}, for some fixed  C,c>0C,c>0?

By definition, the function  SYT\operatorname{{\rm SYT}}  is a restriction of the function ee counting the number of linear extensions to two-dimensional posets corresponding to skew shapes. Therefore, if the answer to the question is positive, the proof will be very challenging (cf. §\S6.11). There is, however, some negative evidence.

Proposition 5.25.

Function  SYT\operatorname{{\rm SYT}}  restricted to straight shapes ((i.e. μ=)\mu=\varnothing), is  not  concise.

Proof.

Recall that  SYT(λ)=χλ(1)|n!\operatorname{{\rm SYT}}(\lambda)=\chi^{\lambda}(1)\hskip 0.85355pt|\hskip 0.85355ptn!  for all  λn\lambda\vdash n. Thus  SYT(λ)\operatorname{{\rm SYT}}(\lambda)  cannot be a prime >n>n. Since  SYT\operatorname{{\rm SYT}}  is almost complete even when restricted to straight shapes, it is not concise. ∎

Note that the argument in the proof does not apply to general skew shapes, as we can have large prime even for the zigzag shapeρk/ρk2\rho_{k}/\rho_{k-2}, where  ρk=(k,k1,,1)\rho_{k}=(k,k-1,\ldots,1), see §\S6.4. Note also that function  SYT\operatorname{{\rm SYT}}  restricted to straight shapes can have unbounded coincidences (beyond conjugation). This was proved in [Cra08] by an elegant construction.

Proposition 5.26.

Function  SYT\operatorname{{\rm SYT}}  restricted to self-conjugate straight shapes ((i.e. μ=\mu=\varnothing and λ=λ)\lambda=\lambda^{\prime}), is  not  almost complete.

Proof.

Observe that  SYT\operatorname{{\rm SYT}}  restricted to self-conjugate straight shapes, has exponential growth:  SYT(λ)2n(1o(1))\operatorname{{\rm SYT}}(\lambda)\geq 2^{n(1-o(1))}  for all  λ=λn\lambda=\lambda^{\prime}\vdash n. Indeed, let  hii(λ)=2ai+1h_{ii}(\lambda)=2a_{i}+1,  1id1\leq i\leq d, be the principal hook length  of λ\lambda, where  dd  is the size of the Durfee square  of λ\lambda. We have:

SYT(λ)(2a1a1)(2adad)22a11a122ad1ad2n2dnd/22n2nnn 2nO(nlogn),\operatorname{{\rm SYT}}(\lambda)\,\geq\,\binom{2a_{1}}{a_{1}}\cdots\binom{2a_{d}}{a_{d}}\,\geq\,\frac{2^{2a_{1}-1}}{\sqrt{a_{1}}}\,\cdots\,\frac{2^{2a_{d}-1}}{\sqrt{a_{d}}}\,\geq\,\frac{2^{n-2d}}{n^{d/2}}\,\geq\,\frac{2^{n-2\sqrt{n}}}{n^{\sqrt{n}}}\,\geq\,2^{n\hskip 0.85355pt-\hskip 0.85355ptO(\sqrt{n}\hskip 0.85355pt\log n)}\,,

since  dnd\leq\sqrt{n}  and  (2kk)22k1/k\binom{2k}{k}\geq 2^{2k-1}/\sqrt{k}  for all  k1k\geq 1. On the other hand, there are at most  np(n)=eO(n)n\hskip 0.85355ptp(n)=e^{O(\sqrt{n})}  possible values of  SYT\operatorname{{\rm SYT}}  on symmetric partitions of size at most nn. The disparity with the growth implies the result. ∎

Example 5.27 (Kronecker coefficients).

As we discussed in §\S4.3, the Kronecker coefficients function  g(λ,μ,ν)g(\lambda,\mu,\nu)  has a parsimonious reduction from  #3SAT. From above, this immediately implies that function g\hskip 0.85355ptg\hskip 0.85355pt is concise.

On the other hand, the symmetric Kronecker coefficients function  gs(λ)=g(λ,λ,λ)g_{s}(\lambda)=g(\lambda,\lambda,\lambda)  is more mysterious. Although  gsg_{s}  is complete (see Remark 5.2), it was proved only recently [PP20, Thm 1.3], that  maxλngs(λ)=expΩ(nα)\max_{\lambda\vdash n}\hskip 0.85355ptg_{s}(\lambda)\hskip 0.85355pt=\hskip 0.85355pt\exp\Omega(n^{\alpha}), via an explicit construction based on the approach in [IMW17]. The exact asymptotics  maxλngs(λ)=expΘ(nlogn)\max_{\lambda\vdash n}\hskip 0.85355ptg_{s}(\lambda)\hskip 0.85355pt=\hskip 0.85355pt\exp\Theta(n\log n)  was given in [PP22] by a nonconstructive argument.

Conjecture 5.28.

Symmetric Kronecker coefficient function gsg_{s} is concise.

Example 5.29 (Contingency tables).

In notation of Example 5.4, denote by  n:=|a|=|b|n:=|\textbf{{a}}|=|\textbf{{b}}|  the size  of the contingency table, where  ar\textbf{{a}}\in\mathbb{N}^{r}  and  bs\textbf{{b}}\in\mathbb{N}^{s}. Recall that  CT  is complete.

Conjecture 5.30.

The function  CT  is concise, i.e. for every  k>0k>0  there exist vectors  a,b\textbf{\em{a}},\textbf{\em{b}}  of size nn, such that  CT(a,b)=k\text{\em CT}(\textbf{\em{a}},\textbf{\em{b}})=k  and  nC(logk)cn\leq C\hskip 0.85355pt(\log k)^{c}, for some fixed  C,c>0C,c>0.

Recall that

CT(a,b)(1+rsn)n(1+nrs)rs 4n,\text{CT}(\textbf{{a}},\textbf{{b}})\,\leq\,\left(1+\frac{rs}{n}\right)^{n}\left(1+\frac{n}{rs}\right)^{rs}\,\leq\,4^{n}\hskip 1.70709pt,

where the second inequality is under assumption  rsnrs\leq n, see [PP20, Thm 1.1]. Moreover, this upper bound is tight up to lower order terms. On the other hand, the number of pairs  (a,b)r+s(\textbf{{a}},\textbf{{b}})\in\mathbb{N}^{r+s}  is  nO(r+s)n^{O(r+s)}, suggesting that  c2c\geq 2  in the conjecture. We refer to [Bar17, BLP23], for an overview of the known bounds on  CT(a,b)\text{CT}(\textbf{{a}},\textbf{{b}}), and further references. See also §\S6.6 and §\S6.12, for two closely related problems.

5.4. Back to coincidence problems

Let us first summarize what we know. Recall the notion of the coincidence problem  Cf  defined in (4.1). When  fFPf\in\textup{{{FP}}}, we trivially have  CfP{}_{f}\in\textup{{{P}}}. The examples include the number of standard Young tableaux of skew shapes, spanning trees in graphs, pattern containment of a fixed pattern, and linear extensions of width two posets.

When  ff  has a parsimonious reduction from  #3SAT, the complexity of  Cf  is given by Proposition 4.1. The examples are given in Theorem 1.3. When  ff  is  not  known to be either in  FP  nor  #P-hard, none of our tools apply. The examples include the number of contingency tables, the Littlewood–Richardson coefficients, and symmetric Kronecker coefficients (additional examples are given in Section 6).

Finally, as the examples above show, there are many cases when  ff  is #P-complete, but the corresponding decision problem is in P. The examples are given in the Main Theorem 1.4 which we are now ready to prove.

Proof of Main Theorem 1.4.

For  (1)(1),  (2)(2)  and (3)(3), recall that these counting functions are concise, and that computing them is #P-complete (see above). The rest of the proof follows verbatim the proof of Theorem 4.2.

For  (4)(4), the problem   #LE  is #P-complete, the counting function  ee  concise, but the proof is a little less straightforward. Indeed, the proof in [KS21] does not produce a polynomial time construction of the width two poset  QQ  with  n=O(logkloglogk)n=O(\log k\hskip 0.85355pt\log\log k)  elements and  e(Q)=ke(Q)=k. Even the first step requires factoring of kk which not known to be in polynomial time (cf. Remark 5.31).

The way to get around the issue is to guess  the width two poset QQ of size  nn  with  e(Q)=ke(Q)=k. Such poset exists by [KS21], and  {e(Q)=?k}\{e(Q)=^{?}k\}  can be verified in polynomial time  O((logk)c)O\big{(}(\log k)^{c}\big{)}  since  QQ  has width two and size  n=O(logkloglogk)n=O(\log k\hskip 0.85355pt\log\log k). Denote

Ve:={e(Q)=?k}andCe:={e(P)=?e(Q)},\text{\sc V}_{e}\hskip 1.70709pt:=\hskip 1.70709pt\{e(Q)=^{?}k\}\quad\text{and}\quad\text{\sc C}_{e}\hskip 1.70709pt:=\hskip 1.70709pt\{e(P)=^{?}e(Q)\},

the verification and coincidence problems for the numbers of linear extensions. Note that a version of (4.3) still holds in this case:

(5.3) PH P#PNPVeNP{e(P)=?e(Q)},{e(Q)=?k}NPCe,\textup{{{PH}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{P}}}^{\textup{{{\#P}}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc V}_{e}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\{e(P)=^{?}e(Q)\}\rangle,\hskip 1.70709pt\langle\{e(Q)=^{?}k\}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc C}_{e}\rangle},

where the first inclusion is Toda’s theorem again, in the third inclusion we guess both QQ and kk, and in the fourth inclusion we use  {e(Q)=?k}P\{e(Q)=^{?}k\}\in\textup{{{P}}}. From this point, proceed verbatim the proof of Theorem 4.2.

For  (5)(5), we use the proof of Proposition 5.20 to guess a forest  FF  with  n=O(logkloglogk)n=O(\log k\log\log k)  vertices and  ma(F)=k{\rm ma}(F)=k  matchings. Since we are using the approach in [KS21] again, we similarly conclude that this cannot be used to construct an explicit forest  FF  in polynomial time. On the other hand, note from the proof of the proposition, that give the forest  FF  and the order of removed vertices, the function  ma(F){\rm ma}(F)  can be computed by induction in polynomial time.

We now follow the approach in  (4)(4). Proceed by guessing the desired forest  FF  and the order of vertices to be removed, verify in polynomial time that this forest has exactly  kk  matchings, and obtain the analogue of (5.3). Recall also that the problem of computing the number of matchings is #P-complete for bipartite planar graphs (see above). From this point, proceed as in  (4)(4)  above.

For  (6)(6), recall that computing  b(M)b(M)  is  #P-complete [Sno12]. Recall also that graphic matroids are rational, so  the number of spanning trees  st{\rm st}  is a restriction of the counting function bb. We need only a weaker version [Sto22, Cor. 6.2.1], which states that for all  kk  sufficiently large, there is a simple graph  GG  on nn vertices with exactly  kk  spanning trees, and  n<C(logk)2n<C\hskip 0.85355pt(\log k)^{2}. Examining the proof of this result gives a polynomial time construction, but even without a careful examination the verification problem  {st(G)=?k}\{{\rm st}(G)=^{?}k\}  is clearly in  P  by the matrix-tree theorem. Thus we obtain the analogue of (5.3), and can proceed as above. ∎

Remark 5.31.

The authors of [KS21] reported to us555Personal communication (July, 2023). that there is a way to avoid integer factoring by separating primes into small:  pC(logk)3/2p\leq C(\log k)^{3/2}, and large:  p>C(logk)3/2p>C(\log k)^{3/2}. The former can be found exhaustively, while the latter can treated probabilistically in one swoop. This gives a probabilistic polynomial time algorithm for generating a poset with few elements and desired number of linear extensions. We should mention that the best deterministic version of this problem is known to give only n=O(k)n=O(\sqrt{k}), see [Ten09].

5.5. Most general version

Now that we treated many examples, we are ready to state the most general complexity result which can be used as a black box. For that we need a new definition which combines the notions of concise functions and the polynomial time complexity.

Let  𝒳=𝒳n\mathcal{X}=\sqcup\hskip 0.85355pt\mathcal{X}_{n}  be the set of combinatorial objects (see §\S5.1), and let  𝒴=𝒴n\mathcal{Y}=\sqcup\hskip 0.85355pt\mathcal{Y}_{n}  be a subset, such that  𝒴n𝒳n\mathcal{Y}_{n}\subseteq\mathcal{X}_{n} . For a counting function  f:𝒳f:\mathcal{X}\to\mathbb{N}  denote  𝒯f,𝒴={f(y)y𝒴}\mathcal{T}_{f,\mathcal{Y}}=\{f(y)\mid y\in\mathcal{Y}\}. Similarly, denote by

(5.4) Vf,𝒴:={f(y)=?ky𝒴nandkf(𝒳n)}\text{\sc V}_{f,\mathcal{Y}}\,:=\,\big{\{}\hskip 0.85355ptf(y)\hskip 0.85355pt=^{?}\hskip 0.85355ptk\hskip 1.70709pt\mid\hskip 1.70709pty\in\mathcal{Y}_{n}\ \ \text{and}\ \ k\in f(\mathcal{X}_{n})\hskip 0.85355pt\big{\}}

the restricted verification problem.

Definition 5.32.

A counting function  f:𝒳f:\mathcal{X}\to\mathbb{N}  is called recognizable  if there exists   𝒴=𝒴n,\mathcal{Y}=\sqcup\hskip 0.85355pt\mathcal{Y}_{n}\hskip 0.85355pt,𝒴n𝒳n\mathcal{Y}_{n}\subseteq\mathcal{X}_{n} , such that  𝒯f,𝒳=𝒯f,𝒴\mathcal{T}_{f,\mathcal{X}}=\mathcal{T}_{f,\mathcal{Y}}  and  Vf,𝒴PH\text{\sc V}_{f,\mathcal{Y}}\hskip 0.85355pt\in\hskip 0.85355pt\textup{{{PH}}}.

Note that if  ff  is almost complete on 𝒳\mathcal{X}, then it is almost complete on 𝒴\mathcal{Y}. In this case, we can replace the  “kf(𝒳n)k\in f(\mathcal{X}_{n})”  assumption with “for all sufficiently large kk”. In order for the problem to be in PH, the input of yy has to have a polynomial in the bit-length of kk. Thus, almost complete recognizable functions must also be concise.

Proposition 5.33.

Let  ff  be a counting function which is both #P-complete and recognizable. Then  CfPHPH=Σmp\in\textup{{{PH}}}\ \Rightarrow\ \textup{{{PH}}}=\Sigma^{{\textup{p}}}_{m}  for some mm.

Proof.

Since  ff is recognizable, we have  Vf,𝒴PH{}_{f,\mathcal{Y}}\in\textup{{{PH}}}  for some  𝒴𝒳\mathcal{Y}\subseteq\mathcal{X}. Then  Vf,𝒴Σp{}_{f,\mathcal{Y}}\in\Sigma^{{\textup{p}}}_{\ell}  for some  1\ell\geq 1. We have:

PH P#PNPVfNP{f(x)=?f(y)},{f(y)=?k}NPCf,Vf,𝒴NPCf,ΣpΣ+1Cf,\textup{{{PH}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{P}}}^{\textup{{{\#P}}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc V}_{f}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\{f(x)=^{?}f(y)\}\rangle,\hskip 1.70709pt\langle\{f(y)=^{?}k\}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc C}_{f}\rangle,\hskip 1.70709pt\langle\text{\sc V}_{f,\mathcal{Y}}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\textup{{{NP}}}^{\langle\text{\sc C}_{f}\rangle,\hskip 1.70709pt\Sigma^{{\textup{p}}}_{\ell}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\Sigma_{\ell+1}^{\langle\text{\sc C}_{f}\rangle},

where in the first inclusion we use Toda’s theorem again, in the third inclusion we guess both yy and kk, and in the fourth inclusion we use the definition of recognizable functions.

Now suppose  CfPH{}_{f}\in\textup{{{PH}}}. Then   CfΣrp{}_{f}\in\hskip 0.85355pt\Sigma^{{\textup{p}}}_{r}  for some rr. Then we have:

(5.5) PH Σ+1CfΣ+1ΣrpΣr++1p,\textup{{{PH}}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\Sigma_{\ell+1}^{\langle\text{\sc C}_{f}\rangle}\hskip 1.70709pt\subseteq\hskip 1.70709pt\Sigma_{\ell+1}^{\Sigma^{{\textup{p}}}_{r}}\hskip 1.70709pt\subseteq\hskip 1.70709pt\Sigma^{{\textup{p}}}_{r+\ell+1}\,,

as desired. ∎

Remark 5.34.

For example, suppose ff is concise and the instance  y𝒴ny\in\mathcal{Y}_{n}  with given  f(y)=kf(y)=k  can be obtained by a probabilistic polynomial time algorithm as in Remark 5.31. Since  BPPΣ2p\textup{{{BPP}}}\subseteq\Sigma^{{\textup{p}}}_{2}, Proposition 5.33 applies in this case.

Similarly, suppose Conjecture 5.17 has a nonconstructive proof, where the family of posets 𝒴\mathcal{Y} with some parameters are introduced, and the almost completeness is proved by a probabilistic method. Recall that  #LE  restricted to height two posets is #P-complete. By Proposition 5.16, the restriction of function  ee  to height two posets is concise. Therefore, if the number of random bits used is at most polynomial, we can again apply Proposition 5.33.

6. Final remarks and further open problems

6.1. On terminology

There is a gulf of a difference between the notions of “equality”, “equality conditions” and “coincidence”. Deciding equality  is typically of the form  {f(x)=?g(x)x}\{f(x)=^{?}g(x)\ \forall x\}, whether two functions are equal for all xx. Deciding the equality conditions  is typically of the form  {f(x)=?g(x)}\{f(x)=^{?}g(x)\}, whether two functions are equal for a given xx. Usually this refer to the inequality  f(x)g(x)f(x)\geqslant g(x)  which holds for all xx (cf. §\S6.19). Finally, deciding the coincidence  is typically of the form  {f(x)=?f(y)}\{f(x)=^{?}f(y)\}, whether the same function has equal values at givenx,yx,y. Despite superficial similarities, these properties have very distinctive flavors from the computational complexity point of view. We refer to [Wig19] for the general background and philosophy.

6.2. Two types of #P-complete problems

There are so many known #P-complete problems, it can be difficult to trace down the literature to see if they are proved by a parsimonious reduction from  #3SAT (often, they are not). Helpfully, [GJ79] addresses this for every reduction, but [Pak22, §\S13] does not, for example. Clearly, a parsimonious reduction from  #3SAT  is not possible if the corresponding decision problem is in P. This applies to the numbers of linear extensions, order ideals, independent sets, matchings and bases of matroids considered in Section 5, since the corresponding decision problems are trivially in P.

6.3. Bases in matroids

Recall that every graphic matroid is binary, so binary matroids  is a more natural class to be used in part (6)(6) of Main Theorem 1.4. For the number of bases of binary matroids, the #P-completeness result was claimed by Vertigan and is widely quoted (see e.g. [Wel93]), despite not appearing in print (cf. [Sno12] and [Pak22, §\S13.5.8]). The complete proof was obtained most recently by Knapp and Noble [KN23, Thm 50], after this paper was already written.

Let us mention that number of bases of matroids is  #P-complete for other concise presentations. Notable classes include sparse paving matroids [Jer06] (proved via a parsimonious reduction to the number of Hamiltonian cycles), and bicircular matroids [GN06] (proved via a non-parsimonious reduction to the permanent).

Conjecture 6.1.

Function  bb  restricted to bicircular matroids is concise.

6.4. Linear extensions of height two posets

In the context of Conjecture 5.17, denote by  ee^{\prime}  the restriction of function  ee^{\prime}  to posets height two. The conjecture claims that ee^{\prime} is almost complete, even though the numerical evidence points in the opposite direction:

(6.1) 7,9,10,11,13,15,17,18,19,21,22,23,26,27,28,29,31,32,𝒯e.7,\hskip 1.70709pt9,\hskip 1.70709pt10,\hskip 1.70709pt11,\hskip 1.70709pt13,\hskip 1.70709pt15,\hskip 1.70709pt17,\hskip 1.70709pt18,\hskip 1.70709pt19,\hskip 1.70709pt21,\hskip 1.70709pt22,\hskip 1.70709pt23,\hskip 1.70709pt26,\hskip 1.70709pt27,\hskip 1.70709pt28,\hskip 1.70709pt29,\hskip 1.70709pt31,\hskip 1.70709pt32,\,\ldots\,\notin\hskip 1.70709pt\mathcal{T}_{e^{\prime}}\hskip 1.70709pt.

However, the number of height two posets on nn elements is asymptotically  expΘ(n2)\exp\Theta(n^{2}), thus overwhelming the numbers of linear extensions (cf. Example 5.23). In fact, almost all posets have bounded height [KR75]. Additionally,  𝒯e\mathcal{T}_{e^{\prime}}  can contain large primes, e.g.  E6=61𝒯e(6)E_{6}=61\in\mathcal{T}_{e^{\prime}}(6)  and

E38= 23489580527043108252017828576198947741𝒯e(38),E_{38}\,=\,23489580527043108252017828576198947741\,\in\hskip 1.70709pt\mathcal{T}_{e^{\prime}}(38),

see [OEIS, A092838]. Here  En:=e(Zn)E_{n}:=e(Z_{n})  is the Euler number (also called secant number), defined as the number of linear extensions of the zigzag posetx1x2x3x4x_{1}\prec x_{2}\succ x_{3}\prec x_{4}\succ\ldots\hskip 1.70709pt, see e.g. [Sta10] and [OEIS, A000111]. Clearly, ZnZ_{n} has height two. Thus we believe that positive evidence outweighs the negative, favoring the conjecture.666Extensive computer computations by David Soukup (personal communication, August 2023), show that  1593350922239999𝒯e1593350922239999\notin\hskip 0.85355pt\mathcal{T}_{e^{\prime}}  suggesting that (6.1) might not be an instance of the strong law of small numbers [Guy88].

6.5. Hamiltonian paths in tournaments

An orientation of all edges of a complete graph  KnK_{n}  is called a tournament. Denote by  h(T)h(T)  the number of (directed) Hamiltonian paths in TT. Following Szele (1943), there is large literature on the maximal value of hh over all tournaments with nn vertices, see an overview in [KO12, §\S4.2]. Rédei (1934) proved that  h(T)h(T)  always odd, see e.g. [Moon68, §\S9] and [Ber85, §10.2\S 10.2]. Curiously, it is not known if  (h1)/2(h-1)/2  is almost complete:

Conjecture 6.2.

Function h\hskip 0.85355pth  takes all odd values except 77 and 2121.

This conjecture was made by the MathOverflow user bof and Gordon Royle (March 2016).777mathoverflow.net/q/232751 Royle reported that the conjecture holds for all odd integers up to 80557 (ibid).

6.6. Degree sequences

Let  d=(d1,,dn)n\textbf{{d}}=(d_{1},\ldots,d_{n})\in\mathbb{N}^{n}, such that  0d1,,dnn10\leq d_{1},\ldots,d_{n}\leq n-1. Denote by  c(d)c(\textbf{{d}})  the number of simple graphs on nn vertices  V={v1,,vn}V=\{v_{1},\ldots,v_{n}\}, with  deg(vi)=di\deg(v_{i})=d_{i}  for all  1in1\leq i\leq n. It was conjectured in [Pak22, §\S13.5.4], that computing c()\hskip 0.85355ptc(\cdot)  is #P-complete. Note that the decision problem  {c(d)>?0}\{c(\textbf{{d}})>^{?}0\}  is in P by the Erdős–Gallai theorem, see [EG61, SH91]. We refer to [Wor18] for a recent survey of the large literature on the subject.

Conjecture 6.3.

Function c\hskip 0.85355ptc  is concise, i.e. for every  k>0k>0  there exists  d=(d1,,dn)n\textbf{{d}}=(d_{1},\ldots,d_{n})\in\mathbb{N}^{n}, such that  c(d)=kc(\textbf{{d}})=k  and  nC(logk)cn\leq C(\log k)^{c}, for some fixed  C,c>0C,c>0.

Note that graphs with a given degree sequences can also be viewed as binary (0/1)(0/1) symmetric contingency tables, making this conjecture a variation on Conjecture 5.30.

6.7. Plane triangulations

Let  X={p1,,pn}2X=\{p_{1},\ldots,p_{n}\}\in\mathbb{Q}^{2}  be the set of nn points in general positions, i.e. where no three points lie on a line. Denote by  t(X)t(X)  the number of triangulations of the convex hull of  XX  which contain all vertices in XX. It was conjectured in [DRS10, Exc. 8.16], that computing  tt  is  #P-complete (see also [Epp20] for a closely related problem). It is known that  tt  has exponential growth (see e.g. [DRS10, §\S3.3]):

C12.3nt(X)C243nfor some C1,C2>0.C_{1}\cdot 2.3^{n}\hskip 1.70709pt\leq\hskip 1.70709ptt(X)\hskip 1.70709pt\leq\hskip 1.70709ptC_{2}\cdot 43^{n}\quad\text{for some \ $C_{1},C_{2}>0$.}
Conjecture 6.4.

Function t\hskip 0.85355ptt  is almost complete.

By Proposition 5.14, the conjecture implies that  tt is concise. One reason to believe the conjecture is the exponential number of topological triangulations on nn vertices, given by Tutte’s formula (see e.g. [DRS10, Thm 3.3.5]), combined with Fáry’s theorem  that all topological triangulations are realizable as plane triangulations (see e.g. [PA95, Thm 8.2]). Thus, there is no growth discrepancy as in the proof of Proposition 5.26.

6.8. Trees in planar graphs

Let  G=(V,E)G=(V,E)  be a simple planar graph, and let  tr(G)tr(G)  be the number of trees in GG (of all sizes). Trivially, the number of trees in an empty graph with nn vertices is nn, so  trtr  is almost complete. On the other hand, Jerrum showed that computing  trtr  is #P-complete [Jer94], but the proof uses a non-parsimonious reduction to  #2SAT. This suggests the following:

Conjecture 6.5.

Function tr\hskip 0.85355pttr  is concise.

Note that the number of planar graphs on nn vertices is asymptotically  Cn7/2γnC\hskip 0.85355ptn^{-7/2}\hskip 0.85355pt\gamma^{n}  where  γ27.23\gamma\approx 27.23, see [GN09], which is large enough to make the conjecture plausible.

6.9. Critical group

In [GK20, p. 19], Glass and Kaplan ask about the smallest number of vertices  n=n(H)n=n({\text{\rm H}})  of a graph with a given critical groupH (also known as the sandpile group), that can be defined by the Smith normal form  of the graph Laplacian. Recall that for every graph GG, the size of its critical group is the number of spanning trees in GG :  |HG|=st(G)|{\text{\rm H}}_{G}|={\rm st}(G). Thus this problem refines the problem in the Example 5.21. In particular, Stong’s theorem [Sto22] implies that  n(H)=o((logk)3/2)n({\text{\rm H}})=o((\log k)^{3/2})  for all square-freek=|H|k=|{\text{\rm H}}|, and suggests that the same bound holds for all H. We refer to [CP18, Ch. 6] and [Kli19, Ch. 4] for the background and further references on the critical group.

6.10. Determinants

In their study of hypergraphs with positive discrepancy, Cherkashin and Petrov [CP19] considered the following problem. Fix a parameter  δ\delta\in\mathbb{N}. Denote by  𝒯δ(n)\mathcal{T}_{\delta}(n)  the set of all possible determinants of  n×nn\times n  matrices with entries in  {0,1,,δ}\{0,1,\ldots,\delta\}. They show that

𝒯4(n){0,1,,cn}for somec>1.\mathcal{T}_{4}(n)\,\supset\,\big{\{}\hskip 0.85355pt0,1,\ldots,c^{n}\hskip 0.85355pt\big{\}}\quad\text{for some}\ \ c>1.

This result shows that the determinant is a concise  GapP  function. A stronger result is proved by Shah [Shah22] in connection with random binary (0/1)(0/1) matrices:

𝒯1(n){0,1,,α2nn}for someα>0,\mathcal{T}_{1}(n)\,\supset\,\big{\{}\hskip 0.85355pt0,1,\ldots,\hskip 0.85355pt\tfrac{\alpha\hskip 0.85355pt2^{n}}{n}\hskip 1.70709pt\big{\}}\quad\text{for some}\ \ \alpha>0,

(cf. [OEIS, A013588]).

In [CP19, §\S4], the authors ask whether these results can be strengthened to

(6.2) 𝒯δ(n){0,1,,cnlogn}for someδ1andc>1.\mathcal{T}_{\delta}(n)\,\supset\,\big{\{}\hskip 0.85355pt0,1,\ldots,\hskip 0.85355ptc^{n\log n}\hskip 0.85355pt\big{\}}\quad\text{for some}\ \ \delta\geq 1\ \text{and}\ c>1.

Denote  βn:=max𝒯1(n)\beta_{n}:=\max\mathcal{T}_{1}(n). Finding  βn\beta_{n}  is the famous Hadamard maximal determinant problem, see e.g.  [BC72, BEHO21], [OEIS, A003432] and references therein. It was shown by Hadamard (1893) that  βnnn/2\beta_{n}\leq n^{n/2}. In a different direction, it is known that  βn>Cenlogncn\beta_{n}>C\hskip 0.85355pte^{n\log n-cn}  for an infinite set of integer nn and fixed c,C>0c,C>0, see [BEHO21, Cor. 27]. This lends a (weak) partial evidence towards (6.2).

Finally, recall that by the matrix-tree theorem, the number of spanning trees is a determinant of the Laplacian. Thus, since the graphs in [Sto22] have degrees at most six (cf. Example 5.21), Stong’s result gives a (6.2)–style inclusion for matrices with entries in  {1,0,1,,6}\{-1,0,1,\ldots,6\}, but with a weaker upper limit  cn2/3c^{n^{2/3}}. By contrast, Azarija and Škrekovski conjecture in [AS13] that one can take the upper bound to be  eω(n)e^{\omega(n)}  for general simple graphs. Of course, this bound cannot be achieved on graphs with bounded (average) degree, since they have at most exponential number of spanning trees, see e.g. [Gri76].

6.11. Reduced factorizations

For a permutation  σSn\sigma\in S_{n}, a reduced factorization  is a product  σ=(i1,i1+1)(i,i+1)\sigma=(i_{1},i_{1}+1)\cdots(i_{\ell},i_{\ell}+1)  of  =inv(σ)\ell=\operatorname{{\rm inv}}(\sigma)  adjacent transpositions. Denote by  red(σ)\mathrm{red}(\sigma)  the number of reduced factorizations of σ\sigma. It was conjectured in [DP18] that computing  red\mathrm{red}  is #P-complete. Observe that function  red\mathrm{red}  is almost complete, since

red(2,1,n,3,4,,n1)=n2forn3.\mathrm{red}(2,1,n,3,4,\ldots,n-1)\,=\,n-2\quad\,\text{for}\ \ \hskip 1.70709ptn\hskip 0.85355pt\geq\hskip 0.85355pt3.
Conjecture 6.6.

Function  red\mathrm{red}  is concise.

Let  σSn\sigma\in S_{n}  be a 321321-avoiding  permutation, i.e. suppose that  PC321(σ)=0\mathrm{PC}_{321}(\sigma)=0. We have in this case:  red(σ)=SYT(λ/μ)\mathrm{red}(\sigma)=\operatorname{{\rm SYT}}(\lambda/\mu), where  λ/μ\lambda/\mu  is a skew shape associated to σ\sigma, see [BJS93, pp. 358-359]. Vice versa, for every skew shape  λ/μ\lambda/\mu  of size nn, there is a 321321-avoiding permutation  σSn\sigma\in S_{n}  such that  red(σ)=SYT(λ/μ)\mathrm{red}(\sigma)=\operatorname{{\rm SYT}}(\lambda/\mu). In other words, a positive answer to Question 5.24 implies Conjecture 6.6. For further discussions and applications of this result, see e.g. [Man01, §\S2.6.6] and [MPP19, §\S2.2, §\S5.3].

6.12. Littlewood–Richardson coefficients

Recall that the Littlewood–Richardson coefficients is a complete function (see Example 5.3), and that the maximal LR-coefficient is exponential in nn, see [PPY19]. This suggest the following:

Conjecture 6.7.

The Littlewood–Richardson coefficients  cμνλc^{\lambda}_{\mu\nu}  is a concise function, i.e. for every  k>0k>0  there exists three partitions  λ,μ,ν\lambda,\mu,\nu, such that  cμνλ=kc^{\lambda}_{\mu\nu}=k,  |λ|=|μ|+|ν|=n|\lambda|=|\mu|+|\nu|=n  and  nC(logk)cn\leq C(\log k)^{c}, for some fixed  C,c>0C,c>0.

Using standard reductions, this conjecture follows from Conjecture 5.30 that the numbers of contingency tables is a concise function, see e.g. [PV10]. Let us mention that similarly to the domino tilings, the LR-coefficients have a combinatorial interpretation as the number of Knutson–Tao puzzles  with boundary markings as shown in Figure 6.1, see e.g. [Knu23].

Refer to caption
Figure 6.1. Tiles of the Knutson–Tao puzzles (up to rotation, reflections are not allowed), and an example of a puzzle.

6.13. Integer points

Let  PdP\in\mathbb{R}^{d}  be a rational convex H-polytope, i.e. defined by inequalities over \mathbb{Q}, and let  ζ(P)\zeta(P)  be the number of integer points in PP. One can ask if  ζ\zeta  is concise in full generality? Here we are assuming that the input is in unary. Note that  ζ\zeta  is a counting function which generalizes the number of contingency tables (Example 5.4) and LR-coefficients (Example 5.3), but we don’t know if either of these is concise (see Conjectures 5.30 and 6.7).

There is a sneaky way to establish that  ζ\zeta  is concise, by noting that in a special case this becomes one of the known  #P-complete problems in unary, with a parsimonious reduction from the  0/1 PERMANENT (see e.g. [DO04, Thm 1.2])  or the  #3SAT (see e.g. [IMW17, §\S5.1] and references therein). Recall that both of these problems have concise counting functions. In both cases, Proposition 5.12 implies that  ζ\zeta  is concise. We omit the details.

6.14. Domino tilings of simply connected regions

Recall that our proof of Theorem 3.1 uses regions that are not simply connected. Note also that Nadeau’s construction in Remark 3.2 is simply connected.

Conjecture 6.8.

Theorem 3.1 holds for simply connected regions.

If true, the proof would have to be much more elaborate than our proof for general regions. Indeed, note that domino tilings of simply connected regions have additional properties given by the height functions  (see e.g. [PST16, Rém04, Thu90]).

Additionally, the generalized Temperley’s bijection [KPW00] relates the number of domino tilings in simply connected regions with the number of spanning trees in certain grid graphs. Thus, Conjecture 6.8 can be reformulated to say that a certain restriction of  st{\rm st}  is concise, see Example 5.21. This gives another reason to think that the conjecture might be difficult.

On the other hand, it follows from the approach [ÇS18, Sch19] and [CP24] that Conjecture 6.8 follows from Zaremba’s conjecture. Moreover, for the snake regions defined in [ÇS18] this approach gives the following analogue of (5.1); we omit the details.

Theorem 6.9.

There is  c>1c>1, such that  𝒯snake(n){0,1,,cn/(logn)}\mathcal{T}_{\text{\rm snake}}(n)\supseteq\big{\{}0,1,\ldots,c^{n/(\log n)}\big{\}}, for all  n1n\geq 1.

6.15. Tilings with small tiles

The literature on tilings in the plane is much too large to be reviewed here, but let us mention a few relevant highlight. In [MR01], the authors proved NP-completeness for the decision problem and #P-completeness for the counting problem of tilings of general regions for a small set with just two tiles (up to rotation). In [BNRR95],  NP-completeness was proved for tilings with barsk×1k\times 1  and  1×k1\times k, where  k3k\geq 3. This is especially remarkable since for the same set of tiles, tileability of simply connected regions is in P, see [KK93]. Finally, in [PY13b], a finite set of rectangle tiles was given, with an NP-complete decision problem and #P-complete corresponding counting problem for simply connected regions.

6.16. The induction two-step

The inductive arguments in the paper can be used to extend some results of concise functions. This phenomenon can be seen, for example, in Corollary 3.3 which arises naturally from the proof of Theorem 3.1. Similarly, for relatively prime  (a,b)(a,b), the analogue of the corollary holds for linear extensions of width two posets (Example 5.15), and the number of matchings of trees (Example 5.19). This is a corollary of properties of the Stern–Brocot tree (see references in [KS21]). Another notable example of this phenomenon is a beautiful result in [HKNS10], that for every two groups  Γ0\Gamma_{0}  and  Γ1\Gamma_{1}, there is a graph  G=(V,E)G=(V,E)  and an edge  eEe\in E, such that  Aut(G)=Γ0(G)=\Gamma_{0}  and  Aut(Ge)=Γ1(G-e)=\Gamma_{1} .

6.17. Concise functions

As we mentioned earlier, the notion of a concise function is closely related to  #P-completeness. Indeed, for every parsimonious reduction  fgf\mapsto g, if  ff  is concise, then so is gg. Below we give an example of a #P-complete function that is not  concise. Note that if either Conjecture 5.17 or Conjecture 6.1 is false, this would also give such examples.

On the other hand, note that there are several interesting examples of concise functions that are in FP, such as the number of spanning trees in simple graphs and the number of linear extensions of width two posets (Examples 5.21 and 5.15). This suggests that being concise is an interesting property in its own right, independently of complexity considerations.

6.18. Coincidence problems

Note that none of the natural #P-complete counting functions ff that we consider have the corresponding coincidence problems  Cf  in PH (unless PH collapses). Theorems 1.3 and 1.4 prove this in many cases, and some interesting examples remain open (see above). The same applies for the verification problem  Vf. Note that if  ff  is recognizable, then these two problems are PH-equivalent, i.e.  CfPHVfPH\text{\sc C}_{f}\in\textup{{{PH}}}\ \Longleftrightarrow\ \text{\sc V}_{f}\in\textup{{{PH}}}  by the argument in the proof of Proposition 5.33.

The following elegant example by Greta Panova shows that #P-complete functions can have easy coincidence problems.888Personal communication (August, 2023) In notation of §\S4.4, let  A=(aij)nA=(a_{ij})\in\mathcal{M}_{n}  where n\mathcal{M}_{n}  is the set of n×nn\times n  matrices with entries in  {0,1}\{0,1\}. Denote  g(A):=a11+2a12+4a13++2n21ann+2n2g(A):=a_{11}+2a_{12}+4a_{13}+\ldots+2^{n^{2}-1}a_{nn}+2^{n^{2}}  and let  f:nf:\sqcup\hskip 0.85355pt\mathcal{M}_{n}\to\mathbb{N}  be defined by  f(A):=per(A)+2n2g(A)f(A):=\mathrm{per}(A)+2^{n^{2}}g(A). Observe that  per(A)n!<2n2\mathrm{per}(A)\leq n!<2^{n^{2}}. Since  f(A)=per(A)mod2n2f(A)=\mathrm{per}(A)\mod 2^{n^{2}}, we conclude that computing f\hskip 0.85355ptf\hskip 0.85355pt is  #P-complete. On the other hand, we have  f(A)=f(A)f(A)=f(A^{\prime})  if and only if  A=AA=A^{\prime}, so  CfP{}_{f}\in\textup{{{P}}}.

Curiously, function  ff  is not  almost complete since the first  (n2+1)(n^{2}+1)  bits determine the last  n2n^{2}  bits, but is concise since  n2=O(logf(A))n^{2}=O\big{(}\log f(A)\big{)}  for all  AnA\in\mathcal{M}_{n} . Function  ff  is not recognizable, however, since the corresponding restricted verification problem is exactly  VPER{}_{\text{PER}}  which is not in PH (unless PH collapses). Now take  n:=n×[n]\mathcal{M}_{n}^{\prime}:=\mathcal{M}_{n}\times[n]. Define  f(A,):=f(A)f^{\prime}(A,\ell):=f(A)  if  =0\ell=0, and  f(A,):=1f^{\prime}(A,\ell):=\ell-1  if  >0\ell>0. Then function  ff^{\prime}  is complete, and thus not  concise. Of course, computing  ff^{\prime}  remains  #P-complete.

This leaves us with more questions than answers. Is  CPERcoNP{}_{\text{PER}}\in\textup{{{coNP}}}-hard? Is it  C=P-complete under Turing reductions? What about other natural coincidence problems that we discuss in this paper? Does it make sense to consider a complexity class of  #P-complete recognizable functions as in Proposition 5.33?

6.19. Equality conditions

The equality conditions for many combinatorial inequalities, have been studied widely in the area, especially in order theory. We refer to [Gra83, Win86] for dated surveys, to [Pak22] for a recent overview. Computationally, the equality conditions are decision problems in  C=P. This paper grew out of our efforts to understand the Alexandrov–Fenchel inequality for mixed volumes [CP23].


Acknowledgements

We are grateful to Karim Adiprasito, Sasha Barvinok, Richard Brualdi, Jesús De Loera, Darij Grinberg, Christian Ikenmeyer, Mark Jerrum, Nathan Kaplan, Greg Kuperberg, Alejandro Morales, Steven Noble, Greta Panova, Fedya Petrov, Colleen Robichaux, Rikhav Shah, David Soukup, Richard Stanley and Bridget Tenner, for useful discussions and helpful remarks. Special thanks to Noah Kravitz, Ashwin Sah and Richard Stong for their help explaining [KS21] and [Sto22].

This paper was finished when both authors were visiting the American Institute of Mathematics at their new location in Pasadena, CA. We are grateful to AIM for the hospitality. The first author was partially supported by the Simons Foundation. Both authors were partially supported by the NSF.


References

  • [AB09] Sanjeev Arora and Boaz Barak, Computational complexity. A modern approach, Cambridge Univ. Press, Cambridge, 2009, 579 pp.
  • [AS13] Jernej Azarija and Riste Škrekovski, Euler’s idoneal numbers and an inequality concerning minimal graphs with a prescribed number of spanning trees, Math. Bohem. 138 (2013), 121–131.
  • [Bar17] Alexander Barvinok, Counting integer points in higher-dimensional polytopes, in Convexity and concentration, Springer, New York, 2017, 585–612.
  • [BNRR95] Danièle Beauquier, Maurice Nivat, Éric Remila and Mike Robson, Tiling figures of the plane with two bars, Comput. Geom. 5 (1995), 1–25.
  • [Ber85] Claude Berge, Graphs (second ed.), North-Holland, Amsterdam, 1985, 413 pp.
  • [BJS93] Sara C. Billey, William Jockusch and Richard P. Stanley, Some combinatorial properties of Schubert polynomials, J. Algebraic Combin. 2 (1993), 345–374.
  • [BBL98] Prosenjit Bose, Jonathan F. Buss and Anna Lubiw, Pattern matching for permutations, Inform. Process. Lett. 65 (1998), 277–283.
  • [BLP23] Petter Brändén, Jonathan Leake and Igor Pak, Lower bounds for contingency tables via Lorentzian polynomials, Israel J. Math. 253 (2023), 43–90.
  • [BC72] Joel Brenner and Larry Cummings, The Hadamard maximum determinant problem, Amer. Math. Monthly 79 (1972), 626–630.
  • [BW91] Graham Brightwell and Peter Winkler, Counting linear extensions, Order 8 (1991), 225–247.
  • [BEHO21] Patrick Browne, Ronan Egan, and Fintan Hegarty and Padraig Ó Catháin, A survey of the Hadamard maximal determinant problem, Electron. J. Combin. 28 (2021), no. 4, Paper No. 4.41, 35 pp.
  • [BN65] Richard A. Brualdi and Morris Newman, Some theorems on the permanent, J. Res. Nat. Bur. Standards, Sect. B 69B (1965), 159–163.
  • [ÇS18]  Ilke Çanakçı  and Ralf Schiffler, Cluster algebras and continued fractions, Compos. Math. 154 (2018), 565–593.
  • [CP23] Swee Hong Chan and Igor Pak, Equality cases of the Alexandrov–Fenchel inequality are not in the polynomial hierarchy, preprint (2023), 35 pp.;  arXiv:2309.05764; extended abstract to appear in Proc. 56th STOC, ACM, New York, 2024.
  • [CP24] Swee Hong Chan and Igor Pak, Linear extensions and continued fractions, preprint (2024), 14 pp.;  arXiv:2401.09723.
  • [CP19] Danila Cherkashin and Fëdor Petrov, On small nn-uniform hypergraphs with positive discrepancy, J. Combin. Theory, Ser. B 139 (2019), 353–359.
  • [CP18] Scott Corry and David Perkinson, Divisors and sandpiles. An introduction to chip-firing, AMS, Providence, RI, 2018, 325 pp.
  • [Cra08] David A. Craven, Symmetric group character degrees and hook numbers, Proc. LMS 96 (2008), 26–50.
  • [DL92] Paul Dagum and Michael Luby, Approximating the permanent of graphs with large factors, Theoret. Comput. Sci. 102 (1992), 283–305.
  • [DO04] Jesús A. De Loera and Shmuel Onn, The complexity of three-way statistical tables, SIAM J. Comput. 33 (2004), 819–836.
  • [DRS10] Jesús A. De Loera, Jörg Rambau and Francisco Santos, Triangulations, Springer, Berlin, 2010, 535 pp.
  • [DP18] Samuel Dittmer and Igor Pak, Counting linear extensions of restricted posets, preprint (2018), 34 pp.; arXiv:1802.06312.
  • [Epp20] David Eppstein, Counting polygon triangulations is hard, Discrete Comput. Geom. 64 (2020), 1210–1234.
  • [EG61] Paul Erdős and Tibor Gallai, Graphs with prescribed degrees of vertices (in Hungarian), Mat. Lapok 11 (1961), 264–274.
  • [GJ78] Michael R. Garey and David S. Johnson, “Strong” NP-completeness results: motivation, examples, and implications, J. ACM 25 (1978), 499–508.
  • [GJ79] Michael R. Garey and David S. Johnson, Computers and Intractability: A Guide to the Theory of NP-completeness, Freeman, San Francisco, CA, 1979, 338 pp.
  • [For97] Lance Fortnow, Counting complexity, in Complexity theory retrospective II, Springer, New York, 1997, 81–107.
  • [GN06] Omer Giménez and Marc Noy, On the complexity of computing the Tutte polynomial of bicircular matroids, Combin. Probab. Comput. 15 (2006), 385–395.
  • [GN09] Omer Giménez and Marc Noy, Asymptotic enumeration and limit laws of planar graphs, J. Amer. Math. Soc. 22 (2009), 309–329.
  • [GK20] Darren Glass and Nathan Kaplan, Chip-firing games and critical groups, in A project-based guide to undergraduate research in mathematics, Springer, Cham, 2020, 107–152.
  • [Gol08] Oded Goldreich, Computational complexity. A conceptual perspective, Cambridge Univ. Press, Cambridge, UK, 2008, 606 pp.
  • [Gra83] Ronald L. Graham, Applications of the FKG inequality and its relatives, in Mathematical programming: the state of the art, Springer, Berlin, 1983, 115–131.
  • [Gre93] Frederic Green, On the power of deterministic reductions to  C=P, Math. Systems Theory 26 (1993), 215–233.
  • [Gri14] Rostislav Grigorchuk, Milnor’s problem on the growth of groups and its consequences, in Frontiers in complex dynamics, Princeton Univ. Press, Princeton, NJ, 2014, 705–773.
  • [Gri76] Geoffrey R. Grimmett, An upper bound for the number of spanning trees of a graph, Discrete Math. 16 (1976), 323–324.
  • [Gro93] Helmut Groemer, Stability of geometric inequalities, in Handbook of convex geometry, Vol. A, North-Holland, Amsterdam, 1993, 125–150.
  • [GT18] Alexander E. Guterman and Konstantin A. Taranin, On the values of the permanent of (0,1)(0,1)-matrices, Linear Algebra Appl. 552 (2018), 256–276.
  • [Guy88] Richard K. Guy, The strong law of small numbers, Amer. Math. Monthly 95 (1988), 697–712.
  • [HKNS10] Stephen G. Hartke, Hannah Kolb, Jared Nishikawa and Derrick Stolee, Automorphism groups of a graph and a vertex-deleted subgraph, Electron. J. Combin. 17 (2010), no. 1, RP 134, 8 pp.
  • [IMW17] Christian Ikenmeyer, Ketan D. Mulmuley and Michael Walter, On vanishing of Kronecker coefficients, Comp. Complexity 26 (2017), 949–992.
  • [IP22] Christian Ikenmeyer and Igor Pak, What is in #P and what is not?, preprint (2022), 82 pp.; extended abstract in Proc. 63rd FOCS (2022), 860–871;  arXiv:2204.13149.
  • [IPP22] Christian Ikenmeyer, Igor Pak and Greta Panova, Positivity of the symmetric group characters is as hard as the polynomial time hierarchy, preprint (2022), 15 pp.; extended abstract in Proc. 34th SODA (2023), 3573–3586;  arXiv:2207.05423.
  • [Jer87] Mark Jerrum, Two-dimensional monomer-dimer systems are computationally intractable. J. Statist. Phys. 48 (1987), 121–134. Erratum, 59 (1990), 1087–1088.
  • [Jer94] Mark Jerrum, Counting trees in a graph is #P-complete, Inform. Process. Lett. 51 (1994), 111–116.
  • [Jer06] Mark Jerrum, Two remarks concerning balanced matroids, Combinatorica 26 (2006), 733–742.
  • [KK93] Claire Kenyon and Richard Kenyon, Comment paver des polygones avec des barres (in French), C. R. Acad. Sci. Paris, Sér. I Math. 316 (1993), no. 12, 1319–1322;  Extended abstract:  Tiling a polygon with rectangles, in Proc. 33rd FOCS (1992), 610–619.
  • [Ken04] Richard Kenyon, An introduction to the dimer model, in ICTP Lect. Notes XVII, Trieste, 2004, 267–304.
  • [KPW00] Richard Kenyon, James G. Propp, and David B.  Wilson, Trees and matchings, Electron. J. Combin. 7 (2000), RP 25, 34 pp.
  • [Kit11] Sergey Kitaev, Patterns in permutations and words, Springer, Heidelberg, 2011, 494 pp.
  • [KR75] Daniel J. Kleitman and Bruce L. Rothschild, Asymptotic enumeration of partial orders on a finite set, Trans. AMS 205 (1975), 205–220.
  • [Kli19] Caroline J. Klivans, The mathematics of chip-firing, CRC Press, Boca Raton, FL, 2019, 295 pp.
  • [KN23] Christopher Knapp and Steven Noble, The complexity of the greedoid Tutte polynomial, preprint (2023), 44 pp.;  arXiv:2309.04537.
  • [Knu98] Donald E. Knuth, The art of computer programming. Vol. 2. Seminumerical algorithms (Third ed.), Addison-Wesley, Reading, MA, 1998, 762 pp.
  • [Knu23] Allen Knutson, Schubert calculus and quiver varieties, Proc. ICM (2022, virtual), Vol. VI, 4582–4605, EMS Press, Berlin, 2023; available at  tinyurl.com/5n99dbja
  • [KS21] Noah Kravitz and Ashwin Sah, Linear extension numbers of nn-element posets, Order 38 (2021), 49–66.
  • [KO12] Daniela Kühn and Deryk Osthus, A survey on Hamilton cycles in directed graphs, European J. Combin. 33 (2012), 750–766.
  • [LP09] László Lovász and Michael D. Plummer, Matching theory, AMS, Providence, RI, 2009.
  • [Mac95] Ian G. Macdonald, Symmetric functions and Hall polynomials (Second ed.), Oxford U. Press, New York, 1995, 475 pp.
  • [Man01] Laurent Manivel, Symmetric functions, Schubert polynomials and degeneracy loci, SMF/AMS, Providence, RI, 2001, 167 pp.
  • [Moon68] John W. Moon, Topics on tournaments, Holt, Rinehart and Winston, New York, 1968, 104 pp.; available at  www.gutenberg.org/ebooks/42833
  • [MM11] Cristopher Moore and Stephan Mertens, The nature of computation, Oxford Univ. Press, Oxford, UK, 2011, 985 pp.
  • [MR01] Cristopher Moore and John M. Robson, Hard tiling problems with simple tiles, Discrete Comput. Geom. 26 (2001), 573–590.
  • [MPP19] Alejandro H. Morales, Igor Pak and Greta Panova, Hook formulas for skew shapes III. Multivariate and product formulas, Algebraic Comb. 2 (2019), 815–861.
  • [PA95] János Pach and Pankaj K. Agarwal, Combinatorial geometry, John Wiley, New York, 1995, 354 pp.
  • [Pak22] Igor Pak, What is a combinatorial interpretation?, in Open Problems in Algebraic Combinatorics, AMS, Providence, RI, to appear, 58 pp.;  arXiv:2209.06142.
  • [PP20] Igor Pak and Greta Panova, Bounds on Kronecker coefficients via contingency tables, Linear Algebra Appl. 602 (2020), 157–178.
  • [PP22] Igor Pak and Greta Panova, Durfee squares, symmetric partitions and bounds on Kronecker coefficients, J. Algebra 629 (2023), 358–380.
  • [PPY19] Igor Pak, Greta Panova and Damir Yeliussizov, On the largest Kronecker and Littlewood–Richardson coefficients, J. Combin. Theory, Ser. A 165 (2019), 44–77.
  • [PST16] Igor Pak, Adam Sheffer and Martin Tassy, Fast domino tileability, Discrete Comput. Geom. 56 (2016), 377–394.
  • [PV10] Igor Pak and Ernesto Vallejo, Reductions of Young tableau bijections, SIAM J. Discrete Math. 24 (2010), 113–145.
  • [PY13a] Igor Pak and Jed Yang, The complexity of generalized domino tilings, Electron. J. Comb. 20 (2013), no. 4, Paper 12, 23 pp.
  • [PY13b] Igor Pak and Jed Yang, Tiling simply connected regions with rectangles, J. Combin. Theory, Ser. A 120 (2013), 1804–1816.
  • [Pan23] Greta Panova, Computational complexity in algebraic combinatorics, to appear in Current Developments in Mathematics 2023, International Press, Boston, MA, 30 pp.;  arXiv:2306.17511.
  • [Pap94] Christos H. Papadimitriou, Computational Complexity, Addison-Wesley, Reading, MA, 1994, 523 pp.
  • [PB83] J. Scott Provan and Michael O. Ball, The complexity of counting cuts and of computing the probability that a graph is connected, SIAM J. Comput. 12 (1983), 777–788.
  • [RT10] Kári Ragnarsson and Bridget E. Tenner, Obtainable sizes of topologies on finite sets, J. Combin. Theory, Ser. A 117 (2010), 138–151.
  • [Ras04] Etienne Rassart, A polynomiality property for Littlewood–Richardson coefficients, J. Combin. Theory, Ser. A 107 (2004), 161–179.
  • [Rém04] Eric Rémila, The lattice structure of the set of domino tilings of a polygon, Theoret. Comput. Sci. 322 (2004), 409–422.
  • [Sch78] Thomas J. Schaefer, The complexity of satisfiability problems, in Proc.  10th STOC, ACM, New York, 1978, 216–226.
  • [Sch19] Rolf Schiffler, Snake graphs, perfect matchings and continued fractions, in Snapshots of modern mathematics from Oberwolfach, 2019, No 1, 10 pp; available at  publications.mfo.de/handle/mfo/1405
  • [Sch90] Uwe Schöning, The power of counting, in Complexity theory retrospective, Springer, New York, 1990, 204–223.
  • [Sed70] Jiří Sedláček, On the minimal graph with a given number of spanning trees, Canad. Math. Bull. 13 (1970), 515–517.
  • [Shah22] Rikhav Shah, Determinants of binary matrices achieve every integral value up to Ω(2n/n)\Omega(2^{n}/n), Linear Algebra Appl. 645 (2022), 229–236.
  • [SH91] Gerard Sierksma and Han Hoogeveen, Seven criteria for integer sequences being graphic, J. Graph Theory 15 (1991), 223–231.
  • [OEIS] Neil J. A. Sloane, The online encyclopedia of integer sequences, oeis.org
  • [Sno12] Michael Snook, Counting bases of representable matroids, Electron. J. Combin. 19 (2012), no. 4, Paper 41, 11 pp.
  • [Sta10] Richard P. Stanley, A survey of alternating permutations, in Combinatorics and graphs, AMS, Providence, RI, 2010, 165--196.
  • [Sta12] Richard P. Stanley, Enumerative Combinatorics, vol. 1 (second ed.) and vol. 2, Cambridge Univ. Press, Cambridge, UK, 2012, 626 pp., and 1999, 581 pp.
  • [Ste14] John Stembridge, Generalized stability of Kronecker coefficients, preprint (2014), 27 pp.;  available at  tinyurl.com/429pkka3
  • [Sto22] Richard Stong, Minimal graphs with a prescribed number of spanning trees, Australas. J. Combin. 82 (2022), 182--196.
  • [Tar91] Jun Tarui, Randomized polynomials, threshold circuits, and the polynomial hierarchy, in Proc. 8th STACS (1991), 238--250.
  • [Ten09] Bridget E. Tenner, Optimizing linear extensions, SIAM J. Discrete Math. 23 (2009), 1450--1454.
  • [Thu90] William P. Thurston, Conway’s tiling groups, Amer. Math. Monthly 97 (1990), 757--773.
  • [Toda91] Seinosuke Toda, PP is as hard as the polynomial-time hierarchy, SIAM J. Comput. 20 (1991), 865--877.
  • [Tro95] William T. Trotter, Partially ordered sets, in Handbook of combinatorics, vol. 1, Elsevier, Amsterdam, 1995, 433--480.
  • [Vad01] Salil P. Vadhan, The complexity of counting in sparse, regular, and planar graphs, SIAM J. Comput. 31 (2001), 398--427.
  • [Val79a] Leslie G. Valiant, The complexity of enumeration and reliability problems, SIAM J. Comp. 8 (1979), 410--421.
  • [Val79b] Leslie G. Valiant, Completeness classes in algebra, in Proc. 11th STOC (1979), 249--261.
  • [Val79c] Leslie G. Valiant, The complexity of computing the permanent, Theoret. Comput. Sci. 8 (1979), 189--201.
  • [Vat15] Vincent Vatter, Permutation classes, in Handbook of enumerative combinatorics, CRC Press, Boca Raton, FL, 2015, 753--833.
  • [Wel93] Dominic J. A. Welsh, Complexity: knots, colourings and counting, Cambridge Univ. Press, Cambridge, UK, 1993, 163 pp.
  • [Wig19] Avi Wigderson, Mathematics and computation, Princeton Univ. Press, Princeton, NJ, 2019, 418 pp.; available at  math.ias.edu/avi/book
  • [Win86] Peter M. Winkler, Correlation and order, in Combinatorics and ordered sets, AMS, Providence, RI, 1986, 151--174.
  • [Wor18] Nicholas Wormald, Asymptotic enumeration of graphs with given degree sequence, in Proc. ICM Rio de Janeiro, Vol. 3, 2018, 3229--3248.