Table of Contents Author Guidelines Submit a Manuscript
Journal of Applied Mathematics
Volume 2013, Article ID 986985, 11 pages
Research Article

A Rough Penalty Genetic Algorithm for Multicast Routing in Mobile Ad Hoc Networks

Department of Information Management, Chung Yuan Christian University, Jhongli City 32023, Taiwan

Received 26 April 2013; Accepted 13 July 2013

Academic Editor: Anyi Chen

Copyright © 2013 Chih-Hao Lin and Chia-Chun Chuang. 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.


Multicast routing is an effective way to transmit messages to multiple hosts in a network. However, it is vulnerable to intermittent connectivity property in mobile ad hoc network (MANET) especially for multimedia applications, which have some quality of service (QoS) requirements. The goal of QoS provisioning is to well organize network resources to satisfy the QoS requirement and achieve good network delivery services. However, there remains a challenge to provide QoS solutions and maintain end-to-end QoS with user mobility. In this paper, a novel penalty adjustment method based on the rough set theory is proposed to deal with path-delay constraints for multicast routing problems in MANETs. We formulate the problem as a constrained optimization problem, where the objective function is to minimize the total cost of the multicast tree subject to QoS constraints. The RPGA is evaluated on three multicast scenarios and compared with two state-of-the-art methods in terms of cost, success rate, and time complexity. The performance analyses show that this approach is a self-adaptive method for penalty adjustment. Remarkably, the method can address a variety of constrained multicast routing problems even though the initial routes do not satisfy all QoS requirements.

1. Introduction

Multicasting is a service method in which a source node can deliver copies of messages to multiple recipients at different locations in a communication network. Multicasting techniques play a critical role in many applications such as video conference, internet games, and web-based learning. In this paper, multicast routing problems mainly focus on finding a minimum Steiner tree and satisfying quality-of-service (QoS) requirements. Unfortunately, the problem of finding a Steiner tree is known to be a NP-complete problem [1], even if links have unit costs.

The multicast tree in mobile ad hoc networks (MANETs) is vulnerable to intermittent connectivity property during the transmission period [2]. Due to the nodal mobility and the dynamic topology of mobile networks, a service provider should find a cost-effective tree for its multicast customers in real time and assure certain QoS requirements [3]. Without the support of communication infrastructure, the challenge of the dynamic routing on a changing topology is how to decide a multicast tree as soon as possible. To enhance the effectiveness and efficiency, this work uses a routing representation scheme to encode the multicast tree. Therefore, this paper models the multicast routing problem with QoS constraints to support multimedia transmission in MANETs [4, 5].

The battery limitation of a mobile node is a critical constraint while developing multicast routing protocols. Genetic algorithm (GA) presents a potential solution for the multiconstrained multicast routing problem [6]. Traditionally, penalty-function methods are the most popular constraint-handling techniques. According to the degree of constraint violation, penalty coefficient should be determined carefully [7]. The static-penalty (SP) method applies a static coefficient for each constraint and then adjusts penalty coefficient manually [8]. To adjust penalty coefficient automatically, the dynamic-penalty (DP) method combines the generation number and a scaling constant to adjust coefficient automatically [9]. Furthermore, the adaptive-penalties (AP) method tries to avoid infeasible solutions by adjusting the penalty coefficient according to the convergent situation [10].

Different from several related researches [2, 4, 5], this paper synthesizes the rough set theory (RST) and penalized techniques as a rough penalty genetic algorithm (RPGA). The emphasis of the proposed RPGA uses a rough-penalty (RP) method to releases/enforces some penalties on inefficient/efficient constraints during evolution. To facilitate the effectiveness of multiple penalties, the RPGA incorporates with a therapeutic crossover to enlarge the genetic diversity in a population and guide to find the optimal solution. The performance of the RPGA is evaluated by three kinds of constrained multicast routing scenarios. Experimental results show that the proposed RPGA not only can find near-optimal solutions but also can obtain robust feasible results for QoS-based multicast routing problems.

