---
abstract: |
    We study the [*$k$-route*]{} generalizations of various cut problems, the most general of which is *$k$-route multicut* ([$k$-[$\mathsf{MC}$]{}]{}) problem, wherein we have $r$ source-sink pairs and the goal is to delete a minimum-cost set of edges to reduce the edge-connectivity of every source-sink pair to below $k$. The $k$-route extensions of multiway cut ([$k$-[$\mathsf{MWC}$]{}]{}), and the minimum $s$-$t$ cut problem ([$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}), are similarly defined. We present various approximation and hardness results for [$k$-[$\mathsf{MC}$]{}]{}, [$k$-[$\mathsf{MWC}$]{}]{}, and [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}that improve the state-of-the-art for these problems in several cases. Our contributions are threefold.

    [$\bullet$]{}[=0ex ]{}

    For [*$k$-route multiway cut*]{}, we devise simple, but surprisingly effective, combinatorial algorithms that yield bicriteria approximation guarantees that markedly improve upon the previous-best guarantees.

    For [*$k$-route multicut*]{}, we design algorithms that improve upon the previous-best approximation factors by roughly an $O(\sqrt{\log r})$-factor, when $k=2$, and for general $k$ and unit costs and any fixed violation of the connectivity threshold $k$. The main technical innovation is the definition of a new, powerful *region growing* lemma that allows us to perform region-growing in a recursive fashion even though the LP solution yields a [*different metric*]{} for each source-sink pair, and [*without incurring an $O(\log^2 r)$ blow-up*]{} in the cost that is inherent in some previous applications of region growing to $k$-route cuts. We obtain the same benefits as [@EvenNRS00] do in their divide-and-conquer algorithms, and thereby obtain an $O(\ln r\ln\ln r)$-approximation to the cost. We also obtain some extensions to $k$-route node-multicut problems.

    We complement these results by showing that the [*$k$-route $s$-$t$ cut*]{} problem is at least as hard to approximate as the [*densest-$k$-subgraph*]{} ([$\mathsf{D}k\mathsf{S}$]{}) problem on uniform hypergraphs. In particular, this implies that one cannot avoid a ${\operatorname{poly}}(k)$-factor if one seeks a unicriterion approximation, without improving the state-of-the-art for [$\mathsf{D}k\mathsf{S}$]{}on graphs, and proving the existence of a family of one-way functions. Previously, only [[[*NP*]{}]{}-hard]{}ness of [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}was known.
author:
- 'Guru Guruganesh[^1]'
- 'Laura Sanita[^2]'
- 'Chaitanya Swamy[^3]'
title: 'Improved Region-Growing and Combinatorial Algorithms for $k$-Route Cut Problems'
---

Introduction
============

The problem of finding minimum size cuts for a given graph has a rich history in the field of combinatorial optimization, with a wide range of applications in logistics, transportation and telecommunication systems. One key problem of interest is that of disconnecting a given set of node pairs in a network by removing edges at minimum cost. Formally, in the *multicut* problem, we are given an undirected graph $G=(V,E)$ with nonnegative edge costs $\{c_e\}_{e\in E}$ and pairs of nodes $(s_1,t_1),\ldots,(s_r,t_r)$ called source-sink pairs or commodities, and we seek a minimum-cost set of edges whose removal disconnects every $s_i$-$t_i$ pair. Two special cases of this problem have by themselves attracted widespread attention: (i) the celebrated *minimum $s$-$t$ cut* problem, which is the special case when $r=1$; and (ii) the *multiway cut* problem [@DahlhausJPSY92], where every pair of nodes from a given set $T{\subseteq}V$ of terminals forms a commodity. These cut problems and their variants have been widely studied in terms of hardness and approximation (see, e.g., [@Vazirani01; @WilliamsonS10]), have numerous direct applications (e.g., identifying bottlenecks in a network), and algorithms for them serve as important primitives in the design of divide-and-conquer algorithms (see, e.g., [@EvenNRS00; @Shmoys97]) and find application in diverse settings such as image segmentation, VLSI design and network routing (see, e.g., [@LeightonR99; @MalikS00; @BhattL84; @Racke02]).

We study a natural generalization of the above cut problems motivated by the fact that in various settings, we are not interested in a complete disconnection of our terminals but rather in reducing their connectivity below a certain threshold. Specifically, in the [*$k$-route multicut*]{} ([$k$-[$\mathsf{MC}$]{}]{}) problem, the input is a multicut instance and an integer $k\geq 1$; the goal is to find a minimum-cost set $F{\subseteq}E$ of edges so that there are at most $(k-1)$ edge-disjoint $s_i$-$t_i$ paths in $(V,E{\ensuremath{\setminus}}F)$ for all $i=1,\ldots,r$. We define the [*$k$-route multiway cut*]{} ([$k$-[$\mathsf{MWC}$]{}]{}), and the [*$k$-route $(s,t)$-cut*]{} ([$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}) problems analogously.

The study of $k$-route cut problems can be motivated from various perspectives. One motivation comes from the fact that $k$-route cuts are dual objects to [*$k$-route flows*]{} [@Kishimoto96], which can be seen as a robust or fault-tolerant version of flows where we seek to send traffic along tuples of $k$ edge-disjoint paths. A $k$-route cut establishes an upper bound on the value (suitably defined) of the maximum $k$-route flow, and can thus be seen as identifying the bottleneck in a network when we seek a certain level of robustness. $k$-route cut problems can also be directly motivated as abstracting the problem of an attacker who seeks to reduce connectivity in a given network while incurring minimum cost. Viewed from this perspective, $k$-route cut problems are closely related to [*network interdiction*]{} problems, which typically consider the complementary objective of minimizing source-sink connectivity subject to a budget constraint on the edge-removal cost [@Phillips93; @Wood93; @Zenklusen10].

The $k$-route cut problems are at least as hard as their 1-route counterparts. Multicut and multiway cut are [[*APX*]{}]{}-hard [@DahlhausJPSY92], with the former not admitting any constant-factor approximation assuming the unique-games conjecture [@ChawlaGR08], and [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}is [[[*NP*]{}]{}-hard]{}; hence, we focus on approximation algorithms. Moreover, as highlighted in [@ChekuriK08; @BarmanC10; @KolmanS11; @ChuzhoyMVZ12; @KolmanS12], $k$-route cut problems turn out to be much more challenging than their 1-route counterparts, especially for non-constant $k$, so (as in [@ChuzhoyMVZ12]) we consider bicriteria approximation guarantees. (This is further justified by our hardness result for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}in Section \[kstcut\].) We say that a solution $F{\subseteq}E$ is an $(\alpha,\beta)$-approximation for the given [$k$-[$\mathsf{MC}$]{}]{}instance if $\sum_{e\in F}c_e$ is at most $\beta$ times the optimal value, and $(V,E\setminus F)$ contains at most $\alpha(k -1)$ edge-disjoint $s_i$-$t_i$ paths for all $i=1,\ldots,r$.

#### Our results.

We develop various approximation and hardness results for [$k$-[$\mathsf{MC}$]{}]{}, [$k$-[$\mathsf{MWC}$]{}]{} and [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}that improve upon the current-best approximation and hardness results in several cases.

In Section \[kmwc\], we consider the $k$-route multiway cut problem. We devise an $\bigl(O(1),O(1)\bigr)$-approximation for [$k$-[$\mathsf{MWC}$]{}]{}with unit costs (Theorem \[unitthm\]), and an $\bigl(O(1),O(\log r)\bigr)$-approximation with general costs (Theorem \[genthm\]), where $r=|T|$. The previous-best guarantees for [$k$-[$\mathsf{MWC}$]{}]{}(for general $k$) are those that follow from the results of Chuzhoy et al. [@ChuzhoyMVZ12] for [$k$-[$\mathsf{MC}$]{}]{}, namely, an $\bigl(O(1),O(\log^{1.5}r)\bigr)$-approximation for unit costs and an $\bigl(O(\log r),O(\log^3 r)\bigr)$-approximation for general costs. Thus, our guarantees constitute a significant improvement in the state-of-the-art for [$k$-[$\mathsf{MWC}$]{}]{}. We also show that the special case where $T=V$, which we call [*$k$-route all-pairs cut*]{}, is [[*APX*]{}]{}-hard for $k\geq 3$ (Appendix \[append-allpairs\]). (For $k=1, 2$, it is easy to see that all-pairs $k$-route cut is polytime solvable.)

In Section \[kmc\], we design algorithms for the $k$-route multicut problem. We achieve approximation ratios of $O(\ln r\ln\ln r)$ for 2-[$\mathsf{MC}$]{}, and $\bigl({\ensuremath{\gamma}},O(\frac{{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\ln r\ln\ln r)\bigr)$ for [$k$-[$\mathsf{MC}$]{}]{}with unit costs. In contrast, Chuzhoy et al. [@ChuzhoyMVZ12] obtain approximation ratios of $O(\log^{1.5}r)$ for 2-[$\mathsf{MC}$]{}, and $\bigl({\ensuremath{\gamma}},O(\frac{\log^{1.5}r}{\min\{1,{\ensuremath{\gamma}}-1\}})\bigr)$ for [$k$-[$\mathsf{MC}$]{}]{}with unit costs. Thus, for any fixed ${\ensuremath{\gamma}}$ (i.e., independent of $k$ and $r$), our results improve upon the previous-best guarantees for these cases in [@ChuzhoyMVZ12] by roughly an $O(\sqrt{\log r})$-factor. (Setting ${\ensuremath{\gamma}}=\frac{k}{k-1}$, our guarantee and the one in [@ChuzhoyMVZ12] become unicriterion approximations that are incomparable.) In contrast to the algorithms in [@ChuzhoyMVZ12], which rely on approximations to suitable variants of sparsest cut, we devise [*rounding*]{} algorithms for a natural LP-relaxation for [$k$-[$\mathsf{MC}$]{}]{}, and our guarantees therefore also translate to integrality-gap results. In Section \[node-extn\], we consider some extensions to $k$-route [*node-multicut*]{} problems.

Complementing the above results, we show in Section \[kstcut\] that [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}is at least as hard as the densest $k$-subgraph ([$\mathsf{D}k\mathsf{S}$]{}) problem: a $\rho$-approximation for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}yields a $(2\rho^{\ensuremath{\lambda}})$-approximation for [$\mathsf{D}k\mathsf{S}$]{}on ${\ensuremath{\lambda}}$-uniform hypergraphs (Theorem \[kmichard\]). The latter problem is hard to approximate within an $n^{{\ensuremath{\epsilon}}_0}$-factor, for some constant ${\ensuremath{\epsilon}}_0$, for all ${\ensuremath{\lambda}}\geq 3$, unless a certain family of one-way functions exists [@Applebaum11]. This implies that obtaining a unicriterion $O\bigl(k^{{\ensuremath{\epsilon}}_0}{\operatorname{polylog}}(n)\bigr)$-approximation (even) for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}for some constant ${\ensuremath{\epsilon}}_0$ would improve the state-of-the-art for the notoriously hard densest $k$-subgraph problem on graphs, and imply the existence of certain one-way functions. Previously, only [[[*NP*]{}]{}-hard]{}ness of [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}was known, as a consequence of the fact that certain [[[*NP*]{}]{}-hard]{}unbalanced graph partitioning problems [@HayrapetyanKPS05; @LiZ13] can be cast as special cases of [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}.

#### Our techniques.

Our algorithms for [$k$-[$\mathsf{MWC}$]{}]{}are combinatorial, and rely on the following simple, but quite useful observation: if $F{\subseteq}E$ is feasible, then ${\ensuremath{{\ensuremath{\overline{G}}}}}=(V,E{\ensuremath{\setminus}}F)$ has a multiway cut with at most $(k-1)(r-1)$ edges (Claim \[kmwc-prop\]). Using this, we show that we can identify a terminal $t_i\in T$, a $t_i$-isolating cut, and a set of edges of cost $O\bigl(\frac{{\ensuremath{\mathit{OPT}}}}{|T|}\bigr)$ whose removal causes the $t_i$-isolating cut to have $O(k-1)$ edges. We include these edges, drop $t_i$ from $T$, and repeat, which naturally yields an $O(\log r)$-approximation in the cost. The improvement for unit costs stems from the stronger property that either the minimum multiway-cut in $G$ has cost $O({\ensuremath{\mathit{OPT}}})$, or there is some $t_i$-isolating cut of value $O(k-1)$; thus, we may now drop terminals incurring [*zero cost*]{}, which results in an improved $O(1)$ cost-approximation.

Interestingly, [@BarmanC10] use a similar approach to obtain an $\bigl(O(1),O(1)\bigr)$-approximation for single-source [$k$-[$\mathsf{MC}$]{}]{}with unit costs and they remark that such an approach is unlikely to work for [$k$-[$\mathsf{MWC}$]{}]{}because there are examples where every pair of terminals is $2(k-1)$-edge connected but the optimal multiway cut value is $\Omega(r)\cdot{\ensuremath{\mathit{OPT}}}$. Thus, a useful insight to emerge from our work is that whereas a $2$-factor violation in the pairwise terminal connectivity does not ensure that the multiway cut value is $O({\ensuremath{\mathit{OPT}}})$, a $(2+{\ensuremath{\varepsilon}})$-factor violation in connectivity does, for any ${\ensuremath{\varepsilon}}>0$.

Our algorithms for [$k$-[$\mathsf{MC}$]{}]{}are based on rounding an optimal solution to a natural LP-relaxation of the problem. This is technically the most sophisticated part of the paper. The main technique that we use is [*region growing*]{}. The idea is to view the LP solution as a metric, grow a suitable ball in this metric and prove a region-growing lemma showing that the cost of the ball-boundary edges can be charged to the ball-volume, where volume measures the contribution to the LP objective from the edges inside the ball. This was introduced by [@LeightonR99; @GargVY96] in the context of the sparsest cut and multicut problems, and Even et al. [@EvenNRS00], building upon the work of Seymour [@Seymour95], extended the technique to obtain improved guarantees for various divide-and-conquer algorithms that involve recursive applications of region growing. However, in contrast with various applications of region growing considered in [@LeightonR99; @GargVY96; @EvenNSS98; @EvenNRS00], the difficulty in the $k$-route multicut problem stems from the fact that an LP-solution yields a [*different metric for each source-sink pair*]{} instead of a single common metric that can be applied in the region-growing process. (In particular, [$k$-[$\mathsf{MC}$]{}]{}does not fall into the divide-and-conquer framework of Even et al. [@EvenNRS00].) Although [@BarmanC10; @KolmanS11; @KolmanS12] adapted the region-growing lemma in [@LeightonR99; @GargVY96] to the 2-route, 3-route, and the $k$-route single-source settings, their approach seems incapable of obtaining any thing better than an $O(\log^2 r)$-approximation—one loses one $\log$-factor due to region growing and another due to recursion—which is worse than the guarantees in [@ChuzhoyMVZ12]. (In fact, [@ChuzhoyMVZ12] abandoned the region-growing approach and used a greedy set-cover strategy to obtain their improvements over [@BarmanC10; @KolmanS11; @KolmanS12].)

Our chief technical novelty is to prove a region-growing lemma (see Lemmas \[reggrow\] and \[cor2\]) applicable to settings with different metrics, that is inspired by, but more general, than the analogous lemma in [@EvenNRS00], and much more sophisticated than the one used in [@BarmanC10; @KolmanS11; @KolmanS12]. This lemma, coupled with a subtle insight about the metrics derived from the LP solution, allows us to obtain the same kind of savings in our recursive region-growing algorithm that Even et al. [@EvenNRS00] obtain (via their region-growing lemma) in their divide-and-conquer algorithms; this yields our improved approximation guarantees. We believe that our region-growing lemma and its application in the context of different metrics are tools of independent interest that will find further application in the study of cut problems.

The hardness proof for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}dovetails the hardness proof in [@ChuzhoyMVZ12] for the [*vertex-connectivity*]{} version of [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{} (where we want to decrease the $s$-$t$ [*vertex connectivity*]{} to below $k$), who reduce from the [*small-set vertex expansion*]{} ([$\mathsf{SSVE}$]{}) problem, which they show is [$\mathsf{D}k\mathsf{S}$]{}-hard. We observe that this reduction immediately implies the same hardness for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}on a [*directed graph*]{}, and combine this with a useful trick from [@ChakrabartyKLN13] that allows us to move from digraphs to undirected graphs. The idea is to take the digraph used in the hardness proof, remove edge directions, and add some extra nodes and expensive edges so that the *residual digraph* obtained after sending a partial $s$-$t$ flow along the expensive edges essentially coincides with the digraph used in the hardness proof.

#### Related work.

Standard (i.e., 1-route) cut problems have been extensively studied; we refer the reader to the textbooks [@AhujaMO93; @Vazirani01; @WilliamsonS10] for more information.

The study of $k$-route flow and $k$-route cut problems was initiated by Kishimoto [@Kishimoto96], and has since received much attention in the theoretical Computer Science community [@BruhnCHKS08; @ChekuriK08; @BarmanC10; @KolmanS11; @KolmanS12; @ChuzhoyMVZ12]. Bruhn et al. [@BruhnCHKS08] gave a $2(k-1)$-approximation for single-source [$k$-[$\mathsf{MC}$]{}]{} with unit costs, whereas [@ChekuriK08; @BarmanC10; @KolmanS11] obtained efficient polylogarithmic approximation results for [$k$-[$\mathsf{MC}$]{}]{}with small values of $k$. Subsequently, Chuzhoy et al. [@ChuzhoyMVZ12] obtained the first non-trivial results for [$k$-[$\mathsf{MC}$]{}]{}with arbitrary $k$ in the form of bicriteria approximation guarantees. Independently, Kolman and Scheideler [@KolmanS12] obtained an $O\bigl(\exp(k){\operatorname{polylog}}(r)\bigr)$-approximation for single-source [$k$-[$\mathsf{MC}$]{}]{}(with general costs). As shown by our hardness result for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}in Section \[kstcut\], the move to bicriteria approximations is necessary unless one incurs a ${\operatorname{poly}}(k)$-factor in the approximation.

As noted earlier, $k$-route cut problems and [*network interdiction*]{} problems (see, e.g., [@Phillips93; @Wood93; @Zenklusen10; @DinitzG13] and the references therein) can be viewed as complementary problems. For instance, in the [*maximum-flow interdiction problem*]{} ([$\mathsf{MFIP}$]{}) we are given [*edge capacities*]{} in addition to edge costs, and the goal is to minimize the maximum $s$-$t$ flow in the graph remaining after removing edges of total cost at most a given budget. [$\mathsf{MFIP}$]{}with unit capacities is thus complementary to [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}, and bicriteria guarantees for one translate to the other. Unit-capacity [$\mathsf{MFIP}$]{}is known to be polytime solvable for planar graphs [@Phillips93; @Zenklusen10]. Dinitz and Gupta [@DinitzG13] propose a general framework for attacking [*packing interdiction*]{} problems. However, their results do not quite apply to [$\mathsf{MFIP}$]{}(since phrasing max-flow in terms of edge-flows destroys the packing property, and phrasing it in terms of path-flows yields an interdiction problem where one removes paths).

