arxiv: v2 [cs.ds] 27 Aug 2016

Transcription

1 ariv: v2 [cs.ds] 27 Aug 2016 Recoverable robust spanning tree problem under interval uncertainty representations Mikita Hradovich, Adam Kasperski, Pawe l Zieliński Faculty of Fundamental Problems of Technology, Wroc law University of Technology, Wroc law, Poland Faculty of Computer Science and Management, Wroc law University of Technology, Wroc law, Poland {mikita.hradovich,adam.kasperski,pawel.zielinski}@pwr.edu.pl Abstract This paper deals with the recoverable robust spanning tree problem under interval uncertainty representations. A strongly polynomial time, combinatorial algorithm for the recoverable spanning tree problem is first constructed. This problem generalizes the incremental spanning tree problem, previously discussed in literature. The algorithm built is then applied to solve the recoverable robust spanning tree problem, under the traditional interval uncertainty representation, in polynomial time. Moreover, the algorithm allows to obtain several approximation results for the recoverable robust spanning tree problem under the Bertsimas and Sim interval uncertainty representation and the interval uncertainty representation with a budget constraint. Keywords: robust optimization; interval data; recovery; spanning tree 1 Introduction Let G = (V,E), V = n, E = m, be an undirected graph and let Φ be the set of all spanning trees of G. In the minimum spanning tree problem, a cost is specified for each edge, and we seek a spanning tree in G of the minimum total cost. This problem is well known and can be solved efficiently by using several polynomial time algorithms (see, e.g., [1, 18]). In this paper, we first study the recoverable spanning tree problem (Rec ST for short). Namely, for each edge e E, we are given a first stage cost C e and a second stage cost c e (recovery stage cost). Given a spanning tree Φ, let Φ k be the set of all spanning trees Φ such that \ k (the recovery set), where k is a fixed integer in [0,n 1], called the recovery parameter. Note that Φ k can be seen as a neighborhood of containing all spanning trees which can be obtained from by exchanging up to k edges. The Rec ST problem can be stated formally as follows: Rec ST : min Φ ( e C e + min Φ k ) c e. (1) e We thus seek a first stage spanning tree Φ and a second stage spanning tree Φ k, so that the total cost of and for C e and c e, respectively, is minimum. Notice that Rec ST 1

2 generalizes the following incremental spanning tree problem, investigated in [8]: Inc ST : min Φ k c e, (2) where Φ is a given spanning tree. So, we wish to find an improved spanning tree with the minimum cost, within a neighborhood of determined by Φ k. Several interesting practical applications of the incremental network optimization were presented in [8]. It is worth pointing out that Inc ST can be seen as the Rec ST problem with a fixed first stage spanning tree, whereas in Rec ST both the first and the second stage trees are unknown. It has been shown in [8] that Inc ST can be solved in strongly polynomial time by applying the Lagrangian relaxation technique. On the other hand, no strongly polynomial time combinatorial algorithm for Rec ST has been known to date. Thus proposing such an algorithm for this problem is one of the main results of this paper. The Rec ST problem, beside being an interesting problem per se, has an important connection with a more general problem. Namely, it is an inner problem in the recoverable robust model with uncertain recovery costs, discussed in [4, 5, 6, 7, 15, 17]. Indeed, the recoverable spanning tree problem can be generalized by considering its robust version. Suppose that the second stage costs c e, e E, are uncertain and let U contain all possible realizations of the second stage costs, called scenarios. We will denote by c S e the second stage cost of edge e E under scenario S U, where S = (c S e) e E is a cost vector. In the recoverable robust spanning tree problem (Rob Rec ST for short), we choose an initial spanning tree in the first stage, with the cost equal to e C e. Then, after scenario S U reveals, can be modified by exchanging at most k edges, obtaining a new spanning tree Φ k. The second stage cost of under scenario S U is equal to e cs e. Our goal is to find a pair of trees and such that \ k, which minimizes the sum of the first and the second stage costs e C e + e cs e in the worst case. The Rob Rec ST problem is defined formally as follows: ( ) Rob Rec ST : min C e +max min c S Φ S U Φ k e. (3) e If C e = 0 for each e E and k = 0, then Rob Rec ST is equivalent to the following min-max spanning tree problem, examined in [2, 13, 12], in which we seek a spanning tree that minimizes the largest cost over all scenarios: Min-Max ST : min max c S e. (4) Φ S U If C e = 0 for each e E and k = n 1, then Rob Rec ST becomes the following adversarial problem [17] in which an adversary wants to find a scenario which leads to the greatest increase in the cost of the minimum spanning tree: Adv ST : max min c S e. (5) S U Φ We now briefly recall the known complexity results on Rob Rec ST. It turns out that its computational complexity highly relies on the way of defining the scenario set U. There are two popular methods of representing U, namely the discrete and interval uncertainty representations. For the discrete uncertainty representation (see, e.g., [13]), scenario set, e e e e 2

