% EJC papers *must* begin with the following two lines. 
\documentclass[12pt]{article}
\usepackage{e-jc}

% Please remove all commands that change parameters such as
%    margins or page sizes.

% Packages amssymb and amsthm are already loaded. 
% We recommend these AMS packages:
\usepackage{amsmath,amssymb}
% We recommend the graphicx package for importing figures:
\usepackage{graphicx}
% Use this command to create hyperlinks (optional and recommended):
\usepackage[colorlinks=true,citecolor=black,linkcolor=black,urlcolor=blue]{hyperref}

% Theorem-like environments that are declared in the style file are:
% theorem, lemma, corollary, proposition, fact, observation, claim,
% definition, example, conjecture, open, problem, question, remark, note

%%%%%%%%%%%%%METADATA%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Give the submission and acceptance dates in the format shown.
% The editors will insert the publication date in the third argument.
\dateline{October 23, 2017}{April 9, 2018}{TBD}

% Give one or more subject codes separated by commas.
% Codes are available from http://www.ams.org/mathscinet/freeTools.html
\MSC{05A15, 05A17}

% Uncomment exactly one of the following copyright statements.  Alternatively,
% you can write your own copyright statement subject to the approval of the journal.
% See https://creativecommons.org/licenses/ for a full explanation of the 
% Creative Commons licenses.
%
% We strongly recommend CC BY-ND or CC BY. Both these licenses allow others
% to freely distribute your work while giving credit to you. The difference between
% them is that CC BY-ND only allows distribution unchanged and in whole, while
% CC BY also allows remixing, tweaking and building upon your work.  
%
%    One author:
%\Copyright{The author.}
%\Copyright{The author. Released under the CC BY license (International 4.0).}
%\Copyright{The author. Released under the CC BY-ND license (International 4.0).}
%    More than one author:
%\Copyright{The authors.}
%\Copyright{The authors. Released under the CC BY license (International 4.0).}
\Copyright{  The authors. Released under the CC BY-ND license (International 4.0).}



\newcommand{\M}{{\rm M}}
\newcommand{\al}{{\alpha}}
\newcommand{\be}{{\beta}}
\newcommand{\ga}{{\gamma}}
\newcommand{\del}{{\delta}}
\newcommand{\epf}{{\hfill{$\square$}}}



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% If needed, include a line break (\\) at an appropriate place in the title.
\title{Symmetries of Shamrocks II: Axial Shamrocks}

% Input author, affiliation, address and support information as follows;
% The address should include the country, but does not have to include
%    the street address. Give at least one email address.

\author{Mihai Ciucu\thanks{Research supported in part by NSF grant DMS-1101670 and DMS-1501052.}\\
\small Department of Mathematics\\[-0.8ex]
\small Indiana University \\[-0.8ex] 
\small Bloomington, Indiana 47405\\
\small\tt mciucu@indiana.edu\\
}

%\author{First Author\thanks{Supported by NASA grant ABC123.}\\
%\small Department of Inconsequential Studies\\[-0.8ex]
%\small Solatido College\\[-0.8ex] 
%\small North Kentucky, U.S.A.\\
%\small\tt first.author@dis.solatido.edu\\
%\and
%Some Second Author \qquad  Some Third Author\\
%\small School of Hard Knocks\\[-0.8ex]
%\small University of Western Nowhere\\[-0.8ex]
%\small Somewhere, Australia\\
%\small\tt \{ssa,sta\}@uwn.edu.au}

\begin{document}

\maketitle

% E-JC papers must include an abstract. The abstract should consist of a
% succinct statement of background followed by a listing of the
% principal new results that are to be found in the paper. The abstract
% should be informative, clear, and as complete as possible. Phrases
% like "we investigate..." or "we study..." should be kept to a minimum
% in favor of "we prove that..."  or "we show that...".  Do not
% include equation numbers, unexpanded citations (such as "[23]"), or
% any other references to things in the paper that are not defined in
% the abstract. The abstract may be distributed without the rest of the
% paper so it must be entirely self-contained.  Try to include all words
% and phrases that someone might search for when looking for your paper.

\begin{abstract}
The first paper in this series presented the enumeration of cyclically symmetric, and cyclically symmetric and transpose complementary lozenge tilings of a hexagon with a shamrock removed from its center. In this article we address the transpose complementary case. The results we prove are in fact more general and allow us to give an extension of the symmetric case of the original hexagonal regions with shamrocks removed from their center, to what we call axial shamrocks. For the latter, the transpose complementary case is the only symmetry class besides the one requiring no symmetries. The enumeration of both of these follows from our results.  
\end{abstract}

\section{Introduction: Axial shamrocks}


Inspired by the beauty of MacMahon's formula \cite{MacM} stating that the number of lozenge tilings of a hexagonal region of side-lengths $x$, $y$, $z$, $x$, $y$, $z$ (in cyclic order) on a triangular lattice is
%
\begin{equation}
\prod_{i=1}^x\prod_{j=1}^y\prod_{k=1}^z \frac{i+j+k-1}{i+j+k-2},
\nonumber
\end{equation}
%
and driven by the desire to generalize it,
in \cite{vf} we considered hexagons on the triangular lattice with shamrock-shaped regions (consisting of a central triangle with three additional triangular lobes attached to its vertices) removed from their center, and proved product formulas giving the number of their lozenge tilings. Motivated by the celebrated study 
%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{axial.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{axial_ftp.eps}}
\hfill
}
\vskip-0.1in
\caption{An axial shamrock (left); regions obtained by factorization theorem (right).}
\vskip-0.1in
\label{faa}
\end{figure}
%
of symmetry classes of plane partitions (see \cite{And}\cite{Sta}\cite{Kup}\cite{Ste}\cite{KKZ} and the survey \cite{Bres} for more recent developments), we took up the problem of enumerating their symmetry classes in \cite{symffa}, where we solved two of the six new counting problems arising this way. In this paper we solve a third symmetry class, the one corresponding to transpose complementary plane partitions. We use Kuo's graphical condensation method (see \cite{KuoOne} and \cite{KuoTwo}) to prove our results. In order for this method to work, we have to extend the family of regions we consider, and we end up solving the more general case of axial shamrock regions, defined below.

%Let $x$, $y$, $m$, $a$, $b$ and $k$ be non-negative integers, and consider the hexagon $H$ of side-lengths $x+m$, $y+a+2b$, $y+m$, $x+a+2b$, $y+m$, $y+a+2b$ (clockwise from bottom) on the triangular lattice. In

Consider an arbitrary vertically-symmetric hexagon on the triangular lattice, and remove from it a shamrock-shaped region that is also vertically symmetric about the same symmetry axis as the hexagon, in such a way that the remaining region has the same number of up- and down-pointing unit triangles. We call the resulting region an {\it axial shamrock} (see the picture on the left in Figure~\ref{faa} for an example). Note that in an axial shamrock the removed shamrock-shaped core --- which we will often call simply the {\it core} --- can be freely translated up and down the symmetry axis; for precisely one position among these is the resulting region an $S$-cored hexagon as defined in \cite{vf}. 

%\def\epsfsize#1#2{0.21#1}
%\topinsert
%\twolinetwo{\mypic{axial.eps}}{\mypic{axial_ftp.eps}}
%\medskip
%\centerline{\smc Figure~{\faa}. {\rm An axial shamrock, and regions obtained by factorization theorem.}}
%\endinsert


Axial shamrocks have just one symmetry class besides the case of requiring no symmetry, corresponding to its tilings being vertically symmetric. The transpose complementary case of the original $S$-cored hexagons considered in \cite{vf} and \cite{symffa} follows as a special case of this, as we explain in the next section.

In Section 2 we present the statements of our main results. In order to prove them, we first need to enumerate the tilings of two specific families of regions, which we address in Sections~3 and 4. We prove our main results in Section 5. Section 6 poses an asymptotics problem and presents some interesting data. In Section 7 we present some concluding remarks.


