---
abstract: 'MacMahon’s theorem on plane partitions yields a simple product formula for tiling number of a hexagon, and Cohn, Larsen and Propp’s theorem provides an explicit enumeration for tilings of a dented semihexagon via semi-strict Gelfand–Tsetlin patterns. In this paper, we prove a natural hybrid of the two theorems for hexagons with an arbitrary set of unit triangles removed along the a horizontal axis. In particular, we show that the ‘shuffling’ of removed unit triangles only changes the tiling number of the region by a simple multiplicative factor. Our main result generalizes a number of known enumerations and asymptotic enumerations of tilings. We also reveal connections of the main result to the study of symmetric functions and $q$-series.'
address:
- 'Department of Mathematics, University of Nebraska – Lincoln, Lincoln, NE 68588, U.S.A.'
- 'Department of Mathematics and Computer Science, Saint Mary’s College, Notre Dame, IN 46556, U.S.A.'
author:
- Tri Lai
- Ranjan Rohatgi
title: 'A shuffling theorem for lozenge tilings of doubly-dented hexagons'
---

[^1]

Introduction
============

MacMahon’s classical theorem [@Mac] on plane partition fitting in a given box is equivalent to the fact that the number of lozenge tilings of a centrally symmetric hexagon $H(a,b,c)$ of side-lengths $a,b,c,a,b,c$ (in this cyclic order) is given by the simple product: $$\label{Maceq}
{\operatorname{PP}}(a,b,c):=\prod_{i=1}^{a}\prod_{j=1}^{b}\prod_{k=1}^{c}\frac{i+j+k-1}{i+j+k-2}.$$

This formula was generalized by Cohn, Larsen and Propp [@CLP Proposition 2.1] when they presented a correspondence between lozenge tilings of a semihexagon with unit triangles removed on the base and semi-strict Gelfand-Tsetlin patterns. In particular, the *dented semihexagon* $S_{a,b}(s_1,s_2,\dots,s_a)$ is the region obtained from the upper half of the symmetric hexagon of side-lengths $b,a,a,b,a,a$ (in clockwise order, starting from the north side) by removing $a$ up-pointing unit triangles along the base at the positions $s_1,s_2,\dotsc,s_a$ from left to right. The number of lozenge tilings of the dented semihexagon is given by $$\label{CLPeq}
{\operatorname{M}}(S_{a,b}(s_1,s_2,\dots,s_a))=\prod_{1\leq i<j \leq a}\frac{s_j-s_i}{j-i},$$ where we use the notation ${\operatorname{M}}(R)$ for the number of lozenge tilings of the region $R$.

In this paper, we consider a hybrid object between MacMahon’s hexagon and Cohn–Larsen–Propp’s dented semihexagon. Our region is a hexagon on the triangular lattice, as in the case of MacMahon’s theorem, with an arbitrary set of unit triangles removed along a horizontal axis, like the dents in Cohn–Larsen–Propp’s theorem (see Fig. \[multiplefernfig\] A). In general, the tiling numbers of such regions are not given by simple product formula. However, we show that their tiling number only changes by a simple multiplicative factor when we shuffle the positions of up- and down-pointing removed triangles (see Theorem \[factorization\] and its generalizations in Theorems \[genfactorization\] and \[qfactorization\]).

Our main theorems imply a number of known tiling enumerations of regions with ‘holes’ (e.g. [@Ciu1; @Twofern; @HoleDent]). Here, a *hole* is a portion removed from a region. We also show that our main theorems can be used to obtain new results in asymptotic enumeration of tilings, including the enumeration of the so-called ‘*doubly–dented hexagon*’ in [@Twofern], the main results of the first author about hexagon with three arrays of triangles removed in [@HoleDent Theorems 2.11 and 2.12] and Ciucu’s main results about ‘*$F$-cored hexagons*’ in [@Ciu1 Theorems 1.1 and 2.1] (see Corollary \[cordual\]).

Shuffling theorems
==================

Let $x,y,n,u,d$ be nonnegative integers, such that $u,d \leq n$. Consider a symmetric hexagon of side-lengths $x+n-u,y+u,y+d,x+n-d,y+d,y+u$ in clockwise order, starting from the north side. We remove $u+d$ arbitrary unit triangles along the lattice line $l$ that contains the west and the east vertices of the hexagon. Assume further that, among these $u+d$ removed triangles, there are $u$ up-pointing ones and $d$ down-pointing ones. Let $U=\{s_1,s_2,\dotsc,s_u\}$ and $D=\{t_1,t_2,\dots, t_d\}$ be, respectively, the sets of positions of the up-pointing and down-pointing removed unit triangles (ordered from left to right), such that $|U\cup D|=n$ (i.e., $U,D\subseteq [x+y+n]$, $U$ and $D$ are not necessarily disjoint). Denote by $H_{x,y}(U; D)$ the resulting region . See Fig. \[multiplefernfig\] A for an example of such a region and Fig. \[multiplefernfig\] B for a sample tiling; we ignore the two horizontal unit “barriers” at the positions 6 and 13 on $l$ at the moment.

We now consider ‘shuffling’ the up- and down-pointing unit triangles in the symmetric difference $U\Delta D$ to obtain new position sets $U'$ and $D'$ for the up-pointing and down-pointing removed triangles. (In particular, $U$ and $U'$ have the same size, and so do $D$ and $D'$.) The following theorem shows that the shuffling of removed triangles only changes the tiling number by a simple multiplicative factor. Moreover, the factor can be written in a similar form to Cohn–Larsen–Propp’s formula (i.e. the product on the right-hand side of Eq. \[CLPeq\]).

![(A) The region $H_{4,3}(2,4,5,8,11;\ 4,9,11,12)$ and (B) a lozenge tiling of it. The black and shaded triangles indicate the unit triangles removed.[]{data-label="multiplefernfig"}](Garciafest3.eps){width="12cm"}

\[factorization\] For nonnegative integers $x,y,n,u,d$ ($u,d\leq n$) and four ordered subsets $U=\{s_1,s_2,\dotsc,s_u\}$, $D=\{t_1,t_2,\dots, t_d\}$, $U'=\{s'_1,s'_2,\dotsc,s'_{u}\}$, and $D'=\{t'_1,t'_2,\dots, t'_{d}\}$ of $[x+y+n]$ such that $U\cup D =U'\cup D'$, and $U\cap D =U'\cap D'$. Then $$\label{maineq1}
  \frac{{\operatorname{M}}(H_{x,y}(U; D))}{{\operatorname{M}}(H_{x,y}(U'; D'))}=\prod_{1\leq i <j\leq u}\frac{s_j-s_i}{s'_j-s'_i}\displaystyle\prod_{1\leq i <j\leq d}\frac{t_j-t_i}{t'_j-t'_i}.$$