A simple combinatorial algorithm for $k$-route multiway cut {#kmwc}
===========================================================

Recall that in the $k$-route multiway cut ([$k$-[$\mathsf{MWC}$]{}]{}) problem, we are given a set $T=\{t_1,\ldots,t_r\}{\subseteq}V$ of terminals and we seek to remove a minimum-cost set of edges so that the edge-connectivity between any two terminals is less than $k$. The case $k=1$ is the multiway cut problem, which is known to be [[*APX*]{}]{}-hard [@DahlhausJPSY92] even with unit edge costs. We devise an $\bigl(O(1),O(1)\bigr)$-approximation for [$k$-[$\mathsf{MWC}$]{}]{}with unit costs, and an $\bigl(O(1),O(\log r)\bigr)$-approximation with general costs. These improve upon the previous-best guarantees (for general $k$) of $\bigl(O(1),O(\log^{1.5}r)\bigr)$ for unit costs, and $O\bigl(O(\log r),O(\log^3 r)\bigr)$ for general costs due to [@ChuzhoyMVZ12]. Remark \[intgap\] shows that our guarantees also translate to integrality-gap bounds for a suitable LP-relaxation.

Let ${\ensuremath{O^*}}$ denote the optimal set of edges, and let $\overline G = (V, E{\ensuremath{\setminus}}{\ensuremath{O^*}})$ be the remainder graph. Let $k'=k-1$. Our algorithms are quite easy to describe and analyze. We first prove a simple claim about ${\ensuremath{{\ensuremath{\overline{G}}}}}$.

\[kmwc-prop\] There is a set $\overline E$ of edges of ${\ensuremath{{\ensuremath{\overline{G}}}}}$ with ${\ensuremath{\left|{\ensuremath{{\ensuremath{\overline{E}}}}}\right|}}\leq k'(r-1)$ such that ${\ensuremath{O^*}}{\cup}\overline{E}$ is a multiway cut in $G$.

Compute a minimum $t_1$-$t_2$ cut $F$ in ${\ensuremath{{\ensuremath{\overline{G}}}}}$, where $F$ is a set of edges. By the definition of $\overline G$, $|F| \leq k'$. Removing $F$ from $\overline G$ creates at least two components. We can now recurse in each connected component, and after computing at most $r-1$ min cuts, each terminal will be in a different connected component.

The idea behind the algorithm for unit costs is the following. Claim \[kmwc-prop\] shows that the optimal multiway cut in $G$ would be a good approximation to [$k$-[$\mathsf{MWC}$]{}]{}if $|\overline{E}|=O\bigl(|{\ensuremath{O^*}}|\bigr)$. Otherwise, there is a multiway cut in $G$ of cost $O(k'r)$, and so there is some terminal (in fact $\Omega(r)$ terminals) that has an isolating cut in $G$ of size $O(k')$; we simply remove this terminal from $T$ and repeat this process.

\[unitthm\] There is a $\bigl({\ensuremath{\gamma}}, \frac{2{\ensuremath{\gamma}}}{{\ensuremath{\gamma}}-2}\bigr)$-approximation algorithm for [$k$-[$\mathsf{MWC}$]{}]{}with unit costs for any ${\ensuremath{\gamma}}>2$.

For all $i$, compute a minimum $t_i$-isolating cut $F_i$. It is well known that, even with non-unit costs, $\sum_{i=1}^r c(F_i)$ is a 2-approximation to the minimum multiway cut [@DahlhausJPSY92]. In particular, $C=\sum_{i=1}^r|F_i|\leq 2|{\ensuremath{O^*}}|+2|{\ensuremath{{\ensuremath{\overline{E}}}}}|$ (by Claim \[kmwc-prop\]). If $C\geq{\ensuremath{\gamma}}k'r$, then we have $|{\ensuremath{O^*}}|\geq\frac{C}{2}-k'r\geq C\bigl(\frac{1}{2}-\frac{1}{{\ensuremath{\gamma}}}\bigr)$, so taking the union of the $F_i$s yields a $\frac{2{\ensuremath{\gamma}}}{{\ensuremath{\gamma}}-2}$-approximation. Otherwise, there is some $t_i$ such that $|F_i|<{\ensuremath{\gamma}}k'$, so we can simply remove $t_i$ from $T$ and decrease $r$, and repeat.

We remark that the number of iterations can be reduced to $\log_2 r$ at the expense of increasing the connectivity to $2{\ensuremath{\gamma}}k'$, since must be at least $r/2$ terminals such that $|F_i|\leq 2{\ensuremath{\gamma}}k'$.

The condition ${\ensuremath{\gamma}}>2$ above is tight. This follows from an example in [@BarmanC10] where every pair of terminals is $2k'$-edge connected but the minimum multiway cut yields an $\Omega(r)$-approximation.

To generalize this algorithm to general edge costs, assume for now that we know ${\ensuremath{\mathit{OPT}}}=c({\ensuremath{O^*}})$. Unlike in the unit edge-cost case where we could make progress by dropping terminals while incurring zero cost, here we will need to incur cost $O\bigl(\frac{{\ensuremath{\mathit{OPT}}}}{r}\bigr)$ to drop a terminal (or incur cost $O({\ensuremath{\mathit{OPT}}})$ to drop $r/2$ terminals). This naturally leads to an $O(\log r)$-approximation in the cost. Let $H_r:=1+\frac{1}{2}+\ldots+\frac{1}{r}=O(\log r)$.

\[genthm\] There is a $\bigl({\ensuremath{\gamma}},\frac{2{\ensuremath{\gamma}}}{{\ensuremath{\gamma}}-2}H_r\bigr)$-approximation algorithm for [$k$-[$\mathsf{MWC}$]{}]{}with general edge costs, for any ${\ensuremath{\gamma}}>2$.

Let $T'$ initialized to $T$ denote the current terminal set and $r'{\ensuremath{\leftarrow}}r$. Let $F$ initialized to ${\ensuremath{\emptyset}}$ denote the set of edges removed. Let ${\ensuremath{\alpha}}=\frac{2}{{\ensuremath{\gamma}}-2}$. While $|T'|>1$, we do the following. Set $c'_e=\min\{c_e,\frac{{\ensuremath{\alpha}}{\ensuremath{\mathit{OPT}}}}{k'r'}\}$ for every edge $e$. Note that the $c'$-cost of the minimum multiway cut is at most $c'({\ensuremath{O^*}}{\cup}\overline{E})\leq{\ensuremath{\mathit{OPT}}}+k'r'\cdot\frac{{\ensuremath{\alpha}}{\ensuremath{\mathit{OPT}}}}{k'r'}$. For every terminal $t\in T'$, compute a minimum $c'$-cost $t$-isolating cut $F_t$. Then, we have $\sum_{t\in T'}c'(F_t)\leq 2(1+{\ensuremath{\alpha}}){\ensuremath{\mathit{OPT}}}$. So there is some $t\in T'$ such that $c'(F_t)\leq\frac{2(1+{\ensuremath{\alpha}}){\ensuremath{\mathit{OPT}}}}{r'}$. The number of edges in $F_t$ with $c_e>\frac{{\ensuremath{\alpha}}{\ensuremath{\mathit{OPT}}}}{k'r'}$ is less than $\frac{2(1+{\ensuremath{\alpha}})}{{\ensuremath{\alpha}}}k'={\ensuremath{\gamma}}k'$. We add edges in $F_t$ with $c_e\leq\frac{{\ensuremath{\alpha}}{\ensuremath{\mathit{OPT}}}}{k'r'}$ to $F$. This incurs cost at most $c'(F_t)\leq\frac{2(1+{\ensuremath{\alpha}}){\ensuremath{\mathit{OPT}}}}{r'}$, and ensures that $t$ is less than ${\ensuremath{\gamma}}k'$ connected to every other terminal in $T'$ in the remaining graph. We now set $T'{\ensuremath{\leftarrow}}T'{\ensuremath{\setminus}}\{t\},\ r'{\ensuremath{\leftarrow}}r'-1$, and repeat the above process.

Clearly, every pair of terminals is at most ${\ensuremath{\gamma}}k'$ connected in $(V,E{\ensuremath{\setminus}}F)$. Also, $c(F)\leq\sum_{r'=r}^1\frac{2(1+{\ensuremath{\alpha}}){\ensuremath{\mathit{OPT}}}}{r'}={\ensuremath{\mathit{OPT}}}\cdot\frac{2{\ensuremath{\gamma}}}{{\ensuremath{\gamma}}-2}\cdot H_{r}$.

Finally, we can eliminate the need for knowing ${\ensuremath{\mathit{OPT}}}$ as follows. Given a guess $C$ of ${\ensuremath{\mathit{OPT}}}$, if at some iteration we have $\sum_{t\in T'}c'(F_t)>2(1+{\ensuremath{\alpha}})C$ then we know that $C<{\ensuremath{\mathit{OPT}}}$; otherwise, we obtain a solution of cost at most $C\cdot\frac{2{\ensuremath{\gamma}}}{{\ensuremath{\gamma}}-2}\cdot H_r$. So we can try powers of $(1+{\ensuremath{\epsilon}})$ to find the smallest $C$ such that the latter case happens; this blows up the approximation in cost by at most a $(1+{\ensuremath{\epsilon}})$-factor.

\[intgap\] The guarantees in Theorems \[unitthm\] and \[genthm\] also translate to integrality-gap bounds for the following LP-relaxation of [$k$-[$\mathsf{MWC}$]{}]{}. Let ${\ensuremath{\mathcal P}}_{ij}$ be the collection of all $t_i$-$t_j$ paths. $$\min \ \ c^Tx \quad \ \ \text{s.t.} \qquad 
\sum_{e\in P}(x_e+y_e)\geq 1 \quad \forall t_i,t_j\in T,\ P\in{\ensuremath{\mathcal P}}_{ij}; \quad \ \ 
\sum_e y_e\leq k'(r-1); \quad \ \ x,y \geq 0. \tag{P'} \label{kmwclp}$$

Claim \[kmwc-prop\] implies that is indeed a valid relaxation of [$k$-[$\mathsf{MWC}$]{}]{}. Let $(x,y)$ be an optimal solution to and ${{\ensuremath{\mathit{OPT}}}_{\text{\ref{kmwclp}}}}$ be its value. Then, for any ${\ensuremath{\lambda}}\geq 0$, for the cost function $c'_e=\min\{c_e,{\ensuremath{\lambda}}\}$, there is a fractional multiway-cut of $c'$-cost at most ${{\ensuremath{\mathit{OPT}}}_{\text{\ref{kmwclp}}}}+{\ensuremath{\lambda}}\sum_ey_e$. Also, if $F_i$ is a minimum $c'$-cost $t_i$-isolating cut then we have $\sum_ic'(F_i)\leq 2\bigl({{\ensuremath{\mathit{OPT}}}_{\text{\ref{kmwclp}}}}+{\ensuremath{\lambda}}\sum_e y_e\bigr)$. (This follows since an optimal solution to the multiway-cut LP is known to be half-integral (see, e.g., [@Vazirani01]); this implies that 2(cost of an optimal solution) is at least $\sum_{i}$ (cost of a minimum $t_i$-isolating cut).) This implies that we can replace $|{\ensuremath{O^*}}|$ and $|{\ensuremath{{\ensuremath{\overline{E}}}}}|$ in the proof of Theorem \[unitthm\] by ${{\ensuremath{\mathit{OPT}}}_{\text{\ref{kmwclp}}}}$ and $\sum_e y_e$, and ${\ensuremath{\mathit{OPT}}}$ in the proof of Theorem \[genthm\] by ${{\ensuremath{\mathit{OPT}}}_{\text{\ref{kmwclp}}}}$, and all the arguments go through.

#### The all-pairs case.

This is the special case of [$k$-[$\mathsf{MWC}$]{}]{}where $T=V$. To our knowledge, this [*$k$-route all-pairs cut*]{} problem has not been explicitly studied before. When $k=1$, the all-pairs problem is trivial as the remainder graph cannot contain any edge. When $k=2$, this problem is still in [[*P*]{}]{}as the remainder graph is a maximum-cost spanning forest. We prove that the problem is [[*APX*]{}]{}-hard for all $k\geq 3$ (see Appendix \[append-allpairs\]), thus resolving the complexity (with respect to polytime solvability) of $k$-route all-pairs cut. The all-pairs problem can also be stated in terms of properties required of the remainder graph. For example, in $3$-route all-pairs cut, we seek a minimum-cost edge set such that the remainder graph does not contain a [*diamond*]{} as a minor. Interestingly, this is equivalent to requiring that the remainder graph be a maximum-weight *cactus*, which is a graph where every edge lies in at most one cycle. As noted above, this problem is [[*APX*]{}]{}-hard. But we observe that this problem admits an $O(1)$-approximation as a consequence of the results of Fiorini et al. [@FioriniJP10]; see Appendix \[append-allpairs\].

A region-growing algorithm for $k$-route multicut {#kmc}
=================================================

We now consider general $k$-route multicut ([$k$-[$\mathsf{MC}$]{}]{}): given source-sink pairs/commodities $(s_1,t_1),\ldots,(s_r,t_r)$, we want to find a minimum-cost set of edges whose removal reduces the $s_i$-$t_i$ edge connectivity to less than $k$ for all $i=1,\ldots,r$. We consider the following LP-relaxation of the problem, which was also considered by Barman and Chawla [@BarmanC10]. Throughout $e$ indexes the edges in $E$, and $i$ indexes the commodities. Let ${\ensuremath{\mathcal P}}_i$ denote the collection of all (simple) $s_i$-$t_i$ paths in $G$. $$\min \ \ \sum_e c_ex_e \quad \ \ \text{s.t.} \quad \ \ 
\sum_{e\in P} (x_e+y^i_e) \geq 1 \quad \forall i, \forall P\in{\ensuremath{\mathcal P}}_i; \quad \ \ 
\sum_e y^i_e \leq k-1 \quad \forall i; \quad \ \ 
x,y \geq 0.
\tag{P} \label{kmcp}$$

Let $\bigl(x,\{y^i\}\bigr)$ denote an optimal solution to , and ${\ensuremath{\mathit{OPT}}}$ be its value. We show that $\bigl(x,\{y^i\}\bigr)$ can be rounded to yield an $O(\ln r\ln\ln r)$-approximation when $k=2$, and a bicriteria $\bigl({\ensuremath{\gamma}},O\bigl(\frac{{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\ln r\ln\ln r\bigr)\bigr)$-approximation for [$k$-[$\mathsf{MC}$]{}]{}with unit edge costs, for any ${\ensuremath{\gamma}}>1$. Notably, our cost-approximation is with respect to ${\ensuremath{\mathit{OPT}}}$, and so they translate to integrality-gap upper bounds for . Our results improve upon (for any fixed ${\ensuremath{\gamma}}$) the previous-best guarantees for these cases in [@ChuzhoyMVZ12] by roughly a $\sqrt{\log r}$-factor.

Region-growing lemmas {#reglemmas}
---------------------

The main tool that we leverage is [*region growing*]{} [@LeightonR99; @GargVY96; @EvenNRS00]. The idea is to view the LP solution as a metric, grow a suitable ball in this metric and prove that the cost of the ball-boundary edges can be charged to the ball-volume, where volume measures the contribution to the LP objective from the edges inside the ball. The main difficulty in applying this idea to is that, unlike most applications of region growing [@LeightonR99; @GargVY96; @EvenNSS98; @EvenNRS00], the LP solution yields a [*different metric*]{} for each commodity. The key technical ingredient and novelty is a new region-growing lemma (see Lemmas \[reggrow\] and \[cor2\]) that is analogous to, but more general, than the one in [@EvenNRS00], and much more sophisticated than the one used in [@BarmanC10; @KolmanS11; @KolmanS12]. Roughly speaking, we prove that given a current set $S$ of nodes, one can construct a ball around any $s_i$ in the $(x+y^i)$-metric such that the cost of the “boundary $x$-edges” can be charged to $(x\text{-volume of $S$})\cdot
\ln\bigl(\text{$x$-volume of $S$}/\text{$x$-volume of the ball}\bigr)\cdot\ln\ln r$ (Lemma \[cor2\]). A subtle insight that helps deal with the complication that different applications of region growing involve different commodities and therefore different metrics is the following. Since the $x$-contribution is [*common*]{} to all $(x+y^i)$-metrics, even though we consider different commodity-metrics we can leverage the above guarantee and obtain the same kind of savings that [@EvenNRS00] obtain in their divide-and-conquer algorithms (see Lemma \[indnlem\]); this leads to our improved approximation guarantees.

We now state our region-growing lemmas in a general form and then apply these to the optimal solution $\bigl(x,\{y^i\}\bigr)$ to obtain various useful corollaries. In Section \[node-extn\], we extend our arguments to prove region-growing lemmas in settings that involve both edge and node lengths. Let $n=|V|$, $m=|E|$. Let $\ell:V\times V\mapsto{\ensuremath{\mathbb R}}_{\geq 0}$ be a metric on $V\times V$. Our algorithm will iteratively focus on certain regions of the graph $G$. Let $S{\subseteq}V$, which is intended to represent the node-set of the current region. Let $F{\subseteq}E$, which is intended to represent the edges that contribute to the volume, and whose cost we care about. Let $\beta\geq 0$. Let $z\in V$, and $\rho\geq 0$.

[$\bullet$]{}[=0ex =0.5ex ]{}

Define $B_\ell(z,\rho):=\{v\in V: \ell_{zv}\leq\rho\}$ to be the ball of radius $\rho$ around $z$.

Let $B^S_\ell(z,\rho):=B_\ell(z,\rho)\cap S$ and ${\ensuremath{\overline{B^S_\ell}}}(z,\rho):=S{\ensuremath{\setminus}}B_\ell(z,\rho)$.

Define the following volumes: $$\begin{aligned}
{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta; z,\rho) & := & \beta
+\sum_{(u,v)\in F: u,v\in B^S_\ell(z,\rho)}c_{uv}\ell_{uv}
+\sum_{\substack{(u,v)\in F: u\in B^S_\ell(z,\rho) \\ v\in{\ensuremath{\overline{B^S_\ell(z,\rho)}}}}}c_{uv}\bigl(\rho-\ell_{zu}\bigr) \\[0.5ex]
{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta; z,\rho) & := & \beta
+\sum_{(u,v)\in F: u,v\in{\ensuremath{\overline{B^S_\ell(z,\rho)}}}}c_{uv}\ell_{uv}
+\sum_{\substack{(u,v)\in F: u\in B^S_\ell(z,\rho) \\ v\in{\ensuremath{\overline{B^S_\ell(z,\rho)}}}}}c_{uv}\bigl(\ell_{zv}-\rho\bigr)\end{aligned}$$

For a subset $T{\subseteq}S$, let ${\ensuremath{\delta}}^S_F(T)$ denote $\{(u,v)\in F: u\in T, v\in S{\ensuremath{\setminus}}T\}$.

Define ${\ensuremath{\partial}}^{S,F}_\ell(z,\rho):={\ensuremath{\delta}}^S_F\bigl(B^S_\ell(z,\rho)\bigr)={\ensuremath{\delta}}^S_F\bigl({\ensuremath{\overline{B^S_\ell}}}(z,\rho)\bigr)$.

When $F=E$, we drop $F$ from the above pieces of notation (e.g., ${\ensuremath{\delta}}^S_E(T)$ is shortened to ${\ensuremath{\delta}}^S(T)$). For $H{\subseteq}E$, we use $\ell(H)$ to denote $\sum_{e\in H}\ell_e$.

\[reggrow\] Let $F{\subseteq}E$, $S{\subseteq}V$, $z\in V$, and $0\leq a<b$. Let $\rho$ be chosen uniformly at random from $[a,b)$. Then, $$\begin{aligned}
{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c\bigl({\ensuremath{\partial}}^{S,F}_\ell(z,\rho)\bigr)}
{{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,\rho)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,b)}{{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,\rho)}\Bigr)}\right]}}
& \leq & \frac{1}{b-a}\cdot\ln\ln\biggl(\frac{e{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,b)}{{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,a)}\biggr) 
\quad\text{and} \label{regbnd1} \\
{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c\bigl({\ensuremath{\partial}}^{S,F}_\ell(z,\rho)\bigr)}
{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,\rho)\ln\Bigl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,\rho)}\Bigr)}\right]}}
& \leq &
\frac{1}{b-a}\cdot\ln\ln\biggl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,b)}\biggr) 
\label{regbnd2}\end{aligned}$$

