Exact Localisations of Feedback Sets

TL;DR

This paper introduces exact algorithms and structural insights for the feedback arc and vertex set problems, enabling more efficient solutions for transforming graphs into acyclic structures, with practical heuristics for improved performance.

Contribution

It presents polynomial-time computable subgraph structures, introduces the concept of essential minors and isolated cycles, and develops dynamic programming methods for exact solutions in special cases.

Findings

Polynomial-time computable subgraphs for cycle analysis

Exact solutions for resolvable graphs in polynomial time

Dynamic programming algorithms with exponential parameters for weighted problems

Abstract

The feedback arc (vertex) set problem, shortened FASP (FVSP), is to transform a given multi digraph $G=(V,E)$ into an acyclic graph by deleting as few arcs (vertices) as possible. Due to the results of Richard M. Karp in 1972 it is one of the classic NP-complete problems. An important contribution of this paper is that the subgraphs $G_{\mathrm{el}}(e)$, $G_{\mathrm{si}}(e)$ of all elementary cycles or simple cycles running through some arc $e \in E$, can be computed in $\mathcal{O}\big(|E|^2\big)$ and $\mathcal{O}(|E|^4)$, respectively. We use this fact and introduce the notion of the essential minor and isolated cycles, which yield a priori problem size reductions and in the special case of so called resolvable graphs an exact solution in $\mathcal{O}(|V||E|^3)$. We show that weighted versions of the FASP and FVSP possess a Bellman decomposition, which yields exact solutions using a dynamic programming technique in times $\mathcal{O}\big(2^{m}|E|^4\log(|V|)\big)$ and $\mathcal{O}\big(2^{n}\Delta(G)^4|V|^4\log(|E|)\big)$, where $m \leq |E|-|V| +1$, $n \leq (\Delta(G)-1)|V|-|E| +1$, respectively. The parameters $m,n$ can be computed in $\mathcal{O}(|E|^3)$, $\mathcal{O}(\Delta(G)^3|V|^3)$, respectively and denote the maximal dimension of the cycle space of all appearing meta graphs, decoding the intersection behavior of the cycles. Consequently, $m,n$ equal zero if all meta graphs are trees. Moreover, we deliver several heuristics and discuss how to control their variation from the optimum. Summarizing, the presented results allow us to suggest a strategy for an implementation of a fast and accurate FASP/FVSP-SOLVER.

Equations

$$X(G):=\underset{e\in E}{⨁}\mathbb{Z}e$$

$$\iota_{ij}=\left\{\begin{array}[]{rl}0&\,,\quad\text{if}\quad{e_{i}}^{+}\not=v_{j}\,\,\text{and}\,\,{e_{i}}^{-}\not=v_{j}\\ 1&\,,\quad\text{if}\quad{e_{i}}^{+}=v_{j}\,\,\text{and}\,\,{e_{i}}^{-}\not=v_{j}\\ -1&\,,\quad\text{if}\quad{e_{i}}^{+}\not=v_{j}\,\,\text{and}\,\,{e_{i}}^{-}=v_{j}\end{array}\right.\,.$$

$${\dim }_{\mathbb{Z}}\Lambda (G)={\dim }_{\mathbb{Z}}(\ker \mathcal{I}(G))=|E|-|V|+\#G,$$

$${\dim }_{{\mathbb{Z}}_2}{\Lambda }^0(G)={\dim }_{{\mathbb{Z}}_2}\ker ({\mathcal{I}}^0(G))=|E|-|V|+\#G,$$

$${\Omega }_{G,\omega }(\epsilon ):=\sum _{e\in \epsilon }\omega (e)$$

$$\Omega (G,\omega )\le \underset{e\in E}{\max }\omega (e)\cdot |E|/2.$$

$$F^{+}(e)$$

$$F^{-}(e)$$

$$x{\sim }_{u,v}y⟺x=y\text{or}x,y\in \{u,v\}.$$

$$\displaystyle E/_{\sim{u,v}}=\Big{\{}(p,q)\in V/_{\sim{u,v}}\times V/_{\sim{u,v}}\,\,\big{|}\,\,$$

$$G/e:=\big{(}V/_{\sim_{e^{+},e^{-}}},(E\setminus F(e))/_{\sim_{e^{+},e^{-}}}\big{)}.$$

$$\big{(}G/_{\Gamma},\omega/_{\Gamma})\big{)}:=\big{(}V/_{\Gamma},E/_{\Gamma},\omega/_{\Gamma})\big{)}$$

$$(G/_{\Phi},\omega/_{\Phi}):=\big{(}V/_{\Phi},E/_{\Phi},\omega/_{\Phi}\big{)}$$

$$(G_k,{\omega }_k):=(G_{k-1}^{\prime }{/}_{\Phi },{\omega }_{k-1}^{\prime }{/}_{\Phi })\text{ with }(G_{k-1}^{\prime },{\omega }_{k-1}^{\prime }):=(G_{k-1}{/}_{\Gamma },{\omega }_{k-1}{/}_{\Gamma }).$$

$$\epsilon \in \mathcal{S}(G,\omega )⟺{\epsilon }_C\in \mathcal{S}(C,\delta ).$$

$${\mathrm{O}}_{\mathrm{el}}(e):=$$

$${\mathrm{O}}_{\mathrm{si}}(e):=$$

$$P(u,f^{-})=\emptyset\quad\text{with respect to}\quad G\setminus\big{(}N_{E}^{-}(f^{-})\setminus\{f\}\big{)}\,.$$

$$I(c):=\left\{e\in\mathcal{E}(c)\,\,\big{|}\,\,c\cap c^{\prime}=\emptyset\,,\,\forall c^{\prime}\in\mathrm{O}_{\mathrm{el}}(G)\setminus\mathrm{O}_{\mathrm{el}}\big{(}F^{+}(e)\big{)}\right\}$$

$$O_{I}(e):=\left\{c\in\mathrm{O}_{\mathrm{el}}(G)\,\,\big{|}\,\,e\in I(c)\right\}\,,\quad O_{I}(\varepsilon)=\bigcup_{e\in\varepsilon}O_{I}(e)$$

$$\left<c^{\prime},c^{\prime\prime}\right>=0\,,\quad\text{for all}\,\,\,c^{\prime}\in O_{I}(e)/I(c),c^{\prime\prime}\in\big{(}O_{\mathrm{el}}(G/I(c)\big{)}\setminus\big{(}O_{I}(e)/I(c)\big{)}\,,$$

$$\Lambda\big{(}G/I(c)\big{)}=\mathrm{span}\big{(}O_{I}(e)/I(c)\big{)}\oplus\mathrm{span}\big{(}O_{\mathrm{el}}(G/I(c)\big{)}\setminus\big{(}O_{I}(e)/I(c)\big{)}\,.$$

$$\delta(e)=\Omega\big{(}G_{I}(e),\delta_{|G_{I}(e)}\big{)}\,,\quad G_{\mathrm{el}}(e)\not=G_{\mathrm{el}}(f)\,,\,\text{whenever}\,\,\,e\not=f\,.$$

$$G_{i+1}:=(V_{i},E_{i})=\mathcal{G}_{o}\big{(}E_{C_{i}}\setminus E_{i}^{*}\big{)}\,,\quad\omega:=\delta_{i|E_{i}^{\prime}}$$

$$\displaystyle\big{(}\omega(e_{3})-\Omega(H_{e_{3},d_{3}},\sigma)\big{)}-$$

$$(\omega (e_3)-3)-$$

$$\big{(}\omega(e_{3})+\omega(d_{1})\big{)}-\big{(}\omega(d_{3})+\omega(e_{1})\big{)}=8-7=1\,.$$

$$f\in E_e\text{and}{\mathrm{O}}_{\mathrm{el}}(f)\ne \varnothing \text{w.r.t.}G\setminus e.$$

$$V_k:=$$

$$E_k:=$$

$$(G_{\mathrm{el}}(h),\omega)\setminus\big{\{}\{f\}\cup\text{inner nodes of}\,\,M_{c}\setminus f\big{\}}$$

$${\sigma }_{G,M,e,f}:E⟶\mathbb{R}$$

Full text

††journal: Preprint version submitted to Theory of Computing Systems

1 Introduction

The feedback arc set problem, shortened FASP, is to delete as less as possible arcs of a graph such that the resulting subgraph is acyclic, i.e., it contains no directed cycle. Another equivalent formulation is to find a linear ordering of the vertices of the graph such that the number of back arcs is minimized. Therefore the problem is also known as maximum acyclic subgraph problem or linear ordering problem. For directed graphs this problem is one of the classic NP-complete problems Karp, (1972). The problem of deleting a smallest subset of vertices to result in an acyclic subgraph is known as feedback vertex set problem (FVSP). The FASP and FVSP are linear time reducible among each other, by keeping the relevant parameters fix as we will assert in section 2.1, alternatively see Even et al., (1998). Therefore algorithmic properties of the two problems are transferable. In particular, the FVSP is also NP-complete. Analogous problems for undirected graphs can be defined. As shown in Karp, (1972) the feedback vertex set problem remains NP-complete, while the feedback arc set problem can be solved efficiently by solving a maximum spanning tree problem. An excellent overview on feedback sets can be found in Bang-Jensen and Gutin, (2008). The problem of finding minimal transversals of directed cuts is closely related to the FASP, see Lucchesi and Younger, (1978). The FASP stays NP-complete for graphs where every node has an in-degree and out-degree of at most three or line digraphs even when every clique has at most size three Gavril, (1977). It is also NP-complete for tournament graphs Alon, (2006). However, there also exist graph classes possessing polynomial time algorithms, e.g., planar directed graphs or more general weakly acyclic digraphs Grötschel et al., (1985), and reducible flow graphs Ramachandran, (1988). The FASP or FVSP has a multitude of applications, e.g., retiming synchronous circuitry Leiserson and Saxe, (1991), circuit testing Kunzmann and Wunderlich, (1990), computational biology and neuroscience Ispolatov and Maslov, (2008), network analysis and operating systems Silberschatz et al., (2008).

1.1 Outline

In section 2 we provide the graph theoretical concepts, which are fundamental for this article. In section 3 we present our main results. The fact that the FASP/FVSP on multi-digraphs can be reduced to simple graphs is asserted in section 4 and the first a priori problem size reduction is deduced. In section 5 we construct an algorithm determing the induced subgraph of all cycles with one arc in common efficently. This knowledge is used in section 6 to introduce the concept of isolated cycles and resolvable graphs and providing an efficent solution of the FASP/FVSP on resolvable graphs. Moreover, the second a priori problem size reduction is given. Afterwards we concentrate on the main result of the article. Namely, that the FASP/FVSP possesses a Bellmann decomposition and present exact solutions using this fact to apply a dynamic programming technique in section 7. In section 8, we discuss how greedy approaches behave with respect to the problems and develope a strategy for a general FASP/FVSP-SOLVER. Finally, we discuss our results and other alternatives in section 9.

2 Preliminaries

Before we can introduce the FASP (FVSP) formally, some basic concepts of graphs and cycles need to be mentioned.

2.1 Graphs and Cycles

A multi-directed graph, or multi-digraph, $G=(V,E)$ consists of a set of vertices $V$ and a multi-set of arcs $E$ containing elements from $V\times V$. A directed graph or digraph is a multi-digraph with a simple arc set $E$, i.e., $E\subseteq V\times V$ and therefore every $e\in E$ occurs exactly once. A digraph is called simple if there are no loops, i.e., $E\cap\mathbb{D}(V\times V)=\emptyset$, where \mathbb{D}(V\times V)=\{(v,v)\in V\times V\,\,\big{|}\,\,v\in V\} denotes the diagonal.

If not stated otherwise throughout the article $G=(V,E)$ denotes a finite, connected, directed and loop-free multi-digraph and $G\setminus e$, $G\setminus v$ denote the graphs which occur by deleting the arc $e$ and possibly occurring isolated vertices or the vertex $v$ and all its adjacent arcs. For $\varepsilon\subseteq E$ and $\nu\subseteq V$ the graphs $G\setminus\varepsilon$, $G\setminus\nu$ are analogously defined. Moreover, $\mathcal{G}(\cdot)$, $\mathcal{E}(\cdot)$, $\mathcal{V}(\cdot)$ denote the induced graph, the set of all arcs, the set of all vertices which are inherited by a set or set system of graphs, arcs or vertices. With $\mathcal{P}(A)$ we denote the power set of a given set $A$.

For an arc $e=(u,v)\in E$ we denote ${e}^{+}=u$ as the tail and ${e}^{-}=v$ as the head of the arc. Two arcs $e$ and $f$ are called consecutive if ${e}^{-}={f}^{+}$ and are called connected if $\{{e}^{-},{e}^{+}\}\cap\{{f}^{-},{f}^{+}\}\not=\emptyset$. A directed path from a vertex $u$ to a vertex $v$ is a sequence of consecutive arcs where $u$ is the tail of the first arc and $v$ is the head of the last. A connected path from a vertex $u$ to a vertex $v$ is a sequence of connected arcs containing $u$ and $v$ as vertices. A digraph is connected if there is a connected path between every pair of its vertices. A weighted digraph $(G,\omega)$ or $(V,E,\omega)$, is a digraph with an additional weight function $\omega:E\longrightarrow\mathbb{R}$, which assigns a (usually positive) weight to each arc. For a given vertex $v\in V$ the sets N^{+}_{V}(v):=\left\{u\in V\,\,\big{|}\,\,(v,u)\in E\right\}\,,N^{-}_{V}(v):=\left\{u\in V\,\,\big{|}\,\,(u,v)\in E\right\}, N^{\pm}_{E}(v):=\left\{e\in E\,\,\big{|}\,\,e^{\pm}=v\right\} shall denote the set of all outgoing or incoming vertices or arcs of $v$ respectively. The indegree (respectively out degree) of a vertex $u$ is given by $\deg^{\pm}(u)=|N^{\pm}_{E}(u)|$ and the degree of a vertex is $\deg(u)=\deg^{-}(u)+\deg^{+}(u)$. $\Delta^{\pm}(G),\Delta(G)$ shall denote the maximal (in/out) degree, respectively.

A directed (connected) cycle of a digraph is a multiset of arcs $\{e_{0},\dots,e_{k}\}$ such that there is a permutation $\phi:\{0,\dots,k\}\longrightarrow\{0,\dots,k\}$ with $e_{\phi(i)}$ and $e_{\phi(i)+1\mod k+1}$ are consecutive (connected), for $1\leq i\leq k$. A loop is a cycle containing only a single arc. A cycle is simple if the set of contained arcs $\{e_{1},\dots,e_{k}\}$ is a simple set, i.e., it visits every arc, it contains, exactly once. A cycle is elementary if each vertex it contains is visited exactly once. We denote with $O(G)$ the set of all directed cycles and with $\mathrm{O}_{\mathrm{el}}(G)$, $\mathrm{O}_{\mathrm{si}}(G)$ the set of all elementary or simple cycles, respectively. Analogously, $\mathrm{O}_{\mathrm{el}}^{0}(G)$ and $\mathrm{O}_{\mathrm{si}}^{0}(G)$ shall denote the set of all connected (and not necessarily directed) elementary and simple cycles respectively. If not stated otherwise in the whole article a cycle is assumed to be directed and elementary. A feedback vertex set (FVS) of $G$ is a set $\nu\in\mathcal{P}(V)$ such that $G\setminus\nu$ is acyclic, i.e., $G\setminus\nu$ contains no directed cycle. A feedback arc set (FAS) of $G$ is a set $\varepsilon\in\mathcal{P}(E)$ such that $G\setminus\varepsilon$ is acyclic.

Definition 1 (line graph, natural hypergraph)

The directed line graph $L(G)=(V_{L},E_{L})$ of a digraph $G$ is a digraph where each vertex represents one of the arcs of $G$ and two vertices are connected by an arc if and only if the corresponding arcs are consecutive. In contrast the natural hypergraph $\mathcal{H}(G)=(\bar{V},\bar{E})$ of $G$ is constructed by identifying the arcs of $G$ with the vertices of $\mathcal{H}(G)$, i.e., $\bar{V}$ is a simple set of vertices such that $|V|=|E|$, where $|E|$ is counted with multiplicities. By fixing the identification $\bar{V}\cong E$ we introduce a directed hyperarc $h_{v}$ for every vertex $v\in V$ by requiring that head and tail coincide with all outgoing and ingoing arcs respectively, i.e., h_{v}=\big{(}N_{E}^{-}(v),N_{E}^{+}(v)\big{)}. Consequently, $\bar{E}\cong V$ and therefore every hyperarc can be labeled by its corresponding vertex. See Figure 1 for an example.

The directed, elementary cycles of the line graph $L(G)$ of $G$ are in 1 to 1 correspondence to the directed, simple cycles of $G$ while the directed simple cycles of $\mathcal{H}(G)$, i.e., directed cycles which run through every hyper arc exactly once are called Berge cycles, Berge, (1989). Note that if $G$ is a simple digraph, i.e., there are no multi arcs, then the set of all Berge cycles of $\mathcal{H}(G)$ are in 1 to 1 correspondence to the set of all elementary cycles of $G$. For a set of hyperarcs $\bar{\varepsilon}\subseteq\bar{E}$ of $\mathcal{H}(G)$ we denote with $\varepsilon\subseteq V$ the corresponding vertices in $G$. We summarize some facts in this regard.

Proposition 1

Let $G=(V,E)$ be a graph.

1. i)

The line graph $L(G)=(V_{L},E_{L})$ can be constructed in \mathcal{O}\big{(}|E|^{2}\big{)}. 2. ii)