3 denoted by U D, contains K explicitly listed scenarios, i.e. U D = {S 1,S 2,...,S K }. In this case, the Rob Rec ST problem is known to be NP-hard for K = 2 and any constant k [11]. Furthermore, it becomes strongly NP-hard and not at all approximable when both K and k are a part of the input [11]. It is worthwhile to mention that Min-Max ST is NP hard even when K = 2 and becomes strongly NP-hard and not approximable within O(log 1 ǫ n) for any ǫ > 0 unless NP DTIME(n polylogn ), when K is a part of input [13, 12]. It admits an FPTAS, when K is a constant [2] and is approximable within O(log 2 n), when K is a part of the input [12]. The Adv ST problem, under scenario set U D, is polynomially solvable, since it boils down to solving K traditional minimum spanning tree problems. For the interval uncertainty representation, which is considered in this paper, one assumes that the second stage cost of each edge e E is known to belong to the closed interval [c e,c e +d e ], where c e is a nominal cost of e E and d e 0 is the maximum deviation of the cost of e from its nominal value. In the traditional case U, denoted by U I, is the Cartesian product of all these intervals [13], i.e. U I = {S = (c S e) e E : c S e [c e,c e +d e ],e E}. (6) In [4] a polynomial algorithm for the recoverable robust matroid basis problem under scenario set U I was constructed, provided that the recovery parameter k is constant. In consequence, Rob Rec ST under U I is also polynomially solvable for constant k. Unfortunately, the algorithm proposed in [4] is exponential in k. Interestingly, the corresponding recoverable robust version of the shortest path problem (Φ is replaced with the set of all s t paths in G) has been proven to be strongly NP-hard and not at all approximable even if k = 2 [5]. It has been recently shown in [10] that Rob Rec ST under U I is polynomially solvable when k is a part of the input. In order to prove this result, a technique called the iterative relaxation of a linear programming formulation, whose framework was described in [14], has been applied This technique, however, does not imply directly a strongly polynomial algorithm for Rob Rec ST, since it requires the solution of a linear program. In [3] apopularand commonly used modification of thescenario set U I hasbeen proposed. The new scenario set, denoted as U1 I(Γ), is a subset of UI such that under each scenario in U1 I(Γ), the costs of at most Γ edges are greater than their nominal values c e, where Γ is assumed to be a fixed integer in [0,m]. Scenario set U1 I (Γ) is formally defined as follows: U I 1 (Γ) = {S = (cs e ) e E : c S e [c e,c e +δ e d e ],δ e {0,1},e E, e Eδ e Γ}. (7) The parameter Γ allows us to model the degree of uncertainty. When Γ = 0, then we get Rec ST (Rob Rec ST with one scenario S = (c e ) e E ). On the other hand, when Γ = m, then we get Rob Rec ST under the traditional interval uncertainty U I. It turns out that the Adv ST problem under U1 I (Γ) is strongly NP-hard (it is equivalent to the problem of finding Γ most vital edges) [17, 16, 9]. Consequently, the more general Rob Rec ST problem is also strongly NP-hard. Interestingly, the corresponding Min-Max ST problem with U1 I(Γ) is polynomially solvable [3]. et another interesting way of defining scenario set, which allows us to control the amount of uncertainty, is called the scenario set with a budget constraint (see, e.g,. [17]). This scenario set, denoted as U2 I (Γ), is defined as follows: U I 2 (Γ) = {S = (cs e ) e E : c S e = c e +δ e,δ e [0,d e ],e E, e Eδ e Γ}, (8) 3

4 whereγ 0isaafixedparameterthat canbeseenasabudgetofanadversary, andrepresents the maximum total increase of the edge costs from their nominal values. Obviously, if Γ is sufficientlylarge, thenu2 I(Γ)reducestothetraditionalinterval uncertaintyrepresentation UI. Thecomputational complexity of Rob Rec STforscenarioset U2 I isstill open. Weonly know that its special cases, namely Min-Max ST and Adv ST, are polynomially solvable [17]. In this paper we will construct a combinatorial algorithm for Rec ST with strongly polynomial running time. We will apply this algorithm for solving Rob Rec ST under scenario set U I in strongly polynomial time. Moreover, we will show how the algorithm for Rec ST can beused to obtain several approximation results for Rob Rec ST, underscenario sets U1 I(Γ) and UI 2 (Γ). This paper is organized as follows. Section 2 contains the main result of this paper a combinatorial algorithm for Rec ST with strongly polynomial running time. Section 3 discusses Rob Rec ST under the interval uncertainty representations U I, U1 I(Γ), and U2 I(Γ). 2 The recoverable spanning tree problem In this section we construct a combinatorial algorithm for Rec ST with strongly polynomial running time. Since = n 1 for each Φ, Rec ST (see (1)) is equivalent to the following mathematical programming problem: min C e + e e s.t. L,, Φ, c e (9) where L = n 1 k. Problem (9) can be expressed as the following MIP model: Opt = min C e x e + e y e (10) e E e Ec s.t. x e = n 1, (11) e E e E(U) x e U 1, U V, (12) y e = n 1, (13) e E e E(U) y e U 1, U V, (14) x e z e 0, e E, (15) y e z e 0, e E, (16) z e L, (17) e E x e,y e,z e 0, integer e E, (18) where E(U) stands for the set of edges that have both endpoints in U V. We first apply the Lagrangian relaxation (see, e.g., [1]) to (10)-(18) by relaxing the cardinality constraint (17) 4

5 with a nonnegative multiplier θ. We also relax the integrality constraints (18). We thus get the following linear program (with the corresponding dual variables which will be used later): φ(θ) = min C e x e + c e y e θ e +θl (19) e E e E e Ez s.t. x e = n 1, [µ], e E e E(U) y e = n 1, e E e E(U) x e ( U 1), U V, [w U ], [ν], y e ( U 1), U V, [v U ], x e z e 0, e E, [α e ], y e z e 0, e E, [β e ], x e,y e,z e 0, e E. For any θ 0, the Lagrangian function φ(θ) is a lower bound on Opt. It is well-known that φ(θ) is concave and piecewise linear. By the optimality test (see, e.g., [1]), we obtain the following theorem: Theorem 1. Let (x e,y e,z e ) e E be an optimal solution to (19) for some θ 0, feasible to (11)-(18) and satisfying the complementary slackness condition θ( e E z e L) = 0. Then (x e,y e,z e ) e E is optimal to (10)-(18). Let (,),, Φ, be a pair of spanning trees of G (a pair for short). This pair corresponds to a feasible 0 1 solution to (19), defined as follows: x e = 1 for e, y e = 1 for e, and z e = 1 for e ; the values of the remaining variables are set to 0. From now on, by a pair (,) we also mean a feasible solution to (19) defined as above. Given a pair (,) with the corresponding solution (x e,y e,z e ) e E, let us define the partition (E,E,E Z,E W ) of the set of the edges E in the following way: E = {e E : x e = 1,y e = 0}, E = {e E : y e = 1,x e = 0}, E Z = {e E : x e = 1,y e = 1} and E W = {e E : x e = 0,y e = 0}. Thus equalities: = E E Z, = E E Z and E Z = hold. Our goal is to establish some sufficient optimality conditions for a given pair (,) in the problem (19). The dual to (19) has the following form: φ D (θ) = max U V s.t. ( U 1)w U +(n 1)µ U V( U 1)v U +(n 1)ν +θl (20) {U V :e E(U)} {U V :e E(U)} w U +µ C e α e, e E, v U +ν c e β e, e E, α e +β e θ, e E, w U,v U 0, U V, α e,β e 0, e E. 5

