Grey Wolf Algorithm and Multi-Objective Model for the Manycast RSA Problem in EONs

: Manycast routing and spectrum assignment (RSA) in elastic optical networks (EONs) has become a hot research ﬁeld. In this paper, the mathematical model and high efﬁcient algorithm to solve this challenging problem in EONs is investigated. First, a multi-objective optimization model, which minimizes network power consumption, the total occupied spectrum, and the maximum index of used frequency spectrum, is established. To handle this multi-objective optimization model, we integrate these three objectives into one by using a weighted sum strategy. To make the population distributed on the search domain uniformly, a uniform design method was developed. Based on this, an improved grey wolf optimization method (IGWO), which was inspired by PSO (Particle Swarm Optimization, PSO) and DE (Differential Evolution, DE), is proposed to solve the maximum model efﬁciently. To demonstrate high performance of the designed algorithm, a series of experiments are conducted using several different experimental scenes. Experimental results indicate that the proposed algorithm can obtain better results than the compared algorithm.


Introduction
The rapid growth of numerous high-rate various applications, such as internet protocol television, video on demand, and cloud computing, requires an efficient networking infrastructure; the future optical network tends to be dynamic, heterogeneous, and unpredictable [1,2]. To tackle this issue, elastic optical networks (EONs) [3] are proposed to realize flexible and efficient spectrum allocation with much finer spectrum granularity. In particular, EONs provide just-right bandwidth for the arriving request connection dynamically, which brings better spectrum assignment flexibility [4,5]. Orthogonal frequency division multiplexing is a multi-carrier modulation technology. It can distribute the high-speed data stream into several orthogonal low-speed subcarriers [6]. The adjacent subcarriers have the spectrum overlapping of a subcarrier bandwidth. This subcarrier is referred to the frequency slot (FS). The elastic optical network can allocate several consecutive frequency slots to each connection request according to the required bandwidth by using orthogonal frequency division multiplexing as a spectrum-efficient modulation technology. Adjacent spectrum bandwidths assigned to two connection requests in the same link should be separated by the guaranteed frequency slots (GFs). Similar to the routing and wavelength assignment problem in wavelength division multiplexing networks [7], a routing and spectrum assignment (RSA) problem exists in the elastic optical network [8]. In order to establish a light-path for the connection request in the elastic optical network, three constraints should be satisfied as follows: (1) Spectrum consistency means that the start frequency slot index on different links of a path must be identical; (2) Spectrum continuity means that we must assign consecutive frequency slots to a specific connection request. That is to say, a large connection request can not be divided into several smaller connection requests; (3) A frequency slot on a link should be assigned to one connection request at most. Generally speaking, the objective of static routing and spectrum assignment is to minimize the the maximum index of the used frequency slots with unlimited resources, and to minimize the ratio of blocking with limited resources. Certainly, there are some other optimization objectives, such as energy consumption, cost, etc.
In this paper, we focus on the manycast routing and spectrum assignment problem in EONs. The main contributions of this paper are as follows. First, to minimize network power consumption, the total occupied spectrum and maximum index of the used frequency spectrum, we established a multi-objectives optimization model. Second, we first integrate these three objectives into one by using a weighted sum strategy to handle this multi-objective optimization model. Then, an improved grey wolf optimization method (IGWO), which was inspired by PSO (Particle Swarm Optimization, PSO) and DE (Differential Evolution, DE) is proposed to solve the maximum model efficiently. Finally, a series of experiments are conducted in several different experimental scenes.
The rest of this paper is organized as follows. Some related works are introduced in Section 2. Section 3 gives the network architecture and the optimization model. To solve the optimization model effectively, we propose an improved grey wolf optimization algorithm in Section 4. To evaluate the algorithm proposed, simulation experiments are conducted, and the experimental results are analyzed in Section 5. The paper is concluded with a summary in Section 6.

