\section{Introduction}

Products are a common way in mathematics of constructing larger objects from
smaller building blocks. For graphs, hypergraphs, and related set systems
several types of products have been investigated, see
\cite{HOS-12,Hammack:2011a} for recent overviews.

In this contribution we will focus on the \emph{Cartesian product} of
\emph{directed hypergraphs} that are the common generalization of both directed
graphs and (undirected) hypergraphs. In particular, we present a fast and
conceptually very simple algorithm to find the decomposition of directed
hypergraphs into \emph{prime} hypergraphs (its so-called prime factors), where a
(hyper)graph is called prime if it cannot be presented as the product of two
nontrivial (hyper)graphs, that is, as the product of two (hyper)graphs with at
least two vertices.

\paragraph{Graphs and the Cartesian Product}
A graph is a tuple $G=(V,E)$ with non-empty set of vertices $V$ and a set of
edges $E$ containing two-element subsets of $V$. If the edges are ordered pairs,
then $G$ is called directed and undirected, otherwise. The Cartesian
\emph{graph} product was introduced by Gert Sabidussi \cite{Sabidussi:60}. As
noted by Szamko{\l}owicz \cite{Szamko:62} also Shapiro introduced a notion of
Cartesian products of graphs in \cite{Shapiro:53}. Sabidussi and independently
V.G. Vizing \cite{Vizing63:CartProd} showed that connected \emph{undirected}
graphs have a representation as the Cartesian product of prime graphs
that is unique up to the order and isomorphisms of the factors. The question
whether one can find the prime factorization of connected undirected graphs in
polynomial time was answered about two decades later by Feigenbaum et al.\
\cite{FHS:85} who presented an $O(|V|^{4.5})$ time algorithm. From then on, a
couple of factorization algorithms for undirected graphs have been developed
\cite{AHI:92,Feder:92,FHS:85,Imrich07:linear,Winkler:87}. The fastest one is due
to Imrich and Peterin \cite{Imrich07:linear} and runs in linear-time
$O(|V|+|E|)$.

For connected \emph{directed} graphs, Feigenbaum showed that the Cartesian
product satisfies the unique prime factorization property \cite{Fei86}.
Additionally, she provided a polynomial-time algorithm to determine the prime
factors which was improved to a linear time approach by Crespelle et al.\
\cite{CTL:13}.

\paragraph{Hypergraphs and the Cartesian Product}
Hypergraphs are a natural generalization of graphs in which ``edges'' may
consist of more than two vertices. More precisely, a hypergraph is a tuple
$H=(V,E)$ with non-empty set of vertices $V$ and a set of hyperedges $E$, where
each $e\in E$ is an ordered pair of non-empty sets of vertices $e=(t(e),h(e))$.
If $t(e)=h(e)$ for all $e\in E$ the hypergraph is called undirected and
directed, otherwise. Products of hypergraphs have been investigated by several
authors since the 1960s
\cite{Berge:Hypergraphs,Black15,Bretto06:Helly,Bretto:13b,Bretto09:HyperCartProd,Bretto:13,Doerfler79:CoversDirectProd,HNO:14,
Imrich67:Mengensysteme,Imrich70:SchwachKartProd,KA:12,KA:15,OstHellmStad11:CartProd,
Sonntag89:HamCart,Zhu:92}. It was shown by Imrich \cite{Imrich67:Mengensysteme}
that connected \emph{undirected} hypergraphs have a unique prime factor
decomposition (PFD) w.r.t.\ to the Cartesian product, up to isomorphism and the
order of the factors. A first polynomial-time factorization algorithm
for undirected hypergraphs was proposed by Bretto et al.\ \cite{Bretto:13}.

Unique prime factorization properties for \emph{directed} hypergraphs
were derived by Ostermeier et al.\ \cite{OstHellmStad11:CartProd}.
However, up to our knowledge, no
algorithm to determine the Cartesian prime factors of a connected directed
hypergraph is established, so-far.

\paragraph{Summary of the Results}
In this contribution, we present an algorithm to compute the PFD of connected
directed hypergraphs in $O(|V||E|r^2)$ time, where the rank $r$ denotes the maximum number of vertices
contained in the hyperedges. In addition, if we assume to have
hypergraphs with bounded rank the algorithm runs in $O(|E|\log^2(|V|))$ time.
Note, as directed hypergraphs are a natural
generalization of undirected hypergraphs, our method
generalizes and significantly improves the time-complexity of the method by Bretto et al.\ \cite{Bretto:13}.
In fact, the algorithm of Bretto et al.\
has time complexity $O(|V||E|\Delta^6r^6)$, where $\Delta$ is the maximum number of hyperedges a vertex is
contained in. Assuming that given hypergraphs have bounded rank $r$ and bounded maximum degree $\Delta$
this algorithm runs therefore in $O(|V||E|)$ time.

We shortly outline our method. Given an arbitrary connected directed hypergraph
$H=(V,E)$ we first compute its so-called 2-section $[H]_2$, that is, roughly
spoken the underlying \emph{undirected} graph of $H$. This allows us to use the
algorithm of Imrich and Peterin \cite{Imrich07:linear} in order to compute the
PFD of $[H]_2$ w.r.t.\ the Cartesian \emph{graph} product. As we will show, this
provides enough information to compute the Cartesian prime factors of the
directed hypergraph $H$. In distinction from the method of Bretto et
al.\ our algorithm is in a sense conceptually simpler, as \emph{(1)} we do not
need the transformation of the hypergraph $H$ into its so-called L2-section and
back, where the L2-section is is an edge-labeled version of the 2-section
$[H]_2$, and \emph{(2)} the test which (collections) of the factors of the
2-section are prime factors of $H$ follows a complete new idea based on
increments of fixed vertex-coordinate positions,
that allows an easy and efficient check to determine the PFD of $H$.

\section{Preliminaries}

\subsection{Basic Definitions}

A \emph{directed hypergraph} $H=(V,E)$ consists of a finite \emph{vertex} set
$V(H):=V$ and a set of \emph{directed hyperedges} or \emph{(hyper)arcs}
$E(H):=E$. Each arc $e\in E$ is an ordered pair of non-empty sets of vertices
$e=(t(e),h(e))$. The sets $t(e)\subseteq V$ and $h(e)\subseteq V$ are called the
\emph{tail} and \emph{head} of $e$, respectively. The set of vertices, that are
contained in an arc will be denoted by $V(e):=t(e)\cup h(e)$. If $t(e)=h(e)$
holds for all $e\in E$, we identify $e$ with $V(e)$, and we call $H=(V, E)$ an
\emph{undirected hypergraph}.
An undirected hypergraph is an
\emph{undirected graph} if $|V(e))|=2$ for all $e\in E$.
The elements of $E$ are called simply edges, if we consider an undirected graph.
The hypergraph with $|V|=1$ and $E=\emptyset$ is denoted by $K_1$ and is
called \emph{trivial}.

Throughout this contribution, we only consider hypergraphs without multiple
hyperedges and thus, being $E$ a usual set, and without loops, that is, $|V(e)|>1$
holds for all $e\in E$. However, we allow to have hyperedges being properly
contained in other ones, i.e., we might have arcs $e,f\in E$ with
$t(e)\subseteq t(f)$ and $h(e)\subseteq h(f)$.