The rest of paper is organized as follows. Section 2 models the multicast routing problems with QoS constraints in MANETs. In Section 3, the operations of the proposed RPGA are described in detail. Section 4 introduces the proposed RP method. Section 5 reports the experimental results, algorithm analyses, and performance comparisons for three test networks. Finally, the paper is summarized in Section 6.

2. Problem Description and Formulation

2.1. Network Modeling

At a certain time period in a MANET, we assume that the service provider knows traffic load from node to node (), QoS requirements/constraints (e.g., delay constraint), bandwidth available of each link, link cost for the traffic to pass through each link, and time delay for the traffic to pass through each OD pair. We consider the multicast routing problem with bandwidth and delay constraints. The communication network is modeled as a connected weighted, directed graph , where is a finite set of network nodes and is the set of network links. The link connects node to node with positive cost function (i.e., ), available bandwidth (i.e., ), and delay function (i.e., ). The number of nodes and links (i.e., the cardinalities of and ) is and , respectively. For each multicast session, messages are routed from a source node to a set of multicast destination group . A multicast tree represents a solution to the multicast routing problem, where and . Thus, this tree is a subgraph of with root and a set of nodes . Let represent a path in the tree, , from source node to a destination node . We have the following definitions.

Definition 1. The cost of multicast tree is the sum of the links’ cost in the tree. The link cost may represent its monetary cost or resource utilization:

Definition 2. The bottleneck bandwidth of path is the minimum value of links’ bandwidth along the path, which represents the residual bandwidth of a communication path:

Definition 3. The delay of path is the sum of links’ delay along the path from to . The link delay may include its nodal processing, queueing, transmission, and propagation delays:

Definition 4. The delay of multicast tree is the maximum value of paths’ delay in the tree, that is,

Figure 1 depicts an example of a network graph, link parameters, and a multicast tree. Parameters along links are triple (cost, delay, bandwidth). In Figure 1, the source node is Node 1 (i.e., ). The destination nodes are nodes 3, 7, and 8 (i.e., ). The Steiner tree consists of three paths , , and . The total cost of the Steiner tree can be calculated by (1), that is, . By using (2), we can calculate the bottleneck bandwidth of each path along the Steiner tree, that is, , , and . We can also apply (3) to calculate , , and . Finally, we can derive the total delay of the Steiner tree by (4).

Figure 1: An example of a network graph, link parameters, and a Steiner tree.
2.2. Problem Definition

The optimal multicast tree depends on the operator objectives, such as network cost, transmission delay, or target QoS [11]. Therefore, solving the multicast routing problem is equivalent to find the optimal distribution tree on the basis of a certain cost function under a given set of constraints. In this paper, the multicast routing in MANETs can be modeled as a combinatorial optimization problem in the following:

The objective function is to minimize the total cost of the multicast tree. The path-delay constraint enforces that the total delay of each OD pair must be smaller than or equal to its delay bound . The minimum bandwidth requirement is denoted as . Therefore, the multicast routing problem is to determine a multicast tree connecting the source node to every destination node such that the cost of this tree is minimum, while the path-delay and bottleneck-bandwidth from the source node to any destination node satisfy the prescribed QoS requirements.

2.3. Routing Table

Since there are so many candidate paths between two nodes in the network graph , traditional GAs may consume considerable computational effort in searching infeasible solutions because genetic operations do not always preserve feasibility. Therefore, to reduce the search space, this work uses the shortest path routing algorithm to precalculate the first shortest paths and record in a routing table. For the network topology in Figure 1, an example of its routing table for OD pair (1,7) will look like Table 1, which includes the first shortest paths from node 1 to node 7 with the route path, total cost, aggregate delay, and bottleneck bandwidth.

Table 1: An example of routing table for OD pair (1, 7) in Figure 1.

3. RPGA for Multicasting Routing Problem

The proposed RPGA adopts a RP method to enhance the searching abilities of original GAs for handling constrained multicasting routing problems. To enhance the exploration ability, the RPGA adopts the RST to enlarge the genetic diversity by releasing inefficient constraints and also enforcing efficient ones when the generation number is odd. The flow chart of RPGA (in Figure 2) consists of several genetic operations, such as initialization, selection, crossover, mutation, replacement, and RP method. To enhance the exploitation ability, the proposed RPGA applies the therapeutic crossover to improve the convergence rate during the evolution.