6 Lemma 1. The dual problem (20) can be rewritten as follows: ( ) φ D (θ) = min (C e α e )+min (c e β e ) +θl. (21) max {α e 0,β e 0:α e+β e θ, e E} Φ e Φ e Proof. Fix some α e and β e such that α e +β e θ for each e E in (20). For these constant values of α e and β e, e E, using the dual to (20), we arrive to min Φ e (C e α e ) + min Φ e (c e β e )+θl and the lemma follows. Lemma 1 allows us to establish the following result: Theorem 2 (Sufficient pair optimality conditions). A pair (, ) is optimal to (19) for a fixed θ 0 if there exist α e 0, β e 0 such that α e +β e = θ for each e E and (i) is a minimum spanning tree for the costs C e α e, is a minimum spanning tree for the costs c e β e, (ii) α e = 0 for each e E, β e = 0 for each e E. Proof. By the primal-dual relation, the inequality φ D (θ) φ(θ) holds. Using (21), we obtain φ D (θ) (C e α e )+ (c e β e )+θl = C e + c e + (C e +c e θ)+θl e e e E e E e E Z = C e + C e + c e + c e θ E Z +θl e E e E Z e E e E Z = e C e + e c e θ E Z +θl = φ(θ). The Weak Duality Theorem implies the optimality of (,) in (19) for a fixed θ 0. Lemma 2. A pair (,), which satisfies the sufficient pair optimality conditions for θ = 0, can be computed in polynomial time. Proof. Let be a minimum spanning tree for the costs C e and be a minimum spanning tree for the costs c e, e E. Since θ = 0, we set α e = 0, β e = 0 for each e E. It is clear that (, ) satisfies the sufficient pair optimality conditions. Assume that (,) satisfies the sufficient pair optimality conditions for some θ 0. If, for this pair, E Z L and θ( E Z L) = 0, then we are done, because by Theorem 1, the pair(,)isoptimalto(10)-(18). Supposethat E Z < L((,)isnotfeasibleto(10)-(18)). We will now show a polynomial time procedure for finding a new pair (, ), which satisfies the sufficient pair optimality conditions and E Z = E Z + 1. This implies a polynomial time algorithm for the problem (10)-(18), since it is enough to start with a pair satisfying the sufficient pair optimality conditions for θ = 0 (see Lemma 2) and repeat the procedure at most L times, i.e. until E Z = L. Given a spanning tree T in G = (V,E) and edge e = {k,l} T, let us denote by P T (e) the unique path in T connecting nodes k and l. It is well known that for any f P T (e), T = T {e}\{f} is also a spanning tree in G. We will say that T is the result of a move on T. 6

7 Consider a pair (, ) that satisfies the sufficient pair optimality conditions for some fixed θ 0. Set C e = C e α e and c e = c e β e for every e E, where α e and β e, e E, are the numbers which satisfy the conditions in Theorem 2. Thus, by Theorem 2(i) and the path optimality conditions (see, e.g., [1]), we get the following conditions which must be satisfied by (,): for every e / C e C f for every f P (e), (22a) for every e / c e c f for every f P (e). (22b) We now build a so-called admissible graph G A = (V A,E A ) in two steps. We first associate with each edge e E a node v e and include it to V A, V A = E. We then add arc (v e,v f ) to E A if e /, f P (e) and Ce = Cf. This arc is called an -arc. We also add arc (v f,v e ) to E A if e /, f P (e) and c e = c f. This arc is called an -arc. We say that v e V A is admissible if e E, or v e is reachable from a node v g V A, such that g E, by a directed path in G A. In the second step we remove from G A all the nodes which are not admissible, together with their incident arcs. An example of an admissible graph is shown in Figure 1. Each node of this admissible graph is reachable from some node v g, g E. Note that the arcs (v e7,v e6 ) and (v e7,v e10 ) are not present in G A, because v e7 is not reachable from any node v g, g E. These arcs have been removed from G A in the second step. Observe that each -arc (v e,v f ) E A represents a move on, namely = {e}\{f} is a spanning tree in G. Similarly, each -arc (v e,v f ) E A represents a move on, namely = {f}\{e} is a spanning tree in G. Notice that the cost, with respect to Ce, of is the same as and the cost, with respect to c e, of is the same as. So, the moves indicated by -arcs and -arcs preserve the optimality of and, respectively. Observe that e / or e, which implies e / E. Also f or f /, which implies f / E. Hence, no arc in E A can start in a node corresponding to an edge in E and no arc in E A can end in a node corresponding to an edge in E. Observe also that (v e,v f ) E A can be both -arc and -arc only if e E and f E. Such a case is shown in Figure 1 (see the arc (v e1,v e2 )). Since each arc (v e,v f ) E A represents a move on or, e and f cannot both belong to E W or E Z. (a) e 6 (2;3) e 1 (2;2) e 3 (2;2) e10 e 7 e 5 (2;2) (2;3) (2;3) e 4 (1;5) e 8 (2;3) e (5;0) 9 e 2 (2;2) v e6 E v e2 E v e1 E (Ce;c ; e) (b) v e8 E W v e3 E Z v e5 E v e10 E Z Figure 1: (a) A pair (,) such that = {e 2,e 3,e 4,e 6,e 10 } and = {e 1,e 3,e 5,e 9,e 10 }. (b) The admissible graph G A for (,). We will consider two cases: E {e E : v e V A } ande {e E : v e V A } =. The first case means that there is a directed path from v e, e E, to a node v f, f E, 7