We abbreviate $c\bigl({\ensuremath{\partial}}^{S,F}_\ell(z,\rho)\bigr)$ to $c(\rho)$, ${\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,\rho)$ to ${\ensuremath{\mathcal V}}(\rho)$ and ${\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,\rho)$ to ${\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)$. Let $x^-$ be a value infinitesimally smaller than $x$. Let $I=\{\ell(s_i,v): v\in V\}$. Note that ${\ensuremath{\mathcal V}}(\rho)$ and ${\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)$ are differentiable at all $\rho\in[a,b){\ensuremath{\setminus}}I$ and for each such $\rho$, we have $\frac{d{\ensuremath{\mathcal V}}(\rho)}{d\rho}=c(\rho)$ and $\frac{d{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)}{d\rho}=-c(\rho)$. Let $a_0=a$, $a_k=b$, and $\{a_1,\ldots,a_{k-1}\}=(a,b)\cap I$ with $a_1<\ldots<a_{k-1}$. Then, $$\begin{aligned}
{2}
(b-a)&\cdot{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c(\rho)}{{\ensuremath{\mathcal V}}(\rho)\ln\bigl(\frac{e{\ensuremath{\mathcal V}}(b)}{{\ensuremath{\mathcal V}}(\rho)}\bigr)}\right]}}
&& =  
\sum_{i=1}^k\int_{a_{i-1}}^{a_i^-}\frac{d{\ensuremath{\mathcal V}}(\rho)}{{\ensuremath{\mathcal V}}(\rho)\ln\bigl(\frac{e{\ensuremath{\mathcal V}}(b)}{{\ensuremath{\mathcal V}}(\rho)}\bigr)}
=
\sum_{i=1}^k\int_{a_{i-1}}^{a_i^-}\frac{d(\ln{\ensuremath{\mathcal V}}(\rho))}{\ln\bigl(e{\ensuremath{\mathcal V}}(b)\bigr)-\ln({\ensuremath{\mathcal V}}(\rho))} \\
& = \sum_{i=1}^k-\ln\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}(b)}{{\ensuremath{\mathcal V}}(\rho)}\Bigr)\biggr|_{a_{i-1}}^{a_i^-} 
&& = \sum_{i=1}^k
\biggl[\ln\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}(b)}{{\ensuremath{\mathcal V}}(a_{i-1})}\Bigr)-\ln\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}(b)}{{\ensuremath{\mathcal V}}(a_i^-)}\Bigr)\biggr]
\leq \ln\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}(b)}{{\ensuremath{\mathcal V}}(a)}\Bigr).\end{aligned}$$ The final inequality follows since ${\ensuremath{\mathcal V}}(\rho)$ increases with $\rho$, and $\ln\ln\bigl(\frac{e{\ensuremath{\mathcal V}}(b)}{{\ensuremath{\mathcal V}}(\rho)}\bigr)$ decreases with $\rho$. Similarly, we obtain that $
(b-a)\cdot{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c(\rho)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)\ln\bigl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)}\bigr)}\right]}}
= \sum_{i=1}^k\ln\ln\Bigl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)}\Bigr)\biggr|_{a_{i-1}}^{a_i^-} 
\leq \ln\ln\Bigl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(b)}\Bigr)$.

\[cor1\] Let $F, H {\subseteq}E$, $S{\subseteq}V$, $z\in V$, and $0\leq a<b$. For any ${\ensuremath{\alpha}}\in(0,1)$, we can efficiently find a radius $\rho_1\in[a,b)$ such that $$\begin{aligned}
c\bigl({\ensuremath{\partial}}^{S,F}_\ell(z,\rho_1)\bigr)
& \leq & \frac{2}{(1-{\ensuremath{\alpha}})(b-a)}
\cdot{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,\rho_1)
\cdot\ln\Bigl(\tfrac{e{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,b)}{{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,\rho_1)}\Bigr)
\cdot\ln\ln\Bigl(\tfrac{e{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,b)}{{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,a)}\Bigr)
\label{fbound1} \\
c\bigl({\ensuremath{\partial}}^{S,F}_\ell(z,\rho_1)\bigr)
& \leq & \frac{2}{(1-{\ensuremath{\alpha}})(b-a)}
\cdot{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,\rho_1)
\cdot\ln\Bigl(\tfrac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,\rho_1)}\Bigr)
\cdot\ln\ln\Bigl(\tfrac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,b)}\Bigr)
\label{fbound2} \\
\bigl|{\ensuremath{\partial}}^{S,H}_\ell(z,\rho_1)\bigr| & < & \frac{\ell(H)}{{\ensuremath{\alpha}}(b-a)}. \label{hbound2}\end{aligned}$$

Suppose we pick $\rho$ uniformly at random from $[a,b)$. Define the following events. $$\begin{aligned}
{\ensuremath{\Omega}}& := & \Bigl\{\rho\in[a,b):
\bigl|{\ensuremath{\partial}}^{S,H}_\ell(z,\rho)\bigr|\geq\frac{\ell(H)}{{\ensuremath{\alpha}}(b-a)}\Bigr\} \\
{\ensuremath{\Omega}}_1 & := & \biggl\{\rho\in[a,b): c\bigl({\ensuremath{\partial}}^{S,F}_\ell(z,\rho)\bigr)>
\frac{2}{(1-{\ensuremath{\alpha}})(b-a)}\cdot 
{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,\rho)\ln\Bigl(\tfrac{e{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,b)}{{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,\rho)}\Bigr)
\cdot\ln\ln\Bigl(\tfrac{e{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,b)}{{\ensuremath{\mathcal V}}^{S,F}_\ell(\beta;z,a)}\Bigr)\biggr\}
\\
{\ensuremath{\Omega}}_2 & := & \biggl\{\rho\in[a,b): c\bigl({\ensuremath{\partial}}^{S,F}_\ell(z,\rho)\bigr)>
\frac{2}{(1-{\ensuremath{\alpha}})(b-a)}\cdot 
{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,\rho)\ln\Bigl(\tfrac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,\rho)}\Bigr)
\cdot\ln\ln\Bigl(\tfrac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_\ell}}}(\beta;z,b)}\Bigr)\biggr\}\end{aligned}$$

For an edge $(u,v)\in E$, $\Pr\bigl[(u,v)\in{\ensuremath{\partial}}^{S,H}_\ell(z,\rho)]\leq\frac{\ell_{uv}}{b-a}$. Hence, ${\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\bigl|{\ensuremath{\partial}}^{S,H}_\ell(z,\rho)\bigr|\right]}}\leq\frac{\ell(H)}{b-a}$, and therefore $\Pr[{\ensuremath{\Omega}}]\leq{\ensuremath{\alpha}}$. By Lemma \[reggrow\] and Markov’s inequality, we have that $\Pr[{\ensuremath{\Omega}}_1],\Pr[{\ensuremath{\Omega}}_2]<(1-{\ensuremath{\alpha}})/2$. Conditioning on ${\ensuremath{\Omega}}^c:=[a,b){\ensuremath{\setminus}}{\ensuremath{\Omega}}$ increases the probability of an event by at most a factor $\frac{1}{1-\Pr[{\ensuremath{\Omega}}]}\leq\frac{1}{1-{\ensuremath{\alpha}}}$, so $\Pr[{\ensuremath{\Omega}}_1|{\ensuremath{\Omega}}^c],\Pr[{\ensuremath{\Omega}}_2|{\ensuremath{\Omega}}^c]<1/2$. Thus, $\Pr[{\ensuremath{\Omega}}^c\cap{\ensuremath{\Omega}}_1^c\cap{\ensuremath{\Omega}}_2^c]>0$.

We argue that ${\ensuremath{\Omega}}^c$, ${\ensuremath{\Omega}}_1^c$, and ${\ensuremath{\Omega}}_2^c$ are all unions of at most $n$ subintervals of $[a,b)$, and we can find these efficiently. Since $\Pr[{\ensuremath{\Omega}}^c\cap{\ensuremath{\Omega}}_1^c\cap{\ensuremath{\Omega}}_2^c]>0$, we can then efficiently find an interval contained in ${\ensuremath{\Omega}}^c\cap{\ensuremath{\Omega}}_1^c\cap{\ensuremath{\Omega}}_2^c$ (in fact, a non-singleton interval), and hence, find $\rho_1\in[a,b]$ satisfying – (in fact, there are infinitely many such $\rho$).

There are at most $n$ distinct sets $B^S_\ell(z,\rho)$ that one may encounter as $\rho$ varies in $[a,b)$. For each such set $A$, there is an interval $[{\ensuremath{\mathit{lb}}},{\ensuremath{\mathit{ub}}})$ such that $A=B^S_\ell(z,\rho)$ for all $\rho\in[{\ensuremath{\mathit{lb}}},{\ensuremath{\mathit{ub}}})$. Note that the right-hand-sides (RHSs) of and are continuous and differentiable in $({\ensuremath{\mathit{lb}}},{\ensuremath{\mathit{ub}}})$, and are monotonic (increasing and decreasing, respectively) functions of $\rho$. We call $[{\ensuremath{\mathit{lb}}},{\ensuremath{\mathit{ub}}})$ a smooth subinterval of $[a,b)$. By definition, the left-hand-sides (LHSs) of – are invariant over a smooth subinterval. Hence, ${\ensuremath{\Omega}}^c$ is the union of some smooth subintervals. Consider a smooth subinterval $[{\ensuremath{\mathit{lb}}},{\ensuremath{\mathit{ub}}})$. By continuity, if some $\rho\in[{\ensuremath{\mathit{lb}}},{\ensuremath{\mathit{ub}}})$ satisfies or , then we can efficiently find the maximal subinterval of $[{\ensuremath{\mathit{lb}}},{\ensuremath{\mathit{ub}}})$ (which may be a singleton interval) such that all $\rho$ in the subinterval satisfy the given bound. Hence, both ${\ensuremath{\Omega}}_1^c$ and ${\ensuremath{\Omega}}_2^c$ are the union of at most $n$ subintervals of $[a,b)$. By trying out the at most $3n$ possible subintervals of ${\ensuremath{\Omega}}_1^c\cup{\ensuremath{\Omega}}_1^c\cup{\ensuremath{\Omega}}_2^c$, we can find some interval contained in ${\ensuremath{\Omega}}^c\cap{\ensuremath{\Omega}}_!^c\cap{\ensuremath{\Omega}}_2^c$ and hence obtain $\rho_1$ satisfying –.

#### Applications of the region-growing lemmas.

To apply the above results to the metrics obtained from $\bigl(x,\{y^i\}\bigr)$, it will be convenient to modify $G$ by subdividing every edge $e$ into $r+1$ edges $e_0,e_1,\ldots,e_r$, and setting $x_{f}=x_e$ for $f=e_0$ and 0 otherwise, and $y^i_f=y^i_e$ if $f=e_i$ and 0 otherwise. We call $e_0$ an $x$-edge, and we call $e_i$, a $y^i$-edge for $i=1,\ldots,r$. Clearly, any solution in $G$ yields a solution in the subdivided graph of the same cost and vice versa, and this holds even for fractional solutions to . In the sequel, we work with the subdivided graph. To keep notation simple, we continue to use $G=(V,E)$ to denote the subdivided graph, and $\bigl(x,\{y^i\}\bigr)$ to denote the above solution in the subdivided graph. Let $F$ be the set of all $x$-edges, and $H^i$ be the set of all $y^i$-edges for all $i=1,\ldots,r$. Consider a commodity $i$. Let $\ell^i$ denote the shortest-path metric of $G$ (i.e., the subdivided graph) induced by the $\{x_e+y^i_e\}$ edge lengths. Set $\beta={\ensuremath{\mathit{OPT}}}/r$. To avoid cumbersome notation, we shorten:

  ---------------------------------------------------------------------------------------------------------------------------------------------------------------------- ----------------------------------------------------------------------------------------------------------------------------------------------
  $B^S_{\ell^i}(z,\rho)$ and ${\ensuremath{\overline{B^S_{\ell^i}}}}(z,\rho)$                                                                                            $B^S_i(z,\rho)$ and ${\ensuremath{\overline{B^S_i}}}(z,\rho)$ respectively
  \[0.75ex\] ${\ensuremath{\mathcal V}}^{S,F}_{\ell^i}(\beta; z,\rho)$ and ${\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,F}_{\ell^i}}}}(\beta; z,\rho)$           ${\ensuremath{\mathcal V}}^{S,x}_i(z,\rho)$ and ${\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho)$ respectively
  \[0.5ex\] ${\ensuremath{\partial}}^{S,F}_{\ell^i}(z,\rho)$, ${\ensuremath{\partial}}^{S,H^i}_{\ell^i}(z,\rho)$, and ${\ensuremath{\partial}}^{S,E}_{\ell^i}(z,\rho)$   ${\ensuremath{\partial}}^{S,x}_i(z,\rho)$, ${\ensuremath{\partial}}^{S,y}_i(z,\rho)$, and ${\ensuremath{\partial}}^S_i(z,\rho)$ respectively
  \[1.25ex\]                                                                                                                                                             
  ---------------------------------------------------------------------------------------------------------------------------------------------------------------------- ----------------------------------------------------------------------------------------------------------------------------------------------

Also, define ${\ensuremath{\mathcal V}}^x(S):=\beta+\sum_{e\in E(S)}c_ex_e$, where $E(S)$ is the set of edges having both endpoints in $S$. Finally, for an integer $q\geq 1$ and a set $A$ of edges let $c^q(A)$ be the cost of all but the $q-1$ most expensive edges of $A$ (so $c^q(A)=0$ if $|A|<q$).

\[cor2\] Let $S{\subseteq}V$, $z\in V$, and $i$ be some commodity. Let ${\ensuremath{\alpha}}\in(0,1)$ and $q={\ensuremath{\left\lceil\frac{k-1}{{\ensuremath{\alpha}}}\right\rceil}}$. We can efficiently find $\rho_1\in[0,1)$ such that $$\begin{aligned}
c^q\bigl({\ensuremath{\partial}}^S_i(z,\rho_1)\bigr) & \leq &
\frac{2}{1-{\ensuremath{\alpha}}}\cdot{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{x}(S)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr) \label{xbound1} \\
c^q\bigl({\ensuremath{\partial}}^S_i(z,\rho_1)\bigr) & \leq &
\frac{2}{1-{\ensuremath{\alpha}}}\cdot{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{x}(S)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr). \label{xbound2}\end{aligned}$$

We apply Corollary \[cor1\] taking $\ell=\ell^i$, $H=H^i$, and $[a,b)=[0,1)$ (and $S$, $z$, ${\ensuremath{\alpha}}$ as given by the statement of the lemma, and $F$ to be the set of $x$-edges). Note that $${\ensuremath{\mathcal V}}^{S,x}_i(z,1)\leq{\ensuremath{\mathcal V}}^x(S), \qquad {\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,0)\leq{\ensuremath{\mathcal V}}^x(S), \qquad
\frac{{\ensuremath{\mathcal V}}^{S,x}_i(z,1)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,0)}\leq r+1, \qquad
\frac{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,0)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,1)}\leq r+1.$$

Thus, we obtain $0\leq\rho_1<1$ such that $c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho_1)\bigr)$ satisfies the bounds given by the RHS of and respectively. Moreover, since $\ell^i(H^i)\leq k-1$, we have that $\bigl|{\ensuremath{\partial}}^{S,y}_i(z,\rho_1)\bigr|<\frac{k-1}{{\ensuremath{\alpha}}}$ due to , and so $\bigl|{\ensuremath{\partial}}^{S,y}_i(z,\rho_1)\bigr|<q$. Finally, since edges not in $F\cup H^i$ have zero $\ell^i$-length, ${\ensuremath{\partial}}^{S,x}_i(z,\rho)$ and ${\ensuremath{\partial}}^{S,y}_i(z,\rho)$ partition ${\ensuremath{\partial}}^S_i(z,\rho)$ for all $\rho$. Therefore, $c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho_1)\bigr)\geq c^q\bigl({\ensuremath{\partial}}^S_i(z,\rho_1)\bigr)$, and the lemma follows.

