Equal Sum Sequences and Imbalance Sets of Tournaments

Abstract: Reid conjectured that any finite set of non-negative integers is the score set of some tournament and Yao gave a non-constructive proof of Reid's conjecture using arithmetic arguments. No constructive proof has been found since. In this paper, we investigate a related problem, namely, which sets of integers are imbalance sets of tournaments. We completely solve the tournament imbalance set problem (TIS) and also estimate the minimal order of a tournament realizing an imbalance set. Our proofs are constructive and provide a pseudo-polynomial time algorithm to realize any imbalance set. Along the way, we generalize the well-known equal sum subsets problem (ESS) to define the equal sum sequences problem (ESSeq) and show it to be NP-complete. We then prove that ESSeq reduces to TIS and so, due to the pseudo-polynomial time complexity, TIS is weakly NP-complete.

\section{Introduction}
A tournament is an orientation of a complete simple graph. In a tournament, the \emph{score} $s_{i}$ of a vertex $v_{i}$ is the number of arcs directed away from that vertex, that is, the outdegree of $v_{i}$. The \emph{score sequence} of a tournament is formed by listing the scores in nondecreasing order.
Let us write $[x_{i}]_{1}^{n}$ to denote a sequence with $n$ terms. Landau \cite{landau1} gave a simple characterization of the score sequences of tournaments.

\begin{theorem}\label{landau}
A sequence $[s_{i}]_{1}^{n}$ of non-negative integers in nondecreasing order is the score sequence of a tournament if and only if for every $I\subseteq\{1,2,\ldots,n\}$,
\begin{equation}\label{land}
\sum_{i\in I}{s_{i}}\geq{\left|I\right|\choose 2},
\end{equation}
\noindent with equality when $\left|I\right|=n$, where $\left|I\right|$ is the cardinality of the set $I$.
\end{theorem}

Several proofs of Landau's theorem have appeared over the years \cite{bang1, brualdi1, griggs1, landau1, thomassen1} and it continues to play a central role in the theory of tournaments and their generalizations. Brualdi and Shen \cite{brualdi2} strengthened Landau's theorem by deriving a set of inequalities that are individually stronger than inequalities (\ref{land}) but are collectively equivalent to these inequalities.

The set of scores of vertices in a tournament is called the \textit{score set} of the tournament. Reid \cite{reid1} conjectured that any finite nonempty set $S$ of non-negative integers is the score set of some tournament. He gave a constructive proof of the conjecture for the cases $\left|S\right|=1,2,3$, while Hager \cite{hager1} settled the cases $\left|S\right|=4,5$. In 1986, Yao announced a nonconstructive proof of Reid's theorem by arithmetic arguments \cite{yao1}. Pirzada and Naikoo \cite{pirzada1} obtained the construction of a tournament with a given score set in the special case when the score increments are increasing. However, so far no constructive proof has been found for Reid's theorem in general.

In a digraph, the \emph{imbalance} of a vertex $v_{i}$ is defined as $t_{i} = d_{i}^{+} -d_{i}^{-}$, where $d_{i}^{+}$ and $d_{i}^{-}$ are respectively the outdegree and indegree of $v_{i}$. The \emph{imbalace sequence} of a digraph is formed by listing the vertex imbalances in nonincreasing order. If $T$ is a tournament with imbalance sequence $[t_{i}]_{1}^{n}$, we say that $T$ \emph{realizes} $[t_{i}]_{1}^{n}$. Mubayi, Will and West \cite{west1} gave necessary and sufficient conditions for a sequence of integers to be the imbalance sequence of a simple digraph.

\begin{theorem}\label{west}
A sequence of integers $[t_{i}]_{1}^{n}$  with $t_{1}\geq\cdots\geq t_{n}$ is an imbalance sequence of a simple digraph if and only if $\sum_{i=1}^{j}t_i\leq j(n-j)$, for $1\leq j\leq n$ with equality when $j=n$.
\end{theorem}

On rearranging the imbalances in nondecreasing order, we obtain the equivalent inequalities $\sum_{i=1}^{j}t_i\geq j(j-n)$, for $1\leq j\leq n$ with equality when $j=n$.

Koh and Ree \cite{koh1} showed that if an additional parity condition is satisfied the sequence $[t_{i}]_{1}^{n}$ can be realized by a tournament. In fact, they proved the result in the more general setting of hypertournaments. The following corollary of Theorem 6 in \cite{koh1} provides a characterization of imbalance sequences of tournaments.

\begin{theorem}\label{seq}
A nonincreasing sequence $[t_{i}]_{1}^{n}$ of integers is the imbalance sequence of a tournament if and only if $n-1, t_{1}, \ldots, t_{n}$ have the same parity and
\begin{equation}\label{charac}
\sum_{i=1}^{j}{t_{i}} \leq j(n-j),
\end{equation}
\noindent for $j = 1, \ldots , n$ with equality when $j = n$.
\end{theorem}

In a digraph, the set of imbalances of the vertices is called its \textit{imbalance set} \cite{pirzada2}. In \cite{pirzada2} the following result regarding the imbalance sets of oriented graphs is proved.

\begin{theorem}\label{pir}
Let $P = \{p_1 , \ldots, p_m \}$ and $Q = \{-q_1 , \ldots, -q_n \}$, where $p_1 < \cdots < p_m$ and $q_1 < \cdots < q_n$ are positive integers. Then there exists an oriented graph with imbalance set $P \cup Q$.
\end{theorem}

Due to the interest in Reid's score set theorem, it is natural to ask if a similar result holds for imbalacne sets of tournaments. Furthermore, since a constructive proof of Reid's theorem has not yet been found, it would be interesting to look for an algorithm that generates a tournament from its imbalance set. In this paper we address both questions. We study the following decision problem and its search version.

\begin{definition}[\textbf{Tournament Imbalance Set Problem (TIS)}]\label{TIS}
Given a set $Z$ of integers, decide if $Z$ is the imbalance set of a tournament.
\end{definition}

In Section \ref{odd case}, we first show that the obvious necessary conditions for the existence of tournament imbalance sets are not sufficient. We then completely characterize the sets of odd integers that are imbalance sets of tounraments. In Section \ref{even case}, we treat the case of even integers, which is more involved. We show that any set of even integers that contains at least one positive and at least one negative integer or only consists of a single element 0, is the imbalnce set of a partial tournament in which each vertex is joined to every other vertex except one. However, not all such sets are imbalance sets of tournaments. This is followed by necessary and sufficient conditions for a set of even integers to be a tournament imbalance set. In Section \ref{algorithm}, we define a new variant of the \textit{equal sum subsets problem} (ESS) called \textit{equal sum sequences problem} (ESSeq). We show that ESSeq is NP-hard and ESSeq reduces to TIS in polynomial time. Furthermore, we propose a pseudo-polynomial time algorithm that determines if a set of integers is a tournament imbalance set and, in addition, generates a tournmanet realizing any such set. Thus TIS is shown to be weakly NP-complete. We also consider extremal cases and determine upper bounds for the minimal order of a tournament realizing an imbalance set.

\section{Characterizing odd imbalance sets}
\label{odd case}

Consider a tournament of order (number of vertices) $n$. Let $v_{i}$ be a vertex with score $s_{i}$ and imbalance $t_{i}$ then $t_{i}=d_{i}^{+}-d_{i}^{-}=s_{i}-(n-1-s_{i})=2s_{i}-(n-1)$, or by rearranging $s_{i} = \frac{n-1+t_{i}}{2}$. Converesly, assume that $v_{i}$ is a vertex of a tournament with $n$ vertices and $t_{i}$ is the imbalance of $v_{i}$. Let $s_{i} = \frac{n-1+t_{i}}{2}$ then $s_{i}$ is the score of $v_{i}$. Thus we have

\begin{lemma}\label{aux}
Let $t_{i}$ be the imbalance of a vertex $v_{i}$ in a tournament. Then $s_{i}$ is the score of $v_{i}$ if and only if
\begin{equation}\label{inter1}
s_{i} = \frac{n-1+t_{i}}{2},
\end{equation}
\noindent where $n$ is the order of the tournament.
\end{lemma}

A tournament is said to be \emph{regular} if all the vertices have the same score \cite{chartrand1}. Clearly, there exists a regular tournament on $n$ vertices with score $s_{i}$ if and only if $n$ is odd and $s_{i}=\frac{n-1}{2}$. Therefore, the imbalance of any vertex $v_{i}$ of a regular tournament is $t_{i}=2s_{i}-(n-1)=0$ and the imbalance set of any regular tournament is $\{0\}$.

The following is a set of obvious necessary conditions for an imbalance set of a tournament.

\begin{theorem}\label{necessary}
If a finite nonempty set $Z$ of integers is the imbalance set of a tournament of order $n$ then all the elements of $Z$ have the same parity as $n-1$ and it either contains at least one positive and at least one negative integer or contains only a single element 0.
\end{theorem}
\begin{proof}
If $Z$ is the imbalance set of a tournament with $n$ vertices then by Theorem \ref{seq}, the elements of $Z$ must have the same parity as $n-1$. Furthermore, either the tournament is regular and $Z=\{0\}$ or $Z$ must contain at least one positive and at least one negative integer so the corresponding imbalance sequence sums to zero.
\end{proof}