8 PSfrag replacements in the admissible graph G A and in the second case no such a path exists. We will show that in the first case it is possible to find a new pair (, ) which satisfies the sufficient pair optimality conditions and E Z = E Z +1. The idea will be to perform a sequence of moves on and, indicated by the arcs on some suitably chosen path from v e, e E, to v f, f E in the admissible graph G A. Let us formally handle this case. Lemma 3. If E {e E : v e V A }, then there exists a pair (, ) with E Z = E Z +1, which satisfies the sufficient pair optimality conditions for θ. Proof. We begin by introducing the notion of a cycle graph G(T) = (V T,A T ), corresponding to a given spanning tree T of graph G = (V,A). We build G(T) as follows: we associate with each edge e E a node v e and include it to V T, E = V T ; then we add arc (v e,v f ) to A T if e T and f P T (e). An example is shown in Figure 2. (a) (b) f 1 e 2 v e1 v e2 f 2 f 3 f 4 e 1 v f1 v f6 v f7 v f5 v f2 v f3 v f4 f 5 e 3 f 6 f 7 v e3 Figure 2: (a) A graph G with a spanning tree T (the solid lines). (b) The cycle graph G(T). Claim 1. Given a spanning tree T of G, let F = {(v e1,v f1 ),(v e2,v f2 ),...,(v el,v fl )} be a subset of arcs of G(T), where all v ei and v fi (resp. e i and f i ), i [l], are distinct. If T = T {e 1 PSfrag,...,ereplacements l } \ {f 1,...,f l } is not a spanning tree, then G(T) contains a subgraph depicted in Figure 3, where {j 1,...,j κ } [l]. v ej1 v fj1 v ej2 v fj2 v fjκ v ejκ v fj3 v ej3 Figure 3: A subgraph of G(T) from Claim 1. Let us illustrate Claim 1 by using the sample graph in Figure 2. Suppose that F = {(v e1,v f5 ),(v e2,v f2 ),(v e3,v f3 )}. Then T = T {e 1,e 2,e 3,}\{f 5,f 2,f 3 } is not a spanning tree and G(T) contains the subgraph composed of the following arcs (see Figure 2): (v e1,v f2 ),(v e2,v f2 ),(v e2,v f3 ),(v e3,v f3 ),(v e3,v f5 ),(v e1,v f5 ). Proof of Claim 1. WeformT byperformingasequenceofmoves consistinginaddingedgese i and removing edges f i P T (e i ), i [l]. Suppose that, at some step, a cycle appears, which is formed by some edges from {e 1,...,e l } and the remaining edges of T (not removed from T). Such a cycle must appear, since otherwise T would be a spanning tree. Let us relabel the edges so that {e 1,...,e s } are on this cycle, i.e. the first s moves consisting in adding 8