Related Works
To detect potential access conflicts and prevent both processes from updating data simultaneously in distributed database systems, the concept of manycast has been proposed first [9,10]. The problem of manycasting over optical burst-switched networks has been investigated [11][12][13]. The main challenge is providing reliability despite random contentions optical burst-switched networks for the problem of manycasting. This research focuses on distributed routing or unicast routing algorithms to provide reliable manycast for dynamic traffic. Some literature is focusing on the problem of manycast routing and wavelength assignment in wavelength division multiplexing networks (WDM) [14][15][16]. However, this research is focusing on the manycast routing and wavelength assignment in wavelength division multiplexing networks. Thus, these algorithms can not work on the elastic optical networks well. There are some different network properties among optical burst-switched networks, wavelength division multiplexing networks, and elastic optical networks; the manycast routing and spectrum assignment problem in elastic optical networks mainly considers efficient network resource utilization and request blocking probability reduction. The literature [16] investigated the manycast routing and spectrum assignment (MRSA) problem in WDM networks. The proposed heuristics observably improved network performance in required wavelengths reduction over realistic networks. However, some factors have not been considered in EONs, such as modulation level and spectrum assignment constraints. Because of unique spectrum flexibility in EONs, it has an essential difference compared to supporting manycast with WDM networks and optical burst-switched networks. Thus far, there have only been a few studies about the MA-RMLSA problem in EONs. The energy-efficiency MA-RMLSA strategy was proposed by green-energy aware destination nodes selection [17]. The proposed algorithm has a high performance on decreasing the energy consumption of the network. While the authors focus on the network energy consumption, it may not have advantages over some other objectives, such as spectrum resource utilization, maximum index of used frequency slots, etc. Impairment-aware manycast routing, modulation level, and spectrum assignment problem in EONs are investigated. Two decomposed MILPs (Mixed Integer Linear Programming) and corresponding heuristic algorithms are proposed to find a light-tree and assign modulation level and spectrum to the given requests, sequentially [18]. The authors in [19] studied an integrated approach to optimally place content replicas across DCs (Data Centers) by concurrently solving the routing and wavelength assignment (RWA) problem for both inter-DC content replication and synchronization traffic following the manycast routing paradigm, and end-user-driven user-to-DC communication following the anycast routing paradigm, with the objective to reduce the overall network capacity usage.

Problem Description and Mathematical Modeling
In this section, problem description and mathematical modeling of the manycast routing and spectrum assignment (RSA) problem in elastic optical networks (EONs) will be given.

Problem Description
We use an undirected graph V = (V, E) to denote a network, where V = {V 1 , V 2 , · · · , V N V } and N V denote the nodes set and the number of the nodes in the network, respectively. E = {l ij |V i , V j ∈ V} denotes the link set, and N E denotes the number of links in the network. If l ij = l ji = 1, there is a link between V i and V j ; otherwise, l ij = l ji = 0. Let f = { f 1 , f 2 , · · · , f u , · · · , f N F } denote the set of available frequency slots (FSs) in each link, and N F be the number of frequency slots. R = {R 1 , R 2 , · · · , R k , · · · , R N R } denotes a set of connection requests, where N R is the number of connection requests, and R k is the k-th connection request. R k can be represented as R k = (s k , D k , B k ), where s k and B k represent the source node and the numbers of frequency slots of R k required.
} is the set of destination node, where N k D is the number of destination node; when N k D = 1, it will be a unicast routing request. In this paper, we assume that all nodes in the network are able to split their incoming connection request to any number of other nodes. This architecture is the same as the scheme introduced in [16].
The manycast routing and spectrum assignment (RSA) problem in elastic optical networks (EONs) can be summarized as: to achieve some objectives, the proper path should be selected for each connection request. Then, the optimal scheme of spectrum assignment should be determined for all of the connection requests.