We would like to emphasize that, in general, the numbers of tilings of two regions on the left-hand side of Eq. \[maineq1\] are *not* given by simple product formulas.

By Cohn–Larson–Propp’s theorem (see Eq. \[CLPeq\]), Eq. \[maineq1\] in Theorem \[factorization\] can be written in terms of tiling numbers as $$\begin{aligned}
\label{geointereq}
  \frac{{\operatorname{M}}(H_{x,y}(U; D))}{{\operatorname{M}}(H_{x,y}(U'; D'))}=\frac{{\operatorname{M}}(S_{u,x+y+n-u}(U)){\operatorname{M}}(S_{d,x+y+n-d}(D))}{{\operatorname{M}}(S_{u,x+y+n-u}(U')){\operatorname{M}}(S_{d,x+y+n-d}(D'))}.\end{aligned}$$ It would be interesting to have a combinatorial explanation for this identity.

Moreover, one readily sees that the two dented semihexagons in the numerator of the right-hand side are obtained by dividing the region $H_{x+y,0}(U; D)$ along the horizontal axis $l$. Similarly, the two dented semihexagons in the denominator are obtained by dividing $H_{x+y,0}(U'; D')$ along the horizontal axis. This means that, identity (\[geointereq\]) is equivalent to $$\label{geointereq2}
{\operatorname{M}}(H_{x,y}(U; D)){\operatorname{M}}(H_{x+y,0}(U'; D'))={\operatorname{M}}(H_{x,y}(U'; D')){\operatorname{M}}(H_{x+y,0}(U; D)).$$ The both sides of (\[geointereq2\]) count pairs of tilings of doubly-dented hexagons. It would be interesting to find a bijective proof for this identity.

We can generalize our Shuffling Theorem \[factorization\] by additionally allowing the unit triangles in the symmetric difference $U \Delta D$ to ‘flip’ (from up-pointing to down-pointing, and vice versa). It is possible, then, that the new position sets of removed up-pointing triangles and removed down-pointing triangles $U'$ and $D'$ may have sizes *different* than those of $U$ and $D$. We also allow the appearance of “*barriers*” along the axis $l$. A barrier is a unit horizontal lattice interval which is not allowed to be contained within a lozenge in a tiling. Assume that we have a set of barriers at the positions $B\subseteq [x+y+n]-U\cup D$ so that vertical lozenges may not appear at the positions in $B$ and that $|B|\leq x$ (see the red barriers in Fig. \[multiplefernfig\]; $B=\{6,13\}$ in this case). We now consider the tilings of $H_{x,y}(U;\ D)$ which are compatible with the set of barriers $B$. Denote by $H_{x,y}(U;D;B)$ the doubly-dented hexagons with such setup of removed unit triangles and barrier.

\[genfactorization\] For nonnegative integers $x,y,n,u,d$ ($u,d\leq n$) and five ordered subsets $U=\{s_1,s_2,\dotsc,s_u\}$, $D=\{t_1,t_2,\dots, t_d\}$, $U'=\{s'_1,s'_2,\dotsc,s'_{u'}\}$, $D'=\{t'_1,t'_2,\dots, t'_{d'}\}$, and $B:=\{k_1,k_2,\dots,k_b\}$ of $[x+y+n]$ such that $U\cup D =U'\cup D'$, $U\cap D =U'\cap D'$, $B\cup (U\cup D)=\emptyset$, and $|B|\leq x$, we always have

An surprising fact is that the right-hand side of Eq. \[genmaineq\] does *not* depend on the barrier set $B$. It would be quite interesting to explain this phenomenon combinatorially.

We note that the above generalized shuffling theorem can be viewed as the enumeration of certain ‘restricted tilings’ in the case when $U\cap D=\emptyset$. Indeed, the tilings of $H_{x,y}(U;\ D;\ B)$ are in bijection with tilings of $H_{x,y+(u+d-n)}(U\setminus D; D\setminus U)$ which contain a fixed vertical lozenge at each position in $U \cap D$ and which do not contain a vertical lozenge at each position in $B$. These types of tiling enumerations were investigated by Fischer [@Fischer] and Fulmek and Krattenthaler in [@FK1; @FK2].

![Assigning weights to lozenges in a doubly-dented hexagon.[]{data-label="weight"}](PNASweighted.eps){width="7cm"}

We can assign to each right-tilting lozenge above $l$ (resp., below $l$) a weight $q^{z}$ (resp, a weight $q^{-z}$), where $(z-\frac{1}{2})\frac{\sqrt{3}}{2}$ is the distance from its center to the axis $l$ (see Figure \[weight\]). The ‘*tiling-generating function*’ of a region $H_{x,y}(U;D; B)$ is the sum of the weights of all tilings of the region, where the *weight* of a tiling is the product of weights of all its constituent lozenges. We denote this tiling-generating function by ${\operatorname{M}}_q(H_{x,y}(U;D;B))$. We have a $q$-analog of Theorem \[genfactorization\]:

\[qfactorization\] With the same assumptions as in Theorem \[genfactorization\], we have $$\begin{aligned}
\label{qmianeq}
  \frac{{\operatorname{M}}_q(H_{x,y}(U;D;B))}{{\operatorname{M}}_q(H_{x,y}(U';D';B))}=q^C \frac{\displaystyle\prod_{1\leq i <j\leq u}\frac{q^{s_j}-q^{s_i}}{q^{j}-q^{i}}\displaystyle\prod_{1\leq i <j\leq d}\frac{q^{t_j}-q^{t_i}}{j-i}{\operatorname{PP}}_q(u,d,y)}{\displaystyle\prod_{1\leq i <j\leq u'}\frac{q^{s'_j}-q^{s'_i}}{q^{j}-q^{i}}\displaystyle\prod_{1\leq i <j\leq d'}\frac{q^{t'_j}-q^{t'_i}}{q^{j}-q^{i}}{\operatorname{PP}}_q(u'd',y)}.
   \end{aligned}$$ where $${\operatorname{PP}}_q(a,b,c):=\prod_{i=1}^{a}\prod_{j=1}^{b}\prod_{k=1}^{c}\frac{1-q^{i+j+k-1}}{1-q^{i+j+k-2}}.$$ and where $$\begin{aligned}
C=(d-x-n)\binom{y+d+1}{2}-(d'-x-n)\binom{y+d'+1}{2}+udy-u'd'y.\end{aligned}$$

We note that MacMahon’s Theorem [@Mac] states that the number of plane partitions fitting in an $a\times b \times c$- box of given by ${\operatorname{PP}}_q(a,b,c)=\prod_{i=1}^{a}\prod_{j=1}^{b}\prod_{k=1}^{c}\frac{1-q^{i+j+k-1}}{1-q^{i+j+k-2}}$. The formula (\[Maceq\]) follows by letting $q \to 1$.

