<html><head><meta name="color-scheme" content="light dark"></head><body><pre style="word-wrap: break-word; white-space: pre-wrap;">\documentclass[12pt,reqno]{article}

\usepackage[usenames]{color}
\usepackage{amssymb}
\usepackage{graphicx}
\usepackage{amscd}
\usepackage{bm}

\usepackage[colorlinks=true,
linkcolor=webgreen,
filecolor=webbrown,
citecolor=webgreen]{hyperref}

\definecolor{webgreen}{rgb}{0,.5,0}
\definecolor{webbrown}{rgb}{.6,0,0}

\usepackage{color}
\usepackage{fullpage}
\usepackage{float}

\usepackage{psfig}
\usepackage{graphics,amsmath,amssymb}
\usepackage{amsthm}
\usepackage{amsfonts}
\usepackage{latexsym}
\usepackage{epsf}

\setlength{\textwidth}{6.5in}
\setlength{\oddsidemargin}{.1in}
\setlength{\evensidemargin}{.1in}
\setlength{\topmargin}{-.1in}
\setlength{\textheight}{8.4in}

\newcommand{\seqnum}[1]{\href{http://oeis.org/#1}{\underline{#1}}}

\begin{document}

\begin{center}
\epsfxsize=4in
\leavevmode\epsffile{logo129.eps}
\end{center}

\theoremstyle{plain}
\newtheorem{theorem}{Theorem}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}

\theoremstyle{definition}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{example}[theorem]{Example}
\newtheorem{conjecture}[theorem]{Conjecture}

\theoremstyle{remark}
\newtheorem{remark}[theorem]{Remark}

\newcommand{\pq}[2]{\genfrac{[}{]}{0pt}{}{#1}{#2}}

\begin{center}
\vskip 1cm{\Large\bf Two Approaches to Normal Order Coefficients}
\vskip 1cm
\large
Richell O. Celeste\\
Institute of Mathematics and Natural Sciences Research Institute\\
University of the Philippines --- Diliman\\
Quezon City, Philippines 1101\\
\href{mailto:ching@math.upd.edu.ph}{\tt ching@math.upd.edu.ph} \\
\ \\
Roberto B. Corcino\\
Mathematics Department\\
Cebu Normal University\\
Cebu City, Philippines 6000\\
\href{mailto:rcorcino@yahoo.com}{\tt rcorcino@yahoo.com} \\
\ \\
Ken Joffaniel M. Gonzales\\
Department of Physical Sciences and Mathematics\\
University of the Philippines --- Manila\\
Manila, Philippines 1004\\
\href{mailto:kmgonzales@upd.edu.ph}{\tt kmgonzales@upd.edu.ph} \\
\ \\

\end{center}

\vskip .2 in

\baselineskip14pt
\begin{abstract} 
We consider the normal ordering coefficients of strings consisting of
the symbols $V,U$ which satisfy the commutation rule $UV-qVU=hV^s$.
These coefficients are studied using two approaches. First, we continue
the study by Varvak, where the coefficients were interpreted as $q$-rook
numbers under the row creation rook model introduced by Goldman and
Haglund. Second, we express the coefficients in terms of a kind of
generalization of some symmetric functions. We derive identities
involving the coefficients including some explicit formulas.
\end{abstract}

\section{Introduction}
Let $V,U$ be operators (or variables) that satisfy the commutation rule $UV-qVU=hV^s$, where $s\in\mathbb N$ and $h,q\in\mathbb R$. For example, if $s=0, h=1, q=1$, then $V,U$, respectively, can be represented by the creation operator and annihilation operator in quantum physics \cite{Bla}, or by the operators $X,\partial_x$ whose action on a monomial $x^n$ are given by $Xx^n=x^{n+1}$ and $\partial_xx^n=nx^{n-1}$. Given a string $w$ consisting of $V$'s and $U$'s, the \emph{normally ordered} form of $w$ is an equivalent operator expressed as a sum $\sum c_{i,j} V^i U^j$. The normally ordered form can be computed using the commutation rule alone, i.e., by replacing all occurrences of $UV$ with $qVU+hV^s$, but this task can be cumbersome especially for long strings. It turns out that the coefficients $c_{i,j}$, called \emph{normal ordering coefficients}, can be computed more efficiently using combinatorial techniques. In the classical case  $s=0, h=1, q=1$, Navon \cite{Nav} showed that the normal ordering coefficients of an arbitrary string are given by rook numbers on a Ferrers board. Varvak \cite{Var} generalized Navon's result for arbitrary $q$ and derived explicit formulas for these coefficients using rook factorization. Blasiak \cite{Bla}, El-Desouky 
et al.~\cite{DesCakMan}, and Mansour et al.~\cite{ManSchSha1,ManSchSha2} also computed explicit formulas using other methods.

In this paper, we study normal ordering coefficients using two approaches. In Section \ref{first}, we study the coefficients as $q$-rook numbers under the row creation rule introduced by Goldman and Haglund \cite{GolHag}. In Section \ref{second}, we study the coefficients by expressing them in terms of some generalization of elementary and complete homogeneous symmetric functions. Lastly, some special cases are given in Section \ref{special}.

\section{First approach: rook numbers}\label{first}
Let $\textbf v=(v_1,v_2,\ldots,v_n), \textbf u=(u_1,u_2,\ldots,u_n)$, and $H_{{\textbf{v}},{\textbf{u}}}=V^{v_n}U^{u_n}\cdots V^{v_2}U^{u_2}V^{v_1}U^{u_1}$. In this section, we obtain explicit formulas for the normal ordering coefficients of $H_{{\textbf{v}},{\textbf{u}}}$ which uses the known rook theoretic interpretation directly. We also give a representation of $V,U$ in terms of linear operators and use it to find another explicit formula which generalizes Varvak's \cite[Corollary 4.2]{Var}.

Following Blasiak \cite{Bla}, we write the string $H_{{\textbf{v}},{\textbf{u}}}$ in the form
\begin{equation}\label {wnorm}
H_{{\textbf{v}},{\textbf{u}}} = \sum_
{k=u_1}^{|\textbf u|} S^{\textbf v, \textbf u}_{s,h;q}[k]\, V^{|\textbf v| - (|\textbf u|-k)(1-s)} U^k\,,
\end{equation}
where $|\textbf u|=u_1+u_2+\cdots+u_n$ and $|\textbf v|=v_1+v_2+\cdots+v_n$.

Varvak \cite{Var} showed that for $h=1,q=1$, the coefficients $S^{\textbf v,\textbf u}_{s,h;q}[k]$ also occur as rook numbers under the rook model introduced by Goldman and Haglund \cite{GolHag} which we now describe. An $s$-\emph{rook placement} on a Ferrers board $B$ is obtained as follows. First, choose the columns where rooks will be placed. The rooks are then placed one by one from the right such that every time a rook is placed in a cell the entire row to its left is divided into $s$ rows. When $s=0$, ``division'' into $s$ rows can be interpreted as cancellation of the entire row lying to the left of a rook. Denote by $\mathcal R_s(B,k)$ the set of all placements of $k$ rooks on $B$. An example of a rook placement where $s=2$ is shown in Figure \ref{rooktab}. The $k$-th $s$-rook number of a board $B$ is then defined as $R_{s}(B(w),k)=|\mathcal R_s(B,k)|$.

\begin{figure}[htbp]
\begin{center}
\epsfxsize=1.5in
\leavevmode\epsffile{imw.eps}
\caption{A placement of 3 rooks.}
\label{rooktab}
\end{center}
\end{figure}

Varvak \cite[Theorem 7.1]{Var} showed, in our notation (and after correcting for the typo pointed out by Mansour et al.~\cite{ManSchSha1}), that
\begin{equation}\label{rooknorm}
H_{{\textbf{v}},{\textbf{u}}} = \sum_ {k=0}^{|\textbf u|-u_1} R_{s,1;1}[B(H_{{\textbf{v}},{\textbf{u}}}),k]\, V^{|\textbf v|-k(1-s)} U^{|\textbf u|-k}.
\end{equation}
Comparing \eqref{rooknorm} with \eqref{wnorm} gives $S^{{\textbf{v}}, {\textbf{u}}}_{s,1;1}[k] =
R_{s,1;1}[B(H_{{\textbf{v}},{\textbf{u}}}),|\textbf u|-k]$.

The connection between rook numbers and normal ordering lies on the fact that a string $w$ determines a unique Ferrers board (where a column or row is allowed to have length zero) and that the placement or non-placement of a rook corresponds to the choice of replacement for each occurrence of $UV$. Specifically, a string $w$ outlines a Ferrers board $B(w)$ whose underside border is the lattice path obtained from $w$ by replacing $V$ with a unit step up and $U$ with a unit step to the right. For example, the Ferrers board in Figure \ref{rooktab} is outlined by $VUVUUVVUVVU$. Placing a rook on the northeast inner corner of $B(w)$ corresponds to replacing the rightmost $UV$ with $hV^s$ while leaving a cell empty corresponds to replacing $UV$ with $qVU$. Under this correspondence, it is now clear how Varvak's result can be extended to arbitrary $q$. Given a rook placement $\phi$, assign the weight $w(\phi)=h^{t'}q^t$ where $t'$ is the number of rooks, and $t$ is the number of cells not containing a rook and not lying above a rook if $s\neq 0$ and $t$ is the number of cells not containing a rook and not lying above or to the left of a rook if $s=0$. For example, the rook placement in Figure \ref{rooktab} has weight $h^3q^{11}$. We define the generalized rook numbers by
\[
R_{s,h;q}[B,k] = \sum_{\phi\in \mathcal R_s(B,k)} w(\phi).
\]
When all three parameters $s,h,q$ are arbitrary, the normally ordered form of the string $H_{\textbf v,\textbf u}$ is then given by
\begin{equation}
H_{{\textbf{v}},{\textbf{u}}} = \sum_ {k=0}^{|\textbf u|-u_1} R_{s,h;q}[B(H_{{\textbf{v}},{\textbf{u}}}),k]\, V^{|\textbf v|-k(1-s)} U^{|\textbf u|-k}.\label{neweq}
\end{equation}
Again, comparing \eqref{neweq} with \eqref{wnorm} gives $S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k] =
R_{s,h;q}[B(H_{{\textbf{v}},{\textbf{u}}}),|\textbf u|-k]$.

Our first result gives the normally ordered form of strings of the form $U^{m}V^{l}$. The following notation will be used. For $n\in\mathbb N$, let $[n]_q=1+q+q^2+\cdots+q^{n-1}$ if $n&gt;0$ and $[0]_q=0$. Let $[n]_q!=[n]_q[n-1]_q\cdots[2]_q[1]_q$ for $n&gt;0$ and $[0]_q!=1$. The $q$-binomial coefficient is defined by $\pq{n}{k}_q=\frac{[n]_q!}{[k]_q![n-k]_q!}$. They satisfy the property (see [6,\,Table 1 and Identity 2.2]) 
\begin{align}
\pq{n}{k}_q = \sum_{0 \leq i_1 \leq i_2 \leq \ldots \leq i_{n-k} \leq k} q^{i_1+i_2+\cdots+i_{n-k}}\,. \label{qbinomp}
\end{align}

\begin{theorem}\label{lemrep} For $s\in\mathbb N$, the normally ordered form of $U^{m}V^{l}$ is given by
\begin{equation}
U^{m}V^{l} = \sum_{j=0}^{m} \left( h^jq^{l(m-j)} \pq{m}{j}_{q^{s-1}} \prod_{i=0}^{j-1} [l+i(s-1)]_q \right)V^{l+j(s-1)} U^{m-j}.\label{repx}
\end{equation}
\begin{proof} The string $U^mV^l$ outlines a rectangular Ferrers board with $m$ columns and $l$ rows. The total weight of all placements of $j$ rooks on this board equals the coefficient of $V^{l+j(s-1)}U^{m-j}$ in the normally ordered form of $U^mV^l$. We now compute the total weight of such rook placements as follows. Choose $j$ columns where rooks will be placed. If the first rook is placed on the cell in the $i$th row, $1\leq i\leq l$, then the cells below the rook will contribute a weight of $q^{i-1}$. As $i$ varies, a total weight of $h(1+q+\cdots+q^{l-1})=h[l]_q$ will be contributed by all possible placements of the first rook. Since the placement of the first rook adds $s-1$ subcells to every cell to its left, the total weight contributed by all possible placement of the second rook is $[l+(s-1)]_q$. Continuing this process with the other columns, we see that the weight contributed by all possible placements of $j$ rooks in the chosen columns is $h^j\prod_{i=0}^{j-1} [l+i(s-1)]_q$, and that this weight is the same for any choice of $j$ columns.

We now consider the weight contributed by the other columns in which no rooks are placed. For such a column, the weight is completely determined by the number of columns to its right that contains a rook, i.e., if there are $t$ columns to its right containing a rook, then the column will assume a weight of $q^{l+t(s-1)}$. Note that $t$ varies from $0$ to $j$ and that for a given placement of $j$ rooks, the weight contributed by all the columns containing no rooks is $q^{lt_0} q^{(l+(s-1))t_1} q^{(l+2(s-1))t_2} \cdots q^{(l+j(s-1))t_j}$ for some $t_0+t_1+\cdots+t_j=m-j$. Summing this up over all such possible collections $\{t_0,t_1,\ldots,t_i\}$, we have
\begin{align*}
\sum_{t_0+t_1+\cdots+t_j=m-j} &amp;q^{(l+0(s-1))t_0}  q^{(l+1(s-1))t_1} q^{(l+2(s-1))t_2} \cdots q^{(l+j(s-1))t_j} \\
&amp;= q^{l(m-j)} \sum_{t_0+t_1+\cdots+t_j=m-j} q^{0(s-1)t_0} q^{1(s-1)t_1} q^{2(s-1)t_2} \cdots q^{j(s-1)t_j}\\
&amp;= q^{l(m-j)} \sum_{0 \leq i_1 \leq i_2 \leq \ldots \leq i_{m-j} \leq j} q^{(s-1)(i_1+i_2+\cdots+i_{m-j})}\\
&amp;= q^{l(m-j)} \pq{m}{j}_{q^{s-1}}\,,
\end{align*}
where the last equality follows from \eqref{qbinomp}. This proves the theorem. \end{proof}\end{theorem}

\begin{corollary}\label{theoremrep} Let $s\in\mathbb N$. The string $H_{{\textbf{v}},{\textbf{u}}}$ may be written as \begin{multline}\label{longsum}
H_{{\textbf{v}},{\textbf{u}}} = \sum_{j_1=0}^{u_2} \sum_{j_2=0}^{u_3} \cdots \sum_{j_{n-1}=0}^{u_n} \prod_{i=1}^{n-1}
h^{j_1+\cdots+j_{n-1}} \Gamma_{q,s} [j_i,v_1+\cdots+v_i+(j_1+\cdots+j_{i-1})(s-1),u_{i+1}]_q \\ V^{v_1+\cdots+v_n + (j_1+\cdots+j_{n-1})(s-1)}
U^{u_1+\cdots+u_n-(j_1+\cdots+j_{n-1})}.
\end{multline}
where
\begin{equation*} \Gamma_{q,s}[j,l,m] = q^{l(m-j)} \pq{m}{j}_{q^{s-1}}
\prod_{i=0}^{j-1} [l+i(s-1)]_q.
\end{equation*}
Hence, the numbers $S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k]$ are given by
\begin{equation}\label{exp22}
S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k] =
h^{|{\textbf{u}}|-k}\sum_{j_1+\cdots+j_{n-1}=u_1+\cdots+u_n-k}~\prod_{i=1}^{n-1 } \Gamma_{q,s}
[j_i,v_1+\cdots+v_i+(j_1+\cdots+j_{i-1})(s-1),u_{i+1}].
\end{equation}
\begin{proof} Identity \eqref{longsum} is proved by repeated application of \eqref{repx} beginning from $U^{u_2}V^{v_1}$. Identity \eqref{exp22} follows by comparing the coefficient of $U^k$ in \eqref{longsum} and \eqref{wnorm}. \end{proof} \end{corollary}

