About this Journal Submit a Manuscript Table of Contents
Journal of Applied Mathematics
Volume 2012 (2012), Article ID 394721, 8 pages
doi:10.1155/2012/394721
Research Article

Polynomial Time Approximation Schemes for the Constrained Minimum Spanning Tree Problem

Department of Computer Science, Taipei Municipal University of Education, No. 1, Ai-Guo West Road, Taipei 10048, Taiwan

Received 28 October 2011; Accepted 18 January 2012

Academic Editor: Rudong Chen

Copyright © 2012 Yen Hung Chen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Let 𝐺 = ( 𝑉 , 𝐸 ) be an undirected graph with a weight function and a cost function on edges. The constrained minimum spanning tree problem is to find a minimum cost spanning tree T in G such that the total weight in T is at most a given bound B. In this paper, we present two polynomial time approximation schemes (PTASs) for the constrained minimum spanning tree problem.

1. Introduction

Motivated by the applications of quality of service (QoS) routing and multicasting, several multiple criteria problems have been studied (see [16] and references therein). Marathe et al. [5] studied a class of bicriteria network design problems that are defined as follows. Given an undirected graph 𝐺 = ( 𝑉 , 𝐸 ) and two independent minimization criteria with a bound on the first criterion, a generic bicriteria network design problem involves the minimization of the second criterion but satisfies the bound on the first criterion among all possible subgraphs from 𝐺 [5]. They considered three different criteria that have total edge cost, diameter, and maximum degree of a graph. The total edge cost of a graph is the sum of the costs of all edges in the subgraph. The diameter of a graph is the maximum distance between any pair of nodes in the subgraph. The degree of a node is the number of nodes adjacent to this node, and the maximum degree of a graph is maximum over the degrees of all nodes in the subgraph. Let 𝐺 = ( 𝑉 , 𝐸 ) be an undirected graph and a positive integer 𝐵 . Two different nonnegative functions, weight and cost, associate on edges in 𝐸 , respectively. The constrained minimum spanning tree (CMST) problem [710] is to find a minimum total cost spanning tree 𝑇 in 𝐺 such that the total weight in 𝑇 is at most 𝐵 . Aggarwal et al. [7] showed that the CMST problem is weakly NP-hard. Then Hong et al. [9] gave a pseudopolynomial time algorithm to solve the CMST problem. Moreover, some approximation algorithms have been proposed to solve the CMST problem [5, 8, 10].