The natural hypergraph $\mathcal{H}(G)=(\bar{V},\bar{E})$ can be constructed in \mathcal{O}\big{(}\Delta(G)|V|\big{)}. 3. iii)

$\nu\subseteq V_{L}$* is a FVS of $L(G)$ if and only if $\nu$ is a FAS of $G$.* 4. iv)

$\bar{\varepsilon}\subseteq\bar{E}$* is a FAS of $\mathcal{H}(G)$, with respect to the notion of Berge cycles, if and only if $\varepsilon$ is a FVS of $G$.*

Proof

Storing $G$ as a adjacency list and following the definitions immediately implies $i)$ and $ii)$. Since any two vertices $e,f$ of the line graph $L(G)$ are adjacent if and only if the corresponding arcs are consecutive in $G$ and any two arcs $h_{u},h_{v}$ of $\mathcal{H}(G)$ are consecutive, i.e., $h_{u}^{-}\cap h_{v}^{+}\not=\emptyset$, if and only if the corresponding vertices $u$ and $v$ are adjacent in $G$, $ii)$ and $iii)$ follow. ∎

Remark 1

Note that by introducing an additional arc $h^{*}_{v}$ between head and tail of every hyperarc $h_{v}=(N_{E}^{-}(v),N_{E}^{+}(v))$ of $\mathcal{H}(G)$, the natural hypergraph becomes a directed graph $G^{*}=(V^{*},E^{*})$ with $|V^{*}|=|E|+2|V|$, $|E^{*}|\leq\Delta(G)|V|+|V|$. The directed cycles of $G^{*}$ are in 1 to 1 correspondence to the Berge cycles of $\mathcal{H}(G)$ and a FAS $\varepsilon$ of $G^{*}$ is in 1 to 1 correspondence to a FAS $\bar{\varepsilon}$ of $\mathcal{H}(G)$ by identifying $\bar{\varepsilon}$ with the additional introduced arcs belonging to the hyperarcs in $\bar{\varepsilon}$ and identifying $\varepsilon$ with the hyperarcs corresponding to the bipartite graphs cutted by $\varepsilon$. If $\gamma:\bar{E}\longrightarrow\mathbb{R}+$ is an arc weight on $\mathcal{H}(G)$ then setting $\gamma^{*}(h)\equiv\gamma(h_{v})$ for all $h\in N_{E}^{-}(v)\cup N_{E}^{+}(v)\cup\{h^{*}_{v}\}$ yields the translated weight.

To give a more algebraic notion of cycles we consider

[TABLE]

the free $\mathbb{Z}$-module generated by $E$. If we choose coordinates, i.e., a numbering for $E$ and $V$ then we can identify $E$ with $\{e_{1},\dots,e_{|E|}\}$, $V$ with $\{v_{1},\dots,v_{|V|}\}$ and $X$ with $\mathbb{Z}^{|E|}$. In this case an element $x\in X(G)$ is a tuple $x=(x_{1},\dots,x_{|E|})$, which can be interpreted as a set of paths through $G$ where $x_{i}\in\mathbb{Z}$ indicates how often we pass the arc $e_{i}$ and the sign of $x_{i}$ determines in which direction this is done. We denote with $\mathcal{I}(G)$ the incidence matrix of $G$ with respect to these identifications, i.e., $\mathcal{I}(G)=(\iota_{ij})_{\genfrac{}{}{0.0pt}{}{1\leq i\leq|V|}{1\leq j\leq|E|}}$ with

[TABLE]

It is a well known fact, see for instance Biggs, (1993), that $x\in X(G)$ is a cycle of $G$ if and only if $\mathcal{I}x=0$, i.e., the submodule of all cycles of $G$ coincides with the set of homogeneous solutions $\Lambda(G)=\ker\mathcal{I}(G)$. In particular, this implies that $\Lambda(G)$ is a free $\mathbb{Z}$-module with

[TABLE]

where $\#G$ denotes the number of connected components of $G$ and therefore equals $1$ by our assumption on $G$. Note that the vector space $\Lambda^{0}(G):=\Lambda(G)/\mathbb{Z}_{2}^{|E|}$ can be understood as the cycle space of connected cycles, given by the kernel of the incidence matrix $\mathcal{I}^{0}(G)$ defined with respect to $\mathbb{Z}_{2}$ coefficients. In this case

[TABLE]

still holds, see again Biggs, (1993).

Remark 2

If $c\in\mathrm{O}_{\mathrm{el}}(G)$ then $1$ is the only non vanishing entry of $c$, i.e., $c\in\{0,1\}^{|E|}$. Moreover, no elementary cycle is subset of another elementary cycle. Indeed assume the opposite and consider two elementary cycles $c,c^{\prime}\in\mathrm{O}_{\mathrm{el}}(G)$ with $c\subseteq c^{\prime}$, then $\mathcal{I}(G)(c^{\prime}-c)=\mathcal{I}(G)c^{\prime}-\mathcal{I}(G)c=0-0=0$. Thus, $c^{\prime}-c\not=0$ is also a positive oriented cycle and therefore we have $c^{\prime}=(c^{\prime}-c)+c\not\in\mathrm{O}_{\mathrm{el}}(G)$. A contradiction!

Not that every simple cycle is given by the union of arc disjoint elementary cycles. The following example illustrates this fact.

Example 1

Consider the graph $G$ in Figure 2. Then one observes that the cycle $c=\{e,f,g,h,i\}$ is a simple, non-elementary cycle while the cycles $\{e,f,i\}$, $\{f,h\}$ and $\{h,g\}$ are elementary cycles.

2.2 The Feedback Arc Set Problem (FASP)

Now we have all ingredients to give a formal definition of the FASP.

Problem 1

Let $G=(V,E)$ be a finite, connected, directed, and loop-free graph, $\omega:E\longrightarrow\mathbb{N}^{+}$ be a weight function. Then the weighted FASP is to find a set of arcs $\varepsilon\in\mathcal{P}(E)$ such that $G\setminus\varepsilon$ is acyclic, i.e., \mathrm{O}_{\mathrm{el}}\big{(}G\setminus\varepsilon\big{)}=\emptyset and

[TABLE]

is minimized. We denote the set of solutions of this problem with $\mathcal{S}(G,\omega)$ and denote with $\Omega(G,\omega):=\Omega_{G,\omega}(\varepsilon)$, $\varepsilon\in\mathcal{S}(G,\omega)$ the optimal weight or feedback length. If $\omega$ is constant, e.g., equal to $1$, then the problem coincides with the unweighted minimal FASP.

Remark 3

The condition on $G$ to be loop-free is not an essential restriction. This is because every loop is contained in any solution of the minimal FASP.

Remark 4

Note that, every cycle $c\in O(G)$ can be generated by elementary cycles $c_{1},\dots,c_{n}\in\mathrm{O}_{\mathrm{el}}(G)$ using only non-negative coefficients. Thus, if $\varepsilon$ is a solution of Problem 1 then certainly O\big{(}G\setminus\varepsilon\big{)}=\emptyset, which implies that our notion of acyclic graphs is consistent for the problem.

Remark 5

Note that for given graph $G=(V,E)$ the smallest subgraph $G_{o}\subseteq G$, which contains all cycles of $G$, i.e., $G^{\prime}=G_{o}$ whenever $G^{\prime}\subseteq G_{o}$ and $\mathrm{O}_{\mathrm{el}}(G)=\mathrm{O}_{\mathrm{el}}(G^{\prime})$ can be constructed in $\mathcal{O}(|E|^{2})$. Indeed the arc set of $G_{o}$ is constructed from $E$ by removing arcs $(u,v)$ if there is no directed path from $v$ to $u$. For every arc this can be done by depth-first search in linear time if $G$ is stored in an adjacency list. Since removing arcs does not generate new cycles it suffices to check every arc once yielding the estimated runtime performance. Certainly, a solution for $G_{o}$ is a solution for $G$. We shortly denote with $\mathcal{G}_{o}(G):=G_{o}$ and with $\mathcal{G}_{o}^{0}(G):=G^{0}_{o}$ the analogous graph appearing by considering connected cycles instead of directed ones. In particular, by the argumentation above, an elementary or directed cycle $c\in\mathrm{O}_{\mathrm{el}}(G)$, $c^{\prime}\in\mathrm{O}_{\mathrm{si}}(G)$ can be found in $\mathcal{O}(|E|^{2})$ or no cycle exist.

Remark 6

If $G=(V,E)$ is a simple graph then denoting with $\lfloor\cdot\rfloor$ the Gauss bracket we observe that at most $\lfloor|E|/2\rfloor$ arcs have to be deleted to obtain a graph where no connected path of length $2$ exists anymore. In particular, the graph is acyclic in this case and therefore

[TABLE]

See also Saab, (2001).

2.3 The Feedback Vertex Problem (FVSP)

Let $G=(V,E)$ be given and $\gamma:V\longrightarrow\mathbb{R}^{+}$ be a vertex weight. The feedback vertex set problem (FVSP) on $(G,\gamma)$ is obtained by replacing the role of arcs by vertices in Problem 1. In regard of Proposition 1, by introducing the hyperarc weight $\bar{\omega}(h_{v}):=\gamma(v)$, we realize that the FVSP is equivalent to the the FASP on the natural hypergraph $\mathcal{H}(G)$ of $G$, with respect to the notion of Berge cycles. Already in Remark 1 we mentioned how to treat this case. Vice versa the FASP on an arc weighted graph $(G,\omega)$ is equivalent to the FVSP on the line graph $L(G)$ of $G$ by introducing the vertex weight $\gamma(v)=\omega(v)$, $v\in V_{L}=E$. Since the described transformations can be done efficently, see Proposition 1, an efficent solution of the FASP or FVSP for an arc and vertex weighted instance $(G,\omega,\gamma)$ yields an efficent solution of the FVSP or FASP for the transformed instances and vice versa. In particular, by summarizing some already known results we obtain:

Theorem 1

The unweighted FASP and FVSP are APX complete.

Proof

Since there is an $L$-reduction of the Minimum Vertex Cover Problem, which is APX complete due to Dinur and Safra, (2004), to the FVSP, see Karp, (1972), the FVSP is APX complete. Proposition 1 shows that the unweighted FASP on $G$ is equivalent to the unweighted FVSP on $L(G)$. Thus, the feedback length of any solution remains unchanged yielding an $L$-reduction of the FASP to the FVSP implying the claim.