We note that the case $s=0, h=1, q=1$ of Corollary \ref{theoremrep} was derived by El Desouky et al.~\cite{DesCakMan} using the Leibniz formula.

\begin{corollary}\label{correp} Let $s\in\mathbb N$ and set $\textbf v=\textbf u=(\underbrace{1,1,\ldots,1}_{n})$. Then the following explicit formula for $S^{\textbf v,\textbf u}_{s,h;q}[k]$ holds
\begin{equation*}
S^{\textbf v,\textbf u}_{s,h;q}[k] = h^{n-k} \sum_{j_1+\cdots+j_{n-1}=n-k}~\prod_{i=1}^{n-1} q^{(i+(j_1+\cdots+j_{i-1})(s-1))(1-j_i)} \pq{i+(j_1+\cdots+j_{i-1})(s-1)}{j_i}_q.
\end{equation*}
\end{corollary}

Varvak's \cite{Var} use of rook factorization to obtain an explicit formula adapts readily in the case of $S^{{\textbf{v}},{\textbf{u}}}_{s,h;q}[k]$ after some modification. We will need the following analogues of the falling factorial and factorial: for $r\in\mathbb R,j\in\mathbb N$, define $ [r]^{(\underline{j})}_{q,1-s} = [r(1-s)]_q [(r-1)(1-s)]_q \cdots [(r-j+1)(1-s)]_q$ and for $n\in\mathbb N$, define $[n]_{q,1-s}!=[n]^{(\underline{n})}_{q,1-s}$.