Given two positive real numbers 𝛼 and 𝛽 , an ( 𝛼 , 𝛽 ) -approximation algorithm for the CMST problem is defined as a polynomial-time algorithm that produces a solution with the total weight at most 𝛼 times the bound 𝐵 and the total cost at most 𝛽 times the total cost of optimal solution for the CMST problem. The approximation ratio is denoted by ( 𝛼 , 𝛽 ) . An ( 𝛼 , 𝛽 ) -approximation algorithm for the CMST problem is called the polynomial time approximation scheme (PTAS) that produces a spanning tree with an approximation ratio of ( 1 , 1 + 𝜀 ) , for any constant 𝜀 > 0 . For any constant 𝜀 > 0 , Marathe et al. [5] gave a ( 1 + 𝜀 , 1 + 1 / 𝜀 ) -approximation algorithm for the CMST problem. Ravi and Goemans [10] presented a ( 2 , 1 ) -approximation algorithm to solve the CMST problem in 𝑂 ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) time and improved the ratio to ( 1 + 𝜀 , 1 ) , where the time complexity is 𝑂 ( | 𝑉 | 𝑂 ( 1 / 𝜀 ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) for any constant 𝜀 > 0 . These algorithms are based on the Lagrangian relaxation. Hassin and Levin [8] employed Ravi and Goemans' [10] algorithm with the matroid intersection strategy [1113] to design a ( 1 , 1 + 5 𝜀 ) -approximation algorithm in 𝑂 ( ( 1 / 𝜀 2 ) ( 1 / 𝜀 ) | 𝑉 | 3 ) time for the CMST problem. For other related bicriteria problems, Hassin [14] gave a fully polynomial time approximation scheme (FPTAS) for the restricted shortest path problem. The purpose of this problem is to find a minimum cost shortest path 𝑃 from a source to a destination subject to the constraint that the total weight in 𝑃 is at most 𝐵 . Then, Lorenz and Raz [15] gave an improved FPTAS whose time complexity is 𝑂 ( | 𝑉 | | 𝐸 | l o g l o g | 𝑉 | + ( | 𝑉 | | 𝐸 | ) / 𝜀 ) for the restricted shortest path problem. Further, Xue et al. [16] also designed an FPTAS with improved time complexity to 𝑂 ( | 𝑉 | | 𝐸 | l o g l o g l o g | 𝑉 | + ( | 𝑉 | | 𝐸 | ) / 𝜖 ) for the restricted shortest path problem. If the graph is directed acyclic, and Ergun et al. [17] provided an improved FPTAS for the restricted shortest path problem with running time 𝑂 ( ( | 𝑉 | | 𝐸 | ) / 𝜀 ) . Chen and Xue [18, 19] used a similar method to design PTASs for the 𝑘 -pair delay constrained minimum cost routing problem and the weight constrained Steiner tree problem in series-parallel graphs. These algorithms are based on the rounding and scaling strategy that is offered by Hassin [14]. In this paper, we present two polynomial time approximation schemes to find ( 1 , 1 + 𝜀 ) -approximation ratio for the CMST problem. Applying the second PTAS to Hassin and Levin's algorithm [8], the approximation ratio can be improved to ( 1 , 1 + 4 𝜀 ) for the CMST problem.

The rest of this paper is organized as follows. In Section 2, we first describe a ( 1 , ( 1 + 𝜀 ) 2 ) -approximation algorithm for the CMST problem and then improve the approximation ratio to ( 1 , 1 + 𝜀 ) . Finally, we make a conclusion in Section 3.

2. PTASs for the Constrained Minimum Spanning Tree Problem

In this section, we first clarify a ( 1 , ( 1 + 𝜀 ) 2 ) -approximation algorithm for the CMST problem. Hassin and Levin [8] first used this algorithm to obtain a ( 1 , 2 ( 1 + 𝜀 ) ) ratio for the CMST problem. Then they applied the Lagrangian relaxation and matroid intersection methods to improve the ratio to ( 1 , 1 + 5 𝜀 ) in 𝑂 ( ( 1 / 𝜀 2 ) ( 1 / 𝜀 ) | 𝑉 | 3 ) ) time. However, they did not describe the details of this algorithm. Then, we improve this approximation ratio into ( 1 , 1 + 𝜀 ) . Note that applying our algorithm, the approximation ratio of Hassin and Levin's algorithm [8] can be improved to ( 1 , 1 + 4 𝜀 ) ratio with the same time complexity. For convenience, we use 𝐶 ( 𝐻 ) (resp., 𝑊 ( 𝐻 ) ) to denote the sum of the costs (resp., weights) of all edges in any subgraph 𝐻 of 𝐺 . Let 𝑇 O P T be the optimal solution of the CMST problem. For any constant 𝜀 > 0 and an integer 𝐵 , Ravi and Goemans' [10] algorithm, denoted by 𝒜 𝑅 ( 𝜀 , 𝐵 ) , produces a spanning tree 𝑇 of the 𝑊 ( 𝑇 ) at most ( 1 + 𝜀 ) 𝐵 and of the 𝐶 ( 𝑇 ) at most 𝐶 ( 𝑇 O P T ). Our algorithm first finds a cost of the lower bound, denoted by L B , and a cost of the upper bound, denoted by U B , for the CMST problem. Then, iteratively shrink this range between L B and U B until U B ( 1 + 𝜀 ) L B . A trivial lower bound L B can be set to 1, and an upper bound U B can be set to | 𝑉 | 𝐶 m a x , in which 𝐶 m a x is the maximum edge cost in 𝐺 . However, the range between the trivial lower bound and upper bound is too large (i.e., not polynomial in input size). We reduce the range by Lorenz and Raz's method [15] such that UB | 𝑉 | LB. First, let 𝑐 1 , 𝑐 2 , , 𝑐 𝑘 be all the distinct costs of edges in 𝐺 . 𝐺 𝑗 = ( 𝑉 , 𝐸 𝑗 ) denotes the subgraph of 𝐺 ( 𝑉 , 𝐸 ) , where 𝐸 𝑗 is the set of edges with costs not greater than 𝑐 𝑗 . It is clear that 𝐺 𝑘 = 𝐺 and 𝐺 𝑗 𝐺 𝑗 + 1 for 1 𝑗 𝑘 1 . For 1 𝑗 𝑘 , we can use Prim's algorithm [20] to find a minimum weight spanning tree 𝑇 𝑗 in 𝐺 𝑗 . Let 𝐽 be the smallest index 𝑗 such that 𝑊 ( 𝑇 𝑗 ) is no more than 𝐵 . Then, finding the value 𝑐 𝐽 of a graph 𝐺 can be done in 𝑂 ( | 𝑉 | 2 l o g | 𝑉 | ) time by binary search and Prim's algorithm [20]. Therefore, we have 𝑐 𝐽 𝐶 ( 𝑇 O P T ) | 𝑉 | 𝑐 𝐽 . Hence, we can let the initial U B be | 𝑉 | 𝑐 𝐽 and L B be 𝑐 𝐽 . If U B ( 1 + 𝜀 ) L B , then 𝐶 ( 𝑇 𝐽 ) U B ( 1 + 𝜀 ) L B ( 1 + 𝜀 ) 𝐶 ( 𝑇 O P T ) . Otherwise, we first swap the cost (resp., weight) to the weight (resp., cost) for each edge in 𝐺 and find a value 𝑀 between L B and U B . Then, apply algorithm 𝒜 𝑅 ( 𝜀 , 𝑀 ) to find a spanning tree 𝑇 in 𝐺 . If the total weight of 𝑇 is no more than 𝐵 , let U B = ( 1 + 𝜀 ) 𝑀 ; otherwise, let L B = 𝑀 . Hence, the range between L B and U B is shrunk. Then, we repeatedly run the algorithm 𝒜 𝑅 ( 𝜀 , 𝑀 ) until U B ( 1 + 𝜀 ) L B . The next lemma shows that we can obtain either a lower bound 𝑀 or an upper bound ( 1 + 𝜀 ) 𝑀 for the CMST problem after running 𝒜 𝑅 ( 𝜀 , 𝑀 ) .