The question is whether these conditions are also sufficient. The answer is `no' as can be seen from the following example.

\begin{example}
Let $Z=\{6, -10\}$. Then $Z$ satisfies the necessary conditions given in Theorem \ref{necessary} and it can potentially be the imbalance set of a tournament with an odd number of vertices. However, any sequence with elements chosen from $Z$ can sum to zero only if it consists of an even number of elements (e.g., $6$, $6$, $6$, $6$, $6$, $-10$, $-10$, $-10$ ). Thus by Theorem \ref{seq}, we cannot construct a tournament imbalance sequence from $Z$ and so $Z$ is not a tournament imbalance set.
\end{example}

Although the conditions given in Theorem \ref{necessary} are not sufficient in general, they are sufficient if $Z$ consists of odd integers. We first show that

\begin{theorem}\label{odd}
Let $X=\{x_{1}, \ldots, x_{l}\}$ and $Y=\{-y_{1},\ldots,-y_{m}\}$ be disjoint nonempty sets of odd integers, where $x_{1}>\cdots> x_{l}$ are positive odd integers and $y_{1}<\cdots< y_{m}$ are also positive odd integers. Let $L=\sum_{i=1}^{l}{x_{i}}$, $M=\sum_{i=1}^{m}{y_{i}}$ and $n=lM+mL$. Then there exists a tournament of order $n$ with imbalance set $X\cup Y$.
\end{theorem}
\begin{proof} We observe that $n$ is even and all the elements of $X\cup Y$ have the same parity as $n-1$. Let $x^{(p)}$ denote that $x$ is appearing in $p$ consecutive terms of a sequence. We use Theorem \ref{seq} to prove that the $n$-term sequence
\[
[t_{i}]_{1}^{n}={x_{1}}^{(M)},\ldots,{x_{l}}^{(M)},{-y_{1}}^{(L)},\ldots,{-y_{m}}^{(L)}
\]
\noindent is the imbalance sequence of a tournament arranged in nonincreasing order. First note that
\[
\sum_{i=1}^{M}{t_{i}}= Mt_{1}=Mx_{1}\leq M((l-1)M+mL)=M(n-M),
\]
\[
\sum_{i=1}^{2M}{t_{i}}= M(t_{1}+t_{2})=M(x_{1}+x_{2})\leq 2M((l-2)M+mL)=2M(n-2M),
\]
\[\ldots \ \ \ \ \ldots \ \ \ \ \ldots \ \ \ \ \ldots \ \ \ \ \ldots
\]
\[
\sum_{i=1}^{lM}{t_{i}}= M\sum_{i=1}^{l}x_{i}=LM\leq lM(mL)=lM(n-lM),
\]
\[
\sum_{i=1}^{lM+L}{t_{i}}= LM-Ly_{1}\leq (lM+L)(m-1)L=(lM+L)(n-lM-L),
\]
\[
\sum_{i=1}^{lM+2L}{t_{i}}= LM-L(y_{1}+y_{2})\leq (lM+L)(m-2)L=(lM+2L)(n-lM-2L),
\]
\[
\ldots \ \ \ \ \ldots \ \ \ \ \ldots \ \ \ \ \ldots \ \ \ \ \ldots
\]
and
\[
\sum_{i=1}^{n}{t_{i}}=M\sum_{i=1}^{l}{x_{i}}-L\sum_{i=1}^{m}{y_{i}}=0=n(n-lM-mL).
\]

So inequality (\ref{charac}) holds for $j=M,2M,\ldots,lM,lM+L,lM+2L,\ldots,lM+mL(=n)$ with equality when $j=n$. Now suppose that for some other value $j=j_{0}$ we have $\sum_{i=1}^{j_{0}}{t_{i}} > j_{0}(n-j_{0})$ and $j_{0}$ is the smallest such integer. But then $t_{j_{0}}>n-2j_{0}+1$ and as $j_{0}\neq M,2M,\ldots,lM,lM+L,lM+2L,\ldots,n$, we have $t_{j_{0}+1}= t_{j_{0}}>n-2j_{0}+1>n-2j_{0}-1=n-2(j_{0}+1)+1$. Thus $\sum_{i=1}^{j_{0}+1}{t_{i}} > (j_{0}+1)(n-j_{0}-1)$, showing that $j_{0}+1\neq M,2M,\ldots,lM,lM+L,lM+2L,\ldots,n$. Continuing in this way leads to a contradiction as we must reach one of $M,2M,\ldots,lM,lM+L,lM+2L,\ldots,n$ in finitely many steps.
\end{proof}

Together Theorems \ref{necessary} and \ref{odd} immediately give the following necessary and sufficient conditions for odd tournament imbalance sets.

\begin{theorem}\label{oddsuff}
A finite nonempty set of odd integers is the imbalance set of a tournament if and only if it contains at least one positive and at least one negative integer.
\end{theorem}

\section{The case of even imbalances}
\label{even case}
Mubayi, Will and West \cite{west1} considered simple digraphs with maximum number of arcs that realize imbalance sequences.

\begin{lemma}\cite{west1}\label{aux2}
Let $D$ be a simple digraph with maximum number of arcs realizing the imbalance sequence $[t_{i}]_{1}^{n}$. Then any vertex in $D$ has at most one non-neighbour and the number of arcs in $D$ equals $\sum_{i=1}^{n}\left\lfloor \frac{n-1+t_{i}}{2}\right\rfloor$.
\end{lemma}

A \textit{partial tournament} is a simple digraph obtained by removing one or more arcs from a tournament \cite{brualdi1}. We say that a partial tournament of order $n$ is a \emph{near tournament} if each vertex is joined to all the other vertices except exactly one. Clearly, every near tournament has even order.

In this section, we characterize the sets of even integers that are imbalance sets of tournaments. Recall that $\{0\}$ is the imbalance set of every regular tournament. Therefore, in the remainder of this section we focus on nonzero sets of even integers. Example 1 shows that not every set of even integers that satisfies the necessary conditions of Theorem \ref{necessary} is the imbalance set of a tournament. We can nevertheless prove that any such set is the imbalance set of a near tournament.

\begin{theorem}\label{even}
Let $X=\{x_{1}, \ldots, x_{l}\}$ and $Y=\{-y_{1},\ldots,-y_{m}\}$ be disjoint nonempty sets of even integers, where $x_{1}>\cdots> x_{l}$ are non-negative even integers and $y_{1}<\cdots< y_{m}$ are positive even integers. Suppose that $X\cup Y\neq \{0\}$. Let $L=\sum_{i=1}^{l}{x_{i}}$, $M=\sum_{i=1}^{m}{y_{i}}$ and $n=lM+mL$. Then there exists a near tournament of order $n$ with imbalance set $X\cup Y$.
\end{theorem}
\begin{proof}
Since $L$ and $M$ are even, $n$ is even and so we cannot construct a tournament of order $n$ with imbalance set $X\cup Y$. By mirroring the proof of Case 1 of Theorem \ref{odd}, we can show that the $n$-term sequence
\[
[t_{i}]_{1}^{n}={x_{1}}^{(M)},\ldots,{x_{l}}^{(M)},{-y_{1}}^{(L)},\ldots,{-y_{m}}^{(L)}
\]
\noindent is the imbalance sequence of a simple digraph. Let $D$ be a realization of $[t_{i}]_{1}^{n}$ with maximum number of arcs. Since all the imbalances $t_{i}$ are even while $n-1$ is odd, by Lemma \ref{aux2} the number of arcs in $D$ is
\[
\sum_{i=1}^{n}\left\lfloor \frac{n-1+t_{i}}{2}\right\rfloor =\sum_{i=1}^{n}\frac{n-2+t_{i}}{2}=\frac{n(n-2)}{2},
\]

\noindent which is $\frac{n}{2}$ less than the number of arcs of a tournament of order $n$. Therefore, Lemma \ref{aux2} implies that every vertex of $D$ has exactly one non-neighbour and $D$ must be a near tournament.
\end{proof}

The following result shows that under certain conditions we can transform $D$ into a tournament by adding a suitable number of vertices.

\begin{theorem}\label{evensuff}
Let $X$, $Y$, $l$, $m$, $L$, $M$ and $n$ be as defined in Theorem \ref{even}.
\item (i) If $0\in X\cup Y$ then there exists a tournament of order $n+1$ with imbalance set $X\cup Y$.
\item (ii) If there exists an $x_{p}\in X$ and  (not necessarily distinct) $-y_{q},-y_{r}\in Y$ such that $x_{p}=y_{q}+y_{r}$ then there exists a tournament of order $n+3$ with imbalance set $X\cup Y$.
\item (iii) If there exists a $-y_{p}\in Y$ and (not necessarily distinct) $x_{q},x_{r}\in X$ such that $y_{p}=x_{q}+x_{r}$ then there exists a tournament of order $n+3$ with imbalance set $X\cup Y$.
\end{theorem}

\begin{proof} Let $T$ be a near tournament realizing the imbalance sequence $[t_{i}]_{1}^{n}$ as defined in the proof of Theorem \ref{even}. For a vertex $v_{i}$ in $T$ let $v'_{i}$ denote the unique non-neighbour of $v_{i}$. Also let $(u,v)$ denote an arc directed from vertex $u$ to vertex $v$. In the three cases we can transform $T$ into a tournament as follows.

\noindent \textit{(i)} Add a vertex $v$ to $T$ in such a way that for every pair of non-adjacent vertices $v_{i}$ and $v'_{i}$ we insert the arcs $(v_{i},v'_{i})$, $(v'_{i},v)$ and $(v,v_{i})$. Thus the imbalance of all the vertices of $T$ is preserved and the new vertex $v$ has imbalance 0. Since every vertex of $T$ has been linked with every other vertex of $T$ as well as the new vertex $v$, the resulting digraph is a tournament.

\begin{figure}[htb]
	\centering
		\includegraphics[scale=1.3]{3-3-2-ink.pdf}
	\caption{Construction of tournament in case (ii). Figure (a) represents step 1-2 and figure (b) step 3.}
	\label{tournamentex}
\end{figure}

\noindent \textit{(ii)} Add three new vertices $u_{1}$, $u_{2}$ and $u_{3}$ to $T$ and insert arcs in the following manner:
\begin{enumerate}
\item Insert $(u_{1},u_{2})$, $(u_{2},u_{3})$ and $(u_{3},u_{1})$.
\item Choose any $\frac{x_{p}}{2}=\frac{y_{q}+y_{r}}{2}$ pairs $\{v_{i},v'_{i}\}$ of non-neighbouring vertices in $T$ and insert $(v_{i},v'_{i})$. Out of these choose any $\frac{y_{q}}{2}$ pairs. For each of these pairs insert the arcs $(u_{1},v_{i})$, $(u_{1},v'_{i})$, $(v_{i},u_{2})$, $(v'_{i},u_{2})$, $(v'_{i},u_{3})$ and $(u_{3},v_{i})$. For the other $\frac{y_{r}}{2}$ pairs insert the arcs $(u_{1},v_{i})$, $(u_{1},v'_{i})$, $(v_{i},u_{3})$, $(v'_{i},u_{3})$, $(v'_{i},u_{2})$ and $(u_{2},v_{i})$.
\item For the remaining $\frac{n-x_{p}}{2}$ pairs $\{v_{i},v'_{i}\}$ of non-neighbours, insert the arcs  $(u_{1},v_{i})$, $(v'_{i},u_{1})$, $(v_{i},v'_{i})$, $(v_{i},u_{2})$, $(u_{2},v'_{i})$, $(u_{3},v_{i})$ and $(v'_{i},u_{3})$.
\end{enumerate}
Since every vertex is joined with every other vertex, the resulting digraph is a tournament. Furthermore, the imbalance of each vertex of $T$ is preserved, while the new vertices $u_{1}$, $u_{2}$ and $u_{3}$ have imbalances $x_{p}$, $-y_{q}$ and $-y_{r}$ respectively.

\noindent \textit{(iii)} The proof is essentially the same as that of case (ii).
\end{proof}

Theorem \ref{evensuff} can be generalized and it is the generalized version that is of interest to us as it leads to the characterization of even imbalance sets. However, we stated and proved Theorem \ref{evensuff} to provide the reader a concrete perspective of what is happening in the more abstract setting of Theorem \ref{evensuffgen}.

\begin{theorem}\label{evensuffgen}
Let $X$, $Y$, $l$, $m$, $L$, $M$ and $n$ be as defined in Theorem \ref{even}. The set $X\cup Y$ is the imbalance set of a tournament if any one of the following conditions is satisfied:

\item (i) $0\in X\cup Y$,

\item(ii) there exist an odd number of (not necessarily distinct) $x_{p_1},\ldots,x_{p_{2r+1}}\in X$ and an even number of (not necessarily distinct) $-y_{q_1},\ldots,-y_{q_{2s}}\in Y$ such that $\sum_{j=1}^{2r+1}x_{p_j}=\sum_{j=1}^{2s}y_{q_j}$,

\item(iii) there exist an odd number of (not necessarily distinct) $-y_{p_1},\ldots,-y_{p_{2r+1}}\in Y$ and an even number of (not necessarily distinct) $x_{q_1},\ldots,x_{q_{2s}}\in X$ such that $\sum_{j=1}^{2r+1}y_{p_j}=\sum_{j=1}^{2s}x_{q_j}$.
\end{theorem}

\begin{proof}
Let $T$ be a near tournament realizing the imbalance sequence $[t_{i}]_{1}^{n}$ as defined in the proof of Theorem \ref{even}.

\noindent \textit{(i)} The proof is exactly the same as part (i) of Theorem \ref{evensuff}.

\noindent \textit{(ii)} Add $2r+2s+1$ new vertices labelled $u_{p_1},\ldots,u_{p_{2r+1}}, u_{q_1}, \ldots,u_{q_{2s}}$ to $T$. Note that in the construction that follows we will relabel them in different ways, such as $u_{1},\ldots,u_{2r+2s+1}$, for the sake of convenience. We insert arcs in $T$ using the following procedure.

\begin{figure}[htb]
	\centering
		\includegraphics[scale=1.2]{3-4-1-ink.pdf}
	\caption{Figure (a) represents steps 1-4 and figure (b) step 5 of \textsc{Add Arcs}.}
	\label{tournamentex}
\end{figure}

\vspace{2mm}
\textsc{Add Arcs}:
\begin{enumerate}
\item Insert arcs so that the newly added vertices induce a regular tournament of order $2r+2s+1$.

\item Choose any $\frac{\sum_{j=1}^{2r+1}x_{p_j}}{2}=\frac{\sum_{j=1}^{2s}y_{q_j}}{2}$ pairs $\{v_{i},v'_{i}\}$ of non-neighbouring vertices in $T$ and order them arbitrarily. For each $i=1,\ldots,\frac{\sum_{j=1}^{2r+1}x_{p_j}}{2}$ insert the arc $(v_{i}, v'_{i})$. Therefore, the imbalances of $v_{i}$ and $v'_{i}$ change by $+1$ and $-1$ respectively, for all $i$.

\item For each $j=1,\ldots,2r+1$ choose $i=\frac{\sum_{h=1}^{p_{_{j-1}}}x_h}{2}+1,\ldots,\frac{\sum_{h=1}^{p_{j}}x_h}{2}$ and insert $x_{p_{j}}$ arcs $(u_{p_j},v_{i})$ and $(u_{p_j},v'_{i})$. This gives $u_{p_{j}}$ the imbalance $x_{p_{j}}$.

\item For each $j=1,\ldots,2s$ choose $i=\frac{\sum_{h=1}^{q_{_{j-1}}}y_h}{2}+1,\ldots,\frac{\sum_{h=1}^{q_{j}}y_h}{2}$ and insert the arcs $(v_{i},u_{q_j})$, $(v'_{i},u_{q_j})$. Thus the imbalance of $u_{q_{j}}$ is $-y_{q_{j}}$. But the imbalances of $v_{i}$ and $v'_{i}$ are still perturbed by $+1$ and $-1$ respectively, for $i=1,\ldots,\frac{\sum_{j=1}^{2r+1}x_{p_j}}{2}$.

\item For every $i=1,\ldots,\frac{\sum_{j=1}^{2r+1}x_{p_j}}{2}$ list the $u$'s that are not already linked with $v_{i}$ and $v'_{i}$. There are exactly $2r+2s-1$ such $u$'s, for each $i$. Label them arbitrarily from $1,\ldots,2r+2s-1$. For $j=1,\ldots,\frac{2r+2s-2}{2}$ insert the arcs $(u_{j},v_{i})$ and $(v'_{i},u_{j})$. For $j=\frac{2r+2s-2}{2}+1,\ldots,2r+2s-2$ insert the arcs $(u_{j},v_{i})$ and $(v'_{i},u_{j})$. Finally, insert the arcs $(u_{2r+2s-1},v_{i})$ and $(v'_{i}, u_{2r+2s-1})$. This preserves all the imbalances.

\item For the remaining $\frac{n-\sum_{j=1}^{2r+1}x_{p_j}}{2}$ pairs $\{v_{i}, v'_{i}\}$ of non-neighbours insert the arc $(v_{i},v'_{i})$. Label all the $u$'s arbitrarily from $1,\ldots,2r+2s+1$. For $j=1,\ldots,\frac{2r+2s}{2}$ insert the arcs $(u_{j},v_{i})$ and $(v'_{i},u_{j})$. For $j=\frac{2r+2s}{2}+1,\ldots,2r+2s$ insert the arcs $(u_{j},v_{i})$ and $(v'_{i},u_{j})$. Finally, insert the arcs $(u_{2r+2s+1},v_{i})$ and $(v'_{i}, u_{2r+2s+1})$. This preserves all the imbalances.

\begin{figure}[htb]
	\centering
		\includegraphics[scale=1.5]{3-4-2-ink.pdf}
	\caption{Inserting the remaining arcs in step 6 of \textsc{Add Arcs}.}
	\label{tournamentex}
\end{figure}

\end{enumerate}
\vspace{2mm}

Since every vertex is joined with every other vertex, the resulting digraph is a tournament. Furthermore, the imbalance of each vertex of $T$ is preserved, while the new vertices $u_{p_1},\ldots,u_{p_{2r+1}}, u_{q_1},\ldots,u_{q_{2s}}$ have imbalances $x_{p_1},\ldots,x_{p_{2r+1}}$, $-y_{q_1},\ldots,-y_{q_{2s}}$ respectively.

\noindent \textit{(iii)} The proof is essentially the same as that of case (ii).
\end{proof}

The reader can easily draw parallels between the proofs of Theorems \ref{evensuff}(ii) and \ref{evensuffgen}(ii). For instance, the first and the last steps of both proofs are essentially achieveing the same target while steps 2-5 of the later are similar to but more complicated than step 2 of the former.

Analyzing the above proof leads to a couple of simpler sufficient conditions for tournament imbalance sets. The first is one is a fairly straightforward consequence of Theorem \ref{evensuff} (i).
\begin{corr}\label{zero}
If $Z$ is the empty set or it contains at least one positive and at least one negative even integer then $Z\cup \{0\}$ is the imbalance set of a tournament.
\end{corr}

The second condition is not as obvious and is more of an arithmetic result than a combinatorial one. First, we note that for any positive integer $p\geq 1$, the set $\{2^{p}, -2^{p}\}$ is not a tournament imbalance set as any zero sum sequence formed by the elements of this set necessarily consists of an even number of elements. However, the following sufficient condition shows that any other set of positive and negative even integers containing a power of $2$ is a tournament imbalance set.

\begin{corr}\label{arithmetic}
Let $Z$ be a finite nonempty set of even integers containing at least one positive and at least one negative integer. Suppose $Z$ contains an element of the form $2^{p}$ or $-2^{p}$, for some positive integer $p\geq 1$, and $Z\neq \{2^{p}, -2^{p}\}$. Then $Z$ is the imbalance set of a tournament.
\end{corr}
\begin{proof} Let us assume that for some positive integer $p\geq 1$, $2^{p}$ is an element of $Z$. Choose any negative element $-y\in Z$. Then $y=r2^{q}$, where $q\geq 1$ is a positive integer and $r\geq 1$ is an odd positive integer such that if $r=1$ then $q\neq p$. (If this is not possible, we can start with $-2^{p}\in Z$ and choose an $x=r2^{q}$ from $Z$ with $q\neq p$.) Without loss of generality, let $p=\max\{p,q\}$. We have
\[\underbrace{2^{p}+\cdots+2^{p}}_{r \textnormal{ terms}}=\underbrace{y+\cdots+y}_{2^{p-q} \textnormal{ terms}},
\]
\noindent and by Theorem \ref{evensuffgen}, $Z$ is a tournament imbalance set.
\end{proof}

After deriving a number of sufficient conditions for tournament imbalance sets of even integers, the natural question is whether the sufficient conditions given in Theorem \ref{evensuffgen} are also necessary. The answer is positive as seen from the following result.

\begin{theorem}\label{evennecess}
Let $Z=X\cup Y$ be a finite nonempty set of even integers, where $X$ is the set of non-negative integers and $Y$ is the set of negative integers in $Z$. Then $Z$ is the imbalance set of a tournament if and only if either $Z=\{0\}$ or both $X$ and $Y$ are nonempty and satisfy one of the conditions (i), (ii) or (iii) of Theorem \ref{evensuffgen}.
\end{theorem}
\begin{proof}
The sufficiency follows from Theorem \ref{evensuffgen}.

To prove the necessity, suppose that $0\notin X\cup Y$ and let $X\cup Y$ be the imbalance set of a tournament of order $k$. This implies that we can form a sequence $[t_{i}]_{1}^{k}$ consisting of an odd number of not necessarily distinct terms from the elements of $X\cup Y$ that sums to zero. Since $k$ is odd, therefore either the number of terms from $X$ is odd or the number of terms from $Y$ is odd, but not both. Thus we have an odd (respectively even) number of terms $x\in X$ and an even (respectively odd) number of terms $-y\in Y$ such that $\sum{x}=\sum{y}$.
\end{proof}

\section{Algorithmic aspects}
\label{algorithm}
The aim of this section is to study the \textit{tournament imbalance set problem} (TIS) and present an algorithm that generates a tournament realizing any tournament imbalance set. We begin by proving a theorem on the lengths of equal sum sequences chosen from two set of non-negative integers that will play a crucial role in developing the algorithm.

\begin{theorem}\label{zerosumsequence}
Let $X$, $Y$, $l$, $m$, $L$, $M$ and $n$ be as defined in Theorem \ref{even}. If $k=p+q$ is the minimum odd number such that there exists a $p$-term sequence from $X$ and a $q$-term sequence from $-Y=\{y:-y\in Y\}$ having the same sum, then $k<n$.
\end{theorem}
\begin{proof}
We observe that $k$ equals the minimal length of a zero sum sequence from $X\cup Y$. We prove the result by induction on $n$. Note that according to the conditions of Theorem \ref{even}, $X\cup Y\neq \{0\}$ and so the minimum possible value of $n$ is 4 that only corresponds to the sets $X=\{2\}$ and $Y=\{-2\}$. Since it is not possible to form a zero sum sequence with odd number of terms from $\{2, -2\}$. The next smallest value of $n$ is $6$ that corresponds to the sets $\{2,0,-2\}$, $\{2,-4\}$ and $\{4,-2\}$. Each of these sets admits a zero sum sequence of length $k=3$ and so $k<n$.

Now we aim to show that the result holds for any $n>6$ by assuming that it holds for all values less than $n$. Let $X$ and $Y$ be any two sets of integers corresponding to $n$ and let $k\geq 3$ be the minimum odd number such that there exists a $k$-term zero sum sequence $a_{1},\ldots, a_{p},-b_{1},\ldots,-b_{q}$, where $k=p+q$, $a_{1}>\cdots>a_{p}\in X$ and $-b_{1}>\cdots> -b_{q}\in Y$. Assume that $k>n$, then $k\geq 5$. The sequence $a_{1}+a_{2},\ldots,a_{p},-b_{1}\ldots,-b_{q-1}-b_{q}$  is a zero sum sequence of minimal odd length $k'=k-2$ corresponding to the sets $X'=X-\{a_{1},a_{2}\}\cup \{a_{1}+a_{2}\}$ and $Y'=Y-\{-b_{q-1},-b_{q}\}\cup\{-b_{q-1}-b_{q}\}$. For the sets $X'$ and $Y'$ we have $n'<n-2$ and so $k'>n'$, contradicting the induction hypothesis.
\end{proof}

Given a set of non-negative integers, we call the search problem of finding two disjoint nonempty subsets that have identical sums the \textit{equal sum subsets problem} (ESS). Several authors \cite{bazgan1,woeginger1} have considered the corresponding decision and optimization problems. Additionally, many variants of ESS have been studied in literature. For instance, if we require subsets to be found from two different sets of positive integers the problem is called \textit{equal sum subsets from two sets} (ESST). It is known that ESS and ESST are weakly NP-hard as they admit pseudo-polynomial time algorithms \cite{cieliebak1,woeginger1}. The best known algorithm for the ESS is the dynamic programming procedure by Bazgan, Santha and Tuza \cite{bazgan1} that runs in $O(\left|I\right|\times{Sum}^{2})$ time, and determines all possible solutions of an ESS instance. Here $\left|I\right|$ and $Sum$ respectively denote the number of elements and the sum of elements of the input set. This procedure can be easily adapted to solve ESST \cite{cieliebak1}. Here we are interested in the following variation of ESS.

\begin{definition}[\textbf{Equal Sum Sequences Problem (ESSeq)}]\label{ESSeq}
Given two sets $X$ and $Y$ of non-negative integers and a positive integer $k$, find two nonempty finite sequences $[x]$ and $[y]$ consisting of elements from $X$ and $Y$ respectively, with each element allowed to repeat at the most $k$ times, such that $\sum x = \sum y$.
\end{definition}

We now study the complexity of ESSeq. Clearly, ESSeq is in class NP. Furthermore, ESST corresponds to the special case $k=1$ of ESSeq. Thus ESSeq is NP-hard. In fact, ESSeq is weakly NP-hard as we can solve any instance $ESSeq(X,Y,k)$ by using the multisets $X^{(k)}$ and $Y^{(k)}$, in which each element is repeated $k$ times, as input for the pseudo-polynomial ESST algorithm. Let us call the resulting algorithm, which finds all possible solutions to an ESSeq instance, \textsc{Equal Seq}. We have shown the following.

\begin{theorem}\label{esseq np}
The ESSeq decision (search) problem is weakly NP-complete (weakly NP-hard).
\end{theorem}

On the other hand, any algorithm that solves the even case of TIS must be able to check the existence of equal sum sequences $[x]_{1}^{a}$ and $[y]_{1}^{b}$ from any given nonempty finite sets $X$ and $Y$ of even integers such that $a$ and $b$ have different parity. Let $hX$ denote the set obtained by multiplying every element of a set $X$ by $h$ and $X+\{h\}$ denote the set obtained by adding a number $h$ to every element of $X$. We can solve an instance $ESSeq(X,Y,k)$ of the ESSeq decision problem by using any TIS algorithm to solve $\left|X\right|+1$ even instances $TIS(2X,2Y)$, $TIS(2(X+\{x\}),2Y)_{x\in X}$ of TIS.

\begin{theorem}\label{np}
The odd case of TIS can be solved in linear time. On the other hand, the even case of TIS is NP-complete. Hence in general TIS is NP-complete.
\end{theorem}

We now present a pseudo-polynomial time algorithm that not only solves TIS but also generates a tournaments realizing any tournament imbalance set. Our algorithm is based on the proofs of Theorems \ref{odd}, \ref{evensuffgen} and \ref{zerosumsequence}. First we form a suitable $n$-term imbalance sequence and then realize it as a tournament. Lemma \ref{aux2} can be used to construct a simple digraph, with maximum possible number of arcs, realizing an imbalance sequence $[t_{i}]_{1}^{n}$. The idea is to start with an arbitrary vertex $v$ having imbalance $t_{i}$ and attach it to $\left\lfloor \frac{n-1+t_{i}}{2}\right\rfloor$ vertices by arcs directed away from $v$. If $t_{i}$ has the same parity as $n-1$ then it is joined with $n-1-\left\lfloor \frac{n-1+t_{i}}{2}\right\rfloor$ other vertices by arcs directed towards $v$. Otherwise, it is joined with $n-2-\left\lfloor \frac{n-1+t_{i}}{2}\right\rfloor$ other vertices by arcs directed towards $v$. Thus $v$ is joined to every vertex except possibly one. These steps are then repeated for every vertex without attaching any new arcs to the preprocessed vertices. We name this $O(n^{2})$ procedure \textsc{Max Realization}.

Now suppose that $Z$ is a finite nonempty set of integers arranged in decreasing order. Form the sets $X=\{z\in Z:z\geq 0\}=\{x_{1}, \ldots, x_{l}\}$ and $Y=\{z\in Z:z<0\}=\{-y_{1}, \ldots, -y_{m}\}$ arranged in decreasing order. Let $L=\sum_{i=1}^{l}{x_{i}}$, $M=\sum_{i=1}^{m}{y_{i}}$ and $n=lM+mL$ as in the earlier proofs. The following algorithm outputs a tournament that realizes $Z$, whenever such a tournament exists.

\begin{algo}[\textsc{Imbalance Set}]\label{generate}
\begin{enumerate}

\item If either $X$ or $Y$ is empty, then $Z$ is not a tournament imbalance set. Stop.

\item If elements of $Z$ have different parity, $Z$ is not a tournament imbalance set. Stop.

\item Form the sequence $[t_{i}]_{1}^{n}={x_{1}}^{(M)},\ldots,{x_{l}}^{(M)},{-y_{1}}^{(L)},\ldots,{-y_{m}}^{(L)}.$

\item Call the procedure \textsc{Max Realization} to realize $[t_{i}]_{1}^{n}$ as a simple digraph $D$ with maximum number of arcs.

\item If elements of $Z$ have odd parity, output $D$. End.

\item If elements of $Z$ have even parity, call \textsc{Equal Seq} with the input $(X^{(n)}, (-Y)^{(n)}, n)$ to find sequences $[x]_{1}^{a}$ and $[y]_{1}^{b}$, with $a$ and $b$ having different parity and $\sum x=\sum y$. If no such sequences exist then $Z$ is not a tournament imbalance set. End.

\item Add $2a+2b+1$ isolated vertices to $D$.

\item Call \textsc{Add Arcs} to add $a+b$ vertices and arcs to $D$ to form a tournament $T$. Return $T$.

\end{enumerate}
\end{algo}

The following result shows that the procedure \textsc{Imbalance Set} runs in pseudo-polynomial time and hence TIS is weakly NP-hard.
\begin{theorem}\label{correct}
Algorithm \ref{generate} is correct and runs in pseudo-polynomial time.
\end{theorem}
\begin{proof}
The correctness follows immediately from Theorems \ref{odd}, \ref{even}, \ref{evensuffgen} and \ref{zerosumsequence}. In particular, Theorem \ref{zerosumsequence} guarantees that in the case when $Z$ is a set of even integers, step 6 of Algorithm \ref{generate} necessarily finds the required sequences if they exist. Now note that the computational complexity of Algorithm \ref{generate} is dominated by steps 4 and 6. Step 4 can be performed in $O(n^{2})=O((lM+mL)^{2})$ time, whereas step 6 takes $O(( n\left| X \right| + n\left| Y \right|)\times(n\sum_{x\in X}x + n\sum_{-y\in Y}y)^{2})=O(n^{3}(l+m)\times(L+M)^{2})=O(n^{3}\left|Z\right|\times (L+M)^{2})$ time. The overall complexity is therefore $O(\left|Z\right|\times n^{5})$, which is pseudo-polynomial since $n$ depends on the numeric value of the input.
\end{proof}

Thus we can use Algorithm \ref{generate} to check if a given set of integers is the imbalance set of a tournament and moreover, to construct a tournament realizing the set, if it exists. We now illustrate Algorithm \ref{generate} by showing how it generates a tournament realizing the imbalance set $\{4,2,-2\}$.

\begin{figure}[htb]
\centering
\includegraphics[scale=0.7]{tournament1.pdf}
\caption{A tournament realizing the imbalance set $\{4, 2, -2\}$ obtained from  Algorithm 1.}
\label{tournamentex}
\end{figure}

\begin{example}
Consider the set $Z=\{4, 2, -2\}$. Since $Z$ satisfies the conditions in the first two steps of Algorithm 1, the algorithm goes to step 3 and forms the sequence $4, 4, 2, 2, -2, -2, -2, -2, -2, -2$. Step 4 calls the procedure \textsc{Max Realization} to output a simple digraph of order 10 realizing $Z$. However, this simple digraph is only a near tournament and not a tournament (see the digraph induced by the black vertices in Figure \ref{tournamentex}). Since the elements of $Z$ have even parity, the algorithm proceeds to step 6 and finds the the sequences $[x_{i}]_{1}^{1}=4$ and $[y_{i}]_{1}^{2}=-2,-2$ with odd and even number of terms respectively, such that $\sum{x_{i}} = -\sum{y_{i}}$. Then step 7 adds 3 new vertices (colored white in Figure \ref{tournamentex}) to the near tournament obtained in step 4. In the end, step 8 adds the arcs (dashed arcs in Figure \ref{tournamentex}) necessary to form a tournament of order 13 in such a way that the imbalances of the old vertices are preserved and the new vertices have imbalances $4$, $-2$, and $-2$. The output of Algorithm 1 is a tournament with imbalance sequence $4, 4, 4, 2, 2, -2, -2, -2, -2, -2, -2, -2, -2$ as shown in Figure \ref{tournamentex}.
\end{example}

The results presented in Sections 2, 3 and 4 can be used to estimate the order of a tournament realizing an imbalance set. Let us denote by $ord(Z)$ the minimal order of a tournament realizing an imbalance set $Z$.

\begin{theorem}\label{extremal}
Let $Z$ be a tournament imbalance set (i.e., it satisfies the conditions of Theorem 2.4 or Theorem 3.6). Define $X$, $Y$, $l$, $m$, $L$, $M$ and $n$ as in the earlier results.
\item(i) If $Z$ consists of odd integers then $ord(Z)\leq n=lM+mL$.
\item(ii) If $Z$ consists of even integers and $0\in Z$ then $ord(Z)\leq n+1=lM+mL+1$.
\item(iii) If $Z$ consists of even integers and $0\notin Z$ then $ord(Z)< 2n=2(lM+mL)$.
\end{theorem}
\begin{proof}
The proof of (i) and (ii) follows from Theorem \ref{odd}, while (iii) follows from Theorem \ref{evensuff}(i). For (iv), observe that the order of the tournament constructed in the proof of Theorem \ref{evensuffgen} is $n+2r+2s+1$. From Theorem \ref{zerosumsequence}, $2r+2s+1< n$. As a result, the constructed tournament can be at the most of order $2n-1$.
\end{proof}

\bigskip
\noindent \textbf{Acknowledgements}

\noindent The author would like to thank Professor K\'{a}roly Bezdek for having many useful discussions on the topic and helping in impproving this manuscript. The author is also thankful to Dr. Shariefuddin Pirzada for drawing his attention to the imbalance set problem. This research was performed at the Center for Computational and Discrete Geometry, Department of Mathematics and Statistics, University of Calgary.


Duality Mapping for Schatten Matrix Norms

Abstract: In this paper, we fully characterize the duality mapping over the space of matrices that are equipped with Schatten norms. Our approach is based on the analysis of the saturation of the Hölder inequality for Schatten norms. We prove in our main result that, for $p\in (1,\infty)$, the duality mapping over the space of real-valued matrices with Schatten-$p$ norm is a continuous and single-valued function and provide an explicit form for its computation. For the special case $p = 1$, the mapping is set-valued; by adding a rank constraint, we show that it can be reduced to a Borel-measurable single-valued function for which we also provide a closed-form expression.

\section{Introduction}
\label{Sec:intro}
In linear algebra and matrix analysis, Schatten norms are a family of spectral  matrix norms that are defined via the singular-value decomposition \cite{bhatia2013matrix}. They have  appeared in many applications    such as image reconstruction \cite{lefki2013Poisson,lefki2013HS}, image denoising \cite{xie2016weighted}, and tensor decomposition \cite{gao2020robust}, to name a few.

Generally, the Schatten-$p$ norm of a matrix  is the $\ell_p$ norm of  its singular values. The family  contains some well-known matrix norms: The Frobenius and the spectral (operator) norms are special cases in the family, with $p=2$ and $p=\infty$, respectively.  The case $p=1$ (trace or nuclear norm) is of particular interest for applications as it can be used to recover low-rank matrices \cite{davenport2016overview}. This is the current paradigm in matrix completion, where  the goal is to recover an unknown matrix given some of its entries \cite{candes2009exact}.   Prominent examples of applications that can be reduced to low-rank matrix-recovery problems are phase retrieval \cite{candes2015phase}, sensor-array processing \cite{davies2012rank},  system identification \cite{fazel2013hankel}, and index coding \cite{asadi2017fast,esfahanizadeh2014matrix}.

In addition to their many applications in data science, Schatten norms have been extensively studied from a theoretical point of view. Various inequalities concerning Schatten norms have been proven \cite{kittaneh1985inequalities,kittaneh1987inequalities2,kittaneh1986inequalities3,kittaneh1986inequalities4,kittaneh1987inequalities5,bourin2006matrix,Hirzallah2010Schatten,moslehian2011schatten,conde2016norm};
 sharp bounds for commutators  in Schatten spaces have been given \cite{wenzel2010impressions,cheng2015schatten}; moreover,  facial structure \cite{so1990facial}, Fr\'echet differentiablity \cite{potapov2014frechet}, and various other aspects \cite{kittaneh1989continuity,bhatia2000cartesian}  have been studied already.

 Our objective in this paper is to investigate the duality mapping in spaces of matrices that are equipped with Schatten norms. The duality mapping is a powerful tool to understand the topological structure of Banach spaces \cite{beurling1962theorem,cioranescu2012geometry}. It has been used to derive powerful characterizations of the solution of variational problems in function spaces \cite{de1976best,unser2019unifying} and also to determine generalized linear inverse operators \cite{liu2007best}.  Here, we prove that the duality mapping over Schatten-$p$ spaces with $p\in (1,+\infty)$ is a single-valued and continuous function which, in fact, highlights the strict convexity of these spaces. For the special case $p=1$, the mapping is set-valued. However, we prove that, by adding a rank constraint, it reduces to a single-valued Borel-measurable function. In both cases, we also derive closed-form expressions that allow one to compute them explicitly.

The paper is organized as follows: In Section \ref{Sec:Prelim}, we present relevant mathematical tools and concepts that are used in this paper. We study the duality mapping of Schatten spaces and propose our main result in Section \ref{Sec:DualMap}. We provide further discussions regarding the introduced mappings in Section \ref{Sec:discuss}.
\section{Preliminaries}
\label{Sec:Prelim}
\subsection{Dual Norms, H\"older Inequality, and Duality Mapping}
\label{Sec:Holder}
Let $V$ be a finite-dimensional vector space that is equipped with an inner-product $\langle \cdot,\cdot \rangle:V \times V \rightarrow \mathbb{R}$ and let $\|\cdot\|_{X}:V\rightarrow \mathbb{R}_{\geq 0}$ be an arbitrary norm on $V$. We then denote by $X$  the space $V$ equipped with $\|\cdot\|_X$. Clearly, $X$ is a Banach space, because all finite-dimensional normed spaces are complete. The dual norm of $X$, denoted by $\|\cdot\|_{X'}: V\rightarrow \mathbb{R}_{\geq 0}$, is defined as
\begin{equation}\label{Eq:DualNorm}
\|{\bf v}\|_{X'} = \sup_{{\bf u} \in V\backslash \{\boldsymbol{0}\} } \frac{\langle {\bf v} , {\bf u} \rangle}{\|{\bf u}\|_{X}},
\end{equation}
for any ${\bf v} \in V$. Following this definition, one would directly obtain the generic duality bound
\begin{equation}\label{Eq:DualityBound}
\langle {\bf v}, {\bf u} \rangle  \leq \|{\bf v}\|_{X'} \|{\bf u}\|_{X},
\end{equation}
for any ${\bf v}, {\bf u} \in V$. Saturation of Inequality \eqref{Eq:DualityBound} is the key concept of dual conjugates that is formulated in the following definition.
\begin{definition}\label{Def:DualConj}
Let $V$ be a finite-dimensional vector space and let $(\|\cdot\|_X,\|\cdot\|_{X'})$ be a pair of dual norms that are defined over $V$. The pair $({\bf u},{\bf v}) \in V\times V$ is said to be a  $(X,X')$-conjugate, if
\begin{itemize}
\item $\langle {\bf v}, {\bf u}\rangle = \|{\bf v}\|_{X'} \|{\bf u}\|_{X}$,
\item $\|{\bf v}\|_{X'} = \|{\bf u}\|_{X}$.
\end{itemize}
For any ${\bf u}\in V$, the set of all elements $ {\bf v} \in V$ such that $({\bf u},{\bf v})$ forms an  $(X,X')$-conjugate is denoted by $\mathcal{J}_{X}({\bf u})\subseteq V$. We refer to the set-valued mapping $\mathcal{J}_{X}: V \rightarrow 2^V$ as the duality mapping. If, for all ${\bf u}\in V$, the set  $\mathcal{J}_{X}({\bf u})$ is a singleton, then we indicate the duality mapping for the $X$-norm via the single-valued  function ${\rm J}_X:V\rightarrow V$ with $\mathcal{J}_{X}({\bf u}) = \{{\rm J}_X({\bf u})\}$.
\end{definition}
It is worth mentioning that,  for any ${\bf u}\in V$, the set  $\mathcal{J}_{X}({\bf u})$  is nonempty .  In
fact, the closed ball $B= \{ {\bf v}\in V : \|{\bf v} \|_{X'}= \|{\bf u} \|_{X} \}$ is   compact   and, hence,  the function ${\bf v}\mapsto  \langle {\bf v} , {\bf u}\rangle$ attains its maximum value at some  ${\bf v}^*\in B$.  Now, following Definition \ref{Def:DualConj}, one readily verifies that  $({\bf u},{\bf v}^*)$ is an  $({X},{X}')$-conjugate.

We conclude this part by providing a classical and illustrative example. Let $V= \mathbb{R}^n$ for some $n\in\mathbb{N}$. For any $p \in [1,+\infty]$, the $\ell_p$-norm of a vector ${\bf u} =(u_i)\in \mathbb{R}^n$ is defined as
\begin{equation}\label{Eq:lp}
\|{\bf u}\|_{p} = \begin{cases} \left(\sum_{i=1}^n |u_i|^p\right)^{\frac{1}{p}}, & p< +\infty \\ \max_{i} |u_i|, & p= +\infty. \end{cases}
\end{equation}
It is widely known that the dual norm of $\ell_p$ is the $\ell_q$-norm, where $(p,q)$ are H\"older conjugates ({\it i.e.},  $1/p+1/q=1$) \cite{rudin1991functional}. This stems from the H\"older inequality which states that
\begin{equation}\label{Eq:HolderVec}
\langle {\bf v}, {\bf u} \rangle \leq \| {\bf u}\|_{p} \|{\bf v}\|_q,
\end{equation}
for all ${\bf u}=(u_i),{\bf v}=(v_i) \in \mathbb{R}^n$. In the sequel, we exclude the trivial cases   ${\bf u}=\boldsymbol{0}$ and ${\bf v}=\boldsymbol{0}$ to avoid unnecessary complexities in our statements.

When $1<p<+\infty$, Inequality \eqref{Eq:HolderVec} is saturated if and only if $u_i v_i \geq 0$ for $i=1,\ldots,n$ and  there exists a constant $c >0$ such that $|{\bf u}|^p = c |{\bf v}|^q$, where $|{\bf u}|^p = (|u_i|^p)$.  This ensures that the duality mapping is single-valued  and also yields the map
\begin{equation}\label{Eq:DualMapLp}
\mathrm{J}_p({\bf u})= {\rm sign}({\bf u}) \frac{|{\bf u}|^{p-1}}{\|{\bf u}\|_p^{p-2}}.
\end{equation}
 For $p=1$, one can verify that the equality happens if and only if,    for any index $i=1,\ldots,n$ with $u_i \neq 0$, one has that
\begin{equation}\label{Eq:HolderEq1}
v_i =   {\rm sign} (u_i)\|{\bf v}\|_{\infty}.
\end{equation}
In other words, the vector ${\bf v}$ should attain its extreme values at places where ${\bf u}$ has nonzero values, with the sign being determined by the corresponding element in ${\bf u}$.

Due to \eqref{Eq:HolderEq1}, the set $\mathcal{J}_1({\bf u})$ is not necessarily a singleton. However,  if we add an additional sparsity constraint, then the mapping becomes single-valued. This leads us to introduce the new notion of {\it sparse duality mapping} in Definition \ref{Def:SparseDual}.
\begin{definition}\label{Def:SparseDual}
Let $V$ be a finite-dimensional vector space and let $s_0: V \rightarrow \mathbb{N}$ be an integer-valued function that acts as a sparsity measure.  Assuming a pair $(\|\cdot\|_X,\|\cdot\|_{X'})$ of dual norms over $V$, we call the pair $({\bf u},{\bf v})\in V\times V$  a sparse $(X,X')$-conjugate  if
\begin{itemize}
\item $({\bf u},{\bf v})$ forms an $(X,X')$-conjugate pair. In other words, ${\bf v} \in \mathcal{J}({\bf u})$.
\item The quantity $s_0 ( {\bf v})$ attains its minimal value over the  set   $\mathcal{J}({\bf u})$.
\end{itemize}
We denote the set of sparse conjugates of ${\bf u}$ by $\mathcal{J}_{X,s_0}({\bf u})$. Whenever $\mathcal{J}_{X,s_0}({\bf u})$ is a singleton for any ${\bf u}\in V$, we refer to the single-valued function ${\rm J}_{X,s_0}: V \rightarrow V$ with $\mathcal{J}_{X,s_0}({\bf u})=\{{\rm J}_{X,s_0}({\bf u})\}$ as the sparse duality mapping.
\end{definition}
Following Definition \ref{Def:SparseDual}, if we use the $\ell_0$-norm as the sparsity measure, that is $s_0 ({\bf u}) = \|{\bf u}\|_0={\rm Card}\left(\{i: u_i\neq 0\}\right)$\footnote{Although this functional  does not satisfy the homogeneity property of a norm, it has been widely referred to as the $\ell_0$-norm.}, then we have the sparse duality mapping
\begin{align}\label{Eq:DualMapSparseL1}
&{\rm J}_{1,0}:\mathbb{R}^n\rightarrow\mathbb{R}^n:{\bf u}=(u_i) \mapsto {\bf v}=(v_i) ={\rm J}_{1,0}({\bf u}), \nonumber \\
&v_i = \begin{cases} {\rm sign}(u_i) \|{\bf u}\|_1, & u_i\neq 0 \\
0, & u_i=0.\end{cases}
\end{align}
Finally, we mention that, for $p= +\infty$,  the  reduced set $\mathcal{J}_{\infty,0}$ is  not single-valued. Indeed, let us define $I_{\max}({\bf u})=\{i: |u_i|=\|{\bf u}\|_{\infty}\}\subseteq \{1,\ldots,n\}$. We readily deduce from \eqref{Eq:HolderEq1} that ${\bf v}=(v_1,\ldots,v_n) \in  \mathcal{J}_{\infty}({\bf u})$ if and only if $v_i = 0$ whenever $i\not \in I_{\max}({\bf u})$ and ${\rm
sign}(v_i) = {\rm sign}(u_i)$ for $i \in I_{\max}({\bf u})$ with $\sum_{i\in I_{\max}({\bf u})} |
v_i|= \|{\bf u}\|_{\infty}$. This shows that $\mathcal{J}_{\infty}({\bf u})$ is a convex set with $\mathcal{J}_{\infty,0}({\bf u})$ being its extreme points, where  $\mathcal{J}_{\infty,0}({\bf u})=\{u_i {\bf e}_i: i\in  I_{\max}({\bf u})\}$.
\subsection{Schatten $p$-Norm}
It is widely known that any  matrix ${\bf A} \in \mathbb{R}^{m\times n}$  can be decomposed as
\begin{equation}\label{Eq:SVD}
{\bf A}= {\bf U} {\bf S} {\bf V}^T,
\end{equation}
where ${\bf U}\in\mathbb{R}^{m\times m}$  and ${\bf V}\in\mathbb{R}^{n\times n}$ are orthogonal matrices and ${\bf S}$ is an $m$ by $n$ rectangular diagonal matrix with nonnegative real entries $\sigma_1 \geq \sigma_2 \geq \cdots \geq \sigma_{\min(m,n)} \geq 0$ sorted in descending order. In the literature, \eqref{Eq:SVD} is known as the singular-value decomposition (SVD)   and the entries $\sigma_i$  are the singular values of ${\bf A}$. In general, the SVD of a matrix $\bf A$ is not unique.  However, the diagonal matrix ${\bf S}$ and, consequently, its entries, are fully  determined  from ${\bf A}$. In other words, the values of $\sigma_i$ are invariant to a specific choice of decomposition. This is why one can refer to the diagonal entries of ${\bf S}$ as the ``singular values'' of ${\bf A}$.

 When ${\bf A}$ is not full rank, one can obtain a reduced version of \eqref{Eq:SVD}. Indeed, if we denote the rank of ${\bf A}$ by $r$, then we have that
\begin{equation}\label{Eq:ReducedSVD}
{\bf A} = {\bf U}_r {\bf S}_r {\bf V}_r^T,
\end{equation}
where $ {\bf U}_r \in \mathbb{R}^{m\times r}$ and  $ {\bf V}_r \in \mathbb{R}^{n\times r}$   are (sub)-orthogonal matrices such that $ {\bf U}_r^T  {\bf U}_r =  {\bf V}_r^T {\bf V}_r =  {\bf I}_r$ and ${\bf S}_r ={\rm diag}(\boldsymbol{\sigma})$ is a diagonal matrix that contains positive singular values $\boldsymbol{\sigma}=(\sigma_1,\ldots,\sigma_r)\in \mathbb{R}^{r}$ of ${\bf A}$.

Finally, for any  $p\in [1,+\infty]$, the Schatten-$p$ norm of ${\bf A}$ is defined as
\begin{equation} \label{Eq:SchattenNorm}
\|{\bf A}\|_{S_p} = \begin{cases} \left(\sum_{i=1}^r \sigma_i^p\right)^{\frac{1}{p}}, & p< +\infty \\ \sigma_1, & p= +\infty. \end{cases}
\end{equation}

\section{Duality Mapping in Schatten Spaces}
\label{Sec:DualMap}
The dual of the Schatten-$p$ norm is the Schatten-$q$ norm, where $q\in [1,\infty]$ is such that  $\frac{1}{p}+\frac{1}{q}=1$ \cite{bhatia2013matrix}.  This is due to the generalized version of H\"older's inequality for Schatten norms, as stated in Proposition \ref{Prop:Holder}. While this is a known result (see, for example, \cite{Lefki2015StructureTensor}), it is also the basis for the present work, which is the reason why we provide a proof in  \ref{App:Holder}.
\begin{proposition}\label{Prop:Holder}
For any pair $(p,q)\in [1,+\infty]^2$ of H\"older conjugates with $\frac{1}{p}+\frac{1}{q}=1$  and any pair of  matrices ${\bf A},{\bf B}\in \mathbb{R}^{m\times n}$, we have that
\begin{equation}\label{Eq:Holder}
\langle {\bf A}, {\bf B}\rangle = \mathrm{Tr}\left( {\bf A}^T {\bf B} \right) \leq \|{\bf A}\|_{S_p} \|{\bf B}\|_{S_q}.
\end{equation}
\end{proposition}
In Proposition \ref{Prop:HolderSat}, we investigate the case where the H\"older inequality is saturated, in the sense that
\begin{equation}\label{Eq:HolderSat}
{\rm Tr}\left( {\bf A}^T {\bf B} \right) = \|{\bf A}\|_{S_p} \|{\bf B}\|_{S_q}.
\end{equation}
This saturation is central to our work, as it is tightly linked to the notion of duality mapping.
\begin{proposition}\label{Prop:HolderSat}
Let $(p,q)$ be a pair of H\"older conjugates and let ${\bf A},{\bf B}\in \mathbb{R}^{m\times n}$ be a pair of nonzero matrices with reduced SVDs of the form
\begin{equation}\label{Eq:RedSVD}
{\bf A}= {\bf U}_r {\rm diag}(\boldsymbol{\sigma}) {\bf V}_r^T,  \quad {\bf B}= \tilde{\bf U}_{\tilde{r}} {\rm diag}(\boldsymbol{\tilde{\sigma}})\tilde{\bf V}_{\tilde{r}}^T.
\end{equation}
\begin{itemize}
\item If $p\in (1,\infty)$, then the H\"older inequality is saturated if and only if we have that
\begin{equation}\label{Eq:FormBgenP}
{\bf B} = c {\bf U}_r{\rm diag}({\rm J}_p(\boldsymbol{\sigma})){\bf V}_r^T
\end{equation}
or, equivalently,
\begin{equation}\label{Eq:FormAgenP}
{\bf A} = c^{-1} {\bf \tilde{U}}_{\tilde{r}}{\rm diag}({\rm J}_q(\tilde{\boldsymbol{\sigma}})){\bf \tilde{V}}_{\tilde{r}}^T,
\end{equation}
where $c= \frac{\|{\bf B}\|_{S_q}}{\|{\bf A}\|_{S_p}}$ and ${\rm J}_p(\cdot)$ and ${\rm J}_q(\cdot)$ are the duality mappings for  the $\ell_p$ and $\ell_q$ norms, respectively (see \eqref{Eq:DualMapLp}).
\item If $p=1$, then a necessary condition for the saturation of the H\"older inequality is that
\begin{equation}
{\rm rank}({\bf A}) \leq r_1 \leq {\rm rank}({\bf B}),
\end{equation}
where $r_1= {\rm Card}\left(\{i: \tilde{\sigma}_i = \tilde{\sigma}_1\}\right)$ is  the multiplicity of the first singular value of ${\bf B}$. Moreover, if we denote the first $r_1$ singular vectors of ${\bf B}$ in \eqref{Eq:RedSVD} by ${\bf \tilde{U}}_1\in \mathbb{R}^{m\times r_1}$ and ${\bf \tilde{V}}_1\in\mathbb{R}^{n\times r_1}$, then the saturation of the H\"older inequality is equivalent to the existence of a symmetric matrix ${\bf X}\in \mathbb{R}^{r_1\times r_1}$ such that
\begin{align}\label{Eq:FormA1}
{\bf A}  = {\bf \tilde{U}}_1 {\bf X} {\bf \tilde{V}}_1^T.
\end{align}
Finally in the rank-equality case  ${\rm rank}({\bf A}) ={\rm rank}({\bf B})$, we have saturation if and only if
\begin{equation}\label{Eq:FormBrankEq}
{\bf B} = c {\bf U}_r{\bf V}_r^T,
\end{equation}
where $c= \|{\bf B}\|_{S_{\infty}}$ and the matrices ${\bf U}_r$ and ${\bf V}_r$ are defined in \eqref{Eq:RedSVD}.
\end{itemize}
\end{proposition}
\begin{remark}\label{Remark}
The reduced SVD is not unique; there are multiple choices for the sub-orthogonal  matrices in \eqref{Eq:RedSVD}. However, the parametric forms  given in Proposition \ref{Prop:HolderSat} do not depend on  a specific decomposition.
\end{remark}
 The proof of Proposition \ref{Prop:HolderSat} can be found in  \ref{App:HolderSat}. We observe that,  in the case $p\in(1,\infty)$, the saturation of H\"older inequality provides a very tight link between the two matrices: If we know one of them, then the other lies in a one-dimensional ray that is parameterized by the constant $c>0$. However, in the special case $p=1$, the identification is not as simple. There again, for a given  matrix ${\bf B}$,  one can fully characterize the set of admissible matrices ${\bf A}$. However, for the reverse direction,  an additional rank-equality  constraint  is essential to reduce the set of admissible matrices ${\bf B}$ to just   one ray.

Inspired from Proposition \ref{Prop:HolderSat}, we now propose our main result  in Theorem \ref{Thm:p}, where  we explicitly characterize the duality mapping for the Schatten $p$-norms. The proof of Theorem \ref{Thm:p} can be found in  \ref{App:p}.
\begin{theorem}\label{Thm:p}
Let $p,q\in [1,+\infty]$ be a pair of H\"older conjugates with $\frac{1}{p}+\frac{1}{q}=1$ and ${\bf A}\in\mathbb{R}^{m\times n}$ a matrix whose reduced SVD is specified in \eqref{Eq:ReducedSVD}.
\begin{itemize}
\item If $1<p<+\infty$, then the single-valued duality mapping ${\rm J}_{S_p}: \mathbb{R}^{m\times n} \rightarrow \mathbb{R}^{m\times n}$ is well-defined and can be expressed as
\begin{equation}
{\rm J}_{S_p}:{\bf A} = {\bf U }_r {\rm diag}(\boldsymbol{\sigma} ) {\bf V}_r^T \mapsto {\bf A}^*= {\bf U}_r{\rm diag}({\rm J}_p(\boldsymbol{\sigma} )){\bf V}_r^T.
\end{equation}
\item If $p=1$ and if we consider the rank function as the sparsity measure in Definition \ref{Def:SparseDual}, then the sparse duality mapping ${\rm J}_{S_1,{\rm rank}}: \mathbb{R}^{m\times n} \rightarrow \mathbb{R}^{m\times n}$ is well-defined (singleton)  and is given as
\begin{equation}
{\rm J}_{S_1,{\rm rank}}:{\bf A} = {\bf U }_r {\rm diag}(\boldsymbol{\sigma}) {\bf V}_r^T \mapsto {\bf A}^*= \|\boldsymbol{\sigma}\|_{1} {\bf U} _r{\bf V}_r^T.
\end{equation}
\item If $p=+\infty$, then the set-valued mapping $\mathcal{J}_{S_\infty}(\cdot)$  can be described as
\begin{equation}\label{Eq:Jinf}
\mathcal{J}_{S_\infty}({\bf A})= \left\{  \sigma_1 {\bf U}_1 {\bf X} {\bf V}_1^T:  {\bf X}\in\mathbb{R}^{r_1\times r_1} \text{ is symmetric and } \|{\bf X}\|_{S_1} =1 \right\},
\end{equation}
where $r_1$ denotes the multiplicity of the first singular value $\sigma_1$ of ${\bf A}$ and ${\bf U}_1,{\bf V}_1$ are singular vectors that correspond to $\sigma_1$ in \eqref{Eq:ReducedSVD}. It is a convex set whose extreme points  are ${\bf E}_{i,j} = \frac{\sigma_1}{2}({\bf u}_i {\bf v}_j^T +{\bf v}_i {\bf u}_j^T )$ for $1\leq i\leq j\leq r_1$. Finally, the set of sparse dual conjugates is the collection of rank-1 elements of $\mathcal{J}_{S_\infty}({\bf A})$ which can be characterized as
\begin{equation}
\mathcal{J}_{S_\infty, {\rm rank}}({\bf A}) = \{ \sigma_1 {\bf U}_1 {\bf p}{\bf p}^T {\bf V}_1^T: {\bf p}\in \mathbb{R}^{r_1} , \|{\bf p}\|_2=1\}.
\end{equation}
\end{itemize}
\end{theorem}

\section{Discussion}\label{Sec:discuss}
Theorem \ref{Thm:p} provides an interesting characterization of the duality mapping in three scenarios: The first case is $1<p<+\infty$ which is the most straightforward one. Theorem \ref{Thm:p} tells us  that the mapping is single-valued and also gives a formula to compute the dual conjugate ${\bf A}^*$ of any matrix ${\bf A}\in\mathbb{R}^{m\times n}$.  We use this result to deduce the continuity of the duality mapping as well as the strict convexity of the Schatten space in this case (see Corollary \ref{Cor:strict}).  In the second case, with $p=1$, the mapping is not single-valued. However, there is a unique element in the set of dual conjugates with the minimal rank (that is equal to the rank of ${\bf A}$) and, hence, we can construct a single-valued sparse duality mapping. Finally, we showed in the third case, characterized by  $p=+\infty$, that neither the set of dual conjugates nor the ones with the minimal rank are   unique. However, we observe in \eqref{Eq:Jinf} that the entries of ${\bf X}$ can be independently chosen, up to symmetry and normalization assumptions. This suggests that the  dimension of $\mathcal{J}_{S_\infty}({\bf A})$ is $d=\left( \frac{r_1(r_1+1)}{2}-1\right)$. Moreover, we show that this convex set has exactly $(d+1)$ extreme points, which is the minimal number for a convex set of dimension $d$. We also observe that the extreme points of $\mathcal{J}_{S_\infty}({\bf A})$ are low-rank. They are indeed a collection of rank-1 and rank-2 matrices.

In Corollary \ref{Cor:strict}, we highlight some consequences of Theorem \ref{Thm:p} concerning the strict convexity of Schatten spaces and the continuity of the duality mapping.
\begin{corollary}\label{Cor:strict}
The Banach space of $m$ by $n$ matrices equipped with the Schatten-$p$ norm is strictly convex, if and only if $p\in (1,+\infty)$. In this case,   the function ${\rm J}_{S_p} : \mathbb{R}^{m\times n} \rightarrow \mathbb{R}^{m\times n}$ is continuous.
\end{corollary}
\begin{proof}
For $p\in (1,+\infty)$, we know from Theorem \ref{Thm:p} that the duality mapping ${\rm J}_{S_p}$ is bijective. Moreover, it is known that all finite-dimensional Banach spaces are reflexive. Now, following \cite{petryshyn1970characterization}, we deduce the strict convexity of the space of $m$ by $n$ matrices with Schatten-$p$ norm.

For $p=1$ and $p=+\infty$, we can readily verify that
\begin{align*}
&\left \|\alpha \begin{pmatrix}1 & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & 0 \end{pmatrix} +(1-\alpha)\begin{pmatrix}0 & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & 1 \end{pmatrix}  \right\|_{S_1} = \left \|\begin{pmatrix} \alpha & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & (1-\alpha) \end{pmatrix} \right\|_{S_1}=1, \\
&\left\|\alpha \begin{pmatrix}1 & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & 0 \end{pmatrix}  +(1-\alpha)\begin{pmatrix}1 & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & 1 \end{pmatrix}  \right\|_{S_\infty} = \left\| \begin{pmatrix}1 & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & (1-\alpha) \end{pmatrix}  \right\|_{S_\infty}=1,
\end{align*}
for all $\alpha \in (0,1)$, which  shows that the Schatten space is not strictly convex for $p=1,+\infty$.

Finally, the Schatten-$p$ norm is  known to be  Fr\'echet differentiable  for $p\in(1,+\infty)$ \cite{potapov2014frechet}. Moreover,  the duality mapping of any  Banach space with Fr\'echet-differentiable norms is guaranteed to be continuous \cite{giles1978geometrical,contreras1994upper}. Combining the two statements, we deduce the continuity of the duality mapping in this case.
\end{proof}
By contrast, the sparse duality mapping ${\rm J}_{S_1,{\rm rank}}(\cdot)$ is not continuous. This is best explained by providing a  counterexample. Specifically, let us consider the sequence of 2 by 2 matrices
$${\bf S}_k=\begin{pmatrix} 1 & 0\\ 0 & \frac{1}{k}\end{pmatrix}, \quad k\in\mathbb{N}.$$
It is clear that ${\bf S}_k \rightarrow {\bf S}_{\infty} =\begin{pmatrix} 1 & 0\\ 0 & 0\end{pmatrix}$. However, we   have that
$$\forall k\in\mathbb{N}: {\rm J}_{S_1,{\rm rank}}({\bf S}_k) = \begin{pmatrix} 1 & 0\\ 0 & 1\end{pmatrix}, \quad \text{while} \quad {\rm J}_{S_1,{\rm rank}}({\bf S}_{\infty}) = \begin{pmatrix} 1 & 0\\ 0 & 0\end{pmatrix},$$
which shows the discontinuity of ${\rm J}_{S_1,{\rm rank}}$ in the space of 2 by 2 matrices. This can be generalized to space of matrices with arbitrary dimensions $m,n\in\mathbb{N}$.

Although ${\rm J}_{S_1,{\rm rank}}$ is not continuous, we  now show that it is Borel-measurable and, hence, that it can be approximated with arbitrary precision by a continuous mapping due to Lusin's theorem \cite{rudin1991functional}.
\begin{proposition}\label{Prop:measure}
For any $m,n\in\mathbb{N}$, the sparse duality mapping ${\rm J}_{S_1,{\rm rank}}$ is  a Borel-measurable matrix-valued function over the space of $m$ by $n$ matrices.
\end{proposition}
Before going into the proof of Proposition \ref{Prop:measure}, we present a preliminary result.
\begin{lemma}\label{Lem}
The set  $\mathcal{R}_{ r}\subseteq \mathbb{R}^{m\times n}$  of $m$ by $n$  matrices of rank $r$ is Borel-measurable.
\end{lemma}
\begin{proof}
First note that
$$\mathcal{R}_{1}=\{ {\bf u}{\bf v}^T: {\bf u}\in\mathbb{R}^m, {\bf v}\in \mathbb{R}^n\}. $$
The set $\mathcal{R}_{1}$ is the image of the continuous mapping $\mathbb{R}^{m}\times \mathbb{R}^n \rightarrow \mathbb{R}^{m\times n}:({\bf u}, {\bf v}) \mapsto  {\bf u}{\bf v}^T$ and, hence, is Borel-measurable.

Now, denote by  $\mathcal{R}_{\leq r}\subseteq \mathbb{R}^{m\times n}$, the set of matrices with rank no more than $r$. Using the identity
$$\mathcal{R}_{\leq r} =  \mathcal{R}_{1} + \cdots + \mathcal{R}_{1} ,\quad \text{(r times)}, $$
we deduce that $\mathcal{R}_{\leq r}$ and, consequently, $\mathcal{R}_{r} =  \mathcal{R}_{\leq r}  \backslash \mathcal{R}_{\leq (r-1)}$ are also Borel-measurable sets.
\end{proof}
\begin{proof}[Proof of Proposition \ref{Prop:measure}]
Consider a Borel-measurable set $\mathcal{B}\subseteq \mathbb{R}^{m\times n}$. We show that $\mathcal{B}_{\rm inv}= {\rm J}_{S_1,{\rm rank}}^{-1}(\mathcal{B})$ is also Borel-measurable. By defining  $\mathcal{B}_{{\rm inv},r} = \mathcal{B}_{\rm inv}\cap \mathcal{R}_{r}$, we can partition   $ \mathcal{B}_{\rm inv}$ as
$$\mathcal{B}_{\rm inv}= \bigcup_{r=1}^{\min (m,n)}  \mathcal{B}_{\rm inv}\cap \mathcal{R}_{r}.$$
Hence, it is sufficient to show that each partition $\mathcal{B}_{\rm inv}\cap \mathcal{R}_{r}$ is Borel-measurable.

Define the set $\mathcal{P}_{r} \subseteq \mathcal{R}_{r}^2$ as
$$\mathcal{P}_{r} = \{ ({\bf A},{\bf B}) \in  \mathcal{R}_{r} \times \mathcal{B}: {\rm Tr}({\bf A}^T{\bf B}) = \|{\bf A}\|_{S_1} \|{\bf B}\|_{S_\infty}, \quad \|{\bf A}\|_{S_1}= \|{\bf B}\|_{S_\infty}\}.$$
The set $\mathcal{P}_{r}$ introduces a relation over $ \mathcal{R}_{r} $ whose domain is $\mathcal{B}_{\rm inv}\cap \mathcal{R}_{r}$. Since the trace and norm are continuous (and, consequently, Borel-measurable) functions and  $\mathcal{R}_{r}\times \mathcal{B}$ is a Borel-measurable set (using Lemma \ref{Lem}), we deduce that the relation induced from $\mathcal{P}_{r}$ is Borel-measurable as well. Finally, we use \cite[Proposition 2.1]{himmelberg1975measurable} to show that its domain is Borel-measurable.
\end{proof}

\section{Conclusion}
In this paper, we studied the duality mapping in finite-dimensional Schatten spaces. Based on a careful investigation of the cases where the  H\"older inequality saturates, we provided an explicit form for this mapping when $p\in (1,+\infty)$. Furthermore,  by adding a rank constraint, we proved that the mapping becomes single-valued for the special case $p=1$. As for $p=+\infty$, we showed that the mapping yields a convex set whose extreme points are low-rank matrices. Finally, we discussed our theorem and studied the continuity of the introduced mappings as well as the strict convexity of the Schatten spaces. A possible future direction of research is to extend the results of this paper to infinite-dimensional Schatten spaces  and even, in full generality, to linear operators over Hilbert spaces.