\[cor3\] Let $S{\subseteq}V$. Suppose that $s_i,t_i\in S$ and there are ${\ensuremath{\gamma}}(k-1)$ edge-disjoint $s_i$-$t_i$ paths internal to $S$, where ${\ensuremath{\gamma}}>1$. Suppose that $c_e=1$ for all edges $e$. We can efficiently find $\rho_1\in[0,1)$ such that $$\begin{aligned}
c\bigl({\ensuremath{\partial}}^S_i(s_i,\rho_1)\bigr)\leq
\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\cdot{\ensuremath{\mathcal V}}^{S,x}_i(s_i,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^x(S)}{{\ensuremath{\mathcal V}}^{S,x}_i(s_i,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr) \\
c\bigl({\ensuremath{\partial}}^S_i(s_i,\rho_1)\bigr)\leq
\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\cdot{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(s_i,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^x(S)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(s_i,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr)\end{aligned}$$

We apply Lemma \[cor2\] with ${\ensuremath{\alpha}}\in(0,1)$, whose value we will fix later, to obtain $\rho_1\in[0,1)$. Note that $B^S_i(s_i,\rho_1)$ is an $s_i$-$t_i$ cut. Let $q={\ensuremath{\left\lceil\frac{k-1}{{\ensuremath{\alpha}}}\right\rceil}}$. For any $s_i$-$t_i$ cut $A{\subseteq}S$, we know that $c\bigl({\ensuremath{\delta}}^S(A)\bigr)\geq{\ensuremath{\gamma}}(k-1)$. Therefore, since we have unit edge costs, $c\bigl({\ensuremath{\delta}}^S(A)\bigr)\leq c^q\bigl({\ensuremath{\delta}}^S(A)\bigr)+q-1
\leq c^q\bigl({\ensuremath{\delta}}^S(A)\bigr)\cdot\frac{{\ensuremath{\gamma}}}{{\ensuremath{\gamma}}-1/{\ensuremath{\alpha}}}$. Plugging in the bounds for $c^q(.)$ from , , we see that the constant factor multiplying the volume terms on the RHS is $\frac{2{\ensuremath{\gamma}}}{(1-{\ensuremath{\alpha}})({\ensuremath{\gamma}}-1/{\ensuremath{\alpha}})}$. This factor is minimized by setting ${\ensuremath{\alpha}}={\ensuremath{\gamma}}^{-1/2}$, which yields the constant factor $\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}$ and completes the proof.

The rounding algorithms and their analyses {#kmcalg}
------------------------------------------

#### The case $k=2$.

The algorithm for $k=2$ follows a similar template as the algorithm in [@BarmanC10] for 2-[$\mathsf{MC}$]{}. However, its analysis resulting in our improved guarantee relies crucially on Lemma \[cor2\] which is derived from our stronger region-growing lemma. The algorithm proceeds as follows. Given a current set $U$ of nodes and a current set of $N$ source-sink pairs, we repeatedly use Lemma \[cor2\] to “carve out” disjoint regions $A_1,\ldots,A_h{\subseteq}U$ and build a set $Z$ of edges until there is no $2$-edge-connected source-sink pair in $U{\ensuremath{\setminus}}(A_1\cup\ldots\cup A_h)$. Given $A_1,\ldots,A_{p-1}$, we obtain $A_p$ as follows. We choose an $s_i$-$t_i$ pair that is at least $2$-edge connected in $S=U{\ensuremath{\setminus}}\bigcup_{q=1}^{p-1}A_q$, and use Lemma \[cor2\] with center $s_i$, ${\ensuremath{\alpha}}=0.5$ and set $S$. Note that $2(k-1)=2=k$. We set $A_p$ to be $B^S_i(s_i,\rho_1)$ or ${\ensuremath{\overline{B^S_i}}}(s_i,\rho_1)$, so as to ensure there are at most $N/2$ source-sink pairs inside $A_p$. We add the edges corresponding to $c^2\bigl({\ensuremath{\delta}}^S(A_p)\bigr)$ to $Z$; Lemma \[cor2\] ensures that the cost of these edges can be bounded in terms of the volume contained in $A_p$. Having obtained $A_1,\ldots,A_h$ this way, we now recurse on each set $A_p$ and the source-sink pairs contained in $A_p$, to obtain edge-sets $Z_1,\ldots,Z_h$. The solution we return is $Z\cup(Z_1\cup\ldots\cup Z_h)$. A more formal description follows.

------------------------------------------------------------------------

[**Algorithm**]{} ${\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}(U, {\ensuremath{\mathcal T}}=\{(s_1,t_1),\ldots,(s_{N},t_{N})\}$)

[**Input:**]{} A subset $U{\subseteq}V$, and a collection ${\ensuremath{\mathcal T}}=\{(s_i,t_i)\}_{i=1}^{N}$ of source-sink pairs, where $s_i,t_i\in U$ for all $i=1,\ldots, N$.

[**Output:**]{} A set $Z{\subseteq}E(U)$ such that $s_i$ and $t_i$ are at most 1-edge-connected in $(U,E(U){\ensuremath{\setminus}}Z)$ for all $i=1,\ldots,N$.

[A.]{}[ =0.5ex =0.25ex ]{}

Set $S=U$, $Z={\ensuremath{\emptyset}}$, ${\ensuremath{\mathcal S}}={\ensuremath{\emptyset}}$, and ${\ensuremath{\mathcal T}}'=\bigl\{(s_i,t_i)\in{\ensuremath{\mathcal T}}: \text{$s_i$ and $t_i$ are at least $2$-edge-connected in $(S,E(S))$}\bigr\}$.

If ${\ensuremath{\mathcal T}}'={\ensuremath{\emptyset}}$, return $Z$.

While ${\ensuremath{\mathcal T}}'\neq{\ensuremath{\emptyset}}$, we do the following.

[A.]{}[=0.25ex =0ex ]{}

Pick some $(s_i,t_i)\in{\ensuremath{\mathcal T}}'$.

Apply Lemma \[cor2\] with $z=s_i$, ${\ensuremath{\alpha}}=0.5$ and the set $S$ to find a radius $0\leq\rho_1<1$ satisfying , .

If $B^S_i(s_i,\rho_1)$ contains at most $N/2$ pairs from ${\ensuremath{\mathcal T}}$ then set $A=B^S_i(s_i,\rho_1)$, else set $A={\ensuremath{\overline{B^S_i}}}(s_i,\rho_1)$.

Set ${\ensuremath{\mathcal S}}{\ensuremath{\leftarrow}}{\ensuremath{\mathcal S}}\cup\{A\}$. Add the edges contributing to $c^2\bigl({\ensuremath{\delta}}^S(A)\bigr)$ (i.e., all edges of ${\ensuremath{\delta}}^S(A)$ except the most-expensive one) to $Z$.

Set $S{\ensuremath{\leftarrow}}S{\ensuremath{\setminus}}A$. Update ${\ensuremath{\mathcal T}}'$ to be the $s_i$-$t_i$ pairs from ${\ensuremath{\mathcal T}}$ that are at least $2$-edge-connected in $(S,E(S))$.

For every set $A\in{\ensuremath{\mathcal S}}$, set $Z{\ensuremath{\leftarrow}}Z\cup{\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}\bigl(A,\{(s_i,t_i)\in{\ensuremath{\mathcal T}}: s_i,t_i\in A\}\bigr)$.

Return $Z$.

The initial call to ${\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}$, which computes the solution we return, is ${\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}\bigl(V,\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$.

------------------------------------------------------------------------

Let $Z:={\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}\bigl(V,\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$. The feasibility of $Z$ follows from the same arguments as in [@BarmanC10]; Lemma \[2feas\] gives a self-contained proof. We focus on showing that $c(Z)\leq O(\ln r\ln\ln r)\cdot{\ensuremath{\mathit{OPT}}}$. Consider the recursion tree associated with the execution of [$2\text{-}$[$\mathsf{MCAlg}$]{}]{}, where each node is labeled with arguments passed in the current invocation of ${\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}$. Define the depth of a subtree of this recursion tree to be the maximum number of edges on a root to leaf path of the subtree. Recall that $\beta=\frac{{\ensuremath{\mathit{OPT}}}}{r}$.

\[indnlem\] Let $d$ be the depth of a subtree of the recursion tree rooted at $\bigl(U{\subseteq}V,{\ensuremath{\mathcal T}}{\subseteq}\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$. Let $Z_U={\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}(U,{\ensuremath{\mathcal T}})$. We have $c(Z_U)\leq 
4\bigl(\beta|{\ensuremath{\mathcal T}}|+{\ensuremath{\mathcal V}}^x(U)\bigr)\ln\bigl(\frac{e^dr{\ensuremath{\mathcal V}}^x(U)}{{\ensuremath{\mathit{OPT}}}}\bigr)\ln\ln\bigl(e(r+1)\bigr)$.

The proof is by induction on $d$. If $d=0$, there is no recursive call; so $Z^U={\ensuremath{\emptyset}}$, which satisfies the stated bound. Otherwise, suppose that we make the recursive calls ${\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}(A_1,{\ensuremath{\mathcal T}}_1),\ldots,{\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}(A_h,{\ensuremath{\mathcal T}}_h)$ in step A4 to obtain edge-sets $Z_1,\ldots,Z_h$ respectively. For $p=1,\ldots,h$, let $S_p$ be the current set $S$ when set $A_p$ was added to ${\ensuremath{\mathcal S}}$ in step A3.4 (so $A_p{\subseteq}S_p$), and let $E_p$ be the edge-set added to $Z$ in this step. Then, $Z_U=\bigcup_{p=1}^h(E_p\cup Z_p)$. By the induction hypothesis, $c(Z_p)\leq 4\Bigl(\beta|{\ensuremath{\mathcal T}}_p|+{\ensuremath{\mathcal V}}^x(A_p)\Bigr)\ln\Bigl(\frac{e^{d-1}r{\ensuremath{\mathcal V}}^x(A_p)}{{\ensuremath{\mathit{OPT}}}}\Bigr)\ln\ln\bigl(e(r+1)\bigr)$. Let ${\ensuremath{\mathrm{vol}}}_p={\ensuremath{\mathcal V}}^{S_p,x}_i(s_i,\rho_1)$ if $A_p=B^{S_p}_i(s_i,\rho_1)$ and ${\ensuremath{\mathrm{vol}}}_p={\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S_p,x}_i}}}(s_i,\rho_1)$ if $A_p={\ensuremath{\overline{B^{S_p}_i}}}(s_i,\rho_1)$. Note that ${\ensuremath{\mathcal V}}^x(A_p)\leq{\ensuremath{\mathrm{vol}}}_p\leq{\ensuremath{\mathcal V}}^x(S_p)\leq{\ensuremath{\mathcal V}}^x(U)$. By Lemma \[cor2\] and the above upper bounds, we have $$\begin{aligned}
c(E_p) & = & c^2\bigl({\ensuremath{\delta}}^{S_p}(A_p)\bigr) \leq 
4{\ensuremath{\mathrm{vol}}}_p\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^x(U)}{{\ensuremath{\mathcal V}}^{x}(A_p)}\Bigr)\ln\ln\bigl(e(r+1)\bigr), 
\qquad \text{and therefore} \notag \\
c(E_p)+c(Z_p) & \leq &
4\Bigl(\beta|{\ensuremath{\mathcal T}}_p|+{\ensuremath{\mathrm{vol}}}_p\Bigr)\ln\Bigl(\frac{e^dr{\ensuremath{\mathcal V}}^x(U)}{{\ensuremath{\mathit{OPT}}}}\Bigr)\ln\ln\bigl(e(r+1)\bigr). 
\label{indstep}\end{aligned}$$ Note that $\sum_{p=1}^h{\ensuremath{\mathrm{vol}}}_p\leq{\ensuremath{\mathcal V}}^x(U)+\beta(h-1)$ and $\sum_{p=1}^h|{\ensuremath{\mathcal T}}_p|\leq N-h$ (since each time we create a child of $(U,{\ensuremath{\mathcal T}})$ we remove at least one new $(s_i,t_i)$ pair from ${\ensuremath{\mathcal T}}$). So adding over all $p=1,\ldots,h$, we obtain that $c(Z_U)\leq 
4\bigl(\beta|{\ensuremath{\mathcal T}}|+{\ensuremath{\mathcal V}}^x(U)\bigr)\ln\bigl(\frac{e^dr{\ensuremath{\mathcal V}}^x(U)}{{\ensuremath{\mathit{OPT}}}}\bigr)\ln\ln\bigl(e(r+1)\bigr)$.

\[2mcthm\] Algorithm [$2\text{-}$[$\mathsf{MCAlg}$]{}]{}returns a feasible solution of cost at most $O(\ln r\ln\ln r)\cdot{\ensuremath{\mathit{OPT}}}$.

Feasibility of $Z$ is shown in Lemma \[2feas\]. Each time we make a recursive call to [$2\text{-}$[$\mathsf{MCAlg}$]{}]{}, the number of source-sink pairs involved decreases by at least a factor of 2, so the depth $d$ of the overall recursion tree is $O(\log_2 r)$. Since $\beta r+{\ensuremath{\mathcal V}}^x(V)=\bigl(2+\frac{1}{r}\bigr){\ensuremath{\mathit{OPT}}}$ and $\frac{r{\ensuremath{\mathcal V}}^x(V)}{{\ensuremath{\mathit{OPT}}}}=r+1$, by Lemma \[indnlem\], this implies that $c(Z)\leq O({\ensuremath{\mathit{OPT}}})\cdot\bigl(\ln(r+1)+O(\log_2 r)\bigr)\ln\ln\bigl(e(r+1)\bigr)$.

\[2feas\] The solution $Z$ returned by [$2\text{-}$[$\mathsf{MCAlg}$]{}]{}is feasible.

Suppose for a contradiction that there is some $s_i$-$t_i$ pair such that there are (at least) 2 edge-disjoint $s_i$-$t_i$ paths $P_1, P_2$ in $(V,E{\ensuremath{\setminus}}Z)$. Consider the recursion tree of [$2\text{-}$[$\mathsf{MCAlg}$]{}]{}, and let $(Y,{\ensuremath{\mathcal T}}_Y)$ be the node furthest from the root such that $P_1,P_2{\subseteq}E(Y)$. (Such a node must exist since the root satisfies this property.)

Note that there is at least one child $(X,.)$ of $(Y,{\ensuremath{\mathcal T}}_Y)$ such that ${\ensuremath{\delta}}^Y_{P_1\cup P_2}(X)\neq{\ensuremath{\emptyset}}$. If not and both $P_1$ and $P_2$ are contained in some child of $(Y,{\ensuremath{\mathcal T}}_Y)$ then this contradicts the definition of $(Y,{\ensuremath{\mathcal T}}_Y)$. Otherwise, we have $P_1,P_2{\subseteq}E(A)$, where $A=Y{\ensuremath{\setminus}}\bigcup_{\text{children $(X,.)$ of $(Y,{\ensuremath{\mathcal T}}_Y)$}}X$. But then we would have processed $A$ in step A3 and created at least one child $(A',.)$ for some $A'{\subseteq}A$.

We claim that if ${\ensuremath{\delta}}^Y_{P_1\cup P_2}(X)\neq{\ensuremath{\emptyset}}$ for a child $(X,.)$ of $(Y,{\ensuremath{\mathcal T}}_Y)$, then $|{\ensuremath{\delta}}^Y_{P_1\cup P_2}(X)|\geq 2$. This is true if both $s_i$ and $t_i$ are in $X$ or if neither of them are in $X$ since then a path crossing $X$ must cross it at least twice. Otherwise, $X$ is an $s_i$-$t_i$ cut, and since $P_1$ and $P_2$ are edge-disjoint $s_i$-$t_i$ paths in $E(Y)$, we again have $|{\ensuremath{\delta}}^Y_{P_1\cup P_2}(X)|\geq 2$. Among all the children $(X,.)$ of $(Y,{\ensuremath{\mathcal T}}_Y)$ such that ${\ensuremath{\delta}}^Y_{P_1\cup P_2}(X)\neq{\ensuremath{\emptyset}}$, let $(W,.)$ be the child that was added to ${\ensuremath{\mathcal S}}$ earliest in step A3.4 of ${\ensuremath{2\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}(Y,{\ensuremath{\mathcal T}}_Y)$; let $S'{\subseteq}Y$ be the current set $S$ when $W$ was added. Then, $P_1\cup P_2{\subseteq}E(S')$, and so $|{\ensuremath{\delta}}^{S'}_{E{\ensuremath{\setminus}}Z}(W)|\geq |{\ensuremath{\delta}}^{S'}_{P_1\cup P_2}(W')|\geq 2$. But this is a contradiction, since we include in $Z$ all but at most one edge of ${\ensuremath{\delta}}^{S'}(W)$.

#### General $k$ and unit costs.

The algorithm, which we denote by [$k\text{-}$[$\mathsf{MCAlg}$]{}]{}, leading to our bicriteria guarantee is quite similar to [$2\text{-}$[$\mathsf{MCAlg}$]{}]{}. The only changes are the following:

[$\bullet$]{}[=0ex ]{}

In steps A1 and A3.5, we set ${\ensuremath{\mathcal T}}'$ to be the $s_i$-$t_i$ pairs from ${\ensuremath{\mathcal T}}$ that are at least ${\ensuremath{\gamma}}(k-1)$-edge-connected in $(S,E(S))$.

In step A3.2, we apply Corollary \[cor3\] with the set $S$ to find the radius $\rho_1\in[0,1)$.

In step A3.4, we add [*all*]{} edges of ${\ensuremath{\delta}}^S(A)$ to $Z$. (Unlike 2-[$\mathsf{MC}$]{}, if we only include the edges contributing to $c^q\bigl({\ensuremath{\delta}}^S(A)\bigr)$ for some suitable $q$, then we cannot necessarily argue that the final solution satisfies the stated connectivity guarantee.)

Of course, in step A4, we now recursively call [$k\text{-}$[$\mathsf{MCAlg}$]{}]{}(with the same arguments).

\[kmcthm\] For any ${\ensuremath{\gamma}}>1$, algorithm [$k\text{-}$[$\mathsf{MCAlg}$]{}]{}returns a solution $Z$ such that $c(Z)\leq O\bigl(\frac{{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\ln r\ln\ln r\bigr)\cdot{\ensuremath{\mathit{OPT}}}$ and every $s_i$-$t_i$ pair is less than ${\ensuremath{\gamma}}(k-1)$-edge-connected in $(V,E{\ensuremath{\setminus}}Z)$.

Thus, taking ${\ensuremath{\gamma}}=\frac{k}{k-1}$, we obtain a feasible solution of cost at most $O\bigl((k-1)^2\ln r\ln\ln r\bigr)\cdot{\ensuremath{\mathit{OPT}}}$.

Let $Z$ be the output of ${\ensuremath{k\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}\bigl(V,\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$. It is clear that $Z$ is feasible: every $s_i$-$t_i$ pair that is at least ${\ensuremath{\gamma}}(k-1)$-edge-connected in $(U,E(U))$, where $(U,.)$ is a node of the recursion tree is either taken care of (i.e., rendered less than ${\ensuremath{\gamma}}(k-1)$-edge-connected) by the edges added to $Z$ in step A3, or, by induction, is taken care of by a recursive call.

Mimicking the proof of Lemma \[indnlem\], and using Corollary \[cor3\] in place of Lemma \[cor2\] in the proof, one can easily show that if $d$ is the depth of the recursion tree rooted $(U,{\ensuremath{\mathcal T}})$ and $Z_U={\ensuremath{k\text{-}}{\ensuremath{\mathsf{MCAlg}}\xspace}}(U,{\ensuremath{\mathcal T}})$, then $$c(Z_U)\leq
\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\Bigl(\beta|{\ensuremath{\mathcal T}}|+{\ensuremath{\mathcal V}}^x(U)\Bigr)\ln\Bigl(\frac{e^dr{\ensuremath{\mathcal V}}^x(U)}{{\ensuremath{\mathit{OPT}}}}\Bigr)\ln\ln\bigl(e(r+1)\bigr).$$ Since the depth of the overall recursion tree is $O(\log_2 r)$, as argued in the proof of Theorem \[2mcthm\], we obtain that $c(Z)\leq O\bigl(\frac{{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\ln r\ln\ln r\bigr)\cdot{\ensuremath{\mathit{OPT}}}$.

Improved hardness result for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{} {#kstcut}
====================================================================

Theorems \[thm:inapproxSSVE\] and \[kmichard\] together prove that [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}is at least as hard as the densest-$k$-subgraph problem ([$\mathsf{D}k\mathsf{S}$]{}) problem. In [$\mathsf{D}k\mathsf{S}$]{}on hypergraphs, we seek a set of $k$ nodes containing the maximum the number of hyperedges. Our hardness result implies that obtaining a unicriterion $O\bigl(k^{{\ensuremath{\epsilon}}_0}{\operatorname{polylog}}(n)\bigr)$-approximation for some constant ${\ensuremath{\epsilon}}_0$ would improve the current-best guarantee for [$\mathsf{D}k\mathsf{S}$]{}on graphs, and imply the existence of a certain family of one-way functions. Our reduction is from *small set vertex expansion* ([$\mathsf{SSVE}$]{}), wherein we have a bipartite graph $G=(U \cup V,E)$ and a parameter $0 < \alpha \leq 1$, and we seek a subset $S \subseteq U$ with $|S| \geq \alpha |U|$ that minimizes the number of neighbors, $\Gamma(S)$. Chuzoy et al. [@ChuzhoyMVZ12] show that [$\mathsf{SSVE}$]{}reduces to the minimization version of [$\mathsf{D}k\mathsf{S}$]{}, [$\mathsf{Min}$[$\mathsf{D}k\mathsf{S}$]{}]{}, wherein we seek a minimum number of nodes that contain at least $k$ hyperedges. They also show that a $\rho$-approximation for [$\mathsf{Min}$[$\mathsf{D}k\mathsf{S}$]{}]{}on ${\ensuremath{\lambda}}$-uniform hypergraphs yields a $(2\rho^{\ensuremath{\lambda}})$-approximation for [$\mathsf{D}k\mathsf{S}$]{}on ${\ensuremath{\lambda}}$-uniform hypergraphs.

For a graph $H=(V_H,E_H)$, subset $S{\subseteq}V_H$, and $v\in V_H{\ensuremath{\setminus}}S$, we use ${\ensuremath{\delta}}_H(S,v)={\ensuremath{\delta}}_H(v,S)$ to denote the edges between $v$ and $S$ in $H$, and ${\ensuremath{\Gamma}}_H(S)$ to denote the set of neighbors of $S$ in $H$. As is standard, we abbreviate ${\ensuremath{\delta}}_H(\{v\},V_H{\ensuremath{\setminus}}\{v\})$ to ${\ensuremath{\delta}}_H(v)$.

\[thm:inapproxSSVE\] For any $\lambda \geq 2$, there is a polytime approximation-preserving reduction that given a [$\mathsf{Min}$[$\mathsf{D}k\mathsf{S}$]{}]{}-instance on a ${\ensuremath{\lambda}}$-uniform hypergraph with $n$ nodes and $m$ edges, creates an [$\mathsf{SSVE}$]{}-instance with $m+n$ nodes and ${\ensuremath{\lambda}}m$ edges. Hence, a $\rho(m,n)$-approximation for [$\mathsf{SSVE}$]{}yields, for $\lambda$-uniform hypergraphs, a $\rho({\ensuremath{\lambda}}m,m+n)$-approximation for [$\mathsf{Min}$[$\mathsf{D}k\mathsf{S}$]{}]{}and a $\bigl(2\bigl(\rho({\ensuremath{\lambda}}m,m+n)\bigr)^\lambda\bigr)$-approximation for [$\mathsf{D}k\mathsf{S}$]{}.

\[kmichard\] There is a polytime approximation-preserving reduction that given an [$\mathsf{SSVE}$]{}-instance with $n$ nodes and $m$ edges, creates a [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}-instance with $O(n^3)$ nodes, $O(mn^2+n^5)$ edges, and Hence, a $\rho(k,m,n)$-approximation for [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}yields a $\rho\bigl(O(mn^2),O(mn^2+n^5),O(n^3)\bigr)$-approximation for [$\mathsf{SSVE}$]{}.

Given an instance $\bigl(G=(U\cup V,E),{\ensuremath{\alpha}}\bigr)$ of [$\mathsf{SSVE}$]{}, we construct the following instance of [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}; see Fig. \[fig:SSVE\]. Let $N=2|U||V|+1$. Below, an infinite-cost edge $(u,v)$ of capacity $b_{uv}$ is simply a shorthand for $b_{uv}$ parallel infinite-cost edges. Also, unless otherwise specified, an edge has unit capacity.

[()]{}[ =0.5ex =0ex ]{}

We replace each vertex $v \in V$ with a clique $K(v)$ of size $N$. All edges in the clique have infinite cost. For each edge $(u,v) \in E$, we add an edge between $u$ and every vertex in $K(v)$.

We add the source $s$ and connect it to all vertices in $U$; we add the sink $t$ and connect it to all vertices in $K(v)$ for every $v\in V$.

Finally, we add a vertex $b$, an edge $(b,t)$ with capacity $|E| \cdot N$, and edges $(b,u)$ with capacity $|\delta_G(u)| \cdot N$ for all $u\in U$. We also add a vertex $a$, an edge $(s,a)$ with capacity $|E| \cdot N$, and for all $v \in V$, we add edges $(a,x)$ for all $x\in K(v)$ with capacity $|{\ensuremath{\delta}}_G(v)|$.

All edges have infinite cost except for the edges between $\bigcup_{v\in V}K(v)$ and $t$, which have unit cost. We set $k = |U|(1- \alpha) + N|E| +1$.

We claim that there exists a solution of value at most $C$ for the [$\mathsf{SSVE}$]{}instance iff there is a solution of value at most $CN$ for the [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}instance. Note that a solution $F$ consisting of unit-cost edges is feasible if the maximum $s$-$t$ flow in the capacitated remainder graph $G'{\ensuremath{\setminus}}F$ has value at most $k-1$. The intuition is that if we send $N|E|$ units of flow along the paths $s-a-x-u-b-t$ for all $(u,v)\in E, x\in K(v)$, then in the residual digraph, all arcs between $U$ and $\bigcup_{v\in V}K(v)$ leave $U$. Given this, one can mimic the arguments in [@ChuzhoyMVZ12] to show the desired claim.

Suppose the [$\mathsf{SSVE}$]{}instance has a solution $S{\subseteq}U$ of value at most $C$. Construct a [$k$-[$(s,t)\text{-}\mathsf{Cut}$]{}]{}-solution by removing the $(v,t)$ edges for all $v\in\Gamma_G(S)$. Clearly the cost of this set is at most $CN$. We now argue feasibility. Consider the $s$-$t$ cut induced by $A=\{s,a\}\cup S\cup\bigcup_{v\in\Gamma_G(S)}K(v)$. Then $|{\ensuremath{\delta}}_{G'}(A)|$ is equal to $$\underbrace{|U{\ensuremath{\setminus}}S|}_{\normalsize{\substack{\text{edges between}\\ \text{$s$ and $U{\ensuremath{\setminus}}S$}}}}+\
\underbrace{\sum_{u\in S}N|{\ensuremath{\delta}}_G(u)|}_{\substack{\text{edges between}\\ \text{$b$ and
      $S$}}}\ +
\underbrace{\sum_{v\in V{\ensuremath{\setminus}}{\ensuremath{\Gamma}}_G(S)}N|{\ensuremath{\delta}}_G(v)|}_{\substack{\text{edges
    between $a$ and}\\ \text{$\bigcup_{v\in V{\ensuremath{\setminus}}{\ensuremath{\Gamma}}_G(s)}K(v)$}}}\ +\ 
\underbrace{N\bigl(\text{\# of edges in $G$ between $U{\ensuremath{\setminus}}S$ and
    ${\ensuremath{\Gamma}}_G(S)$}\bigr)}_{\text{edges between $U{\ensuremath{\setminus}}S$ and  
    $\bigcup_{v\in{\ensuremath{\Gamma}}_G(S)}K(v)$}}.$$ The sum of the last two terms is $\sum_{u\in U{\ensuremath{\setminus}}S}N|{\ensuremath{\delta}}_G(u)|$ and $|U{\ensuremath{\setminus}}S|\leq|U|(1-{\ensuremath{\alpha}})$, so the size of the $s$-$t$ cut is at most $|U|(1-{\ensuremath{\alpha}})+\sum_{u\in U}N|{\ensuremath{\delta}}_G(u)|\leq|U|(1-{\ensuremath{\alpha}})+N|E|\leq k-1$.

For the other direction, suppose $G'$ has a solution $F$ of value at most $CN$. Clearly, $F$ can consist of only unit-cost edges (incident to $t$). We first argue that we may convert $F$ into a structured feasible solution $F'$ of cost at most $CN$ where $|F'\cap{\ensuremath{\delta}}_{G'}(K(v))|\in\{0,N\}$ for all $v\in V$.

Fix $v\in V$. If $|\delta_{G'}(K(v),t){\ensuremath{\setminus}}F|\leq k':=|U|(1- \alpha)$, then we add all edges of ${\ensuremath{\delta}}_{G'}(K(v),t){\ensuremath{\setminus}}F$ to $F$. Now suppose $|\delta_{G'}(K(v),t){\ensuremath{\setminus}}F|>k'$ and let $w_1,\dots,w_{k'+1}$ be vertices in $K(v)$ such that $(w_i,t)\notin F$ for all $i=1,\ldots,k'+1$. We claim that $F{\ensuremath{\setminus}}{\ensuremath{\delta}}_{G'}(K(v),t)$ is also feasible. Suppose to the contrary that we now have $k'+1$ edge-disjoint $s$-$t$-paths. We may assume that each such path contains at most one vertex from $K(v)$ since we can always shortcut the path to $t$. Consider a path $P$ that contains a vertex $w\in K(v)$ where $w\notin\{w_1, \dots, w_{k'+1}\}$. Then we can construct another path $P'$ by switching $w$ with a distinct vertex $w_j$ for some $j\in\{1\dots, k'+1\}$. Note that $P'$ is an $s$-$t$ path that avoids edges in $F$. If we repeat this argument for all such paths, we obtain $k'+1$ edge-disjoint $s$-$t$-paths in $G'\setminus F$, contradicting the feasibility of $F$.

If we perform the above transformation for all $v\in V$, then we obtain a feasible solution $F'$ of cost at most $|F|+|V|k'<(C+1)N$. But by construction $|F'|$ must be an integer multiple of $N$, so $|F'|\leq CN$.

Consider the residual network $\tilde G$ that is obtained from $G'{\ensuremath{\setminus}}F'$ as follows. We first bidirect the edges of $G'{\ensuremath{\setminus}}F'$, giving each resulting arc the same capacity as that of the corresponding edge of $G'$. ${\ensuremath{\tilde G}}$ is the residual network obtained after we send one unit of flow along the path $s$-$a$-$x$-$u$-$b$-$t$ for every edge $(u,v)\in E$ and every $x\in K(v)$. Note that these paths are edge disjoint, so we send $N|E|$ units of flow. By flow theory [@AhujaMO93], we know that the value of the maximum $s$-$t$-flow in $G'{\ensuremath{\setminus}}F'$ is at most $k':=|U|(1-{\ensuremath{\alpha}})$ iff the value of maximum $s$-$t$-flow in ${\ensuremath{\tilde G}}$, which equals the capacity of the minimum $s$-$t$ cut in ${\ensuremath{\tilde G}}$, is at most $k'$. It follows there is an $s$-$t$ cut in ${\ensuremath{\tilde G}}$ of capacity at most $k'$.

Let $A$ be the vertices that are on the $s$-side of this cut. Let $S := U \cap A$. Then $|S|\geq\alpha|U|$, otherwise the cut would have capacity more than $|U|(1- \alpha)$ due to the arcs $(s,u)$ for $u\in U\setminus S$. Consider $v\in{\ensuremath{\Gamma}}_G(S)$. We must have $K(v){\subseteq}A$: if $K(v)\cap A={\ensuremath{\emptyset}}$, then considering $u\in S$ such that $(u,v)\in E$, the cut has capacity at least $N>k'$ due to the edges between $u$ and $K(v)$; otherwise, since $K(v)$ is split between the $s$- and $t$- sides, the cut has capacity at least $N-1>k'$. Finally, ${\ensuremath{\delta}}_{G'}(K(v),t){\subseteq}F'$, otherwise ${\ensuremath{\delta}}_{G'}(K(v),t)\cap F'={\ensuremath{\emptyset}}$, and again the cut has capacity at least $N>k'$. Thus, $|{\ensuremath{\Gamma}}_G(S)|\leq C$, so $S$ is an [$\mathsf{SSVE}$]{}-solution.

Extensions to node-connectivity versions of [$k$-[$\mathsf{MC}$]{}]{} {#node-extn}
=====================================================================

We now consider variants of [$k$-[$\mathsf{MC}$]{}]{}where we seek to delete edges or nodes so as to reduce the [*node connectivity*]{} of each $s_i$-$t_i$ pair to at most $k-1$. Formally, as before, we are given an undirected graph $G=(V,E)$, $r$ source-sink pairs $(s_1,t_1),\ldots,(s_r,t_r)$, and an integer $k\geq 1$. In the [*edge-deletion $k$-route node-multicut*]{} ([$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}) problem, we have nonnnegative edge-costs $\{c_e\}_{e\in E}$ and we seek a minimum-cost set $F{\subseteq}E$ of edges such that the remainder graph ${\ensuremath{{\ensuremath{\overline{G}}}}}=(V,E{\ensuremath{\setminus}}F)$ contains at most $k-1$ node-disjoint $s_i$-$t_i$ paths for all $i=1,\ldots,r$. In the [*node-deletion $k$-route node-multicut*]{} ([$\mathsf{ND}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}) problem, we have nonnegative node costs $\{c_v\}_{v\in V}$ and we seek a minimum-cost set $A{\subseteq}V{\ensuremath{\setminus}}\{s_1,t_1,\ldots,s_r,t_r\}$ of nodes such that the remainder graph ${\ensuremath{{\ensuremath{\overline{G}}}}}=\bigl(V{\ensuremath{\setminus}}A,E(V{\ensuremath{\setminus}}A)\bigr)$ contains at most $k-1$ node-disjoint $s_i$-$t_i$ paths for all $i=1,\ldots,r$.

The LP-relaxations of these $k$-route node-multicut problems induce both edge and node lengths, so to round these we develop region-growing lemmas that also incorporate node lengths. To keep notation simple, instead of proving an overly- general region-growing lemma and obtaining the lemmas required for [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}and [$\mathsf{ND}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}as corollaries, we specifically focus on [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}(Section \[edknmc\]) and [$\mathsf{ND}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}(Section \[ndknmc\]) and prove suitable region-growing lemas. We use these to obtain an $O(\ln r\ln\ln r)$-approximation for [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}with $k=2$, and a bicriteria $\bigl({\ensuremath{\gamma}},O\bigl(\frac{{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\ln r\ln\ln r\bigr)\bigr)$-approximation for [$\mathsf{ND}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}with general $k$ and unit node costs.

Edge-deletion $k$-route node-multicut ([$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}) with $k=2$ {#edknmc}
-----------------------------------------------------------------------------------------------

The LP-relaxation for [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}is as follows. $$\begin{aligned}
{3}
\min & \quad & \sum_e c_ex_e & \tag{P2} \label{edknmcp} \\
\text{s.t.} && \sum_{e\in E(P)} x_e &+\sum_{v\in V(P)}y^i_v \geq 1 \qquad && 
\forall i, \forall P\in{\ensuremath{\mathcal P}}_i \\
&& \sum_v y^i_v & \leq k-1, \quad y^i_{s_i}=y^i_{t_i}=0 \qquad && \forall i \\
&& x, y & \geq 0. \end{aligned}$$

#### Region-growing lemma.

Let $\bigl(x,\{y^i\}\bigr)$ be an optimal solution to , and ${\ensuremath{\mathit{OPT}}}$ be its value. Let $S{\subseteq}V$ represent the node-set of the current region. For $T{\subseteq}S{\subseteq}V$, recall that $E(S)$ is the set of edges with both endpoints in $S$ and ${\ensuremath{\delta}}^S(T)$ is the set of boundary edges of $T$ in $E(S)$. Set $\beta={\ensuremath{\mathit{OPT}}}/r$. As before, define ${\ensuremath{\mathcal V}}^x(S):=\beta+\sum_{e\in E(S)}c_ex_e$. Let $\rho\geq 0$. Let $z\in V$. Fix a commodity $i$.

[$\bullet$]{}[=0ex =0.5ex ]{}

Define ${\ell^i}(u;v)=\min_{P:P\text{ is a $u$-$v$ path}}
\bigl(\sum_{e\in E(P)}x_e+\sum_{w\in V(P):w\neq u}y^i_w\bigr)$, where $E(P)$ and $V(P)$ denote respectively the set of edges and nodes of $P$. Note that ${\ell^i}$ defines an [*asymmetric metric*]{} on $V\times V$.

Define $B_i(z,\rho):=\{v\in V: {\ell^i}(z;v)\leq\rho\}$ to be the ball of radius $\rho$ around $z$. Let $B^S_i(z,\rho):=B_i(z,\rho)\cap S$.

Define the edge-boundary ${\ensuremath{\partial}}^{S,x}_i$, and node-boundary ${\ensuremath{\Gamma}}^{S,y}_i$, of $B^S_i(z,\rho)$ in $S$ as follows. $$\begin{aligned}
{\ensuremath{\partial}}^{S,x}_i(z,\rho) & := & \{(u,v)\in E: u,v\in S,\ {\ell^i}(z;u)\leq\rho,\ {\ell^i}(z;v)-y^i_v>\rho\} \\
{\ensuremath{\Gamma}}^{S,y}_i(z,\rho) & := & \{v\in S: \rho<{\ell^i}(z;v)\leq\rho+y^i_v\}.\end{aligned}$$ Let ${\ensuremath{\overline{B^S_i}}}(z,\rho):=S{\ensuremath{\setminus}}\bigl(B^S_i(z,\rho)\cup{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)\bigr)$.

Define the following volumes: $$\begin{aligned}
{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho) & := & \beta
+\sum_{\substack{(u,v)\in E: u\in B^S_i(z,\rho) \\ v\in B^S_i(z,\rho)\cup{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)}}c_{uv}x_{uv}
+\sum_{\substack{(u,v)\in{\ensuremath{\partial}}^{S,x}_i(z,\rho): \\ u\in B^S_i(z,\rho)}}c_{uv}\bigl(\rho-\ell^i(z;u)\bigr) \\[0.5ex]
{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho) & := & \beta
+\sum_{\substack{(u,v)\in E: u\in{\ensuremath{\overline{B^S_i}}}(z,\rho) \\ v\in{\ensuremath{\overline{B^S_i}}}(z,\rho)\cup{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)}}c_{uv}x_{uv}
+\sum_{\substack{(u,v)\in{\ensuremath{\partial}}^{S,x}_i(z,\rho): \\ u\in B^S_i(z,\rho)}}c_{uv}\bigl(\ell^i(z;v)-\rho-y^i_v\bigr)\end{aligned}$$

\[ednreggrow\] Let $S{\subseteq}V$, $z\in V$, $i$ be some commodity, and $0\leq a<b$. Let $\rho$ be chosen uniformly at random from $[a,b)$. Then, $$\begin{aligned}
{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho)\bigr)}
{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{S,x}_i(z,b)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho)}\Bigr)}\right]}}
& \leq & \frac{1}{b-a}\cdot\ln\ln\biggl(\frac{e{\ensuremath{\mathcal V}}^{S,x}_i(z,b)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,a)}\biggr) 
\quad\text{and} \label{ednbnd1} \\
{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho)\bigr)}
{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho)\ln\Bigl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho)}\Bigr)}\right]}}
& \leq &
\frac{1}{b-a}\cdot\ln\ln\biggl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,b)}\biggr) 
\label{ednbnd2}\end{aligned}$$