Figure 2: Flow chart of the RPGA.
3.1. Encoding and Initialization

In this paper, the encoding method is based on a routing representation for multicast trees. The RPGA maintains a population of chromosomes, which represent a candidate set of Steiner trees for the multicast routing problem. Given a source node and a set of destination nodes , a chromosome can be represented by a string of integers with length . The chromosome is denoted as , where is an integer in interval to represent a route number for the OD pair from to in the routing table. In the example of Figure 1, the second OD pair (i.e., source node 1 to destination node 7) is routed along Path 1-4-5-7. The route number of this path in its routing table (in Table 1) is 2. Thus, we should assign the route number as the value of the second gene, that is, . Therefore, the relationship between chromosome, gene, and routing table for the example in Figure 1 can be illustrated as Figure 3. The RPGA starts with a random population within the gene value interval , no matter whether these multicast trees satisfy the QoS constraints or not. The emphasis is that the RPGA can find the global optimum for constrained problems even though the initial population is infeasible.

Figure 3: Representation of chromosome, genes, and routing table.
3.2. Fitness Function with Self-Adaptive Penalty Adjustment

The fitness value of each chromosome represents the quality of the corresponding multicast tree (i.e., ). However, the penalty adjustment for QoS constraints is difficult to adapt suitably. In this paper, we mix the aspects of the Joines and Houck’s DP method [9] with the RST to find a Steiner tree that satisfies the QoS constraints and minimizes total routing cost as well.

The proposed RP method adjusts penalty terms according to both violation magnitude and evolution time. To solve constrained optimization problems effectively, each individual in generation is evaluated using an expanded objective function (6): where is a “severity” factor, is the total number of constraints, and is the violation magnitude of constraint . This fitness function combines a coefficient with an exponent to increase penalty pressure over time. For constraint in generation , the exponent is a representative penalty multiplier that is initially assigned as 2. And then, the penalty multiplier is tuned iteratively according to the discernible mask and the representative attribute value of superior class . The exponent is defined as Remarkably, the discernible mask can be used to enable by differentiating significant characteristics between classes and . If the th constraint is discernible (i.e., ), the exponent is adjusted by the representative attribute value (); otherwise, the exponent retains the same value as in the previous generation. All these RP coefficients will be introduced in Section 4.

3.3. Selection Operation

A selection operation uses fitness to determine the solution quality and to select high-quality chromosomes for the recombination operation [12]. The RPGA employs a stochastic universal selection to create selection pressure towards the global optimal solution. The measurement of a chromosome’s fitness is its value of the expanded objective function in (6).

3.4. Crossover Operation

The crossover operation represents the mixing of genetic material from two selected parents to produce one child chromosome. The RPGA proposes a therapeutic crossover that incorporates a gene-therapy method with a conventional crossover scheme to enhance the exploitation ability and speed up the convergence rate [13].

Each time the selection operation chooses two crossover parents from the population. The therapeutic crossover gives each gene locus an equal chance of being a crossover point (i.e., belongs to a therapeutic genome where ). The proposed gene-therapy method evaluates the merit of two selected genomes by comparing the changes in the chromosome fitness before and after interchanging the genomes with the other mating chromosome [13].

According to their relative merit, these two genomes combine to generate a new genome for their offspring. Therefore, offspring inherit more genetic material from the superior genome than the inferior one. We depict the pseudo-code of the therapeutic crossover in Pseudocode 1.

Pseudocode 1: Pseudocode of the crossover operation for routing problem.

3.5. Mutation Operation

A mutation operation used in GAs can increase population diversity to enhance its exploration ability [13]. This work uses the bit-flip mutation with a fixed small probability . According to this mutation probability , the mutation operation randomly selects a subset of genes and chooses new paths from its routing table. Thus, the route numbers of these new paths in its routing table replace the original values of selected genes. The resulting chromosome is a new multicast tree and can increase population diversity.

3.6. Replacement and Termination

