The Pseudoforest analogue for the Strong Nine Dragon Tree Conjecture is True

Abstract: We prove that for any positive integers $k$ and $d$, if a graph $G$ has maximum average degree at most $2k + \frac{2d}{d+k+1}$, then $G$ decomposes into $k+1$ pseudoforests $C_{1},\ldots,C_{k+1}$ such that there is an $i$ such that for every connected component $C$ of $C_{i}$, we have that $e(C) \leq d$.

\section{Introduction}

For any graph $G$, a \textit{decomposition} of $G$ is a set of edge disjoint subgraphs of $G$ such that the union of their edges sets is the edge set of the graph.  Graph decompositions are a particularly heavily studied area of graph theory, and one of the most beautiful results about graph decompositions is Nash-Williams characterization of when a graph decomposes into $k$ forests.

\begin{thm}[\cite{Nashwilliams}, Nash-Williams Theorem]
A graph $G$ decomposes into $k$ forests if and only if
\[ \max_{H \subseteq G, v(H) \geq 2} \frac{e(H)}{v(H)-1} \geq k.\]
\end{thm}

Here, $H \subseteq G$ is taken to mean $H$ is a subgraph of $G$. We will refer to $\max_{H \subseteq G, v(H) \geq 2} \frac{e(H)}{v(H)-1}$ as the \textit{fractional arboricity} of $G$.  It is not hard to see that the fractional arboricity of a graph might not be integral. For instance, the fractional arboricity of a cycle on $n$ vertices is $\frac{n}{n-1}$.  By Nash-Williams Theorem, this implies that cycles decompose into two forests, and you cannot decompose a cycle into a single forest.  Additionally, it is easy to see that cycles decompose into a forest and a matching, which is quite a bit stronger than just saying that cycles decompose into two forests. As the fractional arboricity of cycles is much closer to $1$ than to $2$, you might speculate that this is the reason you get this extra structure in the decomposition.  In general,  one might speculate that if the fractional arboricity of a graph is much closer to $k-1$ than to $k$, then not only does $G$ decompose into $k$ forests, but that one of the forests can be assumed to be a matching. Intuitively, this should be believable. If the fractional arboricity is very close to $k-1$ but still greater than $k-1$, then Nash-Williams Theorem says that you need $k$ forests, but just barely. Such a graph would roughly look like the union of $k-1$ forests, and then a few edges left over. It is reasonable to believe that you can force these left over edges form a matching. Even more generally, you could conjecture that there exists an $\varepsilon  \in (0,1)$ such that if the fractional arboricity of $G$ is at most $k + \varepsilon$, then $G$ decomposes into $k+1$ forests, such that one of the forests has maximum degree $d$. The Nine Dragon Tree Conjecture (now proven by Hongbi Jiang and Daqing Yang \cite{ndt}), proposed by Mickeal Montassier, Patrice Ossona de Mendez, Andr Raspaud, and Xuding Zhu in \cite{montassier} does just this.

\begin{thm}[\cite{ndt}, Nine Dragon Tree Theorem]
Let $k$ and $d$ be positive integers. If the fractional arboricity of $G$ is at most $k + \frac{d}{k+d+1}$, then $G$ decomposes into $k+1$ forests such that one of the forests has maximum degree $d$.
\end{thm}

In \cite{montassier}, it was shown that the fractional arboricity bound is sharp for arbitrarily large graphs. Despite this, Montassier et al. proposed a significant strengthening of the Nine Dragon Tree Theorem, aptly named the Strong Nine Dragon Tree Conjecture.

\begin{conj}[\cite{montassier}, Strong Nine Dragon Tree Conjecture]
\label{sndt}
Let $k$ and $d$ be positive integers. If the fractional arboricity of $G$ is at most $k + \frac{d}{k+d+1}$, then $G$ decomposes into $k+1$ forests such that one of the forests has each connected component containing at most $d$ edges.
\end{conj}

This conjecture is wide open. The $d =1$ case was proven by Yang in \cite{Yangmatching}. The $k=1$ and $d=2$ case was proven by Kim, Kostochka, West, Wu, and Zhu \cite{kostochkaetal}.  Recently, the second author resolved the $d \leq k+1$ case \cite{moore2019approximate}. Moreover, it was shown that the conjecture is true if you replace ``$d$ edges" in the conclusion with a function $f(d,k)$ edges.

This paper will focus on pseudoforest decompositions. Recall, a \textit{pseudoforest} is a graph where each connected component contains at most one cycle. All of the above results or conjectures have pseudoforest analogues. The pseudoforest analogue of Nash-Williams Theorem is Hakimi's Theorem.