We abbreviate $c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho)\bigr)$ to $c(\rho)$, ${\ensuremath{\mathcal V}}^{S,x}_i(z,\rho)$ to ${\ensuremath{\mathcal V}}(\rho)$ and ${\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho)$ to ${\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)$. Let $I=\{\ell^i(z;v), \ell^i(z;v)-y^i_v: v\in V\}$. Note that ${\ensuremath{\mathcal V}}(\rho)$ and ${\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)$ are differentiable at all $\rho\in[a,b){\ensuremath{\setminus}}I$ and for each such $\rho$, we have $\frac{d{\ensuremath{\mathcal V}}(\rho)}{d\rho}=c(\rho)$ and $\frac{d{\ensuremath{\overline{{\ensuremath{\mathcal V}}}}}(\rho)}{d\rho}=-c(\rho)$. The proof now follows from exactly the same arguments as in the proof of Lemma \[reggrow\].

\[edncor2\] Let $S{\subseteq}V$, $z\in V$, and $i$ be some commodity. Let ${\ensuremath{\alpha}}\in(0,1)$ and $q={\ensuremath{\left\lceil\frac{k-1}{{\ensuremath{\alpha}}}\right\rceil}}$. We can efficiently find $\rho_1\in[0,1)$ such that $$\begin{aligned}
c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho_1)\bigr) & \leq &
\frac{2}{1-{\ensuremath{\alpha}}}\cdot{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{x}(S)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr) \\ 
c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho_1)\bigr) & \leq &
\frac{2}{1-{\ensuremath{\alpha}}}\cdot{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{x}(S)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr) \\ 
\bigl|{\ensuremath{\Gamma}}^{S,y}_i(z,\rho_1)\bigr| & < & q. \end{aligned}$$

