Abstract: A subset of a convex body $B$ containing the origin in a Euclidean space is {\it parkable in $B$} if it can be translated inside $B$ in such a manner that the translate the origin. We provide characterizations of ellipsoids and of centrally symmetric convex bodies in Euclidean spaces of dimension $\ge 3$ based on the notion of parkability, answering several questions posed by G. Bergman.The techniques used, which are based on characterizations of Hilbert spaces among finite-dimensional Banach spaces in terms of their lattices of subspaces and algebras of endomorphisms, also apply to improve a result of W. Blaschke characterizing ellipsoids in terms of boundaries of illumination.

\section{#1}\setcounter{lemma}{0}}

\title{Parkable convex sets and finite-dimensional Hilbert spaces}
\author{Alexandru Chirvasitu\footnote{University of Washington, \url{chirva@uw.edu}}}

\begin{document}

\date{}

\maketitle

\begin{abstract}
A subset of a convex body $B$ containing the origin in a Euclidean space is {\it parkable in $B$} if it can be translated inside $B$ in such a manner that the translate the origin. We provide characterizations of ellipsoids and of centrally symmetric convex bodies in Euclidean spaces of dimension $\ge 3$ based on the notion of parkability, answering several questions posed by G. Bergman.

The techniques used, which are based on characterizations of Hilbert spaces among finite-dimensional Banach spaces in terms of their lattices of subspaces and algebras of endomorphisms, also apply to improve a result of W. Blaschke characterizing ellipsoids in terms of boundaries of illumination.
\end{abstract}

\noindent {\em Key words: convex body, Blaschke, centrally symmetric, ellipsoid, parkable set, Hilbert space, Banach space}

\vspace{.5cm}

\noindent{MSC 2010: 52A20, 52A21, 47L10, 46C15}

\section*{Introduction}

The present paper answers several questions raised in \cite{berg} regarding convex bodies in euclidean spaces. The setup is based on the following notion; it is introduced by G. Bergman in the process of studying ``efficient'' embeddings of metric spaces into other metric spaces.

All convex sets are understood to be subsets of some ambient Euclidean space $\bR^n$.

\begin{definition'}\label{def.parkable}
  Let $C\subseteq B$ be convex sets with $0\in B$. $C$ is \define{parkable in} $B$ if some translate of $C$ is still contained in $B$ and contains $0$.
\end{definition'}

The questions in \cite{berg} referred to above have to do with characterizing particularly nice convex bodies in $\bR^n$ (i.e. centrally symmetric or ellipsoids) by means of parkability. First recall

\begin{definition'}
  Let $C\subseteq \bR^n$ be a convex subset. A \define{center of symmetry} for $C$ is a point $p\in C$ such that
  \begin{equation*}
    C=2p-C:=\{2p-x\ |\ x\in C\}.
  \end{equation*}
  $C$ \define{has a center of symmetry} if a center of symmetry exists.

  $C$ is \define{centrally symmetric} if $0\in C$ is a center of symmetry for it.
\end{definition'}

\cite[Question 32]{berg} then reads as follows.

\begin{question'}\label{qu.sym}
  Let $C$ be a compact convex subset of $\bR^n$ for some $n>2$, with the property that for every centrally symmetric compact convex subset $B\subset \bR^n$ containing some translate of $C$, that translate is parkable in $B$.

Is it true that $C$ must have a center of symmetry?
\end{question'}

One of the main results is

\begin{theorem'}\label{th.sym}
  The answer to \Cref{qu.sym} is affirmative.
\mbox{}\hfill\ensuremath{\blacksquare}
\end{theorem'}

A more elaborate result has to do with recognizing ellipsoids by means of parkability. Recall

\begin{definition'}
  A \define{convex body} is a compact convex set with non-empty interior which contains $0$.
\end{definition'}

The starting point for the discussion that follows is the following result from \cite{berg}.

\begin{proposition'}\label{pr.ell}
 Let $B\subset \bR^n$ be a convex body for some $n>2$. Then, the
 following properties are successively weaker.
 \begin{enumerate}
 \renewcommand{\labelenumi}{(\roman{enumi})}
   \item $B$ is an ellipsoid centered at $0$.
   \item $B$ is centrally symmetric and its intersection with every hyperplane has a center of symmetry (provided it is not empty).
   \item Every closed convex subset of $B$ is parkable in $B$.
 \end{enumerate}
\end{proposition'}

It is then natural to ask \cite[Question 31]{berg}:

\begin{question'}\label{qu.ell}
  Is either of the implications from \Cref{pr.ell} reversible?
\end{question'}

In this context, the second main result is

\begin{theorem'}\label{th.ell}
  The three properties in \Cref{pr.ell} are equivalent to one another.
  \mbox{}\hfill\ensuremath{\blacksquare}
\end{theorem'}

The rest of this introduction contains a few remarks about the proofs.

\Cref{th.sym,th.ell} do not directly imply one another as such, but they are nevertheless interlinked through the methods used in their proofs. The same auxiliary lemma about plane compact convex sets, for instance, leads both to \Cref{th.sym} and to the fact that condition (iii) in \Cref{pr.ell} implies that $B$ is centrally symmetric (in other words, the partial implication (iii) $\Rightarrow$ (ii)).

Once we show that (iii) implies central symmetry, we can assume $B$ to be centrally symmetric throughout the rest of the proof of \Cref{th.ell}. As such, it is the unit ball of a unique Banach space structure $(\bR^n,\|\cdot\|)$ on $\bR^n$, and now functional-analytic techniques and results can be brought to bear.

More specifically, using the same auxiliary result referred to in passing above we first prove

\begin{proposition'}\label{pr.aux}
  Let $B$ be a convex centrally symmetric body $B$ satisfying condition (iii) from \Cref{pr.ell}. Then, for every linear hyperplane $H\subset \bR^n$, the non-empty intersections
  \begin{equation*}
    (H+x)\cap B,\quad x\in \bR^n
  \end{equation*}
have centers of symmetry. Moreover, these centers are collinear. \mbox{}\hfill\ensuremath{\blacksquare}
\end{proposition'}

Associating to $L$ the unique line through the origin containing the centers of symmetry from the statement of \Cref{pr.aux} extends, it turns out, to an inclusion-reversing involution on the lattice of subspaces of the Banach space $(\bR^n,\|\cdot\|)$ referred to above. This, together with a lattice-theoretic characterization of Hilbert spaces among Banach spaces due to Kakutani et al. leads to the conclusion that $(\bR^n,\|\cdot\|)$ is a Hilbert space. But this is equivalent to its unit ball being an ellipsoid, and the conclusion follows.

There are other results that might be of independent interest, such as an improvement of a result of Blaschke on characterizations of ellipsoids by means of light rays (\cite{bla}).

\section{Preliminaries}\label{se.prel}

\subsection{Convex geometry}\label{subse.conv}

Our main reference on the topic will be \cite{tomo}, to which we refer the reader for basic terminology.

We say that two convex bodies in a Euclidean space $\bR^n$ are {\it mutual translates} if one is the image of the other through some translation of $\bR^n$.

The following result will make an appearance several times in the
sequel; the reader can consult e.g. \cite{Rya15} and the references
therein for background on the result (which is also \cite[Theorem 3.1.3]{tomo}).

\begin{theorem}\label{th.transl}
 Let $K$ and $L$ be two compact convex subsets of $\bR^n$ for $n\ge 3$, and $2\le m\le {n-1}$ a positive integer. If the projections of $K$ and $L$ on every $m$-dimensional linear subspace of $\bR^n$ are mutual translates, then so are $K$ and $L$.
\mbox{}\hfill\ensuremath{\blacksquare}
\end{theorem}

As a consequence, we get (\cite[Corollary 3.1.5]{tomo})

\begin{corollary}\label{cor.centr}
  Let $m$ and $n$ be positive integers as in \Cref{th.transl}, and
  $B\subset \bR^n$ a compact convex subset. Then, $B$ has a center of symmetry
  if and only if its projection on every $m$-dimensional linear
  subspace of $\bR^n$ does.
\end{corollary}
\begin{proof}
  Having a center of symmetry is equivalent to $B$ and $-B$ being
  translates; we can now simply apply \Cref{th.transl} to this pair of
  convex sets.
\end{proof}

\subsection{Banach and Hilbert spaces}\label{subse.ban}

We will assume some basics on Banach and Hilbert spaces and bounded operators thereon, as covered e.g. in the introductory sections of \cite[Chapters 1 and 12]{rud} or in the first three chapters of \cite{con}.

Given a centrally symmetric convex body $B=-B$ in $\bR^n$, we can associate to it the unique Banach space structure on $\bR^n$ making $B$ the unit ball: the norm of $x\in \bR^n$ is defined to be
\begin{equation*}
  \|x\| = \|x\|_B = \inf\{r\ge 0\ |\ x\in rB\}.
\end{equation*}
The correspondence $B\mapsto \|\cdot\|_B$ is a bijection between centrally symmetric convex bodies and Banach space structures on $\bR^n$; the inverse map associates to a Banach space structure $(\bR^n,\|\cdot\|)$ its unit ball.

As \Cref{th.ell} above suggests, one of our main goals will be characterizing ellipsoids among convex bodies. In terms of the bijection $B\leftrightarrow \|\cdot\|_B$ the ellipsoids correspond to the {\it Hilbert space} structures on $\bR^n$, i.e. those Banach space structures whose underlying norm $\|\cdot\|$ arises from an inner product $\langle-,-\rangle$ via the usual formula
\begin{equation*}
  \|x\|^2 = \langle x,x\rangle,\ \forall x\in \bR^n.
\end{equation*}
For this reason, it will be important to have at our disposal results that allow for the recognition of Hilbert spaces among Banach spaces. One such tool is \cite[Theorem 1.1]{istr} (or rather a variant thereof, with real Banach spaces instead of complex ones):

\begin{theorem}\label{th.istr}
  Let $(\bR^n,\|\cdot\|)$ be a Banach space. The norm is induced by an inner product if and only if there exists an operation $T\mapsto T^*$ on the space $M_n(\bR)$ of endomorphisms of $\bR^n$ such that
  \begin{enumerate}
    \item $(T+S)^* = T^*+S^*$;
    \item $(T^*)^* = T$;
    \item $(TS)^* = S^*T^*$;
    \item $\|P\|\le 1$ if $P^2=P=P^*$, where $\|P\|$ is the norm on $M_n(\bR)$ induced by that on $\bR^n$.
  \end{enumerate}
\mbox{}\hfill\ensuremath{\blacksquare}
\end{theorem}

In other words, we have a recognition criterion for Hilbert space norms in terms of involutions $T\mapsto T^*$ on their algebras of endomorphisms.

\section{Universally parkable sets}\label{se.univ}

The aim of this section is to prove \Cref{th.sym} above, answering \Cref{qu.sym} in the affirmative. We recall the statement, after introducing the following notion relevant to the setup of the theorem.

\begin{definition}\label{def.univ_park}
  A compact convex set $C\subset \bR^n$ is {\it universally parkable} if for any centrally symmetric convex subset $B\subset \bR^n$ containing a translate of $C$, that translate is parkable in $B$ in the sense of \Cref{def.parkable}.
\end{definition}

\begin{theorem}\label{th.sym_bis}
  For a positive integer $n\ge 3$, every universally parkable compact
  convex subset $C\subset \bR^n$ has a center of symmetry.
\end{theorem}

We proceed through a series of auxiliary results. First off, we narrow down the class of centrally symmetric convex sets that witness the parkability of $C\subset \bR^n$. For this purpose as well as for use in the sequel, denote by $C_u$ the translate $C+u$ for $u\in \bR^n$.

\begin{lemma}\label{le.Cu}
  A compact convex subset $C\subset \bR^n$ is universally parkable if and only if for every $u\in \bR^n$ the translate $C_u$ is parkable in the convex hull $\co(-C_u\bigcup C_u)$.
\end{lemma}
\begin{proof}
  On the one hand, if $C$ is universally parkable, every $\co(-C_u\bigcup C_u)$ is certainly a centrally symmetric convex subset of $\bR^n$ containing a translate $C_u$, and hence the latter must be parkable in the former.

Conversely, every centrally symmetric compact convex $B\subset \bR^n$ containing a translate $C_u$ contains $\co(-C_u\bigcup C_u)$, so that if $C_u$ is parkable in the latter then it is parkable in the former as well.
\end{proof}

As a consequence, we obtain

\begin{corollary}\label{cor.Cu}
  If $C\subset \bR^n$ is universally parkable, then so are its projections on linear subspaces of $\bR^n$.
\end{corollary}
\begin{proof}
  This follows from the fact that the property from the statement of \Cref{le.Cu} is clearly invariant under taking orthogonal projections.
\end{proof}

In order to fix ideas, we now specialize \Cref{th.sym_bis} to $n=3$.

\begin{lemma}
  If \Cref{th.sym_bis} holds for $n=3$ then it holds in general.
\end{lemma}
\begin{proof}
  Starting with a universally parkable $C\subset \bR^n$ for some $n\ge 3$, \Cref{cor.Cu} ensures that its projections on all $3$-dimensional linear subspaces of $\bR^n$ are universally parkable. If \Cref{th.sym_bis} holds for $\bR^3$ we can conclude that these projections all have centers of symmetry, which according to \Cref{cor.centr} means that so does $C$.
\end{proof}

This now allows us to focus on the case $n=3$. Let $\bS\subset \bR^3$ be an origin-centered sphere whose radius
is large enough to ensure that no translates $C_u$ (and hence no reflections
$-C_u$) contain $0$ as $u$ ranges over $\bS$. The following lemma shows
that when $u\in \bS$ there is a certain rigidity in how one may
translate $C_u$ inside $\co(-C_u\bigcup C_u)$ so as to engulf $0$.

\begin{lemma}\label{le.unique}
  Let $C\subset \bR^3$ be as in the statement of \Cref{th.sym_bis} For
  any $u\in \bS$, the direction of a vector $0\ne v\in \bR^3$ such that
  \begin{equation*}
    0\in C_u+v\subset \co(-C_u\bigcup C_u)
  \end{equation*}
 is uniquely determined.
\end{lemma}
\begin{proof}
  Let $u\in \bS$ and $0\ne v\in \bR^3$ be vectors as in the statement, and let
  $H\subset \bR^3$ be a $2$-dimensional linear such that the
  projection $C_u|H$ does not contain $0\in H$ and is not collinear
  with the origin (i.e. is not a segment whose supporting line
  contains $0$). Note that such planes
  $H$ form an open subset of the Grassmannian $\cG(3,2)$ (because
  their defining property is open), and by our
  choice of radius for $\bS$ this subset of $\cG(3,2)$ is non-empty.

In the plane $H$ the part of the boundary $\partial \co(-C_u|H\bigcup
C_u|H)$ that is not contained in the boundaries of $C_u|H$ and
$-C_u|H$ consists of two open segments forming two opposite edges of a
parallelogram centered at $0$:
  \begin{equation*}
      \begin{tikzpicture}[auto,baseline=(current  bounding
        box.center),scale=1]
        \node (c) at (2.2,.8) {$\scriptstyle C_u|H$};
        \node (-c) at (-2.3,-.8) {$\scriptstyle -C_u|H$};
        \node (aux) at (-1.125,-.425) {};
        \draw (0,0) node[circle, inner sep=1pt, fill=black,
        label={right:{$\scriptstyle 0$}}] (0) {};
        \draw[-] (2,1.5) .. controls (2.5,1.7) and (3.3,1.2) .. (3,.7);
        \draw[-] (3,.7) .. controls (2.7,.2) and (2.625,.25) .. (2.5,.2);
        \draw[-] (2.5,.2) to[bend left=9] (1.5,.4);
        \draw[-] (1.5,.4) .. controls (1,.6) and (1.75,1.4) .. (2,1.5);
        \draw[-] (-2,-1.5) .. controls (-2.5,-1.7) and (-3.3,-1.2) .. (-3,-.7);
        \draw[-] (-3,-.7) .. controls (-2.7,-.2) and (-2.625,-.25) .. (-2.5,-.2);
        \draw[-] (-2.5,-.2) to[bend left=9] (-1.5,-.4);
        \draw[-] (-1.5,-.4) .. controls (-1,-.6) and (-1.75,-1.4) .. (-2,-1.5);
        \draw[-] (2,1.5) -- (-2.5,-.2);
        \draw[-] (-2,-1.5) -- (2.5,.2);
        \draw[->] (0) to node[pos=.3,auto,swap] {$\scriptstyle w$} (aux);
        \draw (-.25,.65) node[circle, inner sep=1pt, fill=black] () {};
        \draw (.25,-.65) node[circle, inner sep=1pt, fill=black] () {};
        \draw[-] (.5,-1.3) -- (-.5,1.3);
        \draw (-.5,1.3) node[circle, inner sep=0pt, fill=black,
        label={right:{$\scriptstyle L$}}] () {};
      \end{tikzpicture}
  \end{equation*}
More formally, as the diagram illustrates, the two segments in question are the tangents to the boundary $\partial\co(-C_u|H\bigcup C_u|H)$ at the points where a hyperplane (i.e. line) $L$ that separates $C_u|H$ and $-C_u|H$ intersects said boundary.

Denoting by $w$ the orthogonal projection of $v$ on $H$, it follows
from the above remark about the boundary of $\co(-C_u|H\bigcup C_u|H)$
that $w$ must be parallel to the two segments and point away from
$C_u|H$. This determines the direction of $w\in H$ uniquely, and since
$w$ was the orthogonal projection of $v$ on any one of the elements of
an open subset of $\cG(3,2)$, we deduce the desired uniqueness of the
direction of $v$.
\end{proof}

\begin{remark}\label{re.unique}
   Note that the proof of \Cref{le.unique} actually shows more than the statement claims: it shows that the unique direction in which $C_u$ can be translated without leaving $\co(-C_u\bigcup C_u)$ is that of $\varphi(u)$.
\end{remark}

This gives us, for every $u\in \bS$, a unique unit vector $\varphi(u)=\frac
v{\|v\|}\in \bS_1$ (unit sphere centered at $0$) such that translation
in its direction will position $C_u$ so that it contains $0$. The map
$\varphi$ is in fact very well-behaved. Recall that an {\it odd} map
between two spheres is one that sends antipodes to antipodes. With this in hand, we have

\begin{lemma}\label{le.cont}
  The map $\varphi:\bS\to \bS_1$ deduced from \Cref{le.unique} as in
  the discussion above is continuous and odd.
\end{lemma}
\begin{proof}
  The fact that $\varphi$ is odd is immediate: if $-\varphi(u)$
  belongs to $C_u$ then $\varphi(u)$ belongs to $-C_u$, meaning that
  translation of $-C_u$ by $-\varphi(u)$ will position the former so
  that it contains $0$.

As for the continuity of $\varphi$, note first that because of our choice of $\bS$ (such that no $C_u$ contains the origin) there is a positive lower bound on the length of vectors $v$ such that $0\in C_u+v$ as $u$ ranges over $\bS$. In conclusion, there is some $t>0$ such that for all $u\in \bS$ we have
\begin{equation}\label{eq:cont}
  C_u+t \varphi(u)\in \co(-C_u\bigcup C_u).
\end{equation}
The condition
\begin{equation*}
  C_u+tv\in \co(-C_u\bigcup C_u)
\end{equation*}
is closed in $(u,v)\in \bS\times \bS_1$, and together with \Cref{re.unique} above, \Cref{eq:cont} shows that the set of pairs $(u,v)$ that satisfy it is precisely the graph of $\varphi:\bS\to \bS_1$. Since we now know that the graph of the map $\varphi$ between compact spaces is closed, we can conclude that it is continuous.
\end{proof}

\begin{corollary}\label{cor.onto}
  The map $\varphi:\bS\to \bS_1$ from \Cref{le.cont} is onto.
\end{corollary}
\begin{proof}
  Indeed, \Cref{le.cont} shows that it is a continuous odd map between $2$-spheres. By one
  version of the Borsuk-Ulam theorem (e.g. the main result of \cite{dold}) it follows that the map is not nullhomotopic and hence, because the complement of a point in the $2$-sphere is contractible,  must be onto.
\end{proof}

\begin{lemma}\label{le.tough}
  For every $u\in \bS$, the orthogonal projection of $C_u$ on $H=\varphi(u)^\perp$ is centrally symmetric.
\end{lemma}
\begin{proof}
We have to show that the projections of $C_u$ and $-C_u$ on $H$ coincide. If they do not, then there is a point $x$ on the boundary $\partial(C_u|H)$ that does not belong to $-C_u|H$. The line in $\bR^3$ through $x$ and orthogonal to $H$ (and hence parallel to $\varphi(u)$) intersects $-C_u\bigcup C_u$ along a segment (possibly degenerate, i.e. a single point) contained in the boundary $\partial C_u$.
  \begin{equation*}
      \begin{tikzpicture}[auto,baseline=(current  bounding
        box.center),scale=1]
        \node (c) at (-.7,1.4) {$\scriptstyle C_u$};
        \node (-c) at (.8,-1.5) {$\scriptstyle -C_u$};
        \draw (0,0) node[circle, inner sep=1pt, fill=black, label={45:{$\scriptstyle 0$}}] () {};
        \draw[-] (-2,2) -- (-2,1.5);
        \draw[-] (-2,2) .. controls (-2,3) and (.5,2) .. (.5,1);
        \draw[-] (-2,1.5) .. controls (-2,1) and (.5,0) .. (.5,1);
        \draw[-] (2,-2) -- (2,-1.5);
        \draw[-] (2,-2) .. controls (2,-3) and (-.5,-2) .. (-.5,-1);
        \draw[-] (2,-1.5) .. controls (2,-1) and (-.5,0) .. (-.5,-1);
        \draw[-] (-3,0) -- (3,0);
        \node (H) at (2.8,.2) {$\scriptstyle H$};
        \draw[->] (2,2) to node[pos=.5] {$\scriptstyle \phi(u)$} (2,1);
        \draw[-] (-2,2.5) -- (-2,-.3);
        \draw (-2,0) node[circle, inner sep=1pt, fill=black, label={135:{$\scriptstyle x$}}] () {};
      \end{tikzpicture}
  \end{equation*}
Translation by $\varphi(u)$ will move one of the endpoints of the segment outside of $\co(-C_u\bigcup C_u)$, contradicting
\begin{equation*}
  C_u+\varphi(u)\subset \co(-C_u\bigcup C_u).
\end{equation*}
It follows from the contradiction that $-C_u|H=C_u|H$, as desired.
\end{proof}

We are now ready to complete the proof of the main result of this section.

\begin{proof_of_symbis}
  \Cref{le.tough,cor.onto} show that the orthogonal projection of $C$ on every hyperplane in $\bR^3$ has a center of symmetry and hence, according to \Cref{cor.centr}, so does $C$.
\end{proof_of_symbis}

\section{Parkability and finite-dimensional Banach spaces}\label{se.park}

Here, we prove \Cref{th.ell}. Recall the statement:

\begin{theorem}\label{th.ell_bis}
For a convex body $B\subset \bR^n$, $n\ge 3$ the following properties
are equivalent.
 \begin{enumerate}
 \renewcommand{\labelenumi}{(\roman{enumi})}
   \item $B$ is an ellipsoid centered at $0$.
   \item $B$ is centrally symmetric and its intersection with every hyperplane has a center of symmetry (provided it is not empty).
   \item Every closed convex subset of $B$ is parkable in $B$.
 \end{enumerate}
\end{theorem}

We already know from \Cref{pr.ell} that (i) implies (ii), which in
turn implies (iii); hence, it suffices to go backwards. We once again
specialize to $n=3$ in order to simplify some of the proofs and
language.

\begin{proposition}\label{pr.n=3}
  If \Cref{th.ell_bis} holds for $n=3$ then it holds in general.
\end{proposition}
\begin{proof}
  Property (iii) is clearly preserved by passing to orthogonal
  projections from $\bR^n$ down to its three-dimensional subspaces. If
  the implication (iii) $\Rightarrow$ (i) holds in $\bR^3$, then we
  know that all orthogonal projections of $B$ on such subspaces are
  ellipsoids. This implies that $B$ itself is an ellipsoid by \cite[Theorem 3.1.7]{tomo}.
\end{proof}

As explained in the introduction, the first priority will be to prove that (iii)
implies the property of being centrally symmetric so as to be able to think of the
convex body as the unit ball of a Banach space structure on
$\bR^n$. The following lemma provides the first step in this direction.

\begin{lemma}\label{le.2impln}
  If the implication (iii) $\Rightarrow$ (centrally symmetric) holds
  in $\bR^2$, then it holds for all $\bR^n$, $n\ge 3$.
\end{lemma}
\begin{proof}
  Let $B\subset \bR^n$, $n\ge 3$ be a convex body satisfying condition
  (iii) of the theorem, and suppose (iii) implies central symmetry in
  the plane.

Clearly, (iii) is preserved upon projecting down to any plane, and
hence, by assumption, all projections of $B$ on planes are centrally
symmetric. The conclusion that $B$ itself is then follows from
\Cref{cor.centr} (or rather a variant thereof whereby the center of
symmetry is in fact the origin).
\end{proof}

\begin{remark}\label{re.hyp}
As seen in \cite[Lemma 29]{berg}, condition (iii) of \Cref{th.ell_bis}
is equivalent to the fact that every hyperplane section of $B$ is
parkable. We will henceforth use this equivalence without further
comment when convenient.
\end{remark}

\begin{proposition}\label{pr.2}
  Let $B\subset \bR^2$ be a convex body such that every intersection
  of $B$ with a line is parkable in $B$. Then, $B$ is centrally
  symmetric.
\end{proposition}
\begin{proof}
Let $v\in \bR^2$ be a non-zero vector, and $H_1$ and $H_2$ the two $B$-supporting lines that are orthogonal to $v$.

{\bf Claim: The convex hull of the segments $H_1\cap B$ and $H_2\cap
  B$ contains the origin.} Indeed, if $p_1\in H_1\cap B$ and $p_2\in
H_2\cap B$ are two points on the boundary of $B$, then by our
hypothesis some translate $\overline{p_1p_2}+w\subset B$ contains the
origin. Because $H_i$ are supporting lines the translates $p_1+w$ and
$p_2+w$ belong to $H_1$ and $H_2$ respectively (i.e. if it is
non-zero, then $w$ is parallel to $H_1$ and $H_2$); this finishes the proof of the claim.
  \begin{equation*}
      \begin{tikzpicture}[auto,baseline=(current  bounding
        box.center),scale=.8]
        \draw (0,-1) node[circle, inner sep=0pt, fill=black, label={225:{$\scriptstyle H_1$}}] () {};
        \draw (-1.2,.2) node[circle, inner sep=1pt, fill=black, label={225:{$\scriptstyle p_1$}}] (p1) {};
        \draw (4,-2) node[circle, inner sep=0pt, fill=black, label={45:{$\scriptstyle H_2$}}] () {};
        \draw (2.5,-.5) node[circle, inner sep=1pt, fill=black, label={45:{$\scriptstyle p_2$}}] (p2) {};
        \draw (.3,.3) node[circle, inner sep=1pt, fill=black, label={45:{$\scriptstyle 0$}}] () {};
        \draw[-] (-2,1) -- (-1,0);
        \draw[-] (1,1) -- (3,-1);
        \draw[-] (-2,1) -- (1,1);
        \draw[-] (-1,0) -- (3,-1);
        \draw[-] (-2,1) .. controls (-2.5,1.5) and (.5,1.5) .. (1,1);
        \draw[-] (-1,0) .. controls (-.5,-.5) and (3.5,-1.5) .. (3,-1);
        \draw[-] (-2.5,1.5) -- (0,-1);
        \draw[-] (0,2) -- (4,-2);
        \draw[->] (4,.5) to node[pos=.5,auto,swap] {$\scriptstyle w$} (3.5,1);
      \end{tikzpicture}
  \end{equation*}
Given the claim, the conclusion now follows from the technical \Cref{le.supports} below.
\end{proof}

\begin{lemma}\label{le.supports}
Let $B\subset \bR^2$ be a convex body. Suppose that for every non-zero $v\in \bR^2$ the two supporting lines $L_1=L_1(v)$ and $L_2=L_2(v)$ of $B$ that are orthogonal to $v$ satisfy
\begin{equation*}
  0\in \co((L_1\cap B)\cup(L_2\cap B)).
\end{equation*}
Then, $B$ is centrally symmetric.
\end{lemma}
\begin{proof}
We have to show that $B$ and $-B$ coincide. Assuming they do not, the condition on $B$ at least ensures that the union $-B\cup B$ is convex. Indeed, otherwise for some $p\in \partial B\cap \partial(-B)$ some $B$-supporting line $L$ through $p$ would not be $-B$-supporting, and hence $v\perp L$ would violate the hypothesis.

In conclusion, $B$ and $-B$ admit common parallel supporting lines $L$ and $-L$ at two points $p$ and $-p$ respectively in the intersection $\partial B\cap \partial(-B)$. We will moreover choose coordinates in $\bR^2$ so that the segment $\overline{(-p)p}$ is horizontal and $L$ is vertical. We will derive a contradiction from the assumption that the portions of $\partial B$ and $\partial(-B)$ lying in the lower half plane are distinct.
  \begin{equation*}
      \begin{tikzpicture}[auto,baseline=(current  bounding
        box.center),scale=1]
        \draw (0,0) node[circle, inner sep=1pt, fill=black, label={90:{$\scriptstyle 0$}}] () {};
        \draw (-1,0) node[circle, inner sep=1pt, fill=black, label={135:{$\scriptstyle -p$}}] () {};
        \draw (1,0) node[circle, inner sep=1pt, fill=black, label={45:{$\scriptstyle p$}}] () {};
        \draw (-1,-2) node[circle, inner sep=0pt, fill=black, label={225:{$\scriptstyle -L$}}] () {};
        \draw (1,-2) node[circle, inner sep=0pt, fill=black, label={315:{$\scriptstyle L$}}] () {};
        \draw (.5,-1) node[circle, inner sep=0pt, fill=black, label={100:{$\scriptstyle B$}}] () {};
        \draw (-.5,-2) node[circle, inner sep=0pt, fill=black, label={85:{$\scriptstyle -B$}}] () {};
        \draw[-] (-1,1) -- (-1,-2);
        \draw[-] (1,1) -- (1,-2);
        \draw[-] (-2,0) -- (2,0);
        \draw[-] (-1,0) .. controls (-1,-.5) and (.1,-1) .. (.5,-1);
        \draw[-] (.5,-1) .. controls (.9,-1) and (1,-.5) .. (1,0);
        \draw[-] (-1,0) .. controls (-1,-.5) and (-.9,-2) .. (-.5,-2);
        \draw[-] (-.5,-2) .. controls (-.1,-2) and (1,-2) .. (1,0);
      \end{tikzpicture}
  \end{equation*}
We can now identify $\overline{(-p)p}$ with a closed interval $[-a,a]$, $a>0$ and the portions of $\partial B$ and $\partial(-B)$ depicted above with graphs of convex functions $f<0$ and $g<0$ respectively defined on $(-a,a)$. Moreover, the assumption $B\ne -B$ translates without loss of generality to $g(x)<f(x)$ for all $x\in (-a,a)$.

In this setup, the meaning of the hypothesis is that whenever a line through the origin intersects the graphs of $f$ and $g$ at smooth points $(x,f(x))$ and $(y,f(y))$ respectively, the tangents to the two graphs at those points are parallel; equivalently, $f'(x)=g'(y)$.

If, in the previous paragraph, we choose $x$ to be negative, then $y<x$. Because $g$ is convex, its derivative is non-decreasing, and hence
\begin{equation*}
  f'(x)=g'(y)\le g'(x)
\end{equation*}
 (always assuming we are working with points where the derivatives exist, which is the case for all but countably many of the elements in the domain $(-a,a)$).

In conclusion, whenever both $f$ and $g$ are differentiable at $x\in (-a,0]$ we have $f'(x)\le g'(x)$. This, together with
\begin{equation*}
  f(0)=\int_{-a}^0f'(x)\ \mathrm{d}x\le g(0)=\int_{-a}^0g'(x)\ \mathrm{d}x
\end{equation*}
contradicts our assumption that $g(0)<f(0)$.
\end{proof}

Using \Cref{pr.2}, we obtain

\begin{corollary}\label{cor.iii_implies_sym}
  Let $B\subset \bR^n$, $n\ge 3$ be a convex body satisfying condition
  (iii) of \Cref{th.ell_bis}. Then, $B$ is centrally symmetric.
\end{corollary}
\begin{proof}
  This follows immediately from \Cref{le.2impln} via \Cref{pr.2}.
\end{proof}

We now know that all convex bodies satisfying the weakest condition (iii) of \Cref{th.ell_bis} are centered at the origin, and as a consequence we henceforth specialize the discussion to such bodies. Next, our goal will be to prove an $n=3$ version of \Cref{pr.aux}, announced above.

\begin{proposition}\label{pr.aux_bis}
  Let $B$ be a convex centrally symmetric body $B\subset \bR^3$ satisfying condition (iii) from \Cref{pr.ell}. Then, for every linear hyperplane $H\subset \bR^3$, the non-empty intersections
  \begin{equation*}
    (H+x)\cap B,\quad x\in \bR^3
  \end{equation*}
have centers of symmetry. Moreover, these centers are collinear.
\end{proposition}

Before embarking on a proof, we need some preparation. As explained in the discussion following \cite[Question 31]{berg}, convex bodies $B$ with parkable hyperplane sections have the following property: for any linear hyperplane $H\subset \bR^3$ there is a line $L\subset \bR^3$ such that $B$ has a supporting translate of $L$ at every point of $H\cap \partial B$ (here, `supporting' is the natural extension to lines of the term `supporting hyperplane'; it simply means that the line intersects only the boundary of $B$).

It is also explained in loc. cit. how this property is in a sense dual to one considered by Blaschke \cite[pp. 157-159]{bla} in relation to illuminating convex bodies. In this latter reference, (smooth) convex bodies are considered with the property that when illuminated with parallel light rays, the boundary curve of the illuminated region is planar. In other words, given a line $L$, the intersections of the $B$-supporting translates of $L$ with $B$ form a planar curve in $\partial B$. This is the precise opposite of the situation in the previous paragraph.

Motivated by these considerations, we label the two dual properties described above.

\begin{definition}\label{def.bla}
  Let $B\subset \bR^3$ be a convex body containing $0$ in its interior.

$B$ has {\it the Blaschke property} (or {\it is Blaschke}, for short) if for any line $L$ there exists a linear hyperplane $H\subset \bR^3$ such that all points in $H\cap \partial B$ lie on some $B$-supporting translate of $L$.

$B$ has {\it the dual Blaschke property} (or {\it is dual Blaschke}) if for any linear hyperplane $H\subset \bR^3$ there is a line $L\in \bR^3$ such that all points in $H\cap \partial B$ lie on some $B$-supporting translate of $L$.
\end{definition}

With these terms in place, the discussion preceding the definition amounts to (\cite[p. 257]{berg})

\begin{proposition}\label{pr.co_bla}
  A centrally symmetric convex body in $\bR^3$ satisfying condition (iii) of \Cref{th.ell_bis} is dual Blaschke.
\mbox{}\hfill\ensuremath{\blacksquare}
\end{proposition}

Our next aim is to prove the dual version of this result.

\begin{proposition}\label{pr.bla}
  A centrally symmetric convex body in $\bR^3$ satisfying condition (iii) of \Cref{th.ell_bis} is Blaschke.
\end{proposition}

We now introduce some tools necessary in the proof of \Cref{pr.bla}. Consider a convex body $B$ as in \Cref{pr.co_bla}. For every unit vector $v\in \bR^3$, denote by $H_v$ the plane orthogonal to $v$. We define the subset $\psi(v)=\psi_B(v)$ of the unit sphere to consist of all unit vectors $w$ such that
\begin{equation*}
  w\cdot v>0 \text{ and the line } \{p+tw\ |\ t\in \bR\} \text{ is }B-\text{supporting for every }p\in \partial B\cap H_v.
\end{equation*}
In other words, it is the set of unit vectors that make acute angles with $v$ and pointing along the directions of those lines $L$ associated to the plane $H_v$ as in \Cref{pr.co_bla,def.bla}.

\Cref{pr.co_bla} is what makes $\psi(v)$ non-empty. Note moreover that since the lines supporting $B$ at points in $\partial B\cap H$ cannot be contained in $H$ (for any plane $H$), $\psi(v)$ is closed in the unit sphere $\bS^2$. Finally, it is also {\it convex} as a subset of the sphere, in the sense that for $w,w'\in \psi(v)$ the rescaled convex combinations
\begin{equation*}
  \frac{tw+(1-t)w'}{\|tw+(1-t)w'\|},\ t\in (0,1)
\end{equation*}
all belong to $\psi(v)$. All in all, we have defined a map
\begin{equation*}
  \psi=\psi_B:\bS^2\to\text{ closed convex subsets of }\bS^2.
\end{equation*}
The map $\psi$ has a number of other properties that are easy to check:

\begin{lemma}\label{le.odd}
  For a convex body $B$ as in \Cref{pr.co_bla}, the map $\psi=\psi_B$ is odd in the sense that
  \begin{equation*}
    \psi(-v)=-\psi(v)
  \end{equation*}
and upper semicontinuous in the sense that its graph
\begin{equation*}
  \{(v,w)\in \bS^2\times \bS^2\ |\ w\in \psi(v)\}
\end{equation*}
is closed in $\bS^2\times \bS^2$.
\end{lemma}
\begin{proof}
  The property of being odd is immediate from the definition of $\psi$: $v$ and $-v$ define the same plane $H_v=H_{-v}$, and hence the same set of lines supporting $B$ at the points of $\partial B\cap H_v$. The upper semicontinuity claim is similarly routine, using the observation that the set of $B$-supporting affine lines is closed in the set of all affine lines of $\bR^3$.
\end{proof}

\begin{lemma}\label{le.psi_onto}
  An odd, upper semicontinuous map
  \begin{equation*}
       \psi=:\bS^2\to\text{ closed convex subsets of }\bS^2
  \end{equation*}
is onto in the sense that every $w\in \bS^2$ belongs to some set $\psi(v)$, $v\in \bS^2$.
\end{lemma}
\begin{proof}
  Assume the contrary, with, say, $w\in \bS^2$ lying outside all of the sets $\psi(v)$ as $v$ ranges over $\bS^2$. Flowing along great circles through $w$ down to the equator in the plane orthogonal to $w$, we can deform $\psi$ into an odd, upper semicontinuous map
  \begin{equation*}
    \psi=:\bS^2\to\text{ closed convex subsets of }\bS^1.
  \end{equation*}
The existence of such a map contradicts the multivalued Borsuk-Ulam theorem as stated e.g. in \cite[Corollary 2.4]{multi_BU}.
\end{proof}

We have now effectively proven \Cref{pr.bla} above.

\begin{proof_of_bla}
Using the notation introduced above, the statement amounts to showing that $\psi_B$ is onto; this is precisely what \Cref{le.psi_onto} does.
\end{proof_of_bla}

This provides us with the necessary tools to attack \Cref{pr.aux_bis}.

\begin{proof_of_aux_bis}
  The statement stipulates a condition that is closed among planes in $\cG(3,2)$, so it suffices to prove that this condition holds for a dense set of planes. Specifically, we will prove it for planes $H$ with the following property:

{\it The two supporting planes of $B$ that are parallel to $H$ each intersect $B$ at a single point.}

Now fix such a plane $H$, and let $p$ and $q=-p$ be the two points where the two $B$-supporting translates of $H$ intersect $B$. Let also $K$ be a planar section $(H+x)\cap B$ with non-empty relative interior. Finally, fix a line $L\in H$.

In the affine plane $H+x$, the convex body $K$ has two supporting lines parallel to $L$; denote these by $L_1$ and $L_2$.

{\bf Claim: The convex hull of $L_i\cap K$ contains the point $r=(H+x)\cap \overline{pq}$.} Given the claim, we can finish the proof of the proposition as follows.

Since the claim holds for any line $L\in H$, we can apply \Cref{le.supports} inside the affine plane $H+x$, with the point $r$ regarded as the origin, to conclude that $K$ has this point as its center of symmetry. In conclusion, we obtained the desired result: all planar sections of $B$ that are parallel to $H$ are centrally symmetric along points on the line containing $p$ and $q$.

It thus remains to prove the claim; this goal will take up the rest of the proof.

According to \Cref{pr.bla}, there is a plane $H'\subset \bR^3$ with the property that every point in $H'\cap \partial B$ is contained in a $B$-supported line parallel to $L$. By our choice of $H$ (so that its two $B$-supporting translates meet $B$ at single points $p$ and $q$ respectively), we have $p,q\in H'$. As a consequence of this, the endpoints of the segment $H'\cap K$ are contained in the $B$-supporting translates $L_1$ and $L_2$; this concludes the proof.
\end{proof_of_aux_bis}

\begin{remark}\label{re.32to31}
  Note that given the plane $H\subset \bR^3$, the line containing the centers of symmetry for the non-empty planar sections $(H+x)\cap B$ is uniquely determined; this provides a map from the Grassmannian of planes $\cG(3,2)$ to the Grassmannian $\cG(3,1)$ of lines in $\bR^3$.
\end{remark}

We now take the first step towards constructing a map $\cG(3,1)\to \cG(3,2)$ in the opposite direction to that of \Cref{re.32to31}.

\begin{lemma}\label{le.aux_bis_dual}
  Suppose we are under the assumptions of \Cref{pr.aux_bis}, and $L$ is the line containing the centers of symmetry of $(H+x)\cap B$ for $x$ ranging over $\bR^3$. Then, the midpoint of every non-empty intersection $(L+y)\cap B$ belongs to $H$.
\end{lemma}
\begin{proof}
  \Cref{pr.aux_bis} can be restated as saying that the linear involution $T$ of $\bR^3$ that acts as the identity along $L$ and as $x\mapsto -x$ on $H$ preserves the convex body $B$. The involution $-T$ also preserves $B$, acts as the identity on $H$, and as $x\mapsto -x$ along $L$. In other words, $-T$ preserves $H$ and reverses every segment $(L+y)\cap B$; this implies the conclusion.
\end{proof}

We can now prove the dual version of \Cref{pr.aux_bis}.

\begin{proposition}\label{pr.aux_bis_dual}
  Suppose we are under the hypotheses of \Cref{pr.aux_bis}. For any line $L\subset \bR^3$, the midpoints of the non-empty intersections $(L+y)\cap B$ are coplanar.
\end{proposition}
\begin{proof}
  \Cref{le.aux_bis_dual} above already proves the statement for those lines $L$ arising as images of planes $H\in \cG(3,2)$ through the map $\cG(3,2)\to \cG(3,1)$ from \Cref{re.32to31}. It is thus sufficient to show that this map is onto. In other words, we have to prove that {\it every} line $L$ contains the symmetry centers of the non-empty planar sections $(H+x)\cap B$ for some $H\in \cG(3,2)$.

Consider the map $\phi=\phi_B:\bS^2\to \bS^2$ defined as follows. For a unit vector $v\in \bS^2$, let $H=H_v$ be the plane orthogonal to $v$, and take $\phi(v)$ be the unit vector making an acute angle with $v$ and pointing along the line $L$ that contains the centers of symmetry of $(H+x)\cap B$, $x\in \bR^3$.

Clearly, $\phi$ is continuous and odd, in the sense that $\phi(-v)=-\phi(v)$. We can now apply \Cref{le.psi_onto} in the simpler case of ordinary (rather than multivalued) maps to conclude via Borsuk-Ulam that $\phi$ is onto.
\end{proof}

\Cref{pr.aux_bis,pr.aux_bis_dual} allow us to introduce the following

\begin{notation}\label{not.prime}
  For a linear subspace $H\subset \bR^3$ define $H'$ to be
  \begin{itemize}
    \item the zero subspace $\{0\}$ if $H=\bR^3$;
    \item the line containing the centers of symmetry of $(H+x)\cap B$ if $H$ is a plane;
    \item the plane containing the midpoints of $(H+y)\cap B$ if $H$ is a line;
    \item $\bR^3$ if $H=\{0\}$.
  \end{itemize}
\end{notation}

The operation $H\mapsto H'$ on the subspace of $\bR^3$ has the following properties.

\begin{lemma}\label{le.inv}
  \begin{enumerate}
  \renewcommand{\labelenumi}{(\arabic{enumi})}
    \item For any two subspaces $L,H$ of $\bR^3$ we have $L\subseteq H\Rightarrow L'\supseteq H'$;
    \item For every subspace $H\subseteq \bR^3$ we have $H''=H$;
    \item For every subspace $H\subseteq \bR^3$ we have $H\cap H'=\{0\}$.
  \end{enumerate}
\end{lemma}
\begin{proof}
  In all three statements, the interesting cases are those where the subspaces in question are 1- or 2-dimensional. For this reason, we treat only these cases in the proof.

{\bf (1)} Here, we may as well assume that $L$ is a line contained in the plane $H$. $H'$ is by definition the line containing the symmetry centers of the sections $(H+x)\cap B$. Every such center is the midpoint of the segment $(L+x)\cap B$, and is thus on $L'$ by the definition of the latter.

{\bf (2)} The identity $H=H''$ for planes is simply a restatement of \Cref{le.aux_bis_dual} above. Applying the operation $\bullet\mapsto \bullet'$ once more we get $H'=H'''$ for planes $H$, and hence $L=L''$ for lines $L$ follows from the observation that every line arises as $H'$ for some plane $H$ (this latter claim is essentially the surjectivity of the map $\phi_B$ from the proof of \Cref{pr.aux_bis_dual}).

{\bf (3)} This is immediate: when $H$ is a plane, for instance, there are non-empty sections $(H+x)\cap B$ for non-zero $x\in \bR^3$, and hence there are points of the line $H'$ that are not contained in $H$. A similar argument applies when $H$ is a line (or alternatively, by part (2), in that case $H=H''$ and we can apply the previous argument to the plane $H'$).
\end{proof}

Properties (1), (2) and (3) in \Cref{le.inv} are, according to \cite[Theorem 1]{KM44}, precisely what is necessary in order to ensure that there is a Hilbert space structure on $\bR^3$ for which $H\mapsto H'$ is the orthogonal complement operation. Such a Hilbert space structure is given by an inner product of the form
\begin{equation*}
  \langle v,w\rangle = v\cdot Aw
\end{equation*}
where $A\in M_3=M_3(\bR)$ is a positive operator and $\cdot$ is the usual dot product. Since the image of an ellipsoid through a positive operator $A:\bR^3\to \bR^3$ is again an ellipsoid, we may as well assume (for the purposes of proving \Cref{th.ell_bis}) that $A$ is the identity matrix and hence $K\mapsto K'$ is the usual orthogonal complement operation in $\bR^3$; we will do this throughout the rest of the section, in order to simplify the discussion.

Now consider the usual transposition operation $T\mapsto T^t$ on $M_3$. It satisfies conditions (1), (2) and (3) of \Cref{th.istr} above. Our goal will be to apply that result to the Banach space $(\bR^3,\|\cdot\|_B)$ induced by the centrally symmetric convex body $B$ as explained in \Cref{subse.ban}. In preparation for that, note that the symmetric idempotents $P\in M_3$ (i.e. those satisfying $P^2=P=P^t$) are precisely those whose range and kernel are orthogonal.

\begin{lemma}\label{le.istr_applies}
  If $\|\cdot\|$ denotes the norm induced on $M_3$ by $\|\cdot\|_B$, then $\|P\|^2\le 1$ for all symmetric idempotents $P\in M_3$
\end{lemma}
\begin{proof}
  If the range of $P$ is three- or zero-dimensional then $P$ is the identity or the zero operator respecticely, so there is nothing to prove. It thus remains to prove the claim when $\dim(\mathrm{Im}(P))$ is $2$ or $1$. We tackle the former case; the latter is entirely analogous.

If $H\in \cG(3,2)$ is the range of $P$, then $P$ is the orthogonal projection on $H$. Since, as explained in the discussion following \Cref{le.inv}, we are assuming that $H'=H^\perp$, the convex body $B$ (i.e. the unit ball of $(\bR^3,\|\cdot\|_B)$) is contained in the orthogonal cylinder based on $H\cap B$. But this means that if $\|x\|_B=1$ (i.e. $x\in \partial B$) then the orthogonal projection $Px$ on $H$ is contained in $B\cap H$, and hence $\|Px\|_B\le 1$. In conclusion,
\begin{equation*}
  \|P\| = \sup_{\|x\|_B\le 1}\|Px\|_B\le 1.
\end{equation*}
This finishes the proof.
\end{proof}

This concludes the preparatory material necessary for the proof of the main result of this section.

\begin{proof_of_ellbis}
  As mentioned after the statement of the theorem, it suffices to prove that (iii) implies that $B$ is an ellipsoid. Moreover, \Cref{pr.n=3} shows that it is enough to work with $n=3$.

\Cref{cor.iii_implies_sym} ensues that $B$ is centrally symmetric, and hence can be regarded as the unit ball of a Banach space $(\bR^3,\|\cdot\|_B)$. \Cref{le.istr_applies,th.istr} now ensure that the norm $\|\cdot\|_B$ is induced by an inner product on $\bR^3$, and hence the unit ball of $\|\cdot\|_B$ is an ellipsoid, as explained in \Cref{subse.ban}.
\end{proof_of_ellbis}

\section{Illuminated bodies and subspace lattice involutions}\label{se.illum}

The techniques employed in \Cref{se.park} will also help in improving on a characterization of ellipsoids via illumination by parallel rays due to Blaschke (\cite[pp. 157-159]{bla}). For the purpose of stating the result briefly, we introduce one more piece of terminology.

\begin{definition}\label{def.wbla}
  A convex body $B$ in $\bR^3$ has {\it the weak Blaschke property} (or {\it is weakly Blaschke}) if for any line $L$ there exists an affine plane $H\subset \bR^3$ such that $H$ intersects the interior of $B$ and all points in $H\cap \partial B$ lie on some $B$-supporting translate of $L$.

$B$ has {\it the weak dual Blaschke property} (or {\it is weakly dual Blaschke}) if for any linear plane $H\subset \bR^3$ there is a translate $H'$ of $H$ and a line $L\in \bR^3$ such that $H'$ intersects the interior of $B$ and all points in $H'\cap \partial B$ lie on some $B$-supporting translate of $L$.
\end{definition}

\begin{remark}\label{re.wbla}
  Note the difference to \Cref{def.bla}: in \Cref{def.wbla} we do not require the planes to pass through a given, fixed point (such as the origin in the case of the former definition). This justifies the adjective `weak' in the definition. In fact, it is easy to see that centrally symmetric weakly (dual) Blaschke convex bodies are automatically (dual) Blaschke.
\end{remark}

\begin{remark}
The weak Blaschke property can be interpreted as saying that if the body is illuminated with parallel rays, then the boundary of the shaded region contains a coplanar curve; hence the title of this section.
\end{remark}

With this language in place, Blaschke's result referred to above states that a smooth, strictly convex convex body in $\bR^3$ that is weakly Blaschke must be an ellipsoid. In this section we remove the smoothness and strict convexity requirements:

\begin{theorem}\label{th.illum}
  A convex $B\subset \bR^3$ that is weakly Blaschke is an ellipsoid.
\end{theorem}

We once more prepare the ground before giving the proof proper. First, we show that the weak Blaschke property entails its dual.

\begin{lemma}\label{le.self-dual}
  A weakly Blaschke convex body is weakly dual Blaschke.
\end{lemma}
\begin{proof}
  This is very similar in spirit to the proof of \Cref{pr.bla} from \Cref{pr.co_bla} via \Cref{le.psi_onto}, using the same multivalued Borsuk-Ulam theorem.

To any unit vector $v\in \bS^2$ associate the subset of the unit $2$-sphere consisting of those vectors $w$ with the property that $w\cdot v>0$ and $w$ is orthogonal to those planes associated as in the definition of the weak Blaschke property to the line containing $v$.

This gives rise, as in the discussion following \Cref{pr.bla}, to a map
\begin{equation*}
       \psi=:\bS^2\to\text{ closed convex subsets of }\bS^2.
\end{equation*}
$\psi$ can be shown to be upper semicontinuous and odd as in \Cref{le.odd}, and hence is onto according to \Cref{le.psi_onto}. This concludes the proof: we have just shown that {\it every} plane in $\bR^3$ is parallel to some affine plane associated to a line as in the definition of the weak Blaschke property.
\end{proof}

Next, we prove the existence of a center of symmetry.

\begin{lemma}\label{le.wbla_symm}
  A weakly Blaschke convex body $B\subset \bR^3$ has a center of symmetry.
\end{lemma}
\begin{proof}
  It suffices to show that the group of affine symmetries of $B$ is infinite. For this purpose, let $H\subset \bR^3$ be a plane  with the property that the two $B$-supporting translates of $H$ intersect $B$ at single point $p$ and $q$.

We now proceed very much as in the proof of \Cref{pr.aux_bis}. We denote by $K$ a planar section $(H+x)\cap B$ of $B$ that has non-empty relative interior. Then choose a line $L\subset H$. Just as in the cited proof, we can show using the weak Blaschke property that if $L_1$ and $L_2$ are the $K$-supporting translates of $L$, then the convex hull of $L_i\cap K$ contains the point $\overline{pq}\cap K$.

\Cref{le.supports} then ensures that $K$ has a center of symmetry at $\overline{pq}\cap K$. Since the $H$-parallel section $K$ of $B$ was arbitrary, we conclude that the affine involution of $\bR^3$ that preserves the points on the line $pq$ and acts as reflection across $pq\cap (H+x)$ along the plane $H+x$ preserves $B$.

The above construction provides us with infinitely many affine symmetries of $B$, one for each choice of plane $H$ with the generic property specified at the beginning.
\end{proof}

\begin{proof_of_illum}
  We now know from \Cref{le.self-dual} that $B$ is both weakly Blaschke and weakly dual Blaschke. Since moreover it has a center of symmetry by \Cref{le.wbla_symm}, we can assume that $B$ is centered at $0$ and hence is both Blaschke and dual Blaschke (cf. \Cref{re.wbla}). Running through the proof of \Cref{th.ell_bis}, this is sufficient to conclude that $B$ is an ellipsoid.
\end{proof_of_illum}

Metric Geometry (math.MG); Functional Analysis (math.FA)


Counting subgroups of fixed order in finite abelian groups

Abstract: We use recurrence relations to derive explicit formulas for counting the number of subgroups of given order (or index) in rank 3 finite abelian p-groups and use these to derive similar formulas in few cases for rank 4. As a consequence, we answer some questions by M. T$\ddot{a}$rn$\ddot{a}$uceanu in \cite{MT} and L. T$\dot{\acute{o}}$th in \cite{LT}. We also use other methods such as the method of fundamental group lattices introduced in \cite{MT} to derive a similar counting function in a special case of arbitrary rank finite abelian p-groups.

\section{Introduction}

The subject of counting various kinds of subgroups of finite abelian groups has a long and rich history. For instance, in the early 1900's, Miller\cite{MI04} determined the number of cyclic subgroups of prime power order in a finite abelian p-group $G$, where p is a prime number. In about the same time, Hilton\cite{HH} found a necessary and sufficient condition for the existence of subgroups of a given order and then gave a general procedure for counting subgroups of order $p^r$ in a finite abelian p-group of order $p^n$, for positive integers $r\leq n.$ For recent work on counting subgroups of finite abelian p-groups, see \cite{GB96}, \cite{GBJW}, \cite[Sections 6 \& 7]{GLP}, \cite{HT}, \cite{AI}, \cite{MI39}, \cite{MT} and \cite{LT}.

As an example, Hilton determined the total number of normal subgroups of index $p^2$ in any finite abelian p-group and noted that giving a general formula for the number of subgroups of order $p^r$ for every value of $r$ would be somewhat complicated. Progress has been made since then and there are now recursive formulas based on works of Stehling\cite{TS} and Birkhoff\cite{GB} that can compute the numbers for any $r.$ It can also be deduced using Hall's polynomials\cite{IGM} that the function that counts the number of subgroups is a polynomial in powers of $p.$ In addition, Shokuev\cite{VNS} finds an exact expression for the number of subgroups of any (possible) order of an arbitrary finite p-group. However, this requires subgroup structure of smaller rank subgroups in order to determine the expression. So far as we know, the explicit form of these counting polynomials has not been determined beyond that for rank 2 finite abelian p-groups. In \cite[Theorem 3]{GBJW01}, Bhowmik and Wu get a unimodality result on the coefficients and determine the leading coefficient of the polynomials in any rank. In this paper, we determine the explicit form of the counting polynomials in rank 3 and for some special values of $r$ in any rank. However, finding the general explicit polynomials in rank 4 or higher will be a tedious and difficult computation and we will illustrate this with examples.

By the fundamental theorem on the structure of finite abelian groups, any finite abelian group is isomorphic to a direct product of prime power order cyclic groups. Therefore, the number of subgroups in the finite abelian group is the product of the numbers of the subgroups contained in these cyclic groups\cite{RS}. This reduces the problem of counting subgroups in finite abelian groups to counting subgroups in finite abelian p-groups.

Let $\mathbb{Z}_{n}$ denote the cyclic group of order n. Given partitions $\lambda,\mu,\nu$ and $G$ an abelian $p$-group of type $\lambda  = (a_d,\dots,a_1),$ i.e., $G  \cong \mathbb{Z}/p^{a_1}\times \cdots \times \mathbb{Z}/p^{a_d},$ for some positive integers $a_1,\dots,a_d$ where $1\leq a_1\leq \dots \leq a_d,$ G has $g^\lambda_{\mu\nu}(p)$ subgroups $K$ such that $K$ has type $\mu$ and $G/K$ has type $\nu,$ where $g^\lambda_{\mu\nu}(X)\in \mathbb{Z}[X]$ is a Hall polynomial \cite{IGM}. Let
\begin{equation}
  \label{no_of_subgrps}
 h_b^\lambda (p) = \sum_{|\nu| = b} \sum_{\mu}g^\lambda_{\mu\nu} (p)
\end{equation}
\noindent
where $\nu = (\nu_l,\dots,\nu_1)$ with $\nu_1\leq \nu_2\leq \dots \leq \nu_l$ and $|\nu| = \nu_1+\dots +\nu_l.$

Then $h_b^\lambda(p) = h_b^{(a_d,\dots,a_1)}(p)$ is the number of subgroups of index $p^b$ (or equivalently of order $p^{m-b}$ where $m=\sum_{i=1}^d a_i$) in the rank $d$ abelian $p$-group $\mathbb{Z}/p^{a_1}\times \cdots \times \mathbb{Z}/p^{a_d}.$ Moreover, $h_b^\lambda(p)$ is a polynomial as it is a sum of Hall polynomials. When the finite abelian $p$-group has rank 2, the number of subgroups of order $p^b$ in $\mathbb{Z}/p^{a_1}\times \mathbb{Z}/p^{a_2}$ is given by \cite[Theorem 3.3]{MT} as follows (for an alternate proof, see the end of Section \ref{sec:convolution}):
\begin{equation}
  \label{rank2ratfns}
  h_b^{(a_2,a_1)}(p)=\begin{cases}
  \frac{p^{b+1}-1}{p-1}  & \text{if } 0\le b \leq a_1\\
  \frac{p^{a_1+1}-1}{p-1}  & \text{if } a_1\le b \leq a_2\\
  \frac{p^{a_1+a_2-b+1}-1}{p-1}  & \text{if } a_2\le b \leq a_1+a_2.\\
\end{cases}
\end{equation}
\noindent
Except for rank 2 finite abelian p-groups and some other specific cases, there was no general direct closed formula (i.e. explicit polynomial depending on the type $\lambda$ and $b$; expressed above as a rational function in each case) that computes $h_b^\lambda(p)$ given $\lambda$ and $b.$ We determine in the next section the polynomials that give the number of subgroups of given order in rank 3 and few cases of rank 4 and one particular case for arbitrary rank finite abelian $p$-group. This answers part of a problem posed by M. T$\ddot{a}$rn$\ddot{a}$uceanu\cite[Problem 5.1]{MT}. As an application of the rational function formulas, we will prove a conjecture posed by Laszlo T${\acute{o}}$th in \cite{LT}.
\section{Rank 3}
G. Birkhoff shows in \cite{GB} that $h_b^{(a_d,\dots,a_1)}$ is symmetric in $b$, that is, $h_b^{(a_d,\dots,a_1)} = h_{a_1+\dots+a_d-b}^{(a_d,\dots,a_1)}.$ So, for instance, in the three cases in equation (\ref{rank2ratfns}) for rank 2, the first case is symmetric with the third case while the second case is symmetric to itself. For example, when $d=2,$ the case $0\le b \leq a_1$ is symmetric to $a_2\le b \leq a_1+a_2$ because $$0\le b \leq a_1$$
$$-a_1\le -b \leq 0$$
$$a_2\leq a_1+a_2-b \leq a_1+a_2.$$
\noindent
Therefore, in order to determine the formula in the case $a_2\le b \leq a_1+a_2$ from the case when $0\le b \leq a_1$, we replace $b$ in the formula for $0\le b \leq a_1$ by $a_1+a_2-b.$ So, for instance in rank 2, it suffices to determine the compact form (as rational functions) of the polynomials in only two cases (first and second cases) and then use symmetry to deduce the remaining case (the last case).

Let $\lambda=(a_3,a_2, a_1)$ and $\lambda^{'}=(a_2, a_1)$. The partition of the interval $[0,a_1+a_2+a_3]$ as in  Hironaka's paper \cite{YH} leads to 10 cases for rank 3. For each of these cases, there is a polynomial $h_b^{(a_1,a_2,a_3)}(p)$ (easily described as a rational function) enumerating the number of subgroups of index $p^b$ (or order $p^{a_1+a_2+a_3-b}$) in $\mathbb{Z}/p^{a_1}\times \mathbb{Z}/p^{a_2}\times \mathbb{Z}/p^{a_3}.$

There are at least four methods of deducing the formulas. The first method is to use a lemma \cite[Lemma 2.3]{YH} that Y. Hironaka deduces from a recurrence relation of T. Stehling in \cite{TS}. This is the method used in the proof of the theorem below as her lemma shows how to express the formulas in higher rank in terms of those of lower rank which we already know and this requires the least amount of computation. The second method was suggested to the first author by Gautam Chinta (the first author's doctoral advisor) while working on a related problem on subgroup growth as such subgroup counting formulas help compute smaller rank cotype zeta functions as in \cite{CKK}. This method involves using a series of contour integrals, the residue theorem and convolution of generating series to arrive at the formulas. The third method is based on an extension of results in the second author's previous paper in \cite{PKASSS}. This will be explicated another time. The last method, called the method of fundamental group lattices, was introduced by M. T$\ddot{a}$rn$\ddot{a}$uceanu in \cite{MT} and we have been able to apply it to derive a compact formula for counting subgroups in case 1 of all ranks. We demo the second method of proof in a special case at the end of the paper.

\begin{theorem}
\label{rank3}
For every $0\leq b \leq a_1+a_2+a_3$, the number $h_b^{(a_3,a_2,a_1)}(p)$ of all subgroups of order $p^b$ (equivalently  of index $p^{a_1+a_2+a_3-b}$) in the finite abelian $p$-group $\mathbb{Z}/p^{a_1}\times \mathbb{Z}/p^{a_2}\times \mathbb{Z}/p^{a_3}$ where $1\leq a_1\leq a_2 \leq a_3,$ is given by one of the following  polynomials expressed as rational functions:\\
Case 1: $0\le b \leq a_1$
\begin{equation*}
   h_b^{(a_3,a_2,a_1)}(p) = \frac{p^{2b+3} - p^{b+2} - p^{b+1} +1}{ (p-1)(p^2-1)}.
\end{equation*}
Case 2:  $a_1\le b \leq a_2$
\begin{equation*}
h_b^{(a_3,a_2,a_1)}(p) = \frac{p^{b+a_1+3}+p^{b+a_1+2}-p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{ (p-1)(p^2-1)}.
\end{equation*}
Case 3: $a_2 < b \le a_3 \leq a_1 + a_2$
\begin{equation*}
h_b^{(a_3,a_2,a_1)}(p) =\frac{(b-a_{2}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(b-a_2)p^{a_1+a_2+1}-p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{(p-1)(p^{2}-1)}.
\end{equation*}
Case 4: $a_2 < b \le a_1 + a_2 \leq a_3$
$$h_b^{(a_3,a_2,a_1)}(p) = \frac{(b-a_{2}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(b-a_2)p^{a_1+a_2+1}-p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{(p-1)(p^{2}-1)}.$$
Case 5: $a_1 + a_2 \le b \le a_3$
$$h_b^{(a_3,a_2,a_1)}(p) = \frac{(a_{1}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-a_{1}p^{a_1+a_2+1}-p^{2a_1+2}-p^{a_{1}+a_{2}+2}-p^{a_{1}+a_{2}+1}+1}{(p-1)(p^{2}-1)}.$$
Case 6: $a_3 < b \leq a_{1}+a_{2} $
\begin{multline*}
  h_b^{(a_3,a_2,a_1)}(p) = \frac{(a_3-a_2+1)p^{a_2+a_1+3}+ 2p^{a_2+a_1+2}-(a_3-a_2-1)p^{a_2+a_1+1}-p^{a_3+a_2+a_1-b+2}}{(p-1)(p^2-1)}\\
  + \frac{-p^{a_3+a_2+a_1-b+1} -p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{(p-1)(p^2-1)}.
\end{multline*}
The remaining cases are symmetric to one of the cases above, i.e., they can be obtained when we replace $b$ by $a_1+a_2+a_3-b.$\\
Case 7: $a_1 + a_2 < a_3  < b \le  a_1 + a_3$
\begin{multline*}
 h_b^{(a_3,a_2,a_1)}(p) = \frac{(a_1+a_3-b+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(a_1+a_3-b)p^{a_1+a_2+1}-p^{2a_1+2}}{(p-1)(p^{2}-1)}\\
 + \frac{-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.
\end{multline*}
Case 8: $ a_3  < a_1 + a_2 < b \le  a_1 + a_3$
\begin{multline*}
h_b^{(a_3,a_2,a_1)}(p) = \frac{(a_1+a_3-b+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(a_1+a_3-b)p^{a_1+a_2+1}-p^{2a_1+2}}{(p-1)(p^{2}-1)}\\ + \frac{-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.
\end{multline*}
Case 9: $a_1 + a_3 \le b \le  a_2 + a_3$
$$h_b^{(a_3,a_2,a_1)}(p) = \frac{p^{2a_1+a_2+a_3-b+3}+p^{2a_1+a_2+a_3-b+2}-p^{2a_1+2}-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.$$
Case 10: $a_2 + a_3 \le b \le  a_1 + a_2 + a_3$
$$h_b^{(a_3,a_2,a_1)}(p) = \frac{p^{2a_1+2a_2+2a_3-2b+3}-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.$$
\end{theorem}
\begin{proofof}{Theorem \ref{rank3}} To simplify notation, $h_b^{(a_3,a_2,a_1)}(p)$ will be denoted by $N_{b}(\lambda)$ as in the notation used in \cite{YH}. By \cite[Lemma 2.3]{YH}, we have the recursive formula
\begin{equation}
    N_{b}(\lambda)=\sum_{i=0}^{b} p^{i}N_{i}(\lambda^{'}) - \sum_{|\lambda| +1-k}^{|\lambda^{'}|} p^{i}N_{i}(\lambda^{'}), 0\leq b\le |\lambda|
\end{equation}
\noindent where $\lambda = (a_d,\dots, a_1), \lambda^{'} = (a_{d-1},\dots, a_1), |\lambda| = a_d+\dots+a_1$ and the second summation appears only when $b > a_d.$ Thus, we compute each of the cases as follows: \\
Case 1: $0\le b \leq a_1$ \\
$N_{b}(\lambda)=\sum_{i=0}^{b} p^{i}N_{i}(\lambda^{'})=\sum_{i=0}^{b} p^{i} \frac{p^{i+1}-1}{p-1}=\sum_{i=0}^{b} \frac{p^{2i+1}-p^{i}}{p-1}=\frac{p^{2b+3}-p^{b+2}-p^{b+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 2: $a_1\le b \leq a_2$ \\
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{b} p^{i}N_{i}(\lambda^{'})=\frac{p^{2a_1+3}-p^{a_1+2}-p^{a_1+1}+1}{(p-1)(p^{2}-1)}+\sum_{i=a_{1}+1}^{b} p^{i} \frac{p^{a_{1}+1}-1}{p-1}$\\
$=\frac{p^{b+a_1+3}+p^{b+a_{1}+2}-p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 3: $a_2 < b \le a_3 \leq a_1 + a_2$\\
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{a_{2}} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{2}+1}^{b} p^{i}N_{i}(\lambda^{'})$\\
$=\frac{p^{a_{2}+a_1+3}+p^{a_{2}+a_{1}+2}-p^{2a_1+2}-p^{a_{2}+2}-p^{a_{2}+1}+1}{(p-1)(p^{2}-1)}+\sum_{i=a_{2}+1}^{b} p^{i}\frac{p^{a_{1}+a_{2}+1-i}-1}{p-1}$\\
$=\frac{(b-a_{2}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(b-a_2)p^{a_1+a_2+1}-p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 4: $a_2 < b \le a_1 + a_2 \leq a_3$ \\
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{a_{2}} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{2}+1}^{b} p^{i}N_{i}(\lambda^{'})$\\
$=\frac{p^{a_{2}+a_1+3}+p^{a_{2}+a_{1}+2}-p^{2a_1+2}-p^{a_{2}+2}-p^{a_{2}+1}+1}{(p-1)(p^{2}-1)}+\sum_{i=a_{2}+1}^{b} p^{i}\frac{p^{a_{1}+a_{2}+1-i}-1}{p-1}$\\
$=\frac{(b-a_{2}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(b-a_2)p^{a_1+a_2+1}-p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 5: $a_1 + a_2 \le b \le a_3$ \\
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{a_{2}} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{2}+1}^{a_{1}+a_{2}} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+a_{2}+1}^{b} p^{i}N_{i}(\lambda^{'})\\ - \sum_{i=a_{1}+a_{2}+a_{3}+1-b}^{a_{1}+a_{2}}p^{i}N_{i}(\lambda^{'})$\\
$=\frac{(a_{1}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-a_{1}p^{a_1+a_2+1}-p^{2a_1+2}-p^{a_{1}+a_{2}+2}-p^{a_{1}+a_{2}+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 6: $a_3 < b \leq a_{1}+a_{2} $ \\
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{a_{2}} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{2}+1}^{a_{3}} p^{i}N_{i}(\lambda^{'})+\sum_{i=(+1}^{b} p^{i}N_{i}(\lambda^{'})\\ - \sum_{i=a_{1}+a_{2}+a3+1-b}^{a_{1}+a_{2}}p^{i}N_{i}(\lambda^{'})$\\
$=\frac{(a_{3}-a_{2}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(a_{3}-a_2)p^{a_1+a_2+1}-p^{2a_1+2}-p^{a_{3}+2}-p^{a_{3}+1}+1}{(p-1)(p^{2}-1)}+\sum_{i=a_{3}+1}^{b} p^{i}\frac{p^{a_{1}+a_{2}+1-i}-1}{p-1}$\\
$-\sum_{i=a_{1}+a_{2}+a_{3}+1-b}^{a_{1}+a_{2}}p^{i}\frac{p^{a_{1}+a_{2}+1-i}-1}{p-1}$\\
$=\frac{(a_3-a_2+1)p^{a_2+a_1+3}+ 2p^{a_2+a_1+2}-(a_3-a_2-1)p^{a_2+a_1+1}-p^{a_3+a_2+a_1-b+2}-p^{a_3+a_2+a_1-b+1}}{(p-1)(p^2-1)} \\
+ \frac{-p^{2a_1+2}-p^{b+2}-p^{b+1}+1}{(p-1)(p^2-1)}. $\\
The remaining cases are symmetric to one of the cases above, i.e., they can be obtained when we replace $b$ by $a_1+a_2+a_3-b.$\\
Case 7: $a_1 + a_2 < a_3  < b \le  a_1 + a_3$ Replace $b$ in case 4 by $a_1+a_2+a_3-b.$\\
$N_{b}(\lambda)= \frac{(a_1+a_3-b+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(a_1+a_3-b)p^{a_1+a_2+1}-p^{2a_1+2}-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 8: $ a_3  < a_1 + a_2 < b \le  a_1 + a_3$ Replace $b$ in case 3 by $a_1+a_2+a_3-b.$\\
$N_{b}(\lambda) = \frac{(a_1+a_3-b+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(a_1+a_3-b)p^{a_1+a_2+1}-p^{2a_1+2}-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 9: $a_1 + a_3 \le b \le  a_2 + a_3$ Replace $b$ in case 2 by $a_1+a_2+a_3-b.$\\
$N_{b}(\lambda) = \frac{p^{2a_1+a_2+a_3-b+3}+p^{2a_1+a_2+a_3-b+2}-p^{2a_1+2}-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.$\\
Case 10: $a_2 + a_3 \le b \le  a_1 + a_2 + a_3$ Replace $b$ in case 1 by $a_1+a_2+a_3-b.$\\
$N_{b}(\lambda) = \frac{p^{2a_1+2a_2+2a_3-2b+3}-p^{a_1+a_2+a_3-b+2}-p^{a_1+a_2+a_3-b+1}+1}{(p-1)(p^{2}-1)}.$
\end{proofof}

\begin{remark} We can deduce all other cases using only case 6 as follows:\\
Case 1: replace $a_3,a_2$ and $a_1$ in  case 6 formula by $b.$\\
Case 2: replace $a_3$ and $a_2$ in  case 6 formula by $b.$\\
Cases 3 and 4: replace $a_3$ in  case 6 formula by $b.$\\
Case 5: replace $a_3$ and $b$  in  case 6 formula  by $a_2+a_1$ \\
Cases 7 and 8: replace $a_3$ in  case 6 formula by $a_1+a_2+a_3-b.$\\
Case 9: replace $a_3$ and $a_2$ in  case 6 formula by $a_1+a_2+a_3-b.$\\
Case 10: replace $a_3,a_2$ and $a_1$ in  case 6 formula by $a_1+a_2+a_3-b.$
\end{remark}
\section{Rank 4}
\subsection{Number of subgroups of given order}
While it was a somewhat lengthy but accomplishable task to determine the cases for rank 3, there does not appear to be a systematic way of determining all the cases for rank 4. We were able to list more than 22 cases some of which we list below. These formulas are also proved in exactly the same way as in Theorem \ref{rank3}. Here $\lambda=(a_4,a_3,a_2, a_1)$ and $\lambda^{'}=(a_3,a_2, a_1).$\\
When $0\le b \leq a_1,$ we get \\
$N_{b}(\lambda)=\sum_{i=0}^{b} p^{i}N_{i}(\lambda^{'})=\sum_{i=0}^{b} p^{i} \frac{p^{2i+3}-p^{i+2}-p^{i+1}+1}{(p-1)(p^{2}-1)}=\sum_{i=0}^{b}  \frac{p^{3i+3}-p^{2i+2}-p^{2i+1}+p^{i}}{(p-1)(p^{2}-1)}$\\
$=\frac{p^{3b+6}-p^{2b+5}-p^{2b+4}-p^{2b+3}+p^{b+3}+p^{b+2}+p^{b+1}-1}{(p-1)(p^{2}-1)(p^{3}-1)}.$\\
When $a_1\le b \leq a_2,$ we get
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{b} p^{i}N_{i}(\lambda^{'})$\\
$=\frac{p^{3a_{1}+6}-p^{2a_{1}+5}-p^{2a_{1}+4}-p^{2a_{1}+3}+p^{a_{1}+3}+p^{a_{1}+2}+p^{a_{1}+1}-1}{(p-1)(p^{2}-1)(p^{3}-1)}+\sum_{i=a_{1}+1}^{b} p^{i} \frac{p^{i+a_1+3}+p^{i+a_{1}+2}-p^{2a_1+2}-p^{i+2}-p^{i+1}+1}{(p-1)(p^{2}-1)}$\\
$=\frac{p^{2b+a_{1}+6}+p^{2b+a_{1}+5}+p^{2b+a_{1}+4}-p^{2a_{1}+b+5}-p^{2a_{1}+b+4}-p^{2a_{1}+b+3}+p^{3a_{1}+3}-p^{2b+5}-p^{2b+4}-p^{2b+3}+p^{b+3}+p^{b+2}+p^{b+1}-1}{(p-1)(p^{2}-1)(p^{3}-1)}.$\\
And as a final example, when $a_2\le b \leq min \{a_3, a_{1}+a_2 \},$ we have \\
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{a_{2}} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{2}+1}^{b} p^{i}N_{i}(\lambda^{'})-\sum_{i=a_{1}+a_{2}+a_{3}+a_{4}+1-b}^{a_{1}+a_{2}+a_{3}}p^{i}N_{i}(\lambda^{'})$\\
$=\frac{p^{2a_{2}+a_{1}+6}+p^{2a_{2}+a_{1}+5}+p^{2a_{2}+a_{1}+4}-p^{2a_{1}+a_{2}+5}-p^{2a_{1}+a_{2}+4}-p^{2a_{1}+a_{2}+3}+p^{3a_{1}+3}}{(p-1)(p^{2}-1)(p^{3}-1)}\\
+\frac{-p^{2a_{2}+5}-p^{2a_{2}+4}-p^{2a_{2}+3}+p^{a_{2}+3}+p^{a_{2}+2}+p^{a_{2}+1}-1}{(p-1)(p^{2}-1)(p^{3}-1)}\\
+\sum_{i=a_{2}+1}^{b} p^{i} \frac{(i-a_{2}+1)p^{a_2+a_1+3}+p^{a_2+a_{1}+2}-(i-a_2)p^{a_1+a_2+1}-p^{2a_1+2}-p^{i+2}-p^{i+1}+1}{(p-1)(p^{2}-1)}$\\
$=\frac{(b+1-a_2)p^{a_{2}+a_{1}+b+6}+(b+1-a_{2})p^{a_{2}+a_{1}+b+5}+(a_{2}-b-1)p^{a_{2}+a_{1}+b+3}+(a_{2}-b-1)p^{a_{1}+a_{2}+b+2}}{(p-1)(p^{2}-1)(p^{3}-1)}\\
+\frac{p^{2a_{2}+a_{1}+4}+p^{a_{1}+2a_{2}+3}+p^{a_{1}+2a_{2}+2}}{(p-1)(p^{2}-1)(p^{3}-1)}$\\
$+\frac{-p^{2a_{1}+b+5}-p^{2a_{1}+b+4}-p^{2a_{1}+b+3}+p^{b+3}+p^{b+2}+p^{b+1}+p^{3a_{1}+3}-p^{2b+5}-p^{2b+4}-p^{2b+3}-1}{(p-1)(p^{2}-1)(p^{3}-1)}.$
\subsection{Total number of subgroups}
Let $\lambda = (a_4,a_3,a_2,a_1)$ and let $N(\lambda)$ denote the total number of subgroups of $\mathbb{Z}/p^{a_1}\times \mathbb{Z}/p^{a_2}\times \mathbb{Z}/p^{a_3}\times \mathbb{Z}/p^{a_4}$ where $1\le a_1\le a_2\le a_3\le a_4,$ and $n= a_1+a_2+a_3+a_4.$ Then
\begin{equation*}
N(\lambda) = \sum_{b=0}^n N_b(\lambda).
\end{equation*}

A conjecture of L. T$\acute{o}$th in \cite[Conjecutre 10]{LT} claims that the degree of the polynomial $N(m,m,m,m)$, i.e. when $\lambda = (m,m,m,m),$ is $4m$ and that it has leading coefficient of 1. We prove that this follows from our formulas in Theorem \ref{rank3}. First, let's state a corollary for the special case of a rank 3 abelian p-group of type $(m,m,m).$ The proof of the corollary follows immediately from setting $a_1=a_2=a_3=m$ in the theorem.

\begin{corollary}\label{mmm}
For every $0\leq b \leq 3m$, the number $h_b^{(m,m,m)}(p)$ of all subgroups of order $p^b$ in the finite abelian $p$-group $\mathbb{Z}/p^{m}\times \mathbb{Z}/p^{m}\times \mathbb{Z}/p^{m}$ where $1\leq m,$ is given by one of the following polynomials expressed as rational functions:\\
Case 1: $0\le b \leq m$
$$h_b^{(m,m,m)}(p) = \frac{p^{2b+3} - p^{b+2} - p^{b+1} +1}{(p-1)(p^2-1)}.$$
Case 2:  $m \le b \leq 2m$
$$h_b^{(m,m,m)}(p) = \frac{p^{2m+3}+p^{2m+2}+p^{2m+1}-p^{3m+2-b}-p^{3m+1-b}-p^{b+2}-p^{b+1}+1}{ (p-1)(p^2-1)}.$$
Case 3: $2m \le b \le 3m$
$$h_b^{(m,m,m)}(p) =\frac{p^{6m+3-2b}-p^{3m+2-b}-p^{3m+1-b}+1}{(p-1)(p^{2}-1)}.$$
\end{corollary}
\begin{theorem}\label{Conj10}
The the total number of subgroups $N(m,m,m,m)$ of the rank 4 finite abelian p-group $\mathbb{Z}/p^{m}\times \mathbb{Z}/p^{m}\times \mathbb{Z}/p^{m}\times \mathbb{Z}/p^{m}$ where $1\le m$ is given by a polynomial expressed as a rational function as below:
\begin{multline*}
 N(m,m,m,m) = \sum_{b=1}^{4m} N_b(m,m,m,m)\\
 = \frac{ {\left(p^{2} + p + 1\right)}^{3} {\left(p^{2} + 1\right)} p^{4 \, m + 2} - {\left({\left(2 \, m + 3\right)} p^{3} - 2 \, m - 1\right)} {\left(p^{3} + p\right)} {\left(p + 1\right)}^{3} p^{3 \, m} }{{\left(p^{2} - 1\right)}^{2} {\left(p^3 - 1\right)}^{2}}\\
 - \frac{4 \, m p^{4} + 4 \, m p^{3} + 7 \, p^{4} + 9 \, p^{3} - 4 \, m p + 6 \, p^{2} - 4 \, m + p - 1}{{\left(p^{2} - 1\right)}^{2} {\left(p^3 - 1\right)}^{2}}.
\end{multline*}
In particular, the degree of the polynomial is $4m$ and its leading coefficient is 1.

\end{theorem}

\begin{proofof}{Theorem \ref{Conj10}}
We break the interval $[0,4m]$ into 4 sub-intervals and use \cite[Lemma 3.2]{YH} and the corollary above in each case to derive a polynomial expression for the number of subgroups of order $p^b$ in a rank 4 abelian p-group of type $(m,m,m,m)$. Summing over these polynomials will result in the required polynomial expression.\\
Case 1: $0\le b\le m$
\begin{align*}
N_b(m,m,m,m) &= \sum_{i=0}^{b} p^{i}N_{i}(m,m,m)=\sum_{i=0}^{b} p^{i} \frac{p^{2i+3} - p^{i+2} - p^{i+1} +1}{(p-1)(p^2-1)}\\
&= \frac{p^{3 \, b + 6} - p^{2 \, b + 5} - p^{2 \, b + 4} - p^{2 \, b + 3} + p^{b + 3} + p^{b + 2} + p^{b + 1} - 1}{{\left(p - 1\right)} {\left(p^2 - 1\right)} {\left(p^3 - 1\right)}}.
\end{align*}
Case 2: $m < b\le 2m$
\begin{align*}
N_b(m,m,m,m) &= \sum_{i=0}^{m} p^{i}N_{i}(m,m,m) + \sum_{i=m+1}^{b} p^{i}N_{i}(m,m,m) - \sum_{i=4m+1-b}^{3m} p^{i}N_{i}(m,m,m)\\
&=  \sum_{i=0}^{m} p^{i}\frac{p^{2i+3} - p^{i+2} - p^{i+1} +1}{(p-1)(p^2-1)}\\ &+ \sum_{i=m+1}^{b} p^{i}\frac{p^{2m+3}+p^{2m+2}+p^{2m+1}-p^{3m+2-i}-p^{3m+1-i}-p^{i+2}-p^{i+1}+1}{ (p-1)(p^2-1)}\\
&- \sum_{i=4m+1-b}^{3m} p^{i}\frac{p^{6m+3-2i}-p^{3m+2-i}-p^{3m+1-i}+1}{(p-1)(p^{2}-1)}\\
&= -\frac{p^{2 \, b + 5} + p^{2 \, b + 4} + p^{2 \, b + 3} - p^{b + 2 \, m + 6} - p^{b + 2 \, m + 5} - 2 \, p^{b + 2 \, m + 4}}{(p-1)(p^2-1)(p^3-1)} \\
&+\frac{ - p^{b + 2 \, m + 3} - p^{b + 2 \, m + 2} - p^{b + 3} - p^{b + 2}}{(p-1)(p^2-1)(p^3-1)} \\
&-\frac{ - p^{b + 1} - p^{-b + 4 \, m + 3} - p^{-b + 4 \, m + 2} - p^{-b + 4 \, m + 1} + p^{3 \, m + 5} + 2 \, p^{3 \, m + 4}}{(p-1)(p^2-1)(p^3-1)}\\
&+ \frac{2 \, p^{3 \, m + 3} + 2 \, p^{3 \, m + 2} + p^{3 \, m + 1} + 1}{(p-1)(p^2-1)(p^3-1)}.
\end{align*}
Case 3: $2m < b \le 3m$
\begin{align*}
N_b(m,m,m,m) &= \sum_{i=0}^{m} p^{i}N_{i}(m,m,m) + \sum_{i=m+1}^{2m} p^{i}N_{i}(m,m,m) + \sum_{i=2m+1}^{b} p^{i}N_{i}(m,m,m)\\ &- \sum_{i=4m+1-b}^{2m} p^{i}N_{i}(m,m,m) -\sum_{i=2m+1}^{3m} p^{i}N_{i}(m,m,m)\\
&=  \sum_{i=0}^{m} p^{i}\frac{p^{2i+3} - p^{i+2} - p^{i+1} +1}{(p-1)(p^2-1)}\\
&+ \sum_{i=m+1}^{2m} p^{i}\frac{p^{2m+3}+p^{2m+2}+p^{2m+1}-p^{3m+2-i}-p^{3m+1-i}-p^{i+2}-p^{i+1}+1}{ (p-1)(p^2-1)}\\
&+ \sum_{i=2m+1}^{b} p^{i}\frac{p^{6m+3-2i}-p^{3m+2-i}-p^{3m+1-i}+1}{(p-1)(p^{2}-1)}\\
&- \sum_{i=4m+1-b}^{2m} p^{i}\frac{p^{2m+3}+p^{2m+2}+p^{2m+1}-p^{3m+2-i}-p^{3m+1-i}-p^{i+2}-p^{i+1}+1}{ (p-1)(p^2-1)}\\ &- \sum_{i=2m+1}^{3m} p^{i}\frac{p^{6m+3-2i}-p^{3m+2-i}-p^{3m+1-i}+1}{(p-1)(p^{2}-1)}\\
&=\frac{p^{b + 3} + p^{b + 2} + p^{b + 1} + p^{-b + 6 \, m + 6} + p^{-b + 6 \, m + 5} + 2 \, p^{-b + 6 \, m + 4}}{(p-1)(p^2-1)(p^3-1)}\\ &+ \frac{p^{-b + 6 \, m + 3} + p^{-b + 6 \, m + 2} + p^{-b + 4 \, m + 3}}{(p-1)(p^2-1)(p^3-1)}\\
&+ \frac{p^{-b + 4 \, m + 2} + p^{-b + 4 \, m + 1} - p^{-2 \, b + 8 \, m + 5} - p^{-2 \, b + 8 \, m + 4} - p^{-2 \, b + 8 \, m + 3} - p^{3 \, m + 5} - 2 \, p^{3 \, m + 4}}{(p-1)(p^2-1)(p^3-1)}\\
&+ \frac{- 2 \, p^{3 \, m + 3} - 2 \, p^{3 \, m + 2} - p^{3 \, m + 1} - 1}{(p-1)(p^2-1)(p^3-1)}.
\end{align*}
Case 4: $3m < b\le 4m$
\begin{align*}
N_b(m,m,m,m) &= \sum_{i=0}^{m} p^{i}N_{i}(m,m,m) + \sum_{i=m+1}^{2m} p^{i}N_{i}(m,m,m) + \sum_{i=2m+1}^{3m} p^{i}N_{i}(m,m,m) \\
&- \sum_{i=4m+1-b}^{m} p^{i}N_{i}(m,m,m) -\sum_{i=m+1}^{2m} p^{i}N_{i}(m,m,m) - \sum_{i=2m+1}^{3m} p^{i}N_{i}(m,m,m)\\
&=  \sum_{i=0}^{m} p^{i}\frac{p^{2i+3} - p^{i+2} - p^{i+1} +1}{(p-1)(p^2-1)}\\
&+ \sum_{i=m+1}^{b} p^{i}\frac{p^{2m+3}+p^{2m+2}+p^{2m+1}-p^{3m+2-i}-p^{3m+1-i}-p^{i+2}-p^{i+1}+1}{ (p-1)(p^2-1)}
\end{align*}
\begin{align*}
&+ \sum_{i=2m+1}^{3m} p^{i}\frac{p^{6m+3-2i}-p^{3m+2-i}-p^{3m+1-i}+1}{(p-1)(p^{2}-1)}\\
&- \sum_{i=4m+1-b}^{m} p^{i}\frac{p^{2i+3} - p^{i+2} - p^{i+1} +1}{(p-1)(p^2-1)}\\
&- \sum_{i=m+1}^{2m} p^{i}\frac{p^{2m+3}+p^{2m+2}+p^{2m+1}-p^{3m+2-i}-p^{3m+1-i}-p^{i+2}-p^{i+1}+1}{ (p-1)(p^2-1)}\\
&- \sum_{i=2m+1}^{3m} p^{i}\frac{p^{6m+3-2i}-p^{3m+2-i}-p^{3m+1-i}+1}{(p-1)(p^{2}-1)}\\
&=\frac{p^{-b + 4 \, m + 3} + p^{-b + 4 \, m + 2} + p^{-b + 4 \, m + 1} - p^{-2 \, b + 8 \, m + 5}}{(p-1)(p^2-1)(p^3-1)}\\
&+ \frac{- p^{-2 \, b + 8 \, m + 4} - p^{-2 \, b + 8 \, m + 3} + p^{-3 \, b + 12 \, m + 6} - 1}{(p-1)(p^2-1)(p^3-1)}.
\end{align*}
Now summing over the 4 cases, we get that
\begin{align*}
N(m,m,m,m)&=\sum_{b=0}^{4m} N_b(m,m,m,m)\\
&= \frac{ {\left(p^{2} + p + 1\right)}^{3} {\left(p^{2} + 1\right)} p^{4 \, m + 2} - {\left({\left(2 \, m + 3\right)} p^{3} - 2 \, m - 1\right)} {\left(p^{3} + p\right)} {\left(p + 1\right)}^{3} p^{3 \, m} }{{\left(p^{2} - 1\right)}^{2} {\left(p^3 - 1\right)}^{2}}\\
&- \frac{4 \, m p^{4} + 4 \, m p^{3} + 7 \, p^{4} + 9 \, p^{3} - 4 \, m p + 6 \, p^{2} - 4 \, m + p - 1}{{\left(p^{2} - 1\right)}^{2} {\left(p^3 - 1\right)}^{2}}.
\end{align*}
Now, we can see that the degree of the polynomial is $4m$ and the leading coefficient is 1, as conjectured by L. T$\acute{o}$th.
\end{proofof}
\begin{remark} In an unpublished paper \cite{CCL}, Chew, Chin, and Lim derive an explicit formula for the number of subgroups of a finite abelian p-group of rank 4. Here we will use the formula to reprove the above theorem and also prove one more conjecture by T$\acute{o}$th in \cite[Conjecture 9]{LT}.
\begin{theorem} (Chew, Chin, and Lim) Let $1 \le w \le x \le y \le z$. The number of subgroups of $\mathbb{Z}/p^{w}\times \mathbb{Z}/p^{x}\times \mathbb{Z}/p^{y}\times \mathbb{Z}/p^{z}$ is
\begin{align*}
N(w,x,y,z) &= \sum_{i=0}^{w-1}\sum_{j=0}^{i} [(w+x+y+z-4i+1)(2i-2j+1)p^{3i+j}\\
&+(w+x+y+z-4i-1)(2i-2j+1)p^{3i+j+1}\\
&+ 2(w+x+y+z-4i-2)(i-j+1)p^{3i+j+2}]\\
&+\sum_{i=0}^{w}\sum_{j=w}^{x-1} (w+j-2i+1)[(x+y+z-3j+1)p^{w+2j+i}\\
\end{align*}
\begin{align*}
&+(x+y+z-3j-1)p^{w+2j+i+1}] \\
&+\sum_{i=0}^{w}\sum_{j=x}^{y} (y+z-2j+1)(w+x+-2i+1)p^{w+x+i+j}.
\end{align*}
\end{theorem}
Since all coefficients are positive, none of the terms will cancel out and the degree of the polynomial in $p$ can be easily determined. By comparing the exponents, we can see that the highest degree appears in the last double sum. Setting $i=w$ and $j=y,$ we get that the leading term of $N(w,x,y,z)$ is $(z-y+1)(x-w+1)p^{2w+x+y},$ with degree $2w+x+y,$ confirming conjecture 9 in \cite{LT}.

Moreover, setting $w=x=y=z=m,$ in the above leading term, we get that the degree of the leading term in $N(m,m,m,m)$ is $4m,$ confirming Conjecture 10 in the same paper.
\end{remark}

\begin{remark} What is new in Theorem \ref{Conj10} is the explicit polynomial. Else, both conjectures 9 and 10 of Toth follow from Theorem 3 in \cite{GBJW01} as they have determined the leading coefficient and degree of the counting polynomial for all ranks.
\end{remark}

\section{A special case of any rank}
The recursive counting of the subgroups used in the proof of Theorem \ref{rank3} based on \cite[Lemma 2.3]{YH} enabled us to deduce a compact rational form for the polynomials in all the cases for ranks 2 and 3. Here we will demo the method of fundamental group lattices in one case for all ranks. We remark that it is difficult to use this method to give an explicit formula in the remaining cases (those that are not symmetric to Case 1) of rank 3 or higher.

Below we extend T$\ddot{a}$rn$\ddot{a}$uceanu's result for counting subgroups as in \cite[Theorem 3.3]{MT} for case 1 to any rank.
\begin{theorem}
\label{funlatt} Let $0\leq b \leq a_1.$ Then
\begin{equation}
h_b^{(a_k,\dots,a_1)}(p) = \prod_{i=2}^k\frac{p^{b+i-1}-1}{p^{i-1}-1}.
\end{equation} Similarly, if $a_2+\dots+a_k\le b \leq a_1+\dots+a_k$, then replace $b$ in the above formula by $a_1+\dots+a_k-b$ to get
\begin{equation}
h_b^{(a_k,\dots,a_1)}(p) = \prod_{i=2}^k\frac{p^{a_1+\dots+a_k-b+i-1}-1}{p^{i-1}-1}.  \end{equation}
\end{theorem}
\begin{proofof}{Theorem \ref{funlatt}} Following the notation in \cite{MT}, let $A= (a_{ij})$ be a solution of $(\ast)$ below corresponding to a subgroup of order $p^{a_1+\dots+a_k-b}$ in $\mathbb{Z}/p^{a_1}\times \mathbb{Z}/p^{a_2}\times \cdots \times \mathbb{Z}/p^{a_k}.$ \\

$$(\ast) \begin{cases} \text{i) } a_{ij} = 0 \text{ for any } i>j\\
\text{ii) } 0\leq a_{1j},a_{2j},\dots,a_{j-1,j} < a_{jj} \text{ for any } j \in \{1,2,\dots, k\}\\
\text{iii) } \text{1) } a_{11}|p^{a_1}\\
\hspace{0.2in}\text{ 2) } a_{22}|(p^{a_2},p^{a_1}\frac{a_{12}}{a_{11}})\\
\hspace{0.2in}\text{ 3) } a_{33}|(p^{a_3},p^{a_2}\frac{a_{23}}{a_{22}},p^{a_1}\frac{\begin{vmatrix}
a_{12} &a_{13} \\
a_{22} & a_{23}
\end{vmatrix}}{a_{22}a_{11}})\\
\vspace{0.1in}\\
\hspace{1in} \longvdots{1.5em}\\
\hspace{0.2in}\text{ k) } a_{kk}|(p^{a_k},p^{a_{k-1}}\frac{a_{k-1,k}}{a_{k-1,k-1}},p^{a_{k-2}}\frac{\begin{vmatrix}
a_{k-2,k-1} &a_{k-2,k} \\
a_{k-1,k-1} & a_{k-1,k}
\end{vmatrix}}{a_{k-1,k-1}a_{k-2,k-2}},\dots,p^{a_1}\frac{\begin{vmatrix}
a_{12} & a_{13} & \cdots & a_{1,k} \\
a_{22} & a_{23} & \cdots & a_{2,k} \\
\textbf{.}  & \textbf{.}  & & \textbf{.}  \\
\textbf{.}  & \textbf{.}  & & \textbf{.}  \\
\textbf{.}  & \textbf{.}  & & \textbf{.}  \\
a_{k-1,2} & a_{k-1,3} & \cdots & a_{k-1,k}
\end{vmatrix}}{a_{k-1,k-1}a_{k-2,k-2}\dots a_{11}}).
\end{cases} $$

Put $a_{11}=p^{i_1}$ where $0\leq i_1\leq a_1.$ By the first remark on pp. 376 of \cite{MT}, the order of a subgroup corresponding to $A= (a_{ij})$ is $$\frac{p^{\sum_{i=1}^k a_i}}{\prod_{i=1}^k a_{ii}} = p^{a_1+\dots+a_k-b}$$
Therefore, $a_{11}\cdots a_{kk}=p^b.$ Since $a_{11}=p^{i_1},$ we have $a_{22}\cdots a_{kk}=p^{b-i_1}.$ Let\\

\hspace{2in} $a_{22} = p^{i_2}$ where $0\leq i_2 \leq  b-i_1,$ \\

\hspace{2in} $a_{33}=p^{i_3}$ where $0\leq i_3\leq b-i_1-i_2,$\\

\hspace{2in} $\dots $ \\

\hspace{2in} $a_{k-1,k-1}=p^{i_{k-1}}$ where $0\leq i_{k-1} \leq b-i_1-\dots - i_{k-2}.$ \\
Then $a_{kk}=p^{b-i_1-\dots - i_{k-1}}.$ Now, the conditions become
$$p^{i_2}|(p^{a_2},p^{a_1-i_1}a_{12}) = p^{a_1-i_1}(p^{a_2-a_1+i_1},a_{12})$$
$$p^{i_3}|(p^{a_3},p^{a_2-i_2}a_{12}, p^{a_1-i_1-i_2}(a_{12}a_{23}-p^{i_2}a_{13})$$
$$\dots$$
$$p^{b-i_1-\dots - i_{k-1}}|(p^{a_k},p^{a_{k-1}}\frac{a_{k-1,k}}{a_{k-1,k-1}},p^{a_{k-2}}\frac{\begin{vmatrix}
a_{k-2,k-1} &a_{k-2,k} \\
a_{k-1,k-1} & a_{k-1,k}
\end{vmatrix}}{a_{k-1,k-1}a_{k-2,k-2}},\dots,p^{a_1}\frac{\begin{vmatrix}
a_{12} & a_{13} & \cdots & a_{1,k} \\
a_{22} & a_{23} & \cdots & a_{2,k} \\
\textbf{.}  & \textbf{.}  & & \textbf{.}  \\
\textbf{.}  & \textbf{.}  & & \textbf{.}  \\
\textbf{.}  & \textbf{.}  & & \textbf{.}  \\
a_{k-1,2} & a_{k-1,3} & \cdots & a_{k-1,k}
\end{vmatrix}}{a_{k-1,k-1}a_{k-2,k-2}\dots a_{11}}).$$
These conditions are satisfied by all
$$a_{12}< a_{22}=p^{i_2},$$
$$a_{13},a_{23} < a_{33}=p^{i_3},$$
$$\cdots$$
$$a_{1,k-1},\dots,a_{k-2,k-1}< a_{k-1,k-1}=p^{i_{k-1}}$$
$$a_{1,k},\dots,a_{k-1,k}< a_{k,k}=p^{b-i_1-\cdots-i_{k-1}}.$$
So, we get $$p^{i_2}\cdot (p^{i_3})^2 \cdots (p^{i_{k-1}})^{k-2}\cdot (p^{b-i_1-\cdots-i_{k-1}})^{k-1} = p^{(k-1)b-(k-1)i_1-(k-2)i_2-\cdots -2i_{k-2}-i_{k-1}}$$ distinct solutions of $(\ast)$ for given $a_{11}=p^{i_1},a_{22}=p^{i_2},\cdots,a_{k-1,k-1}=p^{i_{k-1}},$ and  $a_{k,k}=p^{b-i_1-\cdots-i_{k-1}}.$ The number of subgroups of index $p^b$ (or order $p^{a_1+\dots+a_k-b}$) in $\mathbb{Z}/p^{a_1}\times \cdots \times \mathbb{Z}/p^{a_k}$ is then
\begin{equation}
h_b^{(a_k,\dots,a_1)}(p) = \sum_{i_1=0}^{b}\sum_{i_2=0}^{b-i_1}\cdots \sum_{i_{k-1}=0}^{b-i_1-\cdots-i_{k-2}} p^{(k-1)b-(k-1)i_1-(k-2)i_2-\cdots -2i_{k-2}-i_{k-1}}.
\end{equation}

Now we will use induction to prove that
\begin{equation}\label{anyrankclaim}
\sum_{i_1=0}^{b}\sum_{i_2=0}^{b-i_1}\cdots \sum_{i_{k-1}=0}^{b-i_1-\cdots-i_{k-2}} p^{(k-1)b-(k-1)i_1-(k-2)i_2-\cdots -2i_{k-2}-i_{k-1}} = \prod_{i=2}^k\frac{p^{b+i-1}-1}{p^{i-1}-1}.
\end{equation}

By Theorem \ref{rank3} and previous discussion, equation (\ref{anyrankclaim}) is true for when $k=2,3.$ Let us assume that it is true for rank k. Now consider the left hand side of equation (\ref{anyrankclaim})  for rank $k+1,$ that is,
$$\sum_{i_1=0}^{b}\sum_{i_2=0}^{b-i_1}\cdots \sum_{i_{k-1}=0}^{b-i_1-\cdots-i_{k-2}} \sum_{i_{k}=0}^{b-i_1-\cdots-i_{k-2}-i_{k-1}}p^{(k)b-(k)i_1-(k-1)i_2-\cdots -3i_{k-2}-2i_{k-1}-i_{k}}  $$\\
and put $\alpha=b-i_1$. Then \\
$$\sum_{i_1=0}^{b}p^{b-i_1}\sum_{i_2=0}^{\alpha}\cdots \sum_{i_{k-1}=0}^{\alpha-\cdots-i_{k-2}} \sum_{i_{k}=0}^{\alpha-\cdots-i_{k-2}-i_{k-1}}p^{(k-1)\alpha-(k-1)i_2-\cdots -3i_{k-2}-2i_{k-1}-i_{k} }$$
\begin{align*}
&= \sum_{i_1=0}^{b}p^{b-i_1} \prod_{i=3}^{k+1}\frac{p^{b+i-1}-1}{p^{i-1}-1}\\
&= \prod_{i=3}^{k+1}\frac{p^{b+i-1}-1}{p^{i-1}-1}\sum_{i_1=0}^{b}p^{b-i_1}\\
&= \prod_{i=2}^{k+1}\frac{p^{b+i-1}-1}{p^{i-1}-1}.
\end{align*} This completes the proof.
\end{proofof}
\begin{remark} An application of the method used in Theorem \ref{rank3} also shows that we can rederive the above result in a simpler way as follows. Let $\lambda=(a_k,\dots, a_1), \lambda^{'}=(a_{k-1},\dots, a_1),$ and  $0\le b \leq a_1.$ Then a repeated application of the recurrence (\ref{Stehrec}) for $N_b(\lambda)$ shows that
\begin{align*}
N_{b}(\lambda) &= \sum_{i_1=0}^{b} p^{i_1}N_{i_1}(\lambda^{'})\\
 &=\sum_{i_1=0}^{b} p^{i_1}\sum_{i_2=0}^{i_1} p^{i_2}N_{i_2}(\lambda^{''})\\
 &=\dots \\
 &=\sum_{i_1=0}^{b} \sum_{i_2=0}^{i_1} \dots\sum_{i_{k-1}=0}^{i_{k-2}} p^{i_1+\dots +i_{k-1}} N_{i_{k-1}}(a_1)\\
 &=\sum_{i_1=0}^{b} \sum_{i_2=0}^{i_1} \dots\sum_{i_{k-1}=0}^{i_{k-2}} p^{i_1+\dots+i_{k-1}}.
\end{align*}

Now let $j_t = b - \sum_{l=1}^t i_t$ for $1\le t\le k-2.$ Then $i_t = j_{t-1} - j_t$ for $1\le t\le k-2,$ letting $j_0=b.$ As $i_t$ varies between 0 and $j_{t-1},$ $j_t$ will also vary between 0 and $j_{t-1}.$ Also $(k-1)b-(k-1)i_1-(k-2)i_2-\cdots -2i_{k-2}-i_{k-1}= j_1+\dots+j_{k-1}.$ Hence,
\begin{align*}
\sum_{j_1=0}^{b} \sum_{j_2=0}^{j_1} \dots\sum_{j_{k-1}=0}^{j_{k-2}} p^{j_1+\dots+j_{k-1}}
&= \sum_{i_1=0}^{b} \sum_{i_2=0}^{j_1} \dots\sum_{i_{k-1}=0}^{j_{k-2}} p^{j_1+\dots+j_{k-1}} \\
&=\sum_{i_1=0}^{b}\sum_{i_2=0}^{b-i_1}\cdots \sum_{i_{k-1}=0}^{b-i_1-\cdots-i_{k-2}} p^{(k-1)b-(k-1)i_1-(k-2)i_2-\cdots -2i_{k-2}-i_{k-1}},
\end{align*} which is the same as the left hand side of equation (\ref{anyrankclaim}). Therefore,
\begin{equation}
\sum_{i_1=0}^{b} \sum_{i_2=0}^{i_1} \dots\sum_{i_{k-1}=0}^{i_{k-2}} p^{i_1+\dots+i_{k-1}} =  \prod_{i=2}^k\frac{p^{b+i-1}-1}{p^{i-1}-1}.
\end{equation}

\end{remark}
\section{The method of convolution} \label{sec:convolution}
Here, we will use convolution of generating series and a recurrence relation of T. Stehling\cite{TS} to prove the rational function formula in case 1 of rank 2. This is only to demo the method and the rest of the cases can be handled similarly. To ease the use of the recurrence relation, let's adopt Stehling's notation for type $\alpha = (\alpha_1,\dots,\alpha_d)$ where $ \alpha_1\ge \cdots \alpha_d\ge 0$ which reverses the order in the notation we used in previous sections.

Let $0\le b \leq \alpha_2$ (i.e. case 1 of rank 2), then
\begin{equation}
\label{c1rk2}
h_b^{(\alpha_1,\alpha_2)}(p) = \frac{ p^{b+1} -1}{p-1}.
\end{equation}
The symmetric case (i.e. $\alpha_1\le b \leq \alpha_1+\alpha_2$) is then (replacing $b$ by $\alpha_1+\alpha_2-b$ in the above rational function):
\begin{equation}
h_b^{(\alpha_1,\alpha_2)}(p) = \frac{ p^{\alpha_1+\alpha_2-b+1} -1}{p-1}.
\end{equation}

The idea behind the method is the following. Let $F(x) = \sum a_nx^n$ and $G(x) = \sum b_nx^n$ be two convergent series. Then the series $H(x) = \sum a_nb_nx^n$  can be obtained from $F$ and $G$ as follows. Define the \textbf{convolution} of $F$ and $G$ denoted $F\star G (x)$ as the following contour integral around a small circle about  $0$ where $x$ is a complex number very small in magnitude.

\begin{align*}
F\star G (x) &= \frac{1}{2\pi i} \oint_{y\in B_{\epsilon}(0)} F(y)G(\frac{x}{y})\frac{dy}{y}\\
&= \frac{1}{2\pi i} \oint_{y\in B_{\epsilon}(0)} \sum a_nb_my^n\frac{x^m}{y^m}\frac{dy}{y} \\
&=\sum a_nb_mx^m \frac{1}{2\pi i} \oint_{y\in B_{\epsilon}(0)} y^{n-m-1}dy\\
&= \sum a_nb_nx^n.
\end{align*}
Thus $H(x) = F\star G (x).$ Now, we will prove equation (\ref{c1rk2}) using this method. Let $F_d(x,y)$ be the generating series
\begin{equation}
F_d(x,y) = \sum_{\substack{\alpha = (\alpha_1,\dots,\alpha_d) \\ \alpha_1\ge \cdots \alpha_d\ge 0  \\  0\le r\le \alpha_1+\dots+\alpha_d}} h_r^{(\alpha_1,\dots,\alpha_d)}(p)x^{\alpha}y^r
\end{equation}
where $x = (x_1,\dots, x_d)$ and $x^\alpha = x_1^{\alpha_1}x_2^{\alpha_2}\cdots x_d^{\alpha_d}.$ Let's denote $h_r^{\alpha}(p)$ by $N_{\alpha}^{(d)}(r)$ or simply $N_{\alpha}(r)$ when the rank is understood. By \cite[Corollary]{TS}, we have the recurrence formula
\begin{equation}\label{Stehrec}
    N_{\alpha}(r) = N_{\tilde{\alpha}}(r-1)+p^rN_{\hat{\alpha}}(r)
\end{equation} where $\hat{\alpha} = (\alpha_2,\dots,\alpha_d)$ and $\tilde{\alpha} = \alpha$ with $-1$ added to the $k^{th}$ position where

$$k=\begin{cases}
  1  & \text{if } \alpha_1>\alpha_2\\
  2  & \text{if } \alpha_1=\alpha_2>\alpha_3 \\
  3  & \text{if } \alpha_1=\alpha_2=\alpha_3>\alpha_4\\
  \dots & \dots
\end{cases}$$

Note $N_{\alpha}^{(1)}(r) =1$ if $0\le r\le \alpha.$ So,
\begin{equation}
    F_1(x,y) = \frac{1}{(1-x)(1-xy)}.
\end{equation} For rank 2, let's break down $F_2(x_1,x_2,y)$ as follows:
\begin{align}
\label{F2}
F_2(x_1,x_2,y) &=  \sum_{\substack{ \alpha_1 \ge \alpha_2\ge 0  \\  0\le r\le \alpha_1+\alpha_2}} N_{\alpha_1,\alpha_2}^{(2)}(r)x_1^{\alpha_1}x_2^{\alpha_2}y^r
= F_2^{(0)}+F_2^{(1)}
\end{align} where
\begin{align}
    F_2^{(0)}(x_1,x_2,y) &=  \sum_{\alpha_1 = \alpha_2\ge 0} N_{\alpha_1,\alpha_2}^{(2)}(r)x_1^{\alpha_1}x_2^{\alpha_2}y^r  \\
    F_2^{(1)}(x_1,x_2,y) &=  \sum_{\alpha_1 > \alpha_2\ge 0} N_{\alpha_1,\alpha_2}^{(2)}(r)x_1^{\alpha_1}x_2^{\alpha_2}y^r
\end{align} Let's first compute $F_2^{(0)}(x_1,x_2,y). $
\begin{align*}
    F_2^{(0)}(x_1,x_2,y) &= \sum_{\alpha=0}^\infty \sum_r N_{\alpha,\alpha}^{(2)}(r)(x_1x_2)^{\alpha}y^r  \\
    &= \sum_{\alpha} \sum_r N_{\alpha,\alpha-1}^{(2)}(r-1)(x_1x_2)^{\alpha}y^r + \sum_{\alpha} p^rN_{\alpha}^{(1)}(r)(x_1x_2)^{\alpha}y^r\\
    &= \sum_{\alpha} \sum_r N_{\alpha-1,\alpha-1}^{(2)}(r-2)(x_1x_2)^{\alpha}y^r + \sum_{\alpha} \sum_r p^{r-1}N_{\alpha-1}^{(1)}(r-1)(x_1x_2)^{\alpha}y^r \\
    &+ \sum_{\alpha} p^rN_{\alpha}^{(1)}(r)(x_1x_2)^{\alpha}y^r\\
    &= \sum_{\alpha} \sum_r N_{\alpha,\alpha}^{(2)}(r)(x_1x_2)^{\alpha+1}y^{r+2} + \sum_{\alpha} \sum_r N_{\alpha}^{(1)}(r)(x_1x_2)^{\alpha+1}p^{r}y^{r+1}\\
    &+ \sum_{\alpha} N_{\alpha}^{(1)}(r)(x_1x_2)^{\alpha}p^ry^r\\
    &= x_1x_2y^2F_2^{(0)}(x_1,x_2,y) + x_1x_2yF_1(x_1x_2,py) + F_1(x_1x_2,py).
\end{align*} Solving for $F_2^{(0)},$ results in \begin{equation}
    F_2^{(0)}(x_1,x_2,y) = \frac{1+x_2x_2y}{(1-x_1x_2)(1-px_1x_2y)(1-x_1x_2y^2)}.
\end{equation} Similarly, we have
\begin{align*}
    F_2^{(1)}(x_1,x_2,y) &= \sum_{\alpha_1>\alpha_2\ge 0} \sum_r N_{\alpha}^{(2)}(r)x_1^{\alpha_1}x_2^{\alpha_2}y^r  \\
    &=\sum_{\alpha_1>\alpha_2\ge 0} \left [ N_{\tilde{\alpha}}^{(2)}(r-1)+p^rN_{\hat{\alpha}}^{(1)}(r) \right ]x_1^{\alpha_1}x_2^{\alpha_2}y^r\\
    &=\sum_{\alpha_1>\alpha_2\ge 0} \left [ N_{\alpha_1-1,\alpha_2}^{(2)}(r-1)+p^rN_{\alpha_2}^{(1)}(r)  \right ]x_1^{\alpha_1}x_2^{\alpha_2}y^r\\
    &=\sum_{\alpha_1 \ge \alpha_2\ge 0}  N_{\alpha_1,\alpha_2}^{(2)}(r-1)x_1^{\alpha_1+1}x_2^{\alpha_2}y^{r+1} + \sum_{ \substack{ \alpha_1 > \alpha_2\ge 0\\ 0\le r\le \alpha_2}}  x_1^{\alpha_1}x_2^{\alpha_2}(py)^r\\
    &= x_1y\left [ F_2^{(1)}(x_1,x_2,y) + F_2^{(0)}(x_1,x_2,y)  \right ] + \sum_{ \alpha_1 > \alpha_2\ge 0}  x_1^{\alpha_1}x_2^{\alpha_2}\frac{1-(py)^{\alpha_2+1}}{1-py}.
\end{align*} Solving for $F_2^{(1)}(x_1,x_2,y),$ we have \begin{equation}
    F_2^{(1)}(x_1,x_2,y) = \frac{1+y-x_1y+x_1^2x_2y^2}{(1-x_1)(1-x_1y)(1-x_1x_2)(1-x_1x_2y)(1-qx_1x_2y)}.
\end{equation} Finally, combining the two sums, equation (\ref{F2}) becomes

\begin{equation*}
F_2(x,y) = F_2(x_1,x_2,y) = \sum_{\substack{\alpha = (\alpha_1,\alpha_2) \\ \alpha_1\ge \alpha_2\ge 0  \\  0\le r\le \alpha_1+\alpha_2}} h_r^{(\alpha_1,\alpha_2)}(p)x_1^{\alpha_1}x_2^{\alpha_2}y^r\end{equation*}
\begin{equation}
= \frac{x_1^2x_2y^2 + x_1^2x_2y - x_1x_2y - 1}{(1-x_1)(1-x_1y)(1-x_1x_2)(1-x_1x_2y^2)(px_1x_2y-1)}.
\end{equation}

We will use the idea described before to pick out the rational function expression for case 1 from the general rational function expression of $F_2$ above. The final step will be extracting a general formula for the coefficients of the general term of the power series expansion of the resulting rational function using products of power series and change of variables.

Now consider small circles in the complex plane around 0 where the complex numbers $x_1,x_2,y$ are very small in magnitude.  We evaluate the following repeated contour integrals of the convolution of $F_2$ and $G$ by applying Cauchy's residue theorem to get the generating series for case 1.

Let
\begin{equation}
G(x_1,x_2,y) = \sum_{ \alpha_1\ge \alpha_2\ge r\ge  0 } x_1^{\alpha_1}x_2^{\alpha_2}y^r = \frac{1}{(1-x_1)(1-x_1x_2)(1-x_1x_2y)}.
\end{equation}
Then
\begin{align*}
F_2\star G (x_1,x_2,y) &=\frac{1}{(2\pi i)^3} \oint_{u_3\in B_{\epsilon_3}(0)} \oint_{u_2\in B_{\epsilon_2}(0)} \oint_{u_1\in B_{\epsilon_1}(0)} F_2(u_1,u_2,v)G(\frac{x_1}{u_1},\frac{x_2}{u_2},\frac{y}{v})\frac{du_1}{u_1}\frac{du_2}{u_2}\frac{dv}{v}  \\
&= \sum_{ \alpha_1\ge \alpha_2\ge r\ge  0 } h_r^{(\alpha_1,\alpha_2)}(p)x_1^{\alpha_1}x_2^{\alpha_2}y^r.
\end{align*}

Using the rational functions for $F_2$ and $G$ and evaluating the integrals, we get that
\begin{equation}
\sum_{ \alpha_1\ge \alpha_2\ge r\ge  0 } h_r^{(\alpha_1,\alpha_2)}(p)x_1^{\alpha_1}x_2^{\alpha_2}y^r = \frac{1}{(1-x_1)(1-x_1x_2)(1-x_1x_2y)(1-px_1x_2y)}.
\end{equation}
It remains to extract the coefficient of $x_1^{\alpha_1}x_2^{\alpha_2}y^r$ from the last rational function in order to determine $h_r^{(\alpha_1,\alpha_2)}(p)$ for the case $0\le r \le a_2.$ We do this as follows:
\begin{align*}
&\frac{1}{(1-x_1)(1-x_1x_2)(1-x_1x_2y)(1-px_1x_2y)}\\
&= \big( \sum_{k_1\ge 0} x_1^{k_1}\big)\big( \sum_{k_2\ge 0} (x_1x_2)^{k_2}\big)\big( \sum_{k_3\ge 0} (x_1x_2y)^{k_3}\big)\big( \sum_{k_4\ge 0} (px_1x_2y)^{k_3}\big)\\
&= \sum_{k_1,k_2,k_3,k_4\ge 0} p^{k_4}x_1^{k_1+k_2+k_3+k_4}x_2^{k_2+k_3+k_4}y^{k_3+k_4}.
\end{align*}

Let $r = k_3+k_4, a_2=k_2+k_3+k_4=k_2+r, a_1=k_1+k_2+k_3+k_4=k_1+a_2.$ The last sum above becomes
\begin{align*}
\frac{1}{(1-x_1)(1-x_1x_2)(1-x_1x_2y)(1-px_1x_2y)} &=\sum_{a_1\ge a_2\ge r\ge 0} x_1^{a_1}x_2^{a_2}y^{r}\sum_{k_4=0}^r p^{k_4} \\
&= \sum_{a_1\ge a_2\ge r\ge 0} \frac{p^{r+1}-1}{p-1}x_1^{a_1}x_2^{a_2}y^{r}.
\end{align*}

Finally, by comparing coefficients, we get the required result:
$$h_r^{(\alpha_1,\alpha_2)}(p) = \frac{p^{r+1}-1}{p-1}.$$
\begin{remark}
In comparison, we can quickly prove the formulas for all cases of rank 2 using the recurrence used in the proof of Theorem \ref{rank3} as follows. First let's return to the notation of the previous sections, i.e. $a = (a_d, \dots,a_1)$ with $a_1\leq \dots \le a_d.$

Let $\lambda=(a_2, a_1)$ and $\lambda^{'}=( a_1).$\\
\text{Case 1:}  $0\le b \leq a_1$
$N_{b}(\lambda)=\sum_{i=0}^{b} p^{i}N_{i}(\lambda^{'})=\sum_{i=0}^{b} p^{i}=\frac{p^{b+1}-1}{p-1}.$\\
Note $N_{i}(\lambda^{'})=1$ or 0 according as $ i \leq a_1$ or $ i> a_1 ,$ respectively.\\
\text{Case 2:}  $a_1\le b \leq a_2$
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{b} p^{i}N_{i}(\lambda^{'})=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})=\sum_{i=0}^{a_1} p^{i}=\frac{p^{a_1+1}-1}{p-1}.$\\
\text{Case 3:}  $a_2\le b \leq a_{1}+a_2$ (This case can also be deduced by replacing $b$ in case 1 by $a_1+a_2-b$)
$N_{b}(\lambda)=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{1}+1}^{a_{2}} p^{i}N_{i}(\lambda^{'})+\sum_{i=a_{2}+1}^{b} p^{i}N_{i}(\lambda^{'})-\sum_{i=a_{1}+a_{2}+1-b}^{a_1}p^{i}N_{i}(\lambda^{'})$\\
$=\sum_{i=0}^{a_1} p^{i}N_{i}(\lambda^{'})-\sum_{i=a_{1}+a_{2}+1-b}^{a_1}p^{i}=\frac{p^{a_{1}+a_{2}+1-b}-1}{p-1}.$\\

In summary, we have
$$h_b^{(a_2,a_1)}(p)=\begin{cases}
  \frac{p^{b+1}-1}{p-1}  & \text{if } 0\le b \leq a_2\\
  \frac{p^{a_2+1}-1}{p-1}  & \text{if } a_2\le b \leq a_1\\
  \frac{p^{a_1+a_2-b+1}-1}{p-1}  & \text{if } a_1\le b \leq a_1+a_2.
\end{cases}$$
\end{remark}

\textit{Acknowledgement.} The first named author would like to thank the City University of New York HPCC facilities for letting us verify every formula using the SageMath and Mathematica software on their HPC systems based at the College of Staten Island, New York City.

\bibliographystyle{hplain}

Group Theory (math.GR); Combinatorics (math.CO)