We expect that the theorem above still holds for the weighted versions. However, due to the observations made so far, we will focus our studies on the FASP to increase the understanding of the localisation of feedback sets.

3 Main Results

Though the FASP and FVSP are equivalent problems in graph theory and computer sciences the only exact solutions of the FVSP known to us is are the algorithms of Razgon, (2007) with complexity $\mathcal{O}(1,9977^{|V|}|V|^{\mathcal{O}(1)})$ and Chen et al., (2008) requiring \mathcal{O}\big{(}|E|^{4}4^{\Omega}\Omega^{3}\Omega!\big{)}, where $\Omega$ denotes the feedback length. A detailed comparison to our approach is given in Section 9. For now we just mention our results:

Theorem A

Let $(G,\omega)$ be a graph. Then there is an algorithm testing whether $(G,\omega)$ is resolvable and determing a solution of the weighted FASP on $G$ in case of resolvability in $\mathcal{O}(|V||E|^{3})$.

Though there are infinitely many resolvable graphs not all graphs are resolvable. However, if the graph $(G,\omega)$ is not resolvable, we still can find an exact solution:

Theorem B

Let $(G,\omega)$ be a graph then there is a parameter $m\in\mathbb{N}$, $m\leq\dim_{\mathbb{Z}_{2}}\Lambda^{0}(G)=|E|-|V|+1$, which can be determined in $\mathcal{O}(|E|^{3})$ and an algorithm CUT with run time \mathcal{O}\big{(}2^{m}|E|^{4}\log(|V|)\big{)} solving the weighted FASP.

Due to Proposition 1 the analogous results with respect to the FVSP hold. In particular, we call a vertex weighted graph $(G,\nu)$ resolvable iff its natural hypergraph is resolvable, see section 2.3 again. If we replace every hyperarc of $\mathcal{H}(G)=(\bar{V},\bar{E})$ with its corresponding bipartite graph then by following Remark 1 we have $|V^{*}|=|E|+2|V|$, $|E^{*}|\leq(\Delta(G)+1)|V|$ for the resulting graph $G^{*}=(V^{*},E^{*})$. The translated results therefore become:

Theorem C

Let $(G,\nu)$, $\nu:V\longrightarrow\mathbb{R}^{+}$ be a vertex weighted graph. Then there is an algorithm testing whether $(G,\nu)$ is resolvable and determing a solution of the weighted FVSP on $G$ in case of resolvability in $\mathcal{O}(\Delta(G)^{3}|V|^{3}|E|)$.

In case of non-resolvability we have:

Theorem D

Let $(G,\nu)$ be a graph then there is a parameter $m\in\mathbb{N}$, $m\leq(\Delta(G)-1)|V|-|E|+1$, which can be determined in $\mathcal{O}(\Delta(G)^{3}|V|^{3})$ and an algorithm CUT with run time \mathcal{O}\big{(}2^{m}\Delta(G)^{4}|V|^{4}\log(|E|)\big{)} solving the weighted FVSP.

Note that there are infinitely many instances where $m=0$ and $m\leq|V|$ on homogenous graphs, i.e., if $\Delta(G)\leq(|E|-1)/|V|+2$. Thus, by computing the bound or directly $m$ we can decide whether the algorithm of Razgon, (2007) or our approach will be faster for a given instance and choose the better alternative. Moreover, the feedback length $\Omega$ can not assumed to be constantly bounded. Thus, the algorithm of Chen et al., (2008) actually possesses a complexity of \mathcal{O}\big{(}4^{|V|}|V|^{3}|E|^{4}|V|!\big{)}, which is much slower than our approach. Finally, we want to mention that all theorems are based on the following crucial fact:

Theorem E

Let $G=(V,E)$ be a graph and $e\in E$. Then there exist algorithms which compute:

1. i)

The subgraph $G_{\mathrm{el}}(e)\subseteq G$ induced by all elementary cycles $c\in\mathrm{O}_{\mathrm{el}}(G)$ with $e\in\mathcal{E}(c)$ in $\mathcal{O}(|E|^{2})$. 2. ii)

The subgraph $G_{\mathrm{si}}(e)\subseteq G$ induced by all simple cycles $c\in\mathrm{O}_{\mathrm{si}}(G)$ with $e\in\mathcal{E}(c)$ in $\mathcal{O}(|E|^{4})$.

A proof and a more precise version of the statement is given in Theorem 2.

4 The Essential Minor

In this section we introduce the notion of the essential minor $(C,\delta)$ of given graph $(G,\omega)$, which is a simple, weighted digraph that decodes the topological structure of $G$ in a compact way and is therefore very helpful. Even though there are some crucial differences we want to mention that similar concepts were already introduced in Berge, (2001).

Definition 2 (parallel arcs)

Let $G$ be a graph and $e=(u,v)\in E$ then we denote with

[TABLE]

the sets of all parallel and anti-parallel arcs and set $F(e)=F^{+}(e)\cup F^{-}(e)$.

We recall that for a given set $A$ and an equivalence relation $\sim$ on $A\times A$ the quotient $A/_{\sim}$ is given by the set of all equivalence classes [a]_{\sim}=\left\{x\in A\,\,\big{|}\,\,x\sim a\right\}.

Definition 3 (contracted graph)

Let $G=(V,E)$ be a graph and $u,v\in V$. An equivalence relation $\sim_{u,v}$ on $V$ is defined by

[TABLE]

The equivalence class of $x\in V$ is denoted by $[x]_{\sim{u,v}}$ and $V/_{\sim_{u,v}}$ gives the quotient of $V$ with respect to $\sim_{u,v}$. A multiset $E/_{\sim{u,v}}$ is defined by

[TABLE]

The contracted graph of $G$ with respect to $e\in E$ is defined as the topological minor

[TABLE]

If $e,f\in E$ then one can check easily that by identifying $e,f$ with their images in $G/e$, $G/f$ respectively we have $(G/e)/f=(G/f)/e$. Thus, the definition does not depend on the order of the contracted edges. Hence, if $G^{\prime}\subseteq G$ is a subgraph then $G/G^{\prime}:=G/\mathcal{E}(G^{\prime})$ can be defined by contracting $\mathcal{E}(G^{\prime})$ in arbitrary order.

Definition 4 (essential minor)

Let $G=(V,E,\omega)$ be a positively weighted graph. The equivalence relation $\sim_{\Gamma}$ on $E$ is defined by $e\sim_{\Gamma}f$ if and only if $e=f$ or there exists a directed, branch-point-free path $(w_{0},w_{1},\dots,w_{n},w_{n+1})$ with $\{(w_{0},w_{1}),(w_{n},w_{n+1})\}=\{e,f\}$ and $n\in\mathbb{N}^{+}$, i.e., for $i=1,\dots,n$ it holds that $\deg^{-}(w_{i})=\deg^{+}(w_{i})=1$. We represent an equivalence class $[e]_{\sim\Gamma}$ by an arc $(u,v)$, where $u$ and $v$ coincide with the start and endpoint of the longest directed, branch-point-free path running through $e$. The positively weighted graph

[TABLE]

is defined by

1. i)

V/_{\Gamma}:=\{v\in V\,\,\big{|}\,\,\deg^{-}(v)>1\vee\deg^{+}(v)>1\}, 2. ii)

E/_{\Gamma}:=\left\{(u,v)\in V/_{\Gamma}\times V/_{\Gamma}\,\,\big{|}\,\,\exists\,\,[e]/_{\sim\Gamma}\in E_{/\sim_{\Gamma}}\,\,\text{represented by }(u,v)\right\} 3. iii)

$\omega/_{\Gamma}:E/_{\Gamma}\longrightarrow\mathbb{N}^{+}$ with $\omega/_{\Gamma}(e):=\min_{e^{\prime}\in[e]/_{\sim\Gamma}}\omega(e^{\prime})$.

Let $\sim_{\Phi}$ be an equivalence relation on $E$ with $e\sim_{\Phi}f\Longleftrightarrow e,f\in F^{+}(e).$ For $(G,\omega)$ the positively weighted graph

[TABLE]

is defined by

1. i)

$V/_{\Phi}:=V$. 2. ii)

$E/_{\Phi}:=E/_{\sim_{\Phi}}$, where we identify each equivalence class $[e]_{\sim\Phi}$ with an arbitrary representative of $F^{+}(e)$. 3. iii)

$\omega/_{\Phi}:E/_{\Phi}\longrightarrow\mathbb{N}^{+}\,,\,\,\omega/_{\Phi}(e):=\sum_{e^{\prime}\in[e]_{\sim\Phi}}\omega(e^{\prime})\,.$

Starting with $G_{0}:=G$ for $k>0$ we define

[TABLE]

The weighted graph $(G_{K},\omega_{K})$ with $G_{K}=G_{K+1}$, $K\geq 0$, is called the essential minor of $G$ and is denoted by $(C,\delta):=(V_{C},E_{C},\delta)$.

Example 2

In Figure 3 the construction of the essential minor is illustrated. Furthermore, for a graph $G$ with $D$ diamonds connected in a cycle as in the example we obtain $\deg(G)=3$ and \big{|}\mathrm{O}_{\mathrm{el}}(G)\big{|}=2^{D}=2^{|E|/5}=2^{|V|/4}. In contrast the essential minor $C$ of such a graph satisfies $|\mathrm{O}_{\mathrm{el}}(C)|=1$. Hence, even though the number of cycles in $G$ increases exponentially in $|E|$ and $|V|$ by adding further diamonds, the number of cycles in $C$ remains constant equal to 1 while the weight $\xi$ decodes the number of cycles of the original graph $G$.

Since the following results are quite canonically their simple but technical proofs are given in Appendix A.

Proposition 2

Let $G=(V,E,\omega)$ be a positively weighted graph with essential minor $(C,\omega_{C})$ and let $\varepsilon\in\mathcal{P}(E)$ and $\varepsilon_{C}$ be the image of $\varepsilon$ in $(C,\delta)$. Then

[TABLE]

In particular $\Omega(G,\omega)=\Omega(C,\omega_{C})$.

Proposition 2 states that solving the FASP for the essential minor is equivalent to solving the FASP on the original graph. Even though it is possible that $(C,\delta)=(G,\omega)$ the construction might yields an a priori problem size reduction in many cases as in Example 2.

Proposition 3

Let $G=(V,E,\omega)$ be a finite, connected, directed, weighted multigraph then we can construct $(C,\delta)$ in time $\mathcal{O}(|E|^{2})$. Furthermore, there is an algorithm with run time $\mathcal{O}(|E|^{2})$ which constructs a solution $\varepsilon\in S(G,\omega)$ given a solution $\varepsilon_{C}\in\mathcal{S}(C,\delta)$.

Remark 7

Since for given $\varepsilon_{C}\in\mathcal{S}(C,\delta)$ the construction of some $\varepsilon_{G}\in\mathcal{F}(\varepsilon_{C})$ is easy to compute (Proposition 3) Proposition 2 states that it suffices to solve a the weighted FASP for the essential minor instead of the original graph. As a consequence multigraphs do not need to be considered and the number of elementary cycles of the essential minor can be drastically reduced, see for instance Example 2.

5 Subgraphs of Elementary and Simple Cycles

There are several approaches for generating the set $\mathrm{O}_{\mathrm{el}}(G)$ of all elementary cycles of a graph, see Mateti and Deo, (1976) for an overview. Since there can be an exponential number of cycles in a graph, generating algorithms have an exponential worst case run time. The best algorithms available today are the ones of Tarjan, (1973) and Johnson, (1975) solving the problem in \mathcal{O}\big{(}|\mathrm{O}_{\mathrm{el}}(G)|(|V|+|E|)\big{)}. Of course counting all cycles might be less expansive than generating them. However, by reducing to the Hamiltonian cycle problem, see for instance Arora and Barak, (2009), counting all cycles is a NP-hard problem. For our concerns, and supposedly in many other situations, the generation of all cycles is not necessary, but the knowledge of the arc set of all cycles including a common arc suffices. In the following an algorithm for determing the smallest subgraph $G_{\mathrm{el}}(e)\subseteq G$ which contains all elementary cycles that include the arc $e$ is given.

Definition 5

Let $G=(V,E)$ be a graph and $u,v\in V$ we denote with with $P(u,v),P_{\mathrm{el}}(u,v)$, $P_{\mathrm{si}}(u,v)$ the set of all directed, elementary or simple paths from $u$ to $v$ respectively. For an arc $e\in E$ we let

[TABLE]

be the set of all elementary and simple cycles running through $e$. If $\varepsilon\in\mathcal{P}(E)$ then we set $\mathrm{O}_{\mathrm{el}}(\varepsilon):=\cup_{e\in\varepsilon}\mathrm{O}_{\mathrm{el}}(e)$, $\mathrm{O}_{\mathrm{si}}(\varepsilon):=\cup_{e\in\varepsilon}\mathrm{O}_{\mathrm{si}}(e)$. Moreover, we denote with G(u,v):=\mathcal{G}(P(u,v)),G_{\mathrm{el}}(u,v):=\mathcal{G}\big{(}P_{\mathrm{el}}(u,v)\big{)}, G_{\mathrm{si}}(u,v):=\mathcal{G}\big{(}P_{\mathrm{si}}(u,v)\big{)} the by the corresponding paths induced graphs and with $G(e):=\mathcal{G}(P(u,v))$, G_{\mathrm{el}}(e):=\mathcal{G}\big{(}P_{\mathrm{el}}(e^{-},e^{+})\cup\{e\}\big{)}, G_{\mathrm{si}}(e):=\mathcal{G}\big{(}P_{\mathrm{si}}(e^{-},e^{+})\cup\{e\}\big{)} by the corresponding cycles induced graphs. Moreover, $P^{0}(u,v),\mathrm{P}_{\mathrm{el}}^{0}(u,v)$, $\mathrm{P}_{\mathrm{si}}^{0}(u,v)$, $O(e)$, $O_{\mathrm{el}}^{0}(e)$, $O_{\mathrm{si}}^{0}(u,v)$, $G^{0}(e)$, $G_{\mathrm{el}}^{0}(e)$, $G_{\mathrm{si}}^{0}(e)$ shall denote the connected (and not necessarily directed) analagons of the introduced sets and graphs.