If we choose $\rho$ uniformly at random from $[0,1)$ then ${\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\bigl|{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)\bigr|\right]}}\leq\sum_iy^i_v\leq k-1$. Taking $[a,b)=[0,1)$, the arguments in Corollary \[cor1\] readily generalize to show that we can efficiently find $\rho_1\in[0,1)$ such that ${c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho_1)\bigr)}/
{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)\ln\bigl(\frac{e{\ensuremath{\mathcal V}}^{S,x}_i(z,1)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)}\bigr)}$, and ${c\bigl({\ensuremath{\partial}}^{S,x}_i(z,\rho_1)\bigr)}/
{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)\ln\bigl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,0)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)}\bigr)}$ are at most $\frac{2}{(1-{\ensuremath{\alpha}})}$ times the right-hand-sides of and respectively, and $\bigl|{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)\bigr|<\frac{k-1}{{\ensuremath{\alpha}}}$. The lemma now follows by noting that $${\ensuremath{\mathcal V}}^{S,x}_i(z,1)\leq{\ensuremath{\mathcal V}}^x(S), \qquad {\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,0)\leq{\ensuremath{\mathcal V}}^x(S), \qquad
\frac{{\ensuremath{\mathcal V}}^{S,x}_i(z,1)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,0)}\leq r+1, \qquad
\frac{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,0)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,1)}\leq r+1. 
\tag*{\qedsymbol}$$

#### Algorithm and analysis.

The algorithm and analysis dovetail the one in Section \[kmcalg\] for $k=2$.

------------------------------------------------------------------------

[**Algorithm**]{} ${\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}(U, {\ensuremath{\mathcal T}}=\{(s_1,t_1),\ldots,(s_{N},t_{N})\}$)

[**Input:**]{} A subset $U{\subseteq}V$, and a collection ${\ensuremath{\mathcal T}}=\{(s_i,t_i)\}_{i=1}^{N}$ of source-sink pairs, where $s_i,t_i\in U$ for all $i=1,\ldots, N$.

[**Output:**]{} A set $Z{\subseteq}E(U)$ such that $s_i$ and $t_i$ are at most 1-node-connected in $(U,E(U){\ensuremath{\setminus}}Z)$ for all $i=1,\ldots,N$.

[B.]{}[ =0.5ex =0.25ex ]{}

Set $S=U$, $Z={\ensuremath{\emptyset}}$, ${\ensuremath{\mathcal S}}={\ensuremath{\emptyset}}$, and ${\ensuremath{\mathcal T}}'=\bigl\{(s_i,t_i)\in{\ensuremath{\mathcal T}}: \text{$s_i$ and $t_i$ are at least $2$-node-connected in $(S,E(S))$}\bigr\}$.

If ${\ensuremath{\mathcal T}}'={\ensuremath{\emptyset}}$, return $Z$.

While ${\ensuremath{\mathcal T}}'\neq{\ensuremath{\emptyset}}$, we do the following.

[A.]{}[=0.25ex =0ex ]{}

Pick some $(s_i,t_i)\in{\ensuremath{\mathcal T}}'$.

Apply Corollary \[edncor2\] with $z=s_i$, ${\ensuremath{\alpha}}=0.5$ and the set $S$ (and $k=2$) to find a radius $0\leq\rho_1<1$.

If $B^S_i(s_i,\rho_1)\cup{\ensuremath{\Gamma}}^{S,y}_i(s_i,\rho_1)$ contains at most $N/2$ pairs from ${\ensuremath{\mathcal T}}$ then set $A=B^S_i(s_i,\rho_1)$, else set $A={\ensuremath{\overline{B^S_i}}}(s_i,\rho_1)$.

Set ${\ensuremath{\mathcal S}}{\ensuremath{\leftarrow}}{\ensuremath{\mathcal S}}\cup\{A\cup{\ensuremath{\Gamma}}^{S,y}_i(s_i,\rho_1)\}$. Add the edges in ${\ensuremath{\partial}}^{S,x}_i(s_i,\rho_1)$ to $Z$.

Set $S{\ensuremath{\leftarrow}}S{\ensuremath{\setminus}}A$. Update ${\ensuremath{\mathcal T}}'$ to be the $s_i$-$t_i$ pairs from ${\ensuremath{\mathcal T}}$ that are at least $2$-node-connected in $(S,E(S))$.

For every set $A\in{\ensuremath{\mathcal S}}$, set $Z{\ensuremath{\leftarrow}}Z\cup{\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}\bigl(A,\{(s_i,t_i)\in{\ensuremath{\mathcal T}}: s_i,t_i\in A\}\bigr)$.

Return $Z$.

The initial call to ${\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}$, which computes the solution we return, is ${\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}\bigl(V,\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$.

------------------------------------------------------------------------

Let $Z:={\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}\bigl(V,\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$. Define the depth of a subtree of the recursion tree corresponding to the execution of [$\mathsf{ED}\text{-}2\text{-}$[$\mathsf{NMCAlg}$]{}]{}to be the maximum number of edges on a root to leaf path of the subtree. The following claim will be useful to prove feasibility of $Z$.

\[nodecon\] Let $S, T{\subseteq}V$ with $|S\cap T|\leq 1$. Let $E_S{\subseteq}E(S)$ and $E_T{\subseteq}E(T)$. Let $u,v\in S$ be such that $u$ and $v$ are at most 1-node-connected in $(S,E_S)$. Then, $u$ and $v$ are at most 1-node-connected in $(S\cup T,E_S\cup E_T)$.

If $S\cap T={\ensuremath{\emptyset}}$, this clearly holds. So assume otherwise. Suppose $P_1$, $P_2$ are two simple node-disjoint $u$-$v$ paths in $G'=(S\cup T,E_S\cup E_T)$. At least one of $P_1$ and $P_2$ does not lie completely in $(S,E_S)$; suppose $P_2$ is this path. But since all edges of ${\ensuremath{\delta}}_{G'}(S)$ are incident to a single node, and $P_2$ both exits and leaves $S$, $P_2$ contains a repeated node, which is a contradiction.

\[edn2feas\] The solution $Z$ returned is feasible.

Suppose for a contradiction that there is some $s_i$-$t_i$ pair that is (at least) 2-node-connected in $(V,E{\ensuremath{\setminus}}Z)$. Consider the recursion tree of [$\mathsf{ED}\text{-}2\text{-}$[$\mathsf{NMCAlg}$]{}]{}, and let $(Y,{\ensuremath{\mathcal T}}_Y)$ be the node furthest from the root such that $s_i$ and $t_i$ are at least 2-node-connected in the subgraph $(Y,E(Y))$ induced by $Y$. Suppose that the loop in step B3 of ${\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}(Y,{\ensuremath{\mathcal T}}_Y)$ runs for $h$ iterations. Note that $h\geq 1$ since $s_i$ and $t_i$ are at least 2-node-connected in $(Y,E(Y))$. Let $X_p$ be the set added to ${\ensuremath{\mathcal S}}$ in step B3.4 in the $p$-th iteration of the loop. Let $X_{h+1}{\subseteq}Y$ be the set $S$ at the termination of the loop. Let $S_p=\bigcup_{q=p}^{h+1}X_q$ (so $S_1=Y$). Notice that $|X_p\cap S_{p+1}|\leq 1$, since $X_p\cap S_{p+1}{\subseteq}{\ensuremath{\Gamma}}^{S_p,y}_p(s_p,\rho_p)$, where $s_p$-$t_p$ is the source-sink pair and $\rho_p$ is the radius chosen in the $p$-th iteration, and $|{\ensuremath{\Gamma}}^{S_p,y}_p(s_p,\rho_p)|<2$ by Lemma \[edncor2\].

Let $p$ be the highest index such that $s_i$ and $t_i$ are at least 2-node-connected in $(S_p,E(S_p))$. Note that $p\leq h$, otherwise the loop in step B3 would not have terminated with $S=X_{h+1}$. If $s_i, t_i\in S_{p+1}$, they are at most 1-node-connected in $(S_{p+1},E(S_{p+1}))$. Since $|X_p\cap S_{p+1}|\leq 1$, we have $E(S_p)=E(X_p)\cup E(S_{p+1})$, and by Claim \[nodecon\] it follows that $s_i$ and $t_i$ are at most 1-node-connected in $(S_p,E(S_p))$, which is a contradiction. If $s_i, t_i\in X_p$, they are at most 1-node-connected in $(X_p,E(X_p))$ due to the definition of $(Y,{\ensuremath{\mathcal T}}_Y)$, and so we arrive at the same contradiction. So it must be that $\bigl|\{s_i,t_i\}\cap (X_p{\ensuremath{\setminus}}S_{p+1})\bigr|=1$. But then all $s_i$-$t_i$ paths in $(S_p,E(S_p))$ contain the singleton node in $X_p\cap S_{p+1}$. So $s_i$ and $t_i$ are at most 1-node-connected in $(S_p,E(S_p))$, and we have the same contradiction.

\[ednindnlem\] Let $d$ be the depth of a subtree of the recursion tree rooted at $\bigl(U{\subseteq}V,{\ensuremath{\mathcal T}}{\subseteq}\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$. Let $Z_U={\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}(U,{\ensuremath{\mathcal T}})$. Then $c(Z_U)\leq 
4\bigl(\beta|{\ensuremath{\mathcal T}}|+{\ensuremath{\mathcal V}}^x(U)\bigr)\ln\bigl(\frac{e^dr{\ensuremath{\mathcal V}}^x(U)}{{\ensuremath{\mathit{OPT}}}}\bigr)\ln\ln\bigl(e(r+1)\bigr)$.

When $d=0$, we have $Z_U={\ensuremath{\emptyset}}$, so the statement holds. Suppose in step B3 of ${\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}(U,{\ensuremath{\mathcal T}})$, we add sets $A_1,\ldots,A_h$ to ${\ensuremath{\mathcal S}}$ (where $h\geq 1$), in that order. For $p=1,\ldots,h$, let $S_p$ be the current set $S$ when $A_p$ was added to ${\ensuremath{\mathcal S}}$ in step B3.4, and let $E_p$ be the edge-set added to $Z$ in this step. Let $Z_1,\ldots,Z_h$ be the edge-sets returned by the recursive calls to ${\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}(A_1,{\ensuremath{\mathcal T}}_1),\ldots,{\ensuremath{\mathsf{ED}\text{-}2\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}(A_h,{\ensuremath{\mathcal T}}_h)$ in step B4.

Let ${\ensuremath{\mathrm{vol}}}_p={\ensuremath{\mathcal V}}^{S_p,x}_i(s_i,\rho_1)$ if $A_p=B^{S_p}_i(s_i,\rho_1)\cup{\ensuremath{\Gamma}}^{S_p,y}_i(s_i,\rho_1)$ and ${\ensuremath{\mathrm{vol}}}_p={\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S_p,x}_i}}}(s_i,\rho_1)$ if $A_p={\ensuremath{\overline{B^{S_p}_i}}}(s_i,\rho_1)\cup{\ensuremath{\Gamma}}^{S_p,y}_i(s_i,\rho_1)$. The key thing to note is that we still have ${\ensuremath{\mathcal V}}^x(A_p)\leq{\ensuremath{\mathrm{vol}}}_p\leq{\ensuremath{\mathcal V}}^x(S_p)\leq{\ensuremath{\mathcal V}}^x(U)$ and $\sum_{p=1}^h{\ensuremath{\mathrm{vol}}}_p\leq{\ensuremath{\mathcal V}}^x(U)+\beta(h-1)$. The latter follows since an easy induction argument shows that $\sum_{p=q}^h{\ensuremath{\mathrm{vol}}}_p\leq{\ensuremath{\mathcal V}}^x(S_q)+\beta(h-q)$ for all $q=1,\ldots,h$. Given this, the rest of the proof is identical to that of Lemma \[indnlem\].

Each recursive call to [$\mathsf{ED}\text{-}2\text{-}$[$\mathsf{NMCAlg}$]{}]{}, reduces the number of source-sink pairs involved by a factor of at least 2, so the depth $d$ of the entire recursion tree is $O(\log_2 r)$. So we have shown the following.

\[edn2mcthm\] Algorithm [$\mathsf{ED}\text{-}2\text{-}$[$\mathsf{NMCAlg}$]{}]{}returns a feasible solution of cost at most $O(\ln r\ln\ln r)\cdot{\ensuremath{\mathit{OPT}}}$.

Node-deletion $k$-route node-multicut ([$\mathsf{ND}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}) with unit costs {#ndknmc}
----------------------------------------------------------------------------------------------------

The LP-relaxation for [$\mathsf{ND}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}is as follows. $$\begin{aligned}
{3}
\min & \quad & \sum_v c_vx_v & \tag{P3} \label{ndknmcp} \\
\text{s.t.} && \sum_{v\in V(P)}(x_v &+y^i_v) \geq 1 \qquad && 
\forall i, \forall P\in{\ensuremath{\mathcal P}}_i \\
&& \sum_v y^i_v & \leq k-1, \quad y^i_{s_i}=y^i_{t_i}=0 \qquad && \forall i \\
&& x, y & \geq 0, \qquad x_v=0 \quad && \forall v\in\{s_1,t_1,\ldots,s_r,t_r\}.\end{aligned}$$

#### Region-growing lemma.

Let $\bigl(x,\{y^i\}\bigr)$ be an optimal solution to , and ${\ensuremath{\mathit{OPT}}}$ be its value. As before, let $S{\subseteq}V$ represent the node-set of the current region. Set $\beta={\ensuremath{\mathit{OPT}}}/r$. Let $z\in V$ and $\rho\geq 0$. Fix a commodity $i$.

[$\bullet$]{}[=0ex =0.5ex ]{}

Define $\ell^i(u;v)=\min_{P:P\text{ is a $u$-$v$ path}}\sum_{w\in V(P):w\neq u}(x_w+y^i_w)$, where $V(P)$ is the set of nodes of $P$. As before, $\ell^i$ defines an asymmetric metric on $V\times V$.

Define $B_i(z,\rho):=\{v\in V: \ell^i(z;v)\leq\rho\}$, and $B^S_i(z,\rho):=B_i(z,\rho)\cap S$.

Define the $x$-boundary of $B^S_i(z,\rho)$ in $S$ to be ${\ensuremath{\Gamma}}^{S,x}_i(z,\rho):=\{v\in S: \rho+y^i_v<\ell^i(z;v)\leq\rho+x_v+y^i_v\}$. Define the $y$-boundary of $B^S_i(z,\rho)$ in $S$ to be ${\ensuremath{\Gamma}}^{S,y}_i(z,\rho):=\{v\in S: \rho<\ell^i(z;v)\leq\rho+y^i_v\}$. Note that ${\ensuremath{\Gamma}}^{S,x}_i(z,\rho)$ and ${\ensuremath{\Gamma}}^{S,y}_i(z,\rho)$ partition ${\ensuremath{\Gamma}}^S_i(z,\rho):=
\{v\in S{\ensuremath{\setminus}}B^S_i(z,\rho): \exists u\in B^s_i(z,\rho)\text{ s.t. }(u,v)\in E\}$. Let ${\ensuremath{\overline{B^S_i}}}(z,\rho):=S{\ensuremath{\setminus}}\bigl(B^S_i(z,\rho)\cup{\ensuremath{\Gamma}}^{S}_i(z,\rho)\bigr)$.

Define the following volumes: $$\begin{aligned}
{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho) & := & \beta
+\sum_{u\in B^S_i(z,\rho)\cup{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)}c_{u}x_{u}
+\sum_{u\in{\ensuremath{\Gamma}}^{S,x}_i(z,\rho)}c_{u}\bigl(\rho-(\ell^i(z;u)-x_u-y^i_u)\bigr) \\[0.5ex]
{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho) & := & \beta
+\sum_{u\in{\ensuremath{\overline{B^S_i}}}(z,\rho)\cup{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)}c_{u}x_{u}
+\sum_{u\in{\ensuremath{\Gamma}}^{S,x}_i(z,\rho)}c_u(\ell^i(z;u)-y^i_u-\rho)\end{aligned}$$

The following lemma is analogous to Lemma \[ednreggrow\] and follows from the same reasoning.

\[ndnreggrow\] Let $S{\subseteq}V$, $z\in V$, $i$ be some commodity, and $0\leq a<b$. Let $\rho$ be chosen uniformly at random from $[a,b)$. Then, $$\begin{aligned}
{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c\bigl({\ensuremath{\Gamma}}^{S,x}_i(z,\rho)\bigr)}
{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{S,x}_i(z,b)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho)}\Bigr)}\right]}}
& \leq & \frac{1}{b-a}\cdot\ln\ln\biggl(\frac{e{\ensuremath{\mathcal V}}^{S,x}_i(z,b)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,a)}\biggr) 
\quad\text{and} \\ 
{\ensuremath{{\textstyle \operatorname*{E}_{\rho}}\left[\frac{c\bigl({\ensuremath{\Gamma}}^{S,x}_i(z,\rho)\bigr)}
{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho)\ln\Bigl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho)}\Bigr)}\right]}}
& \leq &
\frac{1}{b-a}\cdot\ln\ln\biggl(\frac{e{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,a)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,b)}\biggr) \end{aligned}$$