We note that our shuffling theorems (Theorems \[factorization\], \[genfactorization\], and \[qfactorization\]) above only for symmetric hexagons with removed unit triangles and barriers (we call them generally the ‘*obstacles*’) running along the horizontal axis $l$ that passes the east and west vertices of the hexagons. One would ask for similar results for general hexagons (not necessarily symmetric) with obstacles on an arbitrary horizontal lattice line (not necessarily passing a vertex of the hexagon). However, the latter case can be implied from our main theorems via ‘*forced lozenges*’. A *forced lozenge* in a region $R$ is a lozenge that appears in every tiling of $R$. In the unweighted case, the removal of forced lozenges does not change the tiling number of the region. We consider a general hexagons of side-lengths $x+n-u, y+u, z+d,x+n-d,y+d,z+u$ (in the clockwise order from the north side), and we would like removing unit triangles and placing barriers along an *arbitrary* horizontal lattice line $l$. There are four possible cases depending on the relative positions of the axis $l$ and the east and the west vertices of the hexagon. Figure \[forced2\] shows that, in all four cases, by removing forced lozenges from a symmetric doubly-dented hexagon as considered in our shuffling theorems, we can obtain general doubly-dented hexagon with an arbitrary position of the axis $l$. This means that one can obtain similar shuffling theorems for general hexagons with obstacles on an arbitrary axis from our shuffling theorems.

![Obtaining general doubly-dented hexagons from symmetric doubly-dented hexagons by removing forced lozenges.[]{data-label="forced2"}](PNASforced2.eps){width="12cm"}

An asymptotic enumeration
=========================

Ciucu and Krattenthaler in [@CK] proved a counterpart of MacMahon’s theorem (Eq. \[Maceq\]) by obtaining the asymptotic tiling number of the exterior of a concave hexagonal contour in which we turn $120^{\circ}$ after drawing each edge (see Fig. \[contour\] B; the tiling number in MacMahon’s theorem is for the interior of a hexagonal contour in which each turn $60^{\circ}$ after drawing each edge as in Fig. \[contour\] A). Ciucu [@Ciu1] later obtained a similar counterpart of Cohn–Larsen–Propp’s theorem corresponding to the exterior of a concave polygon with an arbitrary number sides (see the contour in Fig. \[contour\] C). Recently, the first author generalized the asymptotic result of Ciucu to the union of three polygons [@HoleDent]. In Corollary \[cordual\], we will show a multi-parameter generalization of the latter two asymptotic results.

![Three contours: (A) the contour in MacMahon’s theorem, (B) the contour in [@CK], and (C) the contour in [@Ciu1].[]{data-label="contour"}](PNASContour.eps){width="12cm"}

We now assume that the set of removed unit triangles is partitioned into $k$ separated clusters (i.e. chains of contiguous unit triangles). Denote these clusters by $C_1,C_2,\dotsc,C_k$ and the distances between them by $d_1,d_2,\dotsc,d_{k-1}$ ($d_i >0$), as they appear from left to rights. For the sake of convenience, we assume that $C_1$ is attached to the west vertex of the hexagon, that $C_k$ is attached to the east vertex of the hexagon, and that $C_1$ and $C_k$ may be empty. We use the notation $H_{x,y}(C_1,\dots,C_k; d_1,\dots,d_{k-1})$ for these regions (see Fig. \[fern\] for an example; the black unit triangles indicate the ones removed). For each cluster $C_i$, we use the notations $U_i$ and $D_i$ for the index sets of its up-pointing and down-pointing triangles. For each cluster $C_i$, we can shuffle unit triangles at the positions in $U_i \Delta D_i$ to obtain a new cluster $C'_i$, for $i=1,2,\dots,k$. The index sets of triangles in $C'_i$ are denoted by $U'_i$ and $D'_i$. Assume that $|U_i|=u_i$, $D_i=d_i$, $|U'_i|=u'_i$, $|D_i|=d'_i$ and $|U_i\cup D_i|=|U'_i\cup D'_i|=f_i$. We call $f_i$ the *length* of the cluster $C_i$ (and also the length of $C'_i$).

We now consider the behavior of the tiling number of the region when the side-lengths of the outer hexagon and the distances between two consecutive clusters get large.

\[cordual\]For nonnegative integers $x$ and $y$, $$\begin{aligned}
\label{dualeq}
\lim_{N\to \infty}&\frac{{\operatorname{M}}(H_{Nx,Ny}(C_1,\dots,C_k; Nd_1,\dots,Nd_{k-1}))}{{\operatorname{M}}(H_{Nx,Ny}(C'_1,\dots,C'_k; Nd_1,\dots,Nd_{k-1}))}\notag\\
&=\prod_{i=1}^{k}\frac{s^+(C_i)s^-(C_i)}{s^+(C'_i)s^-(C'_i)}\end{aligned}$$ where $s^+(C_i)={\operatorname{M}}(S_{u_i,f_i-u_i}(U_i))$ and $s^-(C_i)={\operatorname{M}}(S_{d_i,f_i-d_i}(D_i))$ are respectively the tiling numbers of the dented semihexagons whose dents are defined by the up-pointing triangles and down-pointing triangles in the cluster $C_i$, and where $s^+(C'_i)$ and $s^-(C'_i)$ are defined similarly with respect to $C'_i$.

Corollary \[cordual\] can be visualized as in Fig. \[geointer\], for $k=3$. The dented semihexagons corresponding to $s^+(C_i)$ and $s^-(C_i)$ are the upper and lower halves of the ‘numerator hexagon’ in the $i$th fraction on the right-hand side; the dented semihexagons corresponding to $s^+(C'_i)$ and $s^-(C'_i)$ are the upper and lower halves of the ‘denominator hexagon’ in the $i$th fraction, for $i=1,2,3$.

![Illustrating Corollary \[cordual\].[]{data-label="geointer"}](PNASgeointer.eps){width="14cm"}