%\bigskip

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{leftregpp.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{rightregpp.eps}}
\hfill
}
\vskip-0.1in
\caption{The regions $R_{x,y,k}(m,a,b)$ (left) and $R'_{x,y,k}(m,a,b)$ (right) for $x=2$, $y=2$, $k=4$, $m=3$, $a=2$, $b=2$.}
\vskip-0.1in
\label{fba}
\end{figure}
%


\section{Statement of main results}

Assuming that the top and bottom sides of the outer hexagon, as well as the sides of the central and top lobe of the shamrock core have even length, the two ``halves'' obtained by applying the factorization theorem of \cite{FT} to the axial shamrock are as shown on the right in Figure \ref{faa}. This leads us to consider the following two families of regions.



%\def\epsfsize#1#2{0.30#1}
%\topinsert
%\twoline{\mypic{leftregpp.eps}}{\mypic{rightregpp.eps}}
%\medskip
%\centerline{\smc Figure~{\fba}. {\rm The regions $R_{x,y,k}(m,a,b)$ (left) and $R'_{x,y,k}(m,a,b)$ (right)}}
%\centerline{for $x=2$, $y=2$, $k=4$, $m=3$, $a=2$, $b=2$.}
%\endinsert

Let $x$, $y$, $k$, $m$, $a$ and $b$ be non-negative integers. We define $R_{x,y,k}(m,a,b)$ to be the region described on the left in Figure \ref{fba}: $x+m$, $m$, $a$ and $b$ are the lengths of the indicated segments, while $y$ and $k$ are the number of ``bumps'' (zig-zags) in the indicated portions of the boundary. 

The region $R'_{x,y,k}(m,a,b)$ is defined analogously, using the picture on the right in Figure~\ref{fba}. 
%Note that the southwestern side has length $y+1$, and there are $k+1$ bumps on the bottom part of the eastern boundary. 
The shaded ovals indicate, throughout this paper, that the corresponding lozenge positions are weighted by $1/2$.

For both the $R$- and the $R'$-regions, the chosen parametrization ensures that the half-shamrock core is inside (or possibly touches the outer boundary of) the half-hexagon.

The main results of this paper are the enumeration of lozenge tilings of these regions, presented in Theorems \ref{tba} and \ref{tbb} below. They allow us to deduce the symmetry classes of axial shamrocks (see Corollaries \ref{tbcorone} and \ref{tbcortwo}).

The formulas we obtain involve the following expressions, known from previous results in the literature, for the enumeration of lozenge tilings of hexagonal regions with a maxi-

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{Proc.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{Procp.eps}}
\hfill
}
\vskip-0.1in
\caption{The regions $H_s(a,b,c)$ (left) and $H'_s(a,b,c)$ (right) for $a=5$, $b=8$, $c=3$.}
\vskip-0.1in
\label{fbaa}
\end{figure}
%

$\!\!\!\!\!\!\!\!$mal staircase shape removed from one of their corners.

For non-negative integers $a$, $b$ and $c$, define $H_s(a,b,c)$ to be the region obtained from the hexagon of side-lengths $a$, $b$, $c$, $a$, $b$, $c$ (in cyclic order, starting with the northwestern side) by removing the maximal staircase shape that fits in its southeastern corner (see the picture on the left in Figure \ref{fbaa} for an illustration). Let $H'_s(a,b,c)$ be the region obtained from $H_s(a,b,c)$ by weighting the lozenge positions that fit in the folds of its southeastern boundary by $1/2$ (this is illustrated on the right in Figure \ref{fbaa}).
For non-negative integers $a$, $b$ and $c$ with $a\leq b+1$, set\footnote{ Under the indicated restrictions for the integers $a$, $b$ and $c$, it is possible for some of the product limits involved in \eqref{eba} and \eqref{ebaa} to be out of order --- more precisely, for the upper limit to be one unit less than the lower limit. In all such instances the corresponding product is taken to be 1.} 
%
\begin{equation}
P(a,b,c):=\prod_{i=1}^a\prod_{j=1}^{b-a+1}\frac{c+i+j-1}{i+j-1}
\prod_{j=b-a+2}^{b-a+i}\frac{2c+i+j-1}{i+j-1}.
\label{eba}
\end{equation}
%
and for non-negative integers $a$, $b$ and $c$ with $a\leq b$, set
%
\begin{equation}
P'(a,b,c):=\frac{1}{2^a}\prod_{i=1}^a\prod_{j=1}^{b-a}\frac{c+i+j-1}{i+j-1}
\prod_{j=b-a+i}^{b}\frac{2c+i+j-1}{i+j-1}.
\label{ebaa}
\end{equation}
%
The lozenge tilings of the regions $H_s(a,b,c)$  were enumerated by Proctor
%\footnote{ The case $a=b+1$ is readily seen to reduce to the case $a=b$ after removing forced lozenges.} 
\cite{Proc}, who proved that\footnote{ Throughout the paper, $\M(G)$ stands for the number of perfect matchings of the graph $G$, or, in case the edges of $G$ carry weights, it denotes the sum over all perfect matchings $\mu$ of $G$ of the weight of $\mu$ (the latter being the product of the weights of the edges that make up $\mu$). Since lozenge tilings of a region $R$ can be identified with perfect matchings of the planar dual graph of $R$, we also use $\M(R)$ to denote the (possibly weighted) count of lozenge tilings of the region $R$ (in this paper, the only tiling weights that occur are 1 and 1/2; the latter appear precisely when a lozenge occupies a tile position weighted by 1/2; all these locations are indicated in the figures by shaded ovals).}, for $a\leq b$,
%
\begin{equation}
\M(H_s(a,b,c))=P(a,b,c).
\label{ebaaa}
\end{equation}
%
One readily sees that \eqref{ebaaa} holds also for $a=b+1$ (upon removing forced lozenges, the case $a=b+1$ reduces to the case $a=b$).

The lozenge tilings of the regions $H'_s(a,b,c)$ were enumerated by Ciucu and Krattenthaler \cite{cutcor}, who showed that for $a\leq b$
%
\begin{equation}
\M(H'_s(a,b,c))=P'(a,b,c).
\label{ebaaaa}
\end{equation}
% 
The formulas we present for the enumeration of the tilings of the regions of type $R$ and $R'$ defined above involve Pochhammer symbols that may have negative index. Recall that, for any integer $k$ (positive or negative) the Pochhammer symbol $(x)_k$ is defined as
%
\begin{equation}
(a)_k=\left\{\begin{matrix} a(a+1)\cdots(a+k-1), & k\geq0\\
\\
                     \dfrac{1}{(a-1)(a-2)\cdots(a-k)}, & k<0\end{matrix}\right.
\label{ebaaaaa}
\end{equation}
% 
All the Pochhammer symbols that occur in this paper are to be understood as being defined by this formula.

%\proclaim{Theorem \tba} For any non-negative values of the parameters for which it is defined, the number of lozenge tilings of the region $R_{x,y,k}(m,a,b)$ is given by
%$$
%\spreadlines{2\jot}
%\align
%&
%\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
%\M(R_{x,y,k}(m,a,b))
%=
%\frac{P(y,y,x)P(k,k,m)}{P(y,y,m)}
%%\\
%%&
%%\times
%P(y-k-b,y-k-b+2m-1,b)
%\\
%&
%\times
%P(y-k-b-2a-1,y-k-b-2a-1,a) 
%\\
%&
%\times
%\prod_{i=1}^a 
%\frac{\left(m+i+\frac12\right)_{y-k-b-2i} (m+i)_{y-k-b-2i+1}}
%     {\left(b+m+i+\frac12\right)_{y-k-b-2i} (b+m+i)_{y-k-b-2i+1}}
%\\
%&
%\times
%\prod_{i=1}^m
%\frac{(x-i+1)_{y-k} (x+i+k+1)_{y-k}}
%     {(m-i+1)_{y-k} (m+i+k+1)_{y-k}}
%\\
%&
%\times
%\prod_{i=1}^a 
%\frac{(m+i+k+1)_{y-2k-2i} \left(m+i+\frac12\right)_{y-2i+1}}
%     {(x+i+k+1)_{y-2k-2i} \left(x+i+\frac12\right)_{y-2i+1}}
%\\
%&
%\times
%\prod_{i=1}^b 
%\frac{(k-i+2)_{y+2m-2k+2i-2} \left(m+a+i+\frac12\right)_{y-2a-2i+1}}
%     {(x-m+k-i+2)_{y+2m-2k+2i-2} \left(x+a+i+\frac12\right)_{y-2a-2i+1}},
%\tag\ebb
%\endalign
%$$
%where $P(a,b,c)$ is given by (\eba).
%
%\endproclaim

The main results of this paper are the following.

\begin{theorem}
\label{tba} 
For any non-negative integer values of the parameters $x,y,k,m,a,b$, 
%for which it is defined, 
the number of lozenge tilings of the region $R_{x,y,k}(m,a,b)$ is given by
%
\begin{align}
&
%\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
\M(R_{x,y,k}(m,a,b))
=
\frac{P(y+2a+2b+k-1,y+2a+2b+k-1,x+m)}
{P(y+2a+2b+k-1,y+2a+2b+k-1,m)}
\nonumber
\\
&\ \ \ \ \ \ \ \ 
\times
P(k+b-1,k+b-1,m) P(y+2a,y+2a+2m-1,b)
%\\
%&
%\times
P(y-1,y-1,a) 
\nonumber
\\
&\ \ \ \ \ \ \ \ 
\times
\prod_{i=1}^a 
\frac{\left(m+i+\frac12\right)_{y+2a-2i} (m+i)_{y+2a-2i+1}}
     {\left(b+m+i+\frac12\right)_{y+2a-2i} (b+m+i)_{y+2a-2i+1}}
\nonumber
\\
&\ \ \ \ \ \ \ \ 
\times
\prod_{i=1}^m
\frac{(x+m-i+1)_{y+2a+b} (x+m+i+k+b)_{y+2a+b}}
     {(m-i+1)_{y+2a+b} (m+i+k+b)_{y+2a+b}}
\nonumber
\\
&\ \ \ \ \ \ \ \ 
\times
\prod_{i=1}^a 
\frac{(m+i+k+b)_{y+2a-k-2i+1} \left(m+i+\frac12\right)_{y+2a+2b+k-2i}}
     {(x+m+i+k+b)_{y+2a-k-2i+1} \left(x+m+i+\frac12\right)_{y+2a+2b+k-2i}}
\nonumber
\\
&\ \ \ \ \ \ \ \ 
\times
\prod_{i=1}^b 
\frac{(k+b-i+1)_{y+2m+2a-k+2i-1} \left(m+a+i+\frac12\right)_{y+2b+k-2i}}
     {(x+k+b-i+1)_{y+2m+2a-k+2i-1} \left(x+m+a+i+\frac12\right)_{y+2b+k-2i}},
\label{ebb}
\end{align}
%
where $P(a,b,c)$ is given by \eqref{eba} when its arguments are non-negative, and is  defined to equal $1$ otherwise.

\end{theorem}



%!!
%{\bfb Double-check formula for new parametrization!}
%!!
%\bigskip
%
%We did --- it checks!! 
%
%See note dated 3/3/17 in file ciuculaptop2/papers/symff/notes.
%
%Recall (see footnote 2) that for a lattice region $R$ for which the lozenge positions are weighted, we denote by $\M(R)$ the weighted count of its lozenge tilings, i.e. the sum of the weights of its lozenge tilings, where the weight of a tiling is the product of the weights of its tiles (in this paper, the only tiling weights that occur are 1 and 1/2; the latter appear precisely when a lozenge occupies a tile position weighted by 1/2; all these locations are indicated in the figures by shaded ovals).

%The formula for the regions $R'_{x,y,k}(m,a,b)$ is in terms of the following counterpart of~(\eba). For non-negative integers $a$, $b$ and $c$ with $a\leq b$, set
%
%$$
%P'(a,b,c):=\frac{1}{2^a}\prod_{i=1}^a\prod_{j=1}^{b-a}\frac{c+i+j-1}{i+j-1}
%\prod_{j=b-a+i}^{b}\frac{2c+i+j-1}{i+j-1}.
%\tag\ebc
%$$
%

\begin{theorem}
\label{tbb} 
For any non-negative integer values of the parameters $x,y,k,m,a,b$, 
%for which it is defined, 
the number of lozenge tilings of the region $R'_{x,y,k}(m,a,b)$ is given by
%
\begin{align}
&
\!\!\!\!
\M(R'_{x,y,k}(m,a,b))
=
2^{y+2a}
\frac{P'(y+2a+2b+k,y+2a+2b+k,x+m)}
{P'(y+2a+2b+k,y+2a+2b+k,m)}
\nonumber
\\
&\ \ \ \ 
\times
P'(k+b,k+b,m)
P'(y+2a,y+2a+2m,b)
%\\
%&
%\times
P'(y,y,a)
\nonumber
\\
&\ \ \ \ 
\times
\prod_{i=1}^a 
\frac{(m+i)_{y+2a-2i+1}\left(m+i-\frac12\right)_{y+2a-2i+2}}
     {(b+m+i)_{y+2a-2i+1}\left(b+m+i-\frac12\right)_{y+2a-2i+2}}
\nonumber
\\
&\ \ \ \ 
\times
\prod_{i=1}^m
\frac{(x+m-i+1)_{y+2a+b}(x+m+i+k+b)_{y+2a+b}}
     {(m-i+1)_{y+2a+b}(m+i+k+b)_{y+2a+b}}
\nonumber
\\
&\ \ \ \ 
\times
\prod_{i=1}^a 
\frac{(m+i+k+b)_{y+2a-k-2i+1}\left(m+i-\frac12\right)_{y+2a+2b+k-2i+2}}
     {(x+m+i+k+b)_{y+2a-k-2i+1}\left(x+m+i-\frac12\right)_{y+2a+2b+k-2i+2}}
\nonumber
\\
&\ \ \ \ 
\times
\prod_{i=1}^b 
\frac{(k+b-i+1)_{y+2m+2a-k+2i-1}\left(m+a+i-\frac12\right)_{y+2b+k-2i+2}}
     {(x+k+b-i+1)_{y+2m+2a-k+2i-1}\left(x+m+a+i-\frac12\right)_{y+2b+k-2i+2}},
\nonumber
\\
\label{ebd}
\end{align}
%
where $P'(a,b,c)$ is given by \eqref{ebaa}.

\end{theorem}


%Note: We checked this also, and it works!! 
%
%See note dated 3/13/17 in file ciuculaptop2/papers/symff/notes.
%

%rs:=proc(x,y,m,a,b,k); p:=2^(y-k-b)*Procp(k+1,k+1,m)*Procp(y+1,y+1,x)/Procp(y+1,y+1,m)*Procp(y-k-b,y-k-b+2*m,b)
%                                                   *Procp(y-k-b-2*a,y-k-b-2*a,a);
%for i to a do
% p:=p*rf(m+i,y-k-b-2*i+1)*rf(m+i-1/2,y-k-b-2*i+2)/
%     (rf(b+m+i,y-k-b-2*i+1)*rf(b+m+i-1/2,y-k-b-2*i+2))
%od;
%for i to m do
% p:=p*rf(x-i+1,y-k)*rf(x+i+k+1,y-k)/
%     (rf(m-i+1,y-k)*rf(m+i+k+1,y-k))
%od;
%for i to a do
% p:=p*rf(m+i+k+1,y-2*k-2*i)*rf(m+i-1/2,y-2*i+3)/
%     (rf(x+i+k+1,y-2*k-2*i)*rf(x+i-1/2,y-2*i+3))
%od;
%for i to b do
% p:=p*rf(m-m+k-i+2,y+2*m-2*k+2*i-2)*rf(m+a+i-1/2,y-2*a-2*i+3)/
%     (rf(x-m+k-i+2,y+2*m-2*k+2*i-2)*rf(x+a+i-1/2,y-2*a-2*i+3))
%od;
%RETURN(p) end:


We can now present the formulas that give the enumeration of the symmetry classes of axial shamrocks.

For non-negative integers $x$, $y$, $k$, $m$, $a$ and $b$, define the axial shamrock $AS_{x,y,k}(m,a,b)$ to be the region obtained from the hexagon of side lengths $2x+2a+2b$, $y+2m+k$, $y+2a+2b+k$, $2x+2m$, $y+2a+2b+k$, $y+2m+k$ (clockwise from top) by removing from it the shamrock-shaped core indicated in Figure \ref{fbc}, placed symmetrically with respect to the vertical symmetry axis of the hexagon. This definition assumes that the central and top lobes of the core have even side-lengths, and that the top and bottom sides of the outer boundary are also even. For the remaining cases there exist enumeration formulas very similar to the ones we present below, but since it captures all the features of the general case, in the interest of clarity we present the formulas just for the case introduced above. 
%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.60\textwidth]{axialp.eps}}
\hfill
}
\vskip-0.1in
\caption{The axial shamrock $AS_{x,y,k}(m,a,b)$ for $x=5$, $y=2$, $k=4$, $m=3$, $a=2$, $b=2$.}
\vskip-0.1in
\label{fbc}
\end{figure}
%


\begin{corollary} 
\label{tbcorone}
For all non-negative integers $x$, $y$, $k$, $m$, $a$ and $b$, the number of lozenge tilings of the axial shamrock $AS_{x,y,k}(m,a,b)$ is given by
%
\begin{equation}
\M(AS_{x,y,k}(m,a,b))=2^{y+k+b}\M(R_{x,y,k}(m,a,b))\M(R'_{x,y,k}(m,a,b)),
\nonumber
\end{equation}
%
where the quantities on the right hand side above are given by the formulas of Theorems~\ref{tba} and \ref{tbb}. 


\end{corollary}

\begin{proof} The statement follows by applying the factorization theorem of \cite{FT} to the axial shamrock $AS_{x,y,k}(m,a,b)$, and using Theorems~\ref{tba} and \ref{tbb}. \end{proof}


\begin{corollary} 
\label{tbcortwo}
For all non-negative integers $x$, $y$, $k$, $m$, $a$ and $b$, the number of vertically symmetric lozenge tilings of the axial shamrock $AS_{x,y,k}(m,a,b)$ is equal to 
%
\begin{equation}
\M_{|}(AS_{x,y,k}(m,a,b))=\M(R_{x,y,k}(m,a,b)),
\nonumber
\end{equation}
%
where the quantity on the right hand side above is given by the formula of Theorem~\ref{tba}. 


\end{corollary}

\begin{proof} This follows directly from Theorem \ref{tba} and the fact that the vertically symmetric lozenge tilings of the axial shamrock $AS_{x,y,k}(m,a,b)$ are in one-to-one correspondence with the tilings of the region $R_{x,y,k}(m,a,b)$. \end{proof}

\medskip
{\bf Remark 1.} The motivating case of the $S$-cored hexagons of \cite{vf} clearly follows as a special case of Corollary \ref{tbcortwo}.

\medskip

Our proofs are based on Kuo's graphical condensation method \cite{KuoOne}. For convenience, we include below the two forms of it that we will employ.

\begin{theorem} \cite[Kuo]{KuoOne}
\label{tbc}
Let $G=(V_1,V_2,E)$ be a plane bipartite graph in which $|V_1|=|V_2|$. Let vertices $\al$, $\be$, $\ga$ and $\del$ appear cyclically on a face of $G$. If $\al,\ga\in V_1$ and $\be,\del\in V_2$, then
%
\begin{align}
\M(G)\M(G-\{\al,\be,\ga,\del\})=
&\M(G-\{\al,\be\})\M(G-\{\ga,\del\})
\nonumber
\\
+&\M(G-\{\al,\del\})\M(G-\{\be,\ga\}).
\label{ebe}
\end{align}
%

\end{theorem}


%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{perregmon.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{perregmonp.eps}}
\hfill
}
\vskip-0.1in
\caption{Applying Kuo condensation for periscope regions.}
\vskip-0.1in
\label{fca}
\end{figure}
%



\begin{theorem} \cite[Kuo]{KuoOne}
\label{tbd}
Let $G=(V_1,V_2,E)$ be a plane bipartite graph in which $|V_1|=|V_2|+1$. Let vertices $\al$, $\be$, $\ga$ and $\del$ appear cyclically on a face of $G$. If $\al,\be,\ga\in V_1$ and $\del\in V_2$, then
%
\begin{align}
\M(G-\be)\M(G-\{\al,\ga,\del\})=
&\M(G-\al)\M(G-\{\be,\ga,\del\})
\nonumber
\\
+&\M(G-\ga)\M(G-\{\al,\be,\del\}).
\label{ebf}
\end{align}
%

\end{theorem}


%








%* enumeration of axial shamrocks

%\newpage

\section{Periscope regions $\Gamma$ and $\Gamma'$}


For any non-negative integers $x,y,m,a$, define the periscope region $\Gamma_{x,y}(m,a)$ to be the region described in the picture on the left in Figure \ref{fca}. The periscope region $\Gamma'_{x,y}(m,a)$ is obtained from it by augmenting it with $y$ unit rhombi that fit in the folds of the top part of its right boundary, as illustrated in the picture on the right in Figure \ref{fca} (the shaded ovals indicate lozenge positions weighted by $1/2$). 



The purpose of this section is to provide enumeration formulas for the lozenge tilings of the regions $\Gamma_{x,y}(m,a)$ and $\Gamma'_{x,y}(m,a)$.

\begin{theorem}
\label{tca}
For any non-negative integers $x$, $y$, $m$ and $a$, the number of lozenge tilings of the periscope region $\Gamma_{x,y}(m,a)$ is given by
%$$
%\spreadlines{2\jot}
%\align
%&
%\M(\Gamma_{x,y}(m,a))
%=
%P(y-2a-1,y-2a-1,a)P(y,y+m-1,x)
%\\
%&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  
%\times
%\prod_{i=1}^a 
%\frac{\left(\frac{m}{2}+i+\frac12\right)_{y-2i} \left(\frac{m}{2}+i\right)_{y-2i+1}}
%     {\left(x+\frac{m}{2}+i+\frac12\right)_{y-2i} \left(x+\frac{m}{2}+i\right)_{y-2i+1}},
%\tag\eca
%\endalign
%$$
%
\begin{align}
&
\M(\Gamma_{x,y}(m,a))
=
P(y-1,y-1,a)P(y+2a,y+2a+m-1,x)
\nonumber
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  
\times
\prod_{i=1}^a 
\frac{\left(\frac{m}{2}+i+\frac12\right)_{y+2a-2i} \left(\frac{m}{2}+i\right)_{y+2a-2i+1}}
     {\left(x+\frac{m}{2}+i+\frac12\right)_{y+2a-2i} \left(x+\frac{m}{2}+i\right)_{y+2a-2i+1}},
\label{eca}
\end{align}
%
where $P(a,b,c)$ is given by \eqref{eba} when all its arguments are non-negative, and is equal to~$1$ otherwise.

\end{theorem}

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.26\textwidth]{perrec1.eps}}
\hfill
{\includegraphics[width=0.21\textwidth]{perrec2.eps}}
\hfill
}
%\end{figure}
%
%\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.21\textwidth]{perrec3.eps}}
\hfill
{\includegraphics[width=0.21\textwidth]{perrec4.eps}}
\hfill
}
%\end{figure}
%
%\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.21\textwidth]{perrec5.eps}}
\hfill
{\includegraphics[width=0.21\textwidth]{perrec6.eps}}
\hfill
}
%\vskip-0.1in
\caption{The recurrence for the periscope regions $\Gamma_{x,y}(m,a)$.}
\vskip-0.1in
\label{fcb}
\end{figure}
%