9 e i and removing f i create the cycle, i [s]. An example of such a situation for s = 4 is shown in Figure 4. The cycle is formed by the edges e 1,...,e 4 and the paths P v2 v 3, P v4 v 5 and P v1 v 6 in T. Consider the edge e 1 = {v 1,v 2 }. Because T is a spanning tree, P T (e 1 ) P v2 v 3 P T (e 2 ) P T (e 3 ) P v4 v 5 P T (e 4 ) P v1 v 6. Observe that f 1 P T (e 1 ) cannot belong to any of P v2 v 3, P v4 v 5 and P v1 v 6. If it would be contained in one of these paths, then no cycle would be created. Hence, f 1 must belong to P T (e 2 ) P T (e 3 ) P T (e 4 ). The above argument is general and, by using it, we can show that for each i [s], f i P T (e j ) for some j [s]\{i}. (a) f 2 (b) P v2v 3 P T (e 2 ) v f1 v e1 f 1 2 P T (e 1 ) f 1 2 P T (e 2 ) v 2 v 3 e 1 e 2 v 1 v f2 e 3 P v1v 6 v 4 P T (e 3 ) f 3 v e2 f 2 2 P T (e 2 ) f 2 2 P T (e 3 ) v 6 e 4 v 5 P T (e 4 ) f 4 P v4v 5 v f3 v e3 f 3 2 P T (e 3 ) f 3 2 P T (e 2 ) Figure 4: Illustration for the proof of Claim 1. (a) The bold lines represent paths in T (not necessarily disjoint); f 1 P T (e 2 ) P T (e 3 ) P T (e 4 ). (b) The subgraph G (T) with the corresponding cycle. We are now ready to build a subgraph depicted in Figure 3. Consider a subgraph G (T) of the cycle graph G(T) built as follows. The nodes of G (T) are v e1,...,v es,v f1,...,v fs. Observe that G (T) has exactly 2s nodes, since all the edges e 1,...,e s,f 1,...,f s are distinct by the assumption of the claim. For each i [s] we add to G (T) two arcs, namely (v ei,v fi ), f i P T (e i ) and (v ej,v fi ), f i P T (e j ) for j [s] \ {i} (see Figure 4). The resulting graph G (T) is bipartite and has exactly 2s arcs. In consequence G (T) (and thus G(T)) must contain a cycle which is of the form depicted in Figure 3. After this preliminary step, we can now return to the main proof. If E {e E : v e V A }, then, by the construction of the admissible graph, there exists a directed path in G A from a node v e, e E, to a node v f, f E. Let P be a shortest such a path from v e to v f, i.e. a path consisting of the fewest number of arcs, called an augmenting path. We need to consider the following cases: 1. The augmenting path P is of the form: E v e E v f If (v e,v f ) is -arc, then = {e}\{f} is an updated spanning tree of G such that = E Z +1. Furthermore is a minimum spanning tree for the costs C e and the new pair (,) satisfies the sufficient pair optimality conditions (E E, so 9

10 condition (ii) in Theorem 2 is not violated). If (v e,v f ) is -arc, then = {f}\{e} is an updated spanning tree of G such that = E Z +1. Also is a minimum spanning tree for the costs c e and the new pair (, ) satisfies the sufficient pair optimality conditions. An example can be seen in Figure 1. There is a path v e1 v e2 in the admissible graph. The arc (v e1,v e2 ) is both -arc and -arc. We can thus choose one of the two possible moves = {e 1 } \ {e 2 } or = {e 2 } \ {e 1 }, which results in (,) or (,). 2. The augmenting path P is of the form: E E Z E W E Z E W E Z E W E Z E (a) v e1 vf1 ve2 vf2 ve3 vf3 vel vfl vel+1 (b) E v fl Let = {e 1,...,e l }\{f 1,...,f l }. Let = {e 2,...,e l+1 }\{f 1,...f l } for case (a), and = {e 2,...,e l }\{f 1,...f l 1 } for case (b). We now have to show that the resultingpair(, )isapairofspanningtrees. Supposethat isnotaspanningtree. Observe that the -arcs (v e1,v f1 ),...,(v el,v fl ) belong to the cycle graph G(). Thus, by Claim 1, the cycle graph G() must contain a subgraph depicted in Figure 3, where {j 1,...,j κ } [l]. An easy verification shows that all edges e i,f i, i {j 1,...,j κ } must have the same costs with respect to Ce. Indeed, if some costs are different, then there exists an edge exchange which decreases the cost of. This contradicts our assumption that is a minimum spanning tree with respect to Ce. Finally, there must be an arc (v ei,v fi ) in the subgraph such that i < i. Since Ce i = C f, the arc (v i e i,v fi ) is present in the admissible graph G A. This leads to a contradiction with our assumption that P is an augmenting path. Now supposethat is not a spanningtree. We consider only the case (a) since the proof of case (b) is just the same. For a convenience, let us number the nodes v ei on P from i = 0 to l, so that = {e 1,...,e l }\{f 1,...,f l }. The arcs (v e1,v f1 ),...,(v el,v fl ), which correspond to the -arcs (v f1,v e1 ),...,(v fl,v el ) of P, belong to the cycle graph G(). Hence, by Claim 1, G() must contain a subgraph depicted in Figure 3, where {i 1,...,i κ } [l]. The rest of the proof is similar to the proof for. Namely, the edges e i and f i for i {i 1,...,i κ } must have the same costs with respect to c e. Also, there must exist an arc (v ei,v fi ) in the subgraph such that i > i. In consequence, the arc (v fi,v ei ) belongs to the admissible graph, which contradicts the assumption that P is an augmenting path. An example of the case (a) is shown in Figure 5. Thus = {e 1,e 2,e 3,e 4 } \ {f 1,f 2,f 3,f 4 } and = {e 2,e 3,e 4,e 5 } \ {f 1,f 2,f 3,f 4 }. An example of the case (b) is shown in Figure 6. In this example is the same as in the previous case and = {e 2,e 3,e 4 }\{f 1,f 2,f 3 }. It is easy to verify that E Z = = E Z + 1 holds (see also the examples in Figures 5 and 6). The spanning trees and are optimal for the costs C e and c e, respectively. Furthermore, E E and E E, so (, ) satisfies the sufficient pair optimality conditions (the condition (ii) in Theorem 2 is not violated). 10

11 (Ce;c e) (0,2) f 1 (2,0) f 2 (2,0) (2,0) f 3 f 4 (0,1) e 5 (1,0) e 1 e 2 e 3 e 4 (2,0) (0,2) (0,2) (0,2) E E Z E W E Z E W E Z E W E v e1 v f1 v e2 v f2 v e3 v f3 Z E v e4 v f4 v e5 Figure 5: A pair (,) and the corresponding admissible graph for the case 2a. (Ce;c e) (0,2) (1,0) e 1 f 1 (2,0) e 2 f 2 (2,0) e 3 f 3 (2,0) f 4 e 4 (1,2) (0,2) (0,2) (0,2) E Z E W E Z E W E Z E v e1 v W E f1 v e2 v f2 v e3 v f3 v e4 v f4 E Figure 6: A pair (,) and the corresponding admissible graph for the case 2b. 3. The augmenting path P is of the form E E W E Z E W E Z E W E Z E W E (a) v e1 vf1 ve2 vf2 ve3 vf3... vel vfl vel+1 (b) E v fl Let = {f 1,...,f l }\{e 2,...,e l+1 } for the case (a) and = {f 1,...,f l 1 }\ {e 2,...e l } for the case (b). Let = {f 1,...,f l } \ {e 1,...e l }. The proof that and are spanning trees follows by the same arguments as for the symmetric case described in point 2. An example of the case (a) is shown in Figure 7. Thus = {f 1,f 2,f 3,f 4 } \{e 2,e 3,e 4,e 5 } and = {f 1,f 2,f 3,f 4 } \ {e 1,e 2,e 3,e 4 }. An example for the case (b) is shown in Figure 8. The spanning tree is the same as in the previous case and = {f 1,f 2,f 3 } \ {e 2,e 3,e 4 }. The equality E Z = = E Z +1 holds. Also, the trees and are optimal for the costs C e and c e, respectively, E E, E E, so (, ) satisfies the sufficient pair optimality conditions. 11