Assume that $U_i=\{a^{(i)}_1,\dots,a^{(i)}_{u_i}\}$ and $U'_i=\{e^{(i)}_1,\dots,e^{(i)}_{u_i}\}$. Applying Theorem \[factorization\] to the regions $H_{Nx,Ny}(U; D)$ and $H_{Nx,Ny}(U'; D')$, for $U:=\bigcup_iU_i$, $D:=\bigcup_jD_j$, $U':=\bigcup_iU'_i$, and $D':=\bigcup_jD'_j$, we get $$\frac{{\operatorname{M}}(H_{Nx,Ny}(C_1,\dots,C_k; Nd_1,\dots,Nd_{k-1}))}{{\operatorname{M}}(H_{Nx,Ny}(C'_1,\dots,C'_k; Nd_1,\dots,Nd_{k-1}))}=\frac{\Delta(U)\Delta(D)}{\Delta(U')\Delta(D')},$$ where the operation $\Delta$ is defined as $\Delta(S):=\prod_{1\leq i<j \leq n}(s_j-s_i)$ for an ordered set $S=\{s_1<s_2<\cdots<s_n\}$. However, we find it more convenient to observe that $\Delta(S)^2=\prod_{1\leq i\not=j \leq n}|s_j-s_i|$ here.

We have $\frac{\Delta(U)^2}{\Delta(U')^2}=\prod_{i,j}\prod_{p,q}\frac{|a^{(i)}_p-a^{(j)}_q|}{|e^{(i)}_p-e^{(j)}_q|}$, where $p\not=q$ if $i=j$. It is easy to see that if $i\not=j$ the fraction $\frac{|a^{(i)}_p-a^{(j)}_q|}{|e^{(i)}_p-e^{(j)}_q|}$ tends to $1$, as $N$ gets large (for any $p,q$). Thus $\frac{\Delta(U)^2}{\Delta(U')^2}$ tends to $\prod_{i=1}^{k}\prod_{p\not=q}\frac{|a^{(i)}_p-a^{(i)}_q|}{|e^{(i)}_p-e^{(i)}_q|}=\prod_{i=1}^{k}\frac{\Delta(U_i)^2}{\Delta(U'_i)^2}$.

Moreover, $\frac{\Delta(U_i)}{\Delta(U'_i)} = \frac{s^{+}(C_i)}{s^{+}(C'_i)}$. Thus, $\frac{\Delta(U)}{\Delta(U')}$ tends to $\prod_{i=1}^{k}\frac{s^+(C_i)}{s^+(C'_i)}$. Similarly, we see that $\frac{\Delta(D)}{\Delta(D')}$ tends to $\prod_{i=1}^{k}\frac{s^-(C_i)}{s^-(C'_i)}$, completing the proof.

By the same arguments, one would imply a $q$-analog of the above asymptotic result from $q$-Shuffling Theorem \[qfactorization\].

We now consider the special case in which $U\cap D=\emptyset$ (i.e. in each cluster $C_i$ we have $U_i \cap D_i=\emptyset$). Each cluster $C_i$ can be partitioned into maximal intervals of triangles of the same orientation (we call each of these intervals an ‘*up-interval*’ or a ‘*down-interval*’ if it consists of up-pointing triangles or down-pointing triangles, respectively). For each cluster $C_i$, we can remove forced vertical lozenges above each up-interval and below each down-interval composed of two or more unit triangles. We obtain a new region with the same tiling number in which each cluster is replaced by a chain of removed equilateral triangles of alternating orientations (see Fig. \[fern\]; the forced lozenges are colored white). Each such chain of triangles is called a ‘*fern*’ (see e.g. [@Ciu1; @HoleDent]); the side-lengths of triangles in a fern are equal to the lengths of the intervals of unit triangles of the same orientation in the corresponding cluster. Denote by $E_{x,y}(F_1,\dotsc,F_k;d_1,\dots,d_{k-1})$ the corresponding hexagon with ferns removed (the fern $F_i$ corresponds to the cluster $C_i$; and the fern $F'_i$ corresponds to the cluster $C'_i$). By setting $k=3$, $d_1=d_2$ or $d_1=d_2-1$, and specifying that the cluster $C'_i$ has all its up-pointing triangles on the left and all its down-pointing triangles on the right, for $i=1,2,3$ (equivalently, the fern $F'_i$ consists of two triangles, an up-pointing one followed by a down-pointing one), our Corollary \[cordual\] implies the first author’s work in [@HoleDent Theorem 2.11]. Similarly, we can recover the work of Ciucu in [@Ciu1 Theorem 1.1] by further specifying that $C_1=C_3=\emptyset$ (i.e. we actually have only a non-empty fern $F_2$ in the center of the region).

![Obtained a hexagon with ferns removed from the region $H_{x,y}(C_1,\dots,C_k; d_1,\dots,d_{k-1})$.[]{data-label="fern"}](PNASfern.eps){width="7cm"}

Proof of the main theorem
=========================

In general, the lozenges in a region can carry weights. In the weighted case, ${\operatorname{M}}(R)$ denotes the sum of weights of the tilings in $R$, where the *weight* of a tiling is the product of weights of its constituent lozenges. The removal of one or more forced lozenges changes the number of tilings of the region by a factor equal to the reciprocal of the weighted product of the forced lozenges. In particular, if we remove the lozenges $l_1,l_2,\dots,l_k$ from a region $R$ and get a new region $R'$, then $${\operatorname{M}}(R')=\left(\prod_{i=1}^{k}wt(l_i)\right)^{-1}{\operatorname{M}}(R),$$ where $wt(l_i)$ is the weight of the forced lozenge $l_i$.

A *(perfect) matching* of a graph is a collection of disjoint edges that covers all the vertices of the graph. When edges of a graph carry weights, we denote by ${\operatorname{M}}(G)$ the weighted sum of matchings in $G$, where the weight of a matching is the product of weights of its edges. (In the unweighted case, ${\operatorname{M}}(G)$ is exactly the number of matchings of $G$.) The *(planar) dual graph* of the region $R$ on the triangular lattice is the graph whose vertices are the unit triangles in $R$ and whose edges connect precisely two unit triangles sharing an edge. The edges of the dual graph inherit the weights from the lozenges of the corresponding region. There is a natural (weight-preserving) bijection between tilings of a region and matchings of its dual graph.

Our proof is based on the following powerful graphical condensation lemma first introduced by Kuo [@Kuo]:

\[kuothm\] Let $G=(V_1,V_2,E)$ be a (weighted) planar bipartite graph with the two vertex classes $V_1$ and $V_2$ such that $|V_1|=|V_2|+1$. Assume that $u,v,w,s$ are four vertices appearing in this cyclic order around a face of $G$, such that $u,v,w\in V_1$ and $s\in V_2$. Then $$\begin{aligned}
{\operatorname{M}}(G-v)&{\operatorname{M}}(G-\{u,w,s\})={\operatorname{M}}(G-u){\operatorname{M}}(G-\{v,w,s\})+{\operatorname{M}}(G-w){\operatorname{M}}(G-\{u,v,s\}).\end{aligned}$$

If a region admits a lozenge tiling, then it must have the same number of up-pointing and down-pointing unit triangles. We call such a region *balanced*. The following lemma allows us to decompose a region into smaller regions when enumerating tilings in certain situations.

\[RS\] Let $R$ be a balanced region on the triangular lattice. Assume that a balanced sub-region $Q$ of $R$ satisfies the condition that the unit triangles in $Q$ that are adjacent to some unit triangle of $R-Q$ have the same orientation. Then ${\operatorname{M}}(R)={\operatorname{M}}(Q)\, {\operatorname{M}}(R-Q).$

Given $k$ positive integers $\lambda_1\geq \lambda_2\geq \dots \geq \lambda_k$, a *plane partition* of shape $(\lambda_1,\lambda_2,\dots,\lambda_k)$ is an array of non-negative integers

[rccccccccc]{} $n_{1,1}$ &$n_{1,2}$ &$n_{1,3}$ & $\dotsc$ & $\dotsc$ & $\dotsc$ & $n_{1,\lambda_1}$\
$n_{2,1}$ & $n_{2,2}$ & $n_{2,3}$ & $\dotsc$ & $\dotsc$ & $n_{2,\lambda_2}$&\
$\vdots$ & $\vdots$ & $\vdots$ & $\vdots$ & & &\
$n_{k,1}$ & $n_{k,2}$ & $n_{k,3}$ & $\dotsc$ & $n_{k,\lambda_k}$ & &\

so that $n_{i,j}\geq n_{i,j+1}$ and $n_{i,j}\geq n_{i+1,j}$ (i.e. all rows and all columns are weakly decreasing from left to right and from top to bottom, respectively). The sum of all entries of a plane partition $\pi$ is called the *volume* (or the *norm*) of the plane partition, and denoted by $|\pi|$.

A *column-strict plane partition* is a plane partition having columns strictly decreasing. A *column-strict plane partition* is a plane partition having columns strictly decreasing. We now assign to each right-tilting lozenge in the dented semihexagon similarly to the upper half of the doubly dented hexagon. In particular, a eight-tilting lozenge is weighted by $q^{z}$, where the distance between its center and the base of the region is $(z-\frac{1}{2})\frac{\sqrt{3}}{2}$. Denote ${\operatorname{M}}_q(S_{a,b}(s_1,s_2,\dots,s_a))$ the corresponding tiling generating function of the dented semihexagon $S_{a,b}(s_1,s_2,\dots,s_a)$.There is a well-known (weight preserving) bijection between the lozenge tilings of $S_{a,b}(s_1,s_2,\dots,s_a)$ and the column-strict plane partitions of shape $(s_{a}-a,s_{a-1}-a+1,\dotsc,s_1-1)$ with positive entries at most $a$ (see e.g. [@CLP] and [@Car]). The following weighted enumeration of the dented semihexagon follows directly from the bijection and equation (7.105) in [@Stanley pp. 375].

\[qCLP\] For nonnegative $a,b$ and positive integers $1\leq s_1<s_2<\cdots<s_a$ $${\operatorname{M}}_q(S_{a,b}(s_1,s_2,\dots,s_a))=\sum_{\pi}q^{|\pi|}=q^{\sum_{i=1}^a(s_i-i)}\prod_{1\leq i <j \leq q}\frac{q^{s_j}-q^{s_i}}{q^{j}-q^{i}},$$ where the sum after the first equality sign is taken over all column-strict plane partitions $\pi$ of shape $(s_{a}-a,s_{a-1}-a+1,\dotsc,s_1-1)$ with positive entries at most $a$.

Denote by $g_{x,y}(U;D;U';D')$ the right-hand side of (\[qmianeq\])[^2], i.e. $$g_{x,y}(U;D;U';D')=q^C \frac{\displaystyle\prod_{1\leq i <j\leq u}\frac{q^{s_j}-q^{s_i}}{q^{j}-q^{i}}\displaystyle\prod_{1\leq i <j\leq d}\frac{q^{t_j}-q^{t_i}}{j-i}{\operatorname{PP}}_q(u,d,y)}{\displaystyle\prod_{1\leq i <j\leq u'}\frac{q^{s'_j}-q^{s'_i}}{q^{j}-q^{i}}\displaystyle\prod_{1\leq i <j\leq d'}\frac{q^{t'_j}-q^{t'_i}}{q^{j}-q^{i}}{\operatorname{PP}}_q(u'd',y)},$$ we need to show that $$\label{refineeq}
{\operatorname{M}}_q(H_{x,y}(U;D;B))=g_{x,y}(U;D;U';D'){\operatorname{M}}_q(H_{x,y}(U';D';B)).$$

We prove (\[refineeq\]) by induction on $x+y$. The base cases are the situations when $x=b$ and when $y=0$.

If $y=0$, then we apply Region-splitting Lemma \[RS\] to the region $R=H_{x,0}(U;D;B)$ with the subregion $Q$ the portion above the horizontal axis $l$. The subregion $Q$ is congruent with the dented semihexagons $S_{u,x+n-u}(U)$ weighted as in Lemma \[qCLP\], and the $R-Q$, after rotated $180^{\circ}$, is congruent with $S_{d,x+n-d}(r(D))$ weighted similarly with $q$ replaced by $q^{-1}$ (see Fig. \[basecase\] A). Here we use the notation $r(S)$ for the reflection of the index set $S$, i.e. the set $(x+y+n+1)-S=\{(x+y+n+1)-s: s \in S\}.$ We have $$\begin{aligned}
{\operatorname{M}}_q(H_{x,0}(U;D;B))={\operatorname{M}}_q(S_{u,x+n-u}(U)){\operatorname{M}}_{q^{-1}}(S_{d,x+n-d}(r(D))).\end{aligned}$$ Similarly, the tiling number of $R'=H_{x,0}(U';D';B)$ is also written by a product of tiling numbers of two dented semihexagons as $$\begin{aligned}
{\operatorname{M}}_q(H_{x,0}(U';D';B))={\operatorname{M}}_q(S_{u,x+n-u}(U')){\operatorname{M}}_{q^{-1}}(S_{d,x+n-d}(r(D'))),\end{aligned}$$ and (\[refineeq\]) follows from Lemma \[qCLP\].