Lemma 2.1. Given an instance of the CMST problem with a constant 𝜀 > 0 and an integer 𝒞 , one swaps the cost (resp., weight) to the weight (resp., cost) for each edge and then applies algorithm 𝒜 𝑅 ( 𝜀 , 𝒞 ) to find a spanning tree 𝑇 . If 𝑊 ( 𝑇 ) 𝐵 , then there is an upper bound ( 1 + 𝜀 ) 𝒞 . Otherwise, one can obtain a lower bound 𝒞 .

Proof. For convenience, we use 𝑇 𝒞 to denote the minimum weight spanning tree of 𝐺 with total cost at most 𝒞 . After running the algorithm 𝒜 𝑅 ( 𝜀 , 𝒞 ) , we have a spanning tree 𝑇 with 𝑊 ( 𝑇 ) 𝑊 ( 𝑇 𝒞 ) and 𝒞 ( 𝑇 ) ( 1 + 𝜀 ) 𝒞 since the costs and weights of all edges are exchanged. Hence, if 𝑊 ( 𝑇 ) 𝐵 , then 𝒞 ( 𝑇 O P T ) 𝒞 ( 𝑇 ) ( 1 + 𝜀 ) 𝒞 . Otherwise, 𝑊 ( 𝑇 𝒞 ) 𝑊 ( 𝑇 ) > 𝐵 . Then we assume 𝒞 ( 𝑇 O P T ) < 𝐶 for contradiction. Hence, we have 𝑊 ( 𝑇 𝒞 ) 𝑊 ( 𝑇 O P T ) . However, 𝑊 ( 𝑇 O P T ) 𝐵 , and we have 𝑊 ( 𝑇 𝒞 ) 𝐵 . Hence, 𝒞 must be less than or equal to 𝐶 ( 𝑇 O P T ) .

For clarification, we describe the ( 1 , ( 1 + 𝜀 ) 2 ) -approximation algorithm for the CMST problem as follows.