12 (Ce;c e) (0,0) f 1 (0,2) f 2 (0,2) (0,2) f 3 f 4 (1,0) e 5 (2,1) e 1 e 2 e 3 e 4 (0,0) (2,0) (2,0) (2,0) E v E W E Z E W E Z E W E Z E e1 v f1 v e2 v f2 v W E e3 v f3 v e4 v f4 v e5 Figure 7: A pair (,) and the corresponding admissible graph for the case 3a. (Ce;c e) (0,0) e 1 (2,1) f 1 (0,2) e 2 f 2 (0,2) (0,2) f 3 e 3 f 4 e 4 (2,0) (2,0) (2,0) E v E W E Z E W E Z E W E e1 v f1 v e2 v f2 v Z E e3 v f3 v e4 v f4 Figure 8: A pair (,) and the corresponding admissible graph for the case 3b. We now turn to the case E {e E : v e V A } =. Fix δ > 0 (the precise value of δ will be specified later) and set: C e (δ) = C e δ, c e (δ) = c e v e V A, (23a) C e (δ) = C e, c e (δ) = c e δ v e / V A. (23b) Lemma 4. There exists a sufficiently small δ > 0 such that the costs C e (δ) and c e (δ) satisfy the path optimality conditions for and, respectively, i.e: for every e / C e (δ) C f (δ) for every f P (e), (24a) for every e / c e (δ) c f (δ) for every f P (e). (24b) Proof. If C e > C f (resp. c e > c f ), e /,f P (e) (resp. e /,f P (e)), then there is δ > 0, such that after setting the new costs (23) the inequality C e (δ) C f (δ) (resp. c e (δ) c f (δ)) holds. Hence, one can choose a sufficiently small δ > 0 such that after setting 12

13 the new costs (23), all the strong inequalities are not violated. Therefore, for such a chosen δ it remains to show that all originally tight inequalities in (22) are preserved for the new costs. Consider a tight inequality of the form: C e = C f, e /, f P (e). (25) On the contrary, suppose that C e (δ) < C f (δ). This is only possible when C e (δ) = C e δ and C f (δ) = C f. Hence and from the construction of the new costs, we have v f / V A (see (23b)) and v e V A (see (23a)). By (25), we obtain (v e,v f ) E A. Thus v f V A, a contradiction. Consider a tight inequality of the form: c e = c f, e /, f P (e). (26) On the contrary, suppose that c e (δ) < c f (δ). This is only possible when c e (δ) = c e δ and c f (δ) = c f. Thus we deduce that v e / V A and v f V A (see (23)). From (26), it follows that (v f,v e ) E A and so v e V A, a contradiction. We are now ready to give the precise value of δ. We do this by increasing the value of δ until some inequalities, originally not tight in (22), become tight. Namely, let δ > 0 be the smallest value of δ for which an inequality originally not tight becomes tight. Obviously, it occurs when C e δ = C f for e /, f P (e) or c f δ = c e for f /, e P (f). By (23), v e V A and v f / V A. Accordingly, if δ = δ, then at least one arc is added to G A. Observe also that no arc can be removed from G A - the admissibility of the nodes remains unchanged. It follows from the fact that each tight inequality for v e V A and v f V A is still tight. This leads to the following lemma. Lemma 5. If E {e E : v e V A } =, then (,) satisfies the sufficient pair optimality conditions for each θ [θ,θ+δ ]. Proof. Set θ = θ + δ, δ [0,δ ]. Lemma 4 implies that is optimal for C e (δ) and is optimal for c e (δ). From (23) and the definition of the costs C e and c e, it follows that C e (δ) = C e α e and c e (δ) = c e β e, where α e = α e + δ and β e = β e for each v e V A, α e = α e and β e = β e +δ for each v e / V A. Notice that α e+β e = α e +β e +δ = θ+δ = θ for each e E. By (23), c e (δ) = c e for each e E (recall that e E implies v e V A ), and thus β e = 0 for each e E. Since E {e E : v e V A } =, C e (δ) = C e = C e holds for each e E, and so α e = 0 for each e E. We thus have shown that there exist α e,β e 0 such that α e + β e = θ for each e E satisfying the conditions (i) and (ii) in Theorem 2, which completes the proof. We now describe a polynomial procedure that, for a given pair (,) satisfying the sufficient pairoptimality conditionsforsomeθ 0, findsanewpairofspanningtrees (, ), which also satisfies the sufficient pair optimality conditions with E Z = E Z +1. We start by building the admissible graph G A = (V A,E A ) for (,). If this graph contains an augmenting path, then by Lemma 3, we are done. Otherwise, we determine δ and modify the costs by using (23). Lemma 5 shows that (,) satisfies the sufficient pair optimality conditions for θ + δ. For δ some new arcs are added to the admissible graph G A (all the previous arcs must be still present in G A ). Thus G A is updated and we set C e := C e (δ ), c e := c e(δ ) for each e E, and θ := θ+δ. We repeat this until there is an augmenting path in G A = (V A,E A ). Note that such a path must appear after at most m = E iterations, 13