\begin{proof} Apply the version of Kuo condensation stated in Theorem \ref{tbc}, with monomers $\al,\be,\ga,\del$ chosen as indicated\footnote{ More precisely, we view the region $\Gamma_{x,y}(m,a)$ as a planar graph whose vertices are the unit triangles inside it, and whose edges connect two vertices precisely if the corresponding unit triangles share an edge. Then we apply Theorem \ref{tbc} to this graph, with vertices $\al,\be,\ga,\del$ corresponding to the indicated four unit triangles.} on the left in Figure \ref{fca}. We obtain (see Figure \ref{fcb})
%
\begin{align}
&\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
\M(\Gamma_{x,y}(m,a))\M(\Gamma_{x,y-1}(m-1,a))=
\nonumber
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
\M(\Gamma_{x-1,y}(m-1,a))\M(\Gamma_{x+1,y-1}(m-1,a))
\nonumber
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
+
\M(\Gamma_{x,y-1}(m,a))\M(\Gamma_{x,y}(m-1,a)),
\label{ecb}
\end{align}
%
for all integers $x,y,m\geq1$ and $a\geq0$ (so that all the involved regions are defined).

One readily sees that the $\Gamma$-regions can always be tiled by lozenges, so in particular we can deduce from \eqref{eda} that 
%
\begin{align}
\M(\Gamma_{x,y}(m,a))=
\frac{1}{\M(\Gamma_{x,y-1}(m-1,a))}
\left\{\right.
&\M(\Gamma_{x-1,y}(m-1,a))\M(\Gamma_{x+1,y-1}(m-1,a))
\nonumber
\\
+
&\left.\M(\Gamma_{x,y-1}(m,a))\M(\Gamma_{x,y}(m-1,a))\right\},
\label{ecc}
\end{align}
%
for all $x,y,m\geq1$.