The proposed RPGA adopts a replacement-with-elitism method to prevent best solutions from being lost through a selection process. A successive population is produced from three sources: the replacement-with-elitism method selects the best 10% chromosomes to join the new population; the crossover operation recombines 80% of child chromosomes; and the mutation operation constructs other child chromosomes for the next generation. The RPGA will stop when it reaches the predefined maximum iterations.

4. RP Method for Constraint Handling

To address the multicast routing problem with QoS constraints, the challenge is how to optimize the objective function value against its constraint violations. Traditional GAs are ineffective in searching feasible solutions because genetic operations do not always preserve feasibility [14]. Therefore, penalty-function methods are the most popular constraint-handling techniques for constrained optimization [15].

The novel RP method has been proposed in our previous work for numerical constrained problems [13]. This proposed RPGA inspired by the Pawlak’s RST [16, 17] has been proved better than several existing algorithms for solving a variety of numerical optimization problems [13]. The RP method can automatically adjust penalty coefficients during the evolution. Furthermore, the method does not depend on extra functional analyses for its solution space. Therefore, this study aims to effectively extend the original RPGA as a new constraint-handling technique to address the multicast routing problem with QoS constraints in MANETs. During the genetic evolution, the proposed RP method uses the attribute reduction concept to find appropriate penalizing strategies and release some inefficient constraints.

4.1. Flow of RP Method

The pseudocode of the RP method is depicted in Pseudocode 2. The proposed RP method not only penalizes constraint violations to exploit feasible space but also releases ineffective constraints to explore infeasible space. Therefore, the RP method is a self-adaptive approach that can measure infeasibility and can adjust each penalty coefficient automatically.

Pseudocode 2: Pseudocode of the RP method.

4.2. Rough Penalty Classification

This work uses information granulation as a key function for implementing a divide-and-conquer strategy. Elementary information granules are indiscernibility classes of constraint violations. The information system is an information table of attribute values containing rows labeled by objects and columns labeled by attributes [18].

Remark 5. A partition granularity is defined for classifying the magnitude of constraint violations. The design principle is that solution quality increases as its constraint penalty moves closer to zero. Therefore, this study uses a smaller range in near-zero regions than in other regions.

4.3. Rough Decision System

A decision system is an IS with the form in which each individual is treated as an object of a nonempty finite set . Attribute set is a nonempty finite set of conditional attributes, where each penalty multiplier () corresponds to a conditional attribute. The supervised knowledge is expressed by a decision attribute (denoted by ). An information function maps each object to a decision attribute, that is, for , . For a minimization problem, the information function is designed as follows: where

Remark 6. Because the penalty multiplier should be adjusted according to the region of its constraint violation, this work enlarges the penalty multiplier when its violation level increases. In the illustration in Figure 4, the penalty multiplier will be assigned as , , , , , and for constraint regions 1, 2, 3, 4, 5, and 6, where denotes the coefficient and the partition granularity is six.

Figure 4: Penalty multiplier classification.
4.4. Significant Penalty and Attribute Reduction

Based on the concept of attribute reduction, attributes may not be equally important, and some of them can be eliminated from a decision table without degrading information quality. Attribute reduction can be generalized by introducing attribute evaluation, which can express the merit of each attribute in the information table [19].

Remark 7. Decision attribute in determines a partition of object set , where is the object classification with respect to decision attribute . The minimal subset of penalized constraints is applied to distinguish above-average individuals (i.e., their decision attributes are “1”) and below-average ones (i.e., those values are “0”). The representative value of each relevant attribute is assigned as the attribute value with the maximum cardinality in the same class.

5. Computational Experiments

5.1. Test Platform and Parameter Setting

In this paper, the proposed RPGA is evaluated by solving multicast routing problems in MANETs. We use the well-known network generation tool [20] to create an asynchronous network based on the Waxman’s techniques [21] and depicted in Figure 5. The network illustrates a random graph in which 40 nodes are connected, and each node has average of four connections to other nodes (i.e., average degree of a node is 4). Each link has its own cost, delay, and bandwidth information. In order to reduce the complexity of the graph representation, Figure 5 only shows the cost/delay information along one direction link (from a smaller ID node to a larger ID one). For example, traffic along the link from node 20 to node 34 will spend 54 units cost and delay 14 msec. In all cases, the maximum number of iterations is 40, the population size is 20, and the number of elite individuals is 2. For each test problem, 30 independent runs with different seeds are performed using the MATLAB environment.