Note that $P_{\mathrm{el}}(u,v)\subseteq P_{\mathrm{si}}(u,v)$ and therefore $\mathrm{O}_{\mathrm{el}}(e)\subseteq\mathrm{O}_{\mathrm{si}}(e)$. Moreover, the graphs $G(u,v),G(e)$ can be determined in $\mathcal{O}(|E|^{2})$ by applying a depht first search technique similar to Remark 5. In the other cases we observe:

Theorem 2

Let $G=(V,E)$ be a graph and $u,v\in V$, $e\in E$. Then there exist algorithms which compute:

1. i)

The graphs $G_{\mathrm{el}}(u,v)$,$G_{\mathrm{el}}^{0}(u,v)$, $G_{\mathrm{el}}(e)$,$G_{\mathrm{el}}^{0}(e)$* in $\mathcal{O}(|E|^{2})$.* 2. ii)

The graphs $G_{\mathrm{si}}(u,v)$, $G_{\mathrm{si}}^{0}(u,v)$, $G_{\mathrm{si}}(e)$, $G_{\mathrm{si}}^{0}(e)$ in $\mathcal{O}(|E|^{4})$.

Proof

Let $p\in P_{\mathrm{el}}(u,v)$ then no vertex $w\in\mathcal{V}(p)$ is passed twice of $p$. Thus, for every $f\in\mathcal{E}(p)$ there is a path $q\in P(u,f^{-})$ with respect to G\setminus\big{(}N_{E}^{-}(f^{-})\setminus\{f\}\big{)}. Vice versa if $p\in P(u,v)$ is such that $\mathcal{E}(p)\not\subseteq G_{\mathrm{el}}(u,v)$ then has to be a vertex $w\in\mathcal{V}(p)$, which is passed at least twice by $p$ implying that there is $f\in\mathcal{E}(p)$ such that

[TABLE]

Thus, by setting G^{\prime}:=G\setminus\left\{f\in E\,\,\big{|}\,\,f\,\,\text{fulfills \eqref{DELf}}\right\}, every path $p\in P(u,v)$ with $\mathcal{G}(p)\not\subseteq G_{\mathrm{el}}(u,v)$ is interrupted in $G^{\prime}$. Hence, $G^{\prime}(u,v)=\mathcal{G}(P(u,v))$ with respect to $G^{\prime}$ coincides with $G_{\mathrm{el}}(u,v)$. Algorithm 1 formalizes this procedure and runs in $\mathcal{O}(|E|^{2})$ if $G$ is stored in an adjacency list, enabling us to test whether $P(u,v)=\emptyset$ in $\mathcal{O}(|E|)$. The other cases of $i)$ can now be solved by replacing $u,v$ with $e$ and directed paths or cycles with connected ones.

To show $ii)$ we add two arcs $u^{*}=(x,u)$, $v^{*}=(v,y)$, $x,y\not\in V$ denote with $G^{*}$ the resulting graph and consider the line graph $L(G^{*})=(V_{L}^{*},E_{L}^{*})$. We recall that $|V_{L}|=|E|$, $|E_{L}|\leq|E^{2}|$ and apply the fact that the elementary paths of $L(G)$ are in $1$ to $1$ correspondence to the simple paths of $G$ and therefore \mathcal{V}\big{(}G_{\mathrm{el}}(u^{*},v^{*})\big{)}\setminus\{u^{*},v^{*}\}=\mathcal{E}\big{(}G_{\mathrm{si}}(u,v)\big{)}. Hence $ii)$ follows analogue to $i)$. ∎

Remark 8

Note that if $(C,\delta)$ is the essential minor of $(G,\omega)$. Then the treatment of “parallel” paths is avoided by the essential minor construction. Thus, we expect that if $C$ is significant smaller than $G$ the run time performance will increase drastically.

6 Isolated Cycles

Of course the question arises whether a solution of the FASP on $G_{\mathrm{el}}(e)$ can be determined independently of the rest of the graph. The notion of isolated cycles is our starting point of investigations in this manner and as it will turn out it is a very helpful concept of answering this question.

We recall that a Min-s-t-Cut with source $s=u$ and sink $t=v$ is given by a set $\varepsilon\subseteq E$ such that $P(u,v)=\emptyset$ in $G\setminus\varepsilon$ and $\Omega_{G,\omega}(\varepsilon)=\sum_{e\in\varepsilon}\omega(e)$ is minimized.

Lemma 1

Let $(G,\omega)$ be a weighted graph and $e\in E$. Then there is an algorithm, which determines a solution \varepsilon\in\mathcal{S}\big{(}G_{\mathrm{el}}(e),\omega\big{)} of the FASP on $(G_{\mathrm{el}}(e),\omega)$ in \mathcal{O}\big{(}|V||E|\log(|V|)\big{)}, where we slightly abused notion by still denoting $\omega$ for the the restriction of $\omega$ to $G_{\mathrm{el}}(e)$.

Proof

Observe that by interpreting $\omega$ as a capacity function on $G_{\mathrm{el}}(e)$ a solution of the FASP on $\varepsilon$ is given by $\{e\}$ or a Min-s-t-Cut $\varepsilon$ with source $s=e^{-}$ and sink $t=e^{+}$. The option with the smaller weight is chosen. Due to the famous Min-Cut-Max-Flow Theorem a Min-s-t-Cut can be determined by solving a Max-Flow problem with respect to $\omega$ and $s=e^{-}$, $t=e^{+}$. The algorithm of Dinic, (1970) solves the Max-Flow problem for arbitrary weights in time $\mathcal{O}(|V|^{2}|E|)$ and can be speedend up to \mathcal{O}\big{(}|V||E|\log(|V|)\big{)}, Dinitz, (2006) by using the data structure of dynamic trees. Thus the statement is proven. ∎

Remark 9

Note that if $(C,\delta)$ is the essential minor of $(G,\omega)$ then the absence of “parallel” paths might speeds up the time required to determine a Min-s-t-Cut drastically. Moreover, the Max-Flow-Problem is very well understood, yielding many alternatives to the algorithm of Dinic, (1970) and providing faster solutions in special cases, see Dinitz, (2006) and B. Korte, (1981) for an overview.

Definition 6

Let $G$ be a graph and $c\in\mathrm{O}_{\mathrm{el}}(G)$ then we denote with

[TABLE]

the set of isolating arcs of $c$, i.e., if $e\in I(c)$ then $c$ has empty intersection with every cycle $c^{\prime}$ that does not contain $e$ or an parallel arc of $e$. For an arc $e\in E$ or set of arcs $\varepsilon\subseteq E$ we set

[TABLE]

and G_{I}(e):=\mathcal{G}\big{(}O_{I}(F^{+}(e))\big{)}.

Remark 10

Note, that the sets of isolated cycles possess a flat hierarchy in the following sense. If $e,f\in E$, $e\neq f$, with $O_{I}(e)\not=O_{I}(f)$ then $O_{I}(e)\cap O_{I}(f)=\emptyset$. If vice versa $O_{I}(e)=O_{I}(f)$ then by definition we obtain $\mathrm{O}_{\mathrm{el}}(F^{+}(e))=\mathrm{O}_{\mathrm{el}}(F^{+}(f))$.

Remark 11

Let $c\in O_{I}(G)$ be an isolated cycle and $I(c)$ the set of all isolating arcs of $c$. If we contract $I(c)$ then the resulting graph $G/I(c)$ fulfills

[TABLE]

where $\left<\cdot,\cdot\right>$ denotes the standard scalar product on $\mathbb{R}^{n}$, n=\big{|}E/I(c)\big{|}. Thus, by detecting isolated cycles we obtain an orthogonal splitting

[TABLE]

Such a splitting is certainly helpful whenever one wants to find a basis of $\Lambda(G)$, e.g., a minimal cycle basis of $\mathrm{span}\left(O_{I}(e)\right)$ can be extended to a minimal basis of $\Lambda(G)$.

Consider an isolating arc $e$ of a graph $(G,\omega)$ or its essential minor $(C,\delta)$. The isolated cycles $O_{I}(e)$ running through $e$ can be cut either by removing the arc set $\varepsilon_{0}=F^{+}(e)$ or another feedback set $\varepsilon_{1}$ of $G_{I}(e)$. By definition the arc set $\varepsilon_{0}$ cuts at least the cycles in $O_{I}(\varepsilon_{0})$. By Remark 10 the feedback set $\varepsilon_{1}$ cuts only the isolated cycles $O_{I}(\varepsilon_{0})$ or is given by $\varepsilon_{1}=F^{+}(f)$ of another isolating arc $f\in O_{I}(e)$ with $\mathrm{O}_{\mathrm{el}}(\varepsilon_{0})=\mathrm{O}_{\mathrm{el}}(\varepsilon_{1})$. Thus, if the weight of $\varepsilon_{0}$ equals the weight of a solution of the FASP on $G_{I}(e)$ then there is a solution $\varepsilon$ of the FASP on $G$ with $\varepsilon_{0}\subseteq\varepsilon$. The following definition reflects this idea more formally.

Definition 7

Let $(G,\omega)$ be a graph. We define a maximal list of graphs $(G_{0},\omega_{0}),\dots,(G_{k},\omega_{k})$ with $(G_{0},\omega_{0})=(G,\omega)$, $(G_{i},\omega_{i})\neq(G_{i+1},\omega_{i+1})\,,\,\forall\,i\in[1:k-1]$ as follows. Let $(C_{i},\delta_{i})$, with $C_{i}=(V_{C_{i}},E_{C_{i}})$, be the essential minor of $(G_{i},\omega_{i})$ and $E^{*}_{i}\subseteq I(C_{i})$ be a maximal subset of pairwise different isolating arcs of $C_{i}$ such that $\forall e\in E^{*}_{i}:$

[TABLE]

Then the weighted graph $(G_{i+1},\omega_{i+1})$ is given by

[TABLE]

where $\delta_{i+1}:=\delta_{i|E_{i+1}}$ denotes the restriction of $\delta_{i}$ to $E_{i+1}$. If $G_{k+1}=G_{k}$ for some $k\in\mathbb{N}$ then (S,\tau):=\big{(}(V_{S},E_{S}),\tau\big{)}=(C_{k},\delta_{k}), $k\in\mathbb{N}$ is called the resolved graph of $G$, which we shortly denote with (S,\tau)=\big{(}S(G),\tau(\omega)\big{)}. A graph $(G,\omega)$ is called resolvable if and only if $S=\emptyset$.

Example 3

Note that $(C,\delta)$ in Figure 4 is resolvable, while $(G,\omega)$ in Figure 6 is not resolvable, but becomes resolvable for uniform weight $\omega\equiv 1$.

The construction has an immediate consequence.

Theorem 3

Let $(G,\omega)$ be a graph. Then there is an algorithm testing whether $(G,\omega)$ is resolvable and determing a solution of the weighted FASP on $G$ in case of resolvability in $\mathcal{O}(|V||E|^{3})$.

Proof

Due to Proposition 3 the construction of the essential minor $(C_{i},\delta_{i})$ can be achieved in $\mathcal{O}(|E|^{2})$ for every $1\leq i\leq k$. Since checking whether $\mathrm{O}_{\mathrm{el}}(f)\not=\emptyset$ can be done in \mathcal{O}\big{(}|E|\big{)} by storing $G$ in an adjacency list and using depht first search to figure out whether $P(f^{+},f^{-})\not=\emptyset$ the Algorithm 2 computes the set $G_{I}(e)$ in $\mathcal{O}(|E|^{2})$ and therefore computing $G_{I}(e)$ for all arcs requires at most $\mathcal{O}(|E|^{3})$ computation steps. Furthermore, a solution of the FASP on $G_{I}(e)$ can be computed due to Lemma 1 in \mathcal{O}\big{(}|E|^{2}\big{)}. Due to the fact that during the construction of $(S,\tau)$ no parallel arcs appear, we have to recompute the isolated cycles at most $|V|$ times. Thus, $(S,\tau)$ can be determined in \mathcal{O}\big{(}|V||E|^{3}\big{)}. Furthermore, we can use the backtracking procedure of Proposition 3 to compute a solution of the FASP in $\mathcal{O}(|V||E|^{2})$ once $(S,\tau)$ is known. ∎

Observe that Theorem 3 was already stated in section 3 as Theorem A. However, the result leads to the question : What are fast (linear, quadratic time) checkable conditions a graph $(G,\omega)$ has to satisfy to be resolvable. Though, we can easily construct resolvable graphs as $(C,\delta)$ in Figure 4 or modified versions of $(C,\delta)$ by adding additional isolated cycles a characterization of resolvable graphs is still open. A better understanding of the non-resolvable graphs might help to solve that problem. In order to investigate these graphs the class of graphs without isolated cycles at all, seems to be interesting. Therefore, the next result might be a good starting point for further studies.

Proposition 4

The directed clique $D_{3}=(V,E)$ in Figure 5 is the smallest graph with $\mathrm{O}_{\mathrm{el}}(G)\not=\emptyset$ and $O_{I}(G)=\emptyset$, i.e. any non-isomorphic graph $G^{\prime}=(V^{\prime},E^{\prime})$ with $O(G^{\prime})\not=\emptyset$ and $O_{I}(G^{\prime})=\emptyset$ satisfies $|V^{\prime}|+|E^{\prime}|>|V|+|E|$ .

Proof