\begin{theorem}\label{expexp}
Let $s\neq 1$. The coefficients $S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k]$ satisfy the explicit formula
\begin{equation*}
S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k] = \frac{h^{|{\textbf{u}}|-k}}{[k]_{q,1-s}!} \sum_{j=0}^k (-1)^{k-j} q^{\binom{k-j}{2}(1-s)} \pq{k}{j}_{q^{1-s}} \Omega^{{\textbf v},{\textbf u}}_{s;q}[j],
\end{equation*}
where
\begin{equation*}
\Omega^{{\textbf v},{\textbf u}}_{s;q}[j]=\prod_{t=1}^{n} [j-(u_1+u_2+\cdots+u_{t-1})+(v_1+v_2+\cdots+v_{t-1})/(1-s)]_{q,1-s}^{(\underline{u_t})}.
\end{equation*}
\begin{proof} We use a representation of $V,U$ as linear operators whose action on the monomial $t^j$ is given by $Vt^j=t^{j+1}$ and $Ut^j=h[j]_qt^{j+s-1}$. One can verify that these operators satisfy $VU-qVU=hV^s$ and that $U^k t^{n(1-s)} =h^k [n]_{q,1-s}^{(\underline{k})} t^{(n-k)(1-s)}$. We then apply both sides of \eqref{wnorm} to $t^{x(1-s)}$. After letting $t=1$ to the resulting equation and using the property $[x]^{(\underline{k})}_{q,1-s}=[1-s]_q^k [x]_{q^{1-s}}[x-1]_{q^{1-s}} \cdots[x-k+1]_{q^{1-s}}$, we obtain
\begin{equation*}
h^{|\textbf u|} \Omega^{{\textbf v},{\textbf u}}_{s;q}[x]= \sum_{k=u_1}^{|\textbf u|} h^{k} S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k][1-s]_q^k [x]_{q^{1-s}}[x-1]_{q^{1-s}} \cdots[x-k+1]_{q^{1-s}}.
\end{equation*}

