\section{Introduction}

Let $n$ be a positive integer and let $\mathbb{Z}_n$ be the ring of integers modulo $n.$ In \cite{gri}, Ralph P. Grimaldi defined a graph $G(\mathbb{Z}_n)$ based on the elements and units of $\mathbb{Z}_n.$ The vertices of the unit graph $G(\mathbb{Z}_n)$ are the elements of $\mathbb{Z}_n$ and distinct vertices $x$ and $y$ are defined to be adjacent if and only if $x + y$ is a unit of $\mathbb{Z}_n.$ That is, $xy$ is an edge if and only if $x + y$ is a unit in $\mathbb{Z}_n.$ For a positive integer $m,$ it follows that $G(\mathbb{Z}_{2m})$ is a $\phi(2m)$-regular
graph where $\phi$ is the Euler phi function.

Codes generated by the incidence matrices and the adjacency matrices of
various  graphs are discussed in \cite{3, cd1, cd2, anna}.

In what follows, all rings are associative with nonzero identity, denoted
by 1, which is preserved by ring homomorphisms and inherited by subrings.
Also throughout the article, by a graph $G$ we mean a finite undirected graph
without loops or multiple edges (unless otherwise specified).

\section{Preliminaries }
In this section, we study the basic properties of unit graphs and linear codes.

\indent Let $\mathbb{F}_q$ denote the finite field with cardinality $q.$ Then $\mathbb{F}_q^n$ is  a $n$-dimensional vector space over the finite field $\mathbb{F}_q.$  The Hamming weight
$w_{H}(a$) of $a\in \mathbb{F}_q$ is \begin{align*}
	w_{H}(a)=\begin{cases}0&\,\,\text{ for }\,\,a=0\\1&\,\,\text{ for }\,\,
		a\neq 0.
	\end{cases}
\end{align*}
Let  $a,b\in\mathbb{F}_q$, then the Hamming distance of $a$ and $b$ is defined to be $d_H(a,b)=w_{H}(a-b).$
 For a word $x=(x_1, \dots, x_n)\in \mathbb{F}_q^n,$ the Hamming weight $w_{H}(x)$ is defined to be $$w_{H}(x)=\sum\limits_{i=1}^nw_{H}(x_{i}).$$
Let $x=(x_1,\dots, x_n), y=(y_1, \dots, y_n) \in \mathbb{F}_q^n,$ then the  Hamming distance between $x$ and $y$ is defined by
$$d_H(x,y)=\sum\limits_{i=1}^{n}d_H(x_i,y_i)=\sum\limits_{i=1}
^{n}w_H(x_i-y_i).$$

A non-empty subset $C$ of $\mathbb{F}_{q}^{n}$ is said to be a {\it $q$-ary code of length $n.$} An element of the code
$C$ is called a { \it codeword.} A $q$-ary linear code $C$ of length $n$  is a
subspace of the
vector space $\mathbb{F}_{q}^{n}$ over
$\mathbb{F}_{q}.$\\
\indent The minimum Hamming distance of a code $C$ is defined by
$$
d_{H}(C)=\min\limits_{c_{1},c_{2}\in C}\{
d_{H}(c_{1},
c_{2}) \mid c_1\neq c_2\}.$$
\indent The minimum weight of a code $C$ is the smallest among all weights of the non-zero codewords of $C.$ For $q$-ary linear code, we have
$d_{H}(C)=w_{H}(C).$

All the codes here are linear codes and the notation $[n, k, d]_{q}$
will be used for a $q$-ary code of length $n,$ dimension $k$ and minimum
distance $d.$ A linear $[n, k, d]_{q}$  code  is said to be a {\it Maximum Distance Separable(MDS) code} if $d=n -k + 1.$ A generator matrix $B$ for a
linear code $C$ is a matrix whose rows form a basis for the subspace $C$ and $C_{q}(B)$ is a code generated by matrix $B$ over a finite field $\mathbb{F}_{q}.$

 Let $G=(V, E)$  be a graph with vertex set $V,$ edge set $E$ and for any $x, y \in V,\,[x,y]$ is denoted by an edge between $x$ and $y.$ An incidence matrix of
$G$ is a $|V|\times|E|$ matrix $H$ with rows labelled by the vertices and columns by the edges and entries $h_{ij}=1$ if $i$th vertex is adjacent to the $j$th edge and $h_{ij}=0$ otherwise.

Let $R$ be a ring and let $U(R)$ be the set of unit elements of $R.$ The unit graph of $R,$ denoted $G(R),$ is the graph obtained by setting all the elements of $R$ to be the vertices and defining distinct vertices $x$ and $y$ to be
adjacent if and only if $x + y \in U(R).$

Let $G$ be a graph. We denote $V(G)$ as the vetex set of $G$ and  $E(G)$ as the edge set of $G.$ Let $x\in V(G),$ then the degree of $x,$ denoted $deg(x),$ is the number of edges of $G$ incident with $x.$ A graph $G$ is said to be a { \it $r$-regular } if the degree of each vertex of $G$ is equal to $r.$

Let $W, X\subseteq V$ with $W\cap X=\emptyset$ and let  $E(W, X)$ be the set of edges that have one end in $W$ and the other end in $X.$ Write $|E(W, X)|=q(W, X).$

An edge-cut of a connected graph $G$ is the set $S\subseteq E$ such that $G- S=(V, E-S)$ is disconnected.

The edge-connectivity $\lambda(G)$ is the minimum cardinality of an edge-cut. That is, $$\lambda(G)=\min\limits_{\emptyset \neq W\subsetneq V} q(W, V-W).$$
For any connected graph $G,$ we have $\lambda(G)\leq \delta(G)$
where $\delta(G)$ is minimum degree of the graph $G.$
\begin{theorem}\cite{ash}\label{thm1}
Let $R$ be a finite ring. Then the following statements hold for the
unit graph of $R.$
\begin{enumerate}
\item[(a)] If $2\notin U(R),$ then the unit graph $G(R)$ is a $|U(R)|$-regular graph.
\item[(b)] If $2\in U(R),$ then for every $x\in U(R)$ we have $deg(x)= |U(R)|- 1$ and for every $x\in R\setminus U(R)$ we have $deg(x)=|U(R)|.$
\end{enumerate}
\end{theorem}

\begin{theorem}\cite{3} \label{thm2}
 Let $G= (V, E)$ be a connected graph and let  $H$ be a $|V|\times|E|$ incidence matrix for $G.$ Then
		\begin{itemize}
		\item[1.] the code $C_2(H)$ is a $[|E|, |V|- 1, \lambda(G)]_2$ code
			\item[2.]  the code
		$C_p(H)$ is $[|E|, |V|- 1, \lambda(G)]_p$ code for an odd prime $p$  and  $G$ is bipartite.
	\end{itemize}
\end{theorem}
In this paper, we obtain a linear codes from the incidence matrix of the unit graph
 $G(\mathbb{Z}_n)$ over $\mathbb{F}_2$ and we determined the main parameters of the code. In section 3, we discussed codes from incidence matrix of a unit graph $G(\mathbb{Z}_{p_1})$ over the field $\mathbb{F}_2.$ In section 4, we discussed codes from incidence matrix of a unit graph $G(\mathbb{Z}_{2p_1})$ over the finite field $\mathbb{F}_2.$ All the unit graphs considered in this article is a simple and undirected graph.
\section{Linear Codes from the Incidence Matrix of a unit Graph $G(\mathbb{Z}_{p_1})$}
In this section, we  study a binary linear code obtained from the incidence matrix of the unit graph $G(\mathbb{Z}_{p_1})$ where $p_1$ is an odd prime and find its parameters.

 The units of $\mathbb{Z}_{p_1}$  is $U(\mathbb{Z}_{p_1})=\{1, 2, \cdots, p_1-1\}.$
Let $G(\mathbb{Z}_{p_1})$ be the unit graph with vertex set $V=\mathbb{Z}_{p_1}$ and two distinct vertices $x$ and $y$ are adjacent if and only if $x+y\in U(\mathbb{Z}_{p_1}).$
\begin{theorem}\label{thm3}
Let $G(\mathbb{Z}_{p_1})=(V, E)$ be a unit graph. Then the graph is connected with	$|V|=p_1$ and $|E|=\dfrac{(p_1-1)^2}{2}.$
\end{theorem}
\begin{proof}
Let $G(\mathbb{Z}_{p_1})$ be a unit graph. Then by definition, $V=\mathbb{Z}_{p_1}.$ Since all non-zero elements are units, $0$ and $y$ are adjacent for all nonzero $y \in \mathbb{Z}_{p_1}.$ Therefore the graph is a connected graph.
\begin{align*}
 \Big|\Big\{[0, y] \mid y\in \mathbb{Z}_{p_1}\setminus \{0\} \Big\}\Big|&=p_1-1\\
\Big|\big\{[j, y] \mid y\in \mathbb{Z}_{p_1}\setminus \{0, 1, \cdots,j, p_1-j\} \big\}\Big|&=p_1-(j+2) \text{ for } 1\leq j\leq \frac{p_1-1}{2}
\end{align*}
For $\dfrac{p_1+1}{2}\leq j\leq p_1-2,$ we have
$\Big|\big\{[j, y] \mid j<y \big\}\Big|=p_1-(j+1).$
Therefore,  $$|E|=p_1-1+\sum\limits_{j=1}^{\dfrac{p_1-1}{2}}p_1-(j+2)+\sum\limits_{j=1}^{\dfrac{p_1-3}{2}}j=\frac{(p_1-1)^2}{2}.$$ Hence, the unit graph $G(\mathbb{Z}_{p_1})$  is a connected graph with $p_1$ vertices and $\dfrac{(p_1-1)^2}{2}$ edges.
\end{proof}

\begin{theorem}\label{thm4}
	Let $G(\mathbb{Z}_{p_1})$ be the unit graph. Then the edge-connectivity $\lambda(G(\mathbb{Z}_{p_1}))$ of the unit graph $G(\mathbb{Z}_{p_1})$ is $p_1-2.$
\end{theorem}
\begin{proof}
If $W=\{0\}\subset V,$ then $q(W, V-W)=p_1-1.$ If $W=\{a\}\subset V, a\neq 0,$ then $q(W, V-W)=p_1-2.$
If $W\subset V$ with $ |W|>1,$ then $q(W, V-W)>p_1-2.$ Therefore, $\lambda(G(\mathbb{Z}_{p_1}))=p_1-2.$
\end{proof}
As a consequence of above theorems, we have
\begin{theorem}\label{thm5}
Let $p_1$ be an odd prime. Then the code generated by the incidence matrix $H$ of the unit graph $G(\mathbb{Z}_{p_1})$ is a  $C_2(H)=\Big[\dfrac{(p_1-1)^2}{2}, p_1-1, p_1-2\Big]_2$ code over the finite field $\mathbb{F}_2.$
\end{theorem}

\begin{proof}
	Let  $G(\mathbb{Z}_{p_1})$ be a unit graph and let $H$ be the incidence matrix of $G(\mathbb{Z}_{p_1}).$ Since  $G(\mathbb{Z}_{p_1})$ is a connected graph, by Theorem \ref{thm2}, the code $C_2(H)$ is a $[|E|, |V|-1, \lambda(G(\mathbb{Z}_{p_1}))]_2$ code. By Theorem \ref{thm3}, we have  $|E|=\dfrac{(p_1-1)^2}{2}$ and $|V|=p_1.$ By Theorem \ref{thm4}, the edge-connective of the unit graph $G(\mathbb{Z}_{p_1})$ is $p_1-2$ and hence $C_2(H)$ is a $\Big[\dfrac{(p_1-1)^2}{2}, p_1-1, p_1-2\Big]_2$ code.
\end{proof}
\begin{example}
	The unit graph $G(\mathbb{Z}_{5})=(V, E)$ with $|V|=5$ and $|E|=\frac{(5-1)^2}{2}=8$ is
	\begin{center}
		\includegraphics{5.pdf}
	\end{center}
		Then the incidence matrix of the unit graph $G(\mathbb{Z}_{5})$ is
	$$H=\begin{pmatrix}
		1&1&1&1&\vline&0&0&0&0\\
		\hline
		1&0&0&0&\vline&1&1&0&0\\
		0&1&0&0&\vline&1&0&1&0\\
		0&0&1&0&\vline&0&1&0&1\\
		0&0&0&1&\vline&0&0&1&1
	\end{pmatrix}_{5\times 8}$$