Figure 5: A randomly generated network with 40 nodes and average degree four.
5.2. Performance Metrics

In this paper, the performance metrics of solution algorithms consist of the total cost of the obtained multicast tree; the success rate with respect to the QoS constraints; and the required CPU time for computing the multicast routing problems. The success rate () concerns about the percentage of feasible routes with respect to the QoS requirements. We can define the success rate as follows [22]: where is the number of OD pairs that satisfy all QoS constraints and is the total number of OD pairs in this multicast group.

5.3. Algorithm Analyses for Different Parameter Settings

It is well known that the performance of GAs significantly depends on the configuration of its operating parameters. To investigate the impact of various parameter settings in the RP method, this study experiments on two parameters: the partition granularity () and the penalty coefficient (). The test network (in Figure 5) has 40 nodes with average degree of 4. We randomly select 20 nodes (50% of total nodes) as the test multicast group, that is, one source node and 19 destination nodes. All OD pairs have the same delay bound for 60 msec.

Firstly, this study conducts 30 runs to find the appropriate partition granularity () for this problem. When the penalty coefficient is fixed (), the experimental results for different partition granularities changing from , , to are shown in Table 2. We normalize these results relative to that of and show the percentage comparison in Figure 6. Obviously, the RPGA with 6 partitions can achieve better results than other partition settings with respect to all the three performance metrics.

Table 2: The results obtained by different partition granularities () when .
Figure 6: Comparison with different partition granularities (in percentage relative to ).

Secondly, when the partition granularity is given (), the experiments on different penalty coefficients ( = 1.0001, 1.001, 1.01, 1.1, and 10) are tested for 30 runs. The experimental results on the average of cost, success rate and computing time are depicted in Table 3. The normalized percentages relative to are shown in Figure 7. Noticeably, the RPGA with can achieve better results than other settings on all the cost, success rate, and computing time metrics.

Table 3: The results obtained by different penalty coefficients () when .
Figure 7: Comparison with different penalty coefficients (in percentage relative to ).
5.4. Comparison with Other Existing Methods

The performance of the proposed RPGA is compared with two well-known penalty methods, which are the Wang’s penalty (WP) method [22] and the DP method [9] for three kinds of multicast scenario. In the first test scenario, the test network has 40 nodes with average degree of 4, which is the same as Figure 5. The multicast group includes 20 OD pairs. All OD pairs have the same delay bound for 60 msec. We execute each method for 30 independent runs and report the experimental results in Table 4.

Table 4: The results of different penalty methods.

For the mean cost in Table 4, the RPGA can find the best result in average 30 runs. Furthermore, the standard deviation of the routing cost obtained by the RPGA is smaller than that by the DP and similar to that of the WP. That is, the RPGA can reliably find the minimum-cost Steiner tree. Compared with the success rate, the RPGA is the best method with respect to the mean and standard deviation of the success rate, even though the WP has a high probability of converging on infeasible solutions. That is, the proposed RPGA can succeed in finding feasible and minimum-cost solutions. From the CPU-time metric, the comparison can help us to realize how much time complexity is needed to compute the RP method because the RPGA is enhanced from the DP method. In Table 4, the RPGA only spends a little computing effort (about extra 6.7% computing time) to obtain better results than the DP method.

For comparison, all the results are normalized as percentages relative to the results of the RPGA (in Figure 8). Obviously, the performance comparison shows that the RPGA can find the minimum-cost multicast tree effectively (with higher success rate) and efficiently (with lower computing effort) than the other two methods.

Figure 8: Comparison with different penalty methods (in percentage relative to the RPGA).