\[ndncor3\] Let $S{\subseteq}V$. Suppose that $s_i,t_i\in S$ and there are ${\ensuremath{\gamma}}(k-1)$ node-disjoint $s_i$-$t_i$ paths internal to $S$, where ${\ensuremath{\gamma}}>1$. Suppose that $c_v=1$ for all nodes $v$. We can efficiently find $\rho_1\in[0,1)$ such that $$\begin{aligned}
c\bigl({\ensuremath{\Gamma}}^S_i(s_i,\rho_1)\bigr)\leq
\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\cdot{\ensuremath{\mathcal V}}^{S,x}_i(s_i,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^x(S)}{{\ensuremath{\mathcal V}}^{S,x}_i(s_i,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr) \\
c\bigl({\ensuremath{\Gamma}}^S_i(s_i,\rho_1)\bigr)\leq
\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\cdot{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(s_i,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^x(S)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(s_i,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr)\end{aligned}$$

Let ${\ensuremath{\alpha}}\in(0,1)$, whose value we will fix later. Take $[a,b)=[0,1)$. We have $${\ensuremath{\mathcal V}}^{S,x}_i(z,1)\leq{\ensuremath{\mathcal V}}^x(S), \qquad {\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,0)\leq{\ensuremath{\mathcal V}}^x(S), \qquad
\frac{{\ensuremath{\mathcal V}}^{S,x}_i(z,1)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,0)}\leq r+1, \qquad
\frac{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,0)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,1)}\leq r+1.$$ Given this, the arguments in Corollary \[cor1\] readily generalize to show that we can efficiently find $\rho_1\in[0,1)$ such that $$\begin{aligned}
c\bigl({\ensuremath{\Gamma}}^{S,x}_i(s_i,\rho_1)\bigr) & \leq & 
\frac{2}{1-{\ensuremath{\alpha}}}\cdot{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{x}(S)}{{\ensuremath{\mathcal V}}^{S,x}_i(z,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr) \label{ndnxbound1} \\
c\bigl({\ensuremath{\Gamma}}^{S,x}_i(s_i,\rho_1)\bigr) & \leq &
\frac{2}{1-{\ensuremath{\alpha}}}\cdot{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)\ln\Bigl(\frac{e{\ensuremath{\mathcal V}}^{x}(S)}{{\ensuremath{\overline{{\ensuremath{\mathcal V}}^{S,x}_i}}}(z,\rho_1)}\Bigr)\ln\ln\bigl(e(r+1)\bigr) \label{ndnxbound2} \\
\bigl|{\ensuremath{\Gamma}}^{S,y}_i(z,\rho)\bigr| & < & \frac{k-1}{{\ensuremath{\alpha}}}. \label{ndnybound}\end{aligned}$$

Note that $t_i\notin{\ensuremath{\Gamma}}^S_i(s_i,\rho_1)$ since $\rho_1<1$. So removing ${\ensuremath{\Gamma}}^S_i(s_i,\rho_1)$ disconnects $s_i$ and $t_i$, and hence, $|{\ensuremath{\Gamma}}^S_i(s_i,\rho_1)|\geq{\ensuremath{\gamma}}(k-1)$. Therefore, since we have unit node costs and ${\ensuremath{\Gamma}}^{S,x}_i(s_i,\rho_1)$ and ${\ensuremath{\Gamma}}^{S,y}_i(s_i,\rho_1)$ partition ${\ensuremath{\Gamma}}^S_i(s_i,\rho_1)$, we have $c\bigl({\ensuremath{\Gamma}}^S_i(s_i,\rho_1)\bigr)\leq c\bigl({\ensuremath{\Gamma}}^{S,x}_i(s_i,\rho_1)\bigr)\cdot\frac{{\ensuremath{\gamma}}}{{\ensuremath{\gamma}}-1/{\ensuremath{\alpha}}}$. Plugging in the bounds from , , we see that the constant factor multiplying the volume terms on the RHS is minimized by setting ${\ensuremath{\alpha}}={\ensuremath{\gamma}}^{-1/2}$, which yields the constant factor $\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}$.

#### Algorithm and analysis.

The algorithm, [$\mathsf{ND}\text{-}k\text{-}$[$\mathsf{NMCAlg}$]{}]{}, for [$\mathsf{ND}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}is quite similar to [$\mathsf{ED}\text{-}2\text{-}$[$\mathsf{NMCAlg}$]{}]{}. The only changes are the following:

[$\bullet$]{}[=0ex ]{}

In steps B1 and B3.5, we set ${\ensuremath{\mathcal T}}'$ to be the $s_i$-$t_i$ pairs from ${\ensuremath{\mathcal T}}$ that are at least ${\ensuremath{\gamma}}(k-1)$-node-connected in $(S,E(S))$.

In step B3.2, we apply Corollary \[ndncor3\] with the set $S$ to find the radius $\rho_1\in[0,1)$.

In step B3.4, we add $A$ to ${\ensuremath{\mathcal S}}$, and add all nodes of ${\ensuremath{\Gamma}}^S_i(s_i,\rho_1)$ to $Z$.

Of course, in step B4, we now recursively call [$\mathsf{ND}\text{-}k\text{-}$[$\mathsf{NMCAlg}$]{}]{}(with the same arguments).