Any four rows of $H$ are linearly independent and the graph is connected, then the dimension of the code $C_2(H)$ over the finite field $\mathbb{F}_2$ generated by $H$ is $4.$ The minimum distance of the code $C_2(H)$ is $3.$ Hence $C_2(H)$ is an $[n, k, d]_2=[8, 4, 3]_2$ code.
\end{example}

\section{Linear Codes from the Incidence Matrix of a unit Graph $G(\mathbb{Z}_{2p_1})$}
In this section, we  study the linear codes obtained from the incidence matrix of the unit graph $G(\mathbb{Z}_{2p_1})$ where $p_1$ is an odd prime over $\mathbb{F}_2$ and we find the parameters of the code.

The units of $\mathbb{Z}_{2p_1}$  is $U(\mathbb{Z}_{2p_1})=\{k\in \mathbb{Z}_{2p_1} \mid (k, 2p_1)=1\}$ and $|U(\mathbb{Z}_{2p_1})|=p_1-1.$

\begin{theorem}\label{thm4}
	Let $G(\mathbb{Z}_{2p_1})=(V, E)$ be a unit graph. Then 	$|V|=2p_1$ and $|E|=p_1(p_1-1).$
\end{theorem}
\begin{proof}
	Let $G(\mathbb{Z}_{p_1})$ be a unit graph. Then by definition, $V=\mathbb{Z}_{2p_1}.$
	Since $(2, 2p_1)=2\neq 1,$ by Theorem \ref{thm1}, $G(\mathbb{Z}_{2p_1})$ is a $|U(\mathbb{Z}_{2p_1})|=(p_1-1)$-regular graph. That is, degree of every vertex is $p_1-1.$
	Since the number of edges of the $k$-regular graph with $n$ vertices is $\dfrac{nk}{2},$ implies $|E|=p_1(p_1-1).$
\end{proof}

\begin{theorem}\label{thm6}
	Let $p_1$ be an odd  prime. Then the code generated by the incidence matrix $H$ of the unit graph $G(\mathbb{Z}_{2p_1})$ is a  $C_2(H)=\Big[p_1(p_1-1), 2p_1-1, p_1-1\Big]_2$ code over the finite field $\mathbb{F}_2.$
\end{theorem}
\begin{proof}
	Let  $G(\mathbb{Z}_{2p_1})$ be a unit graph and let $H$ be the incidence matrix of $G(\mathbb{Z}_{2p_1}).$ Since  $G(\mathbb{Z}_{2p_1})$ is a connected graph, by Theorem \ref{thm2}, the code $C_2(H)$ is a $[|E|, |V|-1, \lambda(G(\mathbb{Z}_{2p_1}))]_2$ code. Since  $G(\mathbb{Z}_{2p_1})$ is $(p_1-1)$-regular graph, the edge-connectivity $\lambda(G(\mathbb{Z}_{2p_1}))$  of the unit graph $G(\mathbb{Z}_{2p_1})$ is $p_1-1$ and hence the minimum distance of the code generated by $H$ is $p_1 -1.$  Since  $|E|=p_1(p_1-1)$ and $|V|=2p_1,$ the main parameters of the code $C_2(H)$ is  $[p_1(p_1-1), 2p_1-1, p_1-1]_2.$
\end{proof}
\begin{example}
	The unit graph $G(\mathbb{Z}_{6})=(V, E)$ with $|V|=6$ and $|E|=6$ is
	\begin{center}
		\includegraphics{6.pdf}
	\end{center}
	Then the incidence matrix of the unit graph $G(\mathbb{Z}_{6})$ is
	$$H=\begin{pmatrix}
		1&0&0&0&0&0&1\\
		1&1&0&0&0&0&0\\
		0&0&1&1&0&0&0\\
		0&0&0&1&1&0&0\\
		0&1&0&0&0&1&0\\
		0&0&1&0&0&0&1
	\end{pmatrix}_{6\times 6}$$
	Any five rows of $H$ are linearly independent and the graph is connected, then the dimension of the code $C_2(H)$ over the finite field $\mathbb{F}_2$ generated by $H$ is $5.$ The minimum distance of the code $C_2(H)$ is $2.$ Hence $C_2(H)$ is an $[n, k, d]_2=[6, 5, 2]_2$ code. Since $d=n-k+1,$ the code $C_2(H)$ is a MDS code.
\end{example}

\section*{Conclusion}
In this paper, we studied the codes generated by the incidence matrix of the unit graph of the different commutative rings with unity. Also we found the main parameters of the code over finite field $\mathbb{F}_2.$ We have consider only simple and undirected graphs in this article. Finding the covering radius of these codes is the further direction to work.

\section{Introduction}

Let $n$ be a positive integer and let $\mathbb{Z}_n$ be the ring of integers modulo $n.$ In \cite{gri}, Ralph P. Grimaldi defined a graph $G(\mathbb{Z}_n)$ based on the elements and units of $\mathbb{Z}_n.$ The vertices of the unit graph $G(\mathbb{Z}_n)$ are the elements of $\mathbb{Z}_n$ and distinct vertices $x$ and $y$ are defined to be adjacent if and only if $x + y$ is a unit of $\mathbb{Z}_n.$ That is, $xy$ is an edge if and only if $x + y$ is a unit in $\mathbb{Z}_n.$ For a positive integer $m,$ it follows that $G(\mathbb{Z}_{2m})$ is a $\phi(2m)$-regular
graph where $\phi$ is the Euler phi function.

Codes generated by the incidence matrices and the adjacency matrices of
various  graphs are discussed in \cite{3, cd1, cd2, anna}.

In what follows, all rings are associative with nonzero identity, denoted
by 1, which is preserved by ring homomorphisms and inherited by subrings.
Also throughout the article, by a graph $G$ we mean a finite undirected graph
without loops or multiple edges (unless otherwise specified).

\section{Preliminaries }
In this section, we study the basic properties of unit graphs and linear codes.

\indent Let $\mathbb{F}_q$ denote the finite field with cardinality $q.$ Then $\mathbb{F}_q^n$ is  a $n$-dimensional vector space over the finite field $\mathbb{F}_q.$  The Hamming weight
$w_{H}(a$) of $a\in \mathbb{F}_q$ is \begin{align*}
	w_{H}(a)=\begin{cases}0&\,\,\text{ for }\,\,a=0\\1&\,\,\text{ for }\,\,
		a\neq 0.
	\end{cases}
\end{align*}
Let  $a,b\in\mathbb{F}_q$, then the Hamming distance of $a$ and $b$ is defined to be $d_H(a,b)=w_{H}(a-b).$
 For a word $x=(x_1, \dots, x_n)\in \mathbb{F}_q^n,$ the Hamming weight $w_{H}(x)$ is defined to be $$w_{H}(x)=\sum\limits_{i=1}^nw_{H}(x_{i}).$$
Let $x=(x_1,\dots, x_n), y=(y_1, \dots, y_n) \in \mathbb{F}_q^n,$ then the  Hamming distance between $x$ and $y$ is defined by
$$d_H(x,y)=\sum\limits_{i=1}^{n}d_H(x_i,y_i)=\sum\limits_{i=1}
^{n}w_H(x_i-y_i).$$

A non-empty subset $C$ of $\mathbb{F}_{q}^{n}$ is said to be a {\it $q$-ary code of length $n.$} An element of the code
$C$ is called a { \it codeword.} A $q$-ary linear code $C$ of length $n$  is a
subspace of the
vector space $\mathbb{F}_{q}^{n}$ over
$\mathbb{F}_{q}.$\\
\indent The minimum Hamming distance of a code $C$ is defined by
$$
d_{H}(C)=\min\limits_{c_{1},c_{2}\in C}\{
d_{H}(c_{1},
c_{2}) \mid c_1\neq c_2\}.$$
\indent The minimum weight of a code $C$ is the smallest among all weights of the non-zero codewords of $C.$ For $q$-ary linear code, we have
$d_{H}(C)=w_{H}(C).$

All the codes here are linear codes and the notation $[n, k, d]_{q}$
will be used for a $q$-ary code of length $n,$ dimension $k$ and minimum
distance $d.$ A linear $[n, k, d]_{q}$  code  is said to be a {\it Maximum Distance Separable(MDS) code} if $d=n -k + 1.$ A generator matrix $B$ for a
linear code $C$ is a matrix whose rows form a basis for the subspace $C$ and $C_{q}(B)$ is a code generated by matrix $B$ over a finite field $\mathbb{F}_{q}.$

 Let $G=(V, E)$  be a graph with vertex set $V,$ edge set $E$ and for any $x, y \in V,\,[x,y]$ is denoted by an edge between $x$ and $y.$ An incidence matrix of
$G$ is a $|V|\times|E|$ matrix $H$ with rows labelled by the vertices and columns by the edges and entries $h_{ij}=1$ if $i$th vertex is adjacent to the $j$th edge and $h_{ij}=0$ otherwise.

Let $R$ be a ring and let $U(R)$ be the set of unit elements of $R.$ The unit graph of $R,$ denoted $G(R),$ is the graph obtained by setting all the elements of $R$ to be the vertices and defining distinct vertices $x$ and $y$ to be
adjacent if and only if $x + y \in U(R).$

Let $G$ be a graph. We denote $V(G)$ as the vetex set of $G$ and  $E(G)$ as the edge set of $G.$ Let $x\in V(G),$ then the degree of $x,$ denoted $deg(x),$ is the number of edges of $G$ incident with $x.$ A graph $G$ is said to be a { \it $r$-regular } if the degree of each vertex of $G$ is equal to $r.$

Let $W, X\subseteq V$ with $W\cap X=\emptyset$ and let  $E(W, X)$ be the set of edges that have one end in $W$ and the other end in $X.$ Write $|E(W, X)|=q(W, X).$

An edge-cut of a connected graph $G$ is the set $S\subseteq E$ such that $G- S=(V, E-S)$ is disconnected.

The edge-connectivity $\lambda(G)$ is the minimum cardinality of an edge-cut. That is, $$\lambda(G)=\min\limits_{\emptyset \neq W\subsetneq V} q(W, V-W).$$
For any connected graph $G,$ we have $\lambda(G)\leq \delta(G)$
where $\delta(G)$ is minimum degree of the graph $G.$
\begin{theorem}\cite{ash}\label{thm1}
Let $R$ be a finite ring. Then the following statements hold for the
unit graph of $R.$
\begin{enumerate}
\item[(a)] If $2\notin U(R),$ then the unit graph $G(R)$ is a $|U(R)|$-regular graph.
\item[(b)] If $2\in U(R),$ then for every $x\in U(R)$ we have $deg(x)= |U(R)|- 1$ and for every $x\in R\setminus U(R)$ we have $deg(x)=|U(R)|.$
\end{enumerate}
\end{theorem}

\begin{theorem}\cite{3} \label{thm2}
 Let $G= (V, E)$ be a connected graph and let  $H$ be a $|V|\times|E|$ incidence matrix for $G.$ Then
		\begin{itemize}
		\item[1.] the code $C_2(H)$ is a $[|E|, |V|- 1, \lambda(G)]_2$ code
			\item[2.]  the code
		$C_p(H)$ is $[|E|, |V|- 1, \lambda(G)]_p$ code for an odd prime $p$  and  $G$ is bipartite.
	\end{itemize}
\end{theorem}
In this paper, we obtain a linear codes from the incidence matrix of the unit graph
 $G(\mathbb{Z}_n)$ over $\mathbb{F}_2$ and we determined the main parameters of the code. In section 3, we discussed codes from incidence matrix of a unit graph $G(\mathbb{Z}_{p_1})$ over the field $\mathbb{F}_2.$ In section 4, we discussed codes from incidence matrix of a unit graph $G(\mathbb{Z}_{2p_1})$ over the finite field $\mathbb{F}_2.$ All the unit graphs considered in this article is a simple and undirected graph.
\section{Linear Codes from the Incidence Matrix of a unit Graph $G(\mathbb{Z}_{p_1})$}
In this section, we  study a binary linear code obtained from the incidence matrix of the unit graph $G(\mathbb{Z}_{p_1})$ where $p_1$ is an odd prime and find its parameters.

 The units of $\mathbb{Z}_{p_1}$  is $U(\mathbb{Z}_{p_1})=\{1, 2, \cdots, p_1-1\}.$