If $G^{\prime}$ is a graph with $O(G^{\prime})\not=\emptyset$, $O_{I}(G^{\prime})=\emptyset$ then $G^{\prime}$ possesses at least three vertices and \big{|}O(G^{\prime})\big{|}\geq 3 has to hold. We claim that there are at least three linear independent cycles. Indeed if $d_{1},d_{2},d_{3}\in O(G^{\prime})$ with $\lambda d_{1}+\mu d_{2}+\eta d_{3}=0$, $\lambda,\mu,\eta\in\mathbb{Z}$ then due to Remark 2 we know that $d_{1},d_{2},d_{3}\in\{0,1\}^{|E|}$ and no elementary cycle is subset of another. So w.l.o.g. we can assume that $d_{1}=d_{2}+d_{3}$, which contradicts that all cycles are elementary. Thus, $\dim_{\mathbb{Z}}\Lambda(G^{\prime})=|E^{\prime}|-|V^{\prime}|+1\geq\dim_{\mathbb{Z}}O(G^{\prime})\geq 3$. This observation implies that $|E^{\prime}|\geq|E|$ whenever $|V^{\prime}|>3$. Since $D_{3}$ is an directed clique we can identify any smaller graph $G^{\prime}$ with a subgraph of $D_{3}$. It is easy to see that deletion of any arc $e\in E$ produces an isolated cycle in $D_{3}$. For instance if we delete $e_{4}$ then $c_{3}$ will be isolated, if we delete $e_{5}$ then $c_{2}$ will be isolated and so on. Hence $G^{\prime}\cong D_{3}$. ∎

If $G$ is a non-resolvable graph then one can think about different methods to solve the FASP of the resolved graph $(S,\tau)$. One possibility is discussed in the next section.

7 The Bellman Decomposition

In this section we formulate an solution of the FASP based on a dynamic programming technique. Such an approach can be applied to optimization problems whenever there is a decomposition of the problem into subproblems which satisfy the Bellman principle, i.e., every optimal solution consists only of optimal subsolutions. To motivate the following definitions we first consider an example.

Example 4

Consider the graph $(G,\omega)$ in Figure 6. If we want to know, which arc of $c_{3}$ we have to cut for an optimal solution then this depends on the cycles $c_{1},c_{2}$. The benefit of cutting $e_{1}$ instead of $e_{2}$ or $e_{3}$ is that we do not have to cut $c_{1}$ anymore which costs at least $2$. Thus we introduce a new weight $\sigma$, which equals $\omega$ on $E\setminus\{e_{1}\}$ and is set to $\sigma(e_{1})=\omega(e_{1})-2=3$ on $e_{1}$. Since no other cycles than $c_{3}$ are cut by $d_{3}$ the weight of $d_{3}$ remains unchanged. Now we consider $H_{e_{3},d_{3}}=G_{\mathrm{el}}(e_{3})\setminus d_{3}$ and $H_{d_{3},e_{3}}=G_{\mathrm{el}}(d_{3})\setminus e_{3}$ and compute

[TABLE]

The best solution, which contains $e_{3}$ is $\{e_{3},d_{1}\}$ and the best solution containing $d_{3}$ is $\{d_{3},e_{1}\}$ and we observe that

[TABLE]

Thus, the difference of the solutions coincides with the difference of the subproblems in (6) with respect to the new weight $\sigma$.

We need to introduce several concepts to show that this observation remains true in general.

Definition 8 (arc sensitivity)

Let $(G,\omega)$ be a graph and $e,f\in E$, $e\not=f$ and $G_{\mathrm{el}}(e)=(V_{e},E_{e})$ be given. Then we say that $f$ is arc sensitive to $e$ with respect to the FASP, denoted by $f\twoheadrightarrow e$, if and only if

[TABLE]

We denote with \mathcal{N}_{\twoheadrightarrow}(e)=\left\{f\in E\,\,\big{|}\,\,f\twoheadrightarrow e\right\} the set of all arcs, which are sensitive to $e$.

Note that the arcs $f$ of an isolated cycle $c\in O_{I}(e)$ can not be sensitive to $e$. Thus, arc sensitivity detects arcs, which might prevent us from solving the FASP on $G_{\mathrm{el}}(e)$ independently from the rest of the graph. An understanding of these dependencies can be reached by understanding the meta graph of $G$ defined in the following.

Definition 9 (meta graph)

Let $(G,\omega)$ be a graph and $c\in\mathrm{O}_{\mathrm{el}}(G)$. We set $V_{0}=\mathcal{E}(c)$, $E_{0}=W_{0}=\emptyset$ and for $k\geq 1$ we define recursively $W_{k}=\cup_{i=0}^{k}V_{i}$ with

[TABLE]

Stopping the recursion if $K\in\mathbb{N}$ is such that $V_{K}=\emptyset$ we introduce the simple, undirected graph $M_{c}:=(V_{M_{c}},E_{M_{c}})=\bigcup_{k=0}^{K}(V_{k},E_{k})$ as the meta graph of $G$ with respect to $c$. Furthermore, we introduce $\mathcal{C}(c)=\mathcal{G}(\mathrm{O}_{\mathrm{el}}(V_{M_{c}}))\subseteq G$ as the subgraph of all arc sensitive cycles containing $c$.

Lemma 2

Let $G=(V,E)$ be a graph, $c\in\mathrm{O}_{\mathrm{el}}(G)$. Then we can construct the meta graph $M_{c}=(V_{M},E_{M})$ in $\mathcal{O}(|E|^{4})$.

Proof

Storing $G$ in an adjacency list enables us to test whether $\mathrm{O}_{\mathrm{el}}(e)\not=\emptyset$ in $\mathcal{O}(|E|)$ by depht first search. Due to Theorem 2 the graph $G_{\mathrm{el}}(e)$ can be determined in $\mathcal{O}(|E|^{2})$. Thus, testing sensitivity requires $\mathcal{O}(|E|^{2})$. Due to the fact that $V_{k}\cap V_{k-1}=\emptyset$, $1\leq k\leq K$ the construction of $M_{c}$ tests for sensitivity at most $|E|^{2}$ times, which yields the claimed complexity. ∎

Next we define the relative weight $\sigma_{G,c,e,f}$ with respect to some $c\in\mathrm{O}_{\mathrm{el}}(G)$ and $e,f\in\mathcal{E}(c)$. As it will turn out $\sigma_{G,c,e,f}$ decodes which arcs of $c$ can be cutted to obtain a minimal feedback set. For a better clarity we firstly restrict ourselves to the case where the meta graph $M_{c}\setminus f$ is a tree. In this case, $e$ is chosen as the root and $\sigma_{G,c,e,f}$ is given by solving the FASP for every leaf $h$ on

[TABLE]

and subtracting this value from the weight $\omega(h)$ of the predecessor of $h$. Afterwards, we delete all leafs of $M_{c}\setminus f$ and iterate this procedure till $e$ becomes a leaf. More precisely:

Definition 10 (relative weight for meta trees)

Let $(G,\omega)$ be a weighted graph $c\in\mathrm{O}_{\mathrm{el}}(G)$, $e,f\in\mathcal{E}(c)$. Let $M_{c}$ be the meta graph of $G$ with respect to $c$ and assume that the connected subgraph $M_{c,f}$ of $M_{c}\setminus f$, which contains $e$ is a tree. Then we define the relative weight of $G$ with respect to $c,e,f$

[TABLE]

as the output of Algorithm 3 with input \big{(}(G,\omega),M=M_{c},e,f\big{)}.

To define the relative weight in general, we have to consider all spanning trees of $M_{c,f}$ generated by deleting edges, which cut cycles for the first time, seen from $e$. In Example 5 we assert the definition for a special meta graph. The precise definition can be found below, using the following notions:

For any tree $M=(V_{M},E_{M})$ and any vertices $h,e\in V_{M}$ we denote with $p_{h}(M,e)$ the predecessor of $h$ with respect to root $e$. If $M=(V_{M},E_{M})$ is an arbitrary simple, undirected graph and $q\in V_{M}$ then we consider the set D_{k}(M,q):=\left\{w\in V_{M}\,\,\big{|}\,\,d(q,w)=k\right\} of all vertices possessing shortest path distance $k$ with respect to $q$ in $M$. Furthermore, we consider

[TABLE]

K(M,q):=\min\{k\in\mathbb{N}\,\,\big{|}\,\,U_{k}(M,q)\not=\emptyset\} and $U(M,q):=U_{K(M,q)}(M,q)$. In other words: $U(M,q)$ denotes the set of vertices $w$ which cut cycles for the fist time, seen from starting point $q$.

We set

[TABLE]

and denote with $M_{h}$ the connected component of $\widetilde{M}_{h}$ containing $h$. Recursively for $n\geq 1$ and an ordered set $F=\{h_{0},\dots,h_{n}\}$ with $h_{n}\in U(M_{h_{0},\dots,h_{n-1}},q)$ we define $M_{h_{0},\dots,h_{n}}:=(M_{h_{0},\dots,h_{n-1}})_{h{{}_{n}}}$. For $h\in U(M,q)$ we consider

[TABLE]

and set $M^{\leq U(M,q)}$ to be the connected component of $\widetilde{M}^{\leq U(M,q)}$ containing $q$, which is therefore a tree.

Definition 11 (relative weight in general)

Let $(G,\omega)$ be a graph, $c\in\mathrm{O}_{\mathrm{el}}(G),e,f\in\mathcal{E}(c)$ and $M_{c}$ be the meta graph of $G$ with respect to $c$. Let $\sigma_{G,M,e}$ be the output of Algorithm 4 with input \big{(}(G,\omega),M_{c},e,f\big{)}, $M=M_{c,e}$. Then, we define

[TABLE]

as the relative weight of $G$ with respect to $c,e,f$.

Example 5

Let $(G,\omega)$ be a graph, $c\in\mathrm{O}_{\mathrm{el}}(G),e_{0},e_{1}\in\mathcal{E}(c)$ and assume that $M_{c,e_{1}}$ coincides with $M$ from Example 7. We follow Algorithm 4 to compute $\sigma_{G,M,e_{0},e_{1}}:E\longrightarrow\mathbb{R}$. Observe that $U(M,e_{0})=\{f_{0},f_{1}\}$ and $U(M_{f_{1}},e_{0})=\{h_{0},h_{1}\}$. The graph $M_{f_{1},h_{0}}$ is sketched in the next picture and turns out to be a tree. Now we delete all vertices which are closer to $e_{0}$ as $p_{h_{1}}(M,e_{0})$ and obtain the graph $N$. Next we compute the relative weight $\sigma(h_{1}):=\sigma_{N,p_{h_{1}}(M,p_{h_{1}}),e_{1}}(h_{1})$ of $h_{1}$ with respect to $N,p_{h_{1}}$. Analogously, we compute $\sigma(h_{0}):=\sigma_{N,p_{h_{0}}(M,p_{h_{0}}),e_{1}}(h_{0})$ and consider the graph $M^{\leq}=M_{f_{1}}^{\leq U(M_{f_{1}},e_{0})}$, which is sketched in the last picture. Now $M^{\leq}$ is a tree and the predecessor of $f_{1}$ is $e_{0}$. Thus, we can compute

[TABLE]

which finishes the computation of $\sigma_{G,M,e_{0},e_{1}}:E\longrightarrow\mathbb{R}$ by replacing $\omega(f_{0}),\omega(f_{1})$ with these weights, respectively.

Proposition 5

Let $G=(V,E)$ be a graph $c\in\mathrm{O}_{\mathrm{el}}(G)$, $e,f\in\mathcal{E}(c)$ and let the meta graph $M_{c}$ of $G$ with respect to $c$ be given. Denote with

[TABLE]

the $\mathbb{Z}_{2}$-dimension of the cycle space $\Lambda^{0}(M_{c,f})$ of $M_{c,f}$. Then

1. i)

The computation of $\sigma_{G,c,e,f}$ can be realized in \mathcal{O}\big{(}2^{m(c,f)}|V||E|^{2}\log(|V|)\big{)}. 2. ii)

$m(c,f)\leq\dim_{\mathbb{Z}_{2}}\Lambda^{0}(G)=|E|-|V|+\#G$**

Proof

Assume that $M_{c,f}$ is a tree. As already mentioned, due to Dinic, (1970) and Dinitz, (2006) the feedback length $\Omega(G_{f},\sigma)$ can be determined in \mathcal{O}\big{(}|V||E|\log(|V|)\big{)}, and has to be computed at most $|V_{M_{c,f}}|\leq|E|$ times. If $M_{c,f}$ is not a tree then we observe that the most expansive computation step in Algorithm 4 is again the computation of the relative weight with respect to a certain subtree of $M_{c,f}$ (lines 8 and 19 in Algorithm 4). This computation step has to be computed for every pair $f,h\in U(M,q)$ twice, for some $M,q$. In worst case the combination of the pairs is independent, i.e., every other pair still appears once $M_{f}$ and $M_{h}$ are considered. In this case the set of remaining cycles running through the remaining pairs $(h^{\prime},f^{\prime})$ do not contain the edges $h$ and $f$ and can therefore not be generated by the cycles running through $(h,f)$ with respect to $\mathbb{Z}_{2}$-coefficients. Consequently, there are at most $2^{m(c,f)}$ iterations. Together with the argumentation above this yields $i)$.

To show $ii)$ we write $d=\{h_{0},\dots,h_{k}\}\subseteq V_{M_{c,f}}$ as a list of connected meta vertices. Now we choose connected paths $p_{i}\subseteq G$, $0=1\dots,k+1$ connecting $h_{i}^{-}$ with $h_{i+1\mod k+1}^{+}$. Then $c=\{\{h_{0}\}\cup p_{0}\cup\dots,\{h_{k}\}\cup p_{k}\in\Lambda^{0}(G)$ is a connected cycle in $G$. If $d_{1},\dots,d_{m}\in\Lambda^{0}(M_{c,f})$ is a set of $\mathbb{Z}_{2}$-linear independent meta cycles then regardless of choices for the paths $p_{i}$ representing an meta edge the corresponding cycles $c_{1},\dots,c_{m}\in\Lambda^{0}(G)$ are $\mathbb{Z}_{2}$-linear independent in $\Lambda^{0}(G)$. Thus, $m(c,f)$ is bounded by the $\mathbb{Z}_{2}$-dimension of $\Lambda^{0}(G)$, proving $ii)$. ∎

Remark 12