Mathematical Modeling
In this section, we present a mathematical model for manycast routing and spectrum assignment in EONs. The first objective is minimize the total power consumption, and the total power consumption is calculated by [17] where Q ij is a boolean variable; Q ij = 1 if and only if link l ij is used in the network to provision all manycast requests. A ij and P OA denote the number of optical amplifiers in the link l ij and the power consumption of each optical amplifier, respectively. P j is a boolean variable; P j = 1 if and only if node V j is used in the network to provision all manycast requests. P OXC (V j ) is the power consumption of optical cross connect in node V j . Since Equations (2) and (3) satisfy the total power consumption can be normalized as Thus, we have 0 ≤ F 1 ≤ 1, and the first objective function can be expressed as The second objective is minimize the total occupied frequency slots, and the total occupied spectrum is where y ku ij represents a binary valuable; y ku ij = 1 if and only if the frequency slot f u on link l ij is allocated to connection request R k . We have F S ≤ N E × N F , and the total occupied frequency slots can be normalized as Thus, 0 ≤ F 2 ≤ 1, and the second objective function can be represented by The third objective is minimize the maximum index of used frequency slots (MIUFS) in the network; we can express this objective function as where n(F ij ) denotes the maximum index of used frequency slots on link l ij . Since we have To simplify the model, we integrate the three objectives into one to be minimized by the sum weighted strategy as follows: where α 1 , α 2 , and α 3 are three weights to adjust the importance of the three objectives, and we have The objective should be made under some conditions. These conditions constitute the constraints of the problem as follows: Constraint (a): the same spectrum slots are not assigned to two requests. That is, Constraint (b): contiguous frequency slots should be allocated to the connection request, where x ku ij represents a binary valuable, x ku ij = 1 if and only if f u is the first frequency slot on link l ij allocated to connection request R k . Constraint (c): all the destinations must be reached, we can express this constraint as Constraint (d): at least one outgoing traffic should leave the source node.
Constraint (e): it ensures that the source node could not have incoming traffic: Constraint (f): each node, except the source node, can have at most one piece of incoming traffic.
Constraint (g): one node, except the source node, could not have outgoing traffic unless it has incoming traffic, Constraint (h): if a node is not one of the target destinations and has incoming traffic, it should have one or more pieces of outgoing traffic, Based on the objectives and constraints above, we can set up a global constrained optimization model as follows: The problem of manycast routing and spectrum assignment in EONs is the hardest combinatorial optimization problems. The existing algorithms cannot be applied directly, and are necessary to make some improvements or revisions. To solve the global constrained optimization model established, we propose an improved grey wolf optimization method and denote it as IGWO.

Grey Wolf Optimization (GWO)
The Grey Wolf Optimization (GWO) algorithm simulates the leadership hierarchy and hunting mechanism of grey wolves [20] and has been proven to be an effective technique for many hard problems [21][22][23]. However, it is not suitable to directly apply the algorithms mentioned above to the problems of manycast routing and spectrum assignment in EONs, and it is necessary to make some improvements or revisions on them. In this section, we will describe the proposed IWGO detailed.

Encoding Scheme
In the manycast routing and spectrum assignment problem, we should determine the optimal scheme of routing and spectrum assignment. For the spectrum assignment, it is much easier to assign spectra using first fit strategy [6] than using the method with encoding. Thus, it only needs to encode for routing scheme.
Each individual in routing population represents a routing scheme for all the connection requests.
k } denotes the candidate paths set of connection request R k that is calculated by the K-Shortest path algorithm in advance, where N Q k is the number of the candidate paths and Q q k is the q-th path. We assume that y = (y 1 , y 2 , · · · , y N R ) is an individual in path selection population. y k = q if and only if R k occupies the path Q q k .

Population Initialization
In the proposed improved Grey Wolf Optimization (IWGO) algorithm, we use uniform design method to generate the population. To generate points to be uniformly distributed on the experimental domain, a uniform design method was developed. It generates a small number of the uniformly distributed representative points in a domain by using a uniform array U(S, H) = [U i,j ] H×S , where U i,j denotes the level of the j-th factor in the i-th combination with the j-th factor representing the j-th variable and its level being its value [24,25].
To construct uniform design array, many methods are presented-not only simple but also efficient methods are proposed. Firstly, we construct a hypercube over an S-dimensional space: where a i and b i are the lower and upper bounds of the i-th factor (i.e., i-th variable), respectively. Then, a hyper-rectangle is formed between a i and d i as follows: Finally, H uniformly distributed points are selected randomly from C S . Assume that H(d) is the number of points fallen into the hyper-rectangle C(d), and the fraction of points in C(d) is H(d)/H. As the volume of hypercube Hence, we can map these H points in C S to the problem domain with S factors and χ levels uniformly, where H is an odd and H > S. It has been proved that U i,j can be given by [26]: where σ is a constant related to the number of factors S and level χ. The H sample points scattered uniformly in the hypercube can be selected.