We prove \eqref{eca} by induction on $x+y+m$. Suppose that \eqref{eca} holds for all non-negative integers $x,y,m,a$ with $x+y+m<l$, and let $x,y,m\geq1$ and $a\geq0$ be integers with $x+y+m=l$. Express $\M(\Gamma_{x,y}(m,a))$ as the expression on the right hand side of \eqref{ecc}. By the induction hypothesis, the five quantities on the right hand side of \eqref{ecc} are given by formula \eqref{eca}. One readily checks that the resulting formula for $\M(\Gamma_{x,y}(m,a))$ agrees with the one given by \eqref{eca}. This proves the induction step.

The base cases are $x=0$, $y=0$, and $m=0$. As illustrated in Figure \ref{fcc}, the first and third of these reduce, after removing forced lozenges, to enumerating tilings of special cases of the regions $H_s(a,b,c)$ described at the end of Section 2. For $x=0$, the resulting formula visibly agrees with the special case $x=0$ of \eqref{eca}; a straightforward calculation checks this also for $m=0$. For $y=0$, the region $\Gamma_{x,y}(m,a)$ becomes the hexagon of side-lengths $x$, $2a$, $m$, $x$, $2a$, $m$ (in cyclic order), and one readily verifies that \eqref{eca} follows from MacMahon's original enumeration \cite{MacM} of boxed plane partitions (which are well-known to be in one-to-one correspondence with lozenge tilings of centrally symmetric hexagons on the triangular lattice; see e.g. \cite{DT}). This concludes our proof. \end{proof}