14 (C e;c e) (a) θ = 0 e 4 (0,9) e 1 (2,5) e 5 e 3 (8,3) (1,2) e 2 (0,6) e 6 (8,3) e 8 e 7 (2,3) e 9 (3,2) e 10 (5,9) (7,2) e 11 E v e8 E v e11 E v e3 E W v e7 e 4 (0,8) (7,3) e 1 (2,4) e e 5 8 e 3 (7,3) e E (b) 7 (1,3) (1,0) v e8 e 9 e 2 (3,1) E (0,5) e e 10 6 (5,8) v δ = 1;θ = 1 e11 (6,2) E e 11 v e3 E W v e7 E Z v e5 e 4 (0,7) (6,3) e 1 (2,3) e e 5 8 e 3 (6,3) e E (c) 7 (0,3) (1,-1) (0,1) v e8 e 9 e 2 (3,0) E (0,4) e v e 10 6 (5,7) e11 δ = 1;θ = 2 (5,2) E e 11 v e3 E W v e7 E v e10 E W v e1 E Z v e5 E v e6 Figure 9: Sample computations, = {e 2,e 4,e 5,e 6,e 9,e 10 } and = {e 2,e 3,e 5,e 8,e 9,e 11 }. which follows from the fact that at some step a node v e such that e E must appear in G A. Sample computations are shown in Figure 9. We start with the pair (,), where = {e 2,e 4,e 5,e 6,e 9,e 10 } and = {e 2,e 3,e 5,e 8,e 9,e 11 }, which satisfies the sufficient pair optimality conditions for θ = 0 (see Figure 9a). Observe that in this case it is enough to check that is optimal for the costs C e = C e and is optimal for the costs c e = c e, e E. For θ = 0, the admissible graph does not contain any augmenting path. We thus have to modify the costs C e and c e, according to (23). For δ = 1, a new inequality becomes tight and one arc is added to the admissible graph (see Figure 9b). The admissible graph still does not have an augmenting path, so we have to again modify the costs. For δ = 1 some new inequalities become tight and three arcs are added to the admissible graph (see Figure 9c). Now the admissible graph has two augmenting paths (cases 1 and 3a, see the proof of Lemma 3). Choosing one of them, and performing the modification described in the proof of Lemma 3 we get a new pair (, ) with E Z = E Z +1. Let us now estimate the running time of the procedure. The admissible graph has at most m nodes and at most mn arcs. It can be built in O(nm) time. The augmenting path in the admissible graph can be found in O(nm) time by applying the breath first search. Also the number of inequalities which must be analyzed to find δ is O(nm). Since we have to update the cost of each arc of the admissible graph at most m times, until an augmenting path appears, the required time of the procedure is O(m 2 n). We thus get the following result. 14

15 Theorem 3. The Rec ST problem is solvable in O(Lm 2 n) time, where L = n 1 k. 3 The recoverable robust spanning tree problem In this section we are concerned with the Rob Rec ST problem under the interval uncertainty representation, i.e. for the scenario sets U I, U1 I(Γ), and UI 2 (Γ). Using the polynomial algorithm for Rec ST, constructed in Section 2, we will provide a polynomial algorithm for Rob Rec ST under U I and some approximation algorithms for a wide class of Rob Rec ST under U1 I(Γ) and UI 2 (Γ). The idea will be to solve Rec ST for a suitably chosen second stage costs. Let F() = e C e +max min f(,s), S U Φ k where f(,s) = e cs e. It is worth pointing out that under scenario sets UI and U2 I(Γ), the value of F(), for a given spanning tree, can be computed in polynomial time [8, 17]. On the other hand, computing F() under U1 I (Γ) turns out to be strongly NP-hard [17, 9]. Given scenario S = (c S e ) e E, consider the following Rec ST problem: ( ) min C e + min f(,s). (27) Φ Φ k e e Problem (27) is equivalent to the formulation (1) for S = (c e ) e E and it is polynomially solvable, according to the result obtained in Section 2. As in the previous section, we denote by pair (,) a solution to (27), where Φ and Φ k. Given S, we call (,) an optimal pair under S if (,) is an optimal solution to (27). The Rob Rec ST problem with scenario set U I can be rewritten as follows: ( ) ( ) min C e + max min c S Φ S U I Φ k e = min C e + min (c e +d e ). (28) Φ Φ k e Thus (28) is (27) for S = (c e +d e ) e E U I. Hence and from Theorem 3 we immediately get the following theorem: Theorem 4. For scenario set U I, the Rob Rec ST problem is solvable in O((n 1 k)m 2 n) time. We now address Rob Rec ST under U I 1 (Γ) and UI 2 (Γ). Suppose that c e α(c e + d e ) for each e E, where α (0,1] is a given constant. This inequality means that for each edge e E the nominal cost c e is positive and c e +d e is at most 1/α greater than c e. It is reasonable to assume that this condition will be true in many practical applications for not very large value of 1/α. Lemma 6. Suppose that c e α(c e +d e ) for each e E, where α (0,1], and let ( ˆ,Ŷ) be an optimal pair under S = (c e ) e E. Then for the scenario sets U1 I(Γ) and UI 2 (Γ) the inequality F( ˆ) 1 αf() holds for any Φ. Proof. We give the proof only for the scenario set U1 I(Γ). The proof for UI 2 (Γ) is the same. Let Φ. The following inequality is satisfied: F() = e C e + max min f(,s) = C e +f(,s ) e +f( S U1 I(Γ) Φ k e e C,S). e e E 15

16 Clearly, (, ) is a feasible pair to (27) under S. From the definition of ( ˆ,Ŷ) we get F() e ˆ C e +f(ŷ,s) = e ˆ C e + e Ŷ c e e ˆ C e + e Ŷ α(c e +d e ) = e ˆC e +αf(ŷ,s), where S = (c e +d e ) e E. Hence F() ˆC e +α max f(ŷ,s) C e +α max S U1 I e (Γ) e ˆ min S U1 I(Γ) Φ kˆ f(,s) (29) α and the lemma follows. e ˆC e + max S U I 1 (Γ) min Φ kˆ f(,s) = αf( ˆ) The condition c e α(c e +d e ), e E, in Lemma 6, can be weakened and, in consequence, the set of instances to which the approximation ratio of the algorithm applies can be extended. Indeed, from inequality (29) it follows that the bounds of the uncertainty intervals are only required to meet the condition e Ŷ c e α e Ŷ (c e + d e ). This condition can be verified efficiently, since Ŷ can be computed in polynomial time. We now focus on Rob Rec ST for U2 I(Γ). Define D = e E d e and suppose that D > 0 (if D = 0, then the problem is equivalent to Rec ST for the second stage costs c e, e E). Consider scenario S under which c S e = min{c e + d e,c e + Γ de D } for each e E. Obviously, S U2 I(Γ), since e E δ e e E Γde D Γ. The following theorem provides another approximation result for Rob Rec ST with scenario set U2 I(Γ): Lemma 7. Let ( ˆ,Ŷ) be an optimal pair under S. Then the following implications are true for scenario set U I 2 (Γ): (i) If Γ βd, β (0,1], then F( ˆ) 1 βf() for any Φ. (ii) If Γ γf( ˆ), γ [0,1) then F( ˆ) 1 1 γf() for any Φ. Proof. Let Φ. Since S U2 I (Γ), we get F() = e C e + max min f(,s) e + min f(,s S U2 I(Γ) Φ k Φ e C ). (30) k We first prove implication (i). By (30) and the definition of ( ˆ,Ŷ), we obtain F() e ˆ e ˆ C e +f(ŷ,s ) = e ˆ C e + e Ŷ C e + e Ŷ min{c e +d e,c e +βd e } = e ˆ min{c e +d e,c e +Γ d e D } C e + e Ŷ (c e +βd e ) e ˆC e +βf(ŷ,s), where S = (c e +d e ) e E. The rest of the proof is the same as in the proof of Lemma 6. We now prove implication (ii). By (30) and the definition of ( ˆ,Ŷ), we have F() e ˆ C e +f(ŷ,s ) e ˆ C e +f(ŷ,s) ˆC e + max f(ŷ,s) Γ S U2 I (Γ) e ˆC e + max min f(,s) Γ = F( ˆ) Γ. S U2 I e (Γ) Φ kˆ 16