\[ndknmcthm\] For any ${\ensuremath{\gamma}}>1$, algorithm [$\mathsf{ND}\text{-}k\text{-}$[$\mathsf{NMCAlg}$]{}]{}returns a solution $Z$ such that $c(Z)\leq O\bigl(\frac{{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\ln r\ln\ln r\bigr)\cdot{\ensuremath{\mathit{OPT}}}$ and every $s_i$-$t_i$ pair is less than ${\ensuremath{\gamma}}(k-1)$-node-connected in $(V{\ensuremath{\setminus}}Z,E(V{\ensuremath{\setminus}}Z))$.

Thus, taking ${\ensuremath{\gamma}}=\frac{k}{k-1}$, we obtain a feasible solution of cost at most $O\bigl((k-1)^2\ln r\ln\ln r\bigr)\cdot{\ensuremath{\mathit{OPT}}}$.

Let $Z$ be the output of ${\ensuremath{\mathsf{ND}\text{-}k\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}\bigl(V,\{(s_1,t_1),\ldots,(s_r,t_r)\}\bigr)$. It is clear that $Z$ is feasible. Mimicking the proof of Lemma \[indnlem\], and using Corollary \[ndncor3\] in place of Lemma \[cor2\] in the proof, one can easily show that if $d$ is the depth of the recursion tree rooted $(U,{\ensuremath{\mathcal T}})$ and $Z_U={\ensuremath{\mathsf{ND}\text{-}k\text{-}}{\ensuremath{\mathsf{NMCAlg}}\xspace}}(U,{\ensuremath{\mathcal T}})$, then $$c(Z_U)\leq
\frac{2{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\Bigl(\beta|{\ensuremath{\mathcal T}}|+{\ensuremath{\mathcal V}}^x(U)\Bigr)\ln\Bigl(\frac{e^dr{\ensuremath{\mathcal V}}^x(U)}{{\ensuremath{\mathit{OPT}}}}\Bigr)\ln\ln\bigl(e(r+1)\bigr).$$ Since the depth of the recursion tree is $O(\log_2 r)$, we obtain that $c(Z)\leq O\bigl(\frac{{\ensuremath{\gamma}}}{(\sqrt{{\ensuremath{\gamma}}}-1)^2}\ln r\ln\ln r\bigr)\cdot{\ensuremath{\mathit{OPT}}}$.

[10]{}

S. Arora and C. Lund. Hardness of approximations. In D. Hochbaum, editor, [*Approximation Algorithms for NP-hard Problems*]{}, PWS Publishing co., Boston, MA, USA, 1997.

R. Ahuja, T. Magnanti, J. Orlin, Network Flows: Theory, Algorithms, and Applications. 1993.

B. Applebaum, Pseudorandom generators with long stretch and low locality from random local one-way functions. In [*Electronic Colloquium on Computational Complexity*]{} (ECCC), 18:7, 2011.

S. Arora, L. Babai, J. Stern, and Z. Swedyk, The hardness of approximate optima in lattices, codes, and systems of linear equations. , 54:317Ð331, 1997

S. Barman and S. Chawla, Region growing for multi-route cuts. In [*Proceedings of the 21st SODA*]{}, pages 404–418, 2010.

S. Bhatt, and F. T. Leighton, A framework for solving VLSI graph layout problems. 28.2 300-343, 1984.

H. Bruhn, J. Cerny, A. Hall, P. Kolman, and J. Sgall. Single source multiroute flows and cuts on uniform capacity networks. , 4(1):1–20, 2008.

D. Chakrabarty, R. Krishnaswamy, S. Li, and S. Narayanan, Capacitated Network Design on Undirected Graphs. In [*Proceedings of the 16th APPROX*]{}, pages 71–80, 2013.

S. Chawla, A. Gupta, and H. Räcke. Embeddings of negative-type metrics and an improved approximation to generalized sparsest cut. , 4(2), 2008.

C. Chekuri and S. Khanna, Algorithms for 2-Route Cut Problems. , 472-484, 2008.

J. Chuzhoy, Y. Makarychev, A. Viajayaraghavan, and Y. Zhou, Approximation algorithms and hardness of the $k$-route cut problem. In [*Proceedings of the 23rd SODA*]{}, pages 780–799, 2012.

D. Dahlhaus, D. S. Johnson, C. H. Papadimitriou, P. D. Seymour, and M. Yannakakis, The complexity of multiway cuts. , 241-251, 1992

M. Dinitz and A. Gupta. Packing interdiction and partial covering problems. In [*Proceedings of the 16th IPCO*]{}, pages 157–168, 2013.

M. Dinitz, G. Kortsarz, and R. Raz. Label cover instances with large girth and the hardness of approximating basic k-spanner, 290-301, 2012.

G. Even, J. Naor, S. Rao, and B. Schieber, Divide-and-conquer approximation algorithms via spreading metrics. , 47(4):585–616, 2000.

G. Even, J. Naor, B. Schieber, and M. Sudan, Approximating minimum feedback sets and multicuts in directed graphs. , 20:151–174, 1998.

S. Fiorini, G. Joret, U. Pietropaoli, Hitting Diamonds and Growing Cacti. , pages 191-204, 2010.

F. Fomin, D. Lokshtanov, N. Misra, and S. Saurabh, Planar $F$-Deletion: Approximation, Kernelization and Optimal FPT Algorithms. In [*Proceedings of the 53rd FOCS*]{}, pages 470–479, 2012.

N. Garg, V. Vazirani, and M. Yannakakis, Approximate max-flow min-(multi)-cut theorems and their applications. , 25:235–251, 1996.

A. Hayrapetyan, D. Kempe, M. P‡l, and Z. Svitkina, Unbalanced graph cuts. , pages 191–202, 2005.

W. Kishimoto, A method for obtaining the maximum multi-route flow in a network. , 27(4):279-291, 1996.

P. Kolman, and C. Scheideler, Approximate Duality of Multicommodity Multiroute Flows and Cuts: Single Source Case. In [*Proceedings of the 23rd SODA*]{}, pages 800–810, 2012.

P. Kolman, and C. Scheideler, Towards Duality of Multicommodity Multiroute Cuts and Flows: Multilevel Ball-Growing. In [*Proceedings of the 28th STACS*]{}, pages 129–140, 2011.

G. Kortsarz, On the Hardness of Approximating Spanners. , 30, 1999.

B. Laekhanukit, , October 2013.

B. Laekhanukit, Parameters of two-prover-one-round game and the hardness of connectivity problems. In [*Proceedings of the 25th SODA*]{}, pages 1626–1643, 2014.

F. Leighton and S. Rao., Multicommodity max-flow min-cut theorems and their use in designing approximation algorithms. , 46:787-832, 1999.

A. Li, and P. Zhang, Unbalanced graph partitioning. , 53:454-466, 2013.

J. Shi, and J. Malik, Normalized cuts and image segmentation. , 888–905, 2000.

C. A. Phillips. The network inhibition problem. In [*Proceedings of the 25th STOC*]{}, pages 776–785, 1993.

C. Papadimitriou and M. Yannakakis. Optimization, approximation, and complexity classes. , 43:425–440, 1991.

H. Räcke. Minimizing congestion in general networks. In [*Proceedings of the 43rd FOCS*]{}, pages 43–52, 2002.

P. Seymour. Packing directed circuits fractionally. , 15:281–288, 1995.

D. Shmoys Cut problems and their application to divide-and-conquer. In D. Hochbaum, editor, [*Approximation Algorithms for NP-hard Problems*]{}, PWS Publishing co., Boston, MA, USA, 1997.

C. Tovey, A simplified NP-complete satisfiability problem. 8.1: 85-89, 1984.

V. Vazirani, . Springer Verlag, 2001.

D. Williamson and D. Shmoys. , Cambridge University Press, 2010.

R. K. Wood. Deterministic network interdiction. , 17(2):1–18, 1993.

R. Zenklusen. Network flow interdiction on planar graphs. , 158:1441–1455, 2010.

The $k$-route all-pairs cut problem {#append-allpairs}
===================================

\[threeroute\] The $3$-route all-pairs cut problem is APX-hard.

We give an $L$-reduction from vertex cover on bounded-degree graphs, which is known to be [[*APX*]{}]{}-hard [@PY91]. Given a vertex-cover instance $\hat G=(\hat V, \hat E)$, where $\hat G$ has maximum degree $\alpha=O(1)$, we construct an instance $G=(V,E)$ of the $3$-route all-pairs cut problem. In the following, to avoid confusion, we will refer to the elements $(\hat V, \hat E)$ of the vertex-cover instance $\hat G$ as *vertices* and *edges*, and to the elements $(V,E)$ of the constructed $3$-route all-pairs-cut instance as *nodes* and *links*.

Let the vertices in $\hat V$ be numbered $1,2, \dots, |\hat V|$. For every vertex $v \in \hat V$, we introduce a path $P_v$ in $G$ that contains as many links as the degree of $v$. That is, $P_v$ has one link $f_{e^v}$ for every edge $e \in \hat E$ incident to $v$ in $\hat G$. We give infinite cost to such links. Let $a_v$ be the first node of the path $P_v$ and $b_v$ be the last. We add a link $(a_v,b_v)$ of unit cost in $G$. Note that $P_v$ and $(a_v,b_v)$ yields a cycle in $G$ for every $v \in \hat V$. We also connect $a_v$ to $a_{v+1}$ through a cycle formed by 3 links with infinite cost. That is, we introduce $|\hat V| - 1$ triangles connecting all paths. For each edge $e =(u,v)\in \hat E$ we introduce a node ${\ensuremath{{\ensuremath{\sigma}}}}_e$ and we connect ${\ensuremath{{\ensuremath{\sigma}}}}_e$ to the endpoints of $f_{e^u}$ and to the endpoints of $f_{e^v}$, with links of unit cost. We let $G=(V, E)$ be the resulting graph for our $3$-route all-pairs cut instance (see Fig.\[fig:cactus\]).

Let $p^*$ and $c^*$ be the cost of an optimal solution for the vertex-cover instance and the cost of an optimal solution for the $3$-route all-pairs cut instance, respectively. We claim that:

[()]{}[ =0ex =0ex ]{}

If there exists a vertex cover in $\hat G$ of size $p$, then there is a solution for the $3$-route all-pairs cut instance of cost at most $2|\hat E| +p$. Note that this implies that $c^* \leq 2|\hat E| +p^*  \leq (2\alpha + 1) p^*$.

For any feasible solution for the $3$-route all-pairs cut instance of cost at most $2|\hat E| +p$ we can construct a cover of $\hat G$ of size at most $p$.

This implies that we have an $L$-reduction, and shows that if we have a $\beta$-approximation for $3$-route all-pairs cut, then we can obtain a vertex-cover solution of size at most $\beta c^*-2|\hat E|=O(\beta)p^*$, yielding the theorem.

In proving this, a useful observation is that a set $F$ of edges is feasible for $3$-route all-pairs cut problem iff the remainder graph ${\ensuremath{{\ensuremath{\overline{G}}}}}=(V,E{\ensuremath{\setminus}}F)$ has the property that every two simple cycles meet at most at one vertex. Such a graph is called a [*cactus*]{} graph.

For (i), suppose there exists a vertex cover of size $p$. For every $v$ in the cover, we add $(a_v,b_v)$ in $F$. Furthermore, for each edge $e=(u,w) \in \hat E$, we select one vertex between $u$ and $w$ that is in the cover (at least one of them is in the cover by definition), say $u$, and we add to $F$ the links connecting ${\ensuremath{{\ensuremath{\sigma}}}}_e$ to the endpoints of $f_{e^w}$ for the other vertex $w$. It is not difficult to check that $F$ yields a feasible solution for the $3$-route all-pairs cut instance (using the relationship to cactus graphs) of the claimed cost.

For (ii), suppose we have a feasible solution $F$ for the $3$-route all-pairs cut instance, and consider the remainder graph obtained by removing $F$. Clearly, all links of infinite cost are still present. Note that each node ${\ensuremath{{\ensuremath{\sigma}}}}_e$ can have at most two links incident to it in the remainder graph, and both these links must be incident to two nodes of the same path $P_v$ for some $v$. If not, then we would have an infinite-cost link of some triangle that connects the vertices $\{a_v\}_{v \in \hat V}$ that is contained into another cycle other than the triangle, contradicting feasibility of our solution. We first argue that we may assume that each ${\ensuremath{{\ensuremath{\sigma}}}}_e$ has exactly two links incident to it in the remainder graph (and hence, also in $F$). Let $e=(u,w)$ be the edge in $\hat E$ corresponding to ${\ensuremath{{\ensuremath{\sigma}}}}_e$. As argued above, $F$ contains at least one pair of links that connect ${\ensuremath{{\ensuremath{\sigma}}}}_e$ to the nodes of a path, say $P_u$. Suppose that $F$ also contains some links connecting ${\ensuremath{{\ensuremath{\sigma}}}}_e$ to $P_w$. If we remove such links from $F$, we create one additional cycle, containing the link $f_{e^w}$, in the remainder graph. Thus, the new remainder graph is not a cactus iff $f_{e^w}$ is already contained in some cycle in $(V,E\setminus F)$. But there is only one possible cycle in $(V,E\setminus F)$ containing $f_{e^w}$, namely the cycle formed by $P_w$ and $(a_w,b_w)$. This observation implies that $F\cup\{(a_w,b_w)\}\setminus\{\text{the two links connecting ${\ensuremath{{\ensuremath{\sigma}}}}_e$ to $P_w$}\}$ is a feasible solution to our $3$-route all-pairs cut instance. Furthermore, this solution has no greater cost since we are adding at most one link of unit cost, and we removing at least one link of the same cost.

Since the cost of our solution is at most $2|\hat E| + p$, it follows that there are at most $p$ links in $F$ of the form $(a_v, b_v)$. We claim that these vertices $v$ form a cover in $\hat G$. Suppose not. Then there is at least one edge $e \in \hat E$ that is not covered by these vertices. We know that the node ${\ensuremath{{\ensuremath{\sigma}}}}_e$ is connected in ${\ensuremath{{\ensuremath{\overline{G}}}}}$ to the endpoints of the link $f_{e^u}$ for one of the endpoints, say $u$, of $e$. The link $f_{e^u}$ is therefore contained in the cycle formed by $P_u$ and $(a_u,b_u)$, since $(a_u,b_u)$ is not in $F$, and is also contained in the triangle with the node ${\ensuremath{{\ensuremath{\sigma}}}}_e$, which contradicts feasibility of $F$.

\[allk\] The $k$-route all-pairs cut problem is APX-hard for all $k\geq 3$.

The reduction is very similar to the one in the proof of Theorem \[threeroute\]. The only change is that in the graph $G$ created from the given vertex-cover instance, we now have: (a) $k-2$ parallel links between every pair of consecutive nodes of every path $P_v$; and (b) $k-2$ parallel links between $a_j, a_{j+1}$ for all $j=1,\ldots,|\hat V| -1$.

Suppose there exists a vertex cover of size $p$. As before, for each node ${\ensuremath{{\ensuremath{\sigma}}}}_e$ we remove exactly one pair of links incident to ${\ensuremath{{\ensuremath{\sigma}}}}_e$, and in particular we choose the pair of links that connect ${\ensuremath{{\ensuremath{\sigma}}}}_e$ to the endpoints of the edge $f_{e^u}$ if $e=(u,v)$ and $v$ is in the cover, and the other pair otherwise. We also remove all edges of the form $(a_v,b_v)$ for $v$ in the cover. We remove $2|\hat E|+p$ edges in total.

We claim that for every pair of nodes $z,w$ of the remainder graph, we have at most $k-1$ edge-disjoint paths. Every node that is not a node of a path $P_v$ has maximum degree two, and therefore this is clear. If $z$ and $w$ belongs to two different paths $P_v$ and $P_u$ with $u>v$, then every path connecting them must use either the link $(a_v, a_{v +1})$ (and there are at most $k-2$ such links) or the two upper links of the triangle formed with $a_v$ and $a_{v+1}$. Therefore, we can have at most $k-1$ edge-disjoint such paths. Finally, if $z$ and $w$ belong to the same path $P_u$, we have $k-2$ paths given by the infinite cost links, plus at most one additional path that uses either the edge $(a_u,b_u)$ or a sequence of pairs of links incident into the nodes $\{{\ensuremath{{\ensuremath{\sigma}}}}_e\}$ for the edges $e$ that have $u$ as an endpoint in $\hat G$. Note that, by construction, if $(a_u,b_u)$ is still in the graph, then all the pairs of links incident into the nodes $\{{\ensuremath{{\ensuremath{\sigma}}}}_e\}$ for the edges $e$ that have $u$ as an endpoint have been removed, and therefore, once again we get at most $k-1$ edge-disjoint paths.

For the other direction, suppose we have a feasible solution $F$ for the $k$-route all-pairs cut instance, and consider the remainder graph obtained by removing such set of links. Clearly, all links of infinite cost are still present. Once again, each node ${\ensuremath{{\ensuremath{\sigma}}}}_e$ corresponding to an edge $e=(u,v)$ can have at most two links incident to it in the remainder graph, and both these links must be incident to two nodes of the same path. If not, then we would have a pair of nodes ($a_v$ and $a_u$) that are connected by $k$ edge-disjoint paths: $k-1$ given by the infinite-cost links not in the paths $P_u$ and $P_v$, and one which uses the links in the paths $P_u$ and $P_v$ and two links incident to ${\ensuremath{{\ensuremath{\sigma}}}}_e$. Also, as before, we may assume that each ${\ensuremath{{\ensuremath{\sigma}}}}_e$ has exactly two links incident into it in the remainder graph (and hence, in $F$), because otherwise for one endpoint $u$ of $e$ we get that $F \cup\{(a_u,b_u)\} \setminus\{\text{the two links connecting ${\ensuremath{{\ensuremath{\sigma}}}}_e$ to $P_u$}\}$ is a feasible solution to our $k$-route all-pairs cut instance of no larger cost. So if the cost of $F$ is at most $2|\hat E| + p$, it follows that we have at most $p$ links in $F$ of the form $(a_v, b_v)$. We claim that these vertices $v$ form a cover in $\hat G$.

Suppose not. Then there is at least one edge $e  \in \hat E$ that is not covered by such vertices. We know that the node ${\ensuremath{{\ensuremath{\sigma}}}}_e$ is connected in ${\ensuremath{{\ensuremath{\overline{G}}}}}$ to the endpoints of the link $f_{e^u}$ for one of the endpoints, say $u$, of $e$. The endpoints of $f_{e^u}$ are therefore connected by $k-2$ parallel paths using one single link, one path formed by the edge $(a_u,b_u)$ and edges of $P_u \setminus\{f_{e^u}\}$, and one path contained in the triangle with the node ${\ensuremath{{\ensuremath{\sigma}}}}_e$. This contradicts feasibility of $F$.

On the positive side, the $3$-route all-pairs cut problem admits an $O(1)$-approximation. This follows from: (1) the equivalence of $3$-route all-pairs cut and the problem of removing a min-cost set of edges so that the remainder graph is a cactus; (2) the results of Fiorini et al. [@FioriniJP10], who gave an $O(1)$-approximation for the problem of removing a minimum-weight node set so that the remaining graph is a cactus; and (3) the edge-removal version easily reduces to the node-removal version by subdividing edges, and setting the cost of the original vertices to $\infty$ and the cost of each vertex corresponding to an edge to be the cost of the edge.

Recently, Fomin et al. [@FominLMS12] developed an $O(1)$-approximation algorithm for the problem of removing the fewest number of nodes so that the remaining graph excludes a minor from a given list of graphs, at least one of which should be planar. While $k$-route all-pairs cut can be stated as excluding the planar graph with $k$ parallel edges as a minor of the remainder graph, the result of [@FominLMS12] does not directly apply here. This is because our transformation of an edge-weighted instance to a node weighted one introduces non-uniform node weights, whereas the algorithm in [@FominLMS12] is for uniform node weights.

Hardness of the edge-deletion $k$-route node-multicut problem {#knmc}
=============================================================

Recall that in the edge-deletion $k$-route node-multicut ([$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}) problem, we have an undirected graph $G=(V,E)$ with nonnegative edge costs $\{c_e\}_{e\in E}$, and $r$ source-sink pairs $(s_1,t_1),\dots, (s_r,t_r)$ and an integer $k\geq 1$. We seek a minimum-cost set $F{\subseteq}E$ of edges such that the remainder graph ${\ensuremath{{\ensuremath{\overline{G}}}}}=(V,E{\ensuremath{\setminus}}F)$ contains at most $k-1$ node-disjoint $s_i$-$t_i$ paths for all $i=1,\ldots,r$.

Chuzhoy et al. [@ChuzhoyMVZ12] show that [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{} is hard to approximate within a factor $\Omega(k^\epsilon)$. They present a reduction from $3$-SAT$(5)$, which is the variant of $3$-SAT where each variable occurs in at most $5$ clauses, coupled with the parallel-repetition theorem, which is essentially a reduction from (the minimization version) of [*label cover*]{}. However, Laekhanukit [@Laekhanukit14] pointed out some subtle (but fixable) errors in their proof and proposed a correction, but his reduction also suffers from some subtle (again fixable) errors [@Laekhanukit13]. We give a correct proof below via a somewhat simpler reduction than the ones in [@ChuzhoyMVZ12; @Laekhanukit14].

Label cover was first introduced by Arora et al. [@AroraBSS97] and has been subsequently used as a basis for many hardness reductions (see, e.g., [@AroraL97]). Kortsarz [@Kortsarz99] presented a minimization version of label cover (sometimes known as [$\mathsf{MinRep}$]{}) with the same hardness guarantee, that has since found use in various network-design applications (see, e.g., [@DinitzKR12]).

In the [$\mathsf{MinRep}$]{}problem, we are given a bipartite graph $H=(U \cup W, F)$, two sets of labels $L_1$ (for vertices in $U$) and $L_2$ (for vertices in $W$), and a constraint function for each edge $e$ defined as $\pi_e: L_1 \rightarrow L_2$. A *labeling* is given by specifying a set of labels $f(u) \subseteq L_1$ for every vertex $u \in U$ and a set of labels $f(w) \subseteq L_2$ for every vertex $w \in W$. We say that a labeling *covers* an edge $e=uw \in F$ if there exists $a \in f(u)$ and $b \in f(w)$ such that $\pi_e(a) = b$. *Min-Rep* asks for a labeling that covers all the edges while minimizing $\sum_{u \in U} |f(u)| + \sum_{w \in W} |f(w)|$.

\[lchard\] There are constants ${\ensuremath{\epsilon}}_0,{\ensuremath{\delta}}_0>0$ such that there is no polytime algorithm for [$\mathsf{MinRep}$]{}with approximation factor:

[–]{}[=0.5ex =0ex ]{}

$O\bigl(q^{{\ensuremath{\epsilon}}_0}\bigr)$ unless P=NP, where $q=|L_1|+|L_2|$ is the size of the label set;

$O\bigl({\ensuremath{\Delta}}^{{\ensuremath{\delta}}_0}\bigr)$ unless P=NP, where ${\ensuremath{\Delta}}$ is the maximum degree of the underlying graph;

$2^{\log^{1-\epsilon} m}$ for any constant ${\ensuremath{\epsilon}}$, unless NP is contained in deterministic quasipolynomial time, where $m$ is the number of edges.

\[edknmc-hard\] There is a polytime approximation-preserving reduction that given a [$\mathsf{MinRep}$]{}-instance $(H,\pi,L_1,L_2)$ with label-size $q=|L_1|+|L_2|$ and maximum-degree ${\ensuremath{\Delta}}$, constructs an [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}-instance $\bigl(G,\{c_e\},\{s_1,t_1,\ldots,s_r,t_r\},k\bigr)$ with $k=O({\ensuremath{\Delta}}q)$, $r=|E_H|$, and $|E_G|=O(|E_H|q)$.\
Hence, there is no $O\bigl(k^{{\ensuremath{\epsilon}}_0}\bigr)$-approximation for [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}, for some constant ${\ensuremath{\epsilon}}_0>0$, unless P=NP, and no $2^{\log^{1-{\ensuremath{\epsilon}}}|E_G|}$-approximation for any ${\ensuremath{\epsilon}}>0$ unless NP is contained in deterministic quasipolynomial time.

We first describe the construction and then argue the approximation-preservation property.

#### The construction.

For each vertex $u \in U$ and for each label in $a \in L_1$ we introduce two vertices $a^u, \bar a^u$ in $G$ connected by an edge of unit cost. Intuitively, if we select this edge in an [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}solution for the instance we construct, this implies that we are selecting label $a$ for $u$. Similarly, for each $w \in W$ and for each label $b \in L_2$ we introduce two vertices $b^w, \bar b^w$ connected by an edge of unit cost. The above edges will be the only ones having unit cost. All subsequent edges added to this construction will have infinite cost.

Consider an edge $e=(u,w)$ of $H$. For each such edge, we construct the following gadget. We introduce two nodes $s_e,t_e$ that will form a terminal pair in our new instance. $s_e$ and $t_e$ are connected as follows. For each $b \in L_2$, let $L^e_b \subseteq L_1$ be the labels of $L_1$ such that $a \in L^e_b$ implies $\pi_e(a) = b$. Clearly, the sets $L^e_b$, for all $b \in L_2$, form a partition of the labels $L_1$. For each non empty set $L^e_b$ we add a path $P^e_b$ of length 2 with the middle vertex connected to $t_e$. We then add edges to form a cycle $C^e_b$ starting and ending at $s_e$, containing all the edges in $L^e_b$, the edges in the path $P^e_b$ and the edge $b^w \bar b^w$. Finally we split each of these added edges into 2 by introducing a middle vertex. We let $V(C^e_b)$ we the the set of new vertices introduced by this operation, and let $N_e$ be all the vertices participating in this gadget; see Fig. \[fig:1\].

Let $G'=\bigl(\bigcup_{e'\in E_H}N_{e'},\bigcup_{e'\in E_H}E(N_{e'})\bigr)$ be the graph formed by the vertices and edges of all the edge gadgets. For all $e$, we are going to add other edges forming paths (of length 2) between $s_e$ and $t_e$. We add edges $(s_e,v)$, $(v,t_e)$ for all vertices $v\in{\ensuremath{\Gamma}}_{G'}(N_e)$, that is, for all $v\notin N_e$ that are adjacent in some edge-gadget to some node in $N_e$. Note that $N_e\cap N_{e'}={\ensuremath{\emptyset}}$ unless $e$ and $e'$ share an endpoint in $H$, say $u$, in which case, the two gadgets share the vertices $\{a^u, \bar a^u\}$ for all labels $a$ of $u$. Thus, $|{\ensuremath{\Gamma}}_{G'}(N_e)|=O(\Delta q)$.

Define $k_e:=|\{b\in L_2: L^e_b\neq{\ensuremath{\emptyset}}\}|+|{\ensuremath{\Gamma}}_{G'}(N_e)|$. Finally, set $k:= \max_{e} k_e$. For all edges $e$ with $k_e<k$, we add $k-k_e$ new vertices and connect these to $s_e$ and $t_e$ (via $\infty$-cost edges). Let $G$ be the resulting graph. This concludes our construction.

#### Approximation preservation.

We now argue that any feasible solution to the [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}-instance of finite cost yields a feasible solution to the [$\mathsf{MinRep}$]{}-instance of no greater cost, and vice versa. This will complete the proof.

($\Rightarrow$) Let $Z$ be a solution of finite cost for our [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}instance. Consider a node $u\in V_H$ and let $L_u\in\{L_1,L_2\}$ be the label-set of $u$. Set $f(u):= \{a \in L_u: (a^u, \bar a^u)\in Z\}$. Clearly, the cost of the two solutions are the same. We now claim that the resulting labeling is feasible for the label-cover instance. Suppose not, then there is an edge $e\in E_H$ that is not covered. By our construction, this means that for each $b \in L_2$ with $L^e_b\neq{\ensuremath{\emptyset}}$, the subgraph of the remainder subgraph ${\ensuremath{{\ensuremath{\overline{G}}}}}=(V_G,E_G\setminus Z)$ induced by the nodes of the cycle $C^e_b$ and $t_e$ is connected. Each such cycle $C^e_b$, yields therefore one vertex-disjoint path in ${\ensuremath{{\ensuremath{\overline{G}}}}}$ between $s_e$ and $t_e$. Also, all edges incident to $s_e$ and $t_e$ are present in ${\ensuremath{{\ensuremath{\overline{G}}}}}$ (since they have $\infty$ cost), so all length-2 paths in $G$ between $s_e$ and $t_e$ are still present in ${\ensuremath{{\ensuremath{\overline{G}}}}}$. It follows that the vertex connectivity of $s_e$ and $t_e$ is at least $k$, a contradiction.

($\Leftarrow$) For the other direction, given a labeling for the label-cover instance, we construct $Z:=\{(a^u,\bar a^u): u\in V_H,\ a \in f(u)\}$. Clearly, the cost of the two solutions is the same. We claim that $Z$ is a feasible solution to our [$\mathsf{ED}$-[$k$-[$\mathsf{NMC}$]{}]{}]{}instance. Suppose not. Then, for some $e=(u,w)\in E_H$, the $s_e$-$t_e$ vertex connectivity in the remainder graph ${\ensuremath{{\ensuremath{\overline{G}}}}}=(V_G,E_G{\ensuremath{\setminus}}Z)$ is at least $k$. Therefore we can find a set of vertex-disjoint paths $\mathcal P$ between these vertices of size $|\mathcal P|\geq k$. Without loss of generality, we may assume that all the $k-|\{b\in L_2:L^e_b\neq{\ensuremath{\emptyset}}\}|$ length-2 paths between $s_e$, $t_e$ are in $\mathcal P$. If we remove the internal nodes on these length-2 paths from ${\ensuremath{{\ensuremath{\overline{G}}}}}$, the connected component containing $s_e$ in the remaining portion of ${\ensuremath{{\ensuremath{\overline{G}}}}}$ is a subgraph of the gadget $(N_e,E(N_e))$ for edge $e$ (shown in Fig. \[fig:1\]). This means that this subgraph contains at least $|\{b\in L_2:L^e_b\neq{\ensuremath{\emptyset}}\}|$ $s_e$-$t_e$ vertex-disjoint paths. Clearly, this is only possible if, for every label $b$ such that $L^e_b\neq{\ensuremath{\emptyset}}$, either $(b^w,\bar b^w)\notin Z$ or $(a^u,\bar a^u)\notin Z$ for every $a\in L^e_b$. But that means that the edge $e$ is not covered by our labeling, a contradiction.

[^1]: [ggurugan@cs.cmu.edu]{}. Department of Computer Science, Carnegie Mellon University, Pittsburgh, USA. Work done while the author was an undergraduate research assistant at the University of Waterloo under Prof. Sanita.

[^2]: [sanita@uwaterloo.ca]{}. Dept. of Combinatorics and Optimization, Univ. Waterloo, Waterloo, ON N2L 3G1.

[^3]: [cswamy@math.uwaterloo.ca]{}. Dept. of Combinatorics and Optimization, Univ. Waterloo, Waterloo, ON N2L 3G1. Supported in part by NSERC grant 327620-09, an NSERC Discovery Accelerator Supplement Award and an Ontario Early Researcher Award.