%Therefore we can proceed by induction on $x+y+m$. ``Basic'' base cases are $x+y+m=0$ and $x+y+m=1$, and boundary base cases are $x=0$, $y=0$, and $m=0$. The latter cover also the former.

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.27\textwidth]{perxis0.eps}}
\hfill
{\includegraphics[width=0.27\textwidth]{peryis0.eps}}
\hfill
{\includegraphics[width=0.27\textwidth]{permis0.eps}}
\hfill
}
\vskip-0.1in
\caption{The base cases $x=0$ (left), $y=0$ (center)  and $m=0$ (right) for the regions $\Gamma_{x,y}(m,a)$.}
\vskip-0.1in
\label{fcc}
\end{figure}
%


The enumeration of tilings of the periscope regions $\Gamma'_{x,y}(m,a)$ is proved by a completely analogous inductive argument, with the only difference that equation~\eqref{ebaa} is used for the base case instead of \eqref{eba}. We obtain the following result.

\begin{theorem}
\label{tcb}
For any non-negative integers $x$, $y$, $m$ and $a$, the number of lozenge tilings of the periscope region $\Gamma'_{x,y}(m,a)$ is given by
%$$
%\spreadlines{2\jot}
%\align
%&
%\M(\Gamma_{x,y}(m,a))
%=
%P(y-2a-1,y-2a-1,a)P(y,y+m-1,x)
%\\
%&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  
%\times
%\prod_{i=1}^a 
%\frac{\left(\frac{m}{2}+i+\frac12\right)_{y-2i} \left(\frac{m}{2}+i\right)_{y-2i+1}}
%     {\left(x+\frac{m}{2}+i+\frac12\right)_{y-2i} \left(x+\frac{m}{2}+i\right)_{y-2i+1}},
%\tag\eca
%\endalign
%$$
%
\begin{align}
&
\M(\Gamma'_{x,y}(m,a))
=
2^{y}P'(y,y,a)P'(y+2a,y+2a+m,x)
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  
\times
\prod_{i=1}^a 
\frac{\left(\frac{m}{2}+i\right)_{y+2a-2i+1} \left(\frac{m}{2}+i-\frac12\right)_{y+2a-2i+2}}
{\left(x+\frac{m}{2}+i\right)_{y+2a-2i+1} \left(x+\frac{m}{2}+i-\frac12\right)_{y+2a-2i+2}},
\label{eca}
\end{align}
%
where $P'(a,b,c)$ is given by \eqref{ebaa}.

\end{theorem}



%\newpage

\section{Sigma regions $\Sigma$ and $\Sigma'$}

Define the sigma region $\Sigma_{x,y,k}(a,m)$ to be $R_{x,y,k}(m,a,0)$, the $b=0$ specialization of the region $R_{x,y,k}(m,a,b)$ described in Section 2. Similarly, let $\Sigma'_{x,y,k}(a,m)$ be $R'_{x,y,k}(m,a,0)$, the $b=0$ specialization of the region $R'_{x,y,k}(m,a,b)$. In this section we enumerate the tilings of these regions. They will serve as base cases for the inductive proofs of Theorems~\ref{tba} and \ref{tbb}, which we present in the next section.

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{sigmaregmon.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{sigmaregmonp.eps}}
\hfill
}
\vskip-0.1in
\caption{Applying Kuo condensation for the sigma regions.}
\vskip-0.1in
\label{fda}
\end{figure}
%



\begin{theorem}
\label{tda}
For any non-negative integers $x$, $y$, $k$, $m$ and $a$, the number of lozenge tilings of the sigma region $\Sigma_{x,y,k}(m,a)$ is given by
%
\begin{align}
&
\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
\M(\Sigma_{x,y,k}(m,a))
=
P(k-1,k-1,m) 
P(y-1,y-1,a)
\nonumber 
\\
&\ \ \ \ \ \ \ \ 
\times
\frac{P(y+2a+k-1,y+2a+k-1,x+m)}
{P(y+2a+k-1,y+2a+k-1,m)}
\nonumber 
\\
&\ \ \ \ \ \ \ \ 
\times
\prod_{i=1}^m
\frac{(x+m-i+1)_{y+2a} (x+m+i+k)_{y+2a}}
     {(m-i+1)_{y+2a} (m+i+k)_{y+2a}}
\nonumber 
\\
&\ \ \ \ \ \ \ \ 
\times
\prod_{i=1}^a 
\frac{(m+i+k)_{y+2a-k-2i+1} \left(m+i+\frac12\right)_{y+2a+k-2i}}
     {(x+m+i+k)_{y+2a-k-2i+1} \left(x+m+i+\frac12\right)_{y+2a+k-2i}}
\label{eda}
\end{align}
%
%$$
%\spreadlines{2\jot}
%\align
%&
%\!\!\!\!\!
%\M(\Sigma_{x,y,k}(m,a))
%=
%\frac{P(y,y,x)P(k,k,m)}{P(y,y,m)}
%P(y-k-2a-1,y-k-2a-1,a) 
%\\ \ \ 
%&
%\times
%\prod_{i=1}^m
%\frac{(x-i+1)_{y-k} (x+i+k+1)_{y-k}}
%     {(m-i+1)_{y-k} (m+i+k+1)_{y-k}}
%\prod_{i=1}^a 
%\frac{(m+i+k+1)_{y-2k-2i} \left(m+i+\frac12\right)_{y-2i+1}}
%     {(x+i+k+1)_{y-2k-2i} \left(x+i+\frac12\right)_{y-2i+1}},
%\\
%\tag\eda
%\endalign
%$$
where $P(a,b,c)$ is given by \eqref{eba} when all its arguments are non-negative, and is equal to~$1$ otherwise.

\end{theorem}


\begin{proof} By Kuo condensation with monomers picked as in the picture on the left in Figure \ref{fda}, we get (see Figure \ref{fdb})
%
\begin{align}
&\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
\M(\Sigma_{x,y,k}(m,a))\M(\Sigma_{x,y-1,k-1}(m,a))=
\nonumber
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
\M(\Sigma_{x-1,y,k}(m,a))\M(\Sigma_{x+1,y-1,k-1}(m,a))
\nonumber
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
+
\M(\Sigma_{x,y-1,k}(m,a))\M(\Sigma_{x,y,k-1}(m,a)).
\label{edb}
\end{align}
%
We proceed by induction on $x+y+k$. Suppose $x,y,k\geq1$ (so that all the regions involved in \eqref{edb} are defined), and assume that \eqref{eda} holds for all the $\Sigma$-regions whose $x$-, $y$- and $k$-parameters add up to strictly less than $x+y+k$. By \eqref{edb}, $\M(\Sigma_{x,y,k}(m,a))$ is expressed in terms of tiling counts of $\Sigma$-regions whose sums of their $x$-, $y$- and $k$-parameters are strictly less than $x+y+k$, the latter being --- by the induction hypothesis --- given by formula \eqref{eda}. It is routine to verify that the resulting formula for $\M(\Sigma_{x,y,k}(m,a))$ agrees with the one provided by \eqref{eda}. This proves the induction step.

There are three base cases: $x=0$, $y=0$ and $k=0$. (There are two different sources for base cases: since \eqref{edb} relates $\M(\Sigma_{x,y,k}(m,a))$ to $\Sigma$-regions with $(x+y+k)$-sum one or two units less, we have the base cases $x+y+k=0$ and $x+y+k=1$; in addition we have the three base cases $x=0$, $y=0$ and $k=0$, due to the fact that when any of these holds, the regions involved in \eqref{edb} are not all defined. Clearly, the latter three cases cover also the former two.)

When $x=0$, the $\Sigma$-region looks as pictured in Figure \ref{fdc}. It is easy to see that the two half-hexagons cut out by the indicated thick dashed lines
must be internally tiled in every tiling of such a region (indeed, if there were for instance a tiling of the region containing a lozenge that straddles the bottom dashed line, then starting from its bottom edge there would be a path of lozenges in this tiling that ends on a unit segment of the base; but then the $m$ paths of lozenges starting downward from the unit segments of the inner seg-

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.17\textwidth]{sigmarec1.eps}}
\hfill
{\includegraphics[width=0.15\textwidth]{sigmarec2.eps}}
\hfill
}
\centerline{
\hfill
{\includegraphics[width=0.15\textwidth]{sigmarec3.eps}}
\hfill
{\includegraphics[width=0.15\textwidth]{sigmarec4.eps}}
\hfill
}
\centerline{
\hfill
{\includegraphics[width=0.15\textwidth]{sigmarec5.eps}}
\hfill
{\includegraphics[width=0.15\textwidth]{sigmarec6.eps}}
\hfill
}
%\vskip-0.1in
\caption{The recurrence for the sigma regions $\Sigma_{x,y,k}(m,a)$.}
\vskip-0.1in
\label{fdb}
\end{figure}
%