The evolution curves of all three methods on the total cost and path delay are shown in Figures 9(a) and 9(b), respectively. In Figure 9(a), we can observe that both the DP and RPGA methods rapidly converge on the low-cost results after 8 generations; however, the WP method takes 13 generations to slowly converge on a high-cost one. From the path delay aspect in Figure 9(b), all the three methods can satisfy the delay bound. The RPGA, DP, and WP methods take 2, 2, and 4 generations to achieve feasible solutions, respectively. Therefore, the proposed RPGA has a similar convergent effect with the DP. Both the RPGA and DP method outperform the WP with respect to the effectiveness and efficiency abilities.

Figure 9: Evolution curves of three methods on the performance metrics: (a) cost and (b) delay.

In the second test scenario, we randomly generate 8 test networks with the numbers of nodes from 10 to 80 to mimic the stress test for these three methods. In those tests, all delay constraint bounds are 60 msec and the multicast group size is 50% of network nodes. When nodes number increases, the network overhead increases obviously and end-to-end delay increases at the same time. The experimental results depicted in Figure 10 are normalized as percentages relative to the result of the RPGA for comparison. Compared with Figures 10(a) and 10(b), we can observe that RPGA can find the minimum cost of feasible multicast tree, even though the success rates of the DP in the 60-node network and the WP in almost all networks are less than 90%. The WP is the worst method in these three methods on the success rate. However, the computing effort of the RPGA is higher than the other two methods in Figure 10(c). We can also find that the more number of network nodes in the test scenarios, the closer computing time needed for all three penalty methods.

Figure 10: Comparison with three methods for different numbers of network nodes in percentage relative to the RPGA on (a) cost, (b) success rate, and (c) CPU time.

In the third test scenario, we change the delay-bound requirements from 20 msec to 90 msec in a test network, which has 40 nodes and its multicast group size is 20. The comparisons between the success rates and the delay bounds are shown in Figure 11. The success rate of the RPGA is similar to that of the DP and is better than that of the WP, especially, when the delay bounds are lower than 60 msec.

Figure 11: Comparison with three methods on the success rates for different delay bounds.

6. Conclusions and Future Works

The proposed RP method cooperates with GAs for dealing with QoS-based multicast routing problems. The principle of the RP method is that the RPGA releases/enforces some penalties on inefficient/efficient constraints during evolution. Importantly, this approach can find the optimum or near-optimal solutions even though the initial population includes infeasible solutions. The performance of the proposed algorithm was measured using three kinds of test scenarios and compared with two state-of-the-art methods. Experimental results indicate that the proposed RPGA can find near-optimal solutions and outperforms two existing methods for constrained multicast routing problems. The proposed algorithm is also robust in obtaining feasible solutions of all the test functions even though the WP method has smaller success rates for some difficult problems. In conclusion, the performance assessment also demonstrates that the proposed RP method has a remarkable capability to balance the objective function and the constraint violations as an effective and efficient method for solving a variety of QoS-based multicast routing problems.