If $x=b$, consider the subregion $Q$ of $R=H_{b,y}(U;D;B)$ that is obtained from the portion above the axis $l$ by removing all up-pointing unit triangles in $(U\cup D\cup B)^c$ (we note that in this case $|(U\cup D\cup B)^c|=y$). $Q$ is the dented semihexagon $S_{y+u,b+n-u}((U\cup D\cup B)^c\cup U)$ weighted as in Lemma \[qCLP\], and its complement, after removing forced lozenges at the positions in $(U\cup D\cup B)^c$ and rotating $180^{\circ}$, is the dented semihexagon $S_{y+d,b+n-u}((U\cup D\cup B)^c\cup D)$ weighted similarly with $q$ replaced by $q^{-1}$ (see Fig. \[basecase\] B). This way, the tiling number of $R$ is written as the product of tiling numbers of two dented semihexagons. Similarly, the tiling number of $R'=H_{b,y}(U';D';B)$ is also written by a product of tiling numbers of two dented semihexagons. Then (\[refineeq\]) also follows from Lemma \[qCLP\].

For the induction step, we assume that $x>b$, $y>0$, and that (\[refineeq\]) holds for any $H$-type regions whose sum of $x$- and $y$-parameters is strictly less than $x+y$. We will use Kuo condensation in Theorem \[kuothm\] to obtain a recurrence for the tiling generating function on the left-hand side of (\[refineeq\]), and we show that the expression on the right-hand side satisfies the same equation. Then (\[refineeq\]) follows from the induction principle.