$\!\!\!\!\!\!\!\!$ment of length $m$ on the right boundary would not have enough room to end on the unit segments of the base, a contradiction). Since the tiling is obviously forced on the remaining parallelogram, it follows by \eqref{ebaaa} that 


%
\begin{equation}
\M(\Sigma_{0,y,k}(m,a))=P(k-1,k-1,m)P(y-1,y-1,a),
\nonumber
\end{equation}
%
which is readily seen (cf. \eqref{eba})  to agree with the $x=0$ specialization of the right hand side of \eqref{eda}.
%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{sigmaxis0.eps}}
\hfill
}
\vskip-0.1in
\caption{The base case $x=0$ for the sigma regions $\Sigma_{x,y,k}(m,a)$.}
\vskip-0.1in
\label{fdc}
\end{figure}
% 
On the other hand, if $y=0$ or $k=0$, the $\Sigma$-region clearly reduces to a periscope region, and its enumeration follows therefore from Theorem \ref{tca}. It is straightforward to verify that the resulting formulas agree with the $y=0$ and $k=0$ specializations of the right hand side of \eqref{eda}, respectively. This completes the proof.
\end{proof}

The same arguments prove the following result.

\begin{theorem}
\label{tdb}
For any non-negative integers $x$, $y$, $k$, $m$ and $a$, the number of lozenge tilings of the sigma region $\Sigma_{x,y,k}(m,a)$ is given by
%
\begin{align}
&\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
\M(R'_{x,y,k}(m,a,b))
=
2^{y+2a}
P'(k,k,m)
P'(y,y,a)
\frac{P'(y+2a+k,y+2a+k,x+m)}
{P'(y+2a+k,y+2a+k,m)}
\nonumber
\\
&\ \ \ \ 
\times
\prod_{i=1}^m
\frac{(x+m-i+1)_{y+2a}(x+m+i+k)_{y+2a}}
     {(m-i+1)_{y+2a}(m+i+k)_{y+2a}}
\nonumber
\\
&\ \ \ \ 
\times
\prod_{i=1}^a 
\frac{(m+i+k)_{y+2a-k-2i+1}\left(m+i-\frac12\right)_{y+2a+k-2i+2}}
     {(x+m+i+k)_{y+2a-k-2i+1}\left(x+m+i-\frac12\right)_{y+2a+k-2i+2}},
\nonumber
\\
\label{ebd}
\end{align}
%
where $P'(a,b,c)$ is given by \eqref{ebaa}.


\end{theorem}







%xs\newpage

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{leftregmon.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{leftregmonp.eps}}
\hfill
}
\vskip-0.1in
\caption{Choosing the monomers in the regions $R_{x,y,k}(m,a,b)$ and $R'_{x,y,k}(m,a,b)$.}
\vskip-0.1in
\label{fea}
\end{figure}
% 


\section{Proofs of Theorems 2.1--2.2}

We prove our results by induction, relying on Kuo's graphical condensation method \cite{KuoOne} at the induction step, and on the results of Sections 3 and 4 for the base cases.

We present the details for the proof of Theorem \ref{tba}. The proof of Theorem \ref{tbb} is obtained by the same arguments, with exactly the same arguments at the induction step, and minimal modifications for the base cases.

Consider the region obtained from $R_{x,y,k}(m,a,b)$ by reducing by one unit the side-length of the lobe of side $b$, as illustrated in the picture on the left in Figure \ref{fea}, and choose the four monomers in Kuo's graphical condensation (the version of Theorem \ref{tbd}) as indicated in that figure. 

Each of the six regions involved in the resulting recurrence is readily seen to become, after conveniently removing some of the forced lozenges (operation which preserves the number of lozenge tilings), one of our regions of type $R$. More precisely, as can be seen from Figure \ref{feb}, we obtain the recurrence
%
\begin{align}
&\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
\M(R_{x,y,k}(m,a,b))\M(R_{x,y,k}(m,a,b-1))=
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
\M(R_{x,y,k+1}(m,a,b-1))\M(R_{x,y,k-1}(m,a,b))
\\
&\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
+
\M(R_{x+1,y,k}(m,a,b-1))\M(R_{x-1,y,k}(m,a,b)).
\label{eea}
\end{align}
%
We claim that this recurrence is valid for all integers $x,k,b\geq1$ and $y,m,a\geq0$ (clearly necessary conditions for the involved regions to be defined).

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.17\textwidth]{lrecp1.eps}}
\hfill
{\includegraphics[width=0.17\textwidth]{lrecp2.eps}}
\hfill
}
\centerline{
\hfill
{\includegraphics[width=0.17\textwidth]{lrecp3.eps}}
\hfill
{\includegraphics[width=0.17\textwidth]{lrecp4.eps}}
\hfill
}
\centerline{
\hfill
{\includegraphics[width=0.17\textwidth]{lrecp5.eps}}
\hfill
{\includegraphics[width=0.17\textwidth]{lrecp6.eps}}
\hfill
}
%\vskip-0.1in
\caption{The recurrence for the regions $R_{x,y,k}(m,a,b)$.}
\vskip-0.1in
\label{feb}
\end{figure}
%

Indeed, the indicated choices for the monomers (see the picture on the left in Figure~\ref{fea}) are possible (in other words, there is room for the indicated monomers) provided $x+m\geq1$, $x+a+b\geq1$, $b\geq1$ and $b+k\geq1$. In turn, all these inequalities clearly hold if $x,k,b\geq1$.
Furthermore, the pattern of forced tiles indicated in Figure \ref{feb} is valid provided $x,k,b\geq1$. Finally, for instance the fact that the lobe of side length $b$ fits inside the region indicated by the thick contour in the center right picture in Figure \ref{feb} does not need to be checked, being guaranteed by the very fact that the latter is a region of type $R$ (namely, $R_{x,y,k-1}(m,a,b)$; see the observation in the fourth paragraph of Section 2). This completes the proof of~\eqref{eea} for all $y,m,a\geq0$ and $x,k,b\geq1$. The base cases of our induction will be $x=0$, $k=0$ and $b=0$.

%Look at top right region..
%
%The picture tacitly assumes that $x+m\geq1, b\geq1, k+b\geq1$, top side $ = x+a+b\geq1$. These are covered by the base cases $x=0$, $b=0$.
%
%The existence of the other pictures in the form shown (i.e., that forcing occurs in the pattern shown) implies further that $y+2a+2b+k\geq 2$, top left $=y+2m+k\geq1$, that ``pulling in `` SW corner as in center right picture leaves room for the $b$-lobe, and  that ``pulling in `` NW corner as in top right picture leaves room for the $2a$-half-lobe.
%
%All these last four are covered by the $b=0$ and $k=0$ base cases (first, clear; for last three, see page 3, 2/27/17, of earlier printout of draft of this paper).
%
%Everything is thus covered by the base cases $x=0$, $k=0$, $b=0$. The first two of these amount to periscope regions (enumerated Section 3), the last one to sigma regions (enumerated in Section 4). It is readily seen that the formulas we get for these cases agree with those given by Theorem {\fba} (similarly for Theorem {\fbb).

%These will be our base cases.


{\it Proof of Theorem \ref{tba}}. 
%(We prove (\ebb) by induction on $x$ (or $k$?).)
%
%
%\bigskip
Use \eqref{eea} to express $\M(R_{x,y,k}(m,a,b))$ in terms of the other five quantities in that equation. We want to regard the resulting expression as a recurrence  

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{leftregxis0.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{leftregkis0.eps}}
\hfill
}
\vskip-0.1in
\caption{The base cases $x=0$ (left) and $k=0$ (right) for the regions $R_{x,y,k}(m,a,b)$.}
\vskip-0.1in
\label{fec}
\end{figure}
% 