We have observed that the computing effort of the RPGA is higher than that of other two penalty methods in small-scale networks. In the future, we will study the scalability of the proposed RPGA in finding multicast routes for dynamic MANETs with large-scale dimension. Since MANETs allow ubiquitous service access without any fixed infrastructure, developing a distributed algorithm for high mobility environments is also our future work.


  1. M. R. Gareg and D. S. Johnson, Computer and Intractability: A Guide to the Theory of NP-Completeness, W.J. Freeman, New York, NY, USA, 1979.
  2. C.-H. Liu, T.-C. Chiang, and Y.-M. Huang, “A near-optimal multicast scheme for mobile ad hoc networks using a hybrid genetic algorithm,” in Proceedings of the 20th International Conference on Advanced Information Networking and Applications, vol. 1, pp. 465–470, April 2006. View at Publisher · View at Google Scholar · View at Scopus
  3. S. Gangwar, S. Pal, and K. Kumar, “Mobile ad hoc networks: a comparative study of QoS routing protocols,” International Journal of Computer Science & Engineering Technology, vol. 2, no. 1, pp. 771–775, 2012. View at Google Scholar
  4. M. Hamdan and M. E. El-Hawary, “A Novel Genetic Algorithm Searching Approach for Dynamic Constrained Multicast Routing,” in Proceedings of the Canadian Conference on Electrical and Computer Engineering: Toward a Caring and Humane Technology (CCECE '03), pp. 1127–1130, May 2003. View at Scopus
  5. H. T. Tran and R. J. Harris, “Solving QoS multicast routing with genetic algorithms,” in Proceedings of the International Conference on Information Communications and Signal Processing, vol. 3, pp. 1944–1948, 2003.
  6. S. Sumathy and E. Sri Harsha, “Survey of genetic based approach for multicast routing in MANET,” International Journal of Engineering and Technology, vol. 4, no. 6, pp. 474–485, 2013. View at Google Scholar
  7. M. Gen and R. Cheng, “Survey of penalty techniques in genetic algorithms,” in Proceedings of the IEEE International Conference on Evolutionary Computation (ICEC '96), pp. 804–809, May 1996. View at Scopus
  8. A. Homaifar, C. X. Qi, and S. H. Lai, “Constrained optimization via genetic algorithms,” Simulation, vol. 62, no. 4, pp. 242–254, 1994. View at Google Scholar · View at Scopus
  9. J. A. Joines and C. R. Houck, “On the use of non-stationary penalty functions to solve nonlinear constrained optimization problems with GA's,” in Proceedings of the 1st IEEE Conference on Evolutionary Computation, vol. 2, pp. 579–584, June 1994. View at Scopus
  10. A. B. Hadj-Alouane and J. C. Bean, “A genetic algorithm for the multiple-choice integer program,” Operations Research, vol. 45, no. 1, pp. 92–101, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  11. A. T. Haghighat, K. Faez, M. Dehghan, A. Mowlaei, and Y. Ghahremani, “GA-based heuristic algorithms for bandwidth-delay-constrained least-cost multicast routing,” Computer Communications, vol. 27, no. 1, pp. 111–127, 2004. View at Publisher · View at Google Scholar · View at Scopus
  12. M. Zhang, W. Luo, and X. Wang, “Differential evolution with dynamic stochastic selection for constrained optimization,” Information Sciences, vol. 178, no. 15, pp. 3043–3074, 2008. View at Publisher · View at Google Scholar · View at Scopus
  13. C. H. Lin, “A rough penalty genetic algorithm for constrained optimization,” Information Science, vol. 241, pp. 119–137, 2013. View at Google Scholar
  14. D. Powell and M. M. Skolnick, “Using genetic algorithms in engineering design optimization with non-linear constraints,” in Proceedings of the 5th International Conference on Genetic Algorithms, pp. 424–431, 1993.
  15. Z. Michalewicz, “A survey of constraint handling techniques in evolutionary computation methods,” in Proceedings of the 4th Annual Conference on Evolutionary Programming, pp. 135–155, 1995.
  16. Z. Pawlak, “Rough sets,” International Journal of Computer & Information Sciences, vol. 11, no. 5, pp. 341–356, 1982. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  17. Z. Pawlak and A. Skowron, “Rudiments of rough sets,” Information Sciences, vol. 177, no. 1, pp. 3–27, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  18. Z. Pawlak and A. Skowron, “Rough sets and Boolean reasoning,” Information Sciences, vol. 177, no. 1, pp. 41–73, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  19. S. O. Kimbrough, G. J. Koehler, M. Lu, and D. H. Wood, “On a Feasible-Infeasible Two-Population (FI-2Pop) genetic algorithm for constrained optimization: Distance tracing and no free lunch,” European Journal of Operational Research, vol. 190, no. 2, pp. 310–327, 2008. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  20. H. Salama, “The multicast routing simulator,” The Real-Time Communication Project, Version 2, 1997,
  21. B. M. Waxman, “Routing of multipoint connections,” IEEE Journal on Selected Areas in Communications, vol. 6, no. 9, pp. 1617–1622, 1988. View at Publisher · View at Google Scholar · View at Scopus
  22. Z. Wang, B. Shi, and E. Zhao, “Bandwidth-delay-constrained least-cost multicast routing based on heuristic genetic algorithm,” Computer Communications, vol. 24, no. 7-8, pp. 685–692, 2001. View at Publisher · View at Google Scholar · View at Scopus