We apply Kuo condensation to the *dual graph* $G$ of the region $R$ obtained from $H_{x,y}(U;D;B)$ by adding a layer of unit triangle on the top of the hexagon, with the four vertices $u,v,w,s$ as in Fig. \[hexkuo1\] (the region restricted by the bold contour indicates the region $H_{x,y}(U;D;B)$). In particular, the vertices $w$ and $v$ correspond to the up-pointing triangles at the first and the last positions in $(U\cup D\cup B)^c$, the vertex $v$ corresponds to the up-pointing triangle on the northeast corner of the region, and the vertex $s$ corresponds to the down-pointing triangle on the southeast corner of the region.

![Applications of the Region-splitting Lemma \[RS\] in the base cases: (A) $y=0$, (B) $x=b$.[]{data-label="basecase"}](PNASbase2.eps){width="12cm"}

![Applying Kuo condensation to a hexagon.[]{data-label="hexkuo1"}](Pnaskuo2b.eps){width="8cm"}

![Obtaining the recurrence for the tiling numbers.[]{data-label="hexkuo2"}](PNASkuo1b.eps){width="12cm"}

We consider the region corresponding to $G-v$ (i.e., the region obtained from $R$ by removing the $v$-triangle as in Fig. \[hexkuo2\] A). The removal of the $v$-triangle yields several forced lozenges along the top of the region. After removing these forced lozenges (this lozenge removal changes the tiling number of the region by a factor $q^{(y+u+1)(x+n-u-1)}$), we get back the region $H_{x,y}(U;D;B)$. This means that we get $$\label{kuoeqa}
{\operatorname{M}}(G-v)=q^{(y+u+1)(x+n-u-1)}{\operatorname{M}}_q(H_{x,y}(U; D;B)).$$ By considering forced lozenges in the regions corresponding to the graphs $G-\{u,w,s\},$ $G-u,$ $G-\{v,w,s\},$ $G-w,$ and $G-\{u,v,s\}$ (as shown in Fig. \[hexkuo2\] B – F, respectively), we get $$\label{kuoeqb}
{\operatorname{M}}(G-\{u,w,s\})={\operatorname{M}}_q(H_{x-1,y-1}(\alpha \beta U; D;B)),$$ $$\label{kuoeqc}
{\operatorname{M}}(G-u)={\operatorname{M}}_q(H_{x-1,y}(\beta U;D;B)),$$ $$\label{kuoeqd}
{\operatorname{M}}(G-\{v,w,s\})=q^{(y+u+1)(x+n-u-1)}{\operatorname{M}}_q(H_{x,y-1}(\alpha U;D;B)),$$ $$\label{kuoeqe}
{\operatorname{M}}(G-w)={\operatorname{M}}_q(H_{x-1,y}(\alpha U;D;B)),$$ $$\label{kuoeqf}
{\operatorname{M}}(G-\{u,v,s\})=q^{(y+u+1)(x+n-u-1)}{\operatorname{M}}_q(H_{x,y-1}(\beta U;D;B)),$$ where we use the notations $\alpha U$, $\beta U$ and $\alpha \beta U$ for the unions $U\cup\{\alpha \}$, $U\cup\{\beta \}$ and $U\cup\{\alpha,\beta\}$, respectively. Plugging Eqs. \[kuoeqa\]–\[kuoeqf\] into the equation in Kuo’s Lemma \[kuothm\], we get the recurrence (all powers of $q$ cancel out): $$\begin{aligned}
\label{recurrence1}
{\operatorname{M}}_q(H_{x,y}(U;D;B)){\operatorname{M}}_q(H_{x-1,y-1}(\alpha \beta U; D;B))=&{\operatorname{M}}_q(H_{x-1,y}(\beta U;D;B)){\operatorname{M}}_q(H_{x,y-1}(\alpha U;D;B))\notag\\
&+{\operatorname{M}}_q(H_{x-1,y}(\alpha U;D;B)){\operatorname{M}}_q(H_{x,y-1}(\beta U;D;B)).\end{aligned}$$

Denote by $h_{x,y}(U;D;U';D'):=g_{x,y}(U;D;U';D'){\operatorname{M}}(H_{x,y}(U';D';B))$ the expression on the right-hand side of (\[refineeq\]). To finish the proof, we shall show that the $h_{x,y}(U;D;U';D')$ also satisfies recurrence (\[recurrence1\]). Equivalently, we need to verify $$\begin{aligned}
\label{hrecur}
\frac{h_{x-1,y}(\beta U;D; \beta U';D')h_{x,y-1}(\alpha U;D; \alpha U';D'))}{h_{x,y}(U;D;U';D')h_{x-1,y-1}(\alpha\beta U; D; \alpha\beta U'; D')}+\frac{h_{x-1,y}(\alpha U;D;\alpha U';D')h_{x,y-1}(\beta U;D;\beta U';D'))}{h_{x,y}(U;D;U';D')h_{x-1,y-1}(\alpha\beta U; D;\alpha\beta U'; D')}=1.\end{aligned}$$ We claim

$$\label{ratioeq1}
\frac{g_{x-1,y}(\beta U;D; \beta U';D')g_{x,y-1}(\alpha U;D; \alpha U';D'))}{g_{x,y}(U;D;U';D')g_{x-1,y-1}(\alpha\beta U; D; \alpha\beta U'; D')}=1$$

and $$\label{ratioeq2}
\frac{g_{x-1,y}(\alpha U;D;\alpha U';D')g_{x,y-1}(\beta U;D;\beta U';D'))}{g_{x,y}(U;D;U';D')g_{x-1,y-1}(\alpha\beta U; D;\alpha\beta U'; D')}=1.$$