17 If Γ γf( ˆ). Then F() F( ˆ) γf( ˆ) = (1 γ)f( ˆ) and F( ˆ) 1 1 γ F(). Note that the value of F( ˆ) under U2 I (Γ) can be computed in polynomial time [17]. In consequence, the constants β and γ can be efficiently determined for every particular instance of the problem. Clearly, we can assume that d e Γ for each e E, which implies D mγ, where m = E. Hence, we can assume that Γ 1 md for every instance of the problem. We thus get from Lemma 7 (implication (i)) that F( ˆ) mf() for any Φ and the problem is approximable within m. If α, β and γ are the constants from Lemmas 6 and 7, then the following theorem summarizes the approximation results: Theorem 5. Rob Rec ST is approximable within 1 α under scenario set UI 1 (Γ) and it is approximable within min{ 1 β, 1 α, 1 1 γ } under scenario set UI 2 (Γ). Observe that Lemma 6 and Lemma 7 hold of any sets Φ and Φ k (the particular structure of these sets is not exploited). Hence the approximation algorithms can be applied to any problem for which the recoverable version (27) is polynomially solvable. 4 Conclusions In this paper we have studied the recoverable robust spanning tree problem (Rob Rec ST) under various interval uncertainty representations. The main result is the polynomial time combinatorial algorithm for the recoverable spanning tree. We have applied this algorithm for solving Rob Rec ST under the traditional uncertainty representation (see, e.g., [13]) in polynomial time. Moreover, we have used the algorithm for providing several approximation results forrec STwith thescenario set introducedby Bertsimas andsim [3] andthescenario set with a budged constraint (see, e.g,. [17]). There is a number of open questions concerning the considered problem. Perhaps, the most interesting one is to resolve the complexity of the robust problem under the interval uncertainty representation with budget constraint. It is possible that this problem may be solved in polynomial time by some extension of the algorithm constructed in this paper. One can also try to extend the algorithm for the more general recoverable matroid base problem, which has also been shown to be polynomially solvable in [10]. Acknowledgements The second and the third authors were supported by the National Center for Science (Narodowe Centrum Nauki), grant 2013/09/B/ST6/ References [1] R. K. Ahuja, T. L. Magnanti, and J. B. Orlin. Network Flows: theory, algorithms, and applications. Prentice Hall, Englewood Cliffs, New Jersey, [2] H. Aissi, C. Bazgan, and D. Vanderpooten. Approximation of min-max and min-max regret versions of some combinatorial optimization problems. European Journal of Operational Research, 179: ,

18 [3] D. Bertsimas and M. Sim. Robust discrete optimization and network flows. Mathematical Programming, 98:49 71, [4] C. Büsing. Recoverable robustness in combinatorial optimization. PhD thesis, Technical University of Berlin, Berlin, [5] C. Büsing. Recoverable robust shortest path problems. Networks, 59: , [6] C. Büsing, A. M. C. A. Koster, and M. Kutschka. Recoverable robust knapsacks: the discrete scenario case. Optimization Letters, 5: , [7] A. Chassein and M. Goerigk. On the recoverable robust traveling salesman problem. Optimization Letters, DOI /s , [8] O. Şeref, R. K. Ahuja, and J. B. Orlin. Incremental network optimization: theory and algorithms. Operations Research, 57: , [9] N. G. Frederickson and R. Solis-Oba. Increasing the weight of minimum spanning trees. Journal of Algorithms, 33: , [10] M. Hradovich, A. Kasperski, and P. Zieliński. The recoverable robust spanning tree problem with interval costs is polynomially solvable. Optimization Letters, DOI /s x, [11] A. Kasperski, A. Kurpisz, and P. Zieliński. Recoverable robust combinatorial optimization problems. In Operations Research Proceedings 2012, pages , [12] A. Kasperski and P. Zieliński. On the approximability of robust spanning problems. Theoretical Computer Science, 412: , [13] P. Kouvelis and G. u. Robust Discrete Optimization and its Applications. Kluwer Academic Publishers, [14] L. C. Lau, R. Ravi, and M. Singh. Iterative Methods in Combinatorial Optimization. Cambridge University Press, [15] C. Liebchen, M. E. Lübbecke, R. H. Möhring, and S. Stiller. The concept of recoverable robustness, linear programming recovery, and railway applications. In Robust and Online Large-Scale Optimization, volume 5868 of Lecture Notes in Computer Science, pages Springer-Verlag, [16] K. Lin and M. S. Chern. The most vital edges in the minimum spanning tree problem. Information Processing Letters, 45:25 31, [17] E. Nasrabadi and J. B. Orlin. Robust optimization with incremental recourse. CoRR, abs/ , [18] C. H. Papadimitriou and K. Steiglitz. Combinatorial optimization: algorithms and complexity. Dover Publications Inc.,