A \emph{partial hypergraph} or \emph{sub-hypergraph} $H'=(V',E')$ of a
hypergraph $H=(V,E)$, denoted by $H'\subseteq H$, is a hypergraph such that
$V'\subseteq V$ and $E'\subseteq E$. The partial hypergraph $H'=(V',E')$ is
\emph{induced (by $V'$)} if $E' = \{e\in E\mid V(e)\subseteq V'\}$. Induced
hypergraphs will be denoted by $\left\langle V'\right\rangle$.

A \emph{weak path $P$ (joining the vertices $v_0,v_k\in V$)} in a hypergraph
$H=(V,E)$ is a sequence $P=(v_0,e_1,v_1,e_2,\ldots,e_k,v_k)$ of distinct
vertices and arcs of $H$, such that $v_0\in V(e_1)$, $v_k\in V(e_k)$ and $v_j\in
V(e_j)\cap V(e_{j+1})$. A hypergraph $H$ is said to be \emph{weakly connected}
or simply \emph{connected} for short, if any two vertices of $H$ can be joined
by a weak path. A \emph{connected component} of a hypergraph $H$ is a connected
sub-hypergraph $H'\subseteq H$ that is maximal w.r.t.\ inclusion, i.e., there is
no other connected sub-hypergraph $H''\subseteq H$ with $H'\subsetneq H''$.
Usually, we identify connected components $H'=(V',E')$ of $H$ simply by their
vertex set $V'$, since $\left\langle V'\right\rangle = H'$.

A \emph{homomorphism} from $H_1=(V_1, E_1)$ into $H_2=(V_2, E_2)$ is a mapping
$\phi: V_1\rightarrow V_2$ such that $\phi(e)$ is an arc in $H_2$ whenever $e$
is an arc in $H_1$ with the property that $\phi(t(e))=t(\phi(e))$ and
$\phi(h(e))=h(\phi(e))$. A bijective homomorphism $\phi$ whose inverse function
is also a homomorphism is called an \emph{isomorphism}.

The \emph{rank} of a hypergraph $H=(V,E)$ is $r(H)=\max_{e\in E}|V(e)|$.

The \emph{$2$-section $[H]_2$} of a (directed) hypergraph $H=(V,E)$ is the undirected
graph $(V,E')$ with $E'=\left\{xy:=\{x,y\}\subseteq V\mid x\neq y,\,\exists\;
e\in E: x,y\subseteq V(e)\right\}$. In other words, two vertices are linked by an edge
in $[H]_2$ if they belong to the same hyperarc in $H$. Thus, every arc $e\in E$
of $H$ is a \emph{complete graph} in $[H]_2$, i.e., all pairwise different
vertices in $V(e)$ are linked by an edge in $[H]_2$. Complete graphs defined on a vertex
set $V$ will be denoted by $K_{|V|}$.

We will also deal with equivalence relations, for which the following notations
are needed. For an equivalence relations $R$ we write $\varrho \mathrel{\mathsmaller{\mathsmaller{^{\boldsymbol{\sqsubseteq}}}}} R$ to
indicate that $\varrho$ is an equivalence class of $R$. A relation $Q$ is finer
than a relation $R$ while the relation $R$ is coarser than $Q$ if $(e,f)\in Q$
implies $(e,f)\in R$, i.e, $Q\subseteq R$. In other words, for each class
$\varrho$ of $R$ there is a collection $\{ \chi | \chi\subseteq \varrho\}$ of
$Q$-classes, whose union equals $\varrho$. Equivalently, for all $\varphi\mathrel{\mathsmaller{\mathsmaller{^{\boldsymbol{\sqsubseteq}}}}}
Q$ and $\psi\mathrel{\mathsmaller{\mathsmaller{^{\boldsymbol{\sqsubseteq}}}}} R$ we have either $\varphi\subseteq \psi$ or
$\varphi\cap\psi=\emptyset$.

\begin{rem}
If not stated differently, we assume that the hypergraphs considered in this
contribution are connected.
\end{rem}

\subsection{The Cartesian Product, (Pre-)Coordinates and (Pre-)Layers}

Let $H_1$ and $H_2$ be two hypergraphs. The \emph{Cartesian product} $H=H_1\Box
H_2$ has vertex set $V(H)=V(H_1)\times V(H_2)$, that is the Cartesian set
product of the vertex sets of the factors and the arc set
			\begin{align*}
				E(H)=&\big\{ ( \{x\}\times t(f),\{x\}\times h(f)) \mid x\in V(H_1), f\in  E (H_2)\big\} \cup\\
					  & \big\{ (t(e)\times \{y\}\ ,\ h(e)\times \{y\}) \mid e\in  E (H_1), y\in V(H_2)\big\}.
			\end{align*}

Thus, the tuple
$( \{(x_i,y_i), i\in I\} , \{(x'_j,y'_j), j\in J\}) \text{\ with\ }  x_i,x'_j\in V(H_1),\  y_i,y'_j\in V(H_2),\ i\in I, j\in J$
is an arc in $ E (H_1\Box H_2)$ if and only if either
			\begin{align*}
				\text{\ (i)\ }&   ( \{x_i, i\in I\} , \{x'_j, j\in J\}) \text{\ is an arc in\ }  E(H_1) \text{\ and\ } y_i=y'_j \text{\ for all\ } i\in I, j\in J \text{\ or\ }\\
				\text{\ (ii)\ }&  ( \{y_i, i\in I\} , \{y'_j, j\in J\}) \text{\ is an arc in\ }  E(H_2) \text{\ and\ } x_i=x'_j \text{\ for all\ } i\in I, j\in J.
			\end{align*}

The Cartesian product is associative, commutative, and the trivial one-vertex
hypergraph $K_1$ without arcs serves as unit \cite{HOS-12, OstHellmStad11:CartProd}.
Thus, for arbitrary finitely many factors $\{H_i, i\in I\}$ the product $\Box_{i
\in I} H_i$ is well-defined, and each vertex $x\in V(\Box_{i \in I} H_i)$ is
properly ``coordinatized'' by the vector $(x_i)_{i\in I}$ whose entries are the
vertices $x_i$ of the factors $H_i$.

A nontrivial hypergraph $H$ is \emph{prime} with respect to the Cartesian
product if it cannot be represented as the Cartesian product of two nontrivial
hypergraphs. A \emph{prime factor decomposition (PFD)} of $H$ is a
representation as a Cartesian product $H=\Box_{i\in I} H_i$ such that all factors $H_i$,
$i\in I$, are prime and $H_i\not\simeq K_1$. Note, the number $k$ of  prime
factors of $H = ( V , E )$ is bounded by $log(| V |)$, since every Cartesian
product of $k$ non-trivial hypergraphs has at least $2^k$ vertices.

Two important results concerning the Cartesian products of hypergraphs are given now.

\begin{lemma}[\cite{OstHellmStad11:CartProd}]
 The Cartesian product $H=\Box_{i=1}^n H_i$  of directed hypergraphs
	is connected if and only if all of its factors $H_i$ are connected.
\label{lem:Cart_connected}
\end{lemma}

\begin{thm}[\cite{OstHellmStad11:CartProd}]
	Connected (directed) hypergraphs have a unique prime factor decomposition
	with respect to the Cartesian product.
	\label{thm:upfd}
\end{thm}

We will show, that the PFD of a hypergraph $H$ can be obtained
from the PFD of its 2-section $[H]_2$. For this the following lemma is crucial.

\begin{lemma}
If $\Box_{i\in I} H_i$ is an arbitrary factorization of $H$ it holds
that $[H]_2 = \Box_{i\in I} [H_i]_2$.
\label{lem:Cart_2section}
\end{lemma}
\begin{proof}
	Since the Cartesian product is commutative and associative it suffices to
	prove the statement for two factors.
   Assume that $H=H_1\Box H_2$ and every vertex $x$ has coordinates $(x_1,x_2)$.
	Thus, there is an isomorphism
	$\phi:V(H)\to V(H_1\Box H_2)$ via $x\mapsto (x_1,x_2)$.
	We show that $\phi$ is also an isomorphism
	for the graphs $[H]_2$  and $[H_1]_2\Box [H_2]_2$.

	The edge $xy$ is contained in $E([H]_2)$ if and only if
	there is an arc $e\in E(H)$ with $x,y\in V(e) = t(e)\cup h(e)$ if and only if
 	(i)  $x_1=y_1$ and $x_2y_2 \subset V(f) \in E(H_2)$ or
	(ii) $x_2=y_2$ and $x_1y_1 \subset V(f) \in E(H_1)$ if and only if
 	(i)  $x_1=y_1$ and $x_2y_2 \in E([H_2]_2)$ or
	(ii) $x_2=y_2$ and $x_1y_1 \in E([H_1]_2)$ if and only if
	the edge $\phi(x)\phi(y)$ is contained in $E([H_1]_2\Box [H_2]_2)$.
\end{proof}

Now, given the PFD of $H=\Box_{i\in I} H_i$, we can infer that
$[H]_2=\Box_{i\in I} [H_i]_2$.
However, the factors $[H_i]_2$ might not be prime w.r.t. the Cartesian
\emph{graph} product and hence, $[H]_2$ might have more prime factors.
Since the PFD of $[H]_2 = \Box_{j\in J} G_j$ is unique it follows that the 2-section
of the prime factors $H_i$ of $H$ is a combination of the prime
factors $G_j$ of $[H]_2$, that is,
$[H_i]_2 \simeq \Box_{j\in J'} G_j$, $J'\subseteq J$ for all
$i\in I$.

Our algorithm will start with the PFD of $[H]_2=\Box_{j\in J} G_j$
w.r.t. the Cartesian product of undirected graphs and then tries to combine the respective
prime factors of $[H]_2$ to reconstruct the prime factors of $H$. In other words, we need
to find suitable subsets $J'\subseteq J$ so that $[H_i]_2 \simeq \Box_{j\in J'} G_j$
and $H_i$ is a prime factor of $H$.
To this end, we will introduce (pre-)coordinates and (pre)-layers.

\begin{defi}[(Factorization) Coordinatization]
Let $H$ be isomorphic to some product $\Box_{i \in I} H_i$, where
each factor $H_i$ has vertex set $\{1,\dots,l_i\}$.
A \emph{factorization coordinatization} or \emph{coordinatization} for
short, is an isomorphism $\Upsilon$ from $H$ to $\Box_{i \in I} H_i$.
Thus, $\Upsilon$ assigns to a vertex $v \in V(H)$ a vector of coordinates
$(v_i)_{i \in I}$ where $1 \leq v_i \leq l_i$ is a vertex in $V(H_i)$
\label{defi:coord}
\end{defi}

Hence, a coordinatization gives in an explicit way the information of the
underlying product structure of $H$. Hence, to find a factorization of $H$ one
can equivalently ask for a coordinatization of $H$, a fact that we will utilize
in our algorithm.
Note that the coordinatization w.r.t.\ a given
product decomposition is unique up to relabeling the vertices in each factor
$H_i$.

We will also need a notion which is similar to a coordinatization but
is implied by a factorization of the $2$-section $[H]_2$ rather than a decomposition of $H$.

\begin{defi}[Pre-Coordinatization]
Let $H$ be a given hypergraph and assume that $[H]_2$ has a coordinatization
$\Upsilon\colon V([H]_2) \to \mathop{\hbox{\Large{$\times$}}}_{i \in I} \{1,\dots,l_i\}$.
Since $V([H]_2)=V(H)$ we infer that $\Upsilon$ is a bijective
map on $V(H)$ that assigns to each vertex $v\in V(H)$ a unique coordinate-vector
$(v_i)_{i \in I}$ where $1 \leq v_i \leq l_i$.
This map is called \emph{pre-coordinatization} of $H$.
\label{defi:pre-coord}
\end{defi}

For convenience, we will usually omit the function $\Upsilon$
and identify every vertex $v\in V(H)$ with its (pre-)coordinate vector, i.e.,
we will write $v = (v_i)_{i \in I}$ rather than $\Upsilon(v)=(v_i)_{i \in I}$.

\begin{defi}[Layers and Pre-Layers]
Let  $H \simeq \Box_{i \in I} H_i$ with given respective coordinatization
$v = (v_i)_{i \in I} \in V(H)$ and $I' \subseteq I$.
The \emph{$I'$-layer through $v$} (denoted by $H_{I'}^v$) with
respect to this coordinatization is the sub-hypergraph induced by the vertices
$\{u=(u_i)_{i \in I} \mid i \notin I' \implies u_i = v_i \}$, i.e., we fix all
coordinates except those contained in the set $I'$.
Note, $H_{I'}^v \simeq \Box_{i\in I'} H_i$.

Analogously, as layers are defined by means of a coordinatization define
the \emph{pre-layers} by means of a pre-coordinatization.

For simplicity we write
$H_i^v$ instead of $H_{\{i\}}^v$ and $i$-(pre-)layer rather than $\{i\}$-(pre-)layers.
\label{defi:layer}
\end{defi}

For later reference, we need the following observation and lemma.
If $H \simeq \Box_{i \in I} H_i$ and $e$ is an arc of $H$,
then all vertices in $V(e)$ are contained in
the same $i$-layer $H_i^v$ for some $i \in I$ and $v\in V(H)$,
i.e.\ they only differ in the $i$-th coordinate. The same is true
for pre-layers.

\begin{lemma}
Let $H$ be a hypergraph and let $\Upsilon$ be a pre-coordinatization of $H$.
Then every arc $e$ of $H$ contains vertices of exactly one pre-layer w.r.t\ $\Upsilon$, that is,
all vertices in $V(e)$  only differ in the same $i$-th coordinate.
\label{lem:arc-in-prelayer}
\end{lemma}
\begin{proof}
Every hyperarc $e$ forms a complete subgraph $K_{|V(e)|}$
in $[H]_2$. Moreover, complete subgraphs must be contained entirely in one of the
$i$-layers of $[H]_2$, as complete graphs are so-called S-prime graphs, see e.g.
\cite{Bresar:03,Hel-12,HGS-09,Klavzar:02}. Hence, each hyperarc is contained in
one $i$-pre-layer of $H$.
\end{proof}

Note, any isomorphism from $[H]_2$ to $\Box_{i\in I} [H_i]_2$ and thus,
a  pre-coordinatization of $H$, is a coordinatization of $H$ if and
only if $H$ has a factorization $\Box_{i\in I} H_i$.
Lemma \ref{lem:Cart_2section} immediately implies that
every coordinatization of $H$ is also
a pre-coordinatization of $H$, while the converse is not true in general.
On the other hand, we have the following result for so-called increments of coordinates.

\begin{defi}[Increments of Coordinates]
Given a pre-coordinatization $\Upsilon\colon V(H) \to \mathop{\hbox{\Large{$\times$}}}_{i \in I} \{1,\dots,l_i\}$,
of $H$ and a vertex $v=(v_1,\dots,v_i,\dots,v_k) \in V(H)$
we define $\textrm{inc}(v,i)$ (w.r.t.\ $\Upsilon$) as the vertex with
coordinates $(v_1,\dots,v_i+1,\dots,v_k)$ where we set
$v_i+1:=1$ if $v_i+1>l_i$.

For an (ordered) set of vertices $W\subset V$ we define
the (ordered) set $\textrm{inc}(W,i))=\{\textrm{inc}(w,i)\mid w\in W\}$.

Finally, we denote for an arc $e=(t(e),h(e))$ its increment
$(\textrm{inc}(t(e),i), \textrm{inc}(h(e),i))$ by $\textrm{inc}(e,i)$.
\end{defi}

\begin{lemma}
	Let $H =(V,E)$ be a hypergraph and $\Upsilon\colon V\to \mathop{\hbox{\Large{$\times$}}}_{i \in I} \{1,\dots,l_i\}$
	be a pre-coordinatization of $H$.
	If for each arc $e\in E$ (where the vertices of $e$ differ only in the $j$-th coordinate),
	there is an arc $\textrm{inc}(e,i) \in E$ for all $i\neq j$,
	then $\Upsilon$ is a coordinatization of $H$.
	\label{lem:prec-inc}
\end{lemma}
\begin{proof}
	By Lemma	\ref{lem:arc-in-prelayer},
	all vertices within one arc $e\in E$ differ in
	precisely one coordinate.
	Let $e\in E$ be an arbitrary hyperarc and assume the vertices differ in the $j$-th coordinate.

	Let $\mathcal{H}_j$ be the set of $j$-pre-layers contained in $H$.
	Let $i\neq j$ be an arbitrary index $i\in I$. Assume that for each hyperedge $e$
	contained in some $j$-pre-layer $H(1)$ all ``incremental copies''
	$\textrm{inc}(e,i)$ are also contained in $E$, then there is a homomorphism from
	$H(1) = \langle V(H(1))\rangle$ to $H(2)=\langle(inc(V(H(1)),i))\rangle$,
	where $H(2)$ corresponds to some other $j$-pre-layer. Assume that for all such
	``consecutive'' $j$-pre-layer there is a homomorphism from $H(l)$ to
	$H(l+1)$, $1\leq l\leq l_i-1$. By construction, after $l_i-1$ incremental
	steps we arrive at the $l_i$-th $j$-pre-layer $H(l_i)$ and hence,
	$H(1)=\langle(inc(V(H(l_i)),i))\rangle$. If there is an homomorphism from $H(l_i)$
	to $H(1)$, then there is trivially an isomorphism between all
	such $j$-pre-layers $H(1),H(2),\dots,H(l_i)\in\mathcal H_j$.
	Thus, if for all arcs $e\in E$, where the
	vertices of $e$ differ precisely in this $j$-th coordinate, there is a hyperarc
	$\textrm{inc}(e,i) \in E$ for all $i\neq j$, then there isomorphism between all
	$j$-pre-layers contained in $\mathcal H_j$ for this fixed $j\in I$.

	If this is true for \emph{all} arcs $e\in E$, and thus, for
	\emph{all} $i$-pre-layers with $i\in I$, then \emph{all}
	such $i$-pre-layers are isomorphic for each $i\in I$.

	In particular, we can define for vertices $v,w$ and an index $i\in I$
	the map $g_i^{vw}\colon H_j^v \to H_j^w$ which maps every vertex in $H_i^v$ to the
	unique vertex in $H_i^w$ with the same $i$-coordinate.
	By the preceding arguments,
	for each $i\in J$ the map $g_i^{vw}$ is an isomorphism
	between the $i$-pre-layers in $H$ for all $v, w \in V(H)$.

	Finally, assume for contradiction that
	$\Upsilon$ is a not a coordinatization
	and hence,  $\Upsilon$ is not an isomorphism from $H$ to any product
	$\Box_{i \in I} H_i$. Hence, there must be some $i\in I$ such
	that not all $i$-layers are isomorphic by means of
	$g_i^{vw}$, a contradiction.
\end{proof}

Lemma \ref{lem:arc-in-prelayer} allows defining an equivalence relation
$R_\Upsilon$ on the hyperedge set $E(H)$ for a given pre-coordinatization
$\Upsilon\colon V(H) \to \mathop{\hbox{\Large{$\times$}}}_{i \in I} \{1,\dots,l_i\}$ of $H$, as follows:
$(e,f)\in R_\Upsilon$ if $e\in H_i^v$ and $f\in H_i^w$ for some $i\in I$ and
$v,w\in V(H)$. In other words, $e$ and $f$ are in relation $R_\Upsilon$ if they
are both contained in the $i$-pre-layers for the same fixed $i\in I$. Note, in
case that $\Upsilon$ is a coordinatization the relation $R_\Upsilon$ is also known
as product relation, that is, each equivalence of $R_{\Upsilon}$ contains the
hyperedges of all copies of some (not necessarily prime) factor of $H$. In order
to avoid confusion, we sometimes write that $R_{\Upsilon}(H)$ to indicate that
$R_{\Upsilon}$ is defined on the edge set of $H$.

Given two
pre-coordinatizations $\Upsilon_1$ and $\Upsilon_2$, we say that $\Upsilon_1$ is finer
than $\Upsilon_2$, while $\Upsilon_2$ coarser than $\Upsilon_1$ if $R_{\Upsilon_1}$ is
finer than $R_{\Upsilon_2}$.
We can immediately infer the next result.

\begin{lemma}
Let $H$ be a hypergraph and let $\Upsilon_1$ and $\Upsilon_2$ be
pre-coordinatizations of $H$. Then $\Upsilon_1$ is finer than $\Upsilon_2$ if and
only if the factorization of $[H]_2$ corresponding to $\Upsilon_2$ can be obtained
from the factorization $\Box_{i \in I} G_i$ corresponding to $\Upsilon_1$ by
combining some of these factors, i.e., $I$ can be partitioned into
$I_1,\dots,I_l$ so that the pre-coordinatization $\Upsilon_2$ is
an isomporphism between $[H]_2$ and $\Box_{j=1}^l G'_j$ where
$G'_j \simeq  \Box_{i \in I_j} G_i$.
\label{lem:prec-partition}
\end{lemma}
\begin{proof}
 By definition $\Upsilon_1$ is finer than  $\Upsilon_2$ if and only if
 $R_{\Upsilon_1}$ is finer than $R_{\Upsilon_2}$.
 By definition of the Cartesian product and by Theorem \ref{thm:upfd}
 this is the case if and only if every factor in the factorization
 corresponding to $\Upsilon_2$
 is a combination of the factors in the coordinatization
 $\Upsilon_1$.
\end{proof}

\begin{rem}
By Lemma \ref{lem:Cart_2section}, every coordinatization of $H$ is also a
pre-coordinatization of $H$. This implies together with Lemma
\ref{lem:prec-partition} that \emph{every} pre-coordinatization of $H$ can be
achieved by starting with the pre-coordinatization w.r.t.\ the prime
factorization of $[H]_2$ and then combining the corresponding pre-layers of $H$
to obtain the layers w.r.t.\ the prime factorization of $H$. In other words, one
needs to find a partition $\biguplus_{j \in J} I_j$ of the index set $I$ and
then combine all pre-layers corresponding to indices in the same part into one.
\label{rem:combine}
\end{rem}

As we shall see later, in our algorithm we will \emph{only} check increments w.r.t.\ the
pre-coordinatization $\Upsilon'$ coming from the PFD of $[H]_2$.
However, we have to prove that this is indeed sufficient (Theorem \ref{thm:correctness-pfd}).
In order to apply Lemma
\ref{lem:prec-inc} to validate whether we end up with a coordinatization $\Upsilon$ of
$H$ we would need to check increments with respect to this coarser
(pre-)coordinatization $\Upsilon$. Now one might hope that increments with respect to the
coarser pre-coordinatization are automatically increments with respect to the
finer pre-coordinatization $\Upsilon'$ or that at least the coarser pre-coordinates can be
chosen in a suitable way. However, this is not the case as the following example
shows.

Assume that at some point we need to combine $i$- and $j$-pre-layers of sizes
$3$ and $4$ respectively. The resulting $k$-pre-layers with respect to the new
coordinatization will each contain $12$ vertices labeled $1, \ldots , 12$. We now
claim that no matter how we assign the new labels, there is always at least one
increment $\textrm{inc}(\cdot,k)$ which is not an increment $\textrm{inc}(\cdot,i)$ or
$\textrm{inc}(\cdot,j)$. Assume for a contradiction that all increments $\textrm{inc}(\cdot,k)$
were either of the form $\textrm{inc}(\cdot,i)$ or $\textrm{inc}(\cdot,j)$. By applying
$\textrm{inc}(\cdot,k)$ recursively $12$ times to a vertex, we end up at the same vertex
again. This means, that we have applied $\textrm{inc}(\cdot,i)$ a number of times which
must be divisible by $3$ and $\textrm{inc}(\cdot,j)$ a number of times which must be
divisible by $4$. However, no suitable multiples of $3$ and $4$ add up to $12$.

The latter example shows that the single check of increments with respect to the
PFD of $[H]_2$ is not sufficient to invoke Lemma \ref{lem:prec-inc} to conclude
that some coarser pre-coordinatization is indeed a coordinatization. For this
purpose,  we need the following additional lemma.

\begin{lemma}
	Let $H =(V,E)$ be a hypergraph, let $\Upsilon_1,\Upsilon_2$ be
	pre-coordinatizations of $H$ such that $\Upsilon_1$ is finer than $\Upsilon_2$.
	Let $\textrm{inc}_1$ and $\textrm{inc}_2$ the respective increment maps. Assume that for each
	arc $e\in E$ (where the vertices of $e$ differ only in the $j$-th coordinate
	w.r.t.\ $\Upsilon_2$), there is an arc $\textrm{inc}_1(e,i) \in E$ for all $i\notin
	I_j$ (where $I_j$ is defined as in Lemma \ref{lem:prec-partition}).
	Then there is an arc $\textrm{inc}_2(e,k)$
	for all $k \neq j$ and hence $\Upsilon_2$ is a coordinatization of $H$.
	\label{lem:prec-coarsen-inc}
\end{lemma}

\begin{proof}
The coordinates of the vertices in $\textrm{inc}_2(e,k)$ w.r.t.\ $\Upsilon_1$ can be
obtained from those of vertices in $e$ by only changing coordinates outside
$I_j$. This can be achieved by successive applications of $\textrm{inc}_1(\cdot ,i)$ for
$i\notin I_j$. Since we started at an edge (namely $e$) and each of those
applications takes edges to edges we also end at an edge $\textrm{inc}_2(e,k)$. Hence,
Lemma \ref{lem:prec-inc} implies that $\Upsilon_2$ is a coordinatization.
\end{proof}

\section{PFD-algorithm for Directed Hypergraphs}

\subsection{Workflow}

We give here a summary of the workflow of the algorithm to compute the
prime factor of connected directed hypergraphs.
The top-level control structure is summarized in Algorithm \ref{alg:pfd}
\texttt{PFD\_of\_Di-Hypergraphs}\ in which the subroutines
\texttt{Preprocessing} (Alg.\ \ref{alg:pre}) and
\texttt{Combine} (Alg.\ \ref{alg:combine}) are used.

\begin{algorithm}[htbp]
\caption{\texttt{PFD\_of\_Di-Hypergraphs}}
\label{alg:pfd}
\begin{algorithmic}[1]
\small
    \STATE \textbf{INPUT:} A hypergraph $H=(V,E)$;
	 \STATE $H\gets\ $\texttt{Preprocessing}($H$); \COMMENT{Now, $H=(V_{\textrm{lex}},E_{\textrm{lex}})$ and vertex coordinates are known.}	 \label{alg:prepr}
	 \STATE Let $G_{\textrm{aux}} = (\{1,\dots,k\}, \emptyset)$ \COMMENT{Here, $k$ is the number of prime factors of $[H]_2$.};\label{alg:init-Gaux}
	\FOR {each $e\in	 E_\textrm{lex}$} \label{alg:for-start}
			\STATE Let $j\in \{1,\dots,k\}$ be the \emph{(unique)} coordinate where distinct $x,y\in V(e)$ differ; \label{alg:com1}
			\FOR{$i\in I=\{1,\dots,k\}\setminus\{j\}$} \label{alg:for2-begin}
					\IF {$\textrm{inc}(e,i)\notin E_\textrm{lex}$} \label{alg:if}
						\STATE add edge $ij$ to  $G_{\textrm{aux}}$;
					\ENDIF
			\ENDFOR\label{alg:for2-end}
	 \ENDFOR \label{alg:for-end}
		\STATE \texttt{Combine}($H=(V_{\textrm{lex}},E_{\textrm{lex}}), G_\textrm{aux}$); \label{alg:call-combine}
	 \STATE \textbf{OUTPUT:} PFD $\Box_{i=1}^n H_i$ of $H$;
\end{algorithmic}
\end{algorithm}
\begin{algorithm}[htbp]
\caption{\texttt{Preprocessing}}
\label{alg:pre}
\begin{algorithmic}[1]
\small
    \STATE \textbf{INPUT:} A connected hypergraph $H=(V,E)$;
	 \STATE Compute $\Box$-PFD of $[H]_2 = \Box_{i=1}^k G_i$ and vertex-coordinates with the Imrich-Peterin-Algorithm \cite{Imrich07:linear}; \label{alg:2sect-pfd}
	 \STATE Compute the  list  $V_{\textrm{lex}}$ of lexicographic ordered vertices (w.r.t.\ their coordinates);  \label{alg:lexV}
	 \FOR {each $e\in E$} \label{alg:for1-begin}
		\STATE Reorder $t(e)$ and $h(e)$ w.r.t.\ the lexicographic order of the vertices;
	 \ENDFOR \label{alg:for1-end}
	 \STATE Compute the list $E_{\textrm{lex}}$ of lexicographic ordered arcs;
			  w.r.t.\ to the lexicographic ordered sets $t(e)$ and then $h(e)$;
	\label{alg:lexE}
    \STATE \textbf{return} $(V_{\textrm{lex}},E_{\textrm{lex}})$ with respective vertex-coordinates;
\end{algorithmic}
\end{algorithm}
\begin{algorithm}[htbp]
\caption{\texttt{Combine}}
\label{alg:combine}
\begin{algorithmic}[1]
\small
    \STATE \textbf{INPUT:} A hypergraph $H=(V,E)$ with pre-coordinates w.r.t\ the PFD of $[H]_2 $, a graph $G_\textrm{aux}$;
		\STATE Compute connected components $I_1,\dots,I_l$ of $G_{\textrm{aux}}$; \label{alg3:conn}
		\FOR {each $e\in	 E_\textrm{lex}$} \label{alg3:fora}
			\STATE Let $j\in \{1,\dots,k\}$ be the \emph{(unique)} coordinate where distinct $x,y\in V(e)$ differ;
			\STATE Let $I_l$ be the connected component containing vertex $j$;
			\STATE Assign color $l$ to hyperedge $e$;
		\ENDFOR \label{alg3:fore}
		\STATE compute $[H]_2$ where edges $ij$ obtains the unique color of $e$ where $i,j\in V(e)$; \label{alg3:h2}
		\STATE compute coordinates of all vertices in $V$ in $[H]_2$; \COMMENT{This pre-coordinatization is the PFD-coordinatization of $H$} \label{alg3:coord}
	 \STATE \textbf{return} $\Box$-PFD $\Box_{i=1}^n H_i$ of $H$;
\end{algorithmic}
\end{algorithm}

As input of \texttt{PFD\_of\_Di-Hypergraphs}\  a connected hypergraph $H=(V,E)$ is
expected. First of all, subroutine \texttt{Preprocessing} is called. Here, the
PFD of $[H]_2 = \Box_{i=1}^k G_i$ and the respective coordinatization $\Upsilon$
of $[H]_2$ is computed by application of the algorithm of Imrich and Peterin
\cite{Imrich07:linear}.
 Then the
vertices, the vertices within the arcs and the arcs are ordered in
lexicographic order. This helps to achieve the desired time-complexity in later steps.

By definition, coordinatization $\Upsilon$ of $[H]_2$ is a pre-coordinatization of
$H$. By construction of $\Upsilon$ and Lemma \ref{lem:Cart_2section}, the
pre-coordinatization $\Upsilon$ is at least as fine as the coordinatization of $H$
w.r.t.\ its PFD. By Remark \ref{rem:combine} it suffices to find a suitable
partition of $I=\{1,\dots,k\}$ to derive the prime factors of $H$.
To this end, we initialize in Line
\ref{alg:init-Gaux} of Algorithm\ \ref{alg:pfd} the auxiliary graph $G_{\textrm{aux}}$,
where each vertex $i$ represents an element of $I$. The edge set is left empty.
We might later add edges in order trace back which equivalence classes of
$R_\Upsilon$ have to be combined, i.e., all vertices within one connected components of $G_{\textrm{aux}}$
will then be in one class of the respective partition of $I$.

We continue to check in the for-loop in Line \ref{alg:for-start}-\ref{alg:for-end}
of Algorithm\ \ref{alg:pfd}, if for each arc $e\in E_{\textrm{lex}}$ that is contained in
some $j$-layer its ``copies'' are also contained in ``incremental-neighboring''
$j$-layers, i.e., we check if $\textrm{inc}(e,i)\in E_\textrm{lex}$. If this is not the case,
then we add the edge $ij$ to $G_{\textrm{aux}}$.
Finally, we use the information of the connected components $I_1,\dots,I_l$ of
$G_{\textrm{aux}}$ that partition the set $I$ in order to determine the prime factors
of the given hypergraph $H$. To this end, the subroutine \texttt{Combine} is called
and an edge-colored 2-section $[H]_2$ is computed. That is, each edge $e\in E([H]_2)$ that
is contained in the copy of factor $G_i$ with $i\in I_s$ obtains color $s$.
In other words, all prime factors $G_j, j\in I_s$ of $[H]_2$ are combined
to a single factor of $[H]_2$ and the edges in the respective $I_s$-layers obtain color $s$.
W.r.t.\ this coloring it is possible to efficiently determine
new vertex coordinates in  $[H]_2$ which is then a factorization coordinatization
of $H$. We will show, that this leads to a ``finest'' coordinatization of $H$
and hence, to the prime factors $H_1,\dots,H_r$ of $H$.

For illustrative examples see Figures \ref{fig:exmpl1} and \ref{fig:exmpl2}.

\begin{table}[htbp]
\begin{figure}[H]
	\centering{
	\includegraphics[bb=0 0 491 144, scale=.9]{./exmpl1.eps}
	}
	\caption{The prime hypergraph $H$ (left-hand side) with arc set
				$E(H)=\{ {e_1=(\{11,12,13\},\{12,13,14\})},KA:
							{e_2=(\{21,22,23\},\{22,23,24\})},
							e_3=(\{14\},\{12\}),
							e_4=(\{11\},\{21\}),
            			{e_5=(\{12\},\{22\})},
							{e_6=(\{13,23\},\{13,23\})},
							{e_7=(\{14\},\{24\})},
							{e_8=(\{13\},\{23\})}
						 \}$
				admits a non-trivial prime factorization of its 2-section $[H]_2$ (right-hand side)
				into $K_4\Box K_2$. The vertices of $H$ are labeled w.r.t.\ its pre-coordinatization
				given by the coordinatization of $[H]_2$. For the first arc $e_4$ of
				the lexicographic ordered arc set $E_\textrm{lex}  = \{e_4,e_1,e_5,e_8,e_6,e_3,e_7,e_2\}$
				we can observe that the increment $\textrm{inc}(e_4,2) = ((\{12\}, \{22\}) = e_5 \in E_\textrm{lex}$.
				Analogously, $\textrm{inc}(e_1,1)=e_2$, $\textrm{inc}(e_5,2)=e_8$, and $\textrm{inc}(e_8,2)=e_7$
				are all contained in $E_\textrm{lex}$. However, when we arrive at the hyperedge $e_6$
				we obtain that $\textrm{inc}(e_6,2) = (\{14,24\},\{14,24\})$ is not an arc of $E_\textrm{lex}$.
				Hence, applying Algorithm \ref{alg:pfd} would lead to an edge $12$ in $G_\textrm{aux}$, resulting in a connected auxiliary graph
				and, $H$ would be determined as prime.
				A second example for an arc $e\in E_\textrm{lex}$ with $\textrm{inc}(e,i)\notin E_\textrm{lex}$ is the edge $e=e_7$.
				}
	\label{fig:exmpl1}
\end{figure}

\begin{figure}[H]
	\centering{
	\includegraphics[bb=0 0 409 300, scale=.9]{./exmpl2.eps}
	}
	\caption{The non-prime hypergraph $H = H_1 \Box H_2$
				is the product of the directed prime hypergraph $H_1$ in Fig.\ \ref{fig:exmpl1} and
				an undirected hypergraph $H_2$ with two vertices and one hyperedge.
				The vertices of $H$ are labeled w.r.t.\ its pre-coordinatization
				given by the coordinatization of $[H]_2$. The auxiliary graph $G_\textrm{aux}$ is initialized
				as the graph with three vertices and empty edge-set in Algorithm \ref{alg:pfd}.
				While the increment $\textrm{inc}(e,2) = \{(123,223),(123,223)\}$ of the arc $e=\{(113,213),(113,213)\}$ is
				still contained in $E_\textrm{lex}$, the increment $\textrm{inc}(e,3) = \{(124,224),(124,224)\}$ is not.
				Hence, the edge $13$ is added to $G_\textrm{aux}$. Since the increments of all hyperedges of the form
				$\{(i1j,i2j),(i1j,i2j)\}$ are contained in $E_\textrm{lex}$, no further edges will be added
				to $G_\textrm{aux}$. Hence, the sub-hypergraph induced by vertices with identical 2nd coordinate,
				i.e., the $\{1,3\}$-layers, constitute the copies the prime factor $H_1$, while the
				sub-hypergraphs induced by vertices with identical 1st and 3rd coordinate, i.e, the $2$-layers,
				are copies of the prime factor $H_2$.
				}
	\label{fig:exmpl2}
\end{figure}
\end{table}

\subsection{Correctness}

We are now in the position to prove the correctness of the algorithm \texttt{PFD\_of\_Di-Hypergraphs},
summarized in the following theorem.

\begin{thm}
	Algorithm \ref{alg:pfd} is sound and complete.
	\label{thm:correctness-pfd}
\end{thm}
\begin{proof}
	Given a hypergraph $H=(V,E)$. We start with a preprocessing and call in Algorithm\
	\ref{alg:pfd} the Algorithm\ \ref{alg:pre}. Here, the PFD of $[H]_2 = \Box_{i=1}^k
	G_i$ and the respective coordinatization $\Upsilon$ of $[H]_2$ is computed.
	This coordinatization $\Upsilon$ is by definition a
	pre-coordinatization of $H$. Since $[H]_2$ is an undirected graph, it is
	allowed to apply the algorithm of Imrich and Peterin \cite{Imrich07:linear}.
	Finally, the vertices, the vertices within the arcs and the arcs are
	ordered in lexicographic order. The latter task is not important for the
	correctness of the algorithm, but for the time-complexity that we will
	consider later on.

	We are now in Line \ref{alg:init-Gaux} of Algorithm\ \ref{alg:pfd}.
	By construction of $\Upsilon$ and Lemma \ref{lem:Cart_2section}, the pre-coordinatization $\Upsilon$
	is at least as fine as the coordinatization of $H$ w.r.t.\ its PFD. By Lemma
	\ref{lem:prec-partition} it suffices to find a suitable partition of
	$I=\{1,\dots,k\}$. To this end, we initialize in the auxiliary graph
	$G_{\textrm{aux}}$ where each vertex $i$ represents an element of $I$. The edge set
	is left empty. We might later add edges in order trace back which equivalence
	classes have to be combined, i.e., all vertices within one connected
	components of $G_{\textrm{aux}}$ will then be in one class of the respective
	partition of $I$.

	Now consider the for-loop in Line \ref{alg:for-start}-\ref{alg:for-end}. For
	each $e\in E_{\textrm{lex}}$ we check in which coordinates the vertices in $V(e)$
	differ. Since $\Upsilon$ is a pre-coordinatization and by Lemma \ref{lem:arc-in-prelayer},
	this is exactly one	coordinate for each hyperarc. Let $e\in E_{\textrm{lex}}$ be a chosen arc and assume
	that all vertices in $V(e)$ differ in the $j$-th coordinate. Now, it is
	checked if for arc $e$ its ``copies'' are contained in each
	$G_j^w$-layer where $w\in inc(V(e),i)$. If for some arc $e\in E_{\textrm{lex}}$ we
	observe that there is no hyperedge $\textrm{inc}(e,i) = (\textrm{inc}(t(e),i), \textrm{inc}(h(e),i)\in E_\textrm{lex}$ then
	there is no ``copy'' of $e$ in some $j$-pre-layer through $w$ with $w\in
	inc(V(e),i)$. In this case we add the edge $ij$ to $G_\textrm{aux}$ if not already
	set. The latter tasks are repeated for all hyperarcs $e\in E_{\textrm{lex}}$.

	Finally, in Line \ref{alg:call-combine} the Algorithm\ \ref{alg:combine}
	is called. The task of this subroutine is to combine the pre-coordinates
	and thus, the pre-layers in order to determine the layers of the
	final prime-factors of $H$.
   Let $I_1,\dots,I_r$ be the connected components of $G_\textrm{aux}$.
	Clearly, $\Pi=\{I_1,\dots, I_r\}$ is a partition
	of $I$. Let each $I_j$ having $l_j$ elements.
	Lemma \ref{lem:prec-partition} and Remark \ref{rem:combine} imply that
   $\widehat{\Upsilon}\colon V(H) \to \mathop{\hbox{\Large{$\times$}}}_{l=1}^r \{1,\dots,l_r\}$
	is  a pre-coordinatization of	$H$.
	It remains to show that
	\begin{enumerate}
		\item[(1.)]  $\widehat{\Upsilon}$ is a coordinatization and
		\item[(2.)]  $\widehat{\Upsilon}$ is at least as fine as the coordinatization given by the PFD of $H$.
	\end{enumerate}

	\emph{Claim (1.):}
	By construction, all $e\in E_\textrm{lex}$ where the vertices differ in the $i$-th coordinate
	w.r.t.\ $\Upsilon$ are now contained in some $I_s$-layer where $i\in I_s\in \Pi$.
	Moreover, for all $e$ in some $I_s$-layer the increments $\textrm{inc}(e,j)$ with
	$j\neq i$ and $j\notin I_s$ must be contained in $H$, as otherwise we would
	have added the edge $ij$ to $G_\textrm{aux}$ and hence, $j\in I_s$.
	As the latter is true for all $I_s$-layers contained in $\Pi$
	we can apply Lemma \ref{lem:prec-coarsen-inc} and conclude that
	 $\widehat{\Upsilon}$ is a coordinatization of $H$.

	\emph{Claim (2.):}
	Given the pre-coordinatization $\Upsilon$ of $H$.
	By construction of $\Upsilon$ and Lemma \ref{lem:Cart_2section},
   $\Upsilon$
	is at least as fine as the coordinatization of $H$ w.r.t.\ its PFD.
	Thus, there is a partition $\Pi'=\{I'_1,\dots, I'_t\}$ of $I$
	 w.r.t.\ the PFD of $H$.
	It remains to show that if there are two indices $i,j \in I_s\in \Pi$,
   then $i,j$ are also contained in the same class of $\Pi'$.
	If $i,j \in I_s\in \Pi$ then they are in same connected component
	$C_s$ of $G_\textrm{aux}$. Hence, it suffices to consider pairs
	$i,j \in I_s\in \Pi$ that are connected by an edge.
	Assume, for contradiction that $i,j$ are in different classes of $\Pi'$.
	W.l.o.g. let $i\in I'_1$ and $j\in I'_2$.
	Moreover, let $H=H_{I'_1}\Box H_{\cup_{l\geq2} I'_l}$
	Hence, for \emph{all} $l\in I'_1$ and thus, in particular for $l=i$
	it holds that for all arcs $e$ in some $I'_2$-layer there is an arc
	$(\textrm{inc}(t(e),i), \textrm{inc}(h(e),i)\in E_\textrm{lex}$.
	The same holds with the role of $i$ and $j$ switched.
	However, in this case we would not add the edge $ij$ to $G_\textrm{aux}$,
	a contradiction.

	To finish the PFD-computation we have to compute $\widehat\Upsilon$.
	To this end, we compute the 2-section $[H]_2$ with edges $xy$ colored
	with color $j$ whenever $x$ and $y$ are contained in some edge $e$
	that is contained in some $j$-layer of $H$.
	Lemma \ref{lem:Cart_2section} implies that  $\widehat\Upsilon$
   is also 	a coordinatization of $[H]_2$ and hence,
	all edges with same color $j$ in $[H]_2$ are contained in the same
	equivalence class of $R_{\widehat\Upsilon}([H]_2)$.
	In other words,  $(e,f) \in R_{\widehat\Upsilon}(H)$ if and only if
	 $(xy,uv) \in R_{\widehat\Upsilon}([H]_2)$ for all distinct $x,y\in V(e)$
	and distinct $u,v\in V(f)$. By construction, $R_{\widehat\Upsilon}([H]_2)$ is a
	a product relation of $[H]_2$, and thus we can apply  again
	a method proposed the by Imrich and Peterin (cf.\ Theorem 5.1. in \cite{Imrich07:linear}),
	in order to obtain the desired coordinates and hence, $\widehat\Upsilon$.
\end{proof}

\subsection{Time Complexity}

In order to prove the time-complexity results, we first give the following lemma.

\begin{lemma}
Let $H$ be a hypergraph, let $[H]_2 = \Box_{i=1}^k H_i$ be a factorization of its $2$-section into $k$ factors, and let $m$
and $n$ be the number of arcs and vertices of $H$, respectively.
Then for any $l \in \mathbb N_0$ it holds that $k^l \log m = O(n)$.
\label{lem:bound}
\end{lemma}
\begin{proof}
If $m_i$ and $n_i$ are the numbers of arcs and vertices of the factors then
$m_i \leq (2^{n_i})^2$, since we have to consider tail and head independently.
Let $N$ be the maximum number of vertices of a
factor. Then we have
\[
	m
	= \sum_{i=1}^k m_i \prod_{\substack{j=1 \\j \neq i}}^k n_j
	\leq \sum_{i=1}^k 2^{n_i} \prod_{j=1}^k n_j
	\leq k \cdot 2^N \prod_{j=1}^k n_j.
\]
Taking logarithms on both sides of the inequality gives
\[
	\log m
	\leq \log k + N + \sum_{j=1}^k \log(n_i)
	\leq k + N + k\cdot N
	\leq c \cdot k \cdot N
\]
for some suitable constant $c$.

On the other hand by bounding the size of every factor except the biggest one
from below by $2$ we get $n = \prod_{i=1}^k n_i \geq N \cdot 2^{k-1}$.
Clearly, $N \cdot 2^{k-1}\geq c' \cdot N \cdot k^l$ for some suitable constant $c'$
depending on $l$. Together with the estimate for $\log m$ this proves the lemma.
\end{proof}

The next two lemmas are concerned with the time-complexity of the subroutines
\texttt{Preprocessing} and \texttt{Combine}.

\begin{lemma}
	Let $H = (V,E)$ be a connected hypergraph with $|V|=n$, $|E|=m$ and rank $r$.
	Then Algorithm \ref{alg:pre} performs in $O(r^2mn)$ time.
	If we assume that $H$ has
	bounded rank, then  Algorithm \ref{alg:pre} has time-complexity $O(m\log^2(n))$.
	\label{lem:runtime-pre}
\end{lemma}
\begin{proof}
	In Line \ref{alg:2sect-pfd}, the first task is the computation of the 2-section $[H]_2$.
	To this end, we initialize an
	adjacency list  $N[1],\dots,N[n]$ with empty entries,
	 which can be done in $O(n)$ time. We add
	for each arc $e\in E$ and each pair $\{i,j\}\in \binom{V(e)}{2}$ the vertex
	$i$ to $N[j]$ and $j$ to $N[i]$, if these vertices are not already contained
	in the respective adjacency lists. Hence, we must check whether $i\in N[j]$ or not.
	To this end, assume that $N[j]$ is already ordered. Hence we need
	$O(\log(n))$ comparisons to verify if $i\in N[j]$. If this is not the case,
	the vertex $i$ is added to $N[j]$ on the respective position so that $N[j]$ stays sorted.
	Analogously, we add $j$ to $N[i]$, whenever $j$ is not contained in $N[i]$.
	As for each arc $e\in E$ there are at most $\binom{r}{2} = O(r^2)$ pairs
	$\{i,j\}$ and	for each such pair we have $O(\log(n))$ comparisons we end in a time-complexity of
	$O(n+mr^2 \log(n))$ to create the adjacency list
	$N[1],\dots,N[n]$. These lists serve
	than as input for the algorithm of Imrich and Peterin which computes the PFD
	of the 2-section in $O(|E([H]_2)|+n)$ time. Since $[H]_2$  is connected
	and thus, $[H]_2$ has at least $|V|-1$ edges,
	the PFD algorithm runs in fact in $O(|E([H]_2|)=O(mr^2)$ time. Hence, the total
	time complexity of Line \ref{alg:2sect-pfd} of Algorithm\ \ref{alg:pre} is $O(n+ mr^2
	\log(n) + mr^2) = O(mr^2 \log(n))$.

	In what follows, let $k$ be the number of factors of $[H_2]$
	and let each $v\in V$ be identified with its respective (pre-)\
	coordinate vector $(v_1,\dots,v_k)$ computed by the  Imrich-Peterin-Algorithm.
	Note, $k$ is bounded by $\log(n)$.

	In Line \ref{alg:lexV}, the list $V$ of vertices is reordered in lexicographic
	order w.r.t.\ the vertex coordinates, i.e., $v<w$ if there is some
	$i\in\{1,\dots,k\}$ with $v_j\leq w_j$ for all $j\in\{1,\dots,i-1\}$ and
	$v_i<w_i$. This task can be done in $O(n\log(n)k) = O(n\log^2(n))$. Since
	$mr\geq n$ we obtain that $O(n\log^2(n)) = O(mr\log^2(n))$. This new ordered
	vertex list is called $V_\textrm{lex}$.

	We are now concerned with the for-loop in Line \ref{alg:for-start}. For each
	hyperedge $e\in E$ we reorder the vertices of its head and tail w.r.t.\ to the
	order of the vertices in $V_\textrm{lex}$. Each hyperedge contains at most $r$ vertices and
	hence, this task can be done in $O(r\log(r))$ time. Therefore, the entire
	for-loop (Line \ref{alg:for1-begin} - \ref{alg:for1-end}) takes
	$O(mr\log(r))$ time.

	Finally, the arcs are reordered w.r.t. the lexicographic ordered sets $t(e)$
	and $h(e)$. We say $e<f$ if $t(e)<t(f)$ or $t(e)=t(f)$ and $h(e)<h(f)$,
	whereby the tails, resp., heads are compared w.r.t. the lexicographic order
	of their vertices. To determine if $t(e)<t(f)$ or $h(e)<h(f)$ for some arcs
	$e,f\in E$, the at most $2r$ pairs of vertices must be compared, whereby the
	comparison of each such pair can be done in $O(1)$ time, since the vertices
	are already ordered in the tails and heads. The reordering of the arcs need
	than $O(m\log(m))$ comparisons, where each comparison can be done in $O(r)$
	time, by the preceding arguments. Hence, the creation of $E_\textrm{lex}$ takes $O(r
	m\log(m))$ time. By Lemma \ref{lem:bound} this is $O(rmn)$. Moreover, if we
	assume that the rank $r$ is bounded, then $m \leq \sum_{i=1}^r \binom nr \leq rn^r = O(n^r)$. Hence $O(\log(m))
	= O(\log(n^r)) = O(r\log(n))=O(\log(n))$.
	In this case the time complexity for
	determining $E_\textrm{lex}$ is $O(m\log(m)) = O(m\log(n))$.

	Taken together the latter arguments, we end in overall time complexity for
	Algorithm \ref{alg:pre} of $O(r^2mn)$
	and if the rank $r$ is bounded with $O(m\log^2(n))$.
\end{proof}

\begin{lemma}
	Let $H = (V,E)$ be a connected hypergraph with $|V|=n$, $|E|=m$ and rank $r$.
	Moreover, assume that the graph $G_\textrm{aux}$ has $k$ vertices and $m'$ edge with
	$k\leq \log(n)$.
	Then Algorithm \ref{alg:combine} preforms in $O(mnr^2)$ time.
	If we assume that $H$ has
	bounded rank, then  Algorithm \ref{alg:combine} has time-complexity $O(m\log^2(n))$.
	\label{lem:runtime-post}
\end{lemma}
\begin{proof}
   Determining the connected components of $G_\textrm{aux}$ in Line \ref{alg3:conn}
	can be done in $O(k+m')=O(\log(n)+\log^2(n))$ time by application of
   the classical breadth-first search.
	While doing this, we will in addition record in $O(1)$ time for each vertex
	in which connected component it is contained.
	Let $I_1,\dots,I_l$ be the connected components of $G_\textrm{aux}$.

	For each of the $m$ arcs we have to find the indices where the vertices of
	the particular arc differs. To this end, it suffices to take any two vertices
	$x$ and $y$ of $V(e)$ and to compare their $k$ coordinates which takes $O(k)$
	time. Let $j$ be the coordinate where the two vertices differ. We need to
	check in which of the connected components $I_s$ the vertex is contained in,
	which can be done in $O(1)$ time, since we have already recorded for each
	vertex of $G_\textrm{aux}$, in which component it is contained in.
	Now, the color for each arc can be
	recorded in $O(1)$ time.
   Hence, the for-loop (Line \ref{alg3:fora}-\ref{alg3:fore}) has
	overall-time complexity $O(mk) = O(m\log(n))$.

	To compute the 2-section in Line \ref{alg3:h2} with colored edges we
	initialize an extended adjacency list $N[1],\dots,N[n]$ where whenever we add
	some $i\in N[j]$ we also record the respective unique color of $ij$ as a 2nd
	parameter. Recording this parameter can be done in $O(1)$ time, as for each
	arc $e\in E$ it is known which color it has. Hence, we can argue analogously
	as in the proof of Lemma \ref{lem:runtime-pre}, and state that the 2-section
	with additionally colored edges
	can be computed in $O(n+mr^2 \log(n))$ time.

	Finally, the vertex-coordinates in $[H]_2$ can be computed in $O(m')=O(mr^2)$
	time, see Theorem 5.1. in \cite{Imrich07:linear}.

	Hence the overall-time complexity of Algorithm \ref{alg:combine}
	is $O(\log^2(n) + m\log(n) + n+mr^2 \log(n) + mr^2) = O(n+m\log^2(n)r^2)$.
	Since $mr\geq n$ and $\log^2(n) = O(n)$,
	the latter can be expressed as $O(mnr^2)$.
   If we assume in addition that the rank $r$ is bounded we get
	$O(n+m\log^2(n)r^2) =O(mr+m\log^2(n)r^2) = O(m\log^2(n))$.
\end{proof}

We are now in the position to determine the time-complexity of algorithm
\texttt{PFD\_of\_Di-Hypergraphs}.

\begin{thm}
	Let $H = (V,E)$ be a connected hypergraph with $|V|=n$, $|E|=m$ and rank $r$.
	Then Algorithm \ref{alg:pfd} computes the PFD of $H$ in $O(mnr^2)$ time.
	If the rank $r$ is bounded the time-complexity of Algorithm \ref{alg:pfd}
	is $O(m\log^2(n))$.
	\label{thm:runtime-pfd}
\end{thm}

\begin{proof}
	We suppose both the vertices and the hyperarcs of $H$ implemented
	as integers and $E$ implemented as an $m\times 2$ array,
	where each entry $E[e,i]$ contains the list of vertices in $t(e)$
	if $i=1$ and $h(e)$	if $i=2$.
	In Line \ref{alg:prepr} we call \texttt{Preprocessing}($H$) which takes
   $O(r^2mn)$ time and if $r$ is bounded $O(m\log^2(n))$ time
	(Lemma \ref{lem:runtime-pre}).

	In what follows, let $k\leq\log(n)$ be the number of factors of $[H_2]$
	and assume that each $v\in V$ is identified with its respective (pre-)coordinate vector
	$(v_1,\dots,v_k)$.

	In Line \ref{alg:init-Gaux} the auxiliary graph is be initialized.
	In particular, we initialize $G_\textrm{aux}$ as adjacency list, i.e., we create
	empty lists $N[1],\dots,N[k]$ which can be done  in $O(k)$ time.

	We are now concerned with the for-loop in Line \ref{alg:for-start} -
	\ref{alg:for-end}. For each of the $m$ arcs we have to find the
	indices where the vertices of the particular arc differs. To this end,
	any two vertices $x$ and $y$ of $V(e)$ are chosen and their
	$k$ coordinates are compared, which takes $O(k)$ time. The nested for-loop (Line
	\ref{alg:for2-begin} - \ref{alg:for2-end}) is executed for all coordinates $i$
	where the vertices of arc $e$ are identical and it is checked whether $\textrm{inc}(e,i)=(\textrm{inc}(t(e),i),
	\textrm{inc}(h(e),i)$ is contained in $E_\textrm{lex}$ or not. The increment $(\textrm{inc}(t(e),i),
	\textrm{inc}(h(e),i)$ can be computed in $O(r)$ time. Note, the vertices within
	$\textrm{inc}(t(e),i)$ and $\textrm{inc}(h(e),i)$ are still lexicographically ordered as only
	vertex-coordinates are incremented that have been identical for the vertices
	within the arc and thus,their $i$-th positions are
	all still equal after the computation of $\textrm{inc}(e,i)$. We now check
	whether $(\textrm{inc}(t(e),i), \textrm{inc}(h(e),i)\notin E_\textrm{lex}$. Since $E_\textrm{lex}$ is already
	ordered, binary search finds the corresponding arc using at most $O(\log(m))$
	comparisons of arcs and since head and tail of each arc are in lexicographic order
	comparing two arcs takes $O(r)$ time.
	Therefore, the if-condition in Line \ref{alg:if} takes $O(r+r\log(m)) =
	O(r\log(m))$ time.
	In case,	$(\textrm{inc}(t(e),i), \textrm{inc}(h(e),i)\notin E_\textrm{lex}$ we have to add a respective edge
	$ij$ in $G_\textrm{aux}$, if not already set. Hence, to check whether $ij$ exists in
	$G_\textrm{aux}$, we need to validate if $i\in N[j]$. To this end, assume that $N[j]$
	is already ordered. Hence we need $O(\log(k))$ comparisons to verify if $i\in
	N[j]$. If this is not the case $i$ is added to $N[j]$ on the respective
	position so that $N[j]$ stays sorted. Similarly, $j$ is added to $N[i]$
	whenever $i\not\in N[j]$.
	Hence the nested for-loop in (Line \ref{alg:for2-begin} -
	\ref{alg:for2-end}) has time complexity $O(k(r\log(m)+\log(k))) =
	O(kr\log(m))$, since the number of arcs $m$
	is at least as big as the number (non-trivial) factors $k$.
	Take together the latter arguments, the entire for-loop
	in Line \ref{alg:for-start} -	\ref{alg:for-end} has time-complexity
	$O(m(k+kr\log(m))) = O(mkr\log(m))$.
	By Lemma \ref{lem:bound} this is $O(mnr)$.
	Moreover, if we assume that the rank $r$ is bounded, then
	$m\leq n^r$ and hence, $O(\log(m)) = O(\log(n^r)) = O(r\log(n))=O(\log(n))$
	In this case, the time complexity of Line \ref{alg:for-start} -	\ref{alg:for-end}
	is $O(mk\log(m)) = O(mk\log(n)) =O(m\log^2(n))$.

	Finally, we use Algorithm \ref{alg:combine} which performs in $O(mnr^2)$ time and
	if the rank $r$ is bounded it has time-complexity $O(m\log^2(n))$ (Lemma \ref{lem:runtime-post}).

	To summarize, each step of Algorithm \ref{alg:pfd} can be performed in
	$O(mnr^2)$ time and if the rank $r$ is bounded the time-complexity
	is $O(m\log^2(n))$.
\end{proof}

\section*{Acknowledgment}
We thank the organizers of the 8th Slovenian Conference on Graph Theory (2015)
in Kranjska Gora, where the authors participated, met and basically drafted the main ideas
of this paper, while drinking a cold and tasty red Union, or was it
a green La{\v{s}}ko?

We also thank Wilfried Imrich and Iztok Peterin for helpful comments regarding
the time-complexity of our algorithm.

\bibliographystyle{elsarticle-harv}


\section{HESS J1943+213 Classification}
HESS J1943+213 was discovered in very-high-energy (VHE; E $\textgreater$ 100~GeV) gamma rays during
the H.E.S.S. Galactic plane scan~\cite{hess}. Due to its point-like appearance in VHE gamma rays, three possible
source classes were suggested: gamma-ray binary, pulsar wind nebula (PWN), and BL Lac object (blazar).

Assuming the source is a gamma-ray binary, ref.~\cite{hess} used the lack of detection of a massive (O- or B-type) companion star to estimate a distance limit of greater than $\sim$25 kpc. This distance would place the binary well beyond the extent of the Galactic disk and would imply an X-ray luminosity 100-1000 times higher than luminosities of known gamma-ray binaries. Such a distance limit is problematic and this scenario is disfavored. In addition, the point-like appearance in the X-rays and the soft VHE spectrum with a power-law index of $\Gamma$=3.1 $\pm$ 0.3 (in contrast to the softest known PWN index of 2.7) motivated ref.~\cite{hess} to argue against the PWN scenario, leaving the blazar hypothesis. The authors found all observations to be consistent with the blazar scenario: the point-like nature in both X-rays and VHE, the soft VHE spectral index, and a preliminary IR spectrum showing lack of emission lines. Moreover, the hard X-rays observed with INTEGRAL IBIS and \emph{Swift} BAT instruments showed no evidence of a cutoff up to an energy of $\sim$195~keV. If the source is a blazar, it would be categorized as an extreme high-synchrotron-peak BL Lac object (extreme HBL), a class of blazars with the synchrotron peak located at energies \textgreater1 keV~\cite{ehbl}.

Since the discovery publication, the identity of HESS J1943+213 has been the topic of an ongoing debate. 1.6-GHz VLBI observations of the HESS J1943+213 counterpart with the European VLBI Network produced
a detection that was claimed to show extension, with FWHM angular size of 15.7~mas (with 3.5~mas expected size for a point source)~\cite{gabanyi}.
 Based on this measurement, the brightness temperature of the counterpart was estimated to be 7.7$\times$10$^{7}$ K and was used to argue against the blazar scenario, as the expected brightness temperature of HBLs is in the 10$^{8}$--10$^{9}$ K range. In addition, ref.~\cite{gabanyi} employed a 1$'$ feature observed in the 1.4-GHz VLA C-array configuration image to support the PWN hypothesis, with the assertion that the angular size of the feature is consistent with a Crab-like PWN placed at a distance of 17~kpc. On the other hand, ref.~\cite{tanaka} argued in favor of an extreme HBL by constructing a spectral energy distribution and drawing comparisons to a known extreme HBL, 1ES 0347-121. More recently, the case for the extreme blazar scenario has been bolstered with ref.~\cite{peter} observing the near-infrared (K-band) counterpart of HESS J1943+213, and claiming potential detection of an elliptical host galaxy with only 10\% probability the object is a star. Using 5 years of data, they also obtain the first significant detection of the source with Fermi LAT (within 2$\sigma$ positional uncertainty) in the 1--300 GeV energy regime and use it in conjunction with the H.E.S.S. spectrum to set an upper limit on the redshift at z < 0.4.
Finally, ref.~\cite{gab2015} obtained follow-up observations in the 1.5 GHz and 5 GHz bands using the e-Multiple Element Remotely Linked Interferometer Network (e-MERLIN), showing that the source exhibits a flat spectrum between the two bands and claiming a detection of flux density variability in the 1.5 GHz band when compared with the EVN observations of the source.
Though recent evidence is stacking in favor of the blazar hypothesis, there is yet to be a definitive identification of the HESS J1943+213 source class.

Here we present results from VERITAS observations of HESS J1943+213 and discuss implications for the debate regarding the source identity.

\section{VERITAS Observations}
The Very Energetic Radiation Telescope Array System (VERITAS) is an array of four 12-m imaging
atmospheric Cherenkov telescopes, located in Arizona, USA at an elevation of $\sim$1270~m. The camera for
each telescope is composed from 499 photo-multiplier tubes, with a field of view of approximately 3.5$^{\circ}$.
VERITAS is able to reliably reconstruct VHE gamma rays with energies between 85~GeV and \textgreater30~TeV~\cite{holder}, with an energy resolution of 15-25\% and a sensitivity to detect a 1\% Crab Nebula flux source at 5$\sigma$ in 25 hours.

\begin{figure}[h]
\includegraphics[width=120mm]{HESSJ1943_sigmap_fix.pdf}
\centering
\caption{Significance sky map of VERITAS observations of HESS J1943+213. The white circle represents the VERITAS angular resolution 68\% containment at 1 TeV; the cross shows the catalog position of the source.}
\label{sigmap}
\end{figure}

VERITAS observed HESS J1943+213 between May 27, 2014 (MJD 56804) and July 2, 2014 (MJD 56840) with
22.5~hours of total live time after data affected by poor weather conditions are eliminated. Observations took place at elevations between 63$^{\circ}$ and 80$^{\circ}$, leading to 17.5$\sigma$ source detection above 200~GeV.

\begin{figure}[!h]
\includegraphics[width=100mm]{HESSJ1943_spectrum_HESS.pdf}
\centering
\caption{Differential energy spectrum of HESS J1943+213. The spectrum with VERITAS (black points) is between 200~GeV and 2~TeV and is fit to a power-law function. The spectrum from the H.E.S.S. result \cite{hess} is also included for comparison (gray points).}
\label{spectrum}
\end{figure}

\begin{figure}[!h]
\includegraphics[width=100mm]{HESSJ1943_LC.pdf}
\centering
\caption{Light curve of HESS J1943+213 with VERITAS, with flux measured above 200~GeV. The red line is a constant flux fit to the data.}
\label{lc}
\end{figure}

\section{Results}
Figure~\ref{sigmap} shows the significance sky map near HESS J1943+213, including the detected source. The source location is
consistent with the catalog position of HESS J1943+213 (marked with a black cross).
The VERITAS differential energy spectrum of the source is shows in Figure~\ref{spectrum}. The high elevation observations lead to a lower energy threshold and a spectrum that extends down to 200~GeV, compared with 470~GeV from H.E.S.S. observations. The source spectrum is well fit by a power-law function with an index of $\Gamma$=2.82 $\pm$ 0.13 in the energy range 200~GeV--2~TeV. The VERITAS spectrum of HESS J1943+213 appears to be consistent with the spectrum from H.E.S.S. ($\Gamma$=3.1 $\pm$ 0.3), also depicted in Figure~\ref{spectrum}.

The flux of (1.30 $\pm$ 0.20) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ measured with VERITAS above 470 GeV is consistent with the flux of (1.25 $\pm$ 0.20) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ from the H.E.S.S. detection. As VERITAS is able to observe HESS J1943+213 at a much higher elevation than H.E.S.S., the detection rate of the source with VERITAS is $\sim$3.7$\sigma$/$\sqrt{\textrm{hour}}$, compared with $\sim$1.8$\sigma$/$\sqrt{\textrm{hour}}$ with H.E.S.S., allowing VERITAS to test for variability on a factor of four shorter timescales.

Figure \ref{lc} illustrates the light curve of HESS J1943+214 observed with VERITAS, with fluxes measured above 200~GeV (showing only statistical uncertainties). A constant flux line is fit to the data (shown in red) with a $\chi^{2}$ / NDF = 31.34 / 20, corresponding to a p-value of 0.05. With systematic uncertainties in flux included, any hints of variability are not statistically significant.

\section{\emph{Swift}-XRT Observations and Results}

Soft X-ray observations of HESS J1943+213 were obtained with the \emph{Swift} XRT-instrument on 17 June, 19 June, and 21 June, 2014, contemporaneous with VERITAS observations. The data analysis for XRT observations was performed with the standard tools included in HEASoft package Version 6.15.1, specifically with the XRTDAS v3.0.0 tools, while spectral analysis and fitting was performed using Xspec v12.8.1g.

XRT flux spectra were obtained by unfolding the counts spectra with instrument response functions
included in CALDB 1.0.2 and by assuming an absorbed power-law functional form for the intrinsic
spectrum.

The spectrum from each observation is shown in Figure \ref{xrtspectra} and the spectral fit information is included in Table~\ref{tab:xrt}. No significant spectral or flux variability is detected between observations. Moreover, the results are consistent with previous observations from \emph{Swift} XRT and other X-ray instruments in the 0.3--10 keV energy regime \cite{hess}.

\begin{table}[!h]

\begin{tabular}{cccccc}
\hline
\bf{Date} & \bf{Observation ID} & \bf{Exposure} & \bf{Log10[Flux] (2--10 keV)} & \bf{Index} & \bf{$\chi$$^{2}$/NDF}\\
 & & (seconds) & Log10 (erg cm$^{-2}$ s$^{-1}$) & \\
\hline
2014-06-17 & 00033319001 & 967 & -10.74$\pm0.06$ & 1.76$\pm$0.43 & 3.29/8 \\
2014-06-19 & 00033319002 & 769 & -10.66$\pm0.05$  &  1.98$\pm$0.37 & 3.25/7 \\
2014-06-21 & 00033319003 & 1156 & -10.52$\pm0.04$  &  1.65$\pm$0.28 & 11.93/16 \\
\hline
\end{tabular}
\caption{Summary of \emph{Swift}--XRT observations.}
\label{tab:xrt}
\vspace{1 mm}
\end{table}

\begin{figure}[!ht]
\includegraphics[width=100mm]{spectraXRT.pdf}
\centering
\caption{Flux spectra of HESS J1943+213 with \emph{Swift} XRT for the three separate observations.}
\label{xrtspectra}
\end{figure}

\section{Discussion and Outlook}

If HESS J1943+213 is a blazar, the agreement between VHE fluxes measured approximately five years apart with VERITAS and H.E.S.S.  and the lack of variability in the X-rays, including in the recent \emph{Swift} XRT result is surprising, but not unusual. Blazars are known to vary at all energies and at a wide range of timescales~\cite{boettcher}. The detection of variability at 1.5 GHz by ref.~\cite{gab2015} has weakened the lack of variability concern and is a strong piece of evidence in favor of the blazar scenario.

VERITAS can probe timescales that are a factor of four shorter than those available to H.E.S.S and therefore possesses the best VHE dataset available for searches of flux and spectral variability from this source. Planned application of advanced analysis techniques, such as use of boosted decision trees for gamma-hadron separation will provide an even higher sensitivity and allow for an additional factor of two improvement in the minimum variability timescale that can be tested.

The VERITAS-measured spectrum of HESS J1943+213 is in good agreement with the H.E.S.S. spectrum and extends down to 200 GeV. The soft VHE spectral index constituted one of the key pieces of evidence against the PWN hypothesis.

VERITAS will continue observations of HESS J1943+213 and will monitor the source for potential flux variability. In addition, the newly released PASS 8 Fermi-LAT data will allow for an improved detection and a spectrum of the source in GeV gamma rays, providing a significantly better handle on the gamma-ray peak of the HESS J1943+213 spectral energy distribution. Multi-wavelength studies of the source, including studies of its spectral energy distribution will be essential for definitively identifying the source class and investigating its properties.

\section{Acknowledgements}

This research is supported by grants from the U.S. Department of Energy Office of Science, the U.S. National Science Foundation and the Smithsonian Institution, and by NSERC in Canada. We acknowledge the excellent work of the technical support staff at the Fred Lawrence Whipple Observatory and at the collaborating institutions in the construction and operation of the instrument. The VERITAS Collaboration is grateful to Trevor Weekes for his seminal contributions and leadership in the field of VHE gamma-ray astrophysics, which made this study possible.