We have from the definition $$\begin{aligned}
\label{PPeq}
\frac{{\operatorname{PP}}_q(u,d,y)}{{\operatorname{PP}}_q(u',d',y)}=q^{-udy+u'd'y}\frac{\Delta_q([u])\Delta_q([d])}{\Delta_q([u'])\Delta_q([d'])}\frac{\Delta_q([y+u'])\Delta_q([y+d'])}{\Delta_q([y+u])\Delta_q([y+d])},\end{aligned}$$ where, for any ordered index set $S=\{s_1,s_2,\dots,s_k\}$, we define $\Delta_q(S):=\prod_{1\leq i < j\leq k} q^{s_j}-q^{s_i}$. Therefore, we can rewrite the $g$-function as: $$\begin{aligned}
g_{x,y}(U;D;U';D')=q^{(d-x-n)\binom{y+d+1}{2}-(d'-x-n)\binom{y+d'+1}{2}}\frac{\Delta_q(U)\Delta_q(D)}{\Delta_q(U')\Delta_q(D')}\frac{\Delta_q([y+u'])\Delta_q([y+d'])}{\Delta_q([y+u])\Delta_q([y+d])}.\end{aligned}$$ It is easy to see that terms correspond to the faction $\frac{\Delta_q([y+u'])\Delta_q([y+d'])}{\Delta_q([y+u])\Delta_q([y+d])}$ cancel out in (\[ratioeq1\]). The exponents of $q$ also cancel out easily.