$\!\!\!\!\!\!\!\!\!$for $\M(R_{x,y,k}(m,a,b))$. Note that the sum of the $x$-, $k$- and $b$-parameters --- the most natural quantity to consider --- is not strictly less than $x+k+b$ for all five $R$-regions in \eqref{eea} other than $R_{x,y,k}(m,a,b)$. However, for each of these five regions, adding the $x$-, $k$- and {\it twice} the $b$-parameter results in a quantity that is either 1 or 2 units less than $x+k+2b$. We proceed therefore by induction on $x+k+2b$.

%Do induction on $x+k+2b$. That does strictly go down for all last five regions, compared to first region --- by $-2$, $-1$, $-1$, $-1$, and $-1$, respectively (kind of like in Dodgeson condensation!)

Let $x,k,b\geq1$, and assume that \eqref{ebb} holds for all $R$-regions in which the sum between the $x$-, $y$- and twice the $b$-parameter is strictly less than $x+k+2b$. Then \eqref{eea} yields, by the induction hypothesis, an explicit expression for $\M(R_{x,y,k}(m,a,b))$ obtained by applying formula \eqref{ebb} to the other five quantities in \eqref{eea}. A straightforward calculation verifies that this expression agrees with the right hand side of \eqref{ebb}. This proves the induction step.

The base cases of the induction are $x=0$, $k=0$, and $b=0$. For $b=0$, the region $R_{x,y,k}(m,a,b)$ becomes the sigma region $\Sigma_{x,y,k}(m,a)$, and its tiling count follows by Theorem \ref{tda}. The resulting expression visibly agrees with the $b=0$ specialization of the right hand side of \eqref{ebb}.

%Two ``basic'' base cases: $x+k+2b=0$, $x+k+2b=1$, and three ``boundary'' base cases: $x=0$, $k=0$, and $b=0$. Note that $x+k+2b=0$ implies $x=k=0$, and $x+k+2b=1$ implies $b=0$ and either $x=0, k=1$ or $x=1, k=0$. So the boundary base cases suffice.

%\bigskip
%(Tricky this way --- x+k+b does not go down all the time!
%
%\bigskip
%Try: Use original 4 monomers; need ``periscope with a bay'' for base case; deduce formula for that by setting y=0 in T2.1-2.2; enough to have just one additional section --- periscope with bay. However: latter may have same problem in proving it! Still: could get around that problem as above (induction on $x+k+2b$), and then we'd still have the advantage of needing just one additional section, not two. Weigh carefully the pros and cons and decide which way to go.
%
%NOTE: It seems that in proving periscope with a bay, at least if we choose same monomers as above, b=0 will be a base case. But then we again need an additional section!
%
%So perhaps stay with current approach.
%)

%The base case $b=0$ clearly reduces to the sigma regions of Section 4 (and formula for sigma regions there is visibly $b=0$ specialization of (\ebb)).

For $x=0$, the $R$-region looks as illustrated in the picture on the left in Figure \ref{fec}. Using the same argument as in the proof of Theorem \ref{tda}, one readily sees that the half-hexagon cut out by the bottom thick dashed line in the picture on the left in Figure \ref{fec} must be internally tiled. This in turn forces the indicated parallelogram to be tiled internally. It follows that the periscope region delimited by the top thick dashed line is also internally tiled, and therefore
%
\begin{equation}
\M(R_{0,y,k}(m,a,b))=P(b+k-1,b+k-1,m)\M(\Gamma_{b,y+2a}(2m,a)).
\label{eeb}
\end{equation}
%
One readily checks that the formula resulting from \eqref{eeb} using \eqref{eba} and Theorem \ref{tca} agrees with the $x=0$ specialization of the right hand side of \eqref{ebb}.
The base case $k=0$ follows analogously, using that, as seen from the picture on the right in Figure \ref{fec},
%
\begin{equation}
\M(R_{x,y,0}(m,a,b))=P(b-1,b-1,m)\M(\Gamma_{b+x,y+2a}(2m,a)).
\label{eec}
\end{equation}
%
%The base case $k=0$ reduces to a Proc region and a periscope region (see picture on the left in Figure {\fec}). The same is true for the base case $k=0$ (see picture on the right in Figure {\fec}).
This completes the proof. \epf

The proof of Theorem \ref{tbb} is perfectly analogous. One applies Kuo condensation as indicated in the picture on the right in Figure \ref{fea} (note that the choice of the monomers is precisely the same as above). The previous proof of the induction step carries over with no change. The only difference is that at the base cases we now use the counterparts Theorem \ref{tdb}, Theorem \ref{tcb} and \eqref{ebaaaa} of Theorem \ref{tda}, Theorem \ref{tca} and \eqref{ebaaa}, respectively.


\section{The core position that maximizes the number of tilings}

Since the shamrock hole can glide up and down along the vertical symmetry axis, a natural question is for which position is the number of tilings maximum.

We can gain some intuition about this by recalling the electrostatic interpretation of random tilings with gaps worked out in our previous work \cite{sc}\cite{ec}\cite{ov}\cite{gd}\cite{fin}: If gaps are viewed as electrical charges of magnitude given by the difference between the number of up- and down-pointing unit triangles in them, then for large separations between the gaps, the probability of any position of the gaps is proportional to the exponential of the negative 2D Coulomb energy of the system.

We only have one gap in our current situation (the central shamrock), of charge $2m-2a-2b$, but the set-up is such that this gap feels the interaction with the boundary of the dimer system --- more precisely, since it is allowed to glide along the symmetry axis, the interaction with the top and bottom sides of the outer hexagon. Since both these boundaries are constrained lattice lines, for $2m-2a-2b\neq0$ we expect (see \cite{fin} and~\cite{ranglep}) an exponentially strong interaction with the boundary --- overwhelming attraction to the bottom side for $2m-2a-2b>0$, and overwhelming attraction to the top side for $2m-2a-2b<0$.

Therefore, the interesting case to consider is that when the central shamrock has charge zero, i.e. $2m-2a-2b=0$. 

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{ais0p.eps}}
\hfill
{\includegraphics[width=0.40\textwidth]{bis0.eps}}
\hfill
}
\vskip-0.1in
\caption{The case $a=0$ (left) and the case $b=0$ (right).}
\vskip-0.1in
\label{ffa}
\end{figure}
% 

If $a=0$, then the axial shamrock looks as pictured on the left in Figure \ref{ffa}. Since the top lobe has then positive charge (as the $a$-lobe disappears) and the two bottom lobes have negative charge, we expect (and this is confirmed by numerical examples) that the central shamrock will be repelled both by the top and by the bottom sides of the outer hexagon, and thus that the position with maximum number of tilings will be somewhere in the interior. At the other extreme, when $a=m$ and $b=0$, the axial shamrock looks as pictured on the right in Figure \ref{ffa}, and by the same heuristics we expect (and this is again confirmed by numerical examples) that both the top and the bottom side will attract the central shamrock, and that the maximum will occur when the central shamrock touches one of these sides. Furthermore, for $a<<b$ we expect the former kind of behavior, while for $a>>b$ the latter. The interesting question is to describe the point at which one behavior changes into the other.

Numerical examples suggest that there exists a critical ratio $r=1.5...$ so that for large enough outer hexagons, if $a$ and $b$ are fixed\footnote{ Since we are discussing the case of the neutral central shamrock, fixing $a$ and $b$ fixes $m$ as well, as $m=a+b$.} and $b/a<r$ then the first kind of behavior occurs, while if $b/a>r$ then the second kind holds. This could be verified using our explicit formula for $\M(AS_{x,y,k}(m,a,b))$ given in Corollary \ref{tbcorone}, but the calculations are somewhat complicated and we have not carried them out. 