Let $G(\mathbb{Z}_{p_1})$ be the unit graph with vertex set $V=\mathbb{Z}_{p_1}$ and two distinct vertices $x$ and $y$ are adjacent if and only if $x+y\in U(\mathbb{Z}_{p_1}).$
\begin{theorem}\label{thm3}
Let $G(\mathbb{Z}_{p_1})=(V, E)$ be a unit graph. Then the graph is connected with	$|V|=p_1$ and $|E|=\dfrac{(p_1-1)^2}{2}.$
\end{theorem}
\begin{proof}
Let $G(\mathbb{Z}_{p_1})$ be a unit graph. Then by definition, $V=\mathbb{Z}_{p_1}.$ Since all non-zero elements are units, $0$ and $y$ are adjacent for all nonzero $y \in \mathbb{Z}_{p_1}.$ Therefore the graph is a connected graph.
\begin{align*}
 \Big|\Big\{[0, y] \mid y\in \mathbb{Z}_{p_1}\setminus \{0\} \Big\}\Big|&=p_1-1\\
\Big|\big\{[j, y] \mid y\in \mathbb{Z}_{p_1}\setminus \{0, 1, \cdots,j, p_1-j\} \big\}\Big|&=p_1-(j+2) \text{ for } 1\leq j\leq \frac{p_1-1}{2}
\end{align*}
For $\dfrac{p_1+1}{2}\leq j\leq p_1-2,$ we have
$\Big|\big\{[j, y] \mid j<y \big\}\Big|=p_1-(j+1).$
Therefore,  $$|E|=p_1-1+\sum\limits_{j=1}^{\dfrac{p_1-1}{2}}p_1-(j+2)+\sum\limits_{j=1}^{\dfrac{p_1-3}{2}}j=\frac{(p_1-1)^2}{2}.$$ Hence, the unit graph $G(\mathbb{Z}_{p_1})$  is a connected graph with $p_1$ vertices and $\dfrac{(p_1-1)^2}{2}$ edges.
\end{proof}

\begin{theorem}\label{thm4}
	Let $G(\mathbb{Z}_{p_1})$ be the unit graph. Then the edge-connectivity $\lambda(G(\mathbb{Z}_{p_1}))$ of the unit graph $G(\mathbb{Z}_{p_1})$ is $p_1-2.$
\end{theorem}
\begin{proof}
If $W=\{0\}\subset V,$ then $q(W, V-W)=p_1-1.$ If $W=\{a\}\subset V, a\neq 0,$ then $q(W, V-W)=p_1-2.$
If $W\subset V$ with $ |W|>1,$ then $q(W, V-W)>p_1-2.$ Therefore, $\lambda(G(\mathbb{Z}_{p_1}))=p_1-2.$
\end{proof}
As a consequence of above theorems, we have
\begin{theorem}\label{thm5}
Let $p_1$ be an odd prime. Then the code generated by the incidence matrix $H$ of the unit graph $G(\mathbb{Z}_{p_1})$ is a  $C_2(H)=\Big[\dfrac{(p_1-1)^2}{2}, p_1-1, p_1-2\Big]_2$ code over the finite field $\mathbb{F}_2.$
\end{theorem}

\begin{proof}
	Let  $G(\mathbb{Z}_{p_1})$ be a unit graph and let $H$ be the incidence matrix of $G(\mathbb{Z}_{p_1}).$ Since  $G(\mathbb{Z}_{p_1})$ is a connected graph, by Theorem \ref{thm2}, the code $C_2(H)$ is a $[|E|, |V|-1, \lambda(G(\mathbb{Z}_{p_1}))]_2$ code. By Theorem \ref{thm3}, we have  $|E|=\dfrac{(p_1-1)^2}{2}$ and $|V|=p_1.$ By Theorem \ref{thm4}, the edge-connective of the unit graph $G(\mathbb{Z}_{p_1})$ is $p_1-2$ and hence $C_2(H)$ is a $\Big[\dfrac{(p_1-1)^2}{2}, p_1-1, p_1-2\Big]_2$ code.
\end{proof}
\begin{example}
	The unit graph $G(\mathbb{Z}_{5})=(V, E)$ with $|V|=5$ and $|E|=\frac{(5-1)^2}{2}=8$ is
	\begin{center}
		\includegraphics{5.pdf}
	\end{center}
		Then the incidence matrix of the unit graph $G(\mathbb{Z}_{5})$ is
	$$H=\begin{pmatrix}
		1&1&1&1&\vline&0&0&0&0\\
		\hline
		1&0&0&0&\vline&1&1&0&0\\
		0&1&0&0&\vline&1&0&1&0\\
		0&0&1&0&\vline&0&1&0&1\\
		0&0&0&1&\vline&0&0&1&1
	\end{pmatrix}_{5\times 8}$$
Any four rows of $H$ are linearly independent and the graph is connected, then the dimension of the code $C_2(H)$ over the finite field $\mathbb{F}_2$ generated by $H$ is $4.$ The minimum distance of the code $C_2(H)$ is $3.$ Hence $C_2(H)$ is an $[n, k, d]_2=[8, 4, 3]_2$ code.
\end{example}

\section{Linear Codes from the Incidence Matrix of a unit Graph $G(\mathbb{Z}_{2p_1})$}
In this section, we  study the linear codes obtained from the incidence matrix of the unit graph $G(\mathbb{Z}_{2p_1})$ where $p_1$ is an odd prime over $\mathbb{F}_2$ and we find the parameters of the code.

The units of $\mathbb{Z}_{2p_1}$  is $U(\mathbb{Z}_{2p_1})=\{k\in \mathbb{Z}_{2p_1} \mid (k, 2p_1)=1\}$ and $|U(\mathbb{Z}_{2p_1})|=p_1-1.$

\begin{theorem}\label{thm4}
	Let $G(\mathbb{Z}_{2p_1})=(V, E)$ be a unit graph. Then 	$|V|=2p_1$ and $|E|=p_1(p_1-1).$
\end{theorem}
\begin{proof}
	Let $G(\mathbb{Z}_{p_1})$ be a unit graph. Then by definition, $V=\mathbb{Z}_{2p_1}.$
	Since $(2, 2p_1)=2\neq 1,$ by Theorem \ref{thm1}, $G(\mathbb{Z}_{2p_1})$ is a $|U(\mathbb{Z}_{2p_1})|=(p_1-1)$-regular graph. That is, degree of every vertex is $p_1-1.$
	Since the number of edges of the $k$-regular graph with $n$ vertices is $\dfrac{nk}{2},$ implies $|E|=p_1(p_1-1).$
\end{proof}

\begin{theorem}\label{thm6}
	Let $p_1$ be an odd  prime. Then the code generated by the incidence matrix $H$ of the unit graph $G(\mathbb{Z}_{2p_1})$ is a  $C_2(H)=\Big[p_1(p_1-1), 2p_1-1, p_1-1\Big]_2$ code over the finite field $\mathbb{F}_2.$
\end{theorem}
\begin{proof}
	Let  $G(\mathbb{Z}_{2p_1})$ be a unit graph and let $H$ be the incidence matrix of $G(\mathbb{Z}_{2p_1}).$ Since  $G(\mathbb{Z}_{2p_1})$ is a connected graph, by Theorem \ref{thm2}, the code $C_2(H)$ is a $[|E|, |V|-1, \lambda(G(\mathbb{Z}_{2p_1}))]_2$ code. Since  $G(\mathbb{Z}_{2p_1})$ is $(p_1-1)$-regular graph, the edge-connectivity $\lambda(G(\mathbb{Z}_{2p_1}))$  of the unit graph $G(\mathbb{Z}_{2p_1})$ is $p_1-1$ and hence the minimum distance of the code generated by $H$ is $p_1 -1.$  Since  $|E|=p_1(p_1-1)$ and $|V|=2p_1,$ the main parameters of the code $C_2(H)$ is  $[p_1(p_1-1), 2p_1-1, p_1-1]_2.$
\end{proof}
\begin{example}
	The unit graph $G(\mathbb{Z}_{6})=(V, E)$ with $|V|=6$ and $|E|=6$ is
	\begin{center}
		\includegraphics{6.pdf}
	\end{center}
	Then the incidence matrix of the unit graph $G(\mathbb{Z}_{6})$ is
	$$H=\begin{pmatrix}
		1&0&0&0&0&0&1\\
		1&1&0&0&0&0&0\\
		0&0&1&1&0&0&0\\
		0&0&0&1&1&0&0\\
		0&1&0&0&0&1&0\\
		0&0&1&0&0&0&1
	\end{pmatrix}_{6\times 6}$$
	Any five rows of $H$ are linearly independent and the graph is connected, then the dimension of the code $C_2(H)$ over the finite field $\mathbb{F}_2$ generated by $H$ is $5.$ The minimum distance of the code $C_2(H)$ is $2.$ Hence $C_2(H)$ is an $[n, k, d]_2=[6, 5, 2]_2$ code. Since $d=n-k+1,$ the code $C_2(H)$ is a MDS code.
\end{example}

\section*{Conclusion}
In this paper, we studied the codes generated by the incidence matrix of the unit graph of the different commutative rings with unity. Also we found the main parameters of the code over finite field $\mathbb{F}_2.$ We have consider only simple and undirected graphs in this article. Finding the covering radius of these codes is the further direction to work.


\section{Efficient 2D Pose Estimation and Lifting}
\label{sec:cnns}
This section describes the CNN architectures used for accurate bottom-up 2D pose
estimation and our proposed  method for 2D-to-3D body joint lifting and
for handling  miss-detections due to (self-)occlusion or failures.

\subsection{CNN-based 2D Pose Estimation}
We follow recent breakthroughs in multi-person 2D pose estimation
that use a CNN to predict confidence maps $\rho(\cdot)$ for the location of the body
landmarks in the image and part affinity fields $\phi(\cdot)$ for the location and orientation of the limbs~\cite{CPMPaf}.
We analyze different CNN architectures and the impact of their 2D
estimates on the quality of the 3D pose.

Three architectures are considered. The two firsts are the efficient pose machines
based on residual  modules (RPM) and the one based on MobileNets (MPM) introduced in~\cite{Martinez_TCSVT_2019}.
These are lightweight CNNs that refine predictions with a series
of prediction stages and are designed for efficient 2D pose estimation
with real-time performance, see Fig.\ref{fig:2d-CNN}~(a).
Additionally, we consider the Hourglass network architecture~\cite{HourGlass}
which was originally proposed for single person pose estimation.
It comprises a series of UNet-like networks that process image features at
different semantic levels.
We follow the original design  but adapt the output
to predict part affinity fields to match our multi-person scenario
by branching a duplicate of the confidence maps prediction layers
(Fig.\ref{fig:2d-CNN} (b)).

\begin{figure}[t]
\centering
\includegraphics[width=\linewidth]{architectures_2-0-eps-converted-to.pdf}
\vspace{-0.2cm}
\caption{CNN architectures used for  2D pose estimation.
  %

  (a)~Pose Machine architecture implemented by RPM and MPM~\cite{Martinez_TCSVT_2019}.
  %
  (b)~Our extension of the Hourglass network for multi-person 2D pose estimation.
}
\label{fig:2d-CNN}
\vspace{-2mm}
\end{figure}

\subsection{Pose lifting}
\label{sec:lifting}
Given 2D landmark detections, we use their corresponding depth values $Z$
to lift them according to $\bar{x} = Z \cdot K \cdot (x_{img}, y_{img}, 1)^{\top}$,
where $K= diag( 1/f_x, 1/f_y,1)$ is the depth camera matrix.
However, different errors can arise.
For example, a 2D detection might have missing depth value due to sensing failures.
Additionally, as is common in typical HRI scenarios, self and between-person
occlusions will naturally result in missing body detections.

In this paper, in these cases, rather than feeding our regressor with dummy values
which might bias estimations, we propose a simple  recovery method.
First, in case of missing depth values,
we use the mean depth of the points with valid depth information
in the landmark's vicinity.
Second, in case of missed landmark detections, we rely on a 3D pose prior to
infer their expected coordinates.
However, rather than relying on expensive-to-compute prior~\cite{Sigal:IJCV:11}, we follow a
simpler 3D limb prior based on  pairwise relationships between limb vectors.
Following a tree of limbs from the skeleton and taking the spine limb as root (see Fig.~\ref{fig:calibrated-skeleton}(a)),
we consider adjacent limbs, encode their 3D direction and length within a joint Gaussian distribution
$p(l_i, l_{\mathbf{pa}(l_i)})$, and learn the model parameters from training data.
Then, to predict the lifted coordinates $\bar{x_i}$ of a missed landmark,
we consider its associated limb $l_i$ in the skeleton whose other landmark
is already lifted, and compute  the mean of the
conditional Gaussian distribution  $p(l_i|l_{\mathbf{pa}(l_i)})$ of $l_i$
conditioned on its limb parent $\mathbf{pa}(l_i)$ to further compute $\bar{x_i}$.