Moreover, since the position sets of down-pointing triangles do not change, the corresponding $\Delta_q(D)$ and $\Delta_q(D')$ terms canceled out. This makes Eq. (\[ratioeq1\]) become $$\frac{\Delta_q(\beta U)}{\Delta_q(\beta U')}\frac{\Delta_q(\alpha U)}{\Delta_q(\alpha U')} = \frac{\Delta_q(U)}{\Delta_q(U')}\frac{\Delta_q(\alpha\beta U)}{\Delta_q(\alpha\beta U')}.$$ Dividing two sides of the above equation by $\frac{\Delta_q(U)^2}{\Delta_q(U')^2}$, we get both equal to: $$\begin{aligned}
\frac{\prod_{i=1}^{u}|\beta-s_i|_q}{\prod_{i=1}^{u'}|\beta-s'_i|_q}\frac{\prod_{i=1}^{u}|\alpha-s_i|_q}{\prod_{i=1}^{u'}|\alpha-s'_i|_q},\end{aligned}$$ where the *‘$q$-absolute value’* $|x-y|_q$ is define to be $q^{x}-q^{y}$ if $x\geq y$, and is $q^{y}-q^{x}$ if $y>x$. This implies (\[ratioeq1\]). The identity (\[ratioeq2\]) follows from (\[ratioeq1\]) by interchanging the roles of $\alpha$ and $\beta$.

By (\[ratioeq1\]) and (\[ratioeq2\]) in the above claim, we have (\[hrecur\]) simplified as $$\begin{aligned}
\frac{{\operatorname{M}}_q(H_{x-1,y}(\beta U';D';B)){\operatorname{M}}_q(H_{x,y-1}(\alpha U';D';B))}{{\operatorname{M}}_q(H_{x,y}(U';D';B)){\operatorname{M}}_q(H_{x-1,y-1}(\alpha \beta U'; D';B))}+\frac{{\operatorname{M}}_q(H_{x-1,y}(\alpha U';D;B')){\operatorname{M}}_q(H_{x,y-1}(\beta U';D';B))}{{\operatorname{M}}_q(H_{x,y}(U';D';B)){\operatorname{M}}_q(H_{x-1,y-1}(\alpha \beta U'; D';B))}=1,\end{aligned}$$ which is obtained from the application of recurrence (\[recurrence1\]) to the region $H_{x,y}(U';D';B)$. We just verified that the function $h_{x,y}(U;D;U';D')$ (i.e. the expression on the right-hand side of (\[refineeq\])) satisfies (\[recurrence1\]). This finishes the proof.

Generalizations, more applications, and future directions
=========================================================

(1). Recall that when $U\cap D=\emptyset$ and $B=\emptyset$, our region (written in ‘cluster form’) $H_{x,y}(U; D)=H_{x,y}(U; D; \emptyset)=H_{x,y}(C_1,\dotsc,C_k; d_1,\dots,d_{k-1})$ is tiling-equinumerous with a hexagon with ferns removed $E_{x,y}(F_1,\dotsc,F_k;d_1,\dots,d_{k-1})$. Viewed in this way, our Theorem \[genfactorization\] implies a number of known enumerations about regions with ferns removed. The general idea is that we choose suitable position sets $U,D,U',D'$ so that the region $H_{x,y}(U; D)$ in Eq. (\[genmaineq\]) is the one that we want enumerate, and the region $H_{x,y}(U';D')$ is a known region up to removal of forced lozenges. Let us first consider the implication to [@HoleDent Theorem 2.12] as follows. Suppose we specialize our region as follows:

1.  $k=3$,

2.  $C_1$ and $C_3$ contain the same number of unit triangles,

3.  $|d_1-d_2|\leq 1$ (i.e. $C_2$ is centered between $C_1$ and $C_3$, or as close to centered as possible)

4.  all unit triangles in $C'_1$ have the same orientation,

5.  all unit triangles in $C'_3$ have opposite orientation to those in $C'_1$, and

6.  all unit triangles in $C'_2$ are up-pointing.

In this case, the region $H_{x,y}(U; D)$ becomes a hexagon with three ferns removed and $H_{x,y}(U';D')$, after the removal of forced lozenges, becomes a *cored hexagon* as in [@CEKZ] (see the illustration in Figure \[forced\]). Therefore, our Theorem \[genfactorization\] implies the enumeration of tilings of hexagons with three ferns removed in [@HoleDent Theorem 2.12].

By specifying that $C_1=C'_1$ and $C_3=C'_3$, in addition to the six conditions above, we recover Ciucu’s Theorem 2.1 in [@Ciu1] about $F$-cored hexagons (after removing forced lozenges, the region $H_{x,y}(U;D)$ becomes an $F$-cored hexagon and $H_{x,y}(U';D')$, as in the previous case, becomes a cored hexagon).

Similarly, if we require

1.  $k=2$,

2.  all unit triangles in $C'_1$ have the same orientation, and

3.  all unit triangles in $C'_2$ have opposite orientation to those in $C'_1$,

then we get Theorem 1.1 in [@Twofern] about hexagons with two ferns removed, called ‘*doubly–intruded hexagons*’ (the region $H_{x,y}(U;D)$ is now a doubly–intruded hexagon and $H_{x,y}(U';D')$ becomes a hexagon in MacMahon’s theorem, after removing forced lozenges).

![Obtaining cored hexagons from hexagons with three ferns removed.[]{data-label="forced"}](PNASforced.eps){width="12cm"}

(2). By equation (7.105) of [@StanleyB], we have[^3] $$\prod_{1\leq i< j \leq u}\frac{s_j-s_i}{j-i}=s_{\lambda(\{s_1,\dots,s_u\})}(\textbf{1}^{u}),$$ where the partition $\lambda(\{s_1,\dots,s_u\}):=(s_u-u+1,\dotsc,s_2-1,s_1)$. On the other hand, it is not hard to see that we also have $$\begin{aligned}
{\operatorname{M}}\left(H_{x,y}(U;D)\right)
=\sum_{|S|=y}s_{\lambda(U\cup S)}(\textbf{1}^{u+y})s_{\lambda(D \cup S)}(\textbf{1}^{d+y}),\end{aligned}$$ where the sum runs over all $y$-subsets $S$ of $ [x+y+n]-(U\cup D)$. Indeed, this follows from the fact that in each tiling, precisely $y$ of the $x+y$ unit segments on the lattice line from along which we removed the $n$ unit triangles are straddled by vertical lozenges. Hence one can write Eq. \[maineq1\] in Theorem \[factorization\] as $$\label{schureq}
\frac{  \sum_{|S|=y}\textbf{s}_{\lambda(U\cup S)}(\textbf{1}^{u+y})\textbf{s}_{\lambda(D \cup S)}(\textbf{1}^{d+y})}{  \sum_{|S|=y}\textbf{s}_{\lambda(U'\cup S)}(\textbf{1}^{u+y})\textbf{s}_{\lambda(D' \cup S)}(\textbf{1}^{d+y})}= \frac{\textbf{s}_{\lambda(U)}(\textbf{1}^{u})\textbf{s}_{\lambda(D)}(\textbf{1}^{d})}{\textbf{s}_{\lambda(U')}(\textbf{1}^{u})\textbf{s}_{\lambda(D')}(\textbf{1}^{d})},$$ where the sum is taken over all $y$-subsets $S$ of $ [x+y+n]-(U\cup D)$.

We note that $q$-Shuffling Theorem \[qfactorization\] gives further supporting evidence for the existence of a Schur function identity behind Shuffling Theorem \[factorization\], as we have $$\prod_{1\leq i <j\leq u}\frac{1-q^{s_j-s_i}}{1-q^{j-i}}=q^{A}\, \textbf{s}_{\lambda(\{s_1,\dotsc,s_u\})}(q,q^2,q^3,\dots),$$ for some constant $A$. This would imply that Eq. \[schureq\] is still true (up to a $q$-power) when $\textbf{1}^{n}$ is replaced by the sequence $(q,q^2,q^3,\dots,q^{n})$.

It would be interesting to know if the following general sum has a similar simplification: $$\label{schureq2}
\frac{  \sum_{|S|=y}\textbf{s}_{\lambda(U\cup S)}(\textbf{X}^{m+y})\textbf{s}_{\lambda(D \cup S)}(\textbf{X}^{m+(d-u)+y})}{  \sum_{|S|=y}\textbf{s}_{\lambda(U'\cup S)}(\textbf{X}^{m+y})\textbf{s}_{\lambda(D' \cup S)}(\textbf{X}^{m+(d-u)+y})},$$ where $\textbf{X}^{m}$ denotes the sequence of variables $x_1,x_2,\dotsc,x_m$.

(3). Motivated by Stanley’s classical paper [@Stanley] on symmetric plane partitions, we would like to investigate symmetric tilings of $H_{x,y}(U;D)$. There are two natural classes of symmetric tilings: the tilings which are invariant under a reflection over a vertical axis, and those which are invariant under a $180^{\circ}$ rotation (these tilings correspond to the transposed-complementary and self-complementary plane partitions). The shuffling theorems for these symmetry classes will be investigated in separate papers [@shuffling2; @shuffling3].

Acknowledgement {#acknowledgement .unnumbered}
===============

The authors would like to thank Brendon Rhoades for fruitful discussions. The first author (T.L.) was supported in part by Simons Foundation Collaboration Grant \#-585923.

[10]{}

L. Carlitz and R. Stanley, *Branching and partitions*, Proc. Amer. Math. Soc. **53**(1) (1975), 246–249.

M. Ciucu, Another dual of MacMahon’s theorem on plane partitions. 306 (2017):427–450.

M. Ciucu, T. Eisenkölbl, C. Krattenthaler and D. Zare, Enumeration of lozenge tilings of hexagons with a central triangular hole. 95 (2001): 251–334.

M. Ciucu and C. Krattenthaler, A dual of MacMahon’s theorem on plane partitions. 110 (2013):4518–4523.

M. Ciucu, T. Lai, Lozenge tilings of doubly–intruded hexagons. , 167 (2019): 294–339.

H. Cohn, M. Larsen, J. Propp, The shape of a typical boxed plane partition. 4 (1998):137–165.

I. Fischer, Enumeration of rhombus tilings of a hexagon which contain a fixed rhombus in the center. arXiv:math/9906102: 1–43.

M. Fulmek and C. Krattenthaler, The number of rhombus tilings of a symmetric hexagon which contain a fixed rhombus on the symmetry axis, I. 2(1998): 19–40.

M. Fulmek, C. Krattenthaler, The number of rhombus tilings of a symmetric hexagon which contain a fixed rhombus on the symmetry axis, II. 21 (2000): 601–640.

E.H. Kuo, Applications of graphical condensation for enumerating matchings and tilings. 319 (2004): 29–57.

T. Lai, Lozenge tilings of hexagons with central holes and dents. 1803.02792: 1–100.

T. Lai, A $q$-enumeration of lozenge tilings of a hexagon with four adjacent triangles removed from the boundary. 64 (2017): 66–87.

T. Lai, A $q$-enumeration of lozenge tilings of a hexagon with three dents. 82 (2017): 23–57.

T. Lai, A Shuffling Theorem for Reflectively Symmetric Lozenge Tilings, *Preprint: `arXiv:1905.09268`*.

T. Lai, A Shuffling Theorem for Centrally Symmetric Lozenge Tilings, *Preprint: `arXiv:1906.03759`*.

P.A. MacMahon, [*Combinatory Analysis*]{}. (Cambridge Univ. Press, Cambridge, UK) Vol. 2 (1916).

R. Stanley, [*Enumerative Combinatorics*]{}. (Cambridge Univ. Press) Vol. 2 (1999).

R. Stanley, Symmetries of plane partitions. 43 (1986): 103–113.

[^1]: This research was was supported in part by Simons Foundation Collaboration Grant (\# 585923).

[^2]: As the right-hand side of (\[qmianeq\]) does not depend on the index $B$, our $g$-function does not depend on $B$.

[^3]: The notation $\textbf{1}^n$ in the argument of a Schur function stands for $n$ arguments equal to 1.