An interesting variation is to consider the same question, but to scale the central shamrock together with the outer boundary. 

%For $2a+2b=2m$, which height of the shamrock maximizes the number of tilings? (Here all parameters kept constant, except $y$ and $k$ vary so that $y+k=const$)
%
%NOTE: We can study this by taking the ratio of $(y,k)$ case to that of the $(y-1, k+1)$ case!!
%
%If we're lucky and this sequence is unimodal, then at the maximum of this sequence we have the maximum of the number of tilings!!
%
%Look into this. First thing: Get a simple expression for the ratio of the two cases mentioned above, using formulas from Theorems 2.1 and 2.2 --- there should be great simplification!!


\section{Concluding remarks}


Upon looking at the region $R_{x,y,k}(m,a,b)$ on the left in Figure \ref{fba} (or even directly at Figure~\ref{faa}) one may notice the following lack of symmetry: There is a triangular lobe attached to the lower corner of the $\Sigma$-shaped portion of the eastern boundary, but no lobe attached to the upper corner. Surely, one might say, the more general region obtained including one lobe there as well in order to restore this symmetry (and the similar extension of the region $R'_{x,y,k}(m,a,b)$) should also have its number of lozenge tilings given by a simple product formula. This would imply, by the factorization theorem of \cite{FT}, that tilings of regions of the type illustrated in Figure \ref{fga} are also enumerated by such a product formula.

%
\begin{figure}[h]
\centerline{
\hfill
{\includegraphics[width=0.40\textwidth]{axialg.eps}}
\hfill
}
\vskip-0.1in
\caption{A possible generalization of axial shamrocks.}
\vskip-0.1in
\label{fga}
\end{figure}
% 

However, this doesn't seem to be the case. Numerical examples show that relatively large primes occur in the factorization of the number of tilings of such regions, which does not happen in a product formula where all factors are linear in the parameters (this happens even in the most symmetric case, when the outer hexagon is regular, and the removed structure is symmetric about the center). It would be interesting to find a formula for the number of tilings of these more general regions.




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\begin{figure}[!h]
%  \begin{center}
%    % use \includegraphics to import figures 
%    % \includegraphics{filename}
%  \end{center}
%  \caption{\label{fig:InformativeFigure} Here is an informative
%    figure.}
%\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% \bibliographystyle{plain} 
% \bibliography{myBibFile} 
% If you use BibTeX to create a bibliography
% then copy and past the contents of your .bbl file into your .tex file

\begin{thebibliography}{10}


\bibitem{And} G. E. Andrews.  
\newblock Plane partitions (III): The weak Macdonald conjecture.  
\newblock {\em Invent. Math.}, 53:193--225, 1979.


\bibitem{Bres}D. M. Bressoud.  
\newblock Proofs and confirmations --- The story of the alternating sign matrix conjecture.  
\newblock Cambridge University Press, Cambridge, 2009.


\bibitem{FT} M. Ciucu.  
\newblock Enumeration of perfect matchings in graphs with reflective symmetry.  
\newblock {\em J. Combin. Theory Ser. A}, 77:67--97, 1997.


\bibitem{cutcor} M. Ciucu and C. Krattenthaler.  
\newblock Enumeration of lozenge tilings of hexagons with cut off corners.  
\newblock {\em J. Combin. Theory Ser. A}, 100:201--231, 2002.


\bibitem{sc} M. Ciucu.  
\newblock A random tiling model for two dimensional electrostatics.  
\newblock {\em Mem. Amer. Math. Soc. }, 178:1--106, no. 839, 2005.


\bibitem{ec} M. Ciucu.  
\newblock The scaling limit of the correlation of holes on the triangular lattice with periodic boundary conditions.  
\newblock {\em Mem. Amer. Math. Soc.}, 199:1--100, no. 935, 2009.

\bibitem{ov} M. Ciucu.  
\newblock Dimer packings with gaps and electrostatics.  
\newblock {\em Proc. Nat. Acad. Sci. USA}, 105:2766--2772, no. 8, 2008.

\bibitem{vf} M. Ciucu and C. Krattenthaler.  
\newblock A dual of MacMahon's theorem on plane partitions.  
\newblock {\em Proc. Natl. Acad. Sci. USA}, 110:4518--4523, 2013.

\bibitem{gd} M. Ciucu.  
\newblock The interaction of collinear gaps of arbitrary charge in a two dimensional dimer system.  
\newblock {\em Comm. Math. Phys.}, 330:1115--1153, 2014.

\bibitem{fin} M. Ciucu.  
\newblock Macroscopically separated gaps in dimer coverings of Aztec rectangles.  
\newblock {\em Comm. Math. Phys.}, 344:223--274, 2016.

\bibitem{ranglep} M. Ciucu.  
\newblock Lozenge tilings with gaps in a $90^\circ$ wedge domain with mixed boundary conditions.  
\newblock {\em Comm. Math. Phys.}, 334:507--532, 2015.

\bibitem{symffa} M. Ciucu.  
\newblock Symmetries of shamrocks, Part I.  
\newblock {\em J. Combin. Theory Ser. A}, 155:376--397, 2018.

\bibitem{DT} G. David and C. Tomei.  
\newblock The problem of the calissons.  
\newblock {\em Amer\. Math\. Monthly}, 96:429--431, 1989.

%\bibitem{} .  
%\newblock .  
%\newblock {\em }, 96:429--431, 1989.

\bibitem{KKZ} C. Koutschan, M. Kauers and D. Zeilberger.  
\newblock A proof of George Andrews' and David Robbins' $q$-TSPP-conjecture.  
\newblock {\em Proc. Natl. Acad. Sci. USA}, 108:2196--2199, 2011.

\bibitem{KuoOne} E. H. Kuo.  
\newblock Applications of graphical condensation for enumerating matchings and tilings.  
\newblock {\em Theoret. Comput. Sci.}, 319:29--57, 2004.

\bibitem{KuoTwo} E. H. Kuo.  
\newblock Graphical condensation generalizations involving Pfaffians and determinants.  
\newblock arxiv preprint math:CO/06055154, 2006.

\bibitem{Kup} G. Kuperberg.  
\newblock Symmetries of plane partitions and the permanent-de\-ter\-mi\-nant method.  
\newblock {\em J. Combin. Theory Ser. A}, 68:115--151, 1994.

\bibitem{MacM} P. A. MacMahon.  
\newblock Combinatory Analysis, vols. 1--,.  
\newblock Cambridge, 1916, reprinted by Chelsea, New York, 1960..

\bibitem{Proc} R. A. Proctor.  
\newblock Odd symplectic groups.  
\newblock {\em Invent. Math.}, 92:307--332, 1988.

\bibitem{Sta} R. P. Stanley.  
\newblock Symmetries of plane partitions.  
\newblock {\em J. Comb. Theory Ser. A}, 43:103--113, 1986.

\bibitem{Ste} J. R. Stembridge.  
\newblock The enumeration of totally symmetric plane partitions.  
\newblock {\em Adv. in Math.}, 111:227--243, 1995.




%\bibitem{KSU} Masashi Kiyomi, Toshiki Saitoh, and Ryuhei Uehara.
%  \newblock Reconstruction of interval graphs.  \newblock In {\em
%    Computing and combinatorics}, volume 5609 of {\em Lecture Notes in
%    Comput. Sci.}, pages 106--115. Springer, 2009.

%\bibitem{DT} G. David and C. Tomei.  
%\newblock The problem of the calissons.  
%\newblock {\em Amer\. Math\. Monthly}, 96:429--431, 1989.
%
%
%\bibitem{GV} I. M. Gessel and X. Viennot. 
%\newblock Binomial determinants, paths, and hook length formulae.
%\newblock {\em Adv. in Math.}, 58:300--321, 1985.
%
%\bibitem{Kr} G. Kreweras.
%\newblock  Sur une classe de probl\`emes d\'enombrement li\'es au trellis des partitions de entiers.
%\newblock {\em Cahiers du Bur. Univ. de Rech. Op\'er.}, 6:5--105, 1965.  
%
%\bibitem{Na} T. V. Narayana.
%\newblock   A combinatorial problem and its application to probability theory.
%\newblock {\em J. Indian Soc. Agr. Stat.}, 7:169--178, 1955.
%
%\bibitem{Natwo} T. V. Narayana.
%\newblock  Lattice path combinatorics with statistical applications.
%\newblock {\em Mathematical Expositions No. 23}, University of Toronto Press, Toronto, 1979. 


\end{thebibliography}

\end{document}