Note that our approach requires some body landmarks to be detected.
Indeed,  as in our opinion it is unrealistically to attempt determining
the complete 3D pose of the person from a few detected body landmarks, {e.g}\onedot} \def\Eg{{E.g}\onedot the arm,
we assume that at least the spine limb and other two
body landmarks in the trunk (shoulders, heaps) are detected.

\section{Human 3D Pose Estimation}
\label{sec:3dpose}
\begin{figure}[t]
\includegraphics[width=\linewidth]{pose_calibration_and_error-eps-converted-to.pdf}

\vspace{-0.2cm}
\caption{
  (a) Skeleton and limb pairwise relationships;
  %
  (b) illustration of the error introduced by the lifting process of the 2D detected landmarks;
  %
  (c) mean absolute error on each coordinate when using
  the 3D lifted points as the 3D estimation on the ITOP dataset.
}
\vspace{-2mm}
\label{fig:calibrated-skeleton}
\end{figure}
This section presents our residual-pose learning approach to predict (in a camera coordinate frame )
the 3D coordinates of a human  skeleton comprising $\NJoints$ body landmarks.

\subsection{Residual Pose Learning}
Provided the 2D body landmark detections, our lifting step provides
a \emph{rough} estimate of the  3D pose.
Yet, lifted values will exhibit 3D pose estimation errors,
specially since lifted 3D points lie on the depth surface rather
than represent the inner joint (see Fig.~\ref{fig:calibrated-skeleton}).
In this regard, in absence of other sources of errors (missed detections,
occlusion, etc.) we can argue that such estimates differ only
from the true 3D pose by some coordinate offset.
This inspired us to follow a simple yet effective approach to obtain
refined estimates from rough lifted estimates.

Our approach can be set as follows: given a \emph{rough} 3D pose
estimate $\liftedpose\in \mathbb{R}^{\NJoints\times 3}$
obtained from the 2D landmark detection lifting step,
and its \emph{true} corresponding 3D pose $\truepose\in\mathbb{R}^{\NJoints\times 3}$,
the neural regressor $f$ can focus on modelling their residual $\truepose-\liftedpose$ as:
\begin{equation}
f(\liftedpose)+\liftedpose=\truepose.
\label{eq:residualpose}
\end{equation}
The function $f(\liftedpose)$ is the residual to be learned.
Graphically, these residuals represent the vector of coordinate offsets that are necessary
to predict the true 3D pose $\truepose$ (hence a residual pose).
Architecturally speaking, the operation $f(\liftedpose)+\liftedpose$
is performed by a shortcut connection with the identity mapping of $\liftedpose$,
as shown in  Fig.~\ref{fig:residual_pose}.

Additionally, we can augment $\liftedpose$ by incorporating the confidence
of the 2D detections provided by the 2D pose estimation CNN.
This will add an extra dimension for each detected landmark
$\liftedpose\in\mathbb{R}^{J\times 4}$.
In such case the shortcut connection works as a pooling layer that removes  the
extra dimension to match the one of $\truepose$.
We analyze this particular case in Section~\ref{sec:experiments}.

\begin{figure}[t]
\centering
\includegraphics[width=\linewidth]{residual_pose-eps-converted-to.pdf}
\vspace{-0.2cm}
\caption{
  Residual pose learning framework.
  %
  Our neural network regressor receives as input a lifted 3D pose $\liftedpose$.
  %
  Due to the global skip connection, the regressor has to
  predict the residual pose $f(\liftedpose)$ to be added to $\liftedpose$
  to predict the true 3D pose.
  %

  %
  The building block of our neural network regressor is a linear layer
  followed by batch normalization, ReLU activations and dropout, and
  with a skip connection.
}
\vspace{-2mm}
\label{fig:residual_pose}
\end{figure}

\subsection{Neural Network Regressor}
\label{sec:regarch}
We aim to find a simple and efficient network architecture~$f$ that
performs well enough in the regression task.
Fig.\ref{fig:residual_pose} shows a diagram with the basic building
blocks of our architecture.
It is a multi-layer network consisting on a series of fully-connected
layers, each followed by batch normalization, ReLU activations and
dropout layers.
The first layer receives as input the lifted pose $\liftedpose$
and outputs $\NFeatures$ features.
This number of features are kept fixed until the output layer that
generates the residual pose vector in $\mathbb{R}^{\NJoints\times 3}$.
Each of the inner layers have skip connections.
One can normally squeeze as many inner layers $\RegStages$ to make the
regressor deeper.
However, in this paper we set $\RegStages=3$.

\subsection{Pose Learning Loss}
Let $\totalpose=f(\liftedpose)+\liftedpose$ be the 3D pose prediction.
We use the following loss to train our neural network regressor
\begin{equation}
L_{res} = \frac{1}{\NJoints}\sum^{\NJoints}_{i=1} || \totalpose_i - \truepose_i ||_1,
\label{eq:jointloss}
\end{equation}
where $\truepose_i$ is the ground truth of the body landmark $i$ and
$ \totalpose_i$ is the 3D prediction for such landmark.
In our experiments we use the smooth L1 norm as we found out that it
works better than the L2 or plain L1 norms.

\section*{ACKNOWLEDGMENTS}

\mypartitle{Acknowledgments}:
This work was supported by the European Union under the EU Horizon 2020 Research
and Innovation Action MuMMER (MultiModal Mall Entertainment Robot), project ID
688146, as well as the Mexican National Council for Science and Technology~(CONACYT)
under the PhD scholarships~program.

\section{Conclusions} \label{sec:conclusion}

In this paper we addressed the problem of 3D pose estimation from depth images.
We decoupled 2D and 3D pose estimation and predict the 3D pose from lifted
2D detections.
We proposed a residual-pose regression learning to predict the 3D pose
by refining lifted detections.
We introduced a pairwise 3D limb prior to recover from 2D detection failures
and analyse the incorporation of 2D detection confidence in our pipeline.
Despite the simplicity of our approach we achieve competitive results in two
public datasets for single and multi-person pose estimation.
Our method propose
a more efficient alternative for multi-party HRI settings.

Our study opens the way for new research.
One limitation of our model is that it does not consider the skeleton
kinematics in the learning process.
Additionally, body motion modelling can be introduced to introduce
temporal consistency in our 3D predictions.
\section{Introduction}\label{sec:intro}
3D human pose estimation is an essential part of many applications involving
human behavior analysis, like  3D scene understanding, social robotics,
visual surveillance and gaming.
For instance,  in social HRI, the ability to sense the  3D pose
of humans provides to the robot the means to further recognize their activity
or evaluate their interaction engagement.
However, although  3D pose estimation has been a very important  topic of research,
factors like person self occlusions, pose variations, sensing conditions and
low computational budget
increase the challenge of deploying accurate, reliable and efficient
3D pose estimation systems.

\mypartitle{State-of-the-art.}
Early approaches on 3D human pose estimation
detect body landmarks in the image that are then coupled with 3D human
pose priors that account for body kinematics and physical
constraints~\cite{Sigal:IJCV:11,Ramakrishna_ECCV_2012}.
Nowadays, Deep Neural Networks (DNN) have become the mainstream approach,
which lead to the emergence of a large number of methods to address 3D  pose estimation from
color~\cite{Chen_2017_CVPR,martinez_2017_3dbaseline,Tome_CVPR_2017}
and depth images~\cite{DepthInvariant,pavlakos2017volumetric,Moon_2018_CVPR_V2V-PoseNet,ShottonPAMI,Taylor_Viturvian_CVPR_2018}.

From a methodological perspective, methods can nevertheless be grouped into two main
threads: fitting and learning methods.
The former ones extend earlier works, but rather use CNNs to localize 2D body parts,
and then fit a 3D body pose model along with constraints  via an
optimization objective defined in the image domain~\cite{KEEPSMPL_ECCV_2016}
or in the  2D-to-3D joint space~\cite{Li_ICCV_2015,Chen_2017_CVPR,Tome_CVPR_2017}.
Learning based methods take advantage of the recent DNNs
to directly predict and regress the 3D locations of the body parts with
fully connected networks~\cite{Habibie_CVPR_2019,martinez_2017_3dbaseline}.
Although simple, the 3D coordinate regression has proved to be an effective and
efficient  solution.
Moreover, information about pose kinematics can  been incorporated
as an additional limb loss~\cite{Sun2017CompositionalHP},
using a structured prediction layer~\cite{Aksan_2019_ICCV},
or via a re-projection regularizer~\cite{kanazawaHMR_CVPR_18,Wandt_CVPR_2019}.
However, a drawback of these models is that they predict the 3D coordinates
with respect to a root joint that is assumed to be known in advance, or which in
practice needs to be predicted as well.

\begin{figure}[t]
\centering
\includegraphics[width=0.75\linewidth]{pipeline_3-eps-converted-to.pdf}

	\vspace{-0.1cm}
\caption{Proposed decoupled residual pose  approach:
		%
		a) bottom-up multi-person 2D pose detection;
		%
		b) for each detected person, 2D body joints are lifted to the 3D space.
		%
		c) 3D pose estimation using a residual pose regression network.
}
\label{fig:pipeline}
\end{figure}

The depth data modality has also been largely exploited, since compared to
color images it is texture and color invariant,
and helps to remove the ambiguities in scale and shape
by providing direct access to 3D information.
As with color, some methods tend to rely on fitting approaches,
for instance by  identifying one-to-one relationships between cloud points and a
3D mesh via Iterative Closest Point (ICP)~\cite{ICP_ICCV_2011}
or using random forest~\cite{Taylor_Viturvian_CVPR_2018}.
Other approaches model the 3D location distributions of
3D points with respect to their parents in a kinematic tree~\cite{Jung_RTW_CVPR_2015,ShottonPAMI}.
As a typical example, the seminal work of Shotton~\cite{ShottonPAMI}
employed random forests to classify depth pixels into  different human joints
and used weighted voting to estimate their  3D locations.
Deep learning also improved upon these works~\cite{DepthInvariant,Moon_2018_CVPR_V2V-PoseNet,guo2017towards,DepthMultitask}.
In \cite{DepthInvariant}, multi-view human pose estimation is solved by learning a view invariant feature space
and iteratively refining  the 3D coordinates with a Recurrent Neural Network (RNN).
In~\cite{Moon_2018_CVPR_V2V-PoseNet} depth images are transformed into a voxelized representation
and 3D Gaussian likelihoods are predicted for each body joint per voxel using a costly 3D-CNN.
However,  these methods usually work on image crops centered around the person.
As a consequence, to handle the multi-person case,
a person detector is still needed as a first step,
followed by multiple forward passes of a relatively heavy image processing network to estimate
the 3D pose of each detected person, leading to an  increased computational cost.

\mypartitle{Approach and contributions.}
An overview of our approach for accurate and fast multi-person 3D pose estimation
is presented in Fig.~\ref{fig:pipeline}.
Our main idea is to better exploit the depth information and decouple the task in
two main steps:
2D multi-person pose estimation and 3D pose regression.
The motivations are that the first step can benefit from
recent accurate and efficient architectures to achieve this task,
and that the second one can be done efficiently by directly
regressing the 3D pose coordinates from the 2D ones in
two substeps: a simple but effective scheme which lifts the 2D estimates
to 3D using the depth information and pose priors (to handle partial occlusion);
and a novel efficient residual pose 3D regression methods  that works on this set of points.

This makes our approach computationally ligher for multi-person HRI settings
since compared to CNNs applied to image crops for 3D pose prediction,
the cost of our 3D regression scheme is much smaller, and the cost saving is proportional
to the number of person in the scene.
In this context, our contributions can be summarized as:
\begin{compactitem}
  %
\item we investigate an innovative method decoupling the 3D pose estimation task
  into an accurate and efficient CNN-based 2D bottom-up
  multi-person pose estimation method and 3D pose regression;
  %

  %

\item we propose a simple 2D-to-3D lifting scheme which handles 2D body joint miss detections;

  %

\item we introduce  a novel method for 3D pose regression
  from lifted 2D estimates by relying on a residual-pose
  deep-learning architecture;
  %

  %