Algorithm PTAS-CMST
Input: A graph 𝐺 = ( 𝑉 , 𝐸 ) with a weight function and a cost function on edges, an integer bound 𝐵 , and two real numbers 𝜀 > 0 and 𝜌 > 0 .
Output: A spanning tree 𝑇 A P X of 𝐺 subject to 𝑊 ( 𝑇 A P X ) 𝐵 . (1)For each edge in E, swap cost to weight and weight to cost, respectively. end for(2)Let an initial upper bound U B = | 𝑉 | 𝑐 𝐽 and an initial lower bound L B = 𝑐 𝐽 . (3)If   U B ( 1 + 𝜀 ) L B ,  then let 𝑇 A P X = 𝑇 𝐽 and return.Else(3.1)Let 𝐿 = L B and 𝑈 = U B / ( 1 + 𝜌 ) .(3.2)Repeat the following steps until the condition (1) 𝑈 ( 1 + 𝜀 ) 𝐿 or (2) 𝑈 𝐿 + 1 holds. (3.2.1) Let a middle value 𝑀 = 𝑈 𝐿 . (3.2.2) Use algorithm 𝒜 𝑅 ( 𝜌 , 𝑀 ) to find a spanning tree 𝑇 . (3.2.3)/* if 𝑊 ( 𝑇 ) > 𝐵 , then we have a lower bound 𝑀 ; else we have an upper bound ( 1 + 𝜌 ) 𝑀 and 𝐶 ( 𝑇 ) is less than or equal to this upper bound.  */ If   𝑊 ( 𝑇 ) 𝐵 ,  then let 𝑈 = 𝑀 and 𝑇 A P X = 𝑇 ;  else let 𝐿 = 𝑀 .

The approximation ratio of Algorithm PTAS-CMST is shown in the next lemma.

Lemma 2.2. Algorithm PTAS-CMST returns a spanning tree 𝑇 A P X with 𝐶 ( 𝑇 A P X ) ( 1 + 𝜀 ) ( 1 + 𝜌 ) 𝐶 ( 𝑇 O P T ) .

Proof. Let the upper bound 𝑈 𝑖 , the middle value 𝑀 𝑖 , and the lower bound 𝐿 𝑖 be parameters at the beginning of the 𝑖 th iteration of step (3.2) for Algorithm PTAS-CMST. Let 𝑇 𝑖 denote a spanning tree after running the 𝑖 -th iteration of step (3.2.2). It is clear that 𝐿 𝑖 𝑀 𝑖 𝑈 𝑖 . Hence, 𝑈 𝑖 + 1 𝑈 𝑖 and 𝐿 𝑖 + 1 𝐿 𝑖 . Moreover, we have 𝐿 𝑖 + 1 = 𝑈 𝑖 𝐿 𝑖 and 𝑈 𝑖 + 1 = 𝑈 𝑖 when 𝑊 ( 𝑇 𝑖 ) > 𝐵 . Then, we also have 𝑈 𝑖 + 1 = 𝑈 𝑖 𝐿 𝑖 and 𝐿 𝑖 + 1 = 𝐿 𝑖 when 𝑊 ( 𝑇 𝑖 ) 𝐵 . Further, 𝑈 𝑖 + 1 / 𝐿 𝑖 + 1 = 𝑈 𝑖 / 𝑈 𝑖 𝐿 𝑖 , or 𝑈 𝑖 + 1 / 𝐿 𝑖 + 1 = 𝑈 𝑖 𝐿 𝑖 / 𝐿 𝑖 . Therefore, we have 𝑈 𝑖 + 1 / 𝐿 𝑖 + 1 = 𝑈 𝑖 / 𝐿 𝑖 . Note that it successively narrows down the range between the lower bound 𝐿 and the upper bound 𝑈 after each iteration of step (3.2). After running the 𝑖 th iteration of step (3.2.3), we obtain either a lower bound 𝑀 𝑖 or an upper bound ( 1 + 𝜌 ) 𝑀 𝑖 by Lemma 2.1. Let 𝑘 be the number of iterations of step (3.2). After performing this algorithm, we have either 𝑈 𝑘 + 1 ( 1 + 𝜀 ) 𝐿 𝑘 + 1 or 𝑈 𝑘 + 1 𝐿 𝑘 + 1 + 1 . For the former case, we have 𝐶 𝑇 A P X ( 1 + 𝜌 ) 𝑈 𝑘 + 1 ( 1 + 𝜀 ) ( 1 + 𝜌 ) 𝐿 𝑘 + 1 𝑇 ( 1 + 𝜀 ) ( 1 + 𝜌 ) 𝐶 O P T . ( 2 . 1 ) For the latter case, we have 𝐶 𝑇 A P X ( 1 + 𝜌 ) 𝑈 𝑘 + 1 𝐿 ( 1 + 𝜌 ) 𝑘 + 1 𝐶 𝑇 + 1 ( 1 + 𝜌 ) O P T . + 1 ( 2 . 2 ) Let 𝜌 be 𝜀 , and hence Algorithm PTAS-CMST achieves an approximation ratio of ( 1 , ( 1 + 𝜀 ) 2 ) for the CMST problem.

