Some applications of the chromatic polynomials

Abstract: The chromatic polynomials are studied by several authors and have important applications in different frameworks, specially, in graph theory and enumerative combinatorics. The aim of this work is to establish some properties of the coefficients of the chromatic polynomial of a graph. Three applications on restricted Stirling numbers of the second kind are given.

\section{Introduction}
The chromatic polynomial was introduced by Birkhoff \cite{bir1} and studied
later by Whitney \cite{whi1,whi2}, Birkhoff and Lewis \cite{bir2}, Read \cit
{rea} and several other authors.\ The chromatic polynomial of a graph can be
used as a tool to find the number of possible partitions of a finite set
under some particular restraints such that the Stirling numbers of the
second kind \cite{whit}.
For a given graph $G=\left( V,E\right) $ of order $n$ and $\lambda \in
\mathbb{N}
,$ a mapping $f:V\rightarrow \left \{ 1,2,\ldots ,\lambda \right \} $ is
called a $\lambda $-coloring of $G$ if $f(u)\neq f(v)$ whenever the vertices
$u$ \ and $v$ are adjacent in $G$. The $\lambda $-colorings $f$ and $g$ of
G $ are regarded as distinct if $f(x)\neq g(x)$ for some $x$ in $G.$ The
chromatic polynomial $P(G,\lambda )$ counts the number of (proper) $\lambda
-colorings of $G.$ For example, it is known that $P\left( O_{n},\lambda
\right) =\lambda ^{n},$ $P(K_{n},\lambda )=\left( \lambda \right) _{n}$ and
P\left( T_{n},\lambda \right) =\lambda \left( \lambda -1\right) ^{n-1},$
n\geq 1,$ where $O_{n}$ is a graph of order $n$ and without edges, $K_{n}$
is the complete graph of order $n,$ $T_{n}$ is a tree of order $n$ and
\left( \lambda \right) _{n}=\lambda \left( \lambda -1\right) \cdots \left(
\lambda -n+1\right) $ if $n\geq 1$ and $\left( \lambda \right) _{0}=1.$ More
generally, the chromatic polynomial of $G$ can be written as $P(G,\lambda
)=\sum \limits_{i=\chi \left( G\right) }^{n}\alpha _{i}\left( G\right)
\left( \lambda \right) _{i},$ see \cite[Thm. 1.4.1]{don}, where $\alpha
_{i}\left( G\right) $ is the number of ways of partitioning $V$ into $i$
independent sets and $\chi \left( G\right) $ is the chromatic number. For
use later, recall that, if $G_{1}=\left( V_{1},E_{1}\right) $ and
G_{2}=\left( V_{2},E_{2}\right) $ are graphs on disjoint sets of vertices,
their union is defined by the graph $G_{1}\cup G_{2}=\left( V_{1}\cup
V_{2},E_{1}\cup E_{2}\right) $ and $P(G_{1}\cup G_{2},\lambda
)=P(G_{1},\lambda )P(G_{2},\lambda ),$ see for example \cite[Sec.\ 1.2]{don}.
In this paper, for a given graph $H,$ we present some properties for the
families of graphs $O_{n}\cup H,$ $K_{n}\cup H$ and $T_{n}\cup H.$ In the
next section, we give some recurrence relations for the coefficients $\alpha
_{k}\left( O_{n}\cup H\right) ,$ $\alpha _{k}\left( K_{n}\cup H\right) $ and
$\alpha _{k}\left( T_{n}\cup H\right) $ and some results on log-concavity and P\'{o}lya-frequency for sequences  related to these coefficients. In the three last sections, we present three applications on restricted Stirling numbers of the second kind.

\section{Recurrence relations and some consequences}
Let $H$ be any graph of $h$ vertices, $O_{n}$ be the graph of $n$ $\left(
\geq 1\right) $ vertices and no edges and let $O_{0}$ be the graph with no
vertices, $K_{n}$ be the complete graph of $n\  \left( \geq 1\right) $
vertices with $K_{0}$ be a graph with no vertices and $T_{n}$ be a tree of
n\  \left( \geq 1\right) $ vertices with $T_{0}$ be a graph with no vertices.
In this section, we give some recurrence relations for the coefficients
\alpha _{k}\left( O_{n}\cup H\right) ,$ $\alpha _{k}\left( K_{n}\cup
H\right) ,$ $\alpha _{k}\left( T_{n}\cup H\right) $ and some of their consequences.
\begin{theorem}
\label{TT1}Let $n,$ $s,$ $k$ be nonnegative integers with $\ 0\leq s\leq n.$
Then, $\alpha _{k}\left( O_{n}\cup H\right) =0$ if $k<\chi \left( H\right) $
or $k>n+h$ and for $\chi \left( H\right) \leq k\leq n+h$\ we hav
\begin{equation*}
\alpha _{k}\left( O_{n}\cup H\right) =\sum_{j=\chi \left( H\right)
}^{k}{s+j \brace k }_{\!\!j}\alpha _{j}\left( O_{n-s}\cup H\right) .
\end{equation*
In particular, for $s=n,$ we ge
\begin{equation*}
\alpha _{k}\left( O_{n}\cup H\right) =\sum_{j=\chi \left( H\right)
}^{k}{n+j \brace k }_{\!\!j}\alpha _{j}\left( H\right) ,
\end{equation*
and, for $s=1,$ we ge
\begin{equation*}
\alpha _{k}\left( O_{n}\cup H\right) =k\alpha _{k}\left( O_{n-1}\cup
H\right) +\alpha _{k-1}\left( O_{n-1}\cup H\right) ,\  \ n\geq 1,
\end{equation*
where the numbers ${n \brace k }_{\!\!r}$ are the $r$-Stirling numbers of
the second kind.
\end{theorem}

\begin{proof}
From \cite[Sec. 1.2]{don} we have $P\left( O_{n}\cup H,\lambda \right)
=\lambda ^{s}P\left( O_{n-s}\cup H,\lambda \right) ,\  \  \ 0\leq s\leq n.
\newline
Then, since the graph $O_{n}\cup H$ is of order $n+h$ an
\begin{equation*}
\chi \left( O_{n}\cup H\right) =\max \left( \chi \left( O_{n}\right) ,\chi
\left( H\right) \right) =\max \left( 1,\chi \left( H\right) \right) =\chi
\left( H\right) ,
\end{equation*
we may state that $\alpha _{k}\left( O_{n}\cup H\right) =0$ if $k<\chi
\left( H\right) $ or $k>n+h,$ and otherwise, we hav
\begin{equation*}
\underset{k=\chi \left( H\right) }{\overset{n+h}{\sum }}\alpha _{k}\left(
O_{n}\cup H\right) \left( \lambda \right) _{k}=\underset{j=\chi \left(
H\right) }{\overset{n+h}{\sum }}\alpha _{j}\left( O_{n-s}\cup H\right)
\lambda ^{s}\left( \lambda \right) _{j}.
\end{equation*
On using the known identity $\left( z+j\right) ^{s}=\sum_{k=0}^{s}{n+j \brace k+j }_{\!\!j}\left( z\right) _{k},$
see \cite{bro}, we can writ
\begin{equation*}
\lambda ^{s}\left( \lambda \right) _{j}=\left( \lambda \right) _{j}\left(
\lambda -j+j\right) ^{s}=\overset{s}{\underset{k=0}{\sum }}{n+j \brace k+j }_{\!\!j}
\left( \lambda \right) _{j}\left( \lambda -j\right) _{k}
\end{equation*
and since $\left( \lambda \right) _{j}\left( \lambda -j\right) _{k}=\left(
\lambda \right) _{k+j}$ we get $\lambda ^{s}\left( \lambda \right) _{j}
\overset{s}{\underset{k=0}{\sum }}{n+j \brace k+j }_{\!\!j}\left( \lambda
\right) _{k+j}=\overset{s+j}{\underset{k=j}{\sum }}{s+j \brace k }_{\!\!j}\left( \lambda \right) _{k}.$ The
\begin{eqnarray*}
\underset{k=\chi \left( H\right) }{\overset{n+h}{\sum }}\alpha _{k}\left(
O_{n}\cup H\right) \left( \lambda \right) _{k} &=&\underset{j=\chi \left(
H\right) }{\overset{n+h}{\sum }}\alpha _{j}\left( O_{n-s}\cup H\right)
\overset{s+j}{\underset{k=j}{\sum }}{s+j \brace k }_{\!\!j}\left( \lambda
\right) _{k} \\
&=&\overset{s+j}{\underset{k=j}{\sum }}\left( \lambda \right)
_{k}\sum_{j=\chi \left( H\right) }^{k}\alpha _{j}\left( O_{n-s}\cup H\right)
{s+j \brace k }_{\!\!j}.
\end{eqnarray*
So, we get $\alpha _{k}\left( O_{n}\cup H\right) =\sum_{j=\chi \left(
H\right) }^{k}\alpha _{j}\left( O_{n-s}\cup H\right) {s+j \brace k }_{\!\!j}.$
\end{proof}

\begin{theorem}
\label{TT2}Let $n,$ $s,$ $k$ be nonnegative integers with$\ 0\leq s\leq n.$
Then, $\alpha _{k}\left( K_{n}\cup H\right) =0$ if $k<\max \left( n,\chi
\left( H\right) \right) $ or $k>n+h$ and for $\max \left( n,\chi \left(
H\right) \right) \leq k\leq n+h$\ we hav
\begin{equation*}
\alpha _{k}\left( K_{n}\cup H\right) =\underset{j=\max \left( n-s,\chi \left(
H\right) \right) }{\overset{k}{\sum }}\binom{s}{k-j}\left( j+s-n\right)
_{s+j-k}\alpha _{j}\left( K_{n-s}\cup H\right) .
\end{equation*
In particular, for $s=n,$ we ge
\begin{equation*}
\alpha _{k}\left( K_{n}\cup H\right) =\underset{j=\chi \left(
H\right)  }{\overset{k}{\sum }}\binom{n}{k-j}\left( j\right)
_{n+j-k}\alpha _{j}\left( H\right) ,
\end{equation*
and, for $s=1,$ we ge
\begin{equation*}
\alpha _{k}\left( K_{n}\cup H\right) =\left( k-n+1\right) \alpha _{k}\left(
K_{n-1}\cup H\right) +\alpha _{k-1}\left( K_{n-1}\cup H\right) ,\  \ n\geq 1.
\end{equation*}
\end{theorem}

\begin{proof}
From \cite[Sec. 1.2]{don} we have $P\left( K_{n}\cup H,\lambda \right)
=\left( \lambda -n+1\right) P\left( K_{n-1}\cup H,\lambda \right) .$ Hence
\begin{equation*}
P\left( K_{n}\cup H,\lambda \right) =\left( \lambda -n+s\right) _{s}P\left(
K_{n-s}\cup H,\lambda \right) ,\  \ 0\leq s\leq n.
\end{equation*
Then, since the graph $K_{n}\cup H$ is of order $n+h$ an
\begin{equation*}
\chi \left( K_{n}\cup H\right) =\max \left( \chi \left( K_{n}\right) ,\chi
\left( H\right) \right) =\max \left( n,\chi \left( H\right) \right) ,
\end{equation*
we state that $\alpha _{k}\left( K_{n}\cup H\right) =0$ if $k<\max \left(
n,\chi \left( H\right) \right) $ or $k>n+h,$ and otherwise, we hav
\begin{equation*}
\underset{k=\max \left( n,\chi \left( H\right) \right) }{\overset{n+h}{\sum
}\alpha _{k}\left( K_{n}\cup H\right) \left( \lambda \right) _{k}=\underset
j=\max \left( n-s,\chi \left( H\right) \right) }{\overset{n-s+h}{\sum }}\alpha
_{j}\left( K_{n-s}\cup H\right) \left( \lambda -n+s\right) _{s}\left(
\lambda \right) _{j}.
\end{equation*
By the fact that the sequence of polynomials $((\lambda)_n, n\geq0)$ is of binomial type, the expression $\left( \lambda -n+s\right) _{s}\left(\lambda \right) _{j}$ can be written as
\begin{eqnarray*}
\left( \lambda -n+s\right) _{s}\left( \lambda \right) _{j} &=&\left( \lambda
-j+j+s-n\right) _{s}\left( \lambda \right) _{j} \\
&=&\underset{k=0}{\overset{s}{\sum }}\binom{s}{k}\left( j+s-n\right)
_{s-k}\left( \lambda -j\right) _{k}\left( \lambda \right) _{j} \\
&=&\underset{k=0}{\overset{s}{\sum }}\binom{s}{k}\left( j+s-n\right)
_{s-k}\left( \lambda \right) _{k+j} \\
&=&\underset{k=j}{\overset{s+j}{\sum }}\binom{s}{k-j}\left( j+s-n\right)
_{s+j-k}\left( \lambda \right) _{k}.
\end{eqnarray*
So, it results that
\begin{eqnarray*}
&&\underset{k=\max \left( n,\chi \left( H\right) \right) }{\overset{n+h}
\sum }}\alpha _{k}\left( K_{n}\cup H\right) \left( \lambda \right) _{k} \\
&=&\underset{j=\max \left( n-s,\chi \left( H\right) \right) }{\overset{n-s+h}
\sum }}\alpha _{j}\left( K_{n-s}\cup H\right) \underset{k=j}{\overset{s+j}
\sum }}\binom{s}{k-j}\left( j+s-n\right) _{s+j-k}\left( \lambda \right) _{k}
\\
&=&\underset{k=\max \left( n,\chi \left( H\right) \right) }{\overset{n+h}
\sum }}\left( \lambda \right) _{k}\underset{j=\max \left( n-s,\chi \left(
H\right) \right) }{\overset{k}{\sum }}\binom{s}{k-j}\left( j+s-n\right)
_{s+j-k}\alpha _{j}\left( K_{n-s}\cup H\right).
\end{eqnarray*
The identification of the coefficients of $\left( \lambda \right) _{k}$ complete the proof.
\end{proof}

\begin{theorem}
\label{TT3}Let $n,$ $s,$ $k$ be nonnegative integers with$\ 0\leq s\leq n.$
Then
\begin{equation*}
\alpha _{k}\left( T_{n}\cup H\right) =0\text{ \ for \ }k<\max \left(
2-\delta _{\left( n=1\right) }-2\delta _{\left( n=0\right) },\chi \left(
H\right) \right) \text{ \ or \ }k>n+h,
\end{equation*
and, for $\max \left( 2-\delta _{\left( n=1\right) }-2\delta _{\left(
n=0\right) },\chi \left( H\right) \right) \leq k\leq n+h$\ we hav
\begin{equation*}
\alpha _{k}\left( T_{n}\cup H\right) =\underset{j=\max \left( k-s,0\right) }
\overset{k}{\sum }}{s+j-1 \brace k-1 }_{\!\!j-1}\alpha _{j}\left(
T_{n-s}\cup H\right) .
\end{equation*
In particular, for $s=n,$ we ge
\begin{equation*}
\alpha _{k}\left( T_{n}\cup H\right) =\underset{j=1}{\overset{k}{\sum }
{s+j-1 \brace k-1 }_{\!\!j-1}\alpha _{j}\left( H\right)
\end{equation*
and, for $s=1,$ we ge
\begin{equation*}
\alpha _{k}\left( T_{n}\cup H\right) =\left(k-1 \right)\alpha _{k}\left( T_{n-1}\cup
H\right) +\alpha _{k-1}\left( T_{n-1}\cup H\right) ,\  \ n\geq 1.
\end{equation*}
\end{theorem}

\begin{proof}
The graph $T_{n}\cup H$ is of order $n+h$ and we hav
\begin{equation*}
\chi \left( T_{n}\cup H\right) =\max \left( \chi \left( T_{n}\right) ,\chi
\left( H\right) \right) =\max \left( 2-\delta _{\left( n=1\right) }-2\delta
_{\left( n=0\right) },\chi \left( H\right) \right) ,
\end{equation*
we conclude that $\alpha _{k}\left( T_{n}\cup H\right) =0$ if $k<\max \left(
2-\delta _{\left( n=1\right) }-2\delta _{\left( n=0\right) },\chi \left(
H\right) \right) $ or $k>n+h,$ and otherwise, we hav
\begin{equation*}
P\left( T_{n}\cup H,\lambda \right) =P\left( T_{n},\lambda \right) P\left(
H,\lambda \right) =\lambda \left( \lambda -1\right) ^{n-1}P\left( H,\lambda
\right) ,
\end{equation*
which gives $P\left( T_{n}\cup H,\lambda \right) =\left( \lambda -1\right)
^{s}P\left( T_{n-s}\cup H,\lambda \right) ,\  \ 0\leq s\leq n.$ \newline
This relation is equivalent t
\begin{equation*}
\underset{k=0}{\overset{n+h}{\sum }}\alpha _{k}\left( T_{n}\cup H\right)
\left( \lambda \right) _{k}=\underset{j=0}{\overset{n+h-s}{\sum }}\alpha
_{j}\left( T_{n-s}\cup H\right) \left( \lambda -1\right) ^{s}\left( \lambda
\right) _{j}.
\end{equation*
Similarly to $\lambda ^{s}\left( \lambda \right) _{j}$ (see proof of
Theorem \ref{TT1}), we get
\begin{equation*}
\left( \lambda -1\right) ^{s}\left( \lambda \right) _{j}=\underset{k=j}
\overset{s+j}{\sum }}{s+j-1 \brace k-1 }_{\!\!j-1}\left( \lambda \right)
_{k}.
\end{equation*}
Then, the last equality becomes
\begin{eqnarray*}
\underset{k=0}{\overset{n+h}{\sum }}\alpha _{k}\left( T_{n}\cup H\right)
\left( \lambda \right) _{k} &=&\underset{j=0}{\overset{n+h-s}{\sum }}\alpha
_{j}\left( T_{n-s}\cup H\right) \underset{k=j}{\overset{s+j}{\sum }}{s+j-1 \brace k-1 }_{\!\!j-1}\left( \lambda \right) _{k} \\
&=&\underset{k=0}{\overset{n+h}{\sum }}\left( \lambda \right) _{k}\underset
j=\max \left( k-s,0\right) }{\overset{k}{\sum }}{s+j-1 \brace k-1 }_{\!\!j-1}\alpha _{j}\left( T_{n-s}\cup H\right) ,
\end{eqnarray*
which gives thus the result.
\end{proof}

\noindent To give the following proposition, let us recall some definitions
and results on log-concavity, P\'{o}lya-frequency and $q$-log-convexity.
Indeed, let $u_{0},$ $u_{1},$ $u_{2},\ldots ,$ be a sequence of nonnegative
real numbers. It is log-concave (LC) if $u_{i-1}u_{i+1}\leq u_{i}^{2}$ for
all $i>0,$ and, it is called a P\'{o}lya-frequency sequence (or a PF
sequence) if all minors of the matrix $A=(u_{i-j})_{i,j\geq 0}$ have
nonnegative determinants (where $u_{k}=0$ if $k<0$), for more information
see \cite{kar}. A sequence of real polynomials $\left( P_{n}(q),n\geq
0\right) $ is called $q$-log-convex if the polynomial
P_{n}(q)^{2}-P_{n-1}(q)P_{n+1}(q)$ has nonnegative coefficients for all
n\geq 1,$ see \cite{sag1,sag2,zho}. Some known results on such sequences are
given as follows. \newline
Let $\left( T\left( n,k\right) ,\ n,k\geq 0\right) $ be sequence of
nonnegative numbers satisfying the recurrence
\begin{equation*}
T\left( n,k\right) =\left( a_{1}n+a_{2}k+a_{3}\right) T\left( n-1,k\right)
+\left( b_{1}n+b_{2}k+b_{3}\right) T\left( n-1,k-1\right), \ n\geq k\geq 1,
\end{equation*
with $T\left( n,k\right) =0$ unless $0\leq k\leq n,$ $T\left(0,0\right)>0,$ \ $a_{1}\geq
0,$ $a_{1}+a_{2}\geq 0,\ a_{1}+a_{3}\geq 0$ and $b_{1}\geq 0,$
b_{1}+b_{2}\geq 0,$ $b_{1}+b_{2}+b_{3}\geq 0.$ It is shown in \cite[Thm. 2
{kur} that, for each fixed $n,$ the sequence $\left( T\left( n,k\right) ,\
0\leq k\leq n\right) $ is log-concave. If further we have $a_{2}b_{1}\geq
a_{1}b_{2}$ \ and \ $a_{2}\left( b_{1}+b_{2}+b_{3}\right) \geq \left(
a_{1}+a_{3}\right) b_{2},$ this sequence is P\'{o}lya frequency sequence
\cite[Cor. 3]{wan} and further, by setting $T_{n}\left( q\right) =\overset{n
{\underset{k=0}{\sum }}T\left( n,k\right) q^{k},$ if
\begin{equation*}
\left( a_{2}b_{1}-a_{1}b_{2}\right) n+a_{2}b_{2}k+a_{2}b_{3}-a_{3}b_{2}\geq
\text{ for }0<k\leq n,
\end{equation*
then, the sequence of polynomials $\left( T_{n}\left( q\right) ,\ n\geq
0\right) $ is $q$-log-convex \cite[Thm. 4.1]{liu}.

\begin{proposition}
\label{P1}Let $\left( U\left( n,k\right) ,\ n,k\geq 0\right) ,$  $\left( V\left( n,k\right) ,\ n,k\geq 0\right) ,$ and $\left( W\left( n,k\right) ,\ n,k\geq 0\right) ,$ be sequences of nonnegative numbers with
$U\left( n,k\right)= V\left( n,k\right)= W\left( n,k\right)=0$ when  $ k> n$ and for $0\leq k\leq n,$ these sequences are defined by
\begin{eqnarray*}
U\left( n,k\right) & =&\alpha _{k+h}\left( O_{n}\cup H\right) ,\\
V\left( n,k\right)  &=&\alpha _{k+h}\left( T_{n}\cup H\right) ,\\
W\left( n,k\right) & =&\alpha _{k+h}\left( K_{n}\cup H\right)
\end{eqnarray*
and let
\begin{eqnarray*}
U_{n}\left( q\right) &=&\overset{n}{\underset{k=0}{\sum }}U\left( n,k\right) q^{k}, \\
V_{n}\left( q\right) &=&\overset{n}{\underset{k=0}{\sum }}V\left( n,k\right) q^{k}.
\end{eqnarray*
Then, the sequences $\left( U\left( n,k\right) ,\ 0\leq k\leq n\right) $ and
\left( V\left( n,k\right) ,\ 0\leq k\leq n\right) $ are log-concave and
\'{o}lya frequency sequences, the sequences of real polynomials $\left(
U_{n}(q),n\geq 0\right) $ and $\left( V_{n}(q),n\geq 0\right) $ are $q
-log-convex sequences and the sequence $\left( W\left( n,k\right) ,\ 0\leq
k\leq n\right) $ is a P\'{o}lya frequency sequence.
\end{proposition}

\begin{proof}
We have $U\left(0,0\right) =V\left(0,0\right) = W\left(0,0\right) =\alpha _{h}\left(H\right)=1 $ and for $n\geq 1$, Theorems \ref{TT1}, \ref{TT2} and \ref{TT3} imply
\begin{eqnarray*}
U\left(n,k\right) &=&U\left( n-1,k-1\right) +\left( k+h\right) U\left( n-1,k\right) ,\\
V\left(n,k\right) &=&V\left( n-1,k-1\right) +\left( k+h-1\right) V\left( n-1,k\right) , \\
W\left(n,k\right) &=&W\left( n-1,k-1\right) +\left( k-n+h+1\right) W\left(n-1,k\right) .
\end{eqnarray*
So, the log-concavity follows from \cite[Thm. 2]{kur}, P\'{o}lya frequency
follows from \cite[Cor. 3]{wan} and $q$-log-convexity follows from \cite
Thm. 4.1]{liu}.
\end{proof}

\section{Application to the graph $K_{r_{_{1}}}\cup \cdots \cup K_{r_{_{p}}}\cup O_{n}$}

Let $p\geq 1$ be an integer. We consider the graph $G_{n,\mathbf{r
_{p}}=K_{r_{_{1}}}\cup \cdots \cup K_{r_{_{p}}}\cup O_{n}$ of order
n+r_{1}+\cdots +r_{p}$ and chromatic number
\begin{equation*}
\chi \left( G_{n,\mathbf{r}_{p}}\right) =\max \left( \chi \left(
K_{r_{_{1}}}\right) ,\ldots ,\chi \left( K_{r_{_{p}}}\right) ,\chi \left(
O_{n}\right) \right) =\max \left( r_{1},\ldots ,r_{p},1\right) .
\end{equation*
First of all, recall the definition of the $\left( r_{1},\ldots
,r_{p}\right) $-Stirling number of the second kind introduced by Mihoubi et al. in \cite{mih,maa} .

\begin{definition}
Let $R_{1},\ldots ,R_{p}$ be subsets of the set $\left[ n\right] $ with
\left \vert R_{i}\right \vert =r_{i}$ and $R_{i}\cap R_{j}=\varnothing $ for
all $i,j=1,\ldots ,p,$ $i\neq j.$ The $\left( r_{1},\ldots ,r_{p}\right)
-Stirling number of the second kind, $p\geq 1,$ denoted by ${ n \brace k }_{\!\!r_{1},\ldots,r_{p}},$
counts the number of partitions of the set
\left[ n\right] :=\left \{ 1,2,\ldots ,n\right \} $ into $k$ non-empty
subsets such that the elements of each of the $p$ sets $R_{1},\ldots ,R_{p}$
are in distinct subsets.
\end{definition}

\noindent From this definition, one can see easily that these numbers satisf
\begin{align*}
{ n \brace k }_{\!\!r_{1},\ldots,r_{p}}& =0,\text{ \  \  \ }n<r_{1}+\cdots
+r_{p}\text{ or }k<\max \left( r_{1},\ldots ,r_{p},1\right) , \\
{ n \brace k }_{\!\!r_{1},\ldots,r_{p}}& ={n\brace k}_{\!\!r_{p}
\text{ if }r_{1},\ldots ,r_{p-1}\in \left \{ 0,1\right \} , \\
{ n \brace k }_{\!\!r_{1},\ldots,r_{p}}& ={ n \brace k }_{\!\!r_{\sigma
\left( 1\right) },\ldots ,r_{\sigma \left( p\right) }}\text{ for all
permutations }\sigma \text{ on the set }\left \{ 1,\ldots ,p\right \} .
\end{align*
Furthermore, for $r_{1}\leq r_{2}\leq \cdots \leq r_{p},$ the coefficient
\alpha _{k}\left( G_{n,\mathbf{r}_{p}}\right) $ represents the number of
ways of partitioning the set of
n+\left \vert \mathbf{r}_{p}\right \vert $ vertices of $G_{n,\mathbf{r
_{p}} $ into $k$ independent subsets, and by the definition of the graph
G_{n,\mathbf{r}_{p}},$ the elements of each of the $p$ subgraphs
K_{r_{_{1}}},\ldots ,K_{r_{_{p}}}$ must be in distinct subsets. So, this
number is exactly ${n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!\mathbf{r}_{p}}.$
Then, we may state tha
\begin{equation*}
\alpha _{k}\left( G_{n,\mathbf{r}_{p}}\right) ={n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!\mathbf{r}_{p}} .
\end{equation*
In particula
\begin{equation*}
\alpha _{k}\left( K_{r}\cup O_{n}\right) ={n+r \brace k}_{\!\!r} \text{ \
and \ }\alpha _{k}\left( O_{n}\right) ={n\brace k}.
\end{equation*
Now, we present new proofs for some properties of the $\left( r_{1},\ldots
,r_{p}\right) $-Stirling numbers of the second kind and other
results.

\begin{theorem}
\label{T1}For $r_{1}\leq r_{2}\leq \cdots \leq r_{p},$ the polynomial
\left( z+r_{p}\right) _{r_{1}}\cdots \left( z+r_{p}\right) _{r_{p-1}}\left(
z+r_{p}\right) ^{n}$ can be written in the basis $\left \{ \left( \lambda
\right) _{k},\ k=0,1,\ldots,n+\left \vert \mathbf{r}_{p-1}\right \vert \right \} $ as follow
\begin{equation*}
\left( z+r_{p}\right) _{r_{1}}\cdots \left( z+r_{p}\right) _{r_{p-1}}\left(
z+r_{p}\right) ^{n}=\sum_{k=0}^{n+\left \vert \mathbf{r}_{p-1}\right \vert
}{n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!\mathbf{r}_{p}}\left( z\right) _{k},
\end{equation*
where $\mathbf{r}_{p}:=\left( r_{1},\ldots ,r_{p}\right) $ and $\left \vert
\mathbf{r}_{p}\right \vert :=r_{1}+\cdots +r_{p}.$ In particular, we have
\begin{equation*}
\left( z+r\right) ^{n}=\sum_{k=0}^{n}{n+r \brace k+r}_{\!\!r}\left(
z\right) _{k}.
\end{equation*}
\end{theorem}

\begin{proof}
From \cite[Sec. 1.2]{don}, the chromatic polynomial of the graph $G_{n
\mathbf{r}_{p}}$ i
\begin{eqnarray*}
P\left( G_{n,\mathbf{r}_{p}},z+r_{p}\right) &=&P\left( O_{n},z+r_{p}\right)
P\left( K_{r_{1}},z+r_{p}\right) \cdots P\left( K_{r_{p}},z+r_{p}\right) \\
&=&\left( z+r_{p}\right) ^{n}\left( z+r_{p}\right) _{r_{1}}\left(
z+r_{p}\right) _{r_{2}}\cdots \left( z+r_{p}\right) _{r_{p}},
\end{eqnarray*
and, by definition of the chromatic polynomial, we hav
\begin{eqnarray*}
P\left( G_{n,\mathbf{r}_{p}},z+r_{p}\right)& =&\sum_{k=0}^{n+\left \vert
\mathbf{r}_{p-1}\right \vert }{n+\left \vert \mathbf{r
_{p}\right \vert \brace k+r_{p}}_{\!\!\mathbf{r}_{p}}\left( z+r_{p}\right)
_{k+r_{p}} \\ &=&\sum_{k=0}^{n+\left \vert \mathbf{r}_{p-1}\right \vert }
{n+\left \vert \mathbf{r}_{p}\right \vert \brace k+r_{p}}_{\!\!\mathbf{r
_{p}}\left( z+r_{p}\right) _{r_{p}}\left( z\right) _{k}.
\end{eqnarray*
which complete the proof.
\end{proof}

\begin{remark}
\label{R1}From the following identity (see \cite{tom} or \cite[pp. 102]{don
) we hav
\begin{equation*}
\alpha _{k}\left( G_{n,\mathbf{r}_{p}}\right) =\frac{1}{k!
\sum_{j=0}^{k}\left( -1\right) ^{k-j}\binom{k}{j}P\left( G_{n,\mathbf{r
_{p}},j\right) ,
\end{equation*
or equivalently
\begin{eqnarray*}
{n+\left \vert \mathbf{r}_{p}\right \vert \brace k+r_{p}}_{\!\!\mathbf{
}_{p}} &=&\frac{1}{\left( k+r_{p}\right) !}\sum_{j=r_{p}}^{k+r_{p}}\left(
-1\right) ^{k+r_{p}-j}\binom{k+r_{p}}{j}\left( j\right) _{r_{1}}\cdots
\left( j\right) _{r_{p}}j^{n} \\
&=&\frac{1}{\left( k+r_{p}\right) !}\sum_{j=0}^{k}\left( -1\right) ^{k-j
\binom{k+r_{p}}{j+r_{p}}\left( j+r_{p}\right) _{r_{1}}\cdots \left(
j+r_{p}\right) _{r_{p}}\left( j+r_{p}\right) ^{n} \\
&=&\frac{1}{k!}\sum_{j=0}^{k}\left( -1\right) ^{k-j}\binom{k}{j}\left(
j+r_{p}\right) _{r_{1}}\cdots \left( j+r_{p}\right) _{r_{p-1}}\left(
j+r_{p}\right) ^{n}.
\end{eqnarray*
So, the two initial values of these numbers are given b
\begin{eqnarray*}
{n+\left \vert \mathbf{r}_{p}\right \vert \brace r_{p}}_{\!\!\mathbf{r
_{p}} &=&\left( r_{p}\right) ^{n}\left( r_{p}\right) _{r_{1}}\cdots \left(
r_{p}\right) _{r_{p-1}}, \\
{n+\left \vert \mathbf{r}_{p}\right \vert \brace r_{p}+1}_{\!\!\mathbf{
}_{p}} &=&\left( r_{p}+1\right) _{r_{1}}\cdots \left( r_{p}+1\right)
_{r_{p-1}}\left( r_{p}+1\right) ^{n}-\left( r_{p}\right) _{r_{1}}\cdots
\left( r_{p}\right) _{r_{p-1}}\left( r_{p}\right) ^{n}.
\end{eqnarray*}
\end{remark}

\noindent Now, set $K_{r_{_{i}}}\cup H:=K_{r_{_{1}}}\cup \cdots \cup K_{r_{_{p}}}\cup O_{n} \ \left(i=1,\ldots ,p\right)$ in Theorem \ref{TT2} to obtain:

\begin{corollary}
Let $n, \ s$ and $k$ be nonnegative integers with $n\geq \left \vert \mathbf{r}_{p}\right \vert $.
Then, for a fixed $ i$ in the set $\left\{1,\ldots ,p\right\}$,
if $r_{i}\geq 1$  we hav
\begin{equation*}
{n \brace k}_{\!\!\mathbf{r}_{p}}={n-1 \brace k-1}_{\!\!\mathbf{r
_{p}-e_{i}}+\left( k+1-r_{i}\right) {n-1 \brace k}_{\!\!\mathbf{r
_{p}-e_{i}},
\end{equation*
where $\mathbf{e}_{i}$ denote the i-$th$ vector of the canonical basis of
\mathbb{R}^{p}.$
\end{corollary}

\noindent A particular case of the last corollary is when $p=1$ we obtain
the following known identit
\begin{equation*}
{n \brace k}_{\!\!r}={n-1 \brace k-1}_{\!\!r-1}+\left( k+1-r\right){n-1 \brace k}_{\!\!r-1}.
\end{equation*}

\noindent By setting $ H:=K_{r_{_{1}}}\cup \cdots \cup K_{r_{_{p}}}$ in Theorem \ref{TT1} we obtain:

\begin{corollary}
\label{C1}Let $n,\ k$ be integers such that $r_{p}\leq k\leq n$ and $n\geq
\left \vert \mathbf{r}_{p}\right \vert .$ We hav
\begin{equation*}
{n \brace k}_{\!\!\mathbf{r}_{p}}=\sum_{j=r_{p}}^{k}{s+j \brace k}_{\!\!j}
{n-s \brace j}_{\!\!\mathbf{r}_{p}},\  \ 0\leq s\leq n-\left \vert \mathbf{r}_{p}\right \vert .
\end{equation*
In particular, for $s=1,$ these numbers obey to the following recurrence relation
\begin{equation*}
{n \brace k}_{\!\!\mathbf{r}_{p}}={n-1 \brace k-1}_{\!\!\mathbf{r
_{p}}+k{n-1 \brace k}_{\!\!\mathbf{r}_{p}},
\end{equation*
and for $s=n-\left \vert \mathbf{r}_{p}\right \vert $ we ge
\begin{equation*}
{n \brace k}_{\!\!\mathbf{r}_{p}}=\sum_{j=r_{p}}^{k}
{n-\left \vert \mathbf{r}_{p}\right \vert +j \brace k}_{\!\!j}{\left
\vert \mathbf{r}_{p}\right \vert  \brace j}_{\!\!\mathbf{r}_{p}}.
\end{equation*}
\end{corollary}

\noindent In particular, for $p=s=1$ in Corollary \ref{C1}, we obtain the known identit
\begin{equation*}
{n \brace k}_{\!\!r}={n-1 \brace k-1}_{\!\!r}+k{n-1 \brace k}_{\!\!r},
\end{equation*
and for $s=n-\left \vert \mathbf{r}_{p}\right \vert ,$ $p=2$ and $\left(
r_{1},r_{2}\right) =\left( 1,r\right) $ in Corollary \ref{C1}, we obtain
\begin{equation*}
{n \brace k}_{\!\!r}={n \brace k}_{\!\!1,r}=\left( r+1\right)
{n-1 \brace k}_{\!\!r}+{n \brace k}_{\!\!r+1}.
\end{equation*
Now, upon using Proposition \ref{P1}, we may state the following corollary.
\begin{corollary}
\label{C3}For $0\leq k\leq n$ let
\begin{equation*}
U\left( n,k\right) ={n+\left \vert \mathbf{r}_{p}\right \vert
\brace k+\left \vert \mathbf{r}_{p}\right \vert }_{\!\!\mathbf{r}_{p}}.
\end{equation*
Then, the sequence $\left( U\left( n,k\right) ,\ 0\leq k\leq n\right) $ is
log-concave and P\'{o}lya frequency sequence, and, the sequence of
polynomials $\left( U_{n}\left( q\right) ,n\geq 0\right) $ defined by
\begin{equation*}
U_{n}\left( q\right) =\overset{n}{\underset{k=0}{\sum }}{n+\left \vert \mathbf{r}_{p}\right \vert
\brace k+\left \vert \mathbf{r}_{p}\right \vert }_{\!\!\mathbf{r}_{p}}q^{k}
\end{equation*
is a $q$-log-convex sequence.
\end{corollary}

\section{Application to the graph $K_{r_{1},\ldots,r_{p}}\cup O_{n}$}

Let $p\geq 2$ be an integer. For a second application of the chromatic
polynomials, we consider in this section the graph $K_{n,\mathbf{r
_{p}}=K_{r_{1},\ldots ,r_{p}}\cup O_{n}$ of
order $n+r_{1}+\cdots +r_{p}$ and chromatic numbe
\begin{equation*}
\chi \left( K_{n,\mathbf{r}_{p}}\right) =\max \left( \chi \left(
K_{r_{1},\ldots ,r_{p}}\right) ,\chi \left( O_{n}\right) \right) =\max
\left( p,1\right)=p .
\end{equation*
Similarly to the numbers $\left( r_{1},\ldots ,r_{p}\right) $-Stirling
numbers, we present here a new
class of the Stirling numbers having also triangular recurrence relation
\newline
To start, let us giving the following definition.

\begin{definition}
Let $R_{1},\ldots ,R_{p}$ be subsets of the set $\left[ n\right] $ with
\left \vert R_{i}\right \vert =r_{i}$ and $R_{i}\cap R_{j}=\varnothing $ for
all $i,j=1,\ldots ,p,$ $i\neq j.$ The $K\left( r_{1},\ldots ,r_{p}\right)
-Stirling number of the second kind, denoted by $
{n \brace k}_{\!\!K\left( r_{1},\ldots ,r_{p}\right) },$ counts the number of
partitions of the set $\left[ n\right] $ into $k$ non-empty subsets such
that if $x\in R_{i}$ and $y\in R_{j}$ with $i\neq j,$ then each subset do
not containing simultaneously $x$ and $y$.
\end{definition}

\noindent From this definition, we may state the following
\begin{eqnarray*}
{n \brace k}_{\!\!K\left( r_{1},\ldots ,r_{p}\right) } &=&0,\text{ \  \  \
}n<r_{1}+\cdots +r_{p}\text{ or }k<p, \\
{n \brace k}_{\!\!K\left( r_{1},\ldots ,r_{p}\right) } &=&
{n \brace k}_{\!\!r_{1}+\ldots +r_{p}}\text{ \ if }r_{1},\ldots ,r_{p}\in \left \{
0,1\right \} , \\
{n \brace k}_{\!\!K\left( r_{1},\ldots ,r_{p}\right) } &=&
{n \brace k}_{\!\!K\left( r_{\sigma \left( 1\right) },\ldots ,r_{\sigma \left(
p\right) }\right) }\text{ for all permutations }\sigma \text{ on the set
\left \{ 1,\ldots ,p\right \} .
\end{eqnarray*
Furthermore, for $r_{1}\leq r_{2}\leq \cdots \leq r_{p},$ the coefficient
\alpha _{k}\left( K_{n,\mathbf{r}_{p}}\right) $ represents the number of
ways of partitioning the set  of
n+\left \vert \mathbf{r}_{p}\right \vert $ vertices of $K_{n,\mathbf{r
_{p}} $ into $k$ independent sets, and by the definition of the graph $K_{n
\mathbf{r}_{p}},$ any two elements $x$ of the $i$-th block of the subgraph
K_{r_{1},\ldots ,r_{p}}$ and $y$ of the $j$-th block of $K_{r_{1},\ldots
,r_{p}},$ with $i\neq j,$ can't be in the same subset. So, this number is
exactly ${n+r_{1}+\cdots
+r_{p} \brace k}_{\!\!K\left( r_{1},\ldots ,r_{p}\right) }.\ $Then, we may state tha
\begin{equation*}
\alpha _{k}\left( K_{n,\mathbf{r}_{p}}\right) ={n+r_{1}+\cdots
+r_{p} \brace k}_{\!\!K\left( r_{1},\ldots ,r_{p}\right) }:={n+\left
\vert \mathbf{r}_{p}\right \vert  \brace k}_{\!\!K\left( \mathbf{r}_{p}\right) }.
\end{equation*

\noindent Set $ H:=K_{r_{1},\ldots ,r_{p}}$ in Theorem \ref{TT1} to obtain:

\begin{corollary}
\label{C2}Let $n,\ k$ be integers such that $p\leq k\leq n$ and $n\geq
\left \vert \mathbf{r}_{p}\right \vert .$ We hav
\begin{equation*}
{n \brace k}_{\!\!K\left( \mathbf{r}_{p}\right) }=\sum_{j=p}^{k}{s+j \brace k}_{\!\!j}{n-s \brace j}_{\!\!K\left( \mathbf{r}_{p}\right) },\  \
0\leq s\leq n-\left \vert \mathbf{r}_{p}\right \vert .
\end{equation*
For $s=1$ we conclude that these numbers obey to the recurrence relation
\begin{equation*}
{n \brace k}_{\!\!K\left( \mathbf{r}_{p}\right) }=k
{n-1 \brace k}_{\!\!K\left( \mathbf{r}_{p}\right) }+{n-1 \brace k-1}_{\!\!K\left(
\mathbf{r}_{p}\right) },\  \  \ n\geq k\geq 1.
\end{equation*
and for $s=n-\left \vert \mathbf{r}_{p}\right \vert $ in Corollary \ref{C2},
we ge
\begin{equation*}
{n \brace k}_{\!\!K\left( \mathbf{r}_{p}\right) }=\sum_{j=p}^{k}
{n-\left \vert \mathbf{r}_{p}\right \vert +j \brace k}_{\!\!j}
{\left \vert \mathbf{r}_{p}\right \vert  \brace j}_{\!\!K\left( \mathbf{r}_{p}\right) }.
\end{equation*}
\end{corollary}

\noindent In particular, for $s=r_{1}=\cdots =r_{p}=1$ in Corollary \ref{C2}
we obtain the known identity
\begin{equation*}
{n \brace k}_{\!\!p}=k{n-1 \brace k}_{\!\!p}+{n-1 \brace k-1}_{\!\!p}.
\end{equation*}

\begin{remark}
By setting $s=1\ $and $k=p$ or $p+1$ in Corollary \ref{C2} and upon using the identities given in \cite{zou} (see also \cite[Lemma 4.4.2]{don}) by
\begin{eqnarray*}
{\left \vert \mathbf{r}_{p}\right\vert \brace p}_{\!\!K\left( \mathbf{r}_{p}\right) }=1, \ \ \ {\left
\vert \mathbf{r}_{p}\right \vert \brace p+1}_{\!\!K\left( \mathbf{r
_{p}\right) }=\sum_{j=1}^{p}2^{r_{j}-1}-p,
\end{eqnarray*}
we obtain the two initial values of these number
\begin{eqnarray*}
{n \brace p}_{\!\!K\left( \mathbf{r}_{p}\right) } &=&
{n-\left \vert \mathbf{r}_{p}\right \vert +p \brace p}_{\!\!p}, \\
{n \brace p+1}_{\!\!K\left( \mathbf{r}_{p}\right) } &=&
{n-\left \vert \mathbf{r}_{p}\right \vert +p \brace p+1}_{\!\!p}+
{n-\left \vert \mathbf{r}_{p}\right \vert +p+1 \brace p+1}_{\!\!p+1}\left(
\sum_{j=1}^{p}2^{r_{j}-1}-p\right) .
\end{eqnarray*}
\end{remark}

\begin{proposition}
Let
\begin{equation*}
B\left( \lambda ;K_{n,\mathbf{r}_{p}}\right) :=\underset{k\geq 0}{\sum
{n+\left \vert \mathbf{r}_{p}\right \vert  \brace k}_{\!\!K\left(
\mathbf{r}_{p}\right) }\lambda ^{k}.
\end{equation*
Then, we hav
\begin{eqnarray*}
 B\left( \lambda ;K_{n,\mathbf{r}_{p}}\right)
&=&\lambda\exp \left( -\lambda \right) \frac{d}{d\lambda }\left( \exp \left( \lambda \right) B\left(
\lambda ;K_{n-1,\mathbf{r}_{p}}\right) \right) ,\  \ n\geq 1,\\
B\left( \lambda ;K_{0,\mathbf{r}_{p}}\right) &=&B_{r_{1}}\left( \lambda
\right) \cdots B_{r_{p}}\left( \lambda \right),
\end{eqnarray*
where $B_{n}\left(\lambda \right) =\sum_{j=0}^{n} {n \brace j}\lambda ^{j}$ is the classical Bell polynomial.
\end{proposition}

\begin{proof}
From Corollary \ref{C2} we hav
\begin{equation*}
\underset{k\geq 1}{\sum }{n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!K\left( \mathbf{r}_{p}\right)
}\lambda^{k}=\underset{k\geq 1}{\sum }k{n+\left \vert \mathbf{r}_{p}\right \vert-1 \brace k}_{\!\!K\left( \mathbf{r
_{p}\right) }\lambda^{k}+\underset{k\geq 1}{\sum }{n+\left \vert \mathbf{r}_{p}\right \vert-1 \brace k-1}_{\!\!K\left(
\mathbf{r}_{p}\right) }\lambda^{k}
\end{equation*
which gives the first identity. The second one follows from \cite[Lemma 4.4.1]{don}.
\end{proof}

\begin{corollary}
For $n\geq \left \vert \mathbf{r}_{p}\right \vert ,$ the numbers
${n \brace k}_{\!\!K\left( \mathbf{r}_{p}\right) },\ k=p,p+1,\ldots ,n,$ are
log-concave.
\end{corollary}

\begin{proof}
Apply Role's Theorem on the function $f_{n}\left( \lambda \right) :=\exp
\left( \lambda \right) B\left( \lambda ;K_{n,\mathbf{r}_{p}}\right) $ and
use the fact that the roots of the polynomial $B_{n}\left( \lambda \right) $
are real non-positive (see for example \cite {wan}) to conclude (by induction on $n$) that the polynomial
B\left( \lambda ;K_{n,\mathbf{r}_{p}}\right) $ has only real non-positive
roots. After that, apply Newton's inequality \cite[p. 52]{har} to complete
the proof.
\end{proof}
\noindent Similarly to Corollary \ref{C3}, Proposition \ref{P1} states that we have
\begin{corollary}
For $0\leq k\leq n$ let
\begin{equation*}
W\left( n,k\right) ={n+\left \vert \mathbf{r}_{p}\right \vert
\brace k+\left \vert \mathbf{r}_{p}\right \vert }_{\!\!K\left( \mathbf{r}_{p}\right) }.
\end{equation*
Then, the sequence $\left( W\left( n,k\right) ,\ 0\leq k\leq n\right) $ is a P\'{o}lya frequency sequence.
\end{corollary}

\section{Application to the graph $T_{r_{_{1}}}\cup \cdots \cup T_{r_{_{p}}}\cup O_{n}$}

Let $p\geq 1$ be an integer. To give a third application of the chromatic
polynomials, we consider in this section the graph $T_{n,\mathbf{r
_{p}}=T_{r_{_{1}}}\cup \cdots \cup T_{r_{_{p}}}\cup O_{n}.$ Her
\begin{equation*}
\chi \left( T_{n,\mathbf{r}_{p}}\right) =\max \left( \chi \left(
T_{r_{_{1}}}\right) ,\ldots ,\chi \left( T_{r_{_{p}}}\right) ,\chi \left(
O_{n}\right) \right) =\min \left( r_{1}\cdots r_{p},2\right) .
\end{equation*
Similarly to the above numbers, we present here a new class of the Stirling
numbers having also triangular recurrence relation. Then, we may start by
giving the following definition.

\begin{definition}
Let $R_{1},\ldots ,R_{p}$ be subsets of the set $\left[ n\right] $ with
\left \vert R_{i}\right \vert =r_{i}$ and $R_{i}\cap R_{j}=\varnothing $ for
all $i,j=1,\ldots ,p,$ $i\neq j.$ The $T\left( r_{1},\ldots ,r_{p}\right)
-Stirling number of the second kind, denoted by $
{n \brace k}_{\!\!T\left( r_{1},\ldots ,r_{p}\right) },$ counts the number of
partitions of the set $\left[ n\right] $ into $k$ non-empty subsets such
that the minimum element of $R_{i}$ can't be in the same subset with any
other element of $R_{i},$ $i=1,\ldots ,p.$
\end{definition}

\noindent From this definition, we may state the following
\begin{eqnarray*}
{n \brace k}_{\!\!T\left( r_{1},\ldots ,r_{p}\right) } &=&0,\text{ \  \  \
}n<r_{1}+\cdots +r_{p}\text{ or }k<\min \left( r_{1}\cdots r_{p},2\right) ,
\\
{n \brace k}_{\!\!T\left( r_{1},\ldots ,r_{p}\right) } &=&
{n \brace k}\text{ \ if }r_{1},\ldots ,r_{p}\in \left \{ 0,1\right \} , \\
{n \brace k}_{\!\!T\left( r_{1},\ldots ,r_{p}\right) } &=&
{n \brace k}_{\!\!K\left( r_{\sigma \left( 1\right) },\ldots ,r_{\sigma \left(
p\right) }\right) }\text{ for all permutations }\sigma \text{ on the set
\left \{ 1,\ldots ,p\right \} .
\end{eqnarray*
Furthermore, since the trees of same order have the same chromatic polynomial,
it is sufficient to apply such results on the star graphs. So, choice $T_{i}$
to be the star graph with $r_{i}$ vertices ($T_{i}=K_{1,r_{i}-1}$ ) with the universal vertex is the minimum element of
R_{i}$ (the set-vertex of $T_{i}$). Hence, the coefficient $\alpha _{k}\left(
T_{n,\mathbf{r}_{p}}\right) $ represents the number of ways of partitioning
the set of $n+\left \vert \mathbf{r}_{p}\right \vert $ vertices of $T_{n,\mathbf{r}_{p}}$ into $k$
independent sets, and by the definition of the graph $T_{n,\mathbf{r}_{p}},$
the universal vertex of $T_{r_{i}}$ ($i=1,\ldots ,p$) can't be in the same
independent set with any other vertex of $T_{r_{i}}$. So, this number is
exactly ${n+r_{1}+\cdots+r_{p} \brace k}_{\!\!T\left( r_{1},\ldots ,r_{p}\right) }.$ \\ Then, we may state tha
\begin{equation*}
\alpha _{k}\left( T_{n,\mathbf{r}_{p}}\right) ={n+r_{1}+\cdots
+r_{p} \brace k}_{\!\!T\left( r_{1},\ldots ,r_{p}\right) }:={n+\left
\vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!T\left( \mathbf{r}_{p}\right) }.
\end{equation*

\begin{theorem}
For $1\leq r_{1}\leq r_{2}\leq \cdots \leq r_{p},$ the polynomial
z^{n+p}\left( z-1\right) ^{\left \vert \mathbf{r}_{p}\right \vert -p}$ can
be written in the basis $\left \{ \left( \lambda \right) _{k},\ k=0,1,\ldots,n+\left \vert \mathbf{r}_{p}\right \vert
\right \} $ as follows
\begin{equation*}
z^{n+p}\left( z-1\right) ^{\left \vert \mathbf{r}_{p}\right \vert
-p}=\sum_{k=0}^{n+\left \vert \mathbf{r}_{p}\right \vert }
{n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!T\left( \mathbf{r
_{p}\right) }\left( z\right) _{k},
\end{equation*
Furthermore, we hav
\begin{equation*}
{n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!T\left(
\mathbf{r}_{p}\right) }={n+\left \vert \mathbf{r}_{p}\right
\vert \brace k}_{\!\!T\left( \left \vert \mathbf{r}_{p}\right \vert -p+1\right) }.
\end{equation*}
\end{theorem}

\begin{proof}
From \cite[Sec. 1.2]{don}, the chromatic polynomial of the graph $T_{n
\mathbf{r}_{p}}$ is
\begin{eqnarray*}
P\left( T_{n,\mathbf{r}_{p}},z\right) &=&P\left( O_{n},z\right) P\left(
T_{r_{1}},z\right) \cdots P\left( T_{r_{p}},z\right) \\&=&\left( z\right)
^{n+p}\left( z-1\right) ^{\left \vert \mathbf{r}_{p}\right \vert
-p}\\&=&\sum_{i=0}^{n+\left \vert \mathbf{r}_{p}\right \vert }
{n+\left \vert \mathbf{r}_{p}\right \vert \brace i}_{\!\!T\left( \mathbf{r
_{p}\right) }\left( z\right) _{i}.
\end{eqnarray*
This gives the first result. In particular, we have
\begin{equation*}
z^{n+1}\left(z-1\right) ^{r-1}=\sum_{k=0}^{n}{n+r \brace k}_{\!\!T\left( r\right)
}\left( z\right) _{k},
\end{equation*
and by combining this result with the first one we obtain the second result.
\end{proof}

\noindent For the choice $H:=T_{r_{_{1}}}\cup \cdots \cup
T_{r_{_{p}}}$ in Theorem \ref{TT1} we get:

\begin{corollary}
Let $n,\ k$ be integers such that $2\leq k\leq n$ and $n\geq \left \vert
\mathbf{r}_{p}\right \vert .$ Then, for $1\leq r_{1}\leq \cdots
\leq r_{p},$ we hav
\begin{equation*}
{n \brace k}_{\!\!T\left( \mathbf{r}_{p}\right)
}=\sum_{j=2}^{k}{s+j \brace k}_{\!\!j}
{n-s \brace j}_{\!\!T\left( \mathbf{r}_{p}\right) },\  \ 0\leq s\leq n-\left \vert
\mathbf{r}_{p}\right \vert .
\end{equation*
In particular, for $s=1,$ these numbers obey to the following triangular
recurrence relation
\begin{equation*}
{n \brace k}_{\!\!T\left( \mathbf{r}_{p}\right) }=
{n-1 \brace k-1}_{\!\!T\left( \mathbf{r}_{p}\right) }+k{n-1 \brace k}_{\!\!T\left(
\mathbf{r}_{p}\right) },
\end{equation*
and for $s=n-\left \vert \mathbf{r}_{p}\right \vert $ we ge
\begin{equation*}
{n \brace k}_{\!\!T\left( \mathbf{r}_{p}\right)
}=\sum_{j=2}^{k}{n-\left \vert \mathbf{r}_{p}\right \vert
+j \brace k}_{\!\!j}{\left \vert \mathbf{r}_{p}\right \vert
\brace j}_{\!\!T\left( \mathbf{r}_{p}\right) }.
\end{equation*}
\end{corollary}

\noindent For $i=1,\ldots ,p,$ set $\ s=1$ and $T_{r_{i}}\cup H:=T_{r_{_{1}}}\cup
\cdots \cup T_{r_{_{p}}}\cup O_{n}\ $ in Theorem \ref{TT3} to obtain:

\begin{corollary}
Let $n,\ k$ be integers such that $1\leq k\leq n$ and $n\geq \left \vert
\mathbf{r}_{p}\right \vert .$ Then, for a fixed $i$ \ $(i=1,\ldots ,p)$ with $r_{i}\geq 1,$  we hav
\begin{equation*}
{n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!T\left(
\mathbf{r}_{p}\right) } =\left( k-1\right){n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!T\left(
\mathbf{r}_{p}-\mathbf{e}_{i}\right) }+{n+\left \vert \mathbf{r}_{p}\right \vert \brace k-1}_{\!\!T\left(
\mathbf{r}_{p}-\mathbf{e}_{i}\right) } ,\  \ n\geq 1.
\end{equation*}
\end{corollary}

\begin{remark}
For $1\leq r_{1}\leq \cdots \leq r_{p},$ then similarly to Remark
\ref{R1}, we ge
\begin{equation*}
{n+\left \vert \mathbf{r}_{p}\right \vert \brace k}_{\!\!T\left(
\mathbf{r}_{p}\right) }=\frac{1}{k!}\sum_{j=0}^{k}\left( -1\right) ^{k-j
\binom{k}{j}j^{n+p}\left( j-1\right)
^{\left \vert \mathbf{r}_{p}\right \vert -p},
\end{equation*
and since $\chi \left( T_{n,\mathbf{r}_{p}}\right) =2,$ when $r_{1}\cdots r_{p}>1$, then the two initial values of these numbers are given by
\begin{equation*}
{n+\left \vert \mathbf{r}_{p}\right \vert \brace 2}_{\!\!T\left(
\mathbf{r}_{p}\right) }=2 ^{n+p-1}, \ \ {n+\left \vert \mathbf{r}_{p}\right \vert \brace 3}_{\!\!T\left(
\mathbf{r}_{p}\right) }=3 ^{n+p-1}\times2^{\left \vert \mathbf{r}_{p}\right \vert -p-1}-2 ^{n+p-1},
\end{equation*}
and for $r_{1}=\cdots =r_{p}=1$ we get
\begin{equation*}
{n \brace k}_{\!\!T\left(\mathbf{r}_{p}\right) }={n \brace k}.
\end{equation*
\end{remark}
\noindent Similarly to Corollary \ref{C3}, we have:
\begin{corollary}
For $0\leq k\leq n$ let
\begin{equation*}
V\left( n,k\right) ={n+\left \vert \mathbf{r}_{p}\right \vert
\brace k+\left \vert \mathbf{r}_{p}\right \vert }_{\!\!T\left( \mathbf{r}_{p}\right)
}.
\end{equation*
Then, the sequence $\left( V\left( n,k\right) ,\ 0\leq k\leq n\right) $ is
log-concave and P\'{o}lya frequency sequence, and, the sequence of
polynomials $\left( V_{n}\left( q\right) ,n\geq 0\right) $ defined by
\begin{equation*}
V_{n}\left( q\right) =\overset{n}{\underset{k=0}{\sum }}{n+\left \vert \mathbf{r}_{p}\right \vert
\brace k+\left \vert \mathbf{r}_{p}\right \vert }_{\!\!T\left( \mathbf{r}_{p}\right)
}q^{k}
\end{equation*
is a $q$-log-convex sequence.
\end{corollary}


On the densest packing of polycylinders in any dimension

Abstract: Using transversality and a dimension reduction argument, a result of A. Bezdek and W. Kuperberg is applied to polycylinders $\mathbb{D}^2\times \mathbb{R}^n$, showing that the optimal packing density is $\pi/\sqrt{12}$ in any dimension.

\section{Definitions to Include}

\section{Introduction}
G. Fejes T\`oth and W. Kuperberg  \cite{toth1993blichfeldt} describe a method for computing an upper bound for packings by infinite polycylinders  -- objects isometric to $\mathbb{D}^2 \times \mathbb{R}^n$ in $\mathbb{R}^{n+2}$.  This article explicitly computes their bound and then proves the sharp upper bound of $\pi/\sqrt{12}$ for the packing density of infinite polycylinders in any dimension, generalizing a result of A. Bezdek and W. Kuperberg \cite{bezdek1990maximum}.

\subsection{Objects considered}
Open and closed Euclidean unit $n$-balls will be denoted $\mathbb{B}^n$ and $\mathbb{D}^n$ respectively. The closed unit interval is denoted $\mathbb{I}$. A general polycylinder $C$ is a set isometric to $\Pi_{i=1}^{i=m}\lambda_i\mathbb{D}^{k_i}$ in $\mathbb{R}^{ k_1+\dots + k_m}$, where $\lambda_i$ is in $[0,\infty]$.  For this article, the term polycylinder refers to the special case of an infinite polycylinder over a two-dimensional disk of unit radius.

\begin{definition}
A \defv{polycylinder} is a set isometric to\, $\mathbb{D}^2 \times \mathbb{R}^n $ in $ \mathbb{R}^{n+2}$.
\end{definition}

The following are standard definitions.

\begin{definition}
A \defv{polycylinder packing of\, $\mathbb{R}^{n+2}$} is a countable family $$\mathscr{C} = \{C_i\}_{i \in I}$$ of polycylinders $C_i \subset \mathbb{R}^{n+2}$ with mutually disjoint interiors.
\end{definition}

\begin{definition}
The \defv{upper density} $\delta^+ (\mathscr{C})$ of a packing $\mathscr{C}$ of\, $\mathbb{R}^n$ is defined to be
 \[ \delta^+ (\mathscr{C}) = \limsup_{r\rightarrow \infty}\frac{ \Vol(\mathscr{C}\cap r\mathbb{B}^n)}{\Vol(r\mathbb{B}^n)}.\]
  Note that this notion can be generalized further by replacing $\mathbb{B}^n$ with an arbitrary convex body $K.$
 \end{definition}
 \begin{definition}
 The \defv{upper packing density}\, $\delta^+ (C)$ of an object\, $C$ is the supremum of $\delta^+ (\mathscr{C})$ over all packings $\mathscr{C}$ of\, $\mathbb{R}^n$ by \,$C$.
\end{definition}

The definition of density is equivalent to a number of other definitions under some mild assumptions.  For more on density, see for example \cite{brass2005research, radin2004orbits, conway1999recent, federer1969geometric}.

\section{Computing a bound of Fejes T\'oth and Kuperberg }\label{comp}
 In \cite{toth1993blichfeldt}, Fejes T\'oth and Kuperberg describe a method for computing upper bounds on a broad class of objects in any dimension: Blichfeldt-type results for balls \cite{blichfeldt1929minimum, rankin1947closest} extend to results for outer parallel bodies.  Fejes T\'oth and Kuperberg compute upper bounds for the density of packings by cylinders $\mathbb{D}^{n-1} \times t\mathbb{I}$ and also for packings by outer parallel bodies of line segments, i.e. capped cylinders, in $\mathbb{R}^n$ but do not address the case of polycylinders. This section recalls that work and explicitly computes the \cite{toth1993blichfeldt} bound to be $\delta^+(\mathscr{C}) \le 0.941533\dots$ for polycylinder packings in any dimension.

\subsection{Background}
Given a packing of $\mathbb{R}^n$ by congruent objects $\mathscr{C} = \{C_i\}_{i\in I}$, there are a fixed body $C \subset \mathbb{R}^n$ and isometries $\{\phi_i\}_{i\in I}$ of $\mathbb{R}^n$ such that $C_i = \phi_iC$ for all $i$ in $I$.

\begin{definition}
 A function $f:\mathbb{R}^{n} \rightarrow \mathbb{R}^+$ is a \defv{Blichfeldt gauge} for a convex body $C\subset \mathbb{R}^n$ if for any collection of isometries $\Phi=\{\phi_i\}_{i\in I}$ of $\mathbb{R}^n$ where $\mathscr{C} = \{\phi_iC\}_{i\in I}$ is a packing and for all $x$ in $\mathbb{R}^n,$
 $$\sigma_\Phi(f)(x) := \sum_{i\in I} f(\phi_i^{-1}x) \le 1.$$
\end{definition}

Notice that the characteristic function $\bold{1}_C$ of $C$ is a Blichfeldt gauge for $C$. Replacing $\bold{1}_C$ with a more general Blichfeldt gauge $f$ lets one replace the characteristic function of the packing $\bold{1}_{\mathscr{C}}$ with a diffuse version $\sigma_\Phi(f)$. This new function $\sigma_\Phi(f)$ has the same general characteristics as $\bold{1}_{\mathscr{C}}$, is still bounded pointwise by $1$ in the ambient space and is uniformly bounded independent of $\Phi$ in the moduli space of packings. As $f$ may have greater mass than $\bold{1}_C$, this allows one to estimate the volume of the interstices of a packing and thereby bound the packing density.

\begin{example}\label{exmp} Blichfeldt initially uses the radial function $2f_0$ where

\begin{displaymath}
   f_0(r) = \left\{
     \begin{array}{ll}
       \frac{1}{2}(2-r^2) & : 0 \le r \le \sqrt{2} \\
       0 & : r > \sqrt{2}
     \end{array}
   \right.
\end{displaymath}
and showed that $f_0$ is a Blichfeldt gauge. Then, for a packing $\mathscr{C} =\{\phi_i C\}_{i\in I}$ of a cube $t\mathbb{I}^n$ by spheres, the support of $\sigma_\Phi(f)$ is contained in a slightly larger cube $(t+2\sqrt{2}-2)\mathbb{I}^n$. A bound on sphere packing density can then be extracted as follows. From the definition of the Blichfeldt gauge and integrating in spherical coordinates, one finds

$$(t+2\sqrt{2}-2)^n \ge \left\vert I \right\vert\int_{\mathbb{R}^n}\hspace*{-7pt} f_0\, \mathrm{dV} = \frac{\left\vert I \right\vert\Vol(\mathbb{B}^n)2^{\frac{n+2}{2}}}{n+2}.$$

When density is measured relative to a cube,

$$\delta^+_{}(\mathscr{C}) = \frac{\left\vert I \right\vert \Vol(\mathbb{B}^n)}{t^n} \le \frac{n+2}{2^{\frac{n+2}{2}}}\left(1+\frac{2\sqrt{2}-2}{t}\right)^n \hspace*{-7pt}.$$

It is easy to see that the same method works when $t\mathbb{I}^n$ is replaced with $\mathbb{B}^n_{t/2}$.  By passing to the limit, the bound $$\delta^+(\mathscr{C}) \le  \frac{n+2}{2^{\frac{n+2}{2}}}$$ holds for any sphere packing in $\mathbb{R}^n.$
\end{example}

\subsection{Blichfeldt-type bound for polycylinders}

Example \ref{exmp} motivates the following general observations.

\begin{theorem}[Blichfeldt]\label{blich} If $g$ is a Blichfeldt gauge for a body $C$, then $\delta^+(\mathscr{C})\le \Vol(C)/J(g)$ where
$$J(g) = \int_{\mathbb{R}^n}\hspace*{-7pt}g\, \mathrm{dV}.$$
\end{theorem}

\begin{theorem}[Fejes T\'oth--Kuperberg]
If $f(\alpha)$, $\alpha \ge 0$, is a real valued function such that $f(|x|)$ is a Blichfeldt gauge for the unit ball, and $C$ is a convex body with inradius $r(C)$, then for any $\varrho \le r(C)$
$$g(x) = f\left(\frac{d(x,C_{-\varrho})}{\varrho}\right)$$
is a Blichfeldt gauge for $C$, where $C_{-\varrho}$ is the inner parallel body of $C$ at distance $\varrho$.
\end{theorem}

For their more general results, Fejes T\'oth and Kuperberg do not use $f_0$, but rather Blichfeldt's modified version.
\begin{definition} The \defv{modified Blichfeldt gauge} \cite{blichfeldt1929minimum} for\, $\mathbb{D}^n$ is the radial function
\begin{displaymath}
   f_1(r) = \left\{
     \begin{array}{ll}
       1 & : 0 \le r \le 2-\sqrt{2} \\
       \frac{1}{2}(2-r)^2 & : 2 - \sqrt{2}\le r \le 1\\
        \frac{1}{2}(2-r^2) & : 1\le r\le \sqrt{2}\\
       0 & : r > \sqrt{2}
     \end{array}
   \right.
\end{displaymath}
\end{definition}

\begin{definition}[Fejes T\'oth--Kuperberg]For the two-dimensional gauge $f_1$ defined above, $$A_2 := J(f_1)/\Vol(\mathbb{D}^2) = (29-16\sqrt{2})/6.$$
\end{definition}

From the previous theorems and definitions, the results of \cite{toth1993blichfeldt} give a rough estimate for the density of infinite polycylinders as follows.  Consider $C(t) = \mathbb{D}^{n+2} + t\mathbb{I}^n$ in $\mathbb{R}^{n+2}$ and the gauge $g_t(x) =f_1(d(x,C(t)_{-1}))$, where $f_1$ is the modified Blichfeldt gauge and $C(t)_{-1}$ is the inner parallel body at distance $1$, i.e. an $n$-cube of height $t$. From Theorem \ref{blich} an estimate of the integral$ \int_{\mathbb{R}^{n+2}}g_t \, \mathrm{dV}$ gives a density bound. By integrating $g_t$ over $C(t)_{-1} \times \mathbb{R}^2$ and noticing that contribution from the complement $\mathbb{R}^{n+2} \smallsetminus (C(t)_{-1} \times \mathbb{R}^2)$ are of strictly lower order -- it is bounded above by a constant times the $(n-1)$-Hausdorff measure of the boundary $\partial C(t)_{-1} \subset C(t)_{-1}$, it follows that
 $$\delta^+(C(t)) \le \frac{\pi t^n}{\pi A_2 t^n + O(t^{n-1})}.$$
In the limit as $t$ goes to infinity, this gives a bound of $1/A_2 = .941533...$ for infinite polycylinders in any dimension.

\begin{remark}
The Blichfeldt gauge method yields constrained optimization problems of the form
$$\textrm{maximize } J(g)
\textrm{ subject to }\sigma_\Phi(g) \le 1 \textrm{ for all }\Phi \textrm{ where } \Phi(C) \textrm{ is a packing.}$$
\end{remark}

\section{Transversality}\label{trans}
This section introduces the required transversality arguments in affine geometry.
\begin{definition}
A \defv{d-flat} is a d-dimensional affine subspace of\, $\mathbb{R}^n$.
\end{definition}

\begin{definition}
The \defv{parallel dimension} $dim_\parallel \{F,\dots, G\}$ of a collection of flats $\{F, \dots, G\}$ is the dimension of their maximal parallel sub-flats.
\end{definition}

The notion of parallel dimension can be interpreted in several ways, allowing a modest abuse of notation.

For a collection of flats $\{F, \dots, G\}$, consider their tangent cones at infinity $\{F_\infty, \dots, G_\infty\}$.  The parallel dimension of $\{F, \dots, G\}$ is the dimension of the intersection of these tangent cones.  This may be viewed as the limit of a rescaling process $\mathbb{R}^n \rightarrow  r\mathbb{R}^n$ as $r$ tends to $0$, leaving only the scale-invariant information.

For a collection of flats $\{F, \dots, G\}$, consider each flat as a system of linear equations. The corresponding homogeneous equations determine a collection of linear subspaces $\{F_\infty, \dots, G_\infty\}$. The parallel dimension is the dimension of their intersection $F_\infty \,\cap \dots \cap\, G_\infty$.

\begin{definition}
Two disjoint $d$-flats are \defv{parallel} if their parallel dimension is $d$, that is, if every line in one is parallel to a line in the other.
\end{definition}

\begin{definition}
Two disjoint d-flats are \defv{skew} if their parallel dimension is less than $d.$
\end{definition}

\begin{figure}[t]
\centering
   \begin{subfigure}{0.4\linewidth} \centering
     \includegraphics[scale=0.5]{Fig1a.pdf}
     \subcaption*{(a)}
   \end{subfigure}
   \begin{subfigure}{0.4\linewidth} \centering
     \includegraphics[scale=0.5]{Fig1b.pdf}
      \subcaption*{(b)}
       \end{subfigure}
       \caption{\newline Disjoint 2-flats in $\mathbb{R}^4$ with (a) parallel dimension 2 and (b) parallel dimension 1.} \label{fig:twofigs}
\end{figure}

\begin{lemma}
A pair of disjoint $n$-flats in\, $\mathbb{R}^{n+k}$ with $n\ge k$ has parallel dimension strictly greater than $n-k.$
\end{lemma}

\begin{proof}
 By homogeneity of $\mathbb{R}^{n+k}$, let $F=F_\infty.$  As $F_\infty$ and $G$ are disjoint, $G$ contains a non-trivial vector $\bold{v}$ such that $G = G_\infty +\bold{v}$ and $\bold{v}$ is not in $F_\infty + G_\infty.$  It follows that
$$dim (\mathbb{R}^{n+k}) \ge dim(F_\infty + G_\infty + \spanv(\bold{v}) )> dim(F_\infty + G_\infty) $$
$$= dim (F_\infty) + dim (G_\infty) - dim (F_\infty \cap G_\infty).$$
Count dimensions to find $n+k > n +n - dim_\parallel (F_\infty,G_\infty).$
\end{proof}
\begin{corollary}\label{dimcor}
A pair of disjoint $n$-flats in\, $\mathbb{R}^{n+2}$ has parallel dimension at least $n-1$.
\end{corollary}

\section{Results}

\subsection{Pairwise foliations and dimension reduction}

\begin{figure}[t]
\centering
   \begin{subfigure}{0.3\linewidth} \centering
     \includegraphics[scale=0.3]{Fig2.pdf}
     \subcaption*{$t=-\epsilon$}
   \end{subfigure}
   \begin{subfigure}{0.3\linewidth} \centering
     \includegraphics[scale=0.3]{Fig2.pdf}
      \subcaption*{$t=0$}
       \end{subfigure}
        \begin{subfigure}{0.3\linewidth} \centering
     \includegraphics[scale=0.3]{Fig2.pdf}
      \subcaption*{$t=\epsilon$}
       \end{subfigure}
\caption{\newline The leaves of a pair of polycylinders in $\mathbb{R}^4$ indexed by the common parallel direction $t$ are identical.} \label{fig:twofigs}
\end{figure}

\begin{definition}
The \defv{core} $a_i$ of a polycylinder $C_i$ isometric to $\mathbb{D}^2~\times~\mathbb{R}^n$ in $\mathbb{R}^{n+2}$
is the distinguished $n$-flat defining $C_i$ as the set of points at most distance $1$ from $a_i$.
\end{definition}

In a packing $\mathscr{C}$ of\, $\mathbb{R}^{n+2}$ by polycylinders, Corollary \ref{dimcor} shows that, for every pair of polycylinders $C_i$ and $C_j$, one can choose parallel $(n-1)$-dimensional subflats $b_i \subset a_i$ and $b_j \subset a_j$ and define a product foliation $$\mathscr{F}^{b_i,b_j}:\mathbb{R}^{n+2} \rightarrow \mathbb{R}^{n-1} \times \mathbb{R}^3$$ with $\mathbb{R}^3$ leaves that are orthogonal to $b_i$ and to $b_j$. Given a point $x$ in $a_i$, there is a distinguished $\mathbb{R}^3$ leaf $F_x^{b_i,b_j}$ that contains the point $x$. The foliation $\mathscr{F}^{b_i,b_j}$ restricts to foliations of $C_i$ and $C_j$ with right-circular-cylinder leaves.

\subsection{The Dirichlet slice}

\begin{definition}
In a packing $\mathscr{C}$ \hspace*{-5pt} of\, $\mathbb{R}^{n+2}$ by polycylinders, the \defv{Dirichlet cell} $D_i$ associated to a polycylinder $C_i$ is the set of points in $\mathbb{R}^{n+2}$ no further from $C_i$ than from any other polycylinder in $\mathscr{C}$.
\end{definition}
The Dirichlet cells of a packing partition $\mathbb{R}^{n+2},$ because $C_i \subset D_i$ for all polycylinders $C_i$.  To bound the density $\delta^+(\mathscr{C})$, it is enough to fix an $i$ in $I$ and consider the density of $C_i$ in $D_i$. For the Dirichlet cell $D_i,$ there is a slicing as follows.
\begin{definition}
Given a fixed a polycylinder $C_i$ in a packing $\mathscr{C}$ \hspace*{-7pt} of\, $\mathbb{R}^{n+2}$ by polycylinders and a point $x$ on the core $a_i$, the plane $p_x$ is the $2$-flat orthogonal to $a_i$ and containing the point $x$. The \defv{Dirichlet slice} $d_x$ is the intersection of $D_i$ and $p_x.$
\end{definition}

Note that $p_x$ is a sub-flat of $F_x^{b_i,b_j}$ for all $j$ in $I.$

\subsection{Bezdek--Kuperberg bound}

For any point $x$ on the core $a_i$ of a polycylinder $C_i$, the results of Bezdek and Kuperberg \cite{bezdek1990maximum} apply to the Dirichlet slice $d_x$.

\begin{lemma}\label{lem1}
A Dirichlet slice is convex and, if bounded, a parabola-sided polygon.
\end{lemma}
\begin{proof}

Construct the Dirichlet slice $d_x$ as an intersection.  Define $d^j$ to be the set of points in $p_x$ no further from $C_i$ than from $C_j$.  Then the Dirichlet slice $d_x$ is realized as
$$d_x = \{\cap_{j\in I} d^j \}.$$

Each arc of the boundary of $d_x$ in $p_x$ is given by an arc of the boundary of some $d^j$ in $p_x.$ The boundary of $d^j$ in $p_x$ is the set of points in $p_x$ equidistant from $C_i$ and $C_j.$ Since the foliation $\mathscr{F}^{b_i,b_j}$  is a product foliation, the arc of the boundary of $d^j$ in $p_x$ is also the set of points in $p_x$ equidistant from the leaf $C_i \cap F_x^{b_i,b_j}$ of $\mathscr{F}^{b_i,b_j}|_{C_i}$ and the leaf $C_j \cap F_x^{b_i,b_j}$ of $\mathscr{F}^{b_i,b_j}|_{C_j}$. This reduces the analysis to the case of a pair cylinders in $\mathbb{R}^3.$ From \cite{bezdek1990maximum}, it follows that $d^j$ is convex and the boundary of $d_j$ in $p_x$ is a parabola; the intersection of such sets $d^j$ in $p_x$ is convex, and a parabola-sided polygon if bounded.
 \end{proof}

 Let $S_x(r)$ be the circle of radius $r$ in $p_x$ centered at $x$.
\begin{lemma} \label{lem2}The vertices of $d_x$ are not closer to $S_x(1)$ than the vertices of a regular hexagon circumscribed about $S_x(1).$
\end{lemma}
\begin{proof}
A vertex of $d_x$ occurs where three or more polycylinders are equidistant, so the vertex is the center of a $(n+1)$-ball $B$ tangent to three polycylinders.  Thus $B$ is tangent to three disjoint unit $(n+2)$-balls $B_1$, $B_2$, $B_3$.  By projecting into the affine hull of the centers of $B_1$, $B_2$, $B_3$, it is immediate that the radius of $B$ is no less than $2/\sqrt{3} -1.$
\end{proof}

\begin{lemma} \label{lem3}
Let $y$ and $z$ be points on the circle $S_x(2/\sqrt{3})$.  If each of $y$ and $z$ is equidistant from $C_i$ and $C_j$, then the angle $yxz$ is smaller than or equal to $2\arccos (\sqrt{3} -1) = 85.8828\dots^\circ.$
\end{lemma}

\begin{proof} Following \cite{bezdek1990maximum, kusner2013upper}, the existence of a supporting hyperplane of $C_i$ that separates $\Int(C_i)$ from $\Int(C_j)$ suffices.
\end{proof}

In \cite{bezdek1990maximum}, it is shown that planar objects satisfying Lemmas \ref{lem1}, \ref{lem2} and \ref{lem3} have area no less than $\sqrt{12}.$ As the bound holds for all Dirichlet slices, it follows that $\delta^+(\mathbb{D}^2\times \mathbb{R}^n) \le \pi/\sqrt{12}$ in $\mathbb{R}^{n+2}.$  The product of the dense disk packing in the plane with $\mathbb{R}^n$ gives a polycylinder packing in $\mathbb{R}^{n+2}$ that achieves this density.  Combining this with the result of Thue \cite{thue1910densest} for $n=0$ and the result of Bezdek and Kuperberg \cite{bezdek1990maximum} for $n=1$, it follows that

\begin{theorem}
$\delta^+(\mathbb{D}^2 \times \mathbb{R}^n) \le \pi/\sqrt{12} $ for all natural numbers $n.$
\end{theorem}

Metric Geometry (math.MG); Combinatorics (math.CO); Number Theory (math.NT)


On Routh-Steiner Theorem and Generalizations

Abstract: Following Coxeter we use barycentric coordinates in affine geometry to prove theorems on ratios of areas. In particular, we prove a version of Routh-Steiner theorem for parallelograms.

\section{Introduction}

Coxeter in his book \cite[p. 211]{Coxeter}, considered the following theorem
of \textit{affine} geometry;

\begin{theorem}
\label{Routh th}If the sides $BC,CA,AB$ of a triangle $ABC$ are divided at
L,M,N$ $\ $in the respective ratios $\lambda :1,$ $\mu :1,$ $\nu :1,$ the
cevians $AL,BM,CN$ form a triangle whose area is
\[
\frac{(\lambda \mu \nu -1)^{2}}{(\lambda \mu +\lambda +1)(\mu \nu +\mu
+1)(\nu \lambda +\nu +1)}
\]

times that of $ABC.$
\end{theorem}

He emphasized that this result was discovered by Steiner, but simultaneously
cited two references the first was\ Steiner's work \cite[p. 163-168]{Steiner}
and the second was Routh's work: \cite[p. 82]{Routh}. Later in his book he
referred to the result as "Routh's theorem" \cite[p. 219]{Coxeter} admitting
to the contribution of both scientists in revealing the theorem.

Coxeter gave a general proof of this result using \textit{barycentric
coordinates attributed to M\"{o}bius\textit{.} These are \textit{homogeneous}
coordinates $(t_{1},t_{2},t_{3})$, where $t_{1},t_{2},t_{3}$ are masses at
the vertices of a triangle of reference $A_{1}A_{2}A_{3}$. In particular
(1,0,0)$ is $A_{1},$ $(0,1,0)$ is $A_{2}$, $(0,0,1)$ is $A_{3}$ and
(t_{1},t_{2},t_{3})$ corresponds to a point $P$ such that the areas of the
triangles $PA_{2}A_{3},PA_{3}A_{1},PA_{1}A_{2}$ are proportional to the
barycentric coordinates $t_{1},t_{2},t_{3}$ of $P$, respectively (see Fig.
\ref{triangle-fig}). If $t_{1}+t_{2}+t_{3}=1$ then the normalized
barycentric coordinates $(t_{1},t_{2},t_{3})$ are called \textit{areal}
coordinates. In this case the areas of the triangles
PA_{2}A_{3},PA_{3}A_{1},PA_{1}A_{2}$ are $t_{1},t_{2},t_{3}$ times the area
of the whole triangle $A_{1}A_{2}A_{3}$, respectively.

Recently, A. B\'{e}nyi and B. \'{C}urgus \cite{Benyi-Curgus}, proved a
version of Theorem \ref{Routh th} and got a unification of the theorems of
Ceva and Menelaus. In fact they unified two expressions of Routh on area
ratios into one from which they derived both theorems of Ceva and Menelaus
as special cases.

Putting aside the search after "nice expressions", Theorem \ref{Routh th}
implies the following general property:

\begin{theorem}
\label{general division of triangle} If the sides of a triangle
A_{2}A_{3}A_{1}$ are divided at $A_{i,1},A_{i,2},...,A_{i,n-1},1\leq i\leq 3$
in the respective ratios $\lambda _{i,1}:\lambda _{i,2}:...:\lambda
_{i,n-1},1\leq i\leq 3,$ then the ratio of the area of any sub polygon to
the area of the whole triangle depends only on $\left\{ \lambda
_{i,j}\right\} ,$ where $1\leq i\leq 3$ and $1\leq j\leq n-1$ (see Fig. \re
{triangular-network} \ for the case $n=5).$
\end{theorem}

A sub polygon of the triangle is defined as a polygon whose vertices are
points of intersections of the cevians $\left\{ A_{i}A_{i,j}\right\} ,1\leq
i\leq 3,1\leq j\leq n-1$.

To prove the theorem we compute first the barycentric coordinates for each
vertex of the sub polygon; these are points of intersections of the
barycentric equations of the cevians $\left\{ A_{i}A_{i,j}\right\} ,1\leq
i\leq 3,1\leq j\leq n-1$. Next we divide the sub polygon into
non-overlapping triangles (in Fig. \ref{triangular-network} the polygon
MNOPQR$ is divided into $4$ triangles). Then we use the method which Coxeter
gave in his book \cite[p. 219]{Coxeter} for proving Theorem \ref{Routh th}
(which will be illustrated in the next section) to conclude that the ratio
of the area of each triangle to the area of the whole triangle depends only
on $\left\{ \lambda _{i,j}\right\} $ and the result follows.

\section{Patterns}

The search after nice expressions of area ratios succeeds in "symmetric"
divisions. In particular we have the following theorem:

\begin{theorem}
\label{1-lampda-1-th}Suppose the sides of a triangle $A_{2}A_{3}A_{1}$ are
divided at $A_{i,1},A_{i,2},1\leq i\leq 3$ in the respective ratios
1:\lambda :1$ (see Fig. \ref{1-lampda-1}). Let $I,J,K,L,M,N$ be the points
of intersection of the corresponding cevians as shown in the table
\[
\begin{array}{cc}
\text{point} & \text{cevians intersection} \\
I & A_{1}A_{1,1}\cap A_{3}A_{3,2} \\
J & A_{2}A_{2,1}\cap A_{3}A_{3,2} \\
K & A_{2}A_{2,1}\cap A_{1}A_{1,2} \\
L & A_{3}A_{3,1}\cap A_{1}A_{1,2} \\
M & A_{2}A_{2,2}\cap A_{3}A_{3,1} \\
N & A_{1}A_{1,1}\cap A_{2}A_{2,2
\end{array
\]

Then the ratio of the area of the hexagon $IJKLMN$ to the area of the
triangle $A_{2}A_{3}A_{1}$ is
\[
\frac{2\lambda ^{2}}{(3+\lambda )(2\lambda +3)}.
\]
\end{theorem}

\begin{proof}
At first note that such a hexagon exists in any triangle since it exists in
an equilateral triangle, by symmetry of the division ratios, and every other
triangle is affine equivalent to an equilateral triangle. Let $G$ be the
centre of gravity of $A_{2}A_{3}A_{1}.$ Since the barycentric coordinates of
$A_{1},A_{2},A_{3}$ are $(1,0,0),(0,1,0),(0,0,1)$ respectively, the
barycentric coordinates of $G$ are $(\frac{1}{3},\frac{1}{3},\frac{1}{3}).$
By symmetry, the area of the hexagon $IJKLMN$ is $6$ times the area of the
triangle $GMN.$ Therefore, it is sufficient to compute the barycentric
coordinates of $M$ and $N.$

In order to find the barycentric equations of the corresponding cevians, we
exhibit the barycentric coordinates of some of the division points. Now,
A_{1,1}$ divides $A_{2}A_{3}$ in the ratio $1:\lambda +1,$ $A_{2,2}$ divides
$A_{3}A_{1}$ in the ratio $\lambda +1:1$ and $A_{3,1}$divides $A_{1}A_{2}$
in the ratio $1:\lambda +1.$ Thus the barycentric coordinates are given in
the table:
\[
\begin{array}{cc}
\text{point} & \text{barycentric coordinates} \\
A_{1,1} & (1,0,\lambda +1) \\
A_{2,2} & (1,\lambda +1,0) \\
A_{3,1} & (0,\lambda +1,1
\end{array
.
\]

We proceed by computing the equations of the cevians whose points of
intersection are $M$ and $N.$ The cevian $A_{2}A_{2,2}$ has the equation
\[
\left\vert
\begin{array}{ccc}
0 & 0 & 1 \\
1 & \lambda +1 & 0 \\
t_{1} & t_{2} & t_{3
\end{array
\right\vert =0
\
which is equivalent to $-(\lambda +1)t_{1}+t_{2}=0.$ The cevian
A_{1}A_{1,1} $ has the equation
\[
\left\vert
\begin{array}{ccc}
0 & 1 & 0 \\
1 & 0 & \lambda +1 \\
t_{1} & t_{2} & t_{3
\end{array
\right\vert =0
\
which is equivalent to $-(\lambda +1)t_{1}+t_{3}=0,$ and the cevian
A_{3}A_{3,1}$ has the equation
\[
\left\vert
\begin{array}{ccc}
1 & 0 & 0 \\
0 & \lambda +1 & 1 \\
t_{1} & t_{2} & t_{3
\end{array
\right\vert =0
\
which is equivalent to $-t_{2}+(\lambda +1)t_{3}=0.$ Hence, the point $M$
can be computed by the following system of equations
\[
\left\{
\begin{array}{c}
-(\lambda +1)t_{1}+t_{2}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ =0\ \ \ \ \  \\
\ \ \ \ \ \ \ \ \ \ \ \ \ -t_{2}+(\lambda +1)t_{3}=
\end{array
\right. .
\]

Substituting $t_{2}=1$ we get,
\[
M=\left( \frac{1}{\lambda +1},1,\frac{1}{\lambda +1}\right) .
\]

Similarly, the point $N$ is obtained from the following equations;

\[
\left\{
\begin{array}{c}
-(\lambda +1)t_{1}\ \ \ \ \ \ \ \ \ \ \ \ \ \ +t_{3}=0 \\
-(\lambda +1)t_{1}+t_{2}\ \ \ \ \ \ \ \ \ \ \ \ \ \ =
\end{array
\right. .
\]

Substituting $t_{1}=1$ we get,
\[
N=\left( 1,\lambda +1,\lambda +1\right) .
\]

Consequently, the ratio of the area of the triangle $GMN$ to the area of the
triangle $A_{2}A_{3}A_{1}$ is equal to the determinant
\[
\left\vert
\begin{array}{ccc}
1/3 & 1/3 & 1/3 \\
\frac{1}{\lambda +1} & 1 & \frac{1}{\lambda +1} \\
1 & \lambda +1 & \lambda +
\end{array
\right\vert =\frac{\lambda ^{2}}{3(\lambda +1)}
\]

divided by the product of the sums of the rows
\[
\left( 1+\frac{2}{\lambda +1}\right) \left( 2\lambda +3\right) .
\]

Therefore, the ratio of the area of the hexagon $IGKLMN$ to the area of the
triangle $A_{2}A_{3}A_{1}$ i
\[
\frac{\frac{6\lambda ^{2}}{3(\lambda +1)}}{\left( 1+\frac{2}{\lambda +1
\right) \left( 2\lambda +3\right) }=\frac{2\lambda ^{2}}{\left( \lambda
+3\right) \left( 2\lambda +3\right) },
\]

in agreement with the statement of the theorem.
\end{proof}

\subsection{Special cases}

\begin{enumerate}
\item If $\lambda =1$ then the ratio of the area of the hexagon to the area
of the triangle is $\frac{1}{10}.$ This special case is referred to as
Marion Walter's theorem \cite{Walter} which states the following: \textit{If
the trisection points of the sides of any triangle are connected to the
opposite vertices, the resulting hexagon has one-tenth the area of the
original triangle. }

\item If $n$ is odd, $n=2k+1,$ then taking $\lambda =\frac{1}{k}$ implies
that the ratio of the area of the hexagon to the area of the triangle is
\[
\frac{2}{(3k+1)(3k+2)}=\frac{8}{9n^{2}-1}.
\
This special case is referred to as Morgan's theorem \cite{Morgan}, which
was proved by T. Watanabe, R. Hanson and F. D. Nowosielski \cite{Watanabe}
using Routh-Steiner theorem several times.

\item If $n$ is even, then taking $\lambda =n-2$ implies a new result which
was not mentioned in the previous discussion. In this case, we have the
following:

\textit{The ratio of the area of the hexagon to the area of the triangle is}
\[
\frac{2(n-2)^{2}}{(n+1)(2n-1)}.
\]
\end{enumerate}

Note that, if each side of the triangle is divided into $n$ equal parts,
then for odd $n$, $n=2k+1,$ the hexagon is the sub polygon in Theorem \re
{general division of triangle}, formed by the cevians $\left\{
A_{i}A_{k,i},A_{i}A_{k+1,i}\right\} ,1\leq i\leq 3.$ While for even $n$, the
hexagon is formed by the cevians $\left\{
A_{i}A_{1,i},A_{i}A_{n-1,i}\right\} ,1\leq i\leq 3.$

\section{Parallelograms}

We can generalize the Routh-Steiner theorem to parallelograms in the
following manner:

\begin{theorem}
If the sides $BC,CD,DA,AB$ of a parallelogram $ABCD$ are divided at $K,L,M,N$
$\ $in the respective ratios $\kappa :1,\lambda :1,$ $\mu :1,$ $\nu :1,$ the
cevians $AK,BL,CM,DN$ form a quadrilateral whose area is
\begin{eqnarray}
&&\frac{1}{2}\frac{(\frac{\mu }{1+\nu }+\frac{1+\mu }{\lambda })(1+\frac
\kappa }{1+\mu }+\frac{\kappa +1}{\nu (1+\mu )})}{(2+\frac{1}{\lambda }
\frac{1}{1+\mu })(1+\mu +\frac{\nu \mu }{1+\nu })(1+2\kappa +\frac{\kappa +
}{\nu })}+  \nonumber \\
&&\frac{1}{2}\frac{(\frac{\kappa +1}{\nu }+\frac{\kappa }{1+\lambda })(1
\frac{1}{\lambda }+\frac{\kappa }{1+\mu })}{(2+\frac{1}{\lambda }-\frac{1}
1+\mu })(1+\kappa +\frac{\lambda \kappa }{1+\lambda })(1+2\kappa +\frac
\kappa +1}{\nu })}  \label{eq1}
\end{eqnarray}

times that of $ABCD.$
\end{theorem}

\begin{proof}
If the barycentric coordinates of $A,B,C$ are $(0,1,0),(0,0,1),(1,0,0)$
respectively, then the barycentric coordinates of $D$ are $(1,1,-1).$ This
is true since the diagonals in the parallelogram bisect each other and the
barycentric coordinates of the midpoint of $AC$ are $(\frac{1}{2},\frac{1}{2
,0)$ (see Fig. \ref{parallel-fig}). The barycentric coordinates of $K,L,M,N$
are shown in the following table:

\[
\begin{array}{cc}
\text{point} & \text{barycentric coordinates} \\
K & (\kappa ,0,1) \\
L & (\lambda +1,\lambda ,-\lambda ) \\
M & (1,\mu +1,-1) \\
N & (0,1,\nu
\end{array
.
\]

Proceeding as in the proof of Theorem \ref{1-lampda-1-th}, the cevians
AK,BL,CM,DN$ have the following equations
\[
\begin{array}{cc}
\text{cevian} & \text{barycentric equation} \\
AK & -t_{1}+\kappa t_{3}=0 \\
BL & -\lambda t_{1}+(\lambda +1)t_{2}=0 \\
CM & t_{2}+(\mu +1)t_{3}=0 \\
DN & (1+\nu )t_{1}-\nu t_{2}+t_{3}=
\end{array
.
\]

Now, let $X,Y,Z,W$ be the points of intersections of pairs of cevians $BL$
and $CM,$ $CM$ and $DN,$ $DN$ and $AK,$ $AK$ and $BL,$ respectively (see
Fig. \ref{parallel-fig}).

Then the barycentric coordinates of these points are given in the following
table
\[
\begin{array}{cc}
\text{point} & \text{barycentric coordinates} \\
X & (\frac{\lambda +1}{\lambda },1,\frac{-1}{\mu +1}) \\
Y & (-\frac{1+\nu +\nu \mu }{1+\nu },-\mu -1,1) \\
Z & (\kappa ,\frac{1+\kappa +\kappa \nu }{\nu },1) \\
W & (\kappa ,\frac{\lambda \kappa }{\lambda +1},1
\end{array
.
\]

The area of the quadrilateral $XYZW$ equals the sum of the areas of the
triangles $XYZ$ and $ZWX.$ Normalizing the barycentric coordinates of
X,Y,Z,W$ and dividing by $2$ which is the area of the parallelogram $ABCD,$
we get that the ratio of the area of the quadrilateral $XYZW$ to the area of
the parallelogram $ABCD$ equals $\frac{1}{2}\left( r_{1}+r_{2}\right) $
where,
\begin{equation}
r_{1}=-\frac{\left\vert
\begin{array}{ccc}
\frac{\lambda +1}{\lambda } & 1 & \frac{-1}{\mu +1} \\
-\frac{1+\nu +\nu \mu }{1+\nu } & -\mu -1 & 1 \\
\kappa & \frac{1+\kappa +\kappa \nu }{\nu } &
\end{array
\right\vert }{(2+\frac{1}{\lambda }-\frac{1}{1+\mu })(1+\mu +\frac{\nu \mu }
1+\nu })(1+2\kappa +\frac{\kappa +1}{\nu })}  \label{eq2}
\end{equation
and
\begin{equation}
r_{2}=\frac{\left\vert
\begin{array}{ccc}
\kappa & \frac{1+\kappa +\kappa \nu }{\nu } & 1 \\
\kappa & \frac{\lambda \kappa }{\lambda +1} & 1 \\
\frac{\lambda +1}{\lambda } & 1 & \frac{-1}{\mu +1
\end{array
\right\vert }{(2+\frac{1}{\lambda }-\frac{1}{1+\mu })(1+\kappa +\frac
\lambda \kappa }{1+\lambda })(1+2\kappa +\frac{\kappa +1}{\nu })}.
\label{eq3}
\end{equation}

Note that :
\[
\left\vert
\begin{array}{ccc}
\frac{\lambda +1}{\lambda } & 1 & \frac{-1}{\mu +1} \\
-\frac{1+\nu +\nu \mu }{1+\nu } & -\mu -1 & 1 \\
\kappa & \frac{1+\kappa +\kappa \nu }{\nu } &
\end{array
\right\vert =\left\vert
\begin{array}{ccc}
\frac{\lambda +1}{\lambda } & 1 & \frac{-1}{\mu +1} \\
-\frac{1+\nu +\nu \mu }{1+\nu }+(\mu +1)\frac{\lambda +1}{\lambda } & 0 & 0
\\
\kappa & \frac{1+\kappa +\kappa \nu }{\nu } &
\end{array
\right\vert
\]

and
\[
\left\vert
\begin{array}{ccc}
\kappa & \frac{1+\kappa +\kappa \nu }{\nu } & 1 \\
\kappa & \frac{\lambda \kappa }{\lambda +1} & 1 \\
\frac{\lambda +1}{\lambda } & 1 & \frac{-1}{\mu +1
\end{array
\right\vert =\left\vert
\begin{array}{ccc}
0 & \frac{1+\kappa +\kappa \nu }{\nu } & 1 \\
0 & \frac{\lambda \kappa }{\lambda +1} & 1 \\
\frac{\lambda +1}{\lambda }+\frac{\kappa }{\mu +1} & 1 & \frac{-1}{\mu +1
\end{array
\right\vert .
\]

After evaluating the determinants and substituting in (\ref{eq2}) and (\re
{eq3}) the result follows.
\end{proof}

\

In particular, if $\kappa =\lambda =\mu =\nu $ then the expressions in the
denominators of (\ref{eq1}) are equal to:
\[
\frac{\left( 2\lambda ^{2}+2\lambda +1\right) ^{3}}{\lambda ^{2}(\lambda
+1)^{2}},
\]

while the expressions in the numerators\ of (\ref{eq1}) are equal to:
\[
\frac{\left( 2\lambda ^{2}+2\lambda +1\right) ^{2}}{\lambda ^{2}(\lambda
+1)^{2}}.
\]

Hence, we have the following;

\begin{corollary}
If the sides $BC,CD,DA,AB$ of a parallelogram $ABCD$ are divided at $K,L,M,N$
$\ $in the ratio $\lambda :1,$ then the cevians $AK,BL,CM,DN$ form a
quadrilateral whose area is
\[
\frac{1}{2\lambda ^{2}+2\lambda +1}
\
times that of $ABCD.$
\end{corollary}

Substituting $\lambda =\frac{1}{p-1},p\geq 2$, we get the equivalent
expression $\frac{p^{2}-2p+1}{p^{2}+1},$ which was discovered by M. De
Villiers \cite{Villiers}. He did not provide any proof but pointed out the
following: \ "Since a square is affinely equivalent to a parallelogram, the
easiest way to derive and prove this formula is to consider the special case
of a square".

We leave to the reader to see what happens in another particular case of the
theorem: if $\kappa =\mu $ and $\lambda =\nu .$

Finally, the fascinating theorem of Routh-Steiner attracted the interest of
mathematics education researchers, probably because of the use of a dynamic
geometry software to rediscover the theorem, but it still attracts the
interest of mathematics researchers as well.

\bigskip

\qquad