Let $E$ denote the shift operator $EP(x)=P(x+1)$ and $\Delta_{Q}^k$ the $k$-th $Q$-difference operator defined by $\Delta_{Q}^k=(E-1)(E-Q)\cdots (E-Q^{k-1})$. If $P(x)=\sum_{k}p_k [x]_Q [x-1]_Q \cdots [x-k+1]_Q$, then $p_k=\frac{1}{[k]_Q!} \Delta_Q^k P(x) \vert_{x=0}$. By the $q$-binomial theorem, $\Delta_{Q}^k = \sum_{j=0}^k (-1)^{j} Q^{\binom{j}{2}} \pq{k}{j}_{Q} E^{k-j}$. The result then follows by letting $Q=q^{1-s}$, $p_k=h^{k} S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k][1-s]_q^k$ and $P(x)=h^{|\textbf u|} \Omega^{{\textbf v},{\textbf u}}_{s;q}[x]$.
 \end{proof} \end{theorem}

\begin{corollary}\label{secexp} Let $s\in\mathbb N \backslash \{0,1\}$. If $\textbf v=\textbf u=(\underbrace{1,1,\ldots,1}_{n})$, then the numbers $S^{\textbf v,\textbf u}_{s,h;q}[k]$ have the following explicit formula
\begin{equation*}
S^{\textbf v,\textbf u}_{s,h;q}[k] = \frac{h^{n-k}[s]_q^n}{[k]_{q^{1-s}}![1-s]_q^k} \sum_{j=0}^{k} (-1)^{k-j} q^{\binom{k-j}{2}(1-s)}
\pq{k}{j}_{q^{1-s}} \prod_{t=1}^n [(j/s)+t-j-1]_{q^s}.
\end{equation*}
When $s=0$,
\begin{equation*}
S^{\textbf v,\textbf u}_{0,h;q}[k] = \frac{h^{n-k}}{[k]_q!} \sum_{j=0}^{k} (-1)^{k-j} q^{\binom{k-j}{2}} \pq{k}{j}_{q} [j]_q^n.
\end{equation*} \end{corollary}

Theorem \ref{expexp} can be used to derive a Dobinsky formula for the Bell numbers corresponding to $S^{\textbf v,\textbf u}_{s,h;q}[k]$. Define the generalized $q$-Bell polynomials $B^{\textbf v,\textbf u}_{s,h;q}[x]$ and generalized $q$-Bell numbers $B^{\textbf v,\textbf u}_{s,h;q}$ by
\begin{equation*}
B^{\textbf v,\textbf u}_{s,h;q}[x] = \sum_{k=u_1}^{|\textbf u|} S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k] x^{k}, \hspace{0.25in} B^{\textbf v,\textbf u}_{s,h;q}=B^{\textbf v,\textbf u}_{s,h;q}[1].
\end{equation*}
When $s=0,h=1,q=1$  and the associated string $H_{\textbf v,\textbf u}$ outlines a staircase board, i.e., when $\textbf v=\textbf u=(\underbrace{1,1,\ldots,1}_{n})$, the numbers $B^{\textbf v,\textbf u}_{s,h;q}[x]$ and $B^{\textbf v,\textbf u}_{s,h;q}$ reduce to the usual Bell polynomial $B(n;x)$ and Bell number $B(n)$, respectively. The classical Dobinsky formula is given by 
\begin{equation*}
B(n;x) = \frac{1}{e^x} \sum_{j=0}^{\infty} j^n \frac{x^j}{j!}.
\end{equation*}
The Dobinsky formula corresponding to $B^{\textbf v,\textbf u}_{s,h;q}[x]$ is as follows.

\begin{corollary}\label{dobcor} Let $s\in\mathbb N \backslash \{1\}$ and $\Omega^{{\textbf v},{\textbf u}}_{s;q}[j]$ be as in Theorem \ref{expexp}. Then,
\begin{equation*}
B^{\textbf v,\textbf u}_{s,h;q}[x]=\left(\sum_{j=0}^{\infty} h^{|\textbf u|-j} (-1)^j q^{\binom{j}{2}(1-s)}
\frac{x^{j}}{[j]_{q,1-s}!} \right)\left(\sum_{j=0}^{\infty} \Omega^{{\textbf v},{\textbf u}}_{s;q}[j]
\frac{x^{j}}{h^j[j]_{q,1-s}!}\right).
\end{equation*}
\begin{proof} Using Theorem \ref{expexp}, the property $S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k]=0$ when $|\textbf u|&lt;k&lt;u_1$ and the property $\pq{k}{j}_{q^{1-s}}=\frac{[k]_{q,1-s}!}{[j]_{q,1-s}![k-j]_{q,1-s}!}$, respectively, we have
\begin{align*}
B^{\textbf v,\textbf u}_{s,h;q}[hx] &amp;= \sum_{k=u_1}^{|\textbf u|} S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k] h^k x^{k}\\
&amp;= \sum_{k=0}^{\infty} \sum_{j=0}^k \frac{h^{|{\textbf{u}}|}}{[k]_{q,1-s}!} (-1)^{k-j} q^{\binom{k-j}{2}(1-s)} \pq{k}{j}_{q^{1-s}} \Omega^{{\textbf v},{\textbf u}}_{s;q}[j] x^k\\
&amp;= \sum_{k=0}^{\infty} \sum_{j=0}^k \left(\frac{h^{|{\textbf{u}}|}
(-1)^{k-j} q^{\binom{k-j}{2}(1-s)}}{[k-j]_{q,1-s}!}\right)\left( \frac{\Omega^{{\textbf v},{\textbf
u}}_{s;q}[j]}{[j]_{q,1-s}!}\right) x^k.
\end{align*}
The desired result is then obtained after replacing $x$ with $x/h$ and by using the Cauchy product rule.
\end{proof} \end{corollary}