Improved Grey Wolf Optimization (IGWO) for Manycast RSA
In the Grey Wolf Optimization method, each grey wolf denotes an individual. Four types of grey wolves β, γ, δ, and ω denote the optimal individual, suboptimum individual, third-optimum individual, and other individuals. Assume that there are N I grey wolves, and the position of i(i = 1, 2, · · · , N I )-th wolf can be denoted as We can update the position of i-th wolf by where z p = (z 1 i , z 2 i , · · · , z d i , · · · , z D i ) and t denote the position vector of the prey and the current iteration, respectively. A and C denote coefficient vectors, and are calculated as follows: where r 1 , r 2 are two random vectors in [0, 1]; a linearly decreases from 2 to 0 during the course of iterations. The position of other wolves can be updated according to the position of individual β, γ, δ (denoted as x β , x γ and x δ ) To enhance the search ability and increase the convergent speed, an efficient position update method of individual is proposed as follows: where r 3 , r 4 and r 5 are three random vectors in [0, 1]; x i,best denotes the best position of i-th individual in the past. x j (j = i) represents a random individual in current iteration. X d i (t) can be calculated by where µ(1 ≤ µ ≤ t − 1) is a constant. In this position update method of an individual, the position of the other individual and its past position information are used like PSO and DE. Thus, it can enhance the search ability and increase the convergent speed.
In the encoding scheme, each gene in all individuals is a positive integer. However, some real numbers can be obtained with Equation (26). For this situation, we only use the integer portion as the gene of the individual. In addition, some gene value, which is greater than the upper bound, can be obtained. For this situation, we only use the integer portion as the gene of the individual. The gene value is modulus the upper bound to make the individual update to a feasible solution. Through these two methods, an infeasible solution can be modified as a feasible solution.

Framework of the IGWO
To make the proposed improved grey wolf optimization algorithm understood clearly, we give the framework of the proposed algorithm in Algorithm 1. In the algorithm, step 1 is to initialize the population according to the uniform design method. It can improve the search ability of the algorithm. In step 2, fitness is calculated for all the wolves in the population by using the fitness function (objective function is defined as fitness function in this work). Step 4 to step 9 is update the position of all the wolves by using Equation (23), Equation (24), and Equation (26). All the infeasible solutions are converted to feasible solutions in Step 10.
Step 11 is calculate the fitness for all the wolves after the position updated. The position of optimal individual, suboptimum individual, third-optimum individual (x β , x γ , and x δ ) are updated in step 12. for i = 1, 2, · · · , N I do 5 for d = 1, 2, · · · , D do 6 A and C are calculated according to Equations (23) and (24);

Experiments and Analysis
To demonstrate the effectiveness and efficiency of the proposed algorithm, several experiments are conducted, and the results are presented in this section. In Section 5.1, the parameters used in the algorithms are given. Experimental results are presented in Section 5.2. Finally, the experimental analysis is given in Section 5.3.

Parameters Setting
In the experiments, two widely used networks are used as shown in Figures 1 and 2, i.e., NSFNET (National Science Foundation Network) with 14 nodes and 21 links and US Backbone (United States Backbone) with 27 nodes and 44 links, respectively [27,28], respectively, in Figures 1 and 2; each number on the link denotes the distance between adjacent nodes, and the unit of the link distance is Km. We assume that FSs is 12.5 GHz, and transmission distance of BPSK (Binary Phase Shift Keying), QPSK (Quadrature Phase Shift Keying), 8QAM (8 Quadrature Amplitude Modulation), 16QAM (16 Quadrature Amplitude Modulation) are 9600, 4800, 2400 and 1200 km, respectively. Four groups' connection requests are generated, and their numbers are 250, 500, 750 and 1000, respectively. All connection requests in every group satisfy uniform distribution among all nodes in two topologies. In our work, a large number of experiments are conducted for each case. To make the algorithm converge to an optimal solutions, t max = 2000 and µ = 3 are adopted. Generally speaking, when the population size is large, it will require a long computation time. In addition, when the population size is small, it will result in a bad diversity of population. Thus, N I = 100 is selected in many research works. Like existing works, we use N I = 100 in the algorithm. Each connection request requires frequency slots that satisfy uniform distribution in [1,10], and each link has 1000 frequency slots, i.e., N F = 1000.

Experimental Results
To demonstrate the performance of the proposed algorithm, we compare the proposed algorithm IGWO with other three algorithms. The first is EEM, which was proposed in [17]. Another one was proposed in [29], and denoted as RSAGA. RSAGA optimizes the MA-RMLSA problem in routing constitution, modulation level allocation, and spectrum assignment jointly, to enhance the performance of the network. In addition, we also compared IGWO with GWO (Grey Wolf Optimization), which was proposed in [20].
To demonstrate the performance of proposed model and algorithm, we design two experimental scenes. In the first scene, we fixed the number of destination nodes as N D = N V /3 and N D = 2N V /3, i.e., N d = N D /N V = 1/3 and N d = N D /N V = 2/3. That is to say, the number of destination nodes is generated in [N V /6, N V /3] and [N V /3, 2N V /3] randomly. Figure 3 shows the results obtained in NSFNST topology and US Backbone topology when α 1 = 1, α 2 = 0, α 3 = 0. Figure 4 shows the results obtained in NSFNST topology and US Backbone topology when α 1 = 0, α 2 = 1, α 3 = 0. Figure 5 shows the results obtained in NSFNST topology and US Backbone topology when Figure 6 shows the results obtained in NSFNST topology and US Backbone topology when α 1 = α 2 = α 3 = 1/3. In each experiment, the number of connection requests are set as    In the second scene, we fixed the three weights α 1 , α 2 and α 3 as α 1 = 1/3, α 2 1/3 and α 3 = 1/3, i.e., the objective function is min