\item we demonstrate that despite its simplicity,
  our approach achieves very competitive results
  on different public datasets and is suitable for multi-party HRI scenarios.

  %
\end{compactitem}
Models and code will be made  publicly
available\footnote{\color{red}\url{https://github.com/idiap/residual_pose}}.
The paper is organized as follows: Section~\ref{sec:cnns}
introduces our strategy for 2D pose estimation and lifting.
Section~\ref{sec:3dpose} presents our approach and regressor neural network architecture
for residual pose learning.
Experiments are described in Section~\ref{sec:experiments},
and Section~\ref{sec:conclusion} presents our conclusions.

\section{Experiments}
\label{sec:experiments}
We conducted several experiments to evaluate our approach  effectiveness
in single and multi-person scenarios.

\subsection{Depth-image datasets}

\mypartitle{ITOP~\cite{DepthInvariant}}.
This  dataset consists of images in a single person pose estimation setting.
It has 18k and 5k depth images for training and testing, respectively, recorded
with an Asus Xtion camera.
It  was built from 20 subjects performing 15 different actions each.

\mypartitle{CMU-Panoptic~\cite{Joo_2017_TPAMI}}.
It comprises multiple recordings acquired with different sensor devices such as
color and depth cameras (Kinect2).
We consider a subset of the depth recordings from the \emph{Haggling} category.
The setup contains several interacting people with diverse body pose configurations
with respect to the camera and between-person interactions.
For training we selected 15k 3D person instances from the sequence
170407\underline{ }haggling\underline{ }a3 for training.
For testing 1.5k 3D person instances were selected
from the sequence 170407\underline{ }haggling\underline{ }b3.

\subsection{Evaluation metrics}
\mypartitle{Mean average precision (mAP)}.
As standard practice in 3D human pose estimation, we use mean average precision at 10 cm (mAP{@}10cm) to measure the
3D detection performance.
A successful detection is considered when the detected 3D body landmark
falls within a distance less than 10~cm from the ground truth.
We report the average precision~(AP) for individual body landmarks
and to measure the overall performance,
the mean average precision~(mAP) defined as the mean of the APs of all body landmarks.
Larger values are better.

\mypartitle{Mean per joint position error~(MPJPE)}.
It measures the average error in Euclidean distance between the
detected 3D body landmarks and the ground truth.
Lower values are better.
We report MPJPE in  centimeters for each body landmark
and their mean for the overall performance~(mMPJPE).

\mypartitle{Percentage of correct keypoints~(PCK)}.
We use PCK to evaluate the performance of the 2D pose estimation task.
It relies in the precision and recall that result from the percentage
of correct detected keypoints (body landmarks).
We follow the evaluation protocol presented in~\cite{Martinez_IROS_2018}.
For each joint ({e.g}\onedot} \def\Eg{{E.g}\onedot knee), true positives, false positives, and false negatives
are counted using a radius obtained according to the height of the bounding box
(ground truth) containing the person.
Then, the precision and recall rates are calculated by averaging the above
values over a set of varying radius, body landmarks, and dataset samples.

\subsection{Implementation details}

\mypartitle{Image pre-processing}.
We normalize the depth images by linearly scaling the depth sensor values in
$[0,8]$~meter range into the  $[-0.5,0.5]$ range.

\mypartitle{2D CNN architectures and training}.
We keep the performance-efficiency trade-off reported in \cite{Martinez_TCSVT_2019} and
experiment with RPM with 2 stages and MPM with 4 stages.
We configure the Hourglass architecture (HG) to 2 stages as it was shown
that performance saturates at this point~\cite{HourGlass}.

We train the 2D pose estimation CNNs using Adam.
To avoid overfiting due to the low number of depth images in the addressed
datasets, and increase the 2D pose performance, we train the networks for
13 epochs with the large synthetic people dataset introduced
in~\cite{Martinez_IROS_2018}.
Then, the CNNs are finetuned using the real dataset (ITOP or CMU-Panoptic)
for 100 epochs.

\mypartitle{Residual Pose Regressor}.
We train our neural network regressor for 200 epochs using Adam and
minibatches of size 128.
We apply standard normalization to the 3D lifted pose and the 3D ground truth
pose by substracting the mean and dividing by the standard deviation.
We select $1e-3$ as initial learning rate and decrease it by 2 every 20 epochs.

\subsection{Experimental results}

\begin{figure}
\includegraphics[width=0.49\linewidth]{recall_precision.pdf}
\includegraphics[width=0.49\linewidth]{mpjpe_cnns.pdf}

\vspace{-0.2cm}
\caption{Performance of the different CNNs for 2D pose estimation.
  %
  Left: 2D pose estimation performance measured with recall and precision curves.
  %
  Right: resulting 3D estimation pose performance in terms of MPJPE for each
  body part. The lower the better.
}
\vspace{-2mm}
\label{fig:curves}
\end{figure}

\begin{figure*}[t]
\centering
\includegraphics[width=\textwidth]{examples_3dpose_2-eps-converted-to.pdf}
	\vspace{-0.2cm}
	\caption{3D pose estimation examples and their 2D projection
	of our approach on the single person ITOP dataset
	(top row) and the multi-person CMU-Panoptic dataset (bottom row).
	}

\label{fig:3d-skeletons}
\end{figure*}

\begin{table}
\centering

\begin{tabular}{l | c c c }
\Xhline{2\arrayrulewidth}

CNN model   & MPM    & RPM    & HG   \\

\Xhline{2\arrayrulewidth}

FPS         & 84    &  35     & 18  \\

\# Params   & 304.9K &  2.84M & 12.9M\\

\Xhline{2\arrayrulewidth}

F-Score (2D)& 0.96   &  0.96  & \textbf{0.97} \\ \hline

mAP{@}10cm & 85.61  &  85.96  & \textbf{85.97}\\

mMPJPE      & 6.83   & 7.18    & \textbf{6.78}  \\
\Xhline{2\arrayrulewidth}

\end{tabular}

\caption{2D and 3D pose estimation performance obtained
		for the different 2D CNN architectures and their
		computational complexity.
}
\label{tab:efficiency}
\end{table}

\begin{table*}[t]
\centering

\begin{tabular}{l | c c c c c c c c }

\Xhline{2\arrayrulewidth}

         & \multicolumn{8}{c}{ITOP (front-view)} \\

\Xhline{2\arrayrulewidth}

	& \multicolumn{8}{c}{AP{@}10cm} \\

\Xhline{2\arrayrulewidth}

Body part & \cite{Jung_RTW_CVPR_2015} & \cite{DepthInvariant} & \cite{guo2017towards} & \cite{Moon_2018_CVPR_V2V-PoseNet} & \textbf{\ResPose} &\textbf{\ResPose}$^{*}$& \textbf{\ResPose}$^{-}$& \CoordReg     \\

\Xhline{2\arrayrulewidth}
\rowcolor{Gray}
Head      & 97.8 & 98.1 & \textbf{98.7} & \underline{98.29} & 98.27  &  98.13  & 98.33 & 97.8  \\
Neck      & 95.8 & 97.5 & \textbf{99.4} & \underline{99.07} & 98.6   &  98.56  & 98.5 & 98.66  \\
\rowcolor{Gray}
Shoulders & 94.1 & \underline{96.5} & 96.1 & \textbf{97.18} & 95.34  &  95.2   & 92.78 & 95.64 \\
Elbows    & \underline{77.9} & 73.3 & 74.7 & \textbf{80.42} & 76.52  &  75.89  & 74.38 & 74.24 \\
\rowcolor{Gray}
Hands     & \textbf{70.5} & \underline{68.7} & 55.2 & 67.26 & 61.69  &  61.28  & 59.98 & 55.01 \\
Torso     & 93.8 & 85.6 & \underline{98.7} & \textbf{98.73} & 98.56  &  98.64  & 98.62 & 97.57 \\
\rowcolor{Gray}
Hips      & 80.3 & 72   & \underline{91.8} & \textbf{93.23} & 90.07  &  90.31  & 89.4 & 87.09  \\
Knees     & 68.8 & 69   & 89   & \textbf{91.80} & \underline{89.13}  &  88.93  & 88.82 & 88.29 \\
\rowcolor{Gray}
Feet      & 68.4 & 60.8 & 81.1 & \textbf{87.6}  & \underline{84.28}  &  83.52  & 83.66 & 83.99 \\

\Xhline{2\arrayrulewidth}

Mean (mAP)     & 80.5 & 77.4 & 84.9 & \textbf{88.74} & \underline{85.97}  &  85.71  & 84.9 & 84.17  \\

\Xhline{2\arrayrulewidth}

          & \multicolumn{8}{c}{CMU-Panoptic} \\

\Xhline{2\arrayrulewidth}

	     & \multicolumn{4}{c|}{MPJPE (cm)} &  \multicolumn{4}{c}{AP{@}10cm}  \\

\Xhline{2\arrayrulewidth}

Body part & \textbf{\ResPose} & \textbf{\ResPose}$^{*}$ & \textbf{\ResPose}$^{-}$ & \multicolumn{1}{c|}{\CoordReg} &  \textbf{\ResPose} & \textbf{\ResPose}$^{*}$ & \textbf{\ResPose}$^{-}$ & \CoordReg \\

\Xhline{2\arrayrulewidth}
\rowcolor{Gray}
Head      & \textbf{6.59}  & \underline{6.78}  &  10.17  & \multicolumn{1}{c|}{11.17} & \underline{96.4} &  \textbf{96.67}  & 79.47  & 72.33 \\
Neck      & \textbf{7.29}  & \underline{7.45}  &  8.5    & \multicolumn{1}{c|}{11.68} & \textbf{96.53}  & \underline{96.2} & 92.13  & 74.07 \\
\rowcolor{Gray}
Shoulders & \textbf{8.55}  &  \underline{8.66} &  10.96  & \multicolumn{1}{c|}{14.38} & \textbf{87.17}  & \underline{85.6}   & 77.17  &54.33  \\
Elbows    & \underline{14.52} &  \textbf{14.19} &  23.86  & \multicolumn{1}{c|}{20.2}  & \underline{59.17}  &  \textbf{61.97}  & 38.3   & 28.93 \\
\rowcolor{Gray}
Hands     & \textbf{27.85} &  \underline{27.96} &  31.16  & \multicolumn{1}{c|}{26.37} & 16.63  & \underline{17.47}& \textbf{17.77}  & 6.37  \\
Torso     & \underline{9.06}  &  \textbf{8.51}  &  9.92   & \multicolumn{1}{c|}{11.93} & \textbf{93.27}  & \underline{92.67} & 87.6   &67.53  \\
\rowcolor{Gray}
Hips      & \textbf{8.57}  &  \underline{8.67}  &  12.16  & \multicolumn{1}{c|}{12.99} & \textbf{91.97}  & \underline{90.27} & 70.1   &66.1   \\
Knees     & \textbf{9.24}  &  \underline{9.43}  &  14.72  & \multicolumn{1}{c|}{13.96} & \textbf{81.8}   & \underline{80.6} & 58.67  &52.33  \\
\rowcolor{Gray}
Feet      & \underline{11.26} &  \textbf{11.19} &  18.8   & \multicolumn{1}{c|}{15.54} & \textbf{70.77}  & \underline{70.5}   & 52.17  & 48.27 \\

\Xhline{2\arrayrulewidth}

Mean      & \textbf{12.2}  &  \textbf{12.2}  &  16.79  & \multicolumn{1}{c|}{16.11} & \textbf{73.41}  & \underline{73.22} & 59.17  & 48.44 \\

\Xhline{2\arrayrulewidth}

\end{tabular}

\caption{3D pose estimation performance.
		 %
		 Top: mAP of the state-of-the-art on single person pose estimation setting in the
		 ITOP dataset,
		 %
		 Bottom: mAP and mMPJPE for the multi-person pose estimation setting in the
		 CMU-Panoptic dataset.
}
\label{tab:soa-map}
\end{table*}

\mypartitle{2D Pose Network Architectures}.
We evaluate the quality of the 2D pose predictions for the 3D pose estimation
task in the ITOP dataset.
Fig.~\ref{fig:curves} shows the 2D pose estimation performance curves and the
3D pose error in terms of MPJPE for the different CNNs.
Table~\ref{tab:efficiency} summarizes these results with the maximum F-Score
obtained for 2D pose estimation, and the mAP and mMPJPE for
3D pose prediction.
Indeed, providing better 2D pose estimates reflects directly in the 3D
performance.
Overall the HG 2D detections provide the best 3D estimates
achieving the lowest mMPJPE and better mAP.
We select the HG network for the rest of the analysis.

\mypartitle{Computational Requirements}.
Table~\ref{tab:efficiency} reports the number of parameters of each CNN
and the frames per second (FPS) required for the forward pass
in a single Nvidia card GTX 1050.
Note that the FPS is also valid for the multi-person case since the CNNs
predict the pose for each individual in the image in a single forward pass.
Additionally, the neural network regressor requires 12.7M parameters but runs
at 1700 FPS, so its cost, even when applied for multiple person,  is negligible
compared to that of a 2D pose CNN.
Hence our proposed approach can run very
efficiently in real-time in a single GPU.

\mypartitle{Comparison with the state of the art}.
Table~\ref{tab:soa-map} compares the detailed AP scores for each body landmark
of our proposed approach (\textbf{\ResPose}) with the state of the art in
the ITOP dataset.
Overall our residual pose learning approach shows very competitive results
obtaining the second best performance.
The best performing work is~\cite{Moon_2018_CVPR_V2V-PoseNet} that processes
voxelized representations of the 3D space processed with a 3D CNN,
and uses an ensemble of 10 models for the final prediction.
Contrary, our residual pose approach is simpler and efficient.
Example results are shown in Fig.~\ref{fig:3d-skeletons} (top row).

\mypartitle{Multi-person 3D pose estimation}.
Table~\ref{tab:soa-map} reports AP and MPJPE for the
multi-person setting in the CMU-Panoptic dataset.
Naturally the ranges of pose profile, multiple scales, and the quality of sensing
make this setup more challenging than the single person pose setting.
The more affected body landmarks are the hands and elbows with lower AP and
larger MPJPE.
Note these are the elements that are in constant motion and are more affected
by self occlusions,
compared to other elements like the torso and head.
Fig.~\ref{fig:3d-skeletons} shows prediction examples.

\mypartitle{Recovery from 2D Failures}.
We report the results of removing the prior recovery component introduced
in Section~\ref{sec:lifting}.
Table~\ref{tab:soa-map} shows the performance for the single and multi-person
settings (\textbf{\ResPose}$^{-}$).
The performance drops specially for the multi-person scenario.

\mypartitle{2D Landmark Detection Confidence}.
We incorporated the confidence of the 2D detections provided by the CNN
that range in $[0,1]$ in our residual pose learning setting.
When a landmark was recovered by the process of Section~\ref{sec:lifting},
we set a low confidence value of $\sigma=0.1$ to identify them from the rest.
The results are reported in Table~\ref{tab:soa-map} (\textbf{\ResPose}$^{*}$).
The mAP slightly decreases in this case.
However, in the multi-person setting some especific elements
(head, elbows, hands) have slightly better detection rate.

\mypartitle{Coordinate Regression}.
We experimented with 3D coordinate regression using the neural
network architecture introduced in Section~\ref{sec:regarch} and
predict $X,Y,Z$  coordinates of the body landmarks from
lifted 2D detections, dropping the residual pose connection.
Table~\ref{tab:soa-map} compare these results (\CoordReg) with our residual
pose approach.
The performance drops for both single and multi-person settings.
Certainly when people appear roughly in the same position, as
the case in ITOP dataset, 3D coordinate regression presents
a good alternative.
However, our residual pose approach outperforms 3D coordinate
regression in both, single and multi-person settings.


\section*{Introduction}
Multiple zeta values are number series
playing important role in a wide range of subjects,  see e.~g. \cite{B, GM, IKZ}.
Although the arithmetic nature of these numbers is still unknown,
one may try to consider all consequences of natural relations among
these series, that is the algebra of formal multiple zeta values.
One could expect that there no other rational relations
among them, that is the formal algebra of multiple zeta values
is isomorphic to the algebra rationally generated by multiple zeta values.
The long-standing conjecture states that all relations
for the both of them are regularized double shuffle relations (\cite[Conjecture 1]{IKZ}).

In the first section, we introduce an algebra of number series called big zeta values, which generalizes
multiple zeta values. One may consider the corresponding formal algebra.
The main theorem of the paper states that relations among big zeta values imply that any big zeta value
equals a rational linear combination of multiple zeta values. Thus, the algebra of (formal)
big zeta values is another hopefully more convenient form
of the (formal) algebra of multiple zeta values.

In the second section, we describe in some detail relations among big zeta values.
An interesting purely algebraic problem is to find all
relations among formal multiple zeta values inside the formal algebra of big zeta values.
Stuffle relations obviously follow from these relations.
It is harder but possible to prove shuffle relations in terms of number series,
 see  \cite{KMT}.
One may expect that there are no other relations except regularized double shuffle
relations as it is suggested by the conjecture mentioned above.

The third section is devoted to the proof of the main theorem
of the paper, which states that every big zeta value
is a rational linear combination of multiple zeta values. Surprisingly enough,
the proof does not use all relations among big zeta values.
It does not use the invariance
of a big zeta value under a permutation of columns of the matrix defining the value.
Generally,  a basic matrix may become not basic after
such a permutation, except in the case of reflection. The invariance
of a big zeta value under reflection is analogous to the duality theorem (see e.~g. \cite{Dua})
for multiple zeta values.

To prove the main theorem we build an algorithm
which allows us to turn any big zeta value
into a rational linear combination of multiple zeta values.
The algorithm is of independent interest. It may be considered
as a generalization of the stuffle product for the following reason.
The product of two multiple zeta  series is not a multiple zeta  series,
but a big zeta  series corresponding to the direct sum
of matrices of these multiple zeta values.
Applying the algorithm to it
we may expand it back in multiple zeta values. One may see that
the result is the stuffle product.

In the last section,
we connect our construction with integrals
of regular differential
forms with logarithmic singularities at
infinity on the moduli space of stable curves of genus 0 by the standard simplex
in cubical coordinates.
This construction motivated our definition of the algebra of big zeta values.
It turns out that term by term integration of the Taylor series
of such an integral is a rational linear combination of
big zeta series, and conversely, every big zeta series
is a rational linear combination of integrals of such series.

Different choices of coordinates on the integration region
result in different number series representing the same integral.
It gives a bunch of relations among the big zeta series generalizing the duality
relation mentioned above. It would be interesting to derive them
directly from basic relations.

The formal algebra of big zeta values  reminds the formal algebra of periods from  \cite{BCS, M} but is
more manageable. Relations in it mimic relations in the algebra of periods:
 Orlik--Solomon relations (\ref{osr}) correspond to the Stokes theorem
and  harmonic product relations (\ref{hpr}) correspond to  Arnold relations.
It would be interesting to find an analogous algorithm, which using integral relations expresses any integral as above as a rational combination of iterated integrals.
It would prove or disprove conjectures from \cite{BCS, M}.

Being proved that all relations among integrals as above
follows from basic relations among big zeta number series,
one may define purely algebraically the "integral" of regular differential
forms with logarithmic singularities at
infinity on the moduli space of stable curves of genus 0 by the standard simplex,
which takes values in the formal algebra of big zeta values.
This "integral" would have all properties of the usual integral
such as the Stokes theorem and the Fubini theorem.
It would follow that the  Drinfeld associator with coefficients in the formal algebra of  big zeta values defined by this "integral"
obey the standard associator relations.
Thus, the set of relations among formal multiple zeta values inside formal big zeta values algebra, which are conjecturally regularized double shuffle relations,
would imply the Drinfeld associator relations. This is the subject of future research.

{\bf Acknowledgments.} I am grateful to H.~Tsumura and W.~Zudilin for helpful discussions.

\section{Algebra of big  zeta values}

Denote by $e_{ab}$ a positive $A_w$-root,
that is an element $(r_i)\in\mathbb{Z}^w$
defined by
\begin{equation}
r_{i}=
\begin{cases} 1 & \mbox{if}\quad a\le i\le b \\
0 & \mbox{otherwise} \end{cases}
\label{root}
\end{equation}

\begin{definition}
A $d\times w$-matrix $A=(a_{ij})$ is called {\em basic} if it is of rank $d$, has no zero columns and
all its rows are positive $A_w$-roots  (\ref{root}).
\label{def}
\end{definition}
In other words, rows of a basic matrix are coordinates  of positive $A_w$-roots,
which are linearly independent and are not contained in any coordinate hyperplane.

In the following definition the term "formal" means that we do not care
about convergence of the number series.

\begin{definition}
For a basic matrix $A$ the  {\em big zeta  series} is a formal series
\begin{equation}
\mathop{\mathrm{Z}}(A)=\sum_{\substack{n_i\in \mathbb{N}, \\ 1\le i\le d}}\frac{1}{\prod_j(\sum_i a_{ij} n_i)}
\label{bzv}
\end{equation}
The width $w$ of $A$ is called the {\em weight}  and its height
$d$ is called the {\em depth} of the big zeta series.
\end{definition}

Note that the big  zeta series does not depend on the order
of rows of the matrix.
\begin{example}
\label{e1}
$$
A=
  \begin{array}{@{} c @{}}
    \left (
      \begin{array}{ *{11}{c} }
        \overmat{k_d}{1 & \dots & 1} & \overmat{k_{d-1}}{1 & \dots & 1} & \dots & \dots&\overmat{k_1}{1 &\dots &1 } \\
	0 & \dots & 0 & 1 &\dots & 1 & \dots & \dots & \dots&  \dots & 1 \\
	0 & \dots & 0 & 0 &\dots & 0 & 1 & \dots & \dots&  \dots & 1 \\
        \vdots & \vdots & \vdots & \vdots &\vdots & \vdots& \vdots & \vdots & \vdots &  \vdots & \vdots \\
        0 & \dots & 0 & 0 &\dots & 0 & \dots & 0 & 1 &  \dots & 1 \\
      \end{array}
    \right )\\
    \mathstrut
  \end{array}
$$
\medskip

\begin{equation}
\mathop{\mathrm{Z}}(A)=\sum_{\substack{n_i\in \mathbb{N}, \\ 1\le i\le d}}\frac{1}{ {n_1}^{k_d}(n_1+n_2)^{k_{d-1}}\cdots(n_1+\dots+n_d)^{k_1}}
=\mathop{\zeta}(\sq{k_1, \dots, k_d})
\label{mzv}
\end{equation}

This is the multiple zeta series, see e.~g. \cite{IKZ}.
\end{example}

\begin{example}
\bigskip
\medskip
$$A=
  \begin{array}{@{} c @{}}
    \left (
      \begin{array}{ *{9}{c} }
		 \overmat{a+b}{1 & 1 & \dots & 1& \dots& 1} & 0 & \dots  & 0\\
         0& \dots & 0 & \undermat{b+c}{1 & \dots&1 & \dots &1&1}
      \end{array}
    \right )\\
    \mathstrut
  \end{array}
$$
\medskip
\begin{equation}
\mathop{\mathrm{Z}}(A)=\sum_{n,m \in \mathbb{N}}\frac{1}{n^am^c(n+m)^b}
\end{equation}

This is the Tornheim sum, see e.~g. \cite{BZ}.
\end{example}

Space of formal big zeta series is equipped with a product
given by the direct sum of matrices.
Define the {\em formal big zeta values algebra}
as the space rationally generated
by formal big zeta  series with this product modulo natural relations.
These relations are Orlik--Solomon relations (\ref{osr})
and harmonic product relations (\ref{hpr}) below
plus the invariance of the big zeta series under permutations
of rows and columns of the defining matrix.

Example \ref{e1} above shows that formal multiple zeta values lie
in the formal big zeta values algebra. The following theorem
states that they generate the whole algebra as a vector space.

\begin{theorem}
Any formal big zeta value  is a rational linear combination of formal multiple zeta values of the same weight.
\label{main}
\end{theorem}

The proof of Theorem \ref{main} occupies Section \ref{proof} below.

\begin{cor}
The sum of a convergent big zeta series is a linear rational combination
of multiple zeta values of the same weight.
\end{cor}

The following proposition gives a convergence criterion of
a big zeta series. It also follows from results of Section \ref{integral}
below.

\begin{prop}
For a basic matrix $A$  the big zeta series $\mathop{\mathrm{Z}}(A)$
converges iff $A$ does not contain any rows with only one unit.
In other words, it converges iff the corresponding set of positive
$A_w$-roots contains no simple roots.
\label{conv}
\end{prop}
\begin{proof}
Let $A$ contains a row with only one unit.
Consider the subseries, where  all summation parameters are fixed
except the one corresponding to this row.
This subseries is proportional to the harmonic series without some initial interval.
Thus the series diverges.

The big zeta series (\ref{bzv}) is dominated by the series
$$
\sum_{\substack{n_i\in \mathbb{N}, \\ 1\le i\le d}}\frac{1}{\prod_j(\max_i a_{ij} n_i)}
$$
Divide the summation region in components corresponding to the diagonal
stratification of $\mathbb{N}^d$, that is to total orders on the set of rows.
 At each stratum this series
is equal to some multiple zeta series (\ref{mzv}), and $k_1$ of all these
series is more than the minimal number of units in rows of $A$.
If every row of $A$ contains at least two units, then $k_1$
of all these multiple zeta values are more than 1. All such multiple
zeta series are known to be convergent, see e.~g. \cite{IKZ}.
It follows that the big zeta  series converges too.
\end{proof}

\section{Relations}

For  a basic matrix $A$  and a differential operator with constant coefficients
$D\in \mathbb{Q}[\partial/\partial z_1, \dots , \partial/\partial z_d]$
introduce notation:
\begin{equation}
\mathop{\mathrm{Z}}(A, D)=\sum_{\substack{n_i\in \mathbb{N}, \\ 1\le i\le d}}\left(D\frac{1}{\prod_j(\sum_i a_{ij} z_i)}\right)(n_1, \dots, n_d)
\label{diff}
\end{equation}

\begin{prop}
For a basic $d\times w$-matrix $A$  and a homogeneous differential operator  with constant coefficients $D$, series $\mathop{\mathrm{Z}}(A, D)$
is a rational linear combination of series $\mathop{\mathrm{Z}}(A')$ for some $A'$s   of depth $d$ and of weight $w+\deg D$.
\end{prop}

\begin{proof}
If operator $D$ is of degree one, the right hand side of (\ref{diff})
is equal to a rational linear combination of series  $\mathop{\mathrm{Z}}(A')$ ,
where $A'$ are matrices equal to $A$ with one  column doubled.
Then proof proceeds by induction on $\deg D$.
\end{proof}

Let $V$ be a vector space over $\mathbb{Q}$ and $\Delta\subset V^*$ be
a finite set of non-zero linear functions on $V$, which generates the dual space $V^*$.
Denote by $G_\Delta$ the subspace of rational functions on $V$
generated by functions
$$
\frac{1}{\prod_{\alpha\in \kappa}\alpha^{n_\alpha}},
$$
where $\kappa\subset \Delta$ is a subset, which generates $V^*$.
The space $G_\Delta$ is a module over the algebra $\mathbb{Q}[V]$
 of differential operators with constant coefficients on $V$.
\begin{prop}
The $\mathbb{Q}[V]$-module $G_\Delta$ is generated by functions
$$
\frac{1}{\prod_{\alpha\in \kappa}\alpha},
$$
where $\kappa$ ranges over all subsets of $\Delta$, which are bases of $V^*$.
\label{vergne}
\end{prop}
\begin{proof}
The proof is straightforward, see \cite[Lemma 2]{BV}.
\end{proof}

\begin{prop}
For a basic $d\times w$-matrix $A$  and a  differential operator  with constant coefficients $D$,  series $\mathop{\mathrm{Z}}(A, D)$ is a rational linear combination of series $\mathop{\mathrm{Z}}(S, D\cdot D')$ for some $D'$s, where $S$ ranges over all square submatrices of $A$ of rank $d$.
\label{square}
\end{prop}

\begin{proof}
Let $\Delta$ be the set of linear functions given by columns of matrix $A$.
Restrictions imposed on this matrix by Definition \ref{def}
guarantee that conditions of Proposition \ref{vergne}
are satisfied. It follows that
$$
\frac{1}{\prod_j(\sum_i a_{ij} z_i)} = \sum_{S}D'_{S} \frac{1}{\prod_j(\sum_i s_{ij} z_i)}
$$
for some operators  with constant coefficients $D'_S$,
where the summation is taken over all square submatrices $S=(s_{ij})$ of $A$
of rank $d$. Substituting this to  definition (\ref{diff})
we get the statement.
\end{proof}

Let $l_1, \dots,  l_k$ be non-zero linear functions on a vector space $V$
such that $\sum_i\alpha_i l_i=0$ for some numbers $\alpha_i$.
The  equality
\begin{equation}
\sum_{i=1}^k \,\frac{\alpha_i}{l_1\cdots \hat{l}_i\cdots l_k}=0
\label{OS}
\end{equation}
of rational functions on $V$ is referred to as the Orlik--Solomon relation
(\cite[3.5]{OT}).

\begin{prop}[Orlik--Solomon relations]
Let $A$  be a basic matrix  and $v_1, \dots, v_k$ be a set its distinct columns such that
$\sum_i \alpha_i v_i=0$ for some  $0\neq \alpha_i\in \mathbb{Q}$.
Then for a  differential operator  with constant coefficients $D$
\begin{equation}
\sum_{i=1}^k\alpha_i\mathop{\mathrm{Z}}(\hat{A}_i, D)=0,
\label{osr}
\end{equation}
where $\hat{A}_i$ is  matrix $A$ with  column $v_i$ excluded.
\label{os}
\end{prop}

\begin{proof}
Substituting (\ref{OS}) in  definition (\ref{diff}) we get the statement.
\end{proof}

Another family of relations among big zeta series is
given by subdivision of the set of summands of the series into  groups
with  different total orders of the summation arguments.

Introduce some notations. Consider maps of polynomial algebras
\begin{equation}
\begin{gathered}
m_1\colon \mathbb{Q}[z_1, z_2]\ni p(z_1,z_2)\mapsto p(z_1+z_2, z_1)\in \mathbb{Q}[z_1, z_2]\\
m_2 \colon \mathbb{Q}[z_1, z_2]\ni p(z_1,z_2)\mapsto p(z_1+z_2, z_2)\in \mathbb{Q}[z_1, z_2]\\
m_3 \colon  \mathbb{Q}[z]\ni p(z)\mapsto p(z_1+z_2)\in \mathbb{Q}[z_1, z_2]
\end{gathered}
\label{mmm}
\end{equation}
Denote by
\begin{equation}
\begin{gathered}
m_1^*\colon \mathbb{Q}[\partial/\partial z_1, \partial/\partial z_2]\to \mathbb{Q}[\partial/\partial z_1, \partial/\partial z_2]\\
m_2^*\colon \mathbb{Q}[\partial/\partial z_1, \partial/\partial z_2]\to \mathbb{Q}[\partial/\partial z_1, \partial/\partial z_2]\\
m_3^*\colon \mathbb{Q}[\partial/\partial z_1, \partial/\partial z_2]\to \mathbb{Q}[\partial/\partial z]
\end{gathered}
\label{m&ms}
\end{equation}
the maps linear dual to (\ref{mmm}).

For a set of $w-$vectors $\{v_i\}$ and a $d\times w$-matrix $M$ denote by $[v_1; \dots; v_n; M]$
the matrix $M$ with rows $v_i$ added somewhere.
Recall that a big zeta value does not depend on the order of rows.
Thus for a matrix and vectors as above
the value $\mathop{\mathrm{Z}}([v_1; \dots; v_n; M], D)$ is well defined
for a basic matrix $[v_1; \dots; v_n; M]$,
the property of being basic also does not depend
on the order of rows.

\begin{prop}[Harmonic product relations]
Let $A$ be a basic matrix containing rows $e_{ij}$ and $e_{(j+1)k}$ (see (\ref{root})).
Denote by $A'$ the matrix $A$ with these rows excluded.
Then for a differential operator with constant coefficients $D$
acting on the space of rows of  $A$
\begin{equation}
\mathop{\mathrm{Z}}(A, D)= \mathop{\mathrm{Z}}([e_{ik};e_{ij};A'], m_1^*D)+ \mathop{\mathrm{Z}}([e_{ik};e_{(j+1)k};A'], m_2^*D)+\mathop{\mathrm{Z}}([e_{ik};A'], m_3^*D),
\label{hpr}
\end{equation}
where $m^*_i$ acts trivially on the subspace generated by rows of $A'$
and acts as in (\ref{m&ms}) on the subspace generated by the added  rows.
\label{hp}
\end{prop}

\begin{proof}
Denote by $n$ and $m$ the summation parameters corresponding
to rows $e_{ij}$ and $e_{(j+1)k}$
in $\mathop{\mathrm{Z}}([e_{ij};e_{(j+1)k};A'], D)=\mathop{\mathrm{Z}}(A,D)$.
Divide summands of this series in three groups:  ones with $n>m$,
ones with $n<m$ and ones with $n=m$.
One may see, that they correspond to three summands of the right hand side
of   (\ref{hpr}).
\end{proof}

\section{Proof of theorem \ref{main}}
\label{proof}
Denote by $T_d$ the upper triangular $d\times d$-matrix with units
on and above the diagonal.
One may see that for any homogeneous differential operators
with constant coefficients $D_i$ of degree  $i$, the expression
\begin{equation}
\sum_{i=1}^w\mathop{\mathrm{Z}}(T_i, D_{w-i})
\label{form}
\end{equation}
is a rational linear combination of multiple zeta  values (\ref{mzv}) of weight $w$.

On the other hand, by Proposition \ref{square}
any multiple zeta value of weight $w$ may be written in the form (\ref{form}).
It follows that Theorem \ref{main} may be reformulated as follows.

\begin{thmbis}{main}
Any formal big zeta value of weight $w$ is a rational linear combination of the ones of the form
$\mathop{\mathrm{Z}}(T_d, D_{w-d})$ for some  homogeneous differential operators
with constant coefficients $D_i$ of degree  $i$.
\end{thmbis}

\begin{proof}
The proof proceeds by induction on the depth.
For depth $1$ the statement is clear. Suppose that it is proved for depths
less than $d$.

By Proposition \ref{square}, every big zeta value may be written
as a rational linear combination of series of the form $\mathop{\mathrm{Z}}(S, D)$,
where $S$ is a basic $d\times d$-matrix. We need to show that all
these series are linear combinations of ones of the form $\mathop{\mathrm{Z}}(T_d, D')$ modulo
series of a lower depth.  The strategy of the proof is to  replace iteratively
$\mathop{\mathrm{Z}}(S, D)$ with linear combinations of
big zeta series whose matrices are closer and closer to $T_d$.

To define what means "closer", introduce an order on the set of indices of a $d\times d$-matrix as follows
$$
11\,<\,21\,<\,\dots\, <\,d1 \,<\,12\,<\,\dots\,<\, d2\,<\,\dots\, <\,(d-1)d\,< \,dd
$$
We say that two matrices are equal up to place $ij$ if
they have equal entries at places not bigger than $ij$.
We say that two matrices are equal exactly up to the place $ij$
if they are equal up to this place and are not equal up to the next bigger place.
A matrix $A$ is closer to matrix $C$ than matrix $B$ if $A$ is equal to $C$
exactly up to place $ij$, $B$ is equal to $C$
exactly up to place $i'j'$ and $ij>i'j'$.

The proof consists in the iterative application of two procedures.

\begin{proc}
{\em turns a basic matrix into a linear combination of upper triangular
basic matrices.}

Given a basic matrix,
rearrange its rows so that the number of the first nonzero entry in the row
decreases from the top to the bottom.
Consider rows with the first non-zero entries in the row. If there are more than one such rows,
take any two  and apply Proposition \ref{hp} to them. We get three matrices.
The one of the lower rank may be thrown away by the induction assumption.
Take two other matrices.
Rearrange their rows. The number of rows with non-zero entry at the first column
of both matrices is one less than this number of the initial matrix. Repeat this operation until the number of such rows equals one.
Then turn to  rows with the first non-zero element at the second place.
And so on.
\end{proc}

\begin{proc}
{\em turns an upper triangular basic matrix, which is equal to $T_d$
exactly up to place $(i-1)j$, into a linear combination of basic matrices, which are equal to $T_d$
up to   place $ij$.}

Given such a matrix,
let $k$ be the number of the last non-zero entry of its $j$-th row,
"$*$" means an unknown entry:
\renewcommand{\arraystretch}{1.8}
$$
\scalebox{0.75}{
\begin{blockarray}{cccccccccccccc}
 &  & & & &$i$ & & &$j$& &$k$ & & &\\
\begin{block}{c(ccccccccccccc)}
& 1 &1 & \dots& \dots&\dots &\dots&\dots &1 &* &\dots&\dots&\dots & *\\
& 0 &1 &\dots &\dots & \dots& \dots&\dots &1 &* &\dots&\dots &\dots& *\\
&\vdots &\vdots &\vdots &\vdots &\vdots &\vdots&\vdots &\vdots &\vdots &\vdots &\vdots&\vdots&\vdots \\
&0 &\dots &0 &1 &1 &\dots &\dots & 1& *&\dots &\dots&\dots& *\\
$i$&0 &\dots &0 &0 &1 &\dots &1& 0&\dots &\dots &\dots&\dots& 0\\
&0 &\dots &0 &0 &0 &1 &* & \dots&\dots &\dots &\dots&\dots& *\\
&\vdots &\vdots &\vdots &\vdots &\vdots &\vdots&\vdots &\vdots &\vdots&\vdots &\vdots & \vdots&\vdots\\
$j$&0 &\dots &\dots & \dots&\dots & \dots&0 &1 &\dots &1& 0&\dots&0\\
&\dots &\dots &\dots &\dots &\dots &\dots & \dots& \dots&\dots& \dots&\dots&\dots &\dots \\
\end{block}
\end{blockarray}
}
$$
Apply Proposition \ref{hp} to rows $i$ and $j$.
The matrix given by the last term in (\ref{hpr}) is of depth $d-1$ and may be thrown away by the induction
assumption.
Rearrange the rows of two other matrices so that the number of the first nonzero entry
in the row
decreases from the top to the bottom.

 The second term of the right hand side of (\ref{hpr})
gives a matrix, which
is equal to $T_d$ up to place $ij$.
Thus, consider the matrix given by the first  term.

Add to this matrix a column as follows, the added column is marked with an arrow:
\begin{equation}
\scalebox{0.75}{
\begin{blockarray}{cccccccccccccccc}
 &  & & & &$\mathbf\downarrow$&$i+1$ && & $j+1$&& &$k+1$ & & &\\
\begin{block}{c(ccccccccccccccc)}
& 1 &1 & \dots& \dots& 1 &\dots&\dots&\dots &1 &* &\dots&\dots&\dots&\dots & *\\
& 0 &1 &\dots & \dots& \dots&\dots& \dots&\dots &1 &* &\dots&\dots&\dots &\dots& *\\
&\vdots &\vdots &\vdots &\vdots&\vdots &\vdots &\vdots &\vdots &\vdots &\vdots &\vdots&\vdots&\vdots &\vdots&\vdots\\
&0 &\dots &0 &1&1 &1 &\dots&\dots & 1& *&\dots&\dots &\dots&\dots& *\\
$i$&0 &\dots &\dots &0 &1&1 &\dots&\dots & 1&\dots&\dots&1 &0&\dots& 0\\
$i+1$&0 &\dots &\dots &\dots&0 &1 &\dots&1 & 0&\dots&\dots &\dots &\dots&\dots& 0\\
&0 &\dots &\dots &\dots &\dots&0 &1 &\dots& *&\dots&\dots &\dots&\dots&\dots& *\\
&\vdots &\vdots &\vdots&\vdots &\vdots &\vdots&\vdots &\vdots &\vdots&\vdots &\vdots&\vdots &\vdots & \vdots&\vdots\\
$j-1$&0 &\dots &\dots &\dots &\dots& \dots&0&1 &* &\dots &\dots&\dots& \dots&\dots&*\\
$j$&0 &\dots &\dots &\dots &\dots &\dots&\dots & \dots& 0&1& *&\dots&\dots&\dots &* \\
&\dots &\dots &\dots &\dots &\dots &\dots&\dots & \dots&\dots& \dots&\dots& \dots&\dots&\dots &\dots \\
\end{block}
\end{blockarray}
}
\label{m}
\end{equation}
(it was an ambiguity in arranging $i$-th and $(i+1)$-th rows, we choose this variant).
One may see that columns with numbers from $i$ to $j+1$ in this matrix
are linearly dependent, the first one is expressing through others.
Applying Proposition \ref{os} we get that the big zeta series with the matrix (\ref{m})
with the $i$-th column removed (that is the initial matrix)
is a linear combination of matrices, which are (\ref{m})
without some column numbered from $i+1$ to $j+1$.

Consider these matrices. Let the $i'$-th column was removed.
Then in the resulting matrix rows with numbers  $i'$ and $i'+1$
has $0$s before $i'$s place and $1$ at the $i'$s place.
As in Procedure 1, apply Proposition \ref{hp} to these rows and
rearrange the rows of the matrix
so that the number of the first nonzero entry in the row decreases from top to bottom.
One may see that we get a matrix, which is equal to $T_d$ up to place $ij$.
\end{proc}

To prove the statement of the theorem, firstly apply Procedure 1 to matrix $S$. It gives an upper triangular matrix.
Find the first place, where it differs from $T_d$, denote it by $ij$. Apply Procedure 2.
The result is a matrix, which is closer to $T_d$. Apply to it Procedure
1. One may see that it
does not change the first $i$ rows and the first $j$ columns of the matrix.
It follows that it is equal to $T_d$ at least up to place $ij$.
Find the next place where it differs from $T_d$. Again apply Procedure 2 and so on.
\end{proof}

\section{Integral representation}
\label{integral}

Denote by $\M{w+3}$ the moduli space of
$w+3$ distinct points on the complex projective line
considered up to the
action of the M\"obius group.
This group acts freely  on sets of three distinct points.
To introduce coordinates on $\M{w+3}$
fix three of $w+3$ points at $0$, $1$ and $\infty$.
Coordinates $(t_1, \dots, t_w)$ of other points of the configuration
are called {\em simplicial coordinates} on $\M{w+3}$.

\begin{prop}
The algebra of regular differential
forms on $\M{w+3}$ with logarithmic singularities at
infinity
is generated by $1$-forms
\begin{equation}
\omega_{ij}=\frac{dt_i-dt_j} {t_i-t_j} \quad 0\le i< j\le w+1 \quad ij\neq 0w+1,
\label{omegas}
\end{equation}
where $t_i$ are simplicial coordinates, $t_0=1$ and $t_{w+1}=0$.
The only relations among  them
are   Arnold relations:
\begin{equation}
\omega_{ij}\wedge\omega_{jk}+\omega_{jk}\wedge\omega_{ik}+\omega_{ik}\wedge\omega_{ij}=0
\label{arr}
\end{equation}
\label{ar}
\end{prop}
\begin{proof}
See e.~g. \cite[6.1]{B}.
\end{proof}

Introduce cubical coordinates:
\begin{equation}
x_1=t_1 \quad x_2=t_2/t_1\quad \dots \quad x_w=t_w/t_{w-1}
\label{coord}
\end{equation}
The standard simplex
\begin{equation}
\Delta_{w}=\{1>t_1>\dots>t_w>0,\quad t_i\in \mathbb{R} \}
\end{equation}
in cubical coordinates turns into the cube
\begin{equation}
\square=\{0<x_i<1,\quad x_i\in \mathbb{R}\}
\label{cube}
\end{equation}

Let $A=(a_{ij})$ be a basic $d\times w$-matrix such that the series $\mathop{\mathrm{Z}}(A)$
converges, see Proposition \ref{conv}.
Integrating the Taylor series of the integrand term by term we get  a
relation
\begin{equation}
\mathop{\mathrm{Z}}(A)=\int_\square \left(\prod_{j=1}^d \frac{\prod_k x_k^{a_{jk}}}{1-\prod_k x_k^{a_{jk}}} \cdot \prod_{i=1}^{w} \frac{dx_i}{x_i}\right)
\label{int}
\end{equation}
where the integration region is given by (\ref{cube}).

\begin{prop}
The integrand of (\ref{int}) is a regular differential
form on $\M{w+3}$ with logarithmic singularities at
infinity.
\label{reg}
\end{prop}

\begin{proof}
Denote by $l(i)+1$ and $r(i)$ the numbers of the first and the
last unit in the $i$-th row of matrix $A$.
Substituting (\ref{coord}) into the integrand of (\ref{int})
we get the differential form
\begin{equation}
\prod_{j=1}^{d}\frac{t_{r(j)}}{t_{l(j)}-t_{r(j)}}\cdot\prod_{i=1}^{w}\frac{dt_i}{t_i}
\label{dform}
\end{equation}
where $t_0=1$.

Let us prove that this is a linear combination of products of differential forms (\ref{omegas}).

Denote by $G$ an unoriented graph with vertices $V(G)=\{0,1,\dots, w\}$,
with edges $E(G)$ corresponding to rows of $A$,
vertices $i-1$ and $j$ are connected by an edge iff $A$ contains  row $e_{ij}$
(see (\ref{root})). $G$ has no cycles, otherwise rows corresponding to edges
forming the cycle would be linear dependent. Thus, $G$ is a disjoint union of trees.
Make these trees rooted, choosing the vertex labeled by the minimal number
as a root.  Equip $G$ with an orientation by orienting any edge into direction "from the root"
of the tree to which this edge belongs. Call an edge "right" if its source's label
is less than its target's one, otherwise call it  "wrong".

Rewrite (\ref{dform}) as
\begin{equation}
\prod_{\substack{(ij)\in E(G), \\(ij) \mbox{ \small is right}}} \frac{t_j}{t_i-t_j}
\cdot\prod_{\substack{(ij)\in E(G), \\(ij) \mbox{ \small is wrong}}} \left( \frac{t_i}{t_i-t_j}-1\right)\cdot\prod_{i=1}^{w}\frac{dt_i}{t_i}
\end{equation}
where $(ij)$ connects vertices $i$ and $j$, $j>i$.

Expand brackets. Every vertex of $G$ has at most one incoming edge.
It follows that each of the summands has in the numerator
at most one $t_i$ for any $i$,
which cancels with the same term in the denominator.
Thus, each summand is
proportional to a product of differential forms (\ref{omegas}).
\end{proof}

We left to the reader to check that for matrix satisfying conditions
of Proposition \ref{conv} the differential form given by this proposition is convergent on
the standard simplex,  convergence conditions are described in \cite{B, BCS}.

\begin{prop}
Differential forms appearing as integrands in (\ref{int})
for any, not only convergent basic matrix $A$, generate the space of regular differential
forms of top degree on  $\M{w+3}$ with logarithmic singularities at
infinity.
\end{prop}
\begin{proof}
Let
\begin{equation}
\bigwedge_{k=1}^{w} \omega_{i(k)j(k)}=\bigwedge_{k=1}^{w} \mathop{d\, log}(t_{i(k)}-t_{j(k)})
\label{factor}
\end{equation}
 be a top degree form on  $\M{w+3}$ with logarithmic singularities at
infinity. By Proposition \ref{ar}, these forms generate
the space of top degree forms  with logarithmic singularities at
infinity. Consider the subset of indexes $k$ for which $j(k)\neq w+1$. Build a matrix with rows equal to roots
$e_{(i(k)+1)j(k)}$ for this set of indexes. These rows are linearly
independent, otherwise functions $(t_{i(k)}-t_{j(k)})$ would be
linearly dependent, what would follow vanishing of the product of their
$\mathop{d\, log}$'s, this is  a consequence of the Arnlod relations (\ref{arr}).
Thus, the matrix is basic.

Consider the differential form, which is the integrand of (\ref{int})
associated with this matrix. By the proof of Proposition \ref{reg},
these form is equal to the initial one up to some
differential forms of type (\ref{factor})  with a bigger number
factors $\omega_{i(w+1)}$. Applying  decreasing induction
on the number of such factors we get the statement.
\end{proof}

It follows that the space of convergent big zeta value of weight $w$ is the space of periods
of the pair $(\Md{w+3}, \Md{w+3}\setminus\M{w+3})$, notation $\Md{w+3}$
is introduced in \cite{B}.
Theorem \ref{main} states that this space is generated
by multiple zeta values of weight $w$.
It is in good agreement with the fact that all periods of these pair are given by periods of the fundamental group
of a projective line without three points (\cite{Del}),
that is multiple zeta values, what is proved in \cite{B}.
In contrast with this proof, which states that the weights of multiple
zeta values are not bigger than $w$, we get multiple zeta values exactly of weight  $w$.
It is natural because in (\ref{int}) we integrate a differential form
with logarithmic poles.

\bibliographystyle{alpha}