\section{Second approach: generalized symmetric functions}\label{second}
In this section, we derive formulas for $S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k]$ by showing that it can be expressed in terms of a certain generalization of some symmetric functions. Our derivation will involve the following shift in notation. Let $w=H_{{\textbf{v}}, {\textbf{u}}}$. Label the columns of $B(w)$ from the right by $0,1,2,\ldots,|\textbf u|-1$ and let $w_i$ be the length of the column labeled $i$. Let $B_0(w)=\varnothing$ and for $n=1,2,\ldots,|\textbf u| -1 $, let $B_n(w)$ be the Ferrers board whose column lengths are $w_0,w_1,w_2\ldots,w_{n-1}$. Finally, define $S^w_{s,h;q}[n,k]=R_{s,h;q}[B_n(w),n-k]$.

\begin{theorem}Let $n=0,1,2,\ldots,|\textbf u|-1,k\leq n$. The rook numbers $R_{s,h;q}[B_n(w),k]$ satisfy the recurrence relation
\begin{multline}
R_{s,h;q}[B_n(w),k]= q^{w_{n-1}+k(s-1)}R_{s,h;q}[B_{n-1}(w),k] \\ + h[w_{n-1}+(k-1)(s-1)]_q\,R_{s,h;q}[B_{n-1}(w),k-1]\label{rookrec}
\end{multline}
with boundary conditions $R_{s,h;q}[B_n(w),n]=\prod_{i=0}^{n-1} [w_i+(i-1)(s-1)]_q$ and $R_{s,h;q}[B_0(w),0]=R_{s,h;q}[B_n(w),0]=1$.
\begin{proof}
The identity can be proved by induction but a combinatorial proof appears to be more straightforward. Divide the set $\mathcal R(B(w),k)$ into two, with the first set consisting of rook placements without a rook in the $m$-th column and the other set consisting of rook placements with a rook in $m$-th column. The sum of the weights of all rook placements in the first set is $q^{w_{n-1}+k(s-1)}R_{s,h;q}[B_{n-1}(w),k]$ while that of the second set is $h[w_{n-1}+(k-1)(s-1)]_qR_{s,h;q}[B_{n-1}(w),k-1]$.\label{recrook}
\end{proof}
\end{theorem}

From Equation \eqref{rookrec}, we have
\begin{multline}
S^w_{s,h;q}[n,k] = q^{w_{n-1}+(n-k)(s-1)} S^w_{s,h;q}[n-1,k-1] \\ + h[w_{n-1}+(n-k-1)(s-1)]_q S^w_{s,h;q}[n-1,k]. \label{strec}
\end{multline}
Now, multiply both sides of \eqref{strec} by $q^{(s-1)\binom{k+1}{2}-(s-1)\binom{n+1}{2}-\sum_{i=0}^{n-1} w_i}$ and let
\[
s_{n,k}=q^{(s-1)\binom{k+1}{2}-(s-1)\binom{n+1}{2}-\sum_{i=0}^{n-1} w_i}S^w_{s,h;q}[n,k].
\]
Equation \eqref{strec} then becomes
\[
s_{n,k}=s_{n-1,k-1}+ hq^s \left([w_{n-1}+(n-1)(s-1)]_{1/q}-[k(s-1)]_{1/q}\right) s_{n-1,k}.
\]
This shows that the numbers $s_{n,k}$, and hence $S^w_{s,h;q}[n,k]$, can be studied under a more general setting. (A similar observation was made by Mansour
et al.~\cite{ManSchSha2} in the case of the normal ordering coefficients of $(VU)^n$.) Specifically, let $\bm{\alpha}=(\alpha_0,\alpha_1,\ldots)$ and $\textbf w=(\beta_0,\beta_1,\ldots)$ be sequences (or equivalently, \emph{weight functions}) and define the numbers
$A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k}$ by the recurrence relation
\begin{equation}\label{ankrec}
A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k} = A^{{\bm{\alpha}}, {\bm{\beta}}}_{n-1,k-1} + (\alpha_{n-1} + \beta_{k}) A^{{\bm{\alpha}}, {\bm{\beta}}}_{n-1,k},
\end{equation}
with initial conditions $A^{{\bm{\alpha}}, {\bm{\beta}}}_{0,n}=\delta_{0,n}$ and $A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,0} = (\alpha_{n-1} + \beta_0) (\alpha_{n-2} + \beta_0) \cdots (\alpha_{0} + \beta_0)$. When only the value of a weight function at $i=0,1,2,\ldots$ is specified (for instance, $\alpha_i$), it is understood that the corresponding weight function is the same letter in boldface without the subscript.

In the theorem that follows, we denote the matrix whose $(n,k)$ entry is $M_{n,k}$ by $\left[M_{n,k}\right]$. The zero weight function is denoted $\textbf{0}=(0,0,\ldots)$. The matrices are all assumed to be square.

\begin{theorem} There holds the matrix factorization $\left[ A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k}\right] = \left[
A^{{\bm{\alpha}}, {\textbf{0}}}_{n,k} \right] \left[ A^{{\textbf {0}}, {\bm{\beta}}}_{n,k}\right]$. Equivalently, $A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k} = \sum_{j=k}^{n} A^{{\bm{\alpha}},{\textbf {0}}}_{n,j}A^{{\textbf {0}},{\bm{\beta}}}_{j,k}$.
\begin{proof} By \eqref{ankrec},
\begin{align*}
A^{\bm{\alpha}, \textbf{0}}_{n,j}A^{\textbf {0}, \bm{\beta}}_{j,k} &amp;= (A^{\bm{\alpha}, \textbf{0}}_{n-1,j-1} + \alpha_{n-1} A^{\bm{\alpha}, \textbf{0}}_{n-1,j}) (A^{\textbf {0}, \bm{\beta}}_{j-1,k-1} + \beta_{k} A^{\textbf {0}, \bm{\beta}}_{j-1,k}) \\
&amp;= A^{\bm{\alpha},\textbf{0}}_{n-1,j-1} A^{\textbf {0}, \bm{\beta}}_{j-1,k-1} + \alpha_{n-1} A^{\bm{\alpha}, \textbf{0}}_{n-1,j} A^{\textbf {0}, \bm{\beta}}_{j,k} + \beta_k A^{\bm{\alpha}, \textbf{0}}_{n-1,j-1} A^{\textbf {0}, \bm{\beta}}_{j-1,k}
\end{align*}
Hence, the numbers $A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k}$ viewed as the $(n,k)$ entry of the matrix product $\left[
A^{{\bm{\alpha}}, {\textbf{0}}}_{n,k} \right] \left[ A^{{\textbf {0}}, {\bm{\beta}}}_{n,k}\right]$ satisfy the same recursion as the numbers defined by Equation \eqref{ankrec}.
\end{proof}
\label{fact}
\end{theorem}