Experimental Analysis
In the first experimental scene, the experimental results are obtained by the proposed algorithm (IGWO) and three compared algorithms (EEM, RSAGA, and GWO) are shown in Figures 3-6. In Figure 5, the experimental results are obtained in NSFNET and US BackBone when α 1 , α 2 and α 3 are selected as 1, 0, and 0, respectively. Thus, the objective function is to minimize total power consumption. From the experimental results, we can see that the IGWO can obtain better results than the three compared algorithms. The total power consumption obtained by the IGWO is 2.8%-4.9% less than those obtained by EEM, RSAGA, and GWO when the number of connection requests is 0.25 N V (N V − 1). When the number of connection requests is 4 N V (N V − 1), the total power consumption obtained by IGWO is 6.9%-11.6% less than those obtained by EEM, RSAGA, and GWO, respectively. That is to say, the IGWO can obtain a smaller total power consumption and save more power used than EEM, RSAGA, and GWO with the increase of the number of connection requests. The RSAGA algorithm is to minimize the maximum index of used frequency slots, and the EEM algorithm can decrease the energy consumption. Thus, the total power consumption obtained by EEM is less than that obtained by RSAGA. The proposed algorithm IGWO uses uniform design to generate initial population; it can enhance the search ability of the algorithm. Thus, it can obtain the best scheme among the three algorithms. From the experimental results, we can see that the total power consumption obtained when the number of destination node is N DC = 2N V /3 is larger than that obtained when the number of destination node is N DC = N V /3 for the same number of connection requests. With the increase of destination node, it will increase the number of connection requests. Thus, the total power consumption is increased.
When α 1 , α 2 and α 3 are selected as 0, 1, and 0, the objective function is to minimize the total occupied frequency slots. The experimental results obtained in two networks are shown in Figure 4 with the different number of connection requests. From the experimental results, we can see that the IGWO can obtain better results than the three compared algorithms. The total occupied frequency slots obtained by the IGWO is 4.2%-6.2% less than those obtained by EEM, RSAGA, and GWO when the number of connection requests is 0.25 N V (N V − 1). When the number of connection requests is 4 N V (N V − 1), the total power consumption obtained by IGWO is 8.6%-10.7% less than those obtained by EEM, RSAGA, and GWO, respectively. That is to say, the IGWO can obtain a smaller total occupied frequency slots and save more frequency slots used than EEM, RSAGA, and GWO with the increase in the number of connection requests. The proposed algorithm IGWO uses uniform design to generate initial population; it can enhance the search ability of the algorithm. Thus, it can obtain the best scheme among the three algorithms. From the experimental results, we can see that the total occupied frequency slots obtained when the number of destination node is N DC = 2N V /3 is larger than that obtained when the number of destination node is N DC = N V /3 for the same number of connection requests. With the increase of destination node, it will increase the number of connection requests. Thus, the total occupied frequency slots is increased.
In Figure 5, the experimental results are obtained in NSFNET and US BackBone when α 1 , α 2 and α 3 are selected as 0, 0, and 1, respectively. Thus, the objective function minimizes the maximum index of used frequency slots. From the experimental results, we can see that the IGWO can obtain better results than the three compared algorithms. The total occupied frequency slots obtained by the IGWO are 3.8%-5.9% less than those obtained by EEM, RSAGA, and GWO when the number of connection requests is 0.25 N V (N V − 1). When the number of connection requests is 4 N V (N V − 1), the maximum index of used frequency slots obtained by IGWO is 8.1%-11.2% less than those obtained by EEM, RSAGA, and GWO, respectively. That is to say, the IGWO can obtain a smaller maximum index of used frequency slots than EEM, RSAGA, and GWO with the increase of the number of connection requests. The proposed algorithm IGWO uses uniform design to generate initial population and well-designed strategy of position update; it can enhance the search ability of the algorithm. Thus, it can obtain the best scheme among the three algorithms. From the experimental results, we can see that the maximum index of used frequency slots obtained when the number of destination node is N DC = 2N V /3 is larger than that obtained when the number of destination node is N DC = N V /3 for the same number of connection requests. With the increase of destination node, it will increase the number of connection requests. Thus, the maximum index of used frequency slots is increased. Figure 6 shows that the experimental results obtained in NSFNET and US BackBone when α 1 , α 2 and α 3 are selected as 1, 1, and 1, respectively. Similarly, we also can see that the IGWO can obtain better results than EEM, RSAGA, and GWO. The objective function obtained by the IGWO is 3.7%-5.8% less than those obtained by EEM, RSAGA, and GWO when the number of connection requests is 0.25 N V (N V − 1). When the number of connection requests is 4 N V (N V − 1), the objective function obtained by IGWO is 8.2%-11.6% less than those obtained by EEM, RSAGA, and GWO, respectively. The IGWO can obtain a smaller objective function than EEM, RSAGA, and GWO with the increase of the number of connection requests. The proposed algorithm IGWO uses uniform design to generate initial population and a well-designed strategy of position update; it can enhance the search ability of the algorithm. Thus, it can obtain the best scheme among the three algorithms. From the experimental results, we can see the objective function obtained when the number of destination nodes is N DC = 2N V /3, which is larger than that obtained when the number of destination node is N DC = N V /3 for the same number of connection requests. With the increase of destination node, it will increase the number of connection requests. Thus, the objective function is increased.
In the second experimental scene, the experimental results obtained by the proposed algorithm (IGWO) and three compared algorithms (EEM, RSAGA, and GWO) are shown in Figures 7-11. In this experimental scene, we set α 1 , α 2 and α 3 as 1, 1 and 1, respectively. From the experimental results, we can see that IGWO can obtain better results than EEM, RSAGA, and GWO for different connection requests in two networks. In each figure, the objective function is increased with the increase in the number of destination nodes. When the number of destination nodes is N D = 0.2N V , the objective function obtained by IGWO is 3.5%-5.8% less than those obtained by EEM, RSAGA, and GWO, respectively. The objective function obtained by the IGWO is 8.6%-11.7% less than those obtained by EEM, RSAGA, and GWO when the number of connection requests is N D = N V . With the increase of destination node, it will increase the number of connection requests. Thus, the objective function is increased. In addition, EEM minimizes the total power consumption, and RSAGA minimizes maximum index of used frequency slots, so EEM, RSAGA, and GWO cannot be distinguished when α 1 , α 2 and α 3 as 1, 1, and 1.
As shown in the experimental results, we can see that IWGO can obtain better results than that obtained by GWO. In the IGWO, we improved the strategy of the position update method for the individual. The position of another individual and its past position information are used like PSO and DE. Thus, it can enhance the search ability and increase the convergent speed. In addition, the parameter µ is used. It can help to take advantage of the trajectory information of the individuals. When µ = 1, the proposed IGWO degraded to the standard GWO algorithm. Position update method can use the past µ position information when µ ≥ 1. Thus, IGWO is better than GWO for this optimization problem. That is to say, IGWO can obtain a better solution than GWO.

Conclusions
In this paper, we investigate the manycast RSA problem in EONs. A multi-objective optimization model, which minimizes network power consumption, the total occupied spectrum, and maximum index of used frequency spectrum, is established. To solve this multi-objective optimization model, we integrate these three objectives into one by using a weighted sum strategy. Then, an improved grey wolf optimization method (IGWO) is proposed. To demonstrate high performance of the designed algorithm, a series of experiments are conducted in several different experimental scenes. Experimental results indicate that the proposed algorithm can obtain better results than the compared algorithm. According to the experimental results, we can find that the objective function obtained by the proposed algorithm is 3%-12% less than those obtained by compared algorithms in different networks. However, integrating these objectives into one by using the weighted sum method also has some disadvantages. Thus, we will investigate the efficient algorithm based on the multi-objective optimization algorithm, such as MOEA/D, or other swarm intelligent algorithm. In addition, there are other objectives for the manycast RSA Problem in EONs. The multi-objective optimization algorithm will be investigated to obtain the Pareto front. Thus, it can give more decision-making plans to decision-makers.