Note that a graph $G=(V,E)$ with $c\in\mathrm{O}_{\mathrm{el}}(G)$ such that $M_{c}$ coincides with $M$ in Figure 11 can be easily constructed by choosing a starting cycle $c$ and additional cycles $c_{1},c_{2}$ intersecting with $c$ in $f_{0},f_{1}$, respectively. Then we continue this process by follwing $M$ for the ramining cycles. Hence, the set of graphs $G$ with cycle $c$ and meta graphs $M_{c}$ such that $\dim_{\mathbb{Z}_{2}\Lambda^{0}}(M)$ is small, is actually a huge set.

Indeed the relative weights satisfy a Bellman condition, which can be formulated as follows. We recall that $\mathcal{G}_{o}(G)\subseteq G$ denotes the subgraph induced by all cycles of $G$ and state:

Theorem 4

Let $(G,\omega)$ be a weighted graph $c\in\mathrm{O}_{\mathrm{el}}(G)$, $e,f\in\mathcal{E}(c)$ and H_{e,f}=\mathcal{G}_{o}\big{(}G_{\mathrm{el}}(e)\setminus f\big{)},H_{f,e}=\mathcal{G}_{o}\big{(}G_{\mathrm{el}}(f)\setminus e\big{)} then

[TABLE]

where we shorten $\sigma_{e}=\sigma_{G,c,e,f},\sigma_{f}=\sigma_{G,c,f,e}$ and slightly abuse notation by still denoting $\sigma_{e},\sigma_{f},\omega$ for the restriction of the arc weights to the corresponding subgraphs, respectively.

Remark 13

Note that (4) is a quite comfortable way of formulating the Bellman principle, i.e., though we do not know the values of $\Omega(G\setminus e,\omega)$ and $\Omega(G\setminus f,\omega)$ we know that if $e$ maximizes

[TABLE]

for all on $\mathcal{E}(c)$. Then $\omega(e)+\Omega(G\setminus e,\omega)=\Omega(G,\omega)$. Thus, $\{e\}$ can be extended to a global optimal solution. Maybe this relative formulation can be applied also to other problems for which one wants to use a dynamic programming technique. The described observation is also used in the proof of Theorem 4.

In addition to the observation above the following statement is needed to prove Theorem 4.

Lemma 3

Let $G=(V,E)$ be a graph, $c\in\mathrm{O}_{\mathrm{el}}(G)$ and $e,f\in\mathcal{E}(c)$ such that there is $c^{\prime}\in\mathrm{O}_{\mathrm{el}}(e)\setminus\mathrm{O}_{\mathrm{el}}(f)$ and $p\in\mathcal{E}(c^{\prime})$. Then

1. i)

\sigma_{G,c,e,f}(h)=\sigma_{G\setminus f,c^{\prime},e,p}(h)\quad\text{for all}\,\,\,h\in\mathcal{E}\big{(}H^{\prime}_{e,p}\big{)}. 2. ii)

\sigma_{G,c,e,f}(p)=\omega(p)-\Omega\big{(}H^{\prime}_{p,e},\sigma_{G\setminus f,c^{\prime},p,e}\big{)}, where $H^{\prime}_{p,e}$ is understood with respect to $G^{\prime}$.

Proof

To verify $i)$ and $ii)$ one has to follow directly Definitions 9, 10 and 11, which is left to the reader. ∎

Proof (of Theorem 4)

If $\mathrm{O}_{\mathrm{el}}(e)=\mathrm{O}_{\mathrm{el}}(f)$ with respect to $G$ then $\sigma_{G,c,e,f}=\sigma_{G,c,f,e}$ and \Omega(G\setminus e,\omega)\big{)}=\Omega(G\setminus f,\omega)\big{)} and therefore the claim follows. Now we argue by induction on $|\mathrm{O}_{\mathrm{el}}(G)|$. If $|\mathrm{O}_{\mathrm{el}}(G)|=1$ then there is only one totally isolated cycle and therefore $\mathrm{O}_{\mathrm{el}}(e)=\mathrm{O}_{\mathrm{el}}(f)=\{c\}$. Thus, we are in a special case of the situation above and obtain the claim. Now assume that $|\mathrm{O}_{\mathrm{el}}(G)|>1$ and $\mathrm{O}_{\mathrm{el}}(f)\subsetneq\mathrm{O}_{\mathrm{el}}(e)$. We consider $G^{\prime}:=G\setminus f$ and observe that $|\mathrm{O}_{\mathrm{el}}(G^{\prime})|<|\mathrm{O}_{\mathrm{el}}(G)|$. We choose $c^{\prime}\in\mathrm{O}_{\mathrm{el}}(G^{\prime})$ with $e\in\mathcal{E}(c^{\prime})$ and choose $p\in\mathcal{E}(c^{\prime})$ such that

[TABLE]

is maximized on $c^{\prime}$, where $H^{\prime}_{e,p}$ ,$H^{\prime}_{p,e}$ are understood with respect to $G^{\prime}$. Thus, following Remark 4 there holds

[TABLE]

We set $\sigma^{\prime}_{e}:=\sigma_{G^{\prime},c^{\prime},e,p},\sigma^{\prime}_{p}:=\sigma_{G^{\prime},c^{\prime},p,e}$ then by induction and (10) we compute

[TABLE]

On the other side we consider $G^{\prime\prime}=(V^{\prime\prime},E^{\prime\prime}):=H_{e,f}$ with the arc weight

[TABLE]

Now observe that $H^{\prime\prime}_{p,e}=\emptyset$ with respect to $G^{\prime\prime}$ and by Lemma 3 $i)$ we have $\sigma^{\prime\prime}_{e}:=\sigma_{G^{\prime\prime},c^{\prime},e,p}(h)=\sigma_{G^{\prime},c,e,p}(h)=\gamma(h)$ for all $h\in\mathcal{E}(H^{\prime\prime}_{e,p})$, where $H^{\prime\prime}_{e,p}$ is understood with respect to $G^{\prime\prime}$. Moreover, $\gamma(e)=\omega(e)$ and therefore

[TABLE]

Due to Lemma 3 $ii)$ we have that $\sigma_{G,c,e,f}(p)=\omega(p)-\Omega(H^{\prime}_{p,e},\sigma^{\prime}_{p})$, $\sigma_{p}^{\prime}=\sigma_{G^{\prime},c^{\prime},p,e}$. Inserting this fact in (12) gives

[TABLE]

On the other, by (10) we have that (13) is maximized on $c^{\prime}$. Thus, again by induction

[TABLE]

Thus, by combining (13) with (11) and again (13) with (7) we obtain

[TABLE]

If $\mathrm{O}_{\mathrm{el}}(f)\subseteq\mathrm{O}_{\mathrm{el}}(e)$ then $\Omega(H_{f,e},\sigma_{G,c,f,e})=0$ and $\Omega(G\setminus\{e,f\},\omega)=\Omega(G\setminus e,\omega)$. Thus, by (15) this yields the claim. If $\mathrm{O}_{\mathrm{el}}(f)\subsetneq\mathrm{O}_{\mathrm{el}}(e)$ then the analogous of (15) with respect to $f$ yields

[TABLE]

Since $\Omega(G\setminus\{e,f\},\omega)-\Omega(G\setminus\{f,e\},\omega)=0$ this finishes the proof. ∎

We consider the Algorithms 5,6, denote with $\mathrm{output}(A)$ the set of all possible outputs an algorithm $A$ can produce and conclude :

Corollary 1

Let $(G,\omega)$ be a graph then the algorithm CUT is exact and complete with respect to the FASP, i.e,

[TABLE]

while the algorithm CUT & RESOLVE is exact, i.e.,

[TABLE]

*Moreover, there is $m\in\mathbb{N}$ such that CUT * and CUT & RESOLVE possess run times \mathcal{O}\big{(}2^{m}|E|^{4}\log(|V|)\big{)}, where the parameter $m\leq|E|-|V|+1$ can be determined in \mathcal{O}\big{(}|E|^{3})\big{)}.

Indeed the corollary proves Theorem B and by Remark 12 the set of graphs with small $m\in\ N$, $m<<|V|$ is a huge set.

Proof

If $\varepsilon\in\mathcal{S}(G,\omega)$ and $e\in\varepsilon$ then $\varepsilon\setminus\{e\}$ solves the minimal FASP on $G(E\setminus\{e\})$. Thus, the exactness and completeness statements follow directly from Theorems 4,3. For $c\in\mathrm{O}_{\mathrm{el}}(G)$ we set $m(c):=\dim_{\mathbb{Z}_{2}}\Lambda^{0}(M_{c})$. Then $m(c)\geq m(c,f)$ for all $f\in\mathcal{E}(c)$ with $m(c,e)$ and Proposition 5 implies that $m(c)\leq|E|-|V|+1$ holds. Let $c^{\prime}\in\mathcal{C}(c)$ with $c^{\prime}\cap c\not=0$, be any cycle in the component of all arc connected cycles containing $c$, introduced in Definition 9. Then $m(c^{\prime})\leq m(c)$ on $G\setminus e$ for every $e\in\mathcal{E}(c)$. Thus, as long as at least one arc $e\in\mathcal{E}(c)$ was deleted the maximal number of appearing $\mathbb{Z}_{2}$-linear independent meta cycles appearing for the computation of $\sigma_{G\setminus e,c^{\prime},h,k}$, $h,k\in\mathcal{E}(c^{\prime})$ are bounded by $m(c)$. Thus, by setting $G_{0}=G$, choosing a cycle $c_{0}\in\mathrm{O}_{\mathrm{el}}(G_{0})$, determing $M_{c_{0}}=(V_{M_{c_{0}}},E_{M_{c_{0}}})$ and considering $c_{k}\in\mathrm{O}_{\mathrm{el}}(G_{k})$, $G_{k}=G_{k-1}\setminus V_{M_{c_{k-1}}}$, $k\geq 1$ we obtain cycles $c_{0},\dots,c_{n}$, $n\leq|E|$ with $\mathcal{C}(c_{i})\cap\mathcal{C}(c_{j})=\emptyset$ and $\cup_{i=0}^{n}\mathcal{C}(c_{i})=\mathrm{O}_{\mathrm{el}}(G)$. Since \mathcal{E}\big{(}\mathcal{C}(c_{i})\big{)}\cap\mathcal{E}\big{(}\mathcal{C}(c_{j}))=\emptyset the parameter $m:=\max_{i=0,\dots,n}m(c_{i})$ can be determined in \mathcal{O}\big{(}|E|^{3}\big{)} due to Lemma 2.

Since the algorithm CUT computes $\sigma_{G,c,e,h},\sigma_{G,c,h,e}$ for fixed $e\in\mathcal{E}(c)$ and all $h\in\mathcal{E}(c)$ and $|\mathcal{E}(c)|\leq|V|$, cuts the right arc and repeats the computation at most $|E|$ times by observing that $\mathcal{O}(|V|^{2})=\mathcal{O}(|E|)$ the run time of the algorithm CUT can be estimated as claimed. Recall, that due to Theorem 3 the resolved graph can be computed in \mathcal{O}\big{(}|V||E|^{3}\big{)}. Therefore, the analogous argumentation yields the claimed run time for the algorithm CUT & RESOLVE. ∎

Due to the fact that the FASP is NP complete, as expected our approach depends exponentially on some parameter, which in our approach is the number $m$ of liner independent meta cycles. In cases where $m$ is large we have to use another method to solve the FASP or use a heuristic.

8 Valid Greedy Solutions

As for instance shown in Ispolatov and Maslov, (2008) a greedy solution for the FASP needs not to be optimal. We give a criterium on solutions which guarantees optimility. Moreover, we can estimate the failure of every sub optimal solution. Finally, we suggest a heuristic given by a hybrid technique of the already presented approaches.

For given graph $(G,\omega)$ we consider the functions

[TABLE]

Recall, that due to Arora and Barak, (2009) determing $\theta_{G}(e)$ is a NP-hard problem and the results of Tarjan, (1973) and Johnson, (1975) solving the problem in \mathcal{O}\big{(}|\theta_{G}(e)|(|V|+|E|)\big{)}, where $\theta_{G}(e)$ can depend exponentially on $G$. However, Roberts and Kroese, (2007) could establish efficent and close estimations of the number of $s-t$ paths. Since $|\mathrm{O}_{\mathrm{el}}(e)|=|P_{\mathrm{el}}(e^{-},e^{+})|$ the result enables us to determine $\theta_{G}(e)$ efficently, with small failure. In contrast, $\varphi_{G}$ can be determined in $\mathcal{O}(|E|^{2})$ due to Theorem 2.

Definition 12

Let $(G,\omega)$ be a given weighted graph with $G=G_{o}$. We introduce the efficent weights

[TABLE]

and set $\omega_{\max}(G,\omega):=\max_{e\in E}\omega(e)$, $\theta_{\max}(G)=\max_{e\in E}\theta_{G}(e)$, $\varphi_{\max}(G):=\max_{e\in E}\varphi_{G}(e)$, $\xi_{\max}(G,\omega):=\max_{e\in E}\xi(e)$, $\eta_{\max}(G,\omega):=\max_{e\in E}\eta(e)$, and $\mu(G,\omega)=\left\lceil\frac{|\mathrm{O}_{\mathrm{el}}(G)|}{\xi_{\max}(G,\omega)}\right\rceil$, $\upsilon(G,\omega)=\left\lceil\frac{|E|}{\eta_{\max}(G,\omega)}\right\rceil$, where $\left\lceil\cdot\right\rceil$ denotes the Gauss-bracket.

Theorem 5

Let $(G,\omega)$ be a given graph. Then

[TABLE]

Moreover, there are infinitely many weighted graphs $(G,\omega)$ with $\mu(G,\omega)=\Omega(G,\omega)$ or $\upsilon(G,\omega)=\Omega(G,\omega)$.

Proof