The theorem that follows gives different formulations for the numbers $A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k}$ in terms of expressions that are analogous to elementary and complete symmetric functions.

\begin{theorem}\label{symfcn} The following identities hold: 
\begin{align}
A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k} &amp;=\sum_{0\leq i_1 &lt; i_2 &lt; \ldots &lt; i_{n-k} \leq n-1}~\prod_{j=1}^{n-k} (\alpha_{i_j} + \beta_{i_j-j+1}) \label{tdset}\\
A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k} &amp;= \sum_{0\leq i_1 \leq i_2 \leq \ldots \leq i_{n-k} \leq k}~\prod_{j=1}^{n-k} (\alpha_{i_j+j-1} + \beta_{i_j}) \label{tset}\\
A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k} &amp;=\sum_{i_0+ i_1 + i_2 + \cdots + i_k = n-k} ~\prod_{j=0}^{k}~\prod_{l=0}^{i_j-1} (\alpha_{j + l + i_0 + i_1+i_2+ \cdots +i_{j-1}} + \beta_j) .\label{tsetalt}
\end{align}
\begin{proof} It can be shown by partitioning the set of indices that the $RHS$ of \eqref{tdset} and \eqref{tset} satisfy the recursion \eqref{ankrec}. Identity \eqref{tsetalt} is a restatement of \eqref{tset}. \end{proof} \end{theorem}
 
Let $H$ be a set. The $m$-th \emph{elementary symmetric function} $e$ and $m$-th \emph{complete symmetric function} $h$ are defined as follows: $e_{m}(H)$ (resp. $h_{m}(H)$) is the sum of all products of $m$ elements from $H$ taken without (resp. with) replacement. Observe that by \eqref{tdset}, $A^{\bm{\alpha}, \textbf{0}}_{n,k}=e_{n-k}(\{\alpha_0,\alpha_1,\ldots,\alpha_{n-1}\})$ and by \eqref{tset} or \eqref{tsetalt}, $A^{\textbf {0}, \bm{\alpha}}_{n,k}=h_{n-k}(\{\alpha_0,\alpha_1,\ldots,\alpha_{k}\})$. This shows that the numbers $A^{\bm{\alpha}, \textbf{0}}_{n,k}$ and $A^{\textbf {0}, \bm{\alpha}}_{n,k}$ are exactly the $\mathcal{U}$-Stirling numbers studied by Medicis and Leroux \cite{MedLer2} which they defined in terms of $\mathcal{A}$-tableaux. 

We now give an interpretation for $A^{\bm{\alpha}, \bm{\beta}}_{n,k}$ in terms of certain lattice paths which reduces to the $\mathcal{A}$-tableau interpretation for $A^{\bm{\alpha}, \textbf{0}}_{n,k}$ and $A^{\textbf {0}, \bm{\beta}}_{n,k}$. Let $m,m'\in\mathbb N$ and denote by $\mathcal{L}(m,m')$ the set of pairs of lattice paths $(\lambda,\lambda')$ on an $m$ by $m'$ grid consisting of unit left ($L$) and unit up ($U$) steps beginning from the origin (the lower left corner of the grid) such that $\lambda$ contains no two consecutive $L$'s and such that $\lambda'$ is the lattice path obtained from $\lambda$ by replacing every occurrence of $UL$ by $L$. Figure \ref{atab} shows an example of a pair of lattice path in $\mathcal{L}(6,5)$ where $\lambda$ is the lattice path in bold lines and $\lambda'$ is the lattice path in dotted lines.

\begin{figure}[htbp]
\begin{center}
\epsfxsize=1.5in
\leavevmode\epsffile{imy.eps}
\caption{A pair of lattice paths in $\mathcal{L}(6,5)$}
\label{atab}
\end{center}
\end{figure}

Given a weight function $\bm{\gamma}$, we assign to a lattice path $\lambda^*$ with $m^*$ $L$'s the weight $\gamma(\lambda^*)=\prod_{i=1}^{m^*} \gamma_{\lambda^*_i}$ where the $\lambda^*_i$'s are the second coordinates of the $L$'s of $\lambda^*$. For instance, the lattice path in bold lines in Figure \ref{atab} has weight $\gamma_0 \gamma_1 \gamma_2 \gamma_4 \gamma_6$. The next theorem gives the interpretation of $A^{\bm{\alpha}, \bm{\beta}}_{n,k}$ in terms of lattice paths and follows directly from Equation \eqref{tdset} of Theorem \ref{symfcn}.

\begin{theorem} The following interpretation for $A^{\bm{\alpha}, \bm{\beta}}_{n,k}$ holds
\[
A^{\bm{\alpha}, \bm{\beta}}_{n,k} = \sum_{(\lambda,\lambda')\in \mathcal{L}(n-1,n-k)} \left( \bm{\alpha}(\lambda) + \bm{\beta} (\lambda') \right).
\]\label{inter}
\end{theorem}

Going back to $S^w_{s,h;q}[n,k]$, the discussion in the beginning of this section shows that if $\alpha_i=hq^s [w_i+i(s-1)]_{1/q},\beta_i=-hq^s[i(s-1)]_{1/q}$, then
\begin{equation}
q^{(s-1)\binom{k+1}{2}-(s-1)\binom{n+1}{2}-\sum_{i=0}^{n-1} w_i} S^w_{s,h;q}[n,k]= A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k}. \label{conv}
\end{equation}
If $w=H_{{\textbf{v}}, {\textbf{u}}}$ where $B(w)$ has $n=|\textbf{u}|$ columns, then $S^w_{s,h;q}[n,k]=S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k]$. This shows that properties of $S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k]$ can be easily lifted from the properties of $A^{{\bm{\alpha}}, {\bm{\beta}}}_{n,k}$.