Therefore, we have the following theorem.

Theorem 2.3. The CMST problem admits a PTAS. For any constant 𝜀 > 0 , a ( 1 , ( 1 + 𝜀 ) 2 ) -approximation algorithm for the CMST problem can be found in polynomial time.

Proof. By Lemma 2.2, Algorithm PTAS-CMST returns a ( 1 , ( 1 + 𝜀 ) 2 ) -approximation solution for the CMST problem. Next, we analyze the time complexity of Algorithm PTAS-CMST as follows. Let 𝑘 be the number of iterations of step (3.2). We also let 𝐿 1 and 𝑈 1 be the initial lower bound and upper bound, respectively. Since 𝑈 𝑖 + 1 / 𝐿 𝑖 + 1 = 𝑈 𝑖 / 𝐿 𝑖 , we have 𝑈 1 𝐿 1 , 𝑈 2 𝐿 2 = 𝑈 1 𝐿 1 1 / 2 , 𝑈 3 𝐿 3 = 𝑈 2 𝐿 2 1 / 2 = 𝑈 1 𝐿 1 ( 1 / 2 ) 2 , 𝑈 𝑘 𝐿 𝑘 = 𝑈 𝑘 1 𝐿 𝑘 1 1 / 2 𝑈 = = 1 𝐿 1 ( ( 1 / 2 ) 𝑘 1 ) . ( 2 . 3 ) Because ( 𝑈 𝑘 / 𝐿 𝑘 ) > ( 1 + 𝜀 ) and 𝑈 1 / 𝐿 1 ( | 𝑉 | / 1 + 𝜌 ) + 1 , we have 𝑘 l o g 1 / 2 l o g ( | 𝑉 | + 1 + 𝜌 ) / ( 1 + 𝜌 ) ( 1 + 𝜀 ) + 1 . Hence, 𝑘 is 𝑂 ( l o g l o g | 𝑉 | l o g l o g ( 1 + 𝜀 ) ) . Then, step (3.2.2) can be done in 𝑂 ( | 𝑉 | 𝑜 ( 1 / 𝜌 ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) time by Ravi's algorithm [10]. Hence, Algorithm PTAS-CMST is a polynomial time approximation scheme for the CMST problem.

Now, we improve the approximation ratio to ( 1 , 1 + 𝜀 ) for the CMST problem. This algorithm is the modification of the one designed in [17] for the restricted shortest path problem.

Algorithm Modify-PTAS-CMST
Input: A graph 𝐺 = ( 𝑉 , 𝐸 ) with a weight function and a cost function on edges, an integer bound 𝐵 , and two real numbers 𝜀 > 0 and 0 < 𝛾 < 1 .
Output: A spanning tree 𝑇 A P X of 𝐺 subject to 𝑊 ( 𝑇 A P X ) 𝐵 . (1)For each edge in 𝐸 , swap cost to weight and weight to cost, respectively. end for(2)Let an initial upper bound U B = | 𝑉 | 𝑐 𝐽 and an initial lower bound L B = 𝑐 𝐽 . (3)If U B ( 1 + 𝜀 ) L B , then let 𝑇 A P X = 𝑇 𝐽 and return. Else(3.1)Let 𝐿 = L B and 𝑈 = U B . (3.2) Repeat the following steps until the condition (1) 𝑈 ( 1 + 𝜀 ) 𝐿 or (2) 𝑈 𝐿 + 1 holds. (3.2.1) Let a real number 𝜌 = ( 𝑈 / 𝐿 ) 𝛾 1 and a middle value 𝑀 = 𝑈 𝐿 / ( 1 + 𝜌 ) . (3.2.2) Use algorithm A 𝑅 ( 𝜌 , 𝑀 ) to find a spanning tree 𝑇 .(3.2.3)/* if 𝑊 ( 𝑇 ) > 𝐵 , then we have a lower bound 𝑀 ; else we have an upper bound ( 1 + 𝜌 ) 𝑀 and 𝐶 ( 𝑇 ) is less than or equal to this upper bound.  */If   𝑊 ( 𝑇 ) 𝐵 ,  then let 𝑈 = ( 1 + 𝜌 ) 𝑀 and 𝑇 A P X = 𝑇 ;  else let 𝐿 = 𝑀 .

The approximation ratio of Algorithm Modify-PTAS-CMST is shown in the next lemma.

Lemma 2.4. The spanning tree 𝑇 A P X is a (1, 1 + 𝜀 )-approximation solution for the CMST problem.

Proof. Let the upper bound 𝑈 𝑖 , the middle value 𝑀 𝑖 , the lower bound 𝐿 𝑖 , and the real number 𝜌 𝑖 be parameters at the beginning of the 𝑖 th iteration of step (3.2) for Algorithm Modify-PTAS-CMST. Let 𝑇 𝑖 denote a spanning tree after running the 𝑖 -th iteration of step (3.2.2). If 𝑊 ( 𝑇 𝑖 ) 𝐵 , we have 𝑈 𝑖 + 1 = ( 1 + 𝜌 𝑖 ) 𝑀 𝑖 = ( 𝑈 𝑖 1 + 𝛾 𝐿 𝑖 1 𝛾 ) 1 / 2 and 𝐿 𝑖 + 1 = 𝐿 𝑖 . Otherwise, if 𝑊 ( 𝑇 𝑖 ) > 𝐵 , we have 𝐿 𝑖 + 1 = 𝑀 𝑖 = ( 𝑈 𝑖 1 𝛾 𝐿 𝑖 1 + 𝛾 ) 1 / 2 and 𝑈 𝑖 + 1 = 𝑈 𝑖 . Hence, ( 𝑈 𝑖 + 1 / 𝐿 𝑖 + 1 ) = ( 𝑈 𝑖 / 𝐿 𝑖 ) ( ( 𝛾 + 1 ) / 2 ) . Note that it successively narrows down the range between the lower bound 𝐿 and the upper bound 𝑈 after each iteration of step (3.2) since 0 < 𝛾 < 1 . Let 𝑘 be the number of iterations of step (3.2). After performing this algorithm, we have either 𝑈 𝑘 + 1 ( 1 + 𝜀 ) 𝐿 𝑘 + 1 or 𝑈 𝑘 + 1 𝐿 𝑘 + 1 + 1 . For the former case, we have 𝐶 𝑇 A P X 𝑈 𝑘 + 1 ( 1 + 𝜀 ) 𝐿 𝑘 + 1 𝑇 ( 1 + 𝜀 ) 𝐶 O P T . ( 2 . 4 ) For the latter case, we have 𝐶 𝑇 A P X 𝑈 𝑘 + 1 𝐶 𝑇 O P T + 1 . ( 2 . 5 )

Therefore, we have the following theorem.

Theorem 2.5. The CMST problem admits a PTAS. For any constant 𝜀 > 0 , a ( 1 , 1 + 𝜀 ) -approximation algorithm for the CMST problem can be found in polynomial time.

Proof. Clearly, Algorithm Modify-PTAS-CMST returns a ( 1 , 1 + 𝜀 ) -approximation solution for the CMST problem by Lemma 2.4. Next, we analyze the time complexity of Algorithm Modify-PTAS-CMST as follows. Let 𝑘 be the number of iterations of step (3.2). We also let 𝐿 1 and 𝑈 1 be the initial lower bound and upper bound, respectively. Then, we have 𝑈 1 𝐿 1 , 𝑈 2 𝐿 2 = 𝑈 1 𝐿 1 ( 𝛾 + 1 ) / 2 , 𝑈 3 𝐿 3 = 𝑈 2 𝐿 2 ( 𝛾 + 1 ) / 2 = 𝑈 1 𝐿 1 ( ( ( 𝛾 + 1 ) / 2 ) 2 ) , 𝑈 𝑘 𝐿 𝑘 = 𝑈 𝑘 1 𝐿 𝑘 1 ( 𝛾 + 1 ) / 2 𝑈 = = 1 𝐿 1 ( ( ( 𝛾 + 1 ) / 2 ) 𝑘 1 ) ( 2 . 6 ) since 𝑈 𝑖 + 1 / 𝐿 𝑖 + 1 = ( 𝑈 𝑖 / 𝐿 𝑖 ) ( 𝛾 + 1 / 2 ) .
Initially, 𝑈 1 / 𝐿 1 = | 𝑉 | . After running Algorithm Modify-PTAS-CMST, we have ( 𝑈 𝑘 + 1 / 𝐿 𝑘 + 1 ) ( 1 + 𝜀 ) . Step (3.2.2) can be done in 𝑂 ( | 𝑉 | 𝑜 ( 1 / 𝜌 ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) time by Ravi's algorithm [10]. Let 𝐼 be the smallest index 𝑗 such that | 𝑉 | 𝑜 ( 1 / 𝜌 𝑗 ) > | 𝑉 | , for 1 𝑗 𝑘 . For all 1 𝑗 𝐼 1 , the time is 𝐼 1 𝑗 = 1 𝑂 ( | 𝑉 | 𝑜 ( 1 / 𝜌 𝑗 ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) . It is clear that ( 𝑈 𝐼 / 𝐿 𝐼 ) < 2 1 / 𝛾 and ( 𝑈 𝐼 1 / 𝐿 𝐼 1 ) 2 1 / 𝛾 . Hence, we have l o g ( 𝛾 + 1 ) / 2 ( 1 / ( 𝛾 l o g 2 | 𝑉 | ) ) + 1 𝐼 1 > l o g ( 𝛾 + 1 ) / 2 ( 1 / ( 𝛾 l o g 2 | 𝑉 | ) ) since ( 𝑈 𝐼 / 𝐿 𝐼 ) = ( 𝑈 1 / 𝐿 1 ) ( ( ( 𝛾 + 1 ) / 2 ) 𝐼 1 ) . If 𝐼 1 1 , the complexity is at most 𝑂 ( | 𝑉 | ( 1 l o g ( ( 𝛾 + 1 ) / 2 ) 𝛾 l o g ( 𝛾 + 1 / 2 ) ( l o g 2 | 𝑉 | ) ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) , for all 1 𝑗 𝐼 1 . Then, for all 𝐼 𝑗 𝑘 , the time complexity is 𝑘 𝑗 = 𝐼 𝑂 ( | 𝑉 | 𝑜 ( 1 / 𝜌 𝑗 ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) . Since ( 𝑈 𝑘 / 𝐿 𝑘 ) > ( 1 + 𝜀 ) with ( 𝑈 𝑘 / 𝐿 𝑘 ) = ( 𝑈 1 / 𝐿 1 ) ( ( ( 𝛾 + 1 ) / 2 ) ( k 1 ) ) , we have 𝑘 < l o g ( 𝛾 + 1 ) / 2 ( l o g ( 1 + 𝜀 ) / l o g | 𝑉 | ) + 1 . Moreover, ( 1 / 𝜌 𝐼 ) ( 1 / ( 𝜌 𝐼 + 1 ) ) ( 1 / 𝜌 𝑘 ) . Hence we have the time complexity at most 𝑂 ( | 𝑉 | 𝑜 ( 1 / 𝜌 𝑘 ) ( 1 + l o g ( 𝛾 + 1 ) / 2 l o g ( 1 + 𝜀 ) + l o g ( 𝛾 + 1 ) / 2 𝛾 ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) for all 𝐼 𝑗 𝑘 . Because 𝛾 is a constant, the total time complexity is at most 𝑂 ( ( | 𝑉 | l o g l o g | 𝑉 | + | 𝑉 | 1 / ( ( 1 + 𝜀 ) 𝛾 1 ) l o g l o g ( 1 + 𝜀 ) ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) , where 1 / 𝜌 𝑘 1 / ( ( 1 + 𝜀 ) 𝛾 1 ) , 0 < 𝛾 < 1 .

Let 𝜀 = 1 and 𝛾 1 . Use Algorithm Modify-PTAS-CMST to find a feasible solution for the CMST problem with ratio (1,2). The total time complexity of this algorithm is at most 𝑂 ( ( | 𝑉 | 1 / ( 2 𝛾 1 ) + | 𝑉 | l o g l o g | 𝑉 | ) ( | 𝐸 | l o g 2 | 𝑉 | + | 𝑉 | l o g 3 | 𝑉 | ) ) . Then, we can apply Hassin and Levin's [8] algorithm (i.e., the Lagrangian relaxation and matroid intersection), and the approximation ratio can be improved to ( 1 , 1 + 4 𝜀 ) with the same time complexity of Hassin and Levin's algorithm (i.e., 𝑂 ( ( 1 / 𝜀 2 ) ( 1 / 𝜀 ) | 𝑉 | 3 ) ).

3. Conclusion

In this paper, we presented two polynomial time approximation schemes for the CMST problem. Using the second PTAS with the Lagrangian relaxation and matroid intersection, the approximation ratio of Hassin and Levin's algorithm [8] can be improved to the ratio of ( 1 , 1 + 4 𝜀 ) . A open problem left in the paper could involve studying whether there exists a fully polynomial time approximation scheme for the CMST problem.

Acknowledgment

This work was supported in part by the National Science Council of the Republic of China under Contract NSC99-2221-E-133-003.

References

  1. C. H. Chow, “On multicast path finding algorithms,” in Proceedings of the 10th Annual Joint Conference of the IEEE and Communications Societies (IEEE INFOCOM '91), vol. 3, pp. 1274–1283, Bal Harbour, Fla, USA, 1991.
  2. A. Goel, K. G. Ramakrishnan, D. Kataria, and D. Logothetis, “Efficient computation of delay-sensitive routes from one source to all destinations,” in Proceedings of the 20th Annual Joint Conference of the IEEE Computer and Communications Societies, vol. 2, pp. 854–858, Anchorage, Alaska, USA, 2001.
  3. B. Kadaba and J. Jaffe, “Routing to multiple destinations in computer networks,” IEEE Transactions on Communications, vol. 31, pp. 343–351, 1983.
  4. V. P. Kompella, J. C. Pasquale, and G. C. Polyzos, “Multicasting for multimedia applications,” IEEE/ACM Transactions on Networking, pp. 286–292, 1993.
  5. M. V. Marathe, R. Ravi, R. Sundaram, S. S. Ravi, D. J. Rosenkrantz, and H. B. Hunt, III, “Bicriteria network design problems,” Journal of Algorithms, vol. 28, no. 1, pp. 142–171, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  6. Z. Wang, “On the complexity of quality of service routing,” Information Processing Letters, vol. 69, no. 3, pp. 111–114, 1999. View at Publisher · View at Google Scholar · View at MathSciNet
  7. V. Aggarwal, Y. P. Aneja, and K. P. K. Nair, “Minimal spanning tree subject to a side constraint,” Computers and Operations Research, vol. 9, no. 4, pp. 287–296, 1982.
  8. R. Hassin and A. Levin, “An efficient polynomial time approximation scheme for the constrained minimum spanning tree problem using matroid intersection,” SIAM Journal on Computing, vol. 33, no. 2, pp. 261–268, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. S.-P. Hong, S.-J. Chung, and B. K. Park, “A fully polynomial bicriteria approximation scheme for the constrained spanning tree problem,” Operations Research Letters, vol. 32, no. 3, pp. 233–239, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. R. Ravi and M. Goemans, “The constrained minimum spanning tree problem,” in Proceedings of the 5th Scandinavian Workshop on Algorithm Theory, vol. 1097 of Lecture Notes in Computer Science, pp. 66–75, 1996.
  11. P. Chalasani and R. Motwani, “Approximating capacitated routing and delivery problems,” SIAM Journal on Computing, vol. 28, no. 6, pp. 2133–2149, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. W. J. Cook, W. H. Cunningham, W. R. Pulleyblank, and A. Schrijver, Combinatorial optimization, Wiley-Interscience Series in Discrete Mathematics and Optimization, John Wiley & Sons, New York, NY, USA, 1998, A Wiley-Interscience Publication.
  13. E. L. Lawler, Combinatorial Optimization: Networks and Matroids, Holt, New York, NY, USA, 1976.
  14. R. Hassin, “Approximation schemes for the restricted shortest path problem,” Mathematics of Operations Research, vol. 17, no. 1, pp. 36–42, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  15. D. H. Lorenz and D. Raz, “A simple efficient approximation scheme for the restricted shortest path problem,” Operations Research Letters, vol. 28, no. 5, pp. 213–219, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. G. Xue, W. Zhang, J. Tang, and K. Thulasiraman, “Polynomial time approximation algorithms for multi-constrained QoS routing,” IEEE/ACM Transactions on Networking, vol. 16, pp. 656–669, 2008.
  17. F. Ergun, R. Sinha, and L. Zhang, “An improved FPTAS for restricted shortest path,” Information Processing Letters, vol. 83, no. 5, pp. 287–291, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  18. G. Chen and G. Xue, “K-pair delay constrained minimum cost routing in undirected networks,” in Proceedings of the Twelfth Annual ACM-SIAM Symposium on Discrete Algorithms (SODA '01), pp. 230–231, SIAM, Washington, DC, USA,, 2001.
  19. G. Chen and G. Xue, “A PTAS for weight constrained Steiner trees in series-parallel graphs,” Theoretical Computer Science, vol. 304, no. 1–3, pp. 237–247, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  20. T. H. Cormen, C. E. Leiserson, R. L. Rivest, and C. Stein, Introduction to Algorithms, MIT Press, Cambridge, Mass, USA, 2nd edition, 2001.