Let $\varepsilon=\{e_{1},\dots,e_{n}\}\in\mathcal{S}(G,\omega)$ be an arbitrarily ordered solution of the weighted FASP. We set $G_{0}=G$ and G_{i}=(V_{i},E_{i}):=G\big{(}E\setminus\{e_{1},\dots,e_{i-1}\}\big{)}, for $i\geq 1$ and denote with $\omega_{G_{i}}$ the corresponding restriction of $\omega$ to $E_{i}$. Now due to the fact that $\xi_{G_{i},\omega_{i}}(e)\leq\xi_{G,\omega}(e)$ for all $e\in E_{i},1\leq i\leq n$ we obtain

[TABLE]

Since $\Omega(G,\omega)\in\mathbb{N}$ this proves (16). Now let $(G,\omega)$ be a graph with $\omega=1$ and $\mathrm{O}_{\mathrm{el}}(G)=\{c_{0},\dots,c_{n}\}$, which are arranged path like, i.e., $|\mathcal{E}(c_{i}\cap c_{i+1})|=1$ and $|\mathcal{E}(c_{i}\cap c_{j})|=0$ if $j\not=i$. Then one verifies easily that $\mu(G,\omega)=\Omega(G,\omega)$. By replacing $\theta_{G}$ with $\varphi_{G}$ and $\mu(G,\omega)$ with $\upsilon(G,\omega)$ the exact same argumentaion yields the remaining claim. ∎

Remark 14

Note, that of course there are many more graphs with $\mu(G,\omega)=\Omega(G,\omega)$ or $\upsilon(G,\omega)=\Omega(G,\omega)$ then those used in the proof above. Nevertheless, it is hard to give a good condition on a graph such that $\mu(G,\omega)=\Omega(G,\omega)$ or $\upsilon(G,\omega)=\Omega(G,\omega)$ holds. For instance in Ispolatov and Maslov, (2008) an example of a planar graph is given, where this is not the case. Certainly, the lower bounds can be used to improve the performance of a variety of algorithms solving the FASP or to control the quality of a heuristic as the Algorithms 7,8.

We consider the heuristics GREEDY-CUT and GREEDY-CUT & RESOLVE and discuss their properties.

Proposition 6

Let $(G,\omega)$ and $\varepsilon$ be a solution of GREEDY-CUT or GREEDY-CUT & RESOLVE with respect to the effective weigth $\xi$. Then

[TABLE]

If in particular $\omega\equiv 1$ then $|\varepsilon|\leq|E|/2$.

Proof

Certainly, it suffices to prove the first estimate in (17). We show the claim for a solution $\varepsilon=\{e_{1},\dots,e_{n}\}$ of GREEDY-CUT. Assume that $\varepsilon$ is ordered with respect to appearing arcs, set $G_{1}=G$ and G_{i}=(V_{i},E_{i}):=G\big{(}E\setminus\{e_{1},\dots,e_{i-1}\}\big{)}, for $i\geq 2$ and denote with $\omega_{G_{i}}$, $\theta_{G_{i}}$ and $\xi_{G_{i}}$ the corresponding restrictions of $\omega$, $\theta_{G}$, $\xi_{G}$ to $E_{i}$. Then we compute

[TABLE]

Since $\xi_{n}(e_{n})=\xi_{\min}(\varepsilon)$ we use Theorem 5 to compute

[TABLE]

and the claim follows. A proof of the statement for GREEDY-CUT & RESOLVE can be given by an easy adaption of the argument above and is left to the reader. Now let $\omega\equiv 1$ then we argue by induction on $|\varepsilon|$ to show that $|\varepsilon|\leq|E|/2$ for both algorithms. If $|\varepsilon|=1$ then due to the fact that $G$ possesses no loops the claim follows. Now let $|\varepsilon|>1$ we order $\varepsilon=\{e_{1},\dots,e_{n}\}$ with respect to appearance and consider $\varepsilon_{1}:=\{e_{1}\}$, $\varepsilon_{2}=\varepsilon\setminus\{e_{1}\}$ and $G_{1}:=(V_{1},E_{1})=G(\mathrm{O}_{\mathrm{el}}(e_{1}))$, $G_{2}:=(V_{2},E_{2})=G(\mathrm{O}_{\mathrm{el}}(\varepsilon_{2}))$. By induction we have $|\varepsilon_{2}|\leq|E_{2}|/2$. Consider $G/(E_{1}\cap E_{2})$, delete all appearing loops and denote the resulting graph with $G_{1}^{*}$. If $\mathrm{O}_{\mathrm{el}}(G_{1}^{*})=\{c\}$ then all cycles of $\mathrm{O}_{\mathrm{el}}(G)$ are totally isolated and the claim follows by triviality. If $|\mathrm{O}_{\mathrm{el}}(G_{1}^{*})|>1$ then $|\varepsilon_{1}|\leq|E_{1}^{*}|/2$. Since $E_{1}^{*}\cap E_{2}=\emptyset$ this implies that

[TABLE]

as claimed. ∎

Be replacing $\theta_{G}$, with $\varphi_{G}$ and $\xi_{G,\omega}$ with $\eta_{G,\omega}$ in GREEDY-CUT or GREEDY-CUT & RESOLVE the analouge argumentaion yields.

Proposition 7

Let $(G,\omega)$with $G=G_{o}$ and $\varepsilon$ be a solution of GREEDY-CUT or GREEDY-CUT & RESOLVE with respect to the effective weigth $\eta$. Then

[TABLE]

If in particular $\omega\equiv 1$ then $|\varepsilon|\leq|E|/2$.

Example 6

Consider the directed clique $D_{3}$ from Figure 5 with constant weight $w\equiv 1$. Then $D_{3}$ coincides with its resolved graph and regardless of possible choices every candidate $\varepsilon$ the algorithm GREEDY-CUT or GREEDY-CUT & RESOLVE proposes, satisfies $|\varepsilon|=3$. Since $\mu(D_{3})=3$ every candidate is optimal.

Summarizing our results so far the heuristics GREEDY-CUT or GREEDY-CUT & RESOLVE solve the FASP with controlled variance in $\mathcal{O}(|E|^{4})$, due to Theorem 2, in case of effective weight $\eta$ and in $\mathcal{O}(f_{\theta}|E|^{2})$ in case of effective weight $\xi$, where $f_{\theta}$ shall control the computation steps of $\theta_{G}(e)$, $\forall e\in E$. Even if we approximate $\theta(e)$ by the method of Roberts and Kroese, (2007) the resulting algorithm remains an efficient heuristic. However, possibly there is a more accurate method available, given by a hybrid algorithm of the methods introduced in this article. We expect that an implementation of this strategy yields a fast and precise general FASP-SOLVER, which due to section 2 is therfore also a FVSP-SOLVER.

Strategy 1

For given weighted graph $(G,\omega)$

Compute the resolved graph $(S,\tau)$. 2. 2.

Choose a cycle $c\in\mathrm{O}_{\mathrm{el}}(S)$ and compute the meta graph $M_{c}$. 3. 3a.

If the number of meta cycles $m(c)$ is large determine a “good” feedback vertex set $\nu_{M_{c}}$ of $M_{c}$, with respect to the vertex weight

[TABLE]

using one of the known or presented methods. 4. 3b.

Alternatively, compute a maximal spanning tree $T_{M_{c}}=(V_{m_{c}},E_{T_{c}})$ with respect to the arc weight weight

[TABLE]

where $h=[e,f]$, $e,f,\in E$ and set $\varepsilon_{M_{c}}=E_{M_{c}}\setminus E_{T_{c}}$. 5. 4.

Set $G^{*}=G/\nu_{M_{c}}$ and use CUT or CUT & RESOLVE to solve the FASP on on the component $\mathcal{C}(c^{*})$ of arc connected cycles containing $c^{*}=c/\nu_{M_{c}}$. 6. 5.

Choose a new cycle $c^{\prime}$ of the resulting graph and repeat 1.-4. until no such cycle exists. 7. 6.

Use the backtracking procedure of Proposition 3 to compute a feedback arc set $\varepsilon\subseteq E$ of $G$.

The union $\nu_{M}$ of the meta feedback vertex sets of the meta graphs can be interpreted as arcs, which are forbidden to cut in $G$. The resulting feedback arc set $\varepsilon$ will be optimal up to this obstruction ,i.e., we have

[TABLE]

where $M_{c^{*}}^{*}$ denotes the meta graph of $G^{*}$ with respect to $c^{*}=c/\nu_{M}$. Hence, the quality of this heuristic can be evaluated by measuring how good $G^{*}$ approximates $G$. Thus, if $|\nu_{M}|<<|E|$ and the weight of the forbidden arcs is very high, i.e.,

[TABLE]

the arcs of $\nu_{M}$ will probably not be contained in any optimal solution, yielding the correctness of Strategy 1. Additionally, the lower bounds $\mu(G,\omega)$, $\upsilon(G,\omega)$ from section 8 can be used to validate correctness. Analogous controls can be thought of, if we choose the alternative 3b.

9 Discussion

An implementation of the described algorithms is planned to be realized. Certainly, a comparison of real run times with other approaches would be of great interest. So far we compare our results with other theoretical approaches. Due to the immense amount of results during the last decades we restrict our discussion to publications, which do not restrict themselves to very tight graph classes as tournaments Karpinski and Schudy, (2010) or reducible flow graphs Ramachandran, (1988). Finally, we suggest how the approaches of this article might be adapted to related problems.

9.1 The algorithms from I. Razgon and J. Chen et al.

Note that for given graph $(G,\omega,\gamma)$ with arc weight $\omega$ and vertex weight $\gamma$ a brute force method of solving the FASP/FVSP is given by considering every subset $\varepsilon\subseteq E$ or $\nu\subseteq V$ and check whether the graphs $G\setminus\varepsilon$, $G\setminus\nu$ are acyclic, respectively. Since due to Remark 5, checking for acyclicity requires $\mathcal{O}(|E|^{2})$ operations, we can generate a list of all FAS’s or FVS’s possessing length $l_{A}\leq|\mathcal{P}(E)|=2^{|E|}$, $l_{V}\leq|\mathcal{P}(V)|=2^{|V|}$. Choosing the cheapest FAS or FVS yields therefore a brute force algorithm solving the FASP/FVSP in $\mathcal{O}(|2^{|E|}|E|^{2}|)$, $\mathcal{O}(2^{|V|}|E|^{2})$, respectively.

The algorithm of Razgon, (2007) solves the unweighted FVSP on simple graphs in $\mathcal{O}(1.9977^{|V|}|V|^{\mathcal{O}(1)})$. Compared to the brute force algorithm this yields almost no improvement. Therefore, the question occured whether the parametrised version of the FVSP could be solved by an fixed parameter tractable algorithm. Every NP-complete problem can be solved by a fixed parameter tractable algorithm, i.e., by choosing $p$ as the problem size there is an algorithm with complexity $\mathcal{O}(f(p))$, where $f$ is an on the parameter $p$ exponentially depending function. Thus, the term fixed parameter tractable could be misleading. The precise question is whether there exists an algorithm with run time $\mathcal{O}(f(k)|V|^{\mathcal{O}(1)})$ computing a FVS of length less than $k$ or determing that no such set exists. Since the FVSP is NP-complete the function $f$ will be exponentially dependent on $k$ unless $P=NP$. Indeed, the algorithm of Chen et al., (2008) solves the parametrised version of the FVSP in \mathcal{O}\big{(}|E|^{4}4^{k}k^{3}k!\big{)}. Thus, $f(k)=k^{3}4^{k}k!$, increases even worse than exponentially in $k$. Since a small feedback length almost always correlates to small graphs or very special graphs, e.g. tree-like graphs, even improvements of the algorithm won’t be usefull in many applications. Therefore, the article might be seen as an purely theoretical approach answering this question. Indeed, to the best of our knowledge none of the algorithms were used for an implementation of a general FVSP/FASP-SOLVER.

In contrast, the algorithms CUT and CUT & RESOLVE solve the FASP or FVSP on weighted multi-digraphs in \mathcal{O}\big{(}2^{m}|E|^{4}\log(|V|)\big{)} and \mathcal{O}\big{(}2^{n}\Delta(G)^{4}|V|^{4}\log(|E|)\big{)}. The parameters $m$ and $n$ fulfill $m\leq|E|-|V|+1$, $n\leq(\Delta(G)-1)|V|-|E|+1$ and can be computed in $\mathcal{O}(|E|^{3})$, $\mathcal{O}(\Delta(G)^{3}|V|^{3})$, respectively. Thus, in both cases we can efficently control the run time of the exact solutions, which enables us to a priori decide whether the given instance shall be solved exactly or by an heuristic, e.g., Strategy 1. This crucial difference to the other approaches and the fact that Strategy 1 is an heuristic on the meta level and not on the instance itsself, makes us confident that an implementation generates a fast and accurate FASP/FVSP-SOLVER yielding a deep impact on computational an applied sciences.

9.2 The polytop approach from C. Lucchesi et al.

In Lucchesi and Younger, (1978) and Grötschel et al., (1985) a polytope of arc sets is assigned to a given graph. The FASP translates to solve a certain linear optimization over this polytope. In the case of planar or more general * weakly acyclic graphs* the polytope is integral, i.e., it possesses integral corners. Since the optimum will be obtained in at least one of the corners, one can apply the so called ellipsoid method for submodular functions Grötschel et al., (1981) to find the right corner in polynomial time, see also Martí and Reinelt, (2011) for further details. The approach is certainly remarkable though it contains some weaknesses.

The first problem is that though the algorithm runs in polynomial time the degree of the polynomial depends on a variety of parameters and cannot be estimated by hand a priori. Therefore, there are planar or weakly acyclic graphs, which can be solved efficently from a theoretical view point but actually a computer based implementation of the approach cannot ensure to meet a performance behavior applications require.

Secondly, the class of weakly acyclic graphs is not well classified yet. Thus, if we leave the class of planar graphs it is hard to say whether a given graph is weakly acyclic or not. For instance, consider the directed version of $K_{3,3}$ in Figure 8, which is known to be a weakly acyclic graph. Then it is not hard to see that the meta graph $M_{c}$ contains only one cycle. Thus, though $G$ is not resolvable we can efficently solve the FASP on $G$ by applying the algorithm CUT. In fact all examples of weakly acyclic graphs given in Grötschel et al., (1981) turn out to be efficently solvable by CUT. Since the techniques of this article are not sensitive to topological obstructions as planarity we expect that there are instances of the FASP, which are neither planar nor weakly acyclic and even though can be solved efficently by CUT or CUT & RESOLVE. On the other hand, by arranging cycles along several meta cycles it is quite easy to construct a planar graph $G$ with a number $m\in\mathcal{O}(|E|)$ of linear independent meta cycles. Thus, so far none of the approaches can state to solve the “larger” instance class efficently.

However, a deeper understanding of the meta graphs and their topology seem to be the most relevant tasks for further research, which might enable us to classify weakly acyclic graphs and yield a completely new perspective to other questions in graph theory.

9.3 Heuristics

In Saab, (2001) a summarization of heuristic approaches is given and several new ones are introduced. The weak point in all these approaches is that they do not provide an non-empirical control of the variance of the heuristical solution from the optimum. Therefore, it is impossible to guarantee whether a solution is tight to the optimum. In Huang et al., (2013) a good lower bound of the feedback length for Eulerian graphs is given and therefore it would be interesting how our bound behaves on this graph class. In general, if $\varepsilon$ denotes a feedback set the heuristic GREEDY-CUT or GREEDY-CUT & RESOLVE proposes, then by Proposition 6 we have shown that if $\omega\equiv 1$ then

[TABLE]

yielding a controlled variance, as long as $\mu(G)$ can be determined or estimated from below, see Roberts and Kroese, (2007). We conjecture that (18) improves the known estimates given by Berger and Shor, (1990). Furthermore, other heuristics can be improved by the results of this article. For instance, the counter example for the Greedy approach introduced in Ispolatov and Maslov, (2008) is resolvable and therefore RESOLVE & CUT closes this gap. A comparison of Strategy 1, with the approximation of Even et al., (1998), where an approximation ratio in \mathcal{O}\big{(}\log(\Omega)\log(\log(\Omega))\big{)} was established, might be usefull as well.

9.4 New Approach to the Subgraph Homeomorphism Problem

The NP-complete directed subgraph homeomorphism problem studied by Fortune et al., (1980) is to consider two given graphs $G=(V_{G},E_{G})$ and $P=(V_{P},E_{P})$ together with an injective mapping $m:V_{P}\longrightarrow V_{G}$ of vertices of $P$ into the vertices of $G$. Now the problem is given by deciding whether there exists a injective mapping from arcs of $P$ into pairwise node disjoint elementary paths of $G$ such that an arc $f$ with head $h$ and tail $t$ is mapped on an elementary path from $m(t)$ to $m(h)$. For a given graph $G^{*}=(V^{*},E^{*})$, $u,v\in V^{*}$, $f=(p,q)\in E^{*}$ with $\deg(p),\deg(q)\geq 2$ we consider a special instance of the directed subgraph homeomorphism problem by setting $G=G^{*}\setminus f$, $P=(V_{P},E_{P})$ with $V_{P}=\{a,b,c,d\}$, $E_{P}=\{(a,b),(c,d)\}$ and $m(a)=u,m(b)=p,m(c)=q,m(d)=v$. Thus, solving the subgraph homeomorphism problem with respect to these special instances is equivalent to decide whether $f$ is an arc of $G^{*}_{\mathrm{el}}(u,v)$. Hence, in addition to the polynomial time solvable instance classes known from Fortune et al., (1980), e.g. stars, also problem instances as defined above are polynomial time solvable due to Theorem 2. Potentially, our observations can be generalized in regard of this problem.

Appendix A The Essential Minor

This section is used to prove Propositions 2 and 3. To do so we recall that \mathrm{O}_{\mathrm{el}}(e)=\left\{c\in\mathrm{O}_{\mathrm{el}}(G)\,\,\big{|}\,\,e\in\mathcal{E}(c)\right\} and state the following Lemmas.

Lemma A.1

Let $G$ be a graph, $\omega:E\longrightarrow\mathbb{N}^{+}$ be an arc weight, $\varepsilon\in\mathcal{S}(G,\omega)$, and $e\in E$. Then either

[TABLE]

*If in particular, $\varepsilon\cap[e]_{\sim_{\Gamma}}\not=\emptyset$ then $\varepsilon\cap[e]_{\sim_{\Gamma}}$ minimizes $\omega$ on $[e]_{\sim_{\Gamma}}$. *

Proof

Let $\varepsilon\in\mathcal{S}(G,\omega)$ and $e\in\varepsilon$. Since every arc $f\in E$ with with $e\sim_{\Gamma}f$ is connected by a branch point free path with $e$ we have that $\mathrm{O}_{\mathrm{el}}(e)=\mathrm{O}_{\mathrm{el}}(f)$. Thus, at most one arc in $[e]_{\sim_{\Gamma}}$ will be cutted and this arc has to minimize $\omega$ on $[e]_{\sim_{\Gamma}}$. ∎

Lemma A.2

Let $G$ be a positively weighted graph, $\varepsilon\in\mathcal{S}(G,\omega)$, and $e\in E$. Then either

[TABLE]

Proof

Let $\varepsilon\in\mathcal{S}(G,\omega)$ and $e\in\varepsilon$. Assume there is $f\in F^{+}(e)\setminus\varepsilon$ then certainly $\varepsilon\cap F^{-}(e)=F^{-}(e)$ otherwise there would be a two-cycle that is not cutted. Now, let $e_{1},\dots,e_{k}\in\varepsilon\setminus F(e)$, $k\in\mathbb{N}$, be such that $\mathrm{O}_{\mathrm{el}}(f)\cap\mathrm{O}_{\mathrm{el}}(e_{i})\not=\emptyset$ and \mathrm{O}_{\mathrm{el}}\big{(}\{e_{1},\dots,e_{k}\}\big{)}\supseteq\mathrm{O}_{\mathrm{el}}(f). Since the cycles in $\mathrm{O}_{\mathrm{el}}(e)$ and $\mathrm{O}_{\mathrm{el}}(f)$ differ only in a single arc, i.e., $e$ and $f$, it suffices to cut the arcs $F^{-}(e)\cup\{e_{1},\dots,e_{k}\}$ to cut all cycles in $\mathrm{O}_{\mathrm{el}}(e)$, i.e.,

[TABLE]

Since $F^{-}(e)\cup\{e_{1},\dots,e_{k}\}$ is therefore a cheaper possibility than $\varepsilon$ cutting $\mathrm{O}_{\mathrm{el}}(F^{+}(e))$, this contradicts that $\varepsilon\in\mathcal{S}(G,\omega)$ and yields the claim. ∎

Now we state again Proposition 2 and deliver its proof.

Proposition A.1

Let $G=(V,E,\omega)$ be a positively weighted graph with essential minor $(C,\omega_{C})$ and let $\varepsilon\in\mathcal{P}(E)$ and $\varepsilon_{C}$ be the image of $\varepsilon$ in $(C,\delta)$. Then

[TABLE]

In particular $\Omega(G,\omega)=\Omega(C,\omega_{C})$.

Proof

Let $\varepsilon\in\mathcal{P}(E)$ and $\varepsilon_{1}=(\varepsilon/_{\Gamma})/_{\Phi}\subseteq E_{1}$ be the image of $\varepsilon$ in $G_{1}=(G/_{\Gamma})/_{\Phi}$. We recall that $\omega_{1}=(\omega/_{\Gamma})/_{\Phi}$ was defined in Definition 4 and show that

[TABLE]

Assume that $\varepsilon\in\mathcal{S}(G,\omega)$ then by Lemmas A.1, A.2 and the construction of $(G_{1},\omega_{1})$ we obtain $\Omega_{G,\omega}(\varepsilon)=\Omega_{G_{1},\omega_{1}}(\varepsilon_{1})$. Thus, if $\varepsilon_{1}\not\in\mathcal{S}(G_{1},\omega_{1})$ then we choose $\alpha\in\mathcal{S}(G_{1},\omega_{1})$ and a FAS $\varepsilon^{\prime}\subseteq E$ of $G$ such that the equations 19,20 hold and $(\varepsilon^{\prime}/_{\Gamma})/_{\Phi}=\alpha$. Consequently $\Omega_{G,\omega}(\varepsilon^{\prime})=\Omega_{G_{1},\omega_{1}}(\alpha)$ and therefore due to the construction of $(G_{1},\omega_{1})$ we get

[TABLE]

which contradicts that $\varepsilon\in\mathcal{S}(G,\omega)$. Thus, $\varepsilon_{1}\in\mathcal{S}(G_{1},\omega_{1})$.

Vice versa assume that $\varepsilon\subseteq E$ is such that $\varepsilon_{1}\in\mathcal{S}(G_{1},\omega_{1})$. We claim that equations 19,20 are satisfied by $\varepsilon$. Assume the opposite then due to Lemmas A.2 and A.1 we can delete an arc $e\in\varepsilon$ or replace an arc $e\in\varepsilon$ by an arc $f\in[e]_{\sim_{\Gamma}}$ with $\omega(e)>\omega(f)$. If this is not the case then we can delete all arcs $f\in F^{+}(e)\cap\varepsilon$ whenever $e$ is such that $\emptyset\not=F^{+}(e)\cap\varepsilon\not=F^{+}(e)$. If $\varepsilon^{\prime}$ denotes this modified set, then $\varepsilon^{\prime}$ is FAS of $G$ and in all cases

[TABLE]

A contradiction! Hence, the equations 19,20 hold for $\varepsilon$ and therefore the construction of $(G_{1},\omega_{1})$ yields

[TABLE]

Thus, if $\varepsilon\not\in\mathcal{S}(G,\omega)$ then we choose $\beta\in\mathcal{S}(G,\omega)$ and obtain that $\beta_{1}$ is a FAS of $G_{1}$ with $\Omega_{G,\omega}(\beta)=\Omega_{G_{1},\omega_{1}}(\beta)<\Omega(G_{1},\omega_{1})$, which is impossible. Hence $\varepsilon\in\mathcal{S}(G,\omega)$ and the claim follows by iteration of these arguments. ∎

Proposition A.2

Let $G=(V,E,\omega)$ be a finite, connected, directed, weighted multigraph then we can construct $(C,\delta)$ in time $\mathcal{O}(|V||E|^{2})$. Furthermore, there is an algorithm with run time $\mathcal{O}(|E|^{2})$ which constructs a solution $\varepsilon\in S(G,\omega)$ given a solution $\varepsilon_{C}\in\mathcal{S}(C,\delta)$.

Proof

The graph $G/_{\Gamma}$ can be computed in a single iteration over $V$, see Algorithm 9, where $\kappa$ is explained later. Each non branching node $v$ is removed and its two incident arcs $e=(u,v)$ and $f=(v,w)$ are replaced by an arc $(u,w)$ with weight $\min\{\omega(e),\omega(f)\}$. This is possible in time $\mathcal{O}(|V|)$ if the graph is represented as adjacency list where the targets of the outgoing arcs and the origins of the ingoing arcs are stored separately. Thus, the iteration over $V$ yields a runtime of $\mathcal{O}(|V|^{2})$.

To construct $G/_{\sim}\Phi$ we iterate over $E_{\Gamma}$ yielding $G_{1}=(G/_{\Gamma})/_{\Phi}$, see Algorithm 10. For each arc $e$ the weight is updated to $\sum_{e^{\prime}\in F^{+}(e)}\omega(e^{\prime})$ and the parallel arcs $F^{+}(e)\setminus\{e\}$ are purged from the graph. This can be realized in time $\mathcal{O}(|E|+|V|)=\mathcal{O}(|E|)$ with a counting sort prepossessing step if the target nodes in the adjacency list are stored such that equal targets are stored consecutively. Because the construction of $(C,\delta)$ requires at most $|V|$ iteration steps, i.e., if $(G_{K},\omega_{K})=(C,\delta)$ then $K\in\mathcal{O}(|V|)$, the essential minor $(C,\delta)$ can be computed in time $\mathcal{O}(|V|^{3}+|V||E|^{2})\subseteq\mathcal{O}(|V||E|^{2})$.

A simple extension of the algorithms allows to compute the information that is necessary to compute a solution $\varepsilon\in\mathcal{S}(G,\omega)$ once $\varepsilon_{C}\in\mathcal{S}(C,\delta)$ is given. During the application of $\Gamma$ and $\Phi$ we store the set of arcs of $G$ that are part of a solution if the corresponding arc from $G/_{\Gamma}$ and $G/_{\Phi}$, respectively, are in a FAS. That is, an arc that gave the minimum weight of the two arcs in a non branching path or all parallel arcs, respectively, see Lemma A.1 and Lemma A.2 . In Algorithms 9,10 this is realized by $\kappa$ which can be considered as $\kappa:E_{C}\longrightarrow\mathcal{P}(E)$. The mapping is initialized as $\kappa(e)\leftarrow\{e\}$. Storing the arcs as linked list allows to update $\kappa$ in linear time, i.e., the asymptotic run time of Algorithms 10,9 remains unchanged. Note that, $\arg\!\min^{*}$ returns only one arc in the case of equality. Now, replacing each arc $e\in\varepsilon_{C}$ by $\kappa(e)$ yields $\varepsilon$, which due to Proposition 2 is a solution for the FASP on $(G,\omega)$. Thus, the replacement can be realized in time $\mathcal{O}(|E|^{2})$. ∎

Remark 15

Algorithm 9 may be extended to generate all solutions of the FASP for $G$ given all solutions for the FASP on $C$. Therefore the equal weight alternatives in a non branching path need to be stored. The generation of the combinations of the alternatives of different paths yields all solutions. Certainly, then the run time depends exponentially on the number of possible combinations of alternatives.

Acknowledgements.

I want to thank Matthias Bernt for many fruitful discussions and his support for formalizing some algorithms. Moreover, a heartful thank you goes to Peter F. Stadler for the nice time at his bioinformatics institute in Leipzig.