\begin{theorem} Let $n=|\textbf{u}|$. Then the numbers $S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k]$ satisfy the formula
\[
S^{{\textbf{v}}, {\textbf{u}}}_{s,h;q}[k] = f(n,k) h^{n-k} \sum_{i,j,l} (-1)^{j-l} A^{{\bm{\rho}}, {\textbf{0}}}_{n,j} \binom{j}{i} \binom{i}{l} \pq{l}{k}_{q^{1-s}},
\]
where $f_q(n,k)=q^{2(s+1)(n-k)-(s-1)\binom{k+1}{2}+(s-1)\binom{n+1}{2}-\sum_{i=0}^{n-1} w_i} (1-q)^{-2(n-k)}$, $\bm{\rho}$ is the weight function given by $\rho_i=(1/q)^{w_i+i(s-1)}$ and $w_i$ is the length of the $(i+1)$-st column from the rights of the Ferrers board $B(H_{\textbf{v},\textbf{w}})$.
Hence,
\[
H_{{\textbf{v}}, {\textbf{u}}}=\sum_{k=u_1}^{n}  \left( f_q(n,k) \sum_{i,j,l} h^{n-k} (-1)^{j-l} A^{{\bm\rho}, {\textbf{0}}}_{n,j} \binom{j}{i} \binom{i}{l} \pq{l}{k}_{q^{1-s}} \right) V^{|\textbf v| - (|\textbf u|-k)(1-s)} U^k.
\]

\begin{proof} For simplicity, let $h=1$. We can easily recover the more general identity using the fact that $S^w_{s,h;q}[n,k]=h^{n-k}S^w_{s,1;q}[n,k]$. Define the weight functions $\bm\alpha$ and $\bm\beta$ by $\alpha_i=q^s [w_i+i(s-1)]_{1/q},\beta_i=-q^s[i(s-1)]_{1/q}$. We first gather all the facts we will use. Also, let $Q=1/q$  by $\alpha^*_i=Q^{w_i+i(s-1)}$ and $\beta^*_i=Q^{i(s-1)}$. Then we can rewrite the weight functions as $\alpha_i=Q^{-s} \frac{\alpha^*_i-1}{Q-1}$, $\beta_i=-Q^{-s}\frac{\alpha^*_i-1}{Q-1}$.

For an arbitrary weight function $\bm{\gamma}$ and constant $c$, Equations \eqref{tdset} and \eqref{tset} imply that $A^{{c\bm\gamma}, {\textbf{0}}}_{n,k}=c^{n-k} A^{{\bm\gamma}, {\textbf{0}}}_{n,k}$ and $A^{{\textbf {0}}, {c\bm\gamma}}_{n,k}=c^{n-k} A^{{\textbf {0}}, {\bm\gamma}}_{n,k}$, respectively. Furthermore, if $\gamma_i=q^i$, then $A^{{\textbf {0}}, {\bm\gamma}}_{n,k}= \pq{n}{k}_q$ by the recurrence \eqref{ankrec} and the fact that $\pq{n}{k}_q=\pq{n-1}{k-1}_q+q^k \pq{n-1}{k}_q$. Using these facts and by the repeated use of Theorem \ref{fact},
\begin{align*}
A^{{\bm\alpha}, {\bm\beta}}_{n,k} &amp;= (Q^{s}(Q-1))^{-2(n-k)} \sum_i A^{{\bm\alpha}, {\textbf{0}}}_{n,i} A^{{\textbf {0}}, {\bm\beta}}_{i,k}\\
&amp;= (Q^{s}(Q-1))^{-2(n-k)} \sum_i A^{{\bm\alpha}^*, -{\textbf{1}}}_{n,i} A^{- {\textbf{1}}, {\bm\beta}^*}_{i,k}\\
&amp;= (Q^{s}(Q-1))^{-2(n-k)} \sum_{i,j,l} A^{{\bm\alpha}^*, {\textbf{0}}}_{n,j} A^{{\textbf {0}}, -{\textbf{1}}}_{j,i} A^{-{\textbf {1}}, {\textbf{0}}}_{i,l} A^{{\textbf {0}}, {\bm\beta}^*}_{l,k}\\
&amp;= (Q^{s}(Q-1))^{-2(n-k)} \sum_{i,j,l} (-1)^{j-l} A^{{\bm\alpha}^*, {\textbf{0}}}_{n,j} \binom{j}{i} \binom{i}{l} \pq{l}{k}_{Q^{s-1}}.
\end{align*}
The theorem now follows from \eqref{conv} and after letting $\bm\rho=\bm\alpha^*$.
\end{proof}\label{last}
\end{theorem}

We can use the previous theorem to obtain explicit formulas for the rook numbers of some special types of Ferrers boards studied by Goldman and Haglund \cite{GolHag}. If the columns lengths of the Ferrers board $B(H_{{\textbf{v}}, {\textbf{u}}})$ has column lengths $0,c,2c,\ldots,(|\textbf u|-1)c$, we call $B(H_{{\textbf{v}}, {\textbf{u}}})$ a $c$-\emph{jump board}. If $B(H_{{\textbf{v}}, {\textbf{u}}})$ is an $|\textbf u| \times|\textbf u|$ Ferrers board, then we call $B(H_{{\textbf{v}}, {\textbf{u}}})$ an \emph{Abel board}.

\begin{corollary} Let $n=|\textbf u|$ and let $f_q(n,k)$ be as in Theorem \ref{last}. If $w$ outlines a $c$-jump board, then
\[
R_{s,h;q}[B(H_{{\textbf{v}}, {\textbf{u}}}),k]= f_q(n,n-k)h^k \sum_{i,j,l} (-1)^{j-l} q^{\binom{n-j}{2}({c+s-1)}}\pq{n}{j}_{q^{1-s-c}}\binom{j}{i} \binom{i}{l} \pq{l}{n-k}_{q^{1-s}}.
\]
On the other hand, if $w$ outlines a $d \times |\textbf u|$ Ferrers board, then
\[
R_{s,h;q}[B(H_{{\textbf{v}}, {\textbf{u}}}),k] = f_q(n,n-k)h^k \sum_{i,j,l} (-1)^{j-l} q^{-d(k)-\binom{n-j}{2}} \pq{n}{j}_{q^{1-s}}\binom{j}{i} \binom{i}{l} \pq{l}{n-k}_{q^{1-s}}.
\]
The Abel board is the case where $d=n$.
\end{corollary}

\section{Some Special Cases}\label{special}
In this section, we enumerate special cases of the normal ordering coefficients considered in Sections \ref{first} and \ref{second} for some strings, with $h=1, q=1$. We have verified that the OEIS sequences cited below are indeed the coefficients of the given strings by comparing their known recurrence relations with Equation \eqref{strec}. We have also included sequences whose recurrence relations are not given in the OEIS entry, but whose table of values match the ones we have generated. These sequences are marked with an asterisk, and it remains an open question if these sequences indeed occur as normal ordering coefficients. The name or short description of the sequences are also given when available.

\begin{center}
\begin{tabular}{|c|c|c|c|}
\hline \textbf{String} &amp; $s$ &amp; \textbf{OEIS entry} &amp; \textbf{Name or short description} \\
\hline \hline $(VU)^n$ &amp; $-3$ &amp; \seqnum{A265604} ${}^{*}$  &amp; Inverse Bell transform of quartic factorial numbers \\
\hline $(VU)^n$ &amp; $-2$ &amp; \seqnum{A265605} ${}^{*}$  &amp; Inverse Bell transform of the triple factorial numbers \\ 
\hline $(VU)^n$ &amp; $-1$ &amp; \seqnum{A122848} ${}^{*}$  &amp; Exponential Riordan array $(1,x(1+x/2))$ \\
\hline $(VU)^n$ &amp; $0$ &amp; \seqnum{A008277}  &amp; Stirling numbers of the second kind \\
\hline $(VU)^n$ &amp; $1$ &amp; \seqnum{A132393}  &amp; Stirling numbers of the first kind \\
\hline $(VU)^n$ &amp; $2$ &amp; \seqnum{A001497}  &amp; Coefficients of Bessel polynomials \\
\hline $(VU)^n$ &amp; $3$ &amp; \seqnum{A203412} &amp; \\
\hline $(VU)^n$ &amp; $4$ &amp; \seqnum{A265606}${}^{*}$  &amp; Bell transform of the quartic factorial numbers \\
\hline $(V^2U)^{n}$ &amp; $0$ &amp; \seqnum{A271703}   &amp; unsigned Lah numbers \\
\hline $(V^3U)^{n}$ &amp; $0$ &amp; \seqnum{A035342}${}^{*}$   &amp; convolution matrix of 
\seqnum{A001147} \\
\hline $(V^4U)^{n}$ &amp; $0$ &amp; \seqnum{A035469}  &amp; Bell transform of \seqnum{A007559} without column 0\\
\hline $(V^2U)^{n}$ &amp; $1$ &amp; unsigned \seqnum{A039683} &amp; \\
\hline $(V^3U)^{n}$ &amp; $1$ &amp; unsigned \seqnum{A051141}  &amp; \\
\hline $(V^4U)^{n}$ &amp; $1$ &amp; unsigned \seqnum{A051142}  &amp; \\
\hline $(V^2U)^{n}$ &amp; $2$ &amp; \seqnum{A004747}  &amp; Bell transform of \seqnum{A008544} without column 0.\\
\hline
\end{tabular}
\end{center}

\section{Acknowledgment}
The third author was supported in part by the Department of Science and
Technology (DOST) through its ASTHRD Program. The research was also
supported by the grant MAT-15-1-03 from the Natural Sciences Research
Institute (NSRI), University of the Philippines Diliman. The support of
these institutions are gratefully acknowledged.

\begin{thebibliography}{8}
\bibitem{Bla} P.~Blasiak, Combinatorics of boson normal ordering and some applications, Ph.~D. thesis, University of Paris VI and Polish Academy of Sciences,
2005, \url{http://arxiv.org/abs/quant-ph/0507206}.

\bibitem{DesCakMan} B.~S.~El-Desouky, N.~Caki\'c, and T.~Mansour,
Modified approach to generalized Stirling numbers via differential
operators, {\em Appl.~Math.~Lett.} {\bf 23} (2010), 115--120.

\bibitem{GolHag} J.~Goldman and J.~Haglund, Generalized rook
polynomials, {\em J.~Combin.~Theory~Ser.~A} {\bf 91} (2000), 509--530.

\bibitem{ManSchSha1} T.~Mansour, M.~Schork, and M.~Shattuck, On a new
family of generalized Stirling and Bell numbers, {\em
Electron.~J.~Combin.} {\bf 18} (2011), \#P77.

\bibitem{ManSchSha2} T.~Mansour, M.~Schork, and M.~Shattuck, The
generalized Stirling and Bell numbers revisited, {\em J.~Integer~Sequences}
{\bf 15} (2012), 
\href{https://cs.uwaterloo.ca/journals/JIS/VOL15/Schork/schork2.html}{Article 12.8.3}.

\bibitem{MedLer2} A.~M\'edicis and P.~Leroux, Generalized Stirling
numbers, convolution formulae and $p,q$-analogues, {\em
Canad.~J.~Math.} {\bf 47} (1995), 474--499.

\bibitem{Nav} A.~Navon, Combinatorics and fermion algebra, {\em Nuovo
Cimento} {\bf 16} (1973) 324--330.

\bibitem{Var} A.~Varvak, Rook numbers and the normal ordering problem,
{\em J.~Combin.~Theory~Ser.~A} {\bf 112} (2005), 292--307.
\end{thebibliography}

\bigskip
\hrule
\bigskip

\noindent 2010 {\it Mathematics Subject Classification}:
Primary 05A15; Secondary, 11B65, 11B73.

\noindent \emph{Keywords: } Stirling number, normal ordering, rook theory, Bell number, symmetric function.


\bigskip
\hrule
\bigskip

\noindent (Concerned with sequences 
\seqnum{A001147},
\seqnum{A001497},
\seqnum{A004747},
\seqnum{A007559},
\seqnum{A008277},
\seqnum{A008544},
\seqnum{A035342},
\seqnum{A035469},
\seqnum{A039683},
\seqnum{A051141},
\seqnum{A051142},
\seqnum{A122848},
\seqnum{A132393},
\seqnum{A203412},
\seqnum{A265604},
\seqnum{A265605},
\seqnum{A265606}, and
\seqnum{A271703}.)

\bigskip
\hrule
\bigskip

\vspace*{+.1in}
\noindent
Received August 22 2016; revised versions received October 16 2016;
January 3 2017.  Published in
{\it Journal of Integer Sequences}, January 14 2017.

\bigskip
\hrule
\bigskip

\noindent
Return to
\htmladdnormallink{Journal of Integer Sequences home page}{http://www.cs.uwaterloo.ca/journals/JIS/}.
\vskip .1in

\end{document}
</pre></body></html>