A $\frac{3}{2}$-Approximation Algorithm for Tree Augmentation via Chv\'atal-Gomory Cuts
Rounding the Solution
Given a instance $`(G,L,c)`$, the overall algorithm begins by solving the for $`\gamma = \ceil{\tfrac{28M}{\eps^2}}`$. It then computes a solution for $`L^h`$ for the heavily-covered edges $`E^h`$, contracts $`L^h`$ to obtain a new graph $`\overline{G}`$ and decomposes $`\overline{G}`$ and the solution $`x`$ into $`\tfrac{10M}{\eps^2}`$-simple pairs $`(G^i,x^i), i \in [k]`$, as described in Section [sec:decomposition]. In the process, it also computes a set of links $`L^s`$ covering $`E^s`$ as defined in Lemma [lemma:decomposition].
Afterwards, we round each $`(G^i,x^i)`$ individually using the best of the following lemmas – the first one performs well for pairs where cross-links carry much of the cost and is based on the integrality of the odd-cut for instances with only cross- and up-links, while the second performs well for pairs where in-links carry much of the cost and is due to ; it relies on the bundle constraints. The output of the algorithm is then the union of $`L^h`$ and $`L^s`$ as computed by Lemma [lemma:decomposition], with the rounded solutions for all $`(G^i,x^i)`$.
Let $`x`$ be a feasible solution to the odd-cut for $`G`$. Let $`G`$ be rooted at any node $`r \in V[G]`$, then we can compute in polynomial time a solution $`S \subseteq L`$ covering $`E[G]`$ with cost at most
\begin{equation}
c(S) \leq 2c^\tp \xin + c^\tp \xcr.
\end{equation}
Proof. We define a new feasible solution for the odd-cut for $`G`$ by the applying to the following modification to all in-links with mass that are not up-links, all $`\ell \in \supp(x) \cap (\Lin \setminus \Lup)`$. Let $`\ell = uv`$ and $`w := \lca(uv)`$. Then we increase $`x_{uw}`$ and $`x_{wv}`$ by $`x_\ell`$, and set $`x_\ell := 0`$, and we call the resulting solution $`y`$.
By construction, we have $`\supp(x) \cap (\Lin \setminus \Lup) = \emptyset`$ – that means that the support of $`y`$ contains only up- and cross-links. At the same time, we know that $`c^\tp y \leq 2c^\tp \xin + c^\tp\xcr`$ since we doubled the mass of in-links. Notice finally that $`y`$ is feasible for the odd-cut , since the total mass of links covering each edge did not change by the modification.
Now resolve the odd-cut with the set of links restricted to $`\supp(y)`$. We know by Theorem [theorem:exact] that the odd-cut is integer for this set of links, allowing us to obtain an integral solution $`z`$ with $`c^\tp z \leq c^\tp y \leq 2c^\tp \xin + c^\tp\xcr`$. ◻
The idea of the rounding algorithm from is as follows. Given a $`\tfrac{10M}{\eps^2}`$-simple pair $`(H,x)`$ with a $`\beta`$-center $`r`$, we replace all cross-links $`\ell =uv`$ by two up-links $`ur`$ and $`rv`$. After that, we split $`H`$ at $`r`$ into trees $`H_1,\dots, H_m`$ such that $`r`$ is part of all trees. Since we no longer have cross-links, we can split $`x`$ into feasible solutions $`x^1, \dots, x^m`$ for the bundle for $`H_1,\dots,H_m`$ with $`\supp(x^i) \cap \supp(x^j) = \emptyset`$ if $`i \neq j`$. Every sub-tree $`H_1, \dots, H_m`$ of $`H`$ has at most $`\tfrac{10M}{\eps^2}+1`$ leaves (the root causes the additional $`+1`$); it can then be shown that every $`H_j`$ was a $`\ceil{\tfrac{28M}{\eps^2}}`$-bundle in the original graph, allowing us to use a bundle constraint to bound the rounding cost.
lemmaRoundingDavid Let $`x`$ be a feasible solution to the for $`G`$, with $`\gamma = \ceil{\tfrac{28M}{\eps^2}}`$. Given a $`\frac{10M}{\eps^2}`$-simple pair $`(H,x)`$ from the $`\tfrac{4M}{\eps^2}`$-thin decomposition of $`G`$ and $`x`$ that is rooted at a $`\frac{10M}{\eps^2}`$-center $`r`$, we can compute in polynomial time a set of links $`S \subseteq L`$ covering $`E[H]`$ with cost at most
\begin{equation}
c(S) \leq c^\tp \xin + 2c^\tp \xcr + |V^h \cap V[H]|
\end{equation}
with $`V^h`$ being the set of nodes created by the contractions of the links in $`L^h`$.
For completeness, we give a proof of this lemma in Appendix [appendix:rounding]. Together, these two lemmas allow us to proof the following theorem.
Proof. For constant $`\gamma`$, the can be solved in polynomial time due to Lemma [lem:solving_main_lp]. The decomposition can be performed in polynomial time, as can the computation of $`L^h`$ and $`L^s`$ in Lemma [lemma:decomposition] and the rounding schemes of Lemmas [lemma:rounding] and [lemma:david]. Since we cover all $`E[G^i]`$ as well as $`E^h`$ and $`E^s`$, we produce a feasible solution due to Observation [obs:feasibility].
What is left to analyze the approximation ratio of the algorithm. By Lemma [lemma:decomposition], we know that $`c(L^h \cup L^s) \leq O(\eps) \cdot \OPT`$. For a pair $`(H,x)`$ in the decomposition, we know by Lemmas [lemma:rounding] and [lemma:david] that the cost for its covering $`S`$ can be bounded by
\begin{aligned}
c(S) &\leq \min(2c^\tp \xin + c^\tp \xcr,c^\tp \xin + 2c^\tp \xcr + |V^h \cap V[H]|)\\
&\leq \frac{1}{2}\left(2c^\tp \xin + c^\tp \xcr + c^\tp \xin + 2c^\tp \xcr + |V^h \cap V[H]|\right)\\
&\leq \frac{3}{2} c^\tp \xin + \frac{3}{2} c^\tp \xcr + \frac{1}{2} |V^h \cap V[H]|.\\
\end{aligned}
Thus, the cost of the computed coverings $`S_1,\dots,S_k`$ for $`(G^1,x^1), \dots, (G^k,x^k)`$ is bounded by
\begin{aligned}
\sum_{i \in [k]} c(S_i) \leq |V^h| + \sum_{i \in [k]} \frac{3}{2}c^\tp x^i = |V^h| + \frac{3}{2} c^\tp x.
\end{aligned}
Since we have $`|V^h| \leq |L^h| \leq c(L^h) \in O(\eps) \OPT`$, the claim follows (every link costs at least $`1`$, and we need at least one link for a node in $`V^h`$). ◻
Decomposition of solutions
In this section, we will discuss how we use solutions of the to solve a instance approximately. We use the approach of , with two major differences:
-
We use the with $`\gamma = \ceil{\tfrac{28M}{\eps^2}}`$ instead of the bundle as basis for the decomposition and subsequent rounding.
-
The additional $`\set{0,\tfrac{1}{2}}`$-Chvátal-Gomory cuts in the compared to the bundle allow us to round instances with only up-links and cross-links without increasing cost. Since in-links can be split into two up-links, this yields a rounding that increases the cost incurred by in-links by a factor of 2 and leaves the cost of cross-links unchanged. This replaces a rounding based on a reduction to edge cover in that increases the cost of in-links by a factor of $`2\lambda`$ and the cost of cross-links by a factor $`\tfrac{4}{3}\tfrac{\lambda}{\lambda-1}`$ for some $`\lambda > 1`$.
For the decomposition, we will assume that we are given a shadow-complete instance $`(G,L,c)`$ and a fractional solution $`x`$ to the for $`\gamma = \ceil{\tfrac{28M}{\eps^2}}`$. An instance $`(G,L,c)`$ is shadow-complete, if for every link $`\ell \in L`$ all its shadows – links $`\ell'`$ with $`P_{\ell'} \subseteq P_\ell`$ – are also in $`L`$. 1 We assume that the cost of a shadow $`\ell'`$ of $`\ell`$ fulfills $`c_{\ell'} \leq c_\ell`$ – otherwise, we can set $`c_{\ell'} := c_\ell`$. The decomposition begins by computing a solution $`L^h`$ that covers all heavily-covered edges
E^h := \setc{e \in E}{x(\cov(e)) \geq \tfrac{2}{\eps}}
then contracts $`L^h`$ and obtains a new tree $`\overline{G}`$ that only contains edges $`e`$ with $`x(\cov(e)) < \tfrac{2}{\eps}`$. Then, we decompose $`\overline{G}`$ and $`x`$ by splitting repeatedly along $`\alpha`$-thin edges, which are defined as follows.
For an edge $`e = \set{u,v}`$, let $`G^u`$ and $`G^v`$ be the sub-trees of $`\overline{G}`$ that are obtained by deleting $`e`$, with $`G^u`$ containing $`u`$ and $`G^v`$ containing $`v`$, respectively. An edge $`e = \set{u,v}`$ is $`\alpha`$-thin with respect to $`x`$ for an $`\alpha \in \R_{\geq 0}`$ if the costs of links on both sides of the edge are at least $`\alpha`$ (we ignore links covering $`e`$ here):
\begin{equation}
\label{eq:alphathin}
\sum_{\ell \in L, \ell \subseteq V[G^u]} c_\ell x_\ell \geq \alpha, \quad \sum_{\ell \in L, \ell \subseteq V[G^v]} c_\ell x_\ell \geq \alpha.
\end{equation}
If $`\overline{G}`$ has an $`\alpha`$-thin edge $`e = \set{u,v}`$ with respect to $`x`$, we split $`\overline{G}`$ and $`x`$ along $`e`$:
We obtain $`x^u`$ and $`x^v`$ by splitting all links covering $`e`$ in $`x`$ into two parts, the part in $`G^u`$ and the part in $`G^v`$:
x^u_{\ell=pq}\ := \begin{cases}
x_\ell & \text{if } p,q \in V[G^u] \setminus\set{u},\\
0 & \text{if } p \in V[G^v]\text{ or }q \in V[G^v], \\
x_\ell + \displaystyle\sum\limits_{\substack{\ell' \in \cov(e)\\p\in\ell'}} \frac{c_{\ell'}}{c_\ell} x_{\ell'} & \text{if } p \in V[G^u], q = u,
\end{cases} \quad \text{for all } \ell \in L.
$`x^v`$ is defined symmetrically. Notice that $`\tfrac{c_{\ell'}}{c_\ell} \geq 1`$: this term ensures that the cost in the sub-trees does not decrease – this is important for the bundle constraints. Figure 1 depicts this splitting. If $`G^u`$ contains an $`\alpha`$-thin edge with respect to $`x^u`$, we repeat this procedure for $`G^u`$ and $`x^u`$, and similarly for $`G^v`$ and $`x^v`$ if $`G^v`$ contains an $`\alpha`$-thin edge with respect to $`x^v`$. At the end, we obtain a set of pairs $`\setc{(G^i,x^i)}{i \in [k]}`$ such that no $`G^i`$ contains an $`\alpha`$-thin edge with respect to $`x^i`$. We refer to the contraction of $`E^h`$ and the repeated edge-splitting as $`\alpha`$-thin edge decomposition of $`G`$ and $`x`$. Notice that we have $`\supp(x^i) \cap \supp(x^j) = \emptyset`$ for $`i \neq j`$ – this allows and us to round each $`(G^i,x^i)`$ independently from each other later on.
In the next lemma, we list important properties that the decomposition $`(G^i,x^i), i \in [k]`$ has – these properties were shown in for the bundle , but they also apply to the . One of the key properties is that the pairs $`(G^i,x^i)`$ are $`\beta`$-simple. A pair $`(G^i,x^i)`$ is called $`\beta`$-simple for a $`\beta \in \N`$, if there exists a $`\beta`$-center $`v \in V[G^i]`$ such that removal of $`v`$ decomposes $`G^i`$ into a forest of $`t`$ trees $`K_1,\dots, K_t`$ such that for all $`j \in [t]`$:
-
the cost of links in $`K_j`$ is bounded by $`\beta`$:
\begin{equation} \label{eq:beta1} \sum_{\ell \in L, \ell \subseteq V[K_j]} c_\ell x_\ell \leq \beta, \end{equation} -
$`K_j`$ has at most $`\beta`$ leaves.
Notice that the trees $`K_1,\dots,K_t`$ contain neither the $`\beta`$-center $`v`$ nor edges incident to $`v`$.
lemmaLemmaDecomposition Let $`(G,L,c)`$ be a instance with $`c_\ell \leq M`$ for all $`\ell \in L`$ and a constant $`M \in \R_{\geq 1}`$, and let $`x \in \Rnn^L`$ be a feasible solution to the for $`\gamma = \ceil{\tfrac{28M}{\eps^2}}`$. Let $`E^h`$ be the heavily-covered edges of $`G`$ with respect to $`x`$, $`L^h`$ be a set of links that covers $`E^h`$, $`\overline{G}`$ be the graph obtained from $`G`$ by contracting $`L^h`$ and $`(G^i,x^i)`$, $`i \in [k]`$ the $`\tfrac{4M}{\eps^2}`$-thin edge decomposition of $`\overline{G}`$ and $`x`$. Let $`E^s`$ be the set of edges along which we have split in the $`\tfrac{4M}{\eps^2}`$-thin edge decomposition. Then
-
$`x^i`$ is a feasible solution to the odd-cut for $`G^i`$, and it fulfills
\begin{equation} \label{eq:bundle-decomp} \sum_{\ell \in \cov(B)} c_\ell x_\ell \geq \OPT(B) \quad \text{ for all } B \in \mathcal{B}_\gamma^{G^i} \end{equation}where $`\mathcal{B}_\gamma^{G^i}`$ is the collection of all edge sets in $`E[G^i]`$ that are a $`\gamma`$-bundle in $`G`$.
-
$`\supp(x^i) \cap \supp(x^j) = \emptyset`$ for $`i,j \in [k], i \neq j`$,
-
every $`(G^i,x^i)`$, $`i \in [k]`$ is $`\tfrac{10M}{\eps^2}`$-simple and we have
\begin{equation} \label{eq:beta2} \sum_{\ell \in L, \ell \subseteq V[K_j]} c_\ell x_\ell \leq \frac{4M}{\eps^2} \end{equation}for the trees $`K_j, j\in [t]`$ that are created by removing a $`\beta`$-center $`v`$ from $`G^i`$.
-
$`\sum_{i \in [k]} c^\tp x^i \leq \left(1+\eps\right)c^\tp x`$,
-
we can efficiently compute sets of links $`L^h, L^s \subseteq L`$ covering $`E^h, E^s`$, respectively, with $`c(L^h) \leq \eps c^\tp x`$ and $`c(L^s) \leq O(\eps) \sum_{i\in [k]} c^\tp x^i`$.
The proof of this lemma can be found in Appendix [appendix:decomposition], and is based on . Properties 1 and 2 allow us to round each $`x^i`$ individually, and property 3 makes it possible to employ a rounding algorithm from . Properties 4 and 5 ensure that the decomposition costs us only a factor of $`(1+O(\eps))`$, and that the edges not contained in a $`G^i`$ can also be covered at cost $`O(\eps)\OPT`$. Thus, the approximation guarantee is dominated by how well we can round the individual $`x^i`$ solutions.
A Stronger
In this section, we introduce the which strengthens the bundle used in . We use our in the algorithmic framework described in Section [sec:decomposition] to obtain better approximation guarantees. The additional constraints that we need are Chvátal-Gomory cuts obtained from the cut .
Deriving and Separating the Odd-Cut Constraints
Derivation.
Consider an of the form $`\min \{c^\intercal x \mid Ax \geq b,\ x \in \Z^n\}`$, with $`A \in \R^{m \times n}`$ and $`b \in \R^{m}`$. A Chvátal-Gomory cut is a constraint of the form $`\lambda^\intercal A x \geq \lceil \lambda^\intercal b \rceil`$, where the vector of multipliers $`\lambda \in \R^m_{\geq 0}`$ is chosen in such a way that $`\lambda^\intercal A \in \Z^n`$. Clearly, any such constraint is valid for the integer solutions of the . It is well known that the cuts obtained for $`\lambda \in [0,1)^m`$ imply the cuts obtained for larger multipliers, hence one can assume $`\lambda \in [0,1)^m`$ without loss of generality. (A proof of this fact and more background on the Chvátal-Gomory cuts can be found, for instance, in .) In case we restrict further $`\lambda`$ to be in $`\{0,\frac{1}{2}\}^m`$, we obtain a $`\{0,\frac{1}{2}\}`$-Chvátal-Gomory cut . These cuts are a proper specialization of the Chvátal-Gomory cuts, and are precisely the cuts that we use here.
A $`\{0,\frac{1}{2}\}`$-Chvátal-Gomory cut for the cut is any constraint of the form
\begin{equation}
\label{eq:0-1/2-cut} \sum_{e \in E[G]} \lambda_e x(\cov(e)) + \sum_{\ell \in L} \mu_\ell x_\ell \geq \Big\lceil\sum_{e \in E[G]}\lambda_e\Big\rceil
\end{equation}
where $`\lambda \in \{0,\frac{1}{2}\}^{E[G]}`$ and $`\mu \in \{0, \frac{1}{2}\}^{L}`$ are such that the coefficients in the left-hand side are all integral. Notice that for any fixed $`\lambda`$, there is a unique $`\mu`$ that achieves this.
Let $`K := \supp(\lambda) = \{e \in E[G] \mid \lambda = \frac{1}{2}\}`$. Since $`G`$ is a tree, there exists a (not necessarily connected) set $`S \subseteq V[G]`$ such that $`K = \delta_G(S)`$. Notice that the right-hand side of [eq:0-1/2-cut] is $`\lceil |\delta_G(S)|/2 \rceil`$, and that the cut is redundant whenever $`|\delta_G(S)|`$ is even.
Let $`\pi(S)`$ denote the multiset of links $`\ell`$ such that $`P_\ell`$ intersects $`\delta_G(S)`$, in which the multiplicity of $`\ell`$ is defined as $`\lceil \frac{1}{2} |P_\ell \cap \delta_G(S)| \rceil`$ (see Figure 2). Now, assuming that $`|\delta_G(S)|`$ is odd, we can rewrite the constraint [eq:0-1/2-cut] as
\begin{equation}
\label{eq:0-1/2-cut_pi}
x(\pi(S))\geqslant \frac{|\delta_G(S)|+1}2\,.
\end{equation}
Let $`\mathcal{S}`$ denote the collection of all sets $`S`$ such that $`|\delta_G(S)|`$ is odd, and let the be the resulting from the cut after adding the odd-cut constraint [eq:0-1/2-cut_pi] for each set $`S \in \mathcal{S}`$:
\begin{lp}{min}{\sum_{\ell \in L} c_\ell x_\ell}
\lpconstraint{x(\pi(S))}{\geq}{\frac{|\delta_G(S)|+1}{2}}{S \in \mathcal{S}}
\finallpconstraint{x_\ell}{\geq}{0}{\ell\in L}
\end{lp}
We remark that the first set of constraints includes in particular the covering constraint $`x(\cov(e)) \geq 1`$ for each tree edge $`e \in E[G]`$. Indeed, if $`S`$ denotes any of the two connected components arising when $`e`$ is deleted from $`G`$, we have $`\delta_G(S) = \{e\}`$ and thus $`S \in \mathcal{S}`$. The corresponding odd-cut constraint is the covering constraint.
Separation.
A nice fact that is key to our approach is that separation of the $`\{0,\frac{1}{2}\}`$-Chvátal-Gomory cuts arising from an $`\min \{c^\intercal x \mid Ax \geq b,\ x \in \Z^n\}`$ (here we take $`A \in \Z^{m \times n}`$ and $`b \in \Z^m`$) can be done in polynomial time whenever the matroid represented over $`GF(2)`$ by $`(I\,\bar{A})`$ is graphic (or co-graphic), where $`\bar{A} := A \pmod{2}`$ is the parity matrix of $`A`$. This was proved by Caprara and Fischetti , and implies directly that the constraints of the can be separated in polynomial time.
A more direct way to prove this is to rewrite inequality [eq:0-1/2-cut_pi] as the $`T`$-cut constraint 2 $`x(\delta_L(S)) + y(\delta_G(S)) \geq 1`$ after introducing the slack variables $`y_e := x(\cov(e)) - 1 \geq 0`$ for $`e \in E[G]`$, and then solve the separation problem with the algorithm of Padberg and Rao .
To see this, let $`S \in \mathcal{S}`$ and let $`H`$ be the graph obtained from adding all links in $`L`$ to $`G`$. Assuming that $`y_e = x(\cov(e)) - 1`$ for all $`e \in E[G]`$, we can rewrite the odd-cut constraint [eq:0-1/2-cut_pi] as:
\begin{eqnarray*}
x(\pi(S)) & \geq & \frac{|\delta_G(S)| + 1}2\\
\iff \frac12\sum_{e \in \delta_G(S)}x(\cov(e)) + \frac12
x(\delta_L(S)) & \geq & \frac{|\delta_G(S)| + 1}2\\
\iff \sum_{e \in \delta_G(S)} \underbrace{(x(\cov(e)) - 1)}_{=y_e} +
x(\delta_L(S)) & \geq & 1\\
\iff (x,y)(\delta_H(S)) & \geq & 1\,.
\end{eqnarray*}
Above, the second form is the original expression of the odd-cut constraint as a $`\{0,\frac{1}{2}\}`$-Chvátal Gomory cut, see [eq:0-1/2-cut].
Exactness of the Odd-Cut when all the In-links are Up-links
An integer matrix $`M`$ is said to be the incidence matrix of a bidirected graph if $`\sum_i |M_{ij}| \leq 2`$ for every fixed column index $`j`$. A binet matrix is any matrix of the form $`B = S^{-1} R`$ where $`M = (S\, R)`$ is the incidence matrix of a bidirected graph with full row-rank and $`R`$ is a basis of $`M`$. Binet matrices are a generalization of network matrices and were introduced by Appa and Kotnyek .
Appa proved the following result, extending results of Edmonds and Johnson for incidence matrices of bidirected graphs.
For every binet matrix $`A \in \Z^{m \times n}`$ and every vector $`b \in \Z^{m}`$, the integer hull of the polyhedron $`P := \setc{x \in \R^n}{Ax \geq b,\ x \geq 0}`$ is described by its .
Pick any root $`r \in V[G]`$ in tree $`G`$. We are now ready to recall and prove Theorem [theorem:exact], which provides a key property of the that we use in our approach to approximate .
Proof. Let $`(G,L,c)`$ denote the instance and let $`r`$ denote the chosen root. For a directed edge $`e = (u,v)`$, let $`z = z(e) \in \R^{V[G] \setminus \{r\}}`$ denote the truncated incidence vector of $`e`$ (since the row corresponding to root $`r`$ is removed), defined as $`z_a := 1`$ if $`a = u`$, $`z_a := -1`$ if $`a = v`$ and $`z_a = 0`$ otherwise, for $`a \in V[G] \setminus \{r\}`$.
Direct all the edges of $`G`$ away from the root and let $`R \in \R^{(V[G] \setminus \{r\}) \times E[G]}`$ denote the (truncated) incidence matrix of the resulting directed graph $`\overrightarrow{G}`$. Every directed edge $`e = (u,v)`$ of $`\overrightarrow{G}`$ has a corresponding column $`z(e)`$ in $`R`$.
Now, we define a matrix $`S \in \R^{(V[G] \setminus \{r\}) \times L}`$ encoding the links of instance. Each link $`\ell = uv\in L`$ has a corresponding column $`S_{\star \ell}`$ in $`S`$. If $`\ell = uv`$ is an up-link with $`u = \lca(uv)`$, we let $`S_{\star \ell}`$ be the incidence vector $`z(u,v)`$ of the directed edge $`(u,v)`$. If $`\ell = uv`$ is a cross-link, then $`u, v \neq r`$ and the corresponding column $`S_{\star \ell}`$ has $`S_{a \ell} := 1`$ if $`a \in \{u,v\}`$ and $`S_{a \ell} := 0`$ otherwise.
Consider the matrix $`M = (S\, R)`$. Since $`\sum_{a \neq r} |S_{a\ell}| \leq 2`$ for all $`\ell \in L`$ and $`\sum_{a \neq r} |R_{ae}| \leq 2`$ for all $`e \in E[G]`$, we see that $`M`$ is the incidence matrix of a bidirected graph. Moreover, $`M`$ has $`R`$ as a basis. Thus $`B := R^{-1} S`$ is a binet matrix. We claim that $`B \in \R^{(V[G] \setminus \{r\}) \times L} \cong \R^{E[G] \times L}`$ is the constraint matrix $`A`$ of the cut $`\min \{c^\intercal x \mid Ax \geq 1,\ x \geq 0\}`$.
In other words, we claim that $`A = R^{-1}S`$, or equivalently $`R A = S`$. This actually is clear: For every link $`\ell = uv`$, the sum of the columns $`z(e)`$ of $`R`$ that correspond to the tree edges $`e = (u,v)`$ that are in $`P_\ell`$ is the corresponding column of $`S`$.
By Theorem [theorem:binet_zhcgcuts], the cut becomes integral after one round of $`\{0,\frac{1}{2}\}`$-Chvátal Gomory cuts. That is, the is integral. ◻
We point out that for these instances, there is a combinatorial algorithm that finds an optimum solution in strongly polynomial time. This follows from our proof of Theorem [theorem:exact] and an algorithm of Edmonds and Johnson , see also .
The Odd-Cut Bundle
As its name indicates, the contains all of the constraints of the and additionally the bundle constraints for a constant $`\gamma \in \N`$. As before, let $`\mathcal{S}`$ denote the collection of all sets $`S \subseteq V[G]`$ such that $`|\delta_G(S)|`$ is odd, let $`\mathcal{B}_\gamma`$ be the set of all $`\gamma`$-bundles, and let $`\OPT(B)`$ the cost of an integral optimal solution for the instance obtained from contracting all edges not in $`B`$ (with respect to the given costs $`c_\ell`$). Then, the is given by:
\begin{lp}{min}{\sum_{\ell \in L} c_\ell x_\ell}
\lpconstraint{x(\pi(S))}{\geq}{\frac{|\delta_G(S)|+1}{2}}{S \in \mathcal{S}}
\lpconstraint{\sum_{\ell \in \cov(B)} c_\ell x_\ell}{\geq}{\OPT(B)}{B \in \mathcal{B}_\gamma}
\finallpconstraint{x_\ell}{\geq}{0}{\ell\in L}
\end{lp}
From Lemma [lemma:constantleaves] and the discussion above, we obtain the following result.
For any constant $`\gamma \in \N`$, the can be solved in polynomial time.
Proof of Lemma <a href="#lemma:decomposition" data-reference-type=“ref”
data-reference=“lemma:decomposition”>[lemma:decomposition]
In this section, we prove that the properties of the decomposition in Lemma [lemma:decomposition] also hold if we use the as a basis for the decomposition instead of the bundle . The proof ideas are all due to . Notice that we use $`\gamma = \ceil{\frac{28M}{\eps^2}}`$ and $`\beta = \frac{10M}{\eps^2}`$ instead of $`\gamma = \ceil{\frac{56M}{\eps^2}}`$ and $`\beta = \frac{12M}{\eps^2}`$ in ; this is the result of slightly improved estimations. In particular, we use that:
-
Links covering an edge $`e = \set{u,v}`$ can cover at most one leaf in the subtrees $`G^u`$ and $`G^v`$. This argument allows us to reduce $`\beta`$ from $`\frac{12M}{\eps^2}`$ to $`\frac{10M}{\eps^2}`$.
-
The cost of links completely in a sub-tree of a $`\tfrac{10M}{\eps^2}`$-simple pair is actually bounded by $`\tfrac{4M}{\eps^2}`$, even though the pair might not be $`\tfrac{4M}{\eps^2}`$-simple due to the number of leaves of the sub-trees.
-
The previous point allows us bound $`|V^h|`$ by $`\tfrac{6M}{\eps^2}`$; together with the improved $`\beta = \frac{10M}{\eps^2}`$ we obtain $`\gamma = \ceil{\frac{28M}{\eps^2}}`$.
Details can be found in the proof below.
Proof. We now prove the five properties in the order in which they were stated.
-
Feasibility: We start with a solution $`x`$ that is feasible for the for $`G`$, then we contract the links in $`L^h`$ to obtain a new graph $`\overline{G}`$. Since we do not modify the solution during the contraction, $`x`$ fulfills [eq:bundle-decomp] for $`B \in \mathcal{B}^{\overline{G}}_\gamma`$ by definition of $`\mathcal{B}^{\overline{G}}_\gamma`$, because $`\mathcal{B}_\gamma^{G'}`$ is the set of edges in $`G'`$ that are $`\gamma`$-bundles in $`G`$.
The odd-cut constraints for $`E[\overline{G}]`$ are
\begin{equation*} \sum_{\ell \in L} \ceil{\tfrac{1}{2}|P_\ell \cap \delta_{\overline{G}}(S)|} x_\ell \geq \frac{|\delta_{\overline{G}}(S)|+1}{2} \qquad \text{for all } S \in \mathcal{S}_{\overline{G}} \end{equation*}with $`\mathcal{S}_{\overline{G}}`$ being the collection of all sets $`S`$ such that $`|\delta_{\overline{G}}(S)|`$ is odd. Let $`\mathcal{S}_{G}`$ being the collection of all sets $`S`$ such that $`|\delta_{G}(S)|`$ is odd. For a set $`S \in \mathcal{S}_{\overline{G}}`$, let $`S' \in \mathcal{S}_{G}`$ be the set with $`\delta_{\overline{G}}(S) = \delta_G(S')`$. The existence of $`S'`$ is guaranteed since $`\overline{G}`$ is a contraction of the tree $`G`$. Since we do not modify the solution during the contraction, the existence of $`S'`$ for every $`S`$ guarantees that the odd-cut constraints are fulfilled.
Now assume that we have a graph $`H`$ and a solution $`x`$ which is feasible for the for $`G`$, and fulfills [eq:bundle-decomp] for all $`B \in \mathcal{B}^H_\gamma`$. We now show that splitting an edge $`e = \set{u,v} \in E[H]`$ results in two pairs $`(H^u,x^u)`$, $`(H^v,x^v)`$ such that $`x^u`$, $`x^v`$ are feasible for the for $`H^u`$, $`H^v`$ and fulfill [eq:bundle-decomp] for all $`B \in \mathcal{B}^{H^u}_\gamma`$ and $`B \in \mathcal{B}^{H^v}_\gamma`$, respectively.
Recall that splitting an edge $`e = \set{u,v}`$ results in $`x^u`$ (and analogously $`x^v`$) in the following way:
x^u_{\ell=pq}\ := \begin{cases} x_\ell & \text{if } p,q \in V[G^u] \setminus\set{u},\\ 0 & \text{if } p \in V[G^v]\text{ or }q \in V[G^v], \\ x_\ell + \displaystyle\sum\limits_{\substack{\ell' \in \cov(e)\\p\in\ell'}} \frac{c_{\ell'}}{c_\ell} x_{\ell'} & \text{if } p \in V[G^u], q = u, \end{cases} \quad \text{for all } \ell \in L.This definition differs from in the additional $`\frac{c_{\ell'}}{c_\ell}`$ term that ensures that bundle constraints remain feasible in the case that the shadows to which mass is shifted to are cheaper than the links where the mass originates from. Due to this, we have that for every edge $`e \in E[G^u]`$
\sum_{\ell \in \cov(e)} x^u_\ell \geq \sum_{\ell \in \cov(e)} x_\ell, \quad \text{and } \sum_{\ell \in \cov(e)} c_\ell x^u_\ell \geq \sum_{\ell \in \cov(e)} c_\ell x_\ellby definition of the splitting; the same applies to any edge set $`F \subseteq E[G^u]`$:
\sum_{\ell \in \cov(F)} x^u_\ell \geq \sum_{\ell \in \cov(F)} x_\ell, \quad \text{and}\quad \sum_{\ell \in \cov(F)} c_\ell x^u_\ell \geq \sum_{\ell \in \cov(F)} c_\ell x_\ell.This implies that the left hand sides of [eq:bundle-decomp] can only get larger in $`G^u`$, while the right hand sides stay the same.
For the odd-cut constraints, consider a set $`S \subseteq V[G^u]`$ such that $`|\delta_{G^u}(S)|`$ is odd, and let $`\mathcal{S}_{G^u}`$ be the collection of these sets. The important fact is here that the mass of links covering an edge in $`E[G^u]`$ does not decrease. We define $`\alpha_{\ell,S} := \ceil{\tfrac{1}{2}|P_\ell \cap \delta_{G^u}(S)|}`$ for all $`\ell \in L, S \in \mathcal{S}_{G^u}`$, and note that $`\alpha_{\ell,S} = \alpha_{\ell',S}`$ for $`\ell = pu \in L`$, $`\ell' = pq \in \cov(e)`$ yielding
\begin{aligned} &\ \sum_{\ell \in L} \ceil{\tfrac{1}{2}|P_\ell \cap \delta_{G^u}(S)|} x^u_\ell \\[-4mm] =&\ \sum_{\ell \in L, \ell \subseteq V[G^u] \setminus\set{u}} \alpha_{\ell,S}\cdot x_\ell + \sum_{\ell = pu \in L} \alpha_{\ell,S} \left(x_\ell+\sum_{\substack{\ell' \in\cov(e)\\p\in\ell'}} \frac{c_{\ell'}}{c_\ell}x_{\ell'} \right) + \sum_{\ell \in L, \ell \in \cov(e)} \alpha_{\ell,S} \cdot 0\\ \geq&\ \sum_{\ell \in L, \ell \subseteq V[G^u] \setminus\set{u}} \alpha_{\ell,S}\cdot x_\ell + \sum_{\ell = pu \in L} \alpha_{\ell,S} x_\ell + \sum_{\ell = pu \in L} \sum_{\substack{\ell' \in\cov(e)\\p\in\ell'}} \alpha_{\ell',S}\cdot x_{\ell'}\\ \geq&\ \sum_{\ell \in L, \ell \subseteq V[G^u] \setminus\set{u}} \alpha_{\ell,S}\cdot x_\ell + \sum_{\ell = pu \in L} \alpha_{\ell,S} x_\ell + \sum_{p \in V[G^u]} \sum_{\ell = pq \in\cov(e)} \alpha_{\ell,S}\cdot x_{\ell}\\ \geq&\ \sum_{\ell \in L, \ell \subseteq V[G^u] \setminus\set{u}} \alpha_{\ell,S}\cdot x_\ell + \sum_{\ell = pu \in L} \alpha_{\ell,S} x_\ell + \sum_{\ell \in L, \ell \in \cov(e)} \alpha_{\ell,S}\cdot x_{\ell}\\ \geq&\ \frac{|\delta_{G^u}(S)|+1}{2} \qquad \text{for all } S \in \mathcal{S}_{G^u}. \end{aligned}This proves adherence to the odd-cut constraints. Analogous statements apply to $`G^v`$ and $`x^v`$. This argumentation can now be used inductively on the splitting process – this completes the proof of this property.
-
Disjointness: By definition of the edge splitting, we know that $`V[G^i]`$ and $`V[G^j]`$ are disjoint for all $`i \neq j, i,j \in [k]`$. Furthermore, for any link $`\ell \in \supp(x^i)`$, we have $`\ell \subseteq V[G^i]`$, proving this claim.
-
Simplicity: Here, we have to show that every $`(G^i,x^i)`$ is $`\tfrac{10M}{\eps^2}`$-simple, which requires that there is a node $`v`$ for every $`G^i`$ whose removal decomposes $`G^i`$ into trees $`K_1,\dots,K_t`$ with two properties:
-
$`K_j`$ fulfills $`\sum_{\ell \in L, \ell \subseteq V[K_j]} c_\ell x_\ell \leq \frac{10M}{\eps^2}`$ for all $`j \in [t]`$,
-
$`K_j`$ has at most $`\frac{10M}{\eps^2}`$ leaves for all $`j \in [t]`$.
Furthermore, we have to prove that
\begin{equation*} \sum_{\ell \in L, \ell \subseteq V[K_j]} c_\ell x_\ell \leq \frac{4M}{\eps^2} \quad \text{for all } j \in [t]. \end{equation*}Notice that this implies the first property of $`(G^i,x^i)`$ being $`\tfrac{10M}{\eps^2}`$-simple.
Let $`(H,x) := (G^i,x^i)`$ for some $`i \in [k]`$. We know that every edge $`e = \set{u,v} \in E[H]`$ satisfies
\begin{equation} \label{eq:alphaviol} \sum_{\ell \in L, \ell \subseteq V[H^u]} c_\ell x_\ell < \frac{4M}{\eps^2} \quad \text{ and / or } \quad \sum_{\ell \in L, \ell \subseteq V[H^v]} c_\ell x_\ell < \frac{4M}{\eps^2}. \end{equation}We distinguish two cases:
-
In this case, there exists an edge $`e = \set{u,v} \in E[H]`$ with:
\begin{equation*} \label{eq:decompcase1} \sum_{\ell \in L, \ell \subseteq V[H^u]} c_\ell x_\ell < \frac{4M}{\eps^2} \quad \text{ and }\quad \sum_{\ell \in L, \ell \subseteq V[H^v]} c_\ell x_\ell < \frac{4M}{\eps^2}. \end{equation*}Thus, choosing $`u`$ or $`v`$ decomposes $`H`$ into subtrees $`K_1,\dots,K_t`$ with
\begin{equation*} \label{eq:decompcase1b} \sum_{\ell \in L, \ell \subseteq V[K_i]} c_\ell x_\ell < \frac{4M}{\eps^2} \end{equation*}for all $`i \in [t]`$.
-
Otherwise, we have for all edges $`e = \set{u,v} \in E[H]`$, that
\begin{equation} \text{either } \sum_{\ell \in L, \ell \subseteq V[H^u]} c_\ell x_\ell < \frac{4M}{\eps^2} \quad \text{ or } \quad \sum_{\ell \in L, \ell \subseteq V[H^v]} c_\ell x_\ell < \frac{4M}{\eps^2}. \end{equation}We orient every edge $`e = \set{u,v}`$ based on whether $`\sum_{\ell \in L, \ell \subseteq V[H^u]} c_\ell x_\ell < \frac{4M}{\eps^2}`$ holds: if yes, we orient it from $`u`$ to $`v`$, otherwise from $`v`$ to $`u`$. Since $`H`$ is a tree, there is a node $`v \in V[H]`$ with $`\outdeg(v) = 0`$. Removing $`v`$ again decomposes $`H`$ into sub-trees $`K_1,\dots,K_t`$ with
\begin{equation*} \label{eq:decompcase1c} \sum_{\ell \in L, \ell \subseteq V[K_i]} c_\ell x_\ell < \frac{4M}{\eps^2} \end{equation*}for all $`i \in [t]`$.
This completes the first part of the proof. Now we bound the number of leaves for a tree $`K_j, j \in [t]`$.
Now consider a tree $`K_j`$, $`j \in [t]`$, and let $`e`$ be the edge that connected $`K_j`$ to the $`\beta`$-center in $`G^i`$ that was removed. The number of leaves that a tree $`K_j`$ has can be bounded by the following argument. Every link $`\ell = uv`$ has cost at least 1 and can contribute to the covering of at most 2 leaf edges of $`K_j`$ if $`u,v \in V[K_j]`$; if either $`u`$ or $`v`$ is not in $`V[K_j]`$ then $`\ell`$ can cover at most one of $`K_j`$ leaf edge. Finally, if both $`u,v \not\in V[K_j]`$ it cannot cover leaf edges of $`K_j`$ at all.
We get the following bound for the number of leaves $`|\leaves(K_j)|`$ of $`K_j`$ (since $`e \not\in E^h`$):
|\leaves(K_j)| \leq 2\cdot\sum_{\ell \in L, \ell \subseteq V[K_j]} x_\ell + \sum_{\ell \in \cov(e)} x_\ell \leq 2\cdot\frac{4M}{\eps^2} + \frac{2}{\eps} \leq \frac{10M}{\eps^2}.This completes the proof of this property.
-
-
Cost Increase: We have to show that
\sum_{i \in [k]} c^\tp x^i \leq \left(1+\eps\right)c^\tp x.Consider an edge-splitting of a pair $`(H,x) := (G^i,x^i)`$ along an edge $`e = \set{u,v}`$ into $`H^u, H^v, x^u`$ and $`x^v`$. We have $`c^\tp x \leq c^\tp x^u + c^\tp x^v + \sum_{\ell \in \cov(e)} c_\ell x_\ell`$ by definition of the edge-splitting. Thus, each splitting increases the cost by $`\sum_{\ell \in \cov(e)} c_\ell x_\ell < \tfrac{2}{\eps}M`$ due to our link cost bound $`M`$ and the fact that all edges with a mass of more than $`\tfrac{2}{\eps}`$ were contracted. Since we end up with $`k`$ pairs, we have $`k-1`$ splits and the total increase in cost is bounded by $`(k-1)\tfrac{2M}{\eps}`$:
\sum_{i\in [k]} c^\tp x^i \leq c^\tp x + (k-1)\frac{2M}{\eps} \leq c^\tp x + k\frac{2M}{\eps}.On the other hand, we split along $`\tfrac{4M}{\eps^2}`$-thin edges and have therefore
c^\tp x^i = \sum_{\ell \in L, \ell \subseteq V[G^i]} c_\ell x_\ell \geq \frac{4M}{\eps^2} \quad \text{ for all } i \in [k].Combined, we now have for $`\eps \leq 1`$:
\begin{aligned} \sum_{i\in [k]} c^\tp x^i \leq c^\tp x + k\frac{2M}{\eps} &= \left(1+\frac{k\frac{2M}{\eps}}{c^\tp x}\right) c^\tp x\\ &\leq \left(1+\frac{k\frac{2M}{\eps}}{ \sum_{i\in [k]} c^\tp x^i-k\frac{2M}{\eps}}\right) c^\tp x\\ &\leq \left(1+\frac{k\frac{2M}{\eps}}{k\frac{4M}{\eps^2}-k\frac{2M}{\eps}}\right) c^\tp x\\ &= \left(1+\frac{1}{\frac{2}{\eps}-1}\right) c^\tp x = \left(1+\frac{\eps}{2-\eps}\right) c^\tp x \leq \left(1+\eps\right) c^\tp x.\\ \end{aligned}That completes this part of the proof.
-
Remaining Edges: We need to show the existence of link sets $`L^h, L^s`$ that cover $`E^h`$ and $`E^s`$, respectively, with $`c(L^h) \leq \eps c^\tp x`$ and $`c(L^s) \leq O(\eps) \sum_{i\in [k]} c^\tp x^i`$.
Consider the graph $`G'`$ obtained by contracting all edges in $`E[G] \setminus E^h`$. Since $`E^h := \setc{e \in E[G]}{x(\cov(e)) \geq \frac{2}{\eps}}`$ we know that $`y := \tfrac{\eps}{2}x`$ is a feasible solution to the cut for $`G'`$. The cut is known to have an gap of at most 2 due to various 2-approximation algorithms, or ; applying this to $`y`$ yields a set of links $`L^h`$ with
c(L^h) \leq 2 c^\tp y = 2\tfrac{\eps}{2} c^\tp x = \eps c^\tp x.This shows the first part of the claim. For the second part, remember that all edges $`e \in E^s`$ used for splitting were $`\tfrac{4M}{\eps^2}`$-thin edges,
c^\tp x^i = \sum_{\ell \in L, \ell \subseteq V[G^i]} c_\ell x_\ell \geq \frac{4M}{\eps^2} \quad \text{for all } i\in [k].Since $`|E^s| = k-1`$, any inclusion-wise minimal covering of $`E^s`$ costs at most $`(k-1)M`$ and we have
(k-1)M = \frac{(k-1)M}{\sum_{i\in [k]} c^\tp x^i}\sum_{i\in [k]} c^\tp x^i \leq \frac{(k-1)M}{k\tfrac{4M}{\eps^2}}\sum_{i\in [k]} c^\tp x^i \leq \frac{\eps^2}{4}\sum_{i\in [k]} c^\tp x^i \in O(\eps) \sum_{i\in [k]} c^\tp x^i .
◻
Introduction
The tree augmentation problem (weighted or unweighted) is a fundamental and intensively studied problem in the area of network design, see for example the surveys by Khuller and Kortsarz and Nutov . While already in the unweighted case the problem is known to be -hard, the best algorithms for and achieve approximation factors of $`2`$ and $`\tfrac{3}{2}`$ respectively. One of the main open questions about these problems is to improve the quality of approximation algorithms.
Adjiashvili recently managed to push the approximation guarantee for below $`2`$ in case the link costs are bounded by a constant $`M`$, which was the first improvement in over 35 years that did not restrict the structure of the tree or the links. His algorithm is based on an that strengthens the standard for the problem. Letting $`\cov(e)`$ denote the set of links connecting distinct connected components of $`G - e`$, for each tree edge $`e \in E`$, the standard for is
\begin{alignat}
{4}
& \text{min} & & \sum_{\ell \in L} c_\ell &x_\ell & & \label{eq:cut1}\\
& \text{s.\,t.} & \quad & \sum_{\mathclap{\ell \in \cov(e)}} &x_\ell &\geq 1 & \qquad &\text{for all } e\in E,\label{eq:cut2}\\
& & & &x_\ell &\geq 0 & \qquad &\text{for all } \ell \in L.\label{eq:cut3}
\end{alignat}
This is known as the cut .
Our results.
We add to the cut all its , thus obtaining a new for that we call the . The is key to our approach. It has three extremely useful properties:
-
One can solve the efficiently, even though separating is -hard in general .
-
The is compatible with the decomposition approach of to split the given instance and solution into well-structured independent instances together with their own local solutions.
-
The is exact if for a certain choice of root $`r \in V[G]`$, every link $`\ell`$ connects either two different connected components of $`G - r`$ (in which case $`\ell`$ is called a cross-link) or some node of $`G`$ to one of its ancestors (in which case $`\ell`$ is called an up-link).
We prove the last property by establishing that the constraint matrix of the cut is an integral binet matrix. These matrices were introduced by Appa and Kotnyek as a generalization of network matrices. Relying on earlier work by Edmonds and Johnson , Appa proved that the integer hull of polyhedra of the form $`\setc{x}{Ax \geq b,\ x \geq 0}`$ can be described by whenever $`A`$ is an integral binet matrix and $`b`$ is an integer vector. This results in the following theorem.
theoremexact The odd-cut is integral for instances that contain only cross- and up-links.
Although the alone might be sufficient to obtain a $`\tfrac{3}{2}`$-approximation for (or maybe even better approximation), we combine the odd-cut with the bundle constraints from , resulting in the . This last is the one that we use in our algorithm.
We follow the decomposition approach of . After splitting the given instance and its optimum solution into independent rooted instances and corresponding solutions (respectively), at an extra cost of $`\varepsilon \OPT`$, Adjiashvili applies to each one of the local instances two distinct procedures producing feasible solutions, whose cost is bounded in terms of the local solution.
One of the two procedures of produces an integer solution of cost at most $`c^\intercal x^\mathrm{in} + 2 c^\intercal x^\mathrm{cr} + \varepsilon \OPT`$, where $`c^\intercal x^\mathrm{in}`$ is the local cost on in-links (defined as all the links that are not cross-links) and $`c^\intercal x^\mathrm{cr}`$ is the local cost on cross-links. This is the part of the analysis where bundle constraints are used. We keep this procedure as is in our algorithm.
Using Theorem [theorem:exact], we improve the other procedure of to obtain an integer solution of cost at most $`2 c^\intercal x^\mathrm{in} + c^\intercal x^\mathrm{cr}`$. This gives a significant improvement in the approximation factor since combining both procedures, we can construct an integer solution in each of the local instances of cost at most
\begin{align*}
&\ \min\left\{c^\intercal x^\mathrm{in} +
2 c^\intercal x^\mathrm{cr} + \delta,
2 c^\intercal x^\mathrm{in} +
c^\intercal x^\mathrm{cr}\right\}\\
\leq &\
\frac12 \left(c^\intercal x^\mathrm{in} +
2 c^\intercal x^\mathrm{cr} + \delta\right)
+ \frac12 \left(2 c^\intercal x^\mathrm{in} +
c^\intercal x^\mathrm{cr}\right)\\
\leq&\ \frac32 \left(c^\intercal x^\mathrm{in} +
c^\intercal x^\mathrm{cr}\right) + \delta,
\end{align*}
where $`\delta`$ is a small quantity whose sum across the local instances is at most $`\varepsilon \OPT`$. This yields our main result.
theoremOnePointHalf For every fixed $`\eps > 0`$ and $`M \in \R_{\geq 1}`$, there exists an -based polynomial time $`\tfrac{3}{2}+\eps`$-approximation algorithm for with link costs in $`[1,M]`$.
This result is the best known for with link costs bounded by a constant, and a significant improvement over the previously known $`\approx 1.96418+\eps`$-approximation . For , it is the best -based result and matches the results of the -based algorithm of , while only being the $`\eps`$-term worse than the best overall algorithm . Finally, we point out that while the approximation factors are improving, the proofs are actually getting simpler, which we see as another indication of the power behind our approach.
Related work.
Frederickson and JáJá showed that is -hard even if the tree has constant diameter and the link costs are either 1 or 2. For , Cheriyan proved that the problem is -hard even if the links form a cycle on the leaves of the tree. Kortsarz, Krauthgamer and Lee then showed that even is hard, meaning that these problems have no , unless $`\PT = \NP`$.
For , there are several 2-approximations known, the first given by Frederickson and Jájá and then simplified by Khuller and Thurimella . The primal-dual approach of Goemans and the iterative rounding algorithm of Jain also give a 2-approximation for this problem. Until recently, the factor of 2 was best known for general trees. Adjiashvili then managed to give a $`\delta+\eps`$-approximation for the case where all link costs are bound by a constant, for any small $`\eps > 0`$ and $`\delta = \tfrac{8(23+3\sqrt{5})}{121} \approx 1.96418`$.
For , the best known approximation is a $`\tfrac{3}{2}`$-approximation algorithm by Kortsarz and Nutov , which is purely combinatorial. There is also a $`\tfrac{3}{2}+\eps`$-approximation for any $`\eps > 0`$, given by Cheriyan and Gao , which is also combinatorial but whose analysis is based on an relaxation, showing an integrality gap of $`\tfrac{3}{2}+\eps`$ for the . As far as -based algorithms are concerned, the state of the art was a $`\tfrac{7}{4}`$-approximation due to Kortsarz and Nutov until Adjiashvili gave a $`\tfrac{5}{3}+\eps`$-approximation based on the bundle , for any small $`\eps > 0`$.
For special classes of trees, there are better approximation guarantees known. For example, Cohen and Nutov described a $`1+\ln 2 \approx 1.6931`$-approximation for for trees with constant radius. For the special case of where every link ends in two leaves of the tree, Maduel and Nutov gave a $`\tfrac{17}{12}`$-approximation. In case that the tree has radius 3 or 2, they could strengthen the bound to $`\tfrac{11}{8}`$ and $`\tfrac{4}{3}`$, respectively.
In terms of lower bounds on integrality gaps, Cheriyan conjectured that the cut has an integrality gap of $`\tfrac{4}{3}`$ for . However, this conjecture has been refuted by Cheriyan , who proved that the integrality gap of the cut is at least $`\tfrac{3}{2}`$, even for .
Proof of Lemma <a href="#lemma:david" data-reference-type=“ref”
data-reference=“lemma:david”>[lemma:david]
Proof. Consider the following rounding procedure. Given $`(H,x)`$, let $`r`$ be a $`\tfrac{10M}{\eps^2}`$-center of $`(H,x)`$. Let $`H_1,\dots,H_m`$ be the sub-trees of $`H`$ created by removing $`r`$. For a sub-tree $`H_i`$, let $`e_i`$ be the edge that connects $`H_i`$ to $`r`$ in $`H`$. We define new sub-trees $`\overline{H}_i, i \in [m]`$ by adding $`e_i`$ and $`r`$ to $`H_i`$. We now create a new solution $`y`$ that contains no cross-links by splitting all cross-links at $`r`$ into two up-link shadows by defining:
\begin{equation*}
y_{\ell=uv} := \begin{cases} x_\ell & \ell \in \supp(\xin), r \not\in \ell\\ x_\ell + \sum_{\substack{\ell' \in \supp(\xcr)\\v \in \ell'}} x_{\ell'} & u=r,v\in V[H]\setminus\set{r} \\ 0 & \text{else} \end{cases} \quad \text{ for all } \ell \in L.
\end{equation*}
Notice that we have $`c^\tp y \leq c^\tp\xin + 2c^\tp \xcr`$. Since $`y`$ contains no cross-links, it is the union of $`m`$ disjoint solutions, one for each $`\overline{H}_i, i \in [m]`$. We will now focus on rounding the solution $`z`$ for a tree $`\overline{H} \in \setc{\overline{H}_i}{ i \in [m]}`$ with $`z`$ being defined by
\begin{equation*}
z_{\ell} := \begin{cases} y_\ell & \ell \subseteq V[\overline{H}]\\ 0 & \text{else} \end{cases} \quad \text{ for all } \ell \in L.
\end{equation*}
We now make the following claim.
If $`z`$ is a feasible solution to the bundle
\begin{equation*}
\OPT(E[\overline{H}]) \leq c^\tp z + |V^h \cap V[\overline{H}]|.
\end{equation*}
Proof. Let $`e`$ be the edge connecting $`\overline{H}`$ to $`r`$ in $`H`$. Since $`(H,x)`$ is $`\frac{10M}{\eps^2}`$-simple and $`e \not\in E^h`$, we know that $`c^\tp z \leq \tfrac{10M}{\eps^2} + x(\cov(e)) M \leq \frac{12M}{\eps^2}`$. Due to Lemma [lemma:decomposition], property 3 we even have that $`c^\tp z \leq \tfrac{4M}{\eps^2} + x(\cov(e)) M \leq \frac{6M}{\eps^2}`$.
Since $`z`$ fulfills the cut constraints, we know that $`\OPT(E[\overline{H}]) - c^\tp z \leq c^\tp z \leq \frac{6M}{\eps^2}`$. Thus, if $`|V^h \cap V[\overline{H}]| \geq \frac{6M}{\eps^2}`$, we are done. We therefore assume that $`|V^h \cap V[\overline{H}]| < \frac{6M}{\eps^2}`$. In this case, we want to employ the bundle constraints that $`z`$ fulfills, meaning that we have to analyze whether $`\overline{H}`$ is a $`\ceil{\tfrac{28M}{\eps^2}}`$-bundle in $`G`$.
Due to $`(H,x)`$ being $`\frac{10M}{\eps^2}`$-simple, we know that the number of leaves of $`\overline{H}`$ is bounded by $`\tfrac{10M}{\eps^2}+1`$. We now split $`\overline{H}`$ at all nodes with degree at least 3. The number of nodes with degree at least 3 is bounded by the number of leaves. As a consequence, we obtain a decomposition into $`Q_1, \dots, Q_t`$ paths with $`t \leq 2 \cdot \left(\tfrac{10M}{\eps^2}+1\right)`$.
Each $`Q_i`$, $`i\in [t]`$ is the union of paths in $`G`$, separated by components that were contracted into a node of $`V^h`$. If a $`Q_i`$ is the union of $`k`$ paths, there need to be $`k-1`$ nodes in the interior of the path that by construction cannot be part of other $`Q_i`$, since $`Q_i`$ has no degree 3 or higher nodes in its interior. Thus, each node of $`V^h`$ is only capable of splitting a single path in $`\overline{H}`$ into two paths in $`G`$. Thus, we have at most $`2 \cdot \left(\tfrac{10M}{\eps^2}+1\right)`$ paths $`Q_1, \dots, Q_t`$ in $`\overline{G}`$ that are splitted into two paths in $`G`$ at most $`\tfrac{6M}{\eps^2}`$ times. Therefore, $`\overline{H}`$ is the union of at most $`2 \cdot \left(\tfrac{10M}{\eps^2}+1\right) + \tfrac{6M}{\eps^2} \leq \tfrac{28M}{\eps^2}`$ paths in $`G`$ and is therefore a $`\ceil{\tfrac{28M}{\eps^2}}`$-bundle in $`G`$.
Thus, we have by the corresponding bundle-constraint:
c^\tp z = \sum_{\ell \in \cov(\overline{H})} c_\ell z_\ell \geq \sum_{\ell \in \cov(\overline{H})} c_\ell x_\ell \geq \OPT(E[\overline{H}])
which completes the proof of this claim. ◻
Claim [cl:compound] completes the proof, except for the run-time. However, we can use Lemma [lemma:constantleaves] to compute solutions for $`\gamma`$-bundles, and since $`\OPT(E[\overline{H}]) \in O(M)`$, we can use simple enumeration in the case of many contracted nodes. ◻
-
We can assume without loss of generality that the instance is shadow-complete – otherwise, we can add the missing shadows with a cost equal to the cheapest original link of which they are a shadow. Should these links appear in any solution, we can replace them with the original link (which covers even more edges) at no cost. ↩︎
-
Taking $`T`$ to be the set of odd-degree nodes of the tree $`G`$. ↩︎