\begin{thm}[\cite{hakimithm}, Hakimi's Theorem]
A graph $G$ decomposes into $k$ pseudoforests if and only if the maximum average degree of $G$ is at most $2k$.
\end{thm}

Here, the maximum average degree of a graph $G$ is
\[\text{mad}(G) := \max_{H \subseteq G} \frac{2e(H)}{v(H)}.\]

Genghua Fan, Yan Li, Nine Song, and Daqing Yang \cite{pseudoforest} proved the pseudoforest analogue of the Nine Dragon Tree Theorem.

\begin{thm}[\cite{pseudoforest}]
\label{pseudondt}
Let $k$ and $d$ be positive integers. If $\text{mad}(G) \leq 2k + \frac{2d}{d+k+1}$, then $G$ decomposes into $k+1$ pseudoforests, such that one of the pseudoforests has maximum degree $d$.
\end{thm}

Fan et al. also showed the maximum average degree bound is best possible.

\begin{thm}[\cite{pseudoforest}]
\label{sharpness}
For any positive integers $k$ and $d$, there are arbitrarily large graphs $G$ and an edge $e$ where $\text{mad}(G-e) = 2k + \frac{2d}{k+d+1}$, but $G$ does not decompose into $k+1$ pseudoforests where one of the pseudoforests has maximum degree $d$.
\end{thm}

The main result of this paper is that the pseudoforest analogue of the Strong Nine Dragon Tree Conjecture is true.

\begin{thm}
\label{strongpndt}
Let $k$ and $d$ be positive integers. If $\text{mad}(G) \leq 2k + \frac{2d}{d+k+1}$, then $G$ decomposes into $k+1$ pseudoforests $T_{1},\ldots,T_{k},F$, such that each connected component of $F$ contains at most $d$ edges.
\end{thm}

By Theorem \ref{sharpness}, the maximum average degree bound given in Theorem \ref{strongpndt} is sharp. Some of the cases of Theorem \ref{strongpndt} were known. As having maximum degree one is the same as having one of the pseudoforests be a matching, Theorem \ref{pseudondt} implies Theorem \ref{strongpndt} when $d=1$. Interestingly, the proof given in \cite{kostochkaetal} of the Strong Nine Dragon Tree Conjecture when $k = 1$ and $d=2$ implies Theorem \ref{strongpndt} when $k=1$ and $d=2$. Prior to our result, all other cases were open.

As a template for how the proof of Theorem \ref{strongpndt} will proceed, we will give a proof of the non-trivial direction Hakimi's Theorem. The proof we give appears in \cite{pseudoforest} (and is perhaps the first time it appeared in print, we are not aware of any earlier proofs), and the proof of the Nine Dragon Tree Theorem for pseudoforests in some sense follows this proof of Hakimi's Theorem. However, we feel that our proof more faithfully follows this proof of Hakimi's Theorem, and thus leads to the stronger result while having a shorter proof.

Before proceeding, we need some definitions. Given a graph $G$, an \textit{orientation} of $G$ is obtained from $E(G)$ by taking each edge $xy$, and replacing $xy$ with exactly one of the arcs $(x,y)$ or $(y,x)$. To reverse the direction of an arc $(x,y)$ is to replace $(x,y)$ with the arc $(y,x)$. For any vertex $v$, let $d(v)$ denote the degree of $v$, and $d^{+}(v)$ denote the outdegree of $v$. A directed path $P$ from $u$ to $v$ is a path $P$ oriented so that $v$ is the only vertex with no outgoing edge.  The next observation is easy and well known.

\begin{obs}
A graph $G$ is a pseudoforest if and only if $G$ admits an orientation where every vertex has outdegree at most one.
\end{obs}

From this observation, we get an important corollary.

\begin{cor}
\label{easycor}
A graph admits a decomposition into $k$ pseudoforests if and only if it admits an orientation such that every vertex has outdegree at most $k$.
\end{cor}

For a proof of Corollary \ref{easycor} we refer the reader to Corollary $1.2$ and Theorem $1.1$ of \cite{pseudoforest}. We will use Corollary \ref{easycor} repeatedly and implicitly throughout our proofs. With that, we can give a proof of Hakimi's Theorem.

\begin{thm}
Every graph with maximum average degree at most $2k$ decomposes into $k$ pseudoforests.
\end{thm}

\begin{proof}
Suppose towards a contradiction that $G$ has maximum average degree at most $2k$, but $G$ does not decompose into $k$ pseudoforests. Then $G$ does not admit an orientation such that each vertex has outdegree at most $k$.

Consider an orientation $\vec{G}$ of $G$ that minimizes the sum
\[ \rho := \sum_{v \in V(G)} \max\{0, d^{+}(v) -k\}.\]
If this sum is zero, then we have a desired decomposition, a contradiction. Thus there is a vertex $v \in V(G)$ such that $v$ has out degree at least $k+1$. If there is a directed path $P$ from $v$ to $x$ such that $x$ has outdegree at most $k-1$, then we can reverse the directions on all of the arcs on $P$ and obtain a decomposition with smaller $\rho$ value, a contradiction. Consider the subgraph $H$ induced by vertices which are reachable from directed paths of $v$. That is, $x \in V(H)$ if there is a directed path $P$ which starts at $v$ and ends at $x$. Then all vertices in $H$ have outdegree at least $k$, and $v$ has outdegree at least $k+1$. But this implies that the average degree of $H$ is strictly larger than $2k$, a contradiction.
\end{proof}

Now we will give a high level overview of how our proof will proceed. We will take a pseudoforest decomposition $C_{1},\ldots,C_{t},F$ where we will try and bound the size of each connected component in $F$. In the above proof of Hakimi's Theorem, the bad situation was a vertex which had too large outdegree. Now, the bad situation is that there is a component which is too large. In the proof of Hakimi's Theorem, we searched for special paths to augment on from a vertex which had too large outdegree, and in our proof, we will search for paths to augment on from a component which is too large. In the proof of Hakimi's theorem, we
identified a situation where we could augment our decomposition and obtain a better decomposition, namely, directed paths from a vertex with too large outdegree to a vertex with small outdegree. In our proof, we will identify similar situations, namely, finding two components which are small enough to augment our decomposition, or finding a large component which has at least two small components nearby to perform augmentations. Then we will show that when these configurations are removed, either we have a decomposition satisfying Theorem \ref{strongpndt} or our graph actually had too large maximum average degree to begin with.

The paper is structured as follows. In Section \ref{mincounterexample} we  give the necessary definitions on how we will pick our counterexample, and prove basic properties about the counterexample. In Section \ref{flipping}, we describe how we will augment our decomposition in certain situations. In Section \ref{keyarguments}, we show how to use this augmentation strategy to either find an optimal decomposition or show that our graph has too large maximum average degree.

Finally, to fix some notation we will let $v(G) = |V(G)|$ and $e(G) = |E(G)|$. All other undefined graph theory terminology can be found in \cite{BondyAndMurty}, or any other standard graph theory textbook.

\section{Picking the counterexample}
\label{mincounterexample}
In this section we describe how we will pick our counterexample. Fix $k$ and $d$ as positive integers, and suppose that $G$ is a vertex minimal counterexample to Theorem \ref{strongpndt} for the fixed values of $k$ and $d$.

Our first step will be to obtain desirable orientations of $G$. In particular, the orientations we will demand will imply that $G$ decomposes into $k$ pseudoforests each with $v(G)$ edges, and one left over pseudoforest. For this, we use a lemma proved in \cite{pseudoforest} (Lemma $2.1$). Technically, we need a stronger lemma, however the same proof as Lemma $2.1$ will suffice. We give a proof for completeness sake only, there is no new idea needed.

\begin{lemma}[\cite{pseudoforest}]
\label{orientationlemma}
If $G$ is a vertex minimal counterexample to Theorem \ref{strongpndt}, then there exists an orientation of $G$ such that for all $v \in V(G)$, we have $k \leq d^{+}(v) \leq k+1$.
\end{lemma}

\begin{proof}
Suppose no such orientation exists. As $G$ has maximum average degree at most $2k + 2$, by Hakimi's Theorem, $G$ admits an orientation so that every vertex has outdegree at most $k+1$. Orient $G$ so that every vertex has outdegree at most $k+1$, and that the sum
\[\rho := \sum_{v \in V(G)} \max\{0, k - d^{+}(v)\}\]
is minimized. Observe that if $\rho$ is zero, then we have a desirable orientation.

First we claim there is a vertex $v$ with outdegree $k+1$. If not, then all vertices have outdegree at most $k$, and by Hakimi's Theorem $G$ decomposes into $k$ pseudoforests, contradicting that $G$ is a counterexample to Theorem \ref{strongpndt}.

Now we claim there is no directed path $P$ from a vertex $v$ with outdegree $k+1$ to a vertex $u$ with outdegree at most $k-1$. Suppose towards a contradiction that $P$ is such a path. Then reversing the orientation on all of the arcs in $P$ gives a new orientation, where $v$ has outdegree $k$, all internal vertices have the same outdegree, and the outdegree of $u$ increases by one. But this contradicts that we picked our orientation to minimize $\rho$.

Let $S$ be the set of vertices in $G$ with out degree at most $k-1$, and let $S'$ be the set of vertices which have a directed path to a vertex in $S$. Observe that every vertex in $S'$ has outdegree at most $k$. Let $\bar{S'} = V(G) - S'$. Then every edge with one endpoint lying in $S'$ and one endpoint in $\bar{S'}$ is directed from $S'$ to $\bar{S'}$. Observe that $|\bar{S'}| < v(G)$.

As $G$ is a vertex minimal counterexample we can decompose $G[\bar{S'}]$ into $k+1$ pseudoforests such that one of the pseudoforests has each connected component containing at most $d$ edges. Additionally, as every vertex in $S'$ has outdegree at most $k$, by Hakimi's Theorem we can decompose $G[S']$ into $k$ pseudoforests $C_{1},\ldots,C_{k}$. Thus we only need to deal with the edges between $\bar{S}'$ and $S'$. Observe that if $v$ has $t$ arcs $(v,u_{1}),\ldots,(v,u_{t})$ where $u_{i} \in \bar{S'}$ for all $i \in \{1,\ldots,t\}$, then $v$ has outdegree at most $k-t$ in $G[S']$. Thus $v$ has outdegree zero in at least $t$ of the pseudoforests $C_{1},\ldots,C_{k}$. Therefore we can add the arcs $(v,u_{1}),\ldots,(v,u_{t})$ to $t$ of the pseudoforests so that the result is a pseudoforest. As all arcs between $S'$ and $\bar{S}'$ are oriented from $S'$ to $\bar{S}'$, we now have a decomposition of $G$ which satisfies Theorem \ref{strongpndt}. But this contradicts that $G$ is a counterexample to Theorem \ref{strongpndt}.
\end{proof}

 Let $\mathcal{F}$ be the set of orientations of $E(G)$ which satisfy Lemma \ref{orientationlemma}. A useful way of keeping track of our pseudoforest decomposition will be to colour the edges blue and red, where the edges coloured red will induce a pseudoforest. This will be the pseudoforest that we will want to bound the size of each connected component.

\begin{definition}
Suppose we orient $G$ from an orientation in $\mathcal{F}$. Then a \textit{red-blue colouring of $G$} is a (non-proper)-colouring of the edges where for any vertex $v \in V(G)$, we colour $k$ outgoing arcs of $v$ blue, and if after this there is an uncoloured outgoing arc, colour this arc red.
\end{definition}

Note that given an orientation in $\mathcal{F}$, one can generate many different red-blue colourings. As a graph decomposes into $k$ pseudoforests if and only if it admits an orientation where each vertex has outdegree at most $k$, we obtain the following observation.

\begin{obs}
\label{redbluecolouringdecomp}
Given a red-blue colouring of $G$, we can decompose our graph $G$ into $k+1$ pseudoforests such that $k$ of the pseudoforests have all of their edges coloured blue, and the other pseudoforest has all of its edges coloured red.
\end{obs}

Observe that one red-blue colouring can give rise to many different pseudoforest decompositions. Given a pseudoforest decomposition obtained from Observation \ref{redbluecolouringdecomp} we will say a pseudoforest which has all arcs coloured blue is a \textit{blue pseudoforest}, and the pseudoforest with all arcs coloured red is the \textit{red pseudoforest}.

\begin{definition}
Let $f$ be a red-blue colouring of $G$, and let $C_{1},\ldots,C_{k},F$ be a pseudoforest decomposition obtained from $f$ by Observation \ref{redbluecolouringdecomp}. Then we say that \textit{$C_{1},\ldots,C_{k},F$ is a pseudoforest decomposition generated from $f$}. We will always use the convention that $F$ is the red pseudoforest, and $C_{i}$ is a blue pseudoforest.
\end{definition}

 As $G$ is a counterexample, in every pseudoforest decomposition generated from a red-blue colouring, there is a connected component of the red pseudoforest which has more than $d$ edges. We define a \textit{residue function} which simply measures how close a decomposition is to satisfying Theorem \ref{strongpndt}.

\begin{definition}
Let $f$ be a red-blue colouring and $C_{1},\ldots,C_{k},F$ be a pseudoforest decomposition generated by $f$. Let $\mathcal{T}$ be the set of components $F$. Then the \textit{residue function}, denoted $\rho$, is
\[\rho(F) = \sum_{K \in \mathcal{T}} \max\{e(K) -d,0\}.\]
\end{definition}

Using a red-blue colouring, and the resulting pseudoforest decomposition, we define an induced subgraph of $G$ which we will focus our attention on. Intuitively, this subgraph should be thought of as an ``exploration" subgraph similar to how in the proof of Hakimi's theorem we ``explored" from a vertex which had too large outdegree. Here we will ``explore" from a component which is too large.

\begin{definition}
Suppose that $f$ is a red-blue colouring of $G$, and suppose $D= (C_{1},\ldots,C_{k},F)$ is a pseudoforest decomposition generated from $f$. Let $R$ be an component of $F$ such that $e(R) >d$. We define the \textit{exploration subgraph} $H_{f,D,R}$ in the following manner. Let $S \subseteq V(G)$ where $v \in S$ if and only if there exists a path $P = v_{1},\ldots,v_{m}$ such that $v_{m}= v$, $v_{1} \in V(R)$, and either $v_{i}v_{i+1}$ is an arc $(v_{i},v_{i+1})$ coloured blue, or $v_{i}v_{i+1}$ is an arbitrarily directed arc coloured red. Then we let $H_{f,D,R}$ be the graph induced by $S$.
\end{definition}

Given a particular exploration subgraph $H_{f,D,R}$, we say $R$ is the \textit{root component}. We say the \textit{red components} of $H_{f,D,R}$ are the components of $F$ contained in $H_{f,D,R}$.

It might not be clear why we made this particular definition for $H_{f,D,R}$, however the next observation shows that for any exploration subgraph $H_{f,D,R}$, the red edge density must be low. Before stating the observation, we fix some notation. Given a subgraph $K$ of $G$, we will let $E_{b}(K)$ and $E_{r}(K)$ denote the set of edges of $K$ coloured blue and red, respectively. We let $e_{b}(K) = |E_{b}(K)|$ and $e_{r}(K) = |E_{r}(K)|$.

\begin{obs}
\label{finishingobservation}
For any red-blue colouring $f$, any pseudoforest decomposition $D$ generated from $f$, and any choice of root component $R$, the exploration subgraph $H_{f,D,R}$ satisfies
\[\frac{e_{r}(H_{f,D,R})}{v(H_{f,D,R})} \leq \frac{d}{d+k+1}.\]
\end{obs}

\begin{proof}

Suppose towards a contradiction that
\[\frac{e_{r}(H_{f,D,R})}{v(H_{f,D,R})} > \frac{d}{d+k+1}.\]

As $H_{f,D,R}$ is an induced subgraph defined by directed paths, and every vertex $v \in V(G)$ has $k$ outgoing blue edges, each vertex in $H_{f,D,R}$ has $k$ outgoing blue edges. Thus,
\[\frac{e_{b}(H_{f,D,R})}{v(H_{f,D,R})} = k.\]

Then we have

\[\frac{\text{mad}(G)}{2} \geq \frac{e(H_{f,D,R})}{v(H_{f,D,R})} = \frac{e_{r}(H_{f,D,R})}{v(H_{f,D,R})} + \frac{e_{b}(H_{f,D,R})}{v(H_{f,D,R})} > k + \frac{d}{d+k+1}.\]

But this contradicts that $G$ has $\text{mad}(G) \leq 2k +\frac{2d}{d+k+1}$.

\end{proof}

For the entire proof, we will be attempting to show that we can augment a given decomposition in such a way that either we obtain a decomposition satisfying Theorem \ref{strongpndt} or we can find a exploration subgraph $H_{f,D,R}$ which contradicts Observation \ref{finishingobservation}.

As Observation \ref{finishingobservation} allows us to focus only on red edges, it is natural to focus on red components which have small average degree. With this in mind, we define the notion of a \textit{small red component}.

\begin{definition}
Let $C_{1},\ldots,C_{t},F$ be a pseudoforest decomposition generated by a red-blue colouring. Let $K$ be a subgraph of $F$. Then $K$ is a \textit{small red subgraph} if
\[\frac{e_{r}(K)}{v(K)} < \frac{d}{d+k+1}.\]
\end{definition}

In particular, we will be interested in the case when $K$ is connected and a small red component. When $K$ is connected the red subgraph is actually isomorphic to a tree, and we can rewrite the density bound in the definition in a more convenient manner.

\begin{obs}
\label{AcyclicObservation}
Let $K$ be a connected small red subgraph. Then $K$ is acyclic, and further
\[e_{r}(K) < \frac{d}{k+1}.\]
\end{obs}

\begin{proof}
First, suppose that $K$ is not acyclic. Then $K$ contains exactly one cycle. As $K$ is connected, it follows that
\[\frac{e_{r}(K)}{v(K)} = 1.\]
But $\frac{d}{k+d+1} < 1$ as $d$ and $k$ are positive integers, and hence $K$ would not be a small red subgraph. Thus we can assume that $K$ is acyclic, so $e(K) = v(K)-1$. Thus

\[\frac{e(K)}{v(K)} = \frac{e(K)}{e(K)+1} < \frac{d}{d+k+1}.\]
Therefore
\[e(K)(d+k+1) < d(e(K) + 1).\]
Simplifying, we see that this is equivalent to
\[e(K) < \frac{d}{k+1}.\]
\end{proof}

We will want to augment our decomposition, and we will want a notion of progress that our decomposition is improving. Of course, if we reduce the residue function that clearly improves the decomposition. However, this might not always be possible, so we will introduce a notion of a ``legal order" of the red components. This order keeps track of the number of edges in components which are ``close" to the root component, with the idea being that if we can continually perform augmentations to make components ``closer" to the root component have fewer edges without creating any large components, then we eventually reduce the number of edges in the root component, which improves the residue function.   We formalize this in the following manner.

\begin{definition}
  We call an ordering $(R_{1},\ldots,R_{t})$ of the red components of $H_{f,D,R}$ \textit{legal} if all red components are in the ordering, $R_{1}$ is the root component, and for all $j \in \{2,\ldots,t\}$ there exists an integer $i$ with $1 \leq i < j$ such that there is a blue arc $(u,v)$ such that $u \in V(R_{i})$ and $v \in V(R_{j})$.
\end{definition}

Let $(R_{1},\ldots,R_{t})$ be a legal ordering. We will say that $R_{i}$ is a \textit{parent} of $R_{j}$ if $i < j$ and there is a blue arc $(v_{i},v_{j})$ where $v_{i} \in R_{i}$  and $v_{j} \in R_{j}$. In this definition a red component may have many parents. To remedy this, if a red component has multiple parents, we arbitrarily pick one such red component and designate it as the only parent. If $R_{i}$ is the parent of $R_{j}$, then we say that $R_{j}$ is a \textit{child} of $R_{i}$. We say a red component $R_{i}$ is an \textit{ancestor} of $R_{j}$ if we can find a sequence of red components $R_{i_{1}},\ldots,R_{i_{m}}$ such that $R_{i_{1}} = R_{i}$, $R_{j_{m}} = R_{j}$, and $R_{i_{q}}$ is the parent of $R_{i_{q+1}}$ for all $q \in \{1,\ldots,m-1\}$. An important definition is that of vertices determining a legal order.
\begin{definition}
Given a legal order $(R_{1},\ldots,R_{t})$, we say a vertex $v$ \textit{determines the legal order for $R_{j}$} if there is a blue arc $(u,v)$ such that $u \in R_{i}$ and $v \in R_{j}$ such that $i < j$.
\end{definition}
  Observe there may be many vertices which determine the legal order for a given red component. More importantly, for every component which is not the root, there exists a vertex which determines the legal order. We also want to compare two different legal orders.

\begin{definition}
Let $(R_{1},\ldots,R_{t})$ and $(R'_{1},\ldots,R'_{t'})$ be two legal orders. We will say $(R_{1},\ldots,R_{t})$ is \textit{smaller} than $(R'_{1},\ldots,R'_{t})$ if  the sequence $(e(R_{1}),\ldots,e(R_{t}))$ is smaller lexicographically than $(e(R'_{1}),\ldots,e(R'_{t'}))$.

\end{definition}

With this, we will pick our minimal counterexample in the following manner.
First, we pick our counterexample to be minimized with respect to the number of vertices. Then, we pick an orientation in $\mathcal{F}$. Then we take a red-blue colouring $f$ of this orientation. With this red-blue colouring, we pick a a pseudoforest decomposition $D = (C_{1},\ldots,C_{k},F)$ generated by $f$. Then we pick a root component $R$ of $F$. From this root, we pick a legal order $(R_{1},\ldots,R_{t})$ of $H_{f,D,R}$ subject to the following conditions in the following order:

\begin{enumerate}
\item{The number of cycles in $F$ is minimized,}

\item{subject to the above condition, we minimize the residue function $\rho$,}

\item{and subject to both of the above conditions, $(R_{1},\ldots,R_{t})$ is the smallest legal order.}

\end{enumerate}

From here on out, we will assume we are working with a counterexample picked in the manner described. The point of minimizing the number of cycles in $F$ is slightly unintuitive compared to minimizing the residue function and minimizing the legal order. However, we minimize the number of cycles because when we augment we will need to ensure our decomposition is in fact a pseudoforest decomposition, and our augmentations will never create more cycles in $F$. Hence by minimizing the number of cycles in $F$ first, we can easily take care of the cases where cycles occur, which allows us to focus on the more important cases where the components are acyclic.

 It was pointed out to the second author via personal communication by Daqing Yang that one can remove the condition to minimize the number of cycles in $F$ by instead  modifying the definition of legal order to take into account the maximum average degree instead of just the number of edges. This is a nice approach and an argument could be made that this approach is more intuitive. However the proof ends up being fundamentally the same, and the authors prefer this approach.

\begin{figure}
\begin{center}
\begin{tikzpicture}
\node[blackvertex] at (0,0) (v1) {};
\node[blackvertex] at (2,0) (v2) {};
\node[blackvertex] at (4,0) (v3) {};

\node[blackvertex] at (4,-1) (v4) {};
\node[blackvertex] at (2,-1) (v5) {};

\node[blackvertex] at (0,-1) (v6) {};

\node[blackvertex] at (0,-2) (v7) {};

\node[blackvertex] at (2,-2) (v8) {};

\draw[ultra thick,black, ->,>=stealth] (v1)--(v2);
\draw[ultra thick,black, ->, >=stealth] (v2)--(v3);
\draw[ultra thick,black, ->, >=stealth] (v1)--(v6);
\draw[ultra thick,black, ->, >= stealth] (v6)--(v7);
\draw[ultra thick,black, ->, >= stealth] (v4)--(v5);
\draw[ultra thick,black, ->, >= stealth] (v3)--(v4);
\draw[ultra thick,black, ->, >= stealth] (v2)--(v5);
\draw[ultra thick,black, ->, >= stealth] (v5)--(v8);
\draw[ultra thick,black,->, >= stealth] (v8)--(v7);
\draw[ultra thick,black,->, >= stealth] (v7)--(v2);
\draw[ultra thick,black,->, >= stealth] (v4)--(v8);

\begin{scope}[xshift = 6cm]
\node[blackvertex] at (0,0) (v1) {};
\node[blackvertex] at (2,0) (v2) {};
\node[blackvertex] at (4,0) (v3) {};

\node[blackvertex] at (4,-1) (v4) {};
\node[blackvertex] at (2,-1) (v5) {};

\node[blackvertex] at (0,-1) (v6) {};

\node[blackvertex] at (0,-2) (v7) {};

\node[blackvertex] at (2,-2) (v8) {};

\draw[ultra thick,red, ->,>=stealth] (v1)--(v2);
\draw[ultra thick,red, ->, >=stealth] (v2)--(v3);
\draw[ultra thick,blue, ->, >=stealth] (v1)--(v6);
\draw[ultra thick,blue, ->, >= stealth] (v6)--(v7);
\draw[ultra thick,red, ->, >= stealth] (v4)--(v5);
\draw[ultra thick,blue, ->, >= stealth] (v3)--(v4);
\draw[ultra thick,blue, ->, >= stealth] (v2)--(v5);
\draw[ultra thick,blue, ->, >= stealth] (v5)--(v8);
\draw[ultra thick,blue,->, >= stealth] (v8)--(v7);
\draw[ultra thick,blue,->, >= stealth] (v7)--(v2);
\draw[ultra thick,blue,->, >= stealth] (v4)--(v8);

\node[smallwhite] at (-.5,0) (v9) {$R_{1}$};
\node[smallwhite] at (-.5,-1) (v10) {$R_{2}$};
\node[smallwhite] at (-.5,-2) (v11) {$R_{3}$};
\node[smallwhite] at (2.7,-2) (v12) {$R_{5}$};
\node[smallwhite] at (1.55,-1) (v13) {$R_{4}$};
\end{scope}
\end{tikzpicture}
\begin{caption}
\\ In this example we assume $k=1$ and $d=1$. On the left, the orientation is in $\mathcal{F}$. On the right we have one possible red-blue colouring generated by this orientation. Here, the entire graph would be the exploration subgraph, and assuming $R_{1}$ is the root, $(R_{1},R_{2},R_{3},R_{4},R_{5})$ is the smallest legal order. Lastly, the isolated vertices are the small components (and are in fact the only possible small components when $k=1$ and $d=1$)
\end{caption}
\end{center}
\end{figure}

\section{Augmenting the decomposition}
\label{flipping}

In this section we describe a very simple operation which will mostly be how we augment our decomposition. Let $f$ be the red-blue colouring of our counterexample, and let $C_{1},\ldots,C_{k}$ be the blue pseudoforests, and $F$ the red pseudoforest. Let $(R_{1},\ldots,R_{t})$ be the legal ordering picked for our counterexample. As some notation, given a vertex $x \in V(H_{f,D,R})$, we let $R^{x}$ denote the red component of $H_{f,D,R}$ containing $x$.
\begin{definition}
Let $(x,y)$ be a blue arc. Further suppose that  $R^{y}$ is acyclic, and suppose that $e = xv$ is an arbitrarily oriented red arc incident to $x$. To \textit{exchange $e$ and $(x,y)$} is to perform the following procedure. First, take the maximal directed red path in $R^{y}$ starting at $y$, say $Q = y,v_{1},\ldots,v_{l}$ where $(v_{i},v_{i+1})$ is a red arc for $i \in \{1,\ldots,l\}$ and $(y,v_{1})$ is a red arc, and reverse the direction of all arcs of $Q$. Second, change the colour of $(x,y)$ to red and reorient $(x,y)$ to $(y,x)$. Finally, change the colour of $e$ to blue, and if $e$ is oriented $(v,x)$, reorient to $(x,v)$.
\end{definition}

We note that exchanging on an edge $e$ and $(x,y)$ is well-defined. This is because $R^{y}$ is acyclic (and hence a tree), and thus there is a unique maximal directed path in $R^{y}$ which starts at $y$.

\begin{obs}
\label{flipobservation}
Suppose we exchange the edge $e =xv$ and $(x,y)$. Then the resulting orientation is in $\mathcal{F}$, and the resulting colouring is a red-blue colouring of this orientation.
\end{obs}

\begin{proof}

Let us first check the outdegrees of vertices after the exchange. Let $Q = y,v_{1},\ldots,v_{l}$ be the maximal directed red path in $R^{y}$ before exchanging the edges. First suppose that $Q$ is not just $y$. Then all of the internal vertices on this path have the same outdegree after reversing as before. On the other hand, the outdegree of $y$ decreases by $1$, and the outdegree of $v_{l}$ increases by one. As $Q$ is a maximal red path, this implies that the outdegree of $v_{l}$ before reversing the arcs on $Q$ was $k$, and hence after reversing the arcs has outdegree $k+1$. The outdegree of $y$ drops by one after reversing the arcs on $Q$, but we reverse the arc $(x,y)$ to $(y,x)$, and hence the outdegree of $y$ remains the same as before.

If $Q$ is just $y$, then the outdegree of $y$ before exchanging was $k$, and then as we reorient $(x,y)$ to $(y,x)$, the outdegree of $k$ is now $k+1$.

Focusing on $x$ now, if $e$ is oriented from $x$ to $v$, then the outdegree of $x$ is initially $k+1$, and after the exchange ends up being $k$, and otherwise, we reorient $(v,x)$ to $(x,v)$ and so the outdegree of $x$ remains the same as before the exchange. Lastly, the outdegree of $v$ remains the same if $e$ was oriented $(x,v)$, and otherwise the outdegree of $v$ initially was $k+1$, and after reorienting becomes $k$. Thus the resulting orientation is in $\mathcal{F}$.

Now we will see that this new colouring is a red-blue colouring of the orientation. Note after exchanging $e$ and $(x,y)$, $y$ has exactly one outgoing edge coloured red. If $x$ had no outgoing red edge before, it still has no outgoing red edge, and if it did have an outgoing red edge, then the outdegree of $x$ dropped by one, and now $x$ has no outgoing red edge. Finally, if $Q$ was not just $y$, then $v_{l}$ now has no outgoing red edge. It follows that the resulting colouring is in fact a red-blue colouring.
\end{proof}

\begin{figure}
\begin{tikzpicture}
\node[blackvertex] at (0,0) (v1) {};
\node[smallwhite] at (-.3,0) (v5) {$x$};
\node[blackvertex] at (1.5,0) (v2) {};

\draw[ultra thick, ->, red] (v1)--(v2) {};

\node[blackvertex] at (0,-1.25) (v3) {};
\node[smallwhite] at (-.3,-1.25) (v6) {$y$};
\node[blackvertex] at (1.5,-1.25) (v4) {};
\draw[ultra thick, ->, blue] (v1)--(v3);
\draw[ultra thick, <-, red] (v4) --(v3);
\node[smallwhite] at (.75, .3) (v6) {$e$};

\begin{scope}[xshift = 3cm]

\node[blackvertex] at (0,0) (v1) {};
\node[smallwhite] at (-.3,0) (v5) {$x$};
\node[blackvertex] at (1.5,0) (v2) {};

\draw[ultra thick, ->, blue] (v1)--(v2) {};

\node[blackvertex] at (0,-1.25) (v3) {};
\node[smallwhite] at (-.3,-1.25) (v6) {$y$};
\node[blackvertex] at (1.5,-1.25) (v4) {};
\draw[ultra thick, <-, red] (v1)--(v3);
\draw[ultra thick, ->, red] (v4) --(v3);
\node[smallwhite] at (.75, .3) (v6) {$e$};
\end{scope}
\end{tikzpicture}

\begin{caption}
\\ An example of an exchange on an edge $e$ and $(x,y)$.
\end{caption}
\end{figure}

To avoid repetitively mentioning it, we will implicitly make use of Observation \ref{flipobservation}. Now we begin impose some structure on our decomposition. First we make an observation which allows us to effectively ignore parent components with cycles (such components will still exist, but for the purposes of our argument we will not need to worry about them).

\begin{obs}
\label{redcyclesaturationlemma}
Let $(x,y)$ be a blue arc such that $R^{x}$ is distinct from $R^{y}$ and further $R^{y}$ is a tree. Then $x$ does not lie in a cycle of $F$.
\end{obs}

\begin{proof}
Suppose towards a contradiction that $x$ lies in a cycle of $F$. Let $e$ be an edge incident to $x$ which lies in the cycle coloured red. Now exchange $(x,y)$ and $e$. As $(x,y)$ was an arc between two distinct red components, and $e$ was in the cycle coloured red, after performing the exchange, we reduce the number of cycles in $F$ by one. However, this contradicts that we picked our counterexample to have the fewest number of red cycles.
\end{proof}

With this we can show that given two components $K$ and $C$, where $K$ is the parent of $C$, and $C$ is acyclic, that $e_{r}(K) + e_{r}(C) \geq d$.

\begin{lemma}
\label{troublesomedonttroublesome}
Let $R^{x}$ and $R^{y}$ be red components such that $R^{y}$ is the child of $R^{x}$, $R^{y}$ does not contain a cycle, and $(x,y)$ is a blue arc from $x$ to $y$. Then $e_{r}(R^{x}) + e_{r}(R^{y}) \geq d$.
\end{lemma}

\begin{proof}
Suppose towards a contradiction that $e_{r}(R^{x}) + e_{r}(R^{y}) < d$. Hence $e_{r}(R^{x}) < d$. Thus $R^{x}$ is not the root.  Let $w$ be a vertex which determines the legal order for $R^{x}$ ($w$ exists as $R^{x}$ is not the root). By Observation \ref{redcyclesaturationlemma} we know that $x$ does not belong to a cycle of $R^{x}$.

\textbf{Case 1:} $w \neq x$.

Let $e$ be the edge incident to $x$ in $R^{x}$ such that $e$ lies on the path from $x$ to $w$ in $R^{x}$. Then exchange $(x,y)$ and $e$. As $e_{r}(R^{x}) + e_{r}(R^{y}) < d$, all resulting red components have less than $d$ edges, and hence we do not increase the residue function. Furthermore, we claim we can find a smaller legal order. Let $R_{i}$ be the component in the legal order corresponding to $R^{x}$. Then consider the new legal order where the components $R_{1},\ldots,R_{i-1}$ remain in the same position, we replace $R_{i}$ with the new red component containing $w$, and then complete the order arbitrarily. By how we picked $e$, $e(R^{w})$ is strictly smaller than $e(R_{i})$, and hence we have found a smaller legal order, a contradiction.

\textbf{Case 2:} $w =x$.

As $R^{x}$ is not the root component, let $R^{x_{n}}$ be an ancestor of $R^{x}$ such that $e(R^{x_{n}}) \geq 1$ (there is an ancestor with this property, as the root has at least one edge). Let $R^{x_{n}},R^{x_{n-1}},\ldots,R^{x_{1}}$ be a sequence of red components such that for $i \in \{1,\ldots,n-1\}$, $R^{x_{i}}$ is the child of $R^{x_{i+1}}$ and $R^{x}$ is the child of $R^{x_{1}}$. Up to relabelling the vertices,  there is a path $P = x_{n},\ldots,x_{1},x,y$ such that $(x_{i+1},x_{i})$ is a arc coloured blue, and $(x_{1},x)$ is a arc coloured blue. Let $e$ be a red edge incident to $x_{n}$. Now do the following. Colour $(x,y)$ red, and reverse the direction of all arcs in $P$. Colour $e$ blue, and orient $e$ away from $x_{n}$. By the same argument in Observation \ref{flipobservation}, the resulting orientation is in $\mathcal{F}$ and the colouring described is a red-blue colouring. Furthermore as $e_{r}(R^{x}) + e_{r}(R^{y}) < d$, all resulting red components have less than $d$ edges, and hence the residue function did not increase. Finally, we can find a smaller legal order in this orientation, as we simply take the same legal order up to the component containing $x_{n}$, and then complete the remaining order arbitrarily. As the component containing $x_{n}$ has at least one less edge now, this order is a smaller legal order, a contradiction.
\end{proof}

\begin{figure}
\begin{tikzpicture}

\node[blackvertex] at (0,0) (v1) {};
\node[blackvertex] at (2,0) (v2) {};
\node[blackvertex] at (4,0) (v3) {};

\node[blackvertex] at (4,-1) (v4) {};
\node[blackvertex] at (2,-1) (v5) {};

\node[blackvertex] at (0,-1) (v6) {};

\node[blackvertex] at (0,-2) (v7) {};

\node[blackvertex] at (2,-2) (v8) {};

\draw[ultra thick,red, ->,>=stealth] (v1)--(v2);
\draw[ultra thick,red, ->, >=stealth] (v2)--(v3);
\draw[ultra thick,blue, ->, >=stealth] (v1)--(v6);
\draw[ultra thick,blue, ->, >= stealth] (v6)--(v7);
\draw[ultra thick,red, ->, >= stealth] (v4)--(v5);
\draw[ultra thick,blue, ->, >= stealth] (v3)--(v4);
\draw[ultra thick,blue, ->, >= stealth] (v2)--(v5);
\draw[ultra thick,blue, ->, >= stealth] (v5)--(v8);
\draw[ultra thick,blue,->, >= stealth] (v8)--(v7);
\draw[ultra thick,blue,->, >= stealth] (v7)--(v2);
\draw[ultra thick,blue,->, >= stealth] (v4)--(v8);

\node[smallwhite] at (-.5,0) (v9) {$R_{1}$};
\node[smallwhite] at (-1,-1) (v10) {$R_{2}$};
\node[smallwhite] at (-.5,-1) (v14) {$x$};
\node[smallwhite] at(-.5,-2) (v15) {$y$};
\node[smallwhite] at (-1,-2) (v11) {$R_{3}$};
\node[smallwhite] at (2.7,-2) (v12) {$R_{5}$};
\node[smallwhite] at (1.55,-1) (v13) {$R_{4}$};

\begin{scope}[xshift = 6cm]
\node[blackvertex] at (0,0) (v1) {};
\node[blackvertex] at (2,0) (v2) {};
\node[blackvertex] at (4,0) (v3) {};

\node[blackvertex] at (4,-1) (v4) {};
\node[blackvertex] at (2,-1) (v5) {};

\node[blackvertex] at (0,-1) (v6) {};

\node[blackvertex] at (0,-2) (v7) {};

\node[blackvertex] at (2,-2) (v8) {};

\draw[ultra thick,blue, ->,>=stealth] (v1)--(v2);
\draw[ultra thick,red, ->, >=stealth] (v2)--(v3);
\draw[ultra thick,blue, <-, >=stealth] (v1)--(v6);
\draw[ultra thick,red, <-, >= stealth] (v6)--(v7);
\draw[ultra thick,red, ->, >= stealth] (v4)--(v5);
\draw[ultra thick,blue, ->, >= stealth] (v3)--(v4);
\draw[ultra thick,blue, ->, >= stealth] (v2)--(v5);
\draw[ultra thick,blue, ->, >= stealth] (v5)--(v8);
\draw[ultra thick,blue,->, >= stealth] (v8)--(v7);
\draw[ultra thick,blue,->, >= stealth] (v7)--(v2);
\draw[ultra thick,blue,->, >= stealth] (v4)--(v8);

\end{scope}
\end{tikzpicture}
\begin{caption}
\\ An illustration of Case $2$ in Lemma \ref{troublesomedonttroublesome}, where $k=1$ and $d=1$ and $x =w$. A vertex in the root component was the ancestor,  so in this case we reduce the residue function.
\end{caption}
\end{figure}

We note the following important special case of Lemma \ref{troublesomedonttroublesome}, that small red components do not have small red children.

\begin{cor}
\label{smallchildrencorollary}
If $K$ is a small red component, then $K$ does not have any small red children.
\end{cor}

\begin{proof}
Suppose towards a contradiction that $K$ has a small red child $C$. As $K$ is small, then $e_{r}(K) < \frac{d}{k+1}$. Similarly, $e_{r}(C) < \frac{d}{k+1}$. But then $e_{r}(K) + e_{r}(C) < \frac{2d}{k+1} \leq d$, contradicting Lemma \ref{troublesomedonttroublesome}.
\end{proof}

Now we will show that every red component has at most $k$ small red children.

\begin{lemma}
\label{boundingthenumberofsmallchildren}
Let $K$ be a red component, then $K$ has at most $k$ small children.
\end{lemma}

\begin{proof}
Suppose towards a contradiction that $K$ has at least $k+1$ distinct small children. Then by the pigeon-hole principle, there are two distinct small children $T_{1}$ and $T_{2}$ such that there are blue arcs $(x,x')$, $(y,y')$ so that $x \neq y$, $x,y \in V(K)$, $x' \in V(T_{1})$ and $y' \in V(T_{2})$. By Observation \ref{redcyclesaturationlemma} we can assume that neither $x$ nor $y$ lie in a red cycle in $K$. Consider a path $P_{x,y}$ in $K$ from $x$ to $y$ (not directed). Let $e_{x}$ be the edge incident to $x$ in $P_{x,y}$ and $e_{y}$ be the edge incident to $y$ in $P_{x,y}$.

Let $K_{x}$ denote the component of $K-e_{y}$ which contains $x$, and let $K_{y}$ denote the component of $K-e_{x}$ which contains $y$.

\begin{claim}
$e_{r}(K_{x}) \leq e_{r}(T_{2})$ and $e_{r}(K_{y}) \leq e_{r}(T_{1})$.
\end{claim}

\begin{proof}
By symmetry, we will only show that $e_{r}(K_{x}) \leq e_{r}(T_{2})$. So suppose towards a contradiction that $e_{r}(K_{x}) > e_{r}(T_{2})$. Then exchange on $(y,y')$ and $e_{y}$. As $e_{r}(K_{x}) > e(T_{2})$, the residue function does not increase. Observe that if $K$ is the root, then the residue function will decrease, and that will give a contradiction. Thus we may assume that $K$ is not the root. We claim we can find a smaller legal order. If there is a vertex which determines the legal order for $K$ in $K_{x}$, then taking the same legal order up to $K$ and then replacing $K$ with $K_{x}$ gives a smaller legal order. Similarly, if there is no vertex which determines the legal order in $K_{x}$, then because $e_{r}(K_{x}) > e_{r}(T_{2})$, taking the same legal order up to $K$ and replacing $K$ with the component containing $y$ of $K-e_{y}$ after the exchange, and filling in the rest of the order arbitrarily gives a smaller legal order. In both cases, this is a contradiction.
\end{proof}

Note that $e_{r}(K)$ either equals $e_{r}(K_{x}) + e_{r}(K_{y})$ or $e_{r}(K) = e_{r}(K_{x}) + e_{r}(K_{y}) -1$ (the first case occurs if $e_{y} \neq e_{x}$, and the second occurs if  $e_{y} = e_{x}$). Note that
\[e_{r}(K_{x}) + e_{r}(K_{y}) < \frac{d}{k+1} + \frac{d}{k+1} \leq d.\]

Hence, $e_{r}(K) \leq d$.  Hence we can assume that $K$ is not the root. Let $w$ be a vertex which determines the legal order. Without loss of generality, we can assume that $w \in V(K_{x})$. Then exchange on $(y,y')$ and $e_{y}$. We do not increase the residue function as $e_{r}(K_{y}) \leq e_{r}(T_{1}) < \frac{d}{k+1}$. However, we can find a smaller legal order by taking the same legal order up to $K$, and replacing $K$ with $K_{x}$, and completing this order arbitrarily. But this contradicts our choice of legal order, a contradiction.
\end{proof}

We are now in position to prove the theorem.

\section{Bounding the maximum average degree}
\label{keyarguments}

In this section, we give a counting argument to show that our chosen exploration subgraph has too large average degree. We make the following definition for ease of notation.

\begin{definition}
Let $K$ be a red component, and let  $K_{1},\ldots,K_{q}$ be the small red children of $K$. We will let $K_{C}$ denote the subgraph with vertex set $V(K_{C}) = V(K) \cup V(K_{1}) \cup \cdots \cup V(K_{q})$, and contains all red edges from $K,K_{1},\ldots,K_{q}$.
\end{definition}

\begin{lemma}
\label{counting}
Let $K$ be a red component which is not small. Then $K_{C}$ is not small. Further, if $e_{r}(K) > d$, then
\[\frac{e_{r}(K_{C})}{v(K_{C})} > \frac{d}{d+k+1}.\]
\end{lemma}

\begin{proof}
First observe that we can assume that $K$ has a small child. If not, then we can apply Observation \ref{AcyclicObservation}, and deduce the claim.

Suppose that $K$ has small children $K_{1},\ldots,K_{q}$. By Lemma \ref{boundingthenumberofsmallchildren} $q \leq k$.  By Lemma \ref{troublesomedonttroublesome}, we know for every $i \in \{1,\ldots,q\}$ the inequality $e_{r}(K) + e_{r}(K_{i}) \geq d$ holds. As $e_{r}(K_{i}) \geq 0$ for all $i \in \{1,\ldots,q\}$, it follows that $e_{r}(K_{i}) \geq \max\{0,d- e(K)\}$ for all $i \in \{1,\ldots,q\}$.

Then a quick calculation shows
\begin{align*}
\frac{e_{r}(K_{C})}{v(K_{C})} &= \frac{e_{r}(K) +\sum_{i =1}^{q} e_{r}(K_{i})}{v(K) + \sum_{i =1}^{q} v(K_{i})}\\
 & \geq \frac{e_{r}(K) +\sum_{i =1}^{q} \max\{0,d-e_{r}(K)\}}{e_{r}(K) + 1 + q + \sum_{i=1}^{q} \max\{0,d-e_{r}(K)\}}
 \\ & \geq \frac{e_{r}(K) +\sum_{i =1}^{q} \max\{0,d-e_{r}(K)\}}{e_{r}(K) + 1 + k + \sum_{i=1}^{q} \max\{0,d-e_{r}(K)\}}.
\end{align*}

The first equality is simply applying the definition of $K_{C}$. The second inequality uses that $e_{r}(K_{i}) \geq \max\{0,d-e_{r}(K)\}$, and that as $K_{i}$ is small, $K_{i}$ is a tree, so $v(K_{i}) = e_{r}(K_{i}) +1$. Finally, the last inequality is using that $q \leq k$. Now we split this into two cases based on whether or not $\max\{0,d-e_{r}(K)\}$ is $0$ or $d-e_{r}(K)$.

\textbf{Case 1:} $\max\{0,d-e(K)\} = 0$.

If $\max\{0,d-e(K)\} =0$, then $e(K) \geq d$. Thus it follows that,
\begin{align*}
\frac{e(K) +\sum_{i =1}^{q} \max\{0,d-e(K)\}}{e(K) + 1 + k + \sum_{i=1}^{q} \max\{0,d-e(K)\}} &= \frac{e(K)}{e(K) + k+1} \\
&\geq \frac{d}{d+k+1}.
\end{align*}
Further, if $e_{r}(K) >d$, the above inequalities are strict.

\textbf{Case 2:} $\max\{0,d-e(K)\} = d-e(K)$.
Note that in this case $e(K) \leq d$. Calculating we obtain,
\begin{align*}
\frac{e(K) +\sum_{i =1}^{q} \max\{0,d-e(K)\}}{e(K) + 1 + k + \sum_{i=1}^{q} \max\{0,d-e(K)\}} &=
\frac{e(K) +q(d - e(K))}{e(K) + q(d-e(K)) +k+1} \\
&\geq \frac{d - e(K) +e(K)}{e(K) +d - e(K) +k+1} \\
&= \frac{d}{d+k+1}.
\end{align*}

Note that in this case $e_{r}(K) \leq d$.
\end{proof}

Now we finish the proof.  Let $\mathcal{R}$ denote the set of red components which are not small. By Corollary \ref{smallchildrencorollary} it follows that,
\[V(H_{f,D,R}) = \bigcup_{K \in \mathcal{R}} V(K_{C}).\]

This follows since a small component cannot have a small child by  Corollary \ref{smallchildrencorollary}. Therefore it follows that:
\[E_{r}(H_{f,D,R}) = \bigcup_{K \in \mathcal{R}} E(K_{C}).\]

Now we bound the maximum average degree of $H_{f,D,R}$. We have
\begin{align*}
\frac{e_{r}(H_{f,D,R})}{v(H_{f,D,R})} &= \frac{\sum_{K \in \mathcal{R}}e_{r}(K_{C})}{\sum_{K \in \mathcal{R}} v(K_{C})} \\
&> \frac{d}{d+k+1}.
\end{align*}

Here, equality holds in the first line because we chose unique parents for components. The strict inequality follows as $K_{C}$ is not small for any $K \in \mathcal{R}$ by Lemma \ref{counting}, and further the root component satisfies $e(R) >d$. However, this contradicts Observation \ref{finishingobservation}. Theorem \ref{strongpndt} follows.

\begin{ack}
The authors are extremely thankful to Joseph Cheriyan for numerous helpful discussions about the problem. We thank the anonymous referees for helpful comments in improving the presentation.
\end{ack}

\bibliographystyle{plain}


Unique determination of ellipsoids by their dual volumes and the moment problem

Abstract: Gusakova and Zaporozhets conjectured that ellipsoids in $\mathbb R^n$ are uniquely determined (up to an isometry) by their Steiner polynomials. Petrov and Tarasov confirmed this conjecture in $\mathbb R^3$. In this paper we solve the dual problem. We show that any ellipsoid in $\mathbb{R}^n$ centered at the origin is uniquely determined (up to an isometry) by its dual Steiner polynomial. To prove this result we reduce it to a problem of moments. As a by-product we give an alternative proof of the result of Petrov and Tarasov.

\section{Introduction}

The study of behavior of volume under the Minkowski (vector) addition is the main focus of the Brunn-Minkowski theory.
Let $K$ be a convex body in $\mathbb R^n$. Denote by $B_2^n$ the unit Euclidean ball in $\mathbb R^n$. The classical Steiner formula asserts that for every $\epsilon > 0$,
 $$\mathrm{vol}(K+\epsilon B_2^n) =  \sum_{i=0}^n\kappa_{n-i} V_i(K) \epsilon^{n-i},$$
where $\mathrm{vol}$ is the Lebesgue measure on $\mathbb R^n$, the addition $+$ is the Minkowski addition, and $\kappa_{n-i}$ is the volume of $B_2^{n-i}$. The coefficients   $V_i(K)$ are known as the intrinsic volumes of $K$. The geometric interpretation of some of these quantities is the following: $V_n(K)$ is the volume of $K$, $V_{n-1}(K)$ is (a multiple of) the surface area, $V_1(K)$ is (a multiple of) the mean width. A thorough discussion of intrinsic volumes can be found in \cite{Sch}.

  Gusakova and Zaporozhets   asked if ellipsoids are uniquely determined by their intrinsic volumes (up to an isometry). Namely, they conjectured the following.
 \begin{conjecture} \label{conj:main1}
 Let $\ce_1$ and $\ce_2$	be two ellipsoids in $\R^n$ such that $  V_1(\ce_1)= V_1(\ce_2)$, $  V_2(\ce_1)= V_2(\ce_2)$,..., $  V_n(\ce_1)= V_n(\ce_2)$. Then $\ce_1$ and $\ce_2$ are congruent.
\end{conjecture}

 Petrov and Tarasov \cite{PT} confirmed this conjecture in $\mathbb R^3$. For higher dimensions, the problem is still open.

The purpose of this paper is to show that a similar question for dual volumes has a positive answer in all dimensions. Dual volumes were introduced by Lutwak \cite{L} within the framework of the  dual Brunn-Minkowski theory. In this theory the Minkowski addition of convex bodies is replaced by the radial addition of star bodies (see the next section for definitions of these concepts). The dual version of the Steiner formula  asserts that
$$\mathrm{vol}(K\,\tilde{+} \,\epsilon B_2^n) =   \sum_{i=0}^n{ n\choose i} \widetilde V_i(K) \epsilon^{n-i},$$
where $K$ is a star body in $\mathbb R^n$  and $\tilde{+}$ is the radial addition.
The coefficients   $\widetilde V_i(K)$ are called the  dual   volumes. Note that $\widetilde V_n(K)$ is equal to the volume of $K$.
Denoting by $\rho_K$ the radial function of $K$, one can write the dual volumes of $K$ as follows:
\begin{equation}\label{tildeVi}
\widetilde V_i (K) = \frac{1}{n} \int_{S^{n-1}} \rho_K^i(\theta) \, d\theta,
\end{equation}
where the integration is with respect to the spherical Lebesgue measure.

 Note that   while the intrinsic volumes are invariant under translations, the dual volumes depend on the choice of the origin. Both the intrinsic volumes and dual volumes are invariant under orthogonal transformations.

Our main result is the following.

\begin{theorem} \label{thm:main1}
	Let $\ce_1$ and $\ce_2$	be two ellipsoids in $\R^n$, $n\ge 2$, centered at the origin such that $\widetilde V_1(\ce_1)=\widetilde V_1(\ce_2)$, $\widetilde  V_2(\ce_1)= \widetilde V_2(\ce_2)$,..., $\widetilde  V_n(\ce_1)=\widetilde V_n(\ce_2)$. Then $\ce_1$ and $\ce_2$ are congruent.
\end{theorem}

As one can see, the right-hand side of formula (\ref{tildeVi}) makes sense for all real $i$.  This allows to use  (\ref{tildeVi}) as a definition of dual volumes of any order $i$.
In view of this remark, in the statement of Theorem \ref{thm:main1} the collection of the dual volumes $\{ \widetilde V_{i}  \}_{i=1}^n$can be replaced by any $n$-tuple  of the form $\{ \widetilde V_{i_k} \}_{k=1}^n$, where $i_1$,..., $i_n$ are distinct non-zero real numbers from the interval $(-2,n]$. In some cases one can take numbers from a larger interval, as will be discussed later.

To prove Theorem \ref{thm:main1} we reduce it to a   problem of moments. Using this idea we also give an alternative proof of the result of Petrov and Tarasov for the intrinsic volumes of ellipsoids in $\mathbb R^3$. In the final section we investigate Conjecture \ref{conj:main1} in the setting of ellipsoids of revolution.

\section{Preliminaries}
Let $K$ be a convex body in $\mathbb R^n$, that is a convex compact set with non-empty interior. The support function of $K$ is defined by
$$h_K(x) = \max_{y\in K} \langle x,y\rangle,\qquad x\in \mathbb R^n. $$
The first intrinsic volume of $K$ is a multiple of its mean width:
\begin{equation}\label{1stIV}
V_1(K) = \frac{1}{\kappa_{n-1}} \int_{S^{n-1}} h_K(\theta) d\theta,
\end{equation}
where the integration is with respect to the spherical Lebesgue measure.

Let  $\ce$ be an ellipsoid in $\mathbb R^n$ with semi-axes $a_1, \dots, a_n$. Its volume is given by
$$\mathrm{vol}(\ce) = \kappa_n a_1  \cdots  a_n.$$  Below we will always assume that $\ce$ is centered at the origin and its axes coincide with the coordinate axes; i.e.,
\begin{equation}\label{ce}
\ce=\left\{x\in \mathbb R^n: \frac{x_1^2}{a_1^2} + \cdots + \frac{x_n^2}{a_n^2}\le 1\right\}.
\end{equation}
The support function of this ellipsoid $\ce$ is given by
$$h_{\ce} (x) = \left({a_1^2} {x_1^2} + \cdots + {a_n^2}  {x_n^2}\right)^{1/2}, \qquad x\in \mathbb R^n.$$

For a convex body $K\subset \mathbb R^n$ its polar body $K^\circ$ is defined as follows:
$$K^\circ = \{ x\in \mathbb R^n: \langle x,y\rangle \le 1, \forall y \in K\}.$$
For the ellipsoid $\ce$ defined above by (\ref{ce}), its polar body is also an ellipsoid given by
$$\ce^\circ =\left\{x\in \mathbb R^n: {a_1^2} {x_1^2} + \cdots + {a_n^2}  {x_n^2}\le 1\right\}.$$

We say that   a compact set   $K\subset \mathbb R^n$   is star-shaped about the origin $o$ if for every point $x\in K$ each point of the interval $[o,x)$ is an interior point of $K$.
The {\it Minkowski functional} of $K$ is defined by
\begin{align*}
\|x\|_K=\min \{a\ge 0: x \in aK \}, \qquad x\in \mathbb R^n.
\end{align*}
We say that $K$ is a {\it star body} if it is compact, star-shaped about the origin and its Minkowski functional  is a  continuous function on $\mathbb R^n$.

The Minkowski functional of   ellipsoid (\ref{ce}) is given by
$$\|x\|_{\ce} = \left(\frac{x_1^2}{a_1^2} + \cdots + \frac{x_n^2}{a_n^2}\right)^{1/2}, \qquad x \in \mathbb R^n.$$
For the Euclidean ball $B_2^n$, we will denote its Minkowski functional simply by $|x|$.

If $K$ is a convex body containing the origin in its interior then we have the following relation:
$$\|x\|_K = h_{K^\circ}(x), \qquad x\in \mathbb R^n.$$

The {\it radial function} of a star body  $K$ is defined by \begin{align*}
\rho_K(\xi) = \max\{a \ge 0 :a \xi \in K\}, \quad \xi\in S^{n-1}.
\end{align*}
 Clearly,   $\rho_K(\xi) =  \|\xi\|_K^{-1}$ for $\xi\in S^{n-1}$.

If $K$ and $L$ are star bodies, and $\alpha$ and $\beta$ are positive numbers, then the radial sum $\alpha K\,\tilde{+}\, \beta L$ is  the star body with the radial function
$$\rho_{\alpha K\tilde{+} \beta L} = \alpha \rho_{K}+\beta \rho_{L}.$$

Let us now recall some basic facts about the Gamma-function and fractional derivatives. If    $z\in \mathbb C$ with $\Re(z)>0$, then
$$\Gamma (z) = \int_0^\infty t^{z-1} e^{-t}\, dt.$$
The Gamma-function can be extended to the set $ \mathbb C\setminus (-\mathbb N\cup\{0\})$ using the relation $$\Gamma (z+1)=z\Gamma (z).$$

Let $f$ be a continuous integrable function on $\mathbb [0,\infty)$ which is infinitely smooth in some neighborhood of zero.   If $q\in \mathbb C$, $-1<\Re(q)<0$, the fractional derivative of $f$ of order $q$ at zero is defined by
$$f^{(q)}(0) = \frac{1}{\Gamma(-q)}\int_0^\infty t^{-1-q} f(t)\,dt.$$
This formula is in fact valid for all  $q\in \mathbb C$, $\Re(q)<0$, if $ t^{-1-q} f(t)$ decays at infinity sufficiently fast.

For $q$ with  $\Re(q)>0$ that are not positive integers, the fractional derivatives are defined using analytic continuation;  see, e.g., \cite{GS} or \cite{K}.
In particular, if $m$ is a positive integer, then for $m-1<\Re(q)<m$ we have
$$f^{(q)}(0) = \frac{1}{\Gamma(-q)}\int_0^\infty t^{-1-q} \left(f(t)- \sum_{k=0}^{m-1} \frac{f^{(k)}(0)}{k!}t^k\right)\,dt.$$
Fractional derivatives of positive integer orders coincide with ordinary derivatives (up to a sign).

The reader is referred to the books \cite{G}, \cite{K}, \cite{Sch} for other facts on the topics discussed above.

\section{Main results}\label{SectionMR}

It is known that the intrinsic volumes of an ellipsoid and its polar are related by the following formula.
\begin{equation}\label{KZ} V_i(\ce) = \frac{\kappa_{i}}{\kappa_n\kappa_{n-i}}V_n(\ce) V_{n-i}(\ce^\circ);
\end{equation}
see \cite[Prop. 4.8]{KZ}.

Let us prove a similar formula for the dual volumes.

\begin{lemma}
	Let $\ce$ be a centered ellipsoid in $\mathbb R^n$. Then for all real $i$ we have
	 \begin{equation}\label{relation_dual}
	\widetilde{V}_i(\ce) = \frac{1}{\kappa_n} \widetilde{V}_n(\ce)\widetilde{V}_{n-i}(\ce^\circ).
	\end{equation}
\end{lemma}
\begin{proof}
 The formula is trivial when $i=n$, so we will assume that $i\ne n$.
 We have
 	\begin{align*}
 \widetilde{V}_i(\ce) &= \frac{1}{n}\int_{S^{n-1}} \|\theta\|_{\ce}^{-i} \, d\theta  = \frac{n-i}{n(2^{n-i}-1)}\int_{S^{n-1}} \int_1^2 r^{n-i-1} \, dr \, \|\theta\|_{\ce}^{-i} \, d\theta\\
 & = \frac{n-i}{n(2^{n-i}-1)}\int_{2 B_2^{n}\setminus B_2^n}  \|x\|_{\ce}^{-i} \, dx.
 \end{align*}
 Assuming without loss of generality that $$\|x\|_{\ce} = \left(\frac{x_1^2}{a_1^2} + \cdots + \frac{x_n^2}{a_n^2}\right)^{1/2}, \qquad x \in \mathbb R^n,$$
 and making a  change of variables, we get
	\begin{align*}   \widetilde{V}_i(\ce) &=   \frac{n-i}{n(2^{n-i}-1)}\int_{2 B_2^{n}\setminus B_2^n}  \left(\frac{x_1^2}{a_1^2} + \cdots + \frac{x_n^2}{a_n^2}\right)^{-i/2} \, dx\\
	&= a_1\dots a_n \frac{n-i}{n(2^{n-i}-1)}\int_{2 \ce^\circ\setminus \ce^\circ}  \left( {x_1^2}  + \cdots +  {x_n^2} \right)^{-i/2} \, dx\\
		&= a_1\dots a_n \frac{n-i}{n(2^{n-i}-1)} \int_{S^{n-1}}
	 \int_{\|\theta\|_{\ce^\circ}^{-1}}^{2\|\theta\|_{\ce^\circ}^{-1}} r^{n-i-1} \, dr \, d\theta\\
	 	&= a_1\dots a_n  \frac{1}{n}\int_{S^{n-1}}
	 \|\theta\|_{\ce^\circ}^{-n+i} \, d\theta\\
	 &=\frac{1}{\kappa_n} \widetilde{V}_n(\ce)\widetilde{V}_{n-i}(\ce^\circ).
  \end{align*}
\end{proof}

We will now obtain a representation of dual   volumes as certain moments.

\begin{lemma}\label{MainFormula}
	Let $\ce$ be an ellipsoid centered at the origin with semi-axes $a_1$,..., $a_n$. If $i$ is a  real number  (not necessarily an integer) such that $0<i<n$, then
	\begin{equation}  \label{FormulaW_i1}
	\widetilde{V}_i(\ce) = \frac{4 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \int\limits_0^\infty \frac{u^{i-1}}{\sqrt{(1 + \frac{u^2}{a_1^2})\cdots (1 + \frac{u^2}{a_n^2})}}\, du.
	\end{equation}

	If $-2<i<0$, then

		\begin{equation}  \label{FormulaW_i2}
	\widetilde{V}_i(\ce) = \frac{4 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \int\limits_0^\infty u^{i-1}\left( \frac{1}{\sqrt{(1 + \frac{u^2}{a_1^2})\cdots (1 + \frac{u^2}{a_n^2})}}-1\right) \, du.
	\end{equation}

	If $n<i<n+2$, then
	\begin{equation}  \label{FormulaW_i3}	\widetilde{V}_i(\ce) = \frac{4 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \int\limits_0^\infty u^{i-1} \left(\frac{1}{\sqrt{(1 + \frac{u^2}{a_1^2})\cdots (1 + \frac{u^2}{a_n^2})}}-\frac{a_1\cdots a_n}{u^n}\right) \, du  .
	\end{equation}

\end{lemma}
\begin{proof} Applying an orthogonal transformation if needed, we can assume that the Minkowski functional of $\ce$ is given by
	$$\|x\|_{\ce} = \left(\frac{x_1^2}{a_1^2} + \cdots + \frac{x_n^2}{a_n^2}\right)^{1/2}, \qquad x \in \mathbb R^n.$$
Recall that for all $i$ we have
		\begin{equation*}\label{int}
		\widetilde{V}_i(\ce) = \frac{1}{n}\int_{S^{n-1}} \|\theta\|_{\ce}^{-i} \, d\theta.
		\end{equation*}

			First we will consider the case $i\in (0,n)$. Observe that
	\begin{equation}\label{repr}
	\int_{S^{n-1}} \|\theta\|_{\ce}^{-i} d\theta = \frac{2}{ \Gamma\left(\frac{n-i}2\right) } \int_{\R^n} \|x\|_{\ce}^{-i} e^{-|x|^2} \, dx.
	\end{equation}
	To check this, pass to the polar coordinates in the latter integral:
	$$ \int_{\R^n}  \|x\|_{\ce}^{-i} e^{-|x|^2} \, dx =  \int_{S^{n-1}}  \|\theta\|_{\ce}^{-i} \, d\theta \int_0^\infty r^{n-i-1} e^{-r^2} \, dr = \frac12 \Gamma\left(\frac{n-i}2\right)\int_{S^{n-1}} \|\theta\|_{\ce}^{-i}\, d\theta.$$

	Thus, we have
	$$\widetilde{V}_i(\ce) = \frac{2}{n \Gamma\left(\frac{n-i}2\right) } \int_{\R^n}   \|x\|_{\ce}^{-i}   e^{-|x|^2} \, dx.
	$$

Now observe that, for $i>0$, we have
$$\|x\|_{\ce}^{-i}  = \frac{2}{\Gamma\left(\frac{i}{2}\right)} 	\int\limits_0^\infty s^{i-1}  e^{- \|x\|_{\ce}^{2}\, s^2}\, ds.$$
This can be checked by a   change of the variable in the latter integral.

	Therefore,
	\begin{align*}
	\widetilde{V}_i(\ce)& = \frac{4}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) } 	\int\limits_0^\infty  s^{i-1} \int_{\R^n}  e^{-\left(\frac{x_1^2}{a_1^2} + \cdots + \frac{x_n^2}{a_n^2}\right)\, s^2}   e^{-|x|^2} \, dx\, ds \\
	& = \frac{4}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }  	\int\limits_0^\infty s^{i-1} \int\limits_{\R}  e^{-x_1^2 \left(1 + \frac{s^{2}}{a_1^2} \right) }  dx_1 \cdots\int\limits_{\R} e^{-x_n^2 \left(1 + \frac{s^{2}}{a_n^2 } \right)}  dx_n \,  ds \\
		& = \frac{4}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }  	\int\limits_0^\infty s^{i-1} \frac{\sqrt{\pi}}{\sqrt{1+\frac{s^2}{a_1^2} }}    \cdots \frac{\sqrt{\pi}}{\sqrt{1 + \frac{s^2}{a_n^2 }}} \, ds ,
    \end{align*}
    as claimed.

   To treat the case $i<0$, we will use the analytic continuation technique.
   Using the formula for $\widetilde{V}_i(\ce)$, $0<i<n$, proved above and making a change of variables, we get
   $$	\widetilde{V}_i(\ce) = \frac{2 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \int\limits_0^\infty \frac{u^{i/2-1}}{\sqrt{(1 + \frac{u }{a_1^2})\cdots (1 + \frac{u }{a_n^2})}}\, du.$$
    Denoting
    $$f_{\ce}(u) =  \frac{1}{\sqrt{(1 + \frac{u}{a_1^2})\cdots (1 + \frac{u}{a_n^2})}},$$
    we see that
  \begin{equation}\label{frdir}\widetilde{V}_i(\ce) = \frac{2 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)  }  f_{\ce}^{(-i/2)}(0),
  \end{equation}
    where $f_{\ce}^{(-i/2)}(0)$ is the fractional derivative of $f$ of order $-i/2$ at zero.

    The left hand-side of \eqref{frdir} is an analytic function of $i\in \mathbb C$ and the right-hand side is analytic for in $\{i\in \mathbb C : \Re (i) < n\}$. Since they coincide on the interval $i\in (0,n)$, they coincide for all $i<n$. Thus formula \eqref{frdir} is valid for all $i<n$.
    In particular, for $-2<i<0$, we have
    $$	\widetilde{V}_i(\ce) = \frac{2 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \int\limits_0^\infty u^{i/2-1} \left( \frac{1}{\sqrt{(1 + \frac{u }{a_1^2})\cdots (1 + \frac{u }{a_n^2})}}-1\right) \, du.$$
    Replacing $u$ by $u^2$, we get the formula given in the statement of the lemma.

 To deal with the case when $i>n$,
 we will first use \eqref{relation_dual} and then apply  \eqref{frdir} with  $n-i$ instead of $i$ and with $\ce^\circ$ instead of $\ce$:

 $$\widetilde{V}_i(\ce) = a_1\cdots a_n \widetilde{V}_{n-i}(\ce^\circ) = a_1\cdots a_n \frac{2 \pi^{n/2}}{n \Gamma\left(\frac{i}2\right)  }  f_{\ce^\circ}^{((i-n)/2)}(0), $$
 where
 $$f_{\ce^\circ}(u) =  \frac{1}{\sqrt{(1 +  {u}{a_1^2})\cdots (1 +  {u}{a_n^2})}}.$$

    In particular, for $n<i<n+2$, we have
  \begin{align*}	\widetilde{V}_i(\ce) &= \frac{2 \pi^{n/2}a_1\cdots a_n}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \int\limits_0^\infty u^{(n-i)/2-1} \left( \frac{1}{\sqrt{(1 +  {u }{a_1^2})\cdots (1 +  {u }{a_n^2})}}-1\right) \, du\\
   &= \frac{4 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \int\limits_0^\infty u^{i-1} \left(\frac{1}{\sqrt{(1 + \frac{u^2}{a_1^2})\cdots (1 + \frac{u^2}{a_n^2})}}-\frac{a_1\cdots a_n}{u^n}\right) \, du  .
  \end{align*}

\end{proof}

Now we are ready to prove our main result.

{\bf Proof of Theorem \ref{thm:main1}.} Let $\ce_1$ be an ellipsoid with semi-axes $a_1, \ldots, a_n$, and $\ce_2$   an ellipsoid with semi-axes $b_1, \ldots, b_n$. Assume that
$\widetilde V_i(\ce_1)=\widetilde V_i(\ce_2)$, for all $i=1,\ldots,n$.

Note that  $\widetilde  V_n(\ce_1)=\widetilde V_n(\ce_2)$ is the equality  of the volumes of  $\ce_1$ and $\ce_2$, which  gives $a_1\cdots a_n= b_1\cdots b_n$.

By Lemma \ref{MainFormula}, from the equalities of the dual volumes of orders from $1, \dots, n-1$ we get
$$	  \int\limits_0^\infty  \frac{u^{i-1}}{\sqrt{(1 + \frac{u^2}{a_1^2})\cdots (1 + \frac{u^2}{a_n^2})}}\, du =   \int\limits_0^\infty  \frac{u^{i-1}}{\sqrt{(1 +   \frac{u^2}{b_1^2})\cdots (1 + \frac{u^2}{b_n^2})}}  \, du ,$$
for all  $i=1,\ldots,n-1$.

Denoting for brevity
\begin{equation}\label{F}
F(u) =   \frac{1}{\sqrt{(1 + \frac{u^2}{a_1^2})\cdots (1 + \frac{u^2}{a_n^2})}} - \frac{1}{\sqrt{(1 +   \frac{u^2}{b_1^2})\cdots (1 + \frac{u^2}{b_n^2})}},
\end{equation}
we get
$$	  \int\limits_0^\infty  {u^{i-1}} F(u)\, du =0,$$
for all  $i=1,\ldots,n-1$.
This means that for any polynomial $P$ of degree at most $n-2$ we have
$$	  \int\limits_0^\infty P(u) F(u) \, du=0.$$

Now observe that the function $F$
has at most $n-2$ positive real roots, unless it is identically equal to zero. Indeed, solving $F(u)=0$ is equivalent to solving
 $$ \left(1 + \frac{u^2}{a_1^2}\right)\cdots \left(1 + \frac{u^2}{a_n^2}\right)=\left(1 +   \frac{u^2}{b_1^2}\right)\cdots \left(1 + \frac{u^2}{b_n^2}\right) ,$$
 which can be written in the form
 $$
 1+u^2 q(u^2) +\frac{u^{2n}}{a_1^2\cdots a_n^2} = 1 + u^2 r(u^2) + \frac{u^{2n}}{b_1^2\cdots b_n^2},
 $$
 where $q$ and $r$ are polynomials of degree $n-2$.

 Since   $\ce_1$ and $\ce_2$ have equal volumes, the latter equation reduces to $q(u^2) = r(u^2)$, which has at most $n-2$ positive real roots.

 Now choose such a non-zero polynomial $P$ that changes sign exactly at those positive numbers where $F$ changes its sign. This guarantees that the product $PF$ is either everywhere non-negative or non-positive. Since the integral of $PF$ is zero, we conclude that $F$ is identically zero. This means that the sets of numbers $\{a_1, \ldots, a_n\}$ and  $\{b_1, \ldots, b_n\}$ coincide. The theorem is proved.

\qed

We will now explain how to adjust the proof above to the case when dual volumes of other orders are given. Let us take any non-zero numbers $i_k$, $k=1,\dots,n$ from the interval $(-2,n]$. Using either formula \eqref{FormulaW_i1} or \eqref{FormulaW_i2}, the equality  $\widetilde  V_{i_k}(\ce_1)=\widetilde V_{i_k}(\ce_2)$ gives
\begin{equation}\label{eqns}
  \int\limits_0^\infty  {u^{i_k-1}} F(u)\, du =0,
 \end{equation}
if $i_k\ne n$. Here, $F$ is the same function as in \eqref{F}.

If some $i_k$ is equal to $n$ then $F$ has at most $n-2$ positive real roots, and we proceed as above with the help of the following lemma.

\begin{lemma}[Lemma 19, \cite{ENT}] \label{lem:moments}
	Let $f:[0,\infty)\to \mathbb R$ be a continuous function that changes sign at most $N-1$ times in the interval   $(0, \infty)$. If there exist $N$ real numbers $p_1$,..., $p_N$ such that
	\begin{align*}
  \int\limits_0^{\infty} t^{p_k} f(t) dt =0,\quad \textrm{for every} \quad k=1, \ldots, N,
	\end{align*}
	 then $f$ is identically equal to zero.
\end{lemma}

If neither of $i_k$ is equal to $n$, then we can only conclude that $F$ has at most $n-1$ positive real roots, but we have $n$ conditions of the form \eqref{eqns}, so we can still use the above lemma.

To show how to extend Theorem \ref{thm:main1} to a larger set of indices, let us first give a few formulas that follow from the proof of Lemma \ref{MainFormula}.
For $-4<i<-2$ formula \eqref{frdir} yields \pagebreak
\begin{multline*}	\widetilde{V}_i(\ce) = \frac{4 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }   \\
	\times  \int\limits_0^\infty u^{i -1} \left( \frac{1}{\sqrt{(1 + \frac{u^2 }{a_1^2})\cdots (1 + \frac{u^2}{a_n^2})}}-1+\frac12\left(\frac{1 }{a_1^2}+\cdots+\frac{1 }{a_n^2}\right)u^2 \right) \, du.
	\end{multline*}
When $i=-2$, we have
$$	\widetilde{V}_{-2}(\ce) = \frac{ \pi^{n/2}}{n \Gamma\left(\frac{n+2}2\right)  }  \left(\frac{1 }{a_1^2}+\cdots+\frac{1 }{a_n^2}\right) .$$

For $n+2<i<n+4$ we get
\begin{multline*}	\widetilde{V}_i(\ce) = \frac{4 \pi^{n/2}}{n \Gamma\left(\frac{n-i}2\right)\Gamma\left(\frac{i}{2}\right) }  \\ \times \int\limits_0^\infty u^{i-1} \left( \frac{1}{\sqrt{(1 +  \frac{u^2 }{a_1^2})\cdots (1 +  \frac{u^2}{a_n^2})}}-\frac{a_1\cdots a_n}{u^n}+\frac{a_1\cdots a_n\left( {a_1^2}+\cdots+ {a_n^2}\right)}{2u^{n+2}}\right) \, du.
 \end{multline*}
When $i=n+2$, we have
$$	\widetilde{V}_{n+2}(\ce) = \frac{ \pi^{n/2}  a_1\cdots a_n}{n \Gamma\left(\frac{n+2}2\right)  }  \left( {a_1^2}+\cdots+ {a_n^2}\right) .$$

Now observe that in some cases we can take the interval $(-2,n+2)$ instead of the interval $(-2,n]$. If we assume that $\widetilde  V_n(\ce_1)=\widetilde V_n(\ce_2)$ and $\widetilde  V_{i_k}(\ce_1)=\widetilde V_{i_k}(\ce_2)$ for some distinct numbers $i_k\in (-2,n+2)\setminus\{0,n\}$, $k=1,\ldots,n-1$, then along with the formulas \eqref{FormulaW_i1} and \eqref{FormulaW_i2}, we can also use \eqref{FormulaW_i3}, since $\widetilde  V_i(\ce_1)-\widetilde V_i(\ce_2)=0$ with $n<i<n+2$ gives
$$	  \int\limits_0^\infty  {u^{i-1}} F(u)\, du =0,$$ with the same function $F$ as used above.

Using this argument one can proceed even further. For example, if 	$\widetilde{V}_{-2}(\ce_1) = \widetilde{V}_{-2}(\ce_2)$ and $\widetilde{V}_{n}(\ce_1) = \widetilde{V}_{n}(\ce_2)$, then we can take the remaining $n-2$ indices from the interval $(-4,n+2)$. We will have $n-2$ conditions of the form \eqref{eqns}, while the function $F$ will have at most positive $n-3$ roots. To explain the latter, observe that the equation
$$ \left(1 + \frac{u^2}{a_1^2}\right)\cdots \left(1 + \frac{u^2}{a_n^2}\right)=\left(1 +   \frac{u^2}{b_1^2}\right)\cdots \left(1 + \frac{u^2}{b_n^2}\right) $$
  can be written in the form
\begin{multline*}
1+u^2\left( \frac{1}{a_1^2}+\cdots  + \frac{1}{a_n^2}\right)+ u^4 q(u^2) +\frac{u^{2n}}{a_1^2\cdots a_n^2} \\= 1 + u^2\left( \frac{1}{b_1^2}  +\cdots+ \frac{1}{b_n^2}\right)+u^4 r(u^2) + \frac{u^{2n}}{b_1^2\cdots b_n^2},
\end{multline*}
where $q$ and $r$ are polynomials of degree $n-3$.

Using similar ideas, one can get uniqueness results involving $\widetilde{V}_{-4}$, $\widetilde{V}_{n+2}$, and so forth.

\section{Intrinsic volumes of ellipsoids in $\mathbb R^3$}

In this section we will give an alternative proof of the following result of Petrov and Tarasov \cite{PT}.

\begin{theorem}
	Let    $\ce_1$ and $\ce_2$	be two ellipsoids in $\R^3$ such that $  V_1(\ce_1)= V_1(\ce_2)$, $  V_2(\ce_1)= V_2(\ce_2)$,  $  V_3(\ce_1)= V_3(\ce_2)$. Then $\ce_1$ and $\ce_2$ are congruent.
\end{theorem}
\begin{proof}
Let  $\ce$  be an ellipsoid in $\mathbb R^3$. Using formulas \eqref{1stIV}, \eqref{tildeVi}, and \eqref{relation_dual}, we can write its first intrinsic volume as follows:
$$
V_1(\ce ) = \frac{1}{\pi} \int \limits_{S^{2}} h_{\ce }(\theta) d\theta = \frac{1}{\pi} \int \limits_{S^{2}} \|\theta\|_{\ce ^\circ}  d\theta = \frac{3}{\pi}  \widetilde V_{-1}(\ce ^\circ)  = \frac{4}{V_3(\ce)} \widetilde V_{4}(\ce) .
$$
To find the second intrinsic volume of $\ce$, we will use formulas \eqref{KZ}, \eqref{1stIV}, and \eqref{tildeVi} .
\begin{align*}
V_2(\ce)& = \frac{3}{8} V_3(\ce) \cdot V_1(\ce^\circ) =    \frac{3}{8\pi}V_3(\ce) \int_{S^2} h_{\ce^\circ} (\theta) \, d\theta\\
&=  \frac{3}{8\pi}V_3(\ce) \int_{S^2} \|\theta\|_{\ce}   \, d\theta=  \frac{9}{8\pi} V_3(\ce)  \widetilde V_{-1}(\ce).
\end{align*}

Thus the knowledge of $V_1$, $V_2$, $V_3$ is equivalent to the knowledge of $ \widetilde V_{-1}$, $ \widetilde V_4$, $ \widetilde V_3$. The latter information determines the ellipsoid uniquely, as  was shown in the previous section.
\end{proof}

\section{Intrinsic volumes of ellipsoids of revolution in $\R^n$}

In this section we will discuss uniqueness of ellipsoids of revolution with given intrinsic volumes. As we will see below, we do not need to know all intrinsic volumes $V_1, \ldots, V_n$, only a few of them.

Let $K$ be an origin-symmetric convex body. $K$ is called a zonoid if $K$ can be approximated in the Hausdorff metric by finite Minkowski sums of segments.  It is known (see \cite[Thm 3.5.3]{Sch}) that $K$ is a zonoid if and only if
$$
h_{K}(x) = \int \limits_{S^{n-1}} |\langle x , \theta\rangle | d \mu(\theta), \quad x \in \mathbb R^{n},
$$
where $\mu$ is a  non-negative even Borel measure on $S^{n-1}$.   The measure $\mu$ is called the generating measure of $K$.

Let $\ce$ be an ellipsoid in $\mathbb R^n$ centered at the origin. It is easy to see that $\ce$ is a zonoid. Moreover, its generating measure is given by
$$d\mu(\theta) =\frac{\kappa_n}{2 \mathrm{vol}(\ce) \kappa_{n-1}} \|\theta\|_{\ce}^{-n-1} d\theta.$$
Indeed, let $A$ be a  linear transformation such that $\ce = A B_2^n$ . Then
\begin{align*} \int \limits_{S^{n-1}} |\langle x , \theta\rangle |\, \|\theta\|_{\ce}^{-n-1} d\theta & = (n+1) \int \limits_{\ce} |\langle x , y\rangle |\, dy = (n+1) |\det A|  \int \limits_{B_2^n} |\langle x , A y\rangle |\, dy
\\
&= (n+1) |\det A|  \int \limits_{B_2^n} |\langle A^T x ,  y\rangle |\, dy \\
& = (n+1) |\det A| \,|A^T x| \int \limits_{B_2^n} \left|\left\langle \frac{A^T x}{|A^T x|} ,  y\right\rangle \right|\, dy \\
& = (n+1) |\det A| \, |A^T x| C_n,
\end{align*}
where $$C_n=  \int \limits_{B_2^n} \left|\left\langle u ,  y\right\rangle \right|\, dy = \frac{1}{n+1} \int \limits_{S^{n-1}} \left|\left\langle u ,  \theta \right\rangle \right|\, d\theta = \frac{2 \kappa_{n-1}}{n+1} ,$$
for any $u\in S^{n-1}$; see e.g., \cite[Lemma 3.4.5]{Gr}.

Thus,
$$\int \limits_{S^{n-1}} |\langle x , \theta\rangle |\, \|\theta\|_{\ce}^{-n-1} d\theta  = 2 \kappa_{n-1} |\det A|\, |A^T x|  = 2 \frac{\kappa_{n-1}}{\kappa_n} \mathrm{vol}(\ce)  h_{A B_2^n}(x)     =  2  \frac{\kappa_{n-1}}{\kappa_n} \mathrm{vol}(\ce) h_{\ce }(x) ,
$$
as claimed.

Recall (see formulas (5.31) and (5.82) in \cite{Sch}) that   the $k$-th intrinsic volume of a zonoid $K$ with generating measure $\mu$ is given by
\begin{equation}\label{fml:zonoids}
V_k(K) = C(n,k) \int \limits_{S^{n-1}} \ldots \int \limits_{S^{n-1}} D_n(v_1, \ldots, v_n)\, d\mu(v_1)  \cdots d \mu(v_k) \, dv_{k+1} \cdots dv_n,
\end{equation}
where $D_n(v_1, \ldots,v_n)$ is the  volume of the parallelepiped spanned by the vectors $v_1, \ldots, v_n$, and $C(n,k)= \displaystyle\frac{2^k {n\choose k}}{n!\, \kappa_{n-k} \kappa_{n-1}^{n-k}}$.

Therefore, for an ellipsoid $\ce = AB_2^n$ we have
\begin{align*}
V_k(\ce) &= \frac{\kappa^k_n}{2^k \mathrm{vol}(\ce)^k \kappa_{n-1}^k}C(n,k)  \int \limits_{S^{n-1}} \ldots \int \limits_{S^{n-1}} D_n(v_1, \ldots, v_n)\, \|v_1\|_{\ce}^{-n-1} d v_1   \cdots \|v_k\|_{\ce}^{-n-1} d  v_k  \, dv_{k+1} \cdots dv_n\\
 &= \frac{C(n,k)(n+1)^n \kappa^k_n}{2^k \mathrm{vol}(\ce)^k \kappa_{n-1}^k}  \underbrace{\int \limits_{\ce} \ldots \int \limits_{\ce} }_{k} \underbrace{\int \limits_{B_2^n} \ldots \int \limits_{B_2^n}}_{n-k} D_n(v_1, \ldots, v_n)\,   d v_1   \cdots dv_n\\
  &=\frac{C(n,k)(n+1)^n }{2^k  \kappa_{n-1}^k}  \int \limits_{B_2^n} \ldots \int \limits_{B_2^n}  D_n(A v_1, \ldots, A v_k,   v_{k+1}, \ldots  v_n)\,   d v_1   \cdots dv_n\\
   &= \frac{C(n,k) }{2^k  \kappa_{n-1}^k}  \int \limits_{S^{n-1}} \ldots \int \limits_{S^{n-1}}  D_n(A v_1, \ldots, A v_k,  v_{k+1}, \ldots  v_n)\,   dv_1   \cdots dv_n.
\end{align*}

Now assume that $\ce$ is an ellipsoid of revolution of the form
$$\ce=\left\{x\in \mathbb R^n: \frac{x_1^2}{a^2} + \cdots \frac{x_{n-1}^2}{a^2}+ \frac{x_n^2}{b^2}\le 1\right\}.$$
The corresponding matrix $A$ is diagonal, with the entries  $a, \ldots, a, b$ on the main diagonal. It will be convenient to use the following notation: $A = \mathrm{diag} \, \{a, \ldots, a, b\}$.
We have
\begin{align*} D_n(A v_1, \ldots, A v_k,   v_{k+1} \ldots  v_n)& = \mathrm{abs} \left| \begin{array}{llllll}  a v_{1}^{1} &  \ldots & a v_{k}^{1} & v_{k+1}^{1}
& \ldots &  v_{n}^{1}\\
 a v_{1}^{2} &  \ldots & a v_{k}^{2} & v_{k+1}^{2}
& \ldots &   v_{n}^{2}\\
\vdots  & \ldots & \vdots  &\vdots
& \ldots & \vdots \\
 a v_{1}^{n-1}& \ldots & a v_{k}^{n-1} & v_{k+1}^{n-1}
& \ldots &   v_{n}^{n-1}\\
 b v_{1}^{n}& \ldots & b v_{k}^{n} & v_{k+1}^{n}
& \ldots &  v_{n}^{n}
\end{array} \right| \\
& = a^ k \mathrm{abs} \left| \begin{array}{llllll}    v_{1}^{1} &  \ldots &   v_{k}^{1} & v_{k+1}^{1}
& \ldots &  v_{n}^{1}\\
 v_{1}^{2} &  \ldots &  v_{k}^{2} & v_{k+1}^{2}
& \ldots &   v_{n}^{2}\\
\vdots  & \ldots & \vdots  &\vdots
& \ldots & \vdots \\
  v_{1}^{n-1}& \ldots &   v_{k}^{n-1} & v_{k+1}^{n-1}
& \ldots &   v_{n}^{n-1}\\
\frac{b}{a} v_{1}^{n}& \ldots & \frac{b}{a} v_{k}^{n} & v_{k+1}^{n}
& \ldots &  v_{n}^{n}
\end{array} \right| ,
\end{align*}
where abs stands for the absolute value and $v_i^j$ denotes the $j$-th coordinate of vector $v_i$.

Let $Q=\mathrm{diag}\Big\{  \underbrace{\frac{b}{a},\ldots,\frac{b}{a}}_{k}, \underbrace{1, \ldots,1}_{n-k} \Big\}$.
Denoting $v^i= (v_1^i, \ldots v_n^i)$, we see that
$$ D_n(A v_1, \ldots, A v_k,   v_{k+1}, \ldots  v_n) = a^k  D_n( v^1, \ldots,   v^{n-1},   Q  v^n).$$

Therefore, using an integral representation similar to \eqref{repr}, we get
\begin{align*} V_k(\ce)& =  \frac{C(n,k) }{2^k  \kappa_{n-1}^k}  \int \limits_{S^{n-1}} \ldots \int \limits_{S^{n-1}}  D_n(A v_1, \ldots, A v_k,   v_{k+1}, \ldots  v_n)\,   d v_1   \cdots dv_n\\
& = \frac{2^{n-k} C(n,k) }{ \kappa_{n-1}^k \left(\Gamma\left(\frac{n+1}2\right)\right)^n}  \int \limits_{\mathbb R^n} \ldots \int \limits_{\mathbb R^n}  D_n(A v_1, \ldots, A v_k,   v_{k+1}, \ldots  v_n)e^{-|v_1|^2 - \cdots - |v_n|^2}\,   d v_1   \cdots dv_n\\
& = \frac{2^{n-k}   C(n,k) a^k}{ \kappa_{n-1}^k \left(\Gamma\left(\frac{n+1}2\right)\right)^n}  \int \limits_{\mathbb R^n} \ldots \int \limits_{\mathbb R^n}   D_n( v^1, \ldots,   v^{n-1},   Q  v^n) e^{-|v^1|^2 - \cdots - |v^n|^2}\,   d v^1   \cdots dv^n\\
& = \frac{2^{-k} C(n,k)  a^k}{  \kappa_{n-1}^k }  \int \limits_{S^{n-1}} \ldots \int \limits_{S^{n-1}}  D_n( v^1, \ldots,  v^{n-1},   Q  v^n) \,   d v^1   \cdots dv^n\\
& = \displaystyle\frac{{n\choose k} \kappa_{n-1}}{n \, \kappa_{n-k} } a^k\,  V_1(\tilde \ce),
\end{align*}
where $\tilde \ce = QB_2^n$.

Thus, by \eqref{1stIV},
$$V_k(\ce) = \bar C(n,k)  a^k\, \int_{S^{n-1}} \left(\frac{b^2}{a^2} (\theta_1^2+\cdots + \theta_k^2)+ \theta_{k+1}^2+\cdots + \theta_n^2 \right)^{1/2} \, d\theta, $$
where $\bar C(n,k) =  \displaystyle \frac{{n\choose k}}{n\kappa_{n-k}}$.

Now suppose that we are given $V_n(\ce)$ and $V_k(\ce)$, for some $k$, $1\le k\le n-1$. It is easy to see that these two intrinsic volumes do not determine an ellipsoid of revolution uniquely. Indeed, after a rescaling, we can assume that $\ce$ has the same volume as the ball $B_2^n$, i.e., $a^{n-1} b=1$. Thus,
\begin{equation}\label{V_kREV}
V_k(\ce) = \bar C(n,k) \int_{S^{n-1}} \Big( a^{2k-2n} (\theta_1^2+\cdots + \theta_k^2)+ a^{2k}(\theta_{k+1}^2+\cdots + \theta_n^2) \Big)^{1/2} \, d\theta.
\end{equation}

The latter integral is a continuous function of $a\in (0,\infty)$, and it approaches infinity as $a\to 0$ and $a\to \infty$. This shows that there are two non-congruent ellipsoids of revolution with equal $V_n$ and $V_k$. Thus in order to prove a result in the positive direction, we need more intrinsic volumes.

\begin{theorem}
Let $\ce_1$ and $\ce_2$ be two ellipsoids of revolution in $\mathbb R^n$ such that $V_n(\ce_1) = V_n(\ce_2)$,  $V_k(\ce_1) = V_k(\ce_2)$, $V_{n-k}(\ce_1) = V_{n-k}(\ce_2)$ for some $k\ne n/2$, $1\le k\le n-1$. Then $\ce_1$ and $\ce_2$ are congruent.

If $n$ is even, then $V_n$, $V_{n/2}$, $V_k$, for any $k$ different from $n$ and $n/2$, uniquely determine an ellipsoid of revolution (up to an isometry).
\end{theorem}
\begin{proof}
	Without loss of generality we can assume that the volumes of $\ce_1$ and $\ce_2$ are equal to the volume of $B_2^n$. Thus we will deal with two ellipsoids whose semiaxes are $a,\ldots,a, a^{-n+1}$ and $b,\ldots,b, b^{-n+1}$ correspondingly. Our goal is to show that $a=b$.

	Consider the ellipsoids $\tilde \ce_1 = Q_a B_2^n$ and  $\tilde \ce_2 = Q_b B_2^n$, where
	$$ Q_a=\mathrm{diag}\,\{\underbrace{a^{-2k+2n},\ldots,a^{-2k+2n}}_{k}, \underbrace{a^{-2k}, \ldots,a^{2k}}_{n-k}\} $$ and $$ Q_b=\mathrm{diag}\,\{\underbrace{b^{-2k+2n},\ldots,b^{-2k+2n}}_{k}, \underbrace{b^{-2k}, \ldots,b^{-2k}}_{n-k}\}. $$
	By formula \eqref{V_kREV}, we see that
	 $V_k(\ce_1) = V_k(\ce_2)$ is equivalent to 	$\widetilde V_{-1}(\tilde \ce_1) = \widetilde V_{-1}(\tilde \ce_2)$. Replacing $k$ by $n-k$ in \eqref{V_kREV}, we see that $V_{n-k}(\ce_1) = V_{n-k}(\ce_2)$ is equivalent to $\widetilde V_{-1}(\tilde \ce_1^\circ) = \widetilde V_{-1}(\tilde \ce_2^\circ)$, which in turn is equivalent to $\widetilde V_{n+1}(\tilde \ce_1) = \widetilde V_{n+1}(\tilde \ce_2)$, by formula \eqref{relation_dual}. Thus Lemma \ref{MainFormula} gives
	 $$	  \int\limits_0^\infty  {u^{-2}} F(u)\, du =0  \quad\mbox{and} \quad  \int\limits_0^\infty  {u^{n}} F(u)\, du =0,$$
	 where
	 $$F(u)= \frac{1}{(1+a^{2k-2n}u^2)^{\frac{k}{2}}(1+a^{2k}u^2)^{\frac{n-k}{2}}} - \frac{1}{(1+b^{2k-2n}u^2)^{\frac{k}{2}}(1+b^{2k}u^2)^{\frac{n-k}{2}}}.$$
	 To use the ideas of Section \ref{SectionMR} it remains to show that $F$ has at most one positive root, unless it is identically equal to zero.
	 To this end observe that $F$ has the same roots as the function $g$ given by
	 $$g(u) = \ln \left((1+a^{2k-2n}u^2)^{\frac{k}{2}}(1+a^{2k}u^2)^{\frac{n-k}{2}}\right) - \ln \left( (1+b^{2k-2n}u^2)^{\frac{k}{2}}(1+b^{2k}u^2)^{\frac{n-k}{2}}\right).$$
Computing its derivative, we get
	 \begin{align*} g'(u)  =
	 u\left(\frac{ka^{2k-2n}-kb^{2k-2n}}{(1+a^{2k-2n}u^2)(1+b^{2k-2n}u^2)}+\frac{(n-k)a^{2k}-(n-k)b^{2k}}{(1+a^{2k}u^2)(1+b^{2k}u^2)}   \right).
	 \end{align*}
	 One can see that $g'$ has at most two positive real roots. Since $g(0)=0$ and $\lim_{u\to \infty}g(u) = 0$, we conclude that $g$ has at most one positive root.

	Now assume $n$ is even and $V_n(\ce_1) = V_n(\ce_2)=V_n(B_2^n)$, $V_{n/2}(\ce_1) = V_{n/2}(\ce_2)$, $V_k(\ce_1) = V_k(\ce_2)$ for some $k\ne n/2, n$. By Lemma \ref{MainFormula}, $V_{n/2}(\ce_1)$ is (up to a multiplicative constant) equal to
	$$\int\limits_0^\infty  {u^{-2}} \left(1 - \frac{1}{\big((1+a^{-n}u^2) (1+a^{n}u^2)\big)^{\frac{n}{4}}}\right)\, du. $$
	 The function $(1+a^{-n}u^2) (1+a^{n}u^2)$ considered as a function of $a$ is decreasing from $0$ to $1$ and increasing from $1$ to $\infty$ (for every $u$). Thus the integral above is decreasing/increasing on the same intervals with respect to $a$. Moreover, the integral is invariant under the transformation $a\to a^{-1}$.
	 Thus, if $V_{n/2}$ is prescribed a given value, then there are at most  two ellipsoids of revolution with this intrinsic volume and they must be polars of each other. Hence either $\ce_2=\ce_1$ or $\ce_2=\ce_1^\circ$. Assume the latter. By virtue of formula \eqref{KZ}, $V_k(\ce_1) = V_k(\ce_1^\circ)$ yields $V_{n-k}(\ce_1) = V_{n-k}(\ce_1^\circ)$. By what we proved above, this information guarantees that $\ce_1=\ce_1^\circ$. Thus $\ce_1=B_2^n$, and therefore $\ce_2=B_2^n$. So $\ce_1$ and $\ce_2$ cannot be different.
	\end{proof}