Hypergraph-based resource allocation for Device-to-Device underlay H-CRAN network

In the hybrid communication scenario of the Heterogeneous Cloud Radio Access Network and Device-to-Device in 5G, spectrum efficiency promotion and the interference controlling caused by spectrum reuse are still challenges. In this article, a novel resource management method, consisting of power and channel allocation, is proposed to solve this problem. An optimization model to maximum the system throughput and spectrum efficiency of the system, which is constrained by Signal to Interference plus Noise Ratio requirements of all users in diverse layers, is established. To solve the non-convex mixed integer nonlinear optimization problem, the optimization model is decomposed into two sub-problems, which are all solvable quasi-convex power allocation and non-convex channel allocation. The first step is to solve a power allocation problem based on solid geometric programming with the vertex search method. Then, a channel allocation constructed by three-dimensional hypergraph matching is established, and the best result of this problem is obtained by a heuristic greed algorithm based on the bipartite conflict graph and µ-claw search. Finally, the simulation results show that the proposed scheme improves the throughput performance at least 6% over other algorithms.


Introduction
Recently, as more and more mobile devices are included in the wireless network, the demand for the volume of the access point and the system throughput has a rapid growth. 1,2 By 2025, the total number of the wireless devices will surpass 75 billion, and the corresponding requirement of data is over 3.3 ZB. 3,4 In order to meet the requirement of access capacity and system capacity, the Heterogeneous Cloud Radio Access Network (H-CRAN), combined with Device-to-Device (D2D) communication, is introduced into the 5G core network.
In traditional wireless network structure, to guarantee the seamless coverage, the high-power node integrated control and forwarding unit is utilized. However, it is difficult to fulfill the throughput and access point volume requirement at the hot area. Thus, in order to increase the capacity of the wireless access network effectively, the H-CRAN network is introduced, which adds centralized baseband processing unit (BBU) entities and distributed remote radio heads (RRHs) components while retaining the original highpower node. [5][6][7] For meeting the massive data needs of users in hot areas, RRHs as signal forwarding relay are ultra-densely deployed. All BBU units with control plane are grouped together in a centralized location to form a BBU pool aggregating all resources within cloud computing. [8][9][10] Then, BBU pool can not only realize large-scale radio resource allocation but also offer collaborative processing. [11][12][13] BBU controls RRH through fronthaul link while BBU pool communication with traditional high-power node via backhaul links. By this way, the networking benefits are fully utilized without increasing the running cost and the energy consumption of network infrastructure.
D2D communication is defined as user equipments (UEs) in physical proximity that directly communicate with each other without passing through cellular infrastructures, which greatly reduce the burden of cellular infrastructure. 14,15 To efficiently enhance spectral efficiency, in general, D2D communications are considered to utilize the same spectrum as cellular user (CU). Thus, D2D communications have the great potential advantage of high data rate, short delay, low power consumption, and increased network capacity for supporting popular proximity-based applications. 16 Combining D2D communication with H-CRAN provides larger network capacity and better user experience. Due to the limited spectrum resources of cellular networks, the user in high-density mutual communicating RRHs environment with a resource-reused manner as D2D communication, which will bring unpredictable inter-layer interference or co-layer interference. 17,18 This kind of interference will not only affect the quality of service (QoS) of users but also lead to system performance degradation. Therefore, it is of great importance to make reasonable resource allocation.

Related works
To fully boost advantages brought by D2D and H-CRAN, lots of research efforts have been addressed by industry and academia. Based on stochastic geometry framework, Abana et al. 19 analyzed the improvement of performance in terms of coverage and average traffic delivery rate with non-uniform D2D deployment and uniform D2D deployment in H-CRAN. With the same way as Abana et al., 19 the spectral efficiency performance of system is integrated in Liu et al. 20 for D2D communications with downlink coordinated multipoint (CoMP) in C-RANs, whose mode selection scheme is depended on distance. However, Abana et al. 19 and Liu et al. 20 analyze the influence of D2D communication on the performance of H-CRAN network, but ignore the role of D2D performance itself.
Instead, an improved directed-hypergraph-based local altruistic game is proposed in Sun et al., 21 which is aiming to maximize D2D access rate and channel reuse rate as well as minimizing the uplink transmission power in multi-cell scenarios. But it only focuses on D2D's own performance. Expecting to improve the spectrum efficient, the previous studies 22-24 pay more attention on energy efficient. A resource allocation based on distributed non-cooperative game theory is investigated in Zhou et al. 22 Based on D2D in-band and out-band, a novel energy efficiency spectrum sharing mechanism is proposed in Li et al. 23 and Gui et al. 24 However, their optimized object is the edge user of H-CRAN network. All the above-mentioned works take two-layer hybrid network structure into account, while ignoring the overall system performance.
In order to achieve the maximal overall system throughput, the previous studies 25, 26 considered the downlink reuse scenarios and exploited a many-to-one matching channel allocation to optimize underlaying user access rates, where power allocation is forgotten. Different from above, a joint mode selection and channel allocation scheme are researched in Zhen and Sun 27 to maximize the system uplink throughput, by using game theory on distributed coalition formation. Based on Sun et al., 28 a distributed and staged resource allocation mechanism is proposed, which considered the mode selection, channel allocation, as well as the power control. Furthermore, a joint mode selection, user connection, and power allocation dynamic scheme are investigated in Mo et al., 29 which takes the time-varying of traffic and channel into account. However Sun et al. 28 and Mo et al. 29 still considered only D2D users, neglecting the performance improvement of other users.
From the above articles, it is obviously seen that fully utilizing the advantages of D2D communication could effectively improve the coverage, spectrum efficiency, power efficiency, throughput, and other performance of the H-CRAN network. However, most of the research are interested in enhancing the performance of underlaying network such as RRH users or D2D user. Only a few works concerned with the whole system situation. Some of the literature which optimizes global performance only show solicitude for channel allocation and abandons the role of power control. Therefore, a joint power allocation and channel selection scheme is quite essential.

Main contributions
In this article, to maximize overall H-CRAN throughput and guarantee the QoS requirements of each H-CRAN layer, jointly optimizing the performance of overall H-CRAN with power allocation and channel allocation is considered. This resource allocation is formulated into a nonlinear integer programming problem with constraints on maximum power, Signal to Interference plus Noise Ratio (SINR), and channel number. To deal with this optimization problem, a two-step optimization framework is utilized to decompose the original problem into two sub-problems. The main contributions are as follows.
A joint optimization resource allocation is proposed in this article, which not only focuses on the power allocation but also takes the channel allocation into account. To better adapt to the H-CRAN, a semi-centralized scheme is utilized, which is effectively suitable for the distinct processing characteristics of BBU pool. Different from binary mapping model, the 3D mapping-based solid geometric programming and hypergraph is introduced. The solid geometric programming is used to map power and its constraints. The hypergraph is investigated to capture the complex multiple co-channel relationship. Considering the complexity, the non-convex mixed integer programming original problem is decomposed into two sub-problems, the power allocation and channel allocation. Then, based on 3D mapping-based model, the Geometric Vertex Search approach and bipartite conflict graph are used to solve power and channel allocation problem. The simulation results validate that the proposed algorithm has a significant gain over the others in terms of throughput and user access rate. Furthermore, the proposed algorithm is easier to implement based on BBU pool.
The contents of this article are as follows. In section ''System model,'' the system model and optimization problem are explained. In section ''Problem decomposition and algorithm design,'' the solid geometric programming and hypergraph matching model are proposed, correspondingly Geometric Vertex Search approach and 3D Matching Game are elaborated to search the optimal power allocation and channel allocation. Section ''Numerical results'' is numerical results of the proposed algorithm and the performance of system. The conclusion is revealed in section ''Conclusion.''

System model
As illustrated in Figure 1, an H-CRAN with D2D communication system contains M macro-cell users (MUEs), which is M 2 fMUE m m = 1, . . . , M j g. The BBU pools equipped with a high-power base station (e.g. eNB) is the center unit of the network. K D2D pairs are located among the wireless environment, where K 2 fD k k = 1, . . . , K j g. Each D2D pair is composed of a transmitter D kT and a receiver D kR . In addition, R set of RRHs are equipped to provided service, where R 2 fRRH r j, r = 1, . . . , Rg. The R-UEs in each RRH are denoted as N r 2 fRUE r j , j = 1, . . . , N r g. The total number of all RUEs in H-CRAN is denoted as N = P r2R N r , and whose set is represented as N 2 RUE n jn = 1, . . . , N f g . For simplicity, we divide the uplink resources into equal sized sub-channels and exploit L 2 fCh l l j = 1, . . . , Lg to denote the set of all subchannels. To protect the communication requirements of MUE, one MUE occupied one sub-channel (i.e. MUE m occupied one sub-channel l 2 L). Considering the limited spectrum resources, the RRH users and D2D pairs communicate by multiplexing the channel resources of MUE users. Besides, to prevent high interference, one channel is limited to serve at most one RRH user as well as one D2D pair, while one RRH user at one sub-channel shares with at most one D2D pair. For example, D2D-2 and RUE-3 reuse the uplink resource of MUE-3, so the D kR suffer interference from MUE-3 and RUE-3. In the same time, eNB received interference signal from RUE-3 and D kT and RRH-3 received interference signal from RRH and D kT . To prevent excessive interference, each sub-channel is occupied by one MUEs, one RUEs, and one DUEs. Therefore, interference cannot occur between users in the same layer (i.e. no interference exists between D2D pairs as same as RUEs).
For macrocell, the link gain of the mth MUE is denoted as h m , and g N n, m is denoted as the interference gain from the nth RUE and g K k, m is denoted as the interference gain from the kth D2D pair. For RRHs, the link gain of the jth RUE in RRH-r is denoted as h r, j n . For convenience, let h n to denote the link gain of the nth RUE in overall system, for example, h n = h r, j n . Similarly, g M m, n is denoted as the interference gain from the mth MUE, g K k, n is denoted as the interference gain from the kth D2D pair. For kth D2D pair, the link gain is denoted as h k . Let g N n, k represent the interference gain from the nth RUE and g M m, k represent the interference gain from the mth MUE. Owing to channel reciprocity, g N n, m is equal to g M m, n , g N n, k is equal to g K k, n , and g M m, k is equal to g K k, m . The main notation are listed in Table 1.

Problem formulation
The binary variable x M m, l is the allocation status of M-UE lth channel. Correspondingly, the allocation status variables of R-UE and D2D pair are denoted by x N n, l and x K k, l , respectively We assume that the lth channel is shared by at most one M-UE, one R-UE, and one D2D pair simultaneously, which cannot be shared by two or more UEs in the same type. And each M-UE is assigned a channel x N n, l ł 1, Let p K k, f , p N n, f , and p M m, f represent the optimal transmission power variable of D2D pair k, R-UE n and M-UE m on the channel l, respectively. s 2 represents Table 1. Key system notation.

Notation Description
The interference gain of MUE-m from the RUE-n. The interference gain of D2D pair-k from MUE-m.
The interference gain of RUE-n from D2D pair-k.
The interference gain of MUE-m from D2D pair-k. Gaussian noise. Then, the SINR of D2D pair k, R-UE n, and M-UE m on channel l are calculated According to Shannon's theorem, the throughput of D2D pair k, R-UE n, and M-UE m on channel l is obtained To maximize the sum rate of the system constraint by unique condition of each sub-channel, the maximum power, and the SINR of all users, the optimization is depicted as P1 as follows s:t: x N n, l ł 1, where the P max K , P max N , and P max M are the maximum transmitting power of the D2D pairs, R-UEs, and M-UEs, respectively. The g K min , g N min , and g M min denote the SINR thresholds for D2D pairs, R-UEs, and M-UEs, respectively. As constraint (18) shows, each type of user on channel l should satisfy their QoS requirement. Constraints (15)-(17) force the transmit power of D2D pair (R-UE, M-UE) on channel l to be 0 in case x K k, l = 0 (x N n, l = 0, x M m, l = 0) and their power must not exceed the maximum power. Constraints (13) and (14) ensure that each M-UE is assigned with one subchannel and each R-UE is assigned with at most one sub-channel as well as each D2D pair.

Problem decomposition and algorithm design
The optimization problem P1 is non-convex NP-hard due to the binary nonlinear constraints, which cannot be solved by traditional optimization method. To solve this problem, a two-step scheme containing power allocation and channel allocation is proposed in the following section. Moreover, the detailed communication procedure of the proposed scheme is given.

Power allocation based on Geometric Vertex Search
In this section, the optimal power allocation scheme is addressed by Geometric Vertex Search method. In order to maximize the throughput of the D2D pair, the original problem is transformed to P2, in which the R-UE and M-UE have their own QoS requirements. Here, a simplified case containing one D2D pair, one R-UE, and one M-UE is considered in P2 P2 : max By adopting solid geometric approach, the SINR constraints of the power maximum problem are transformed to equation constraints as follows Considering equations (24)-(26), the power p N n, l , p K k, l , and p M m, l span a 3D space. In the 3D space, equations (24)- (26) are the three planes, where g is the slope of the plane and p is the intercept at the p axis. The transparent cube with the edge of black solid line represents the bound of the individual power. The yellow, black, and gray planes represent the SINR constraint conditions of D2D pairs M-UEs and R-UEs, respectively. The intersection of the three planes in 3D space as shown as black multilayer, which is bounded by yellow, black, and gray planes and the three faces of the transparent cube, forms a feasible power region. In order to avoid a computationally expensive exhaustive search, a near optimal solution is proposed.
It is proved that when the maximized sum rate is achieved, at least one of the powers reaches its maximum power. However, the sum rate expression in equation (19) is non-convex with respect to arbitrary combinations of varying powers. Consequently, for arbitrary number of transmitters, the optimal powers may not necessarily lie on the vertices of the power region, leading to a possibly infinite set of points to test. Since the objective function on the boundary and SINR constraint is a quasi convex function, their optimal solution has the same power. Then, the optimal power of throughput can be obtained by finding the optimal solutions of SINR, which lies on the corners or vertices of the power region as shown in Table 2.
Channel allocation based on 3D matching problem Based on the above results of power allocation, the sum throughput R l m, n, k of M-UE m, R-UE n, and D2D k on channel l is obtained. According to the one-to-one correspondence between M-UEs and channels, we can use R m, n, k to represent R l m, n, k . And a three-dimensional throughput matrix R M 3 N 3 K is obtained by all possible combinations of D2D pairs, M-UEs, and R-UEs. To find a subset T M 3 N 3 K 2 R M 3 N 3 K that maximizes throughput is the objective. Similarly, we define a binary variable j m, n, k to indicate the channel allocation state of M-UE m, R-UE n, and D2D k The original problem P1 transformed into a problem P3. Due to the formulated NP-hard problem is too complex to solve this problem via classical optimization approaches To address it, first P3 is modeled as a maximum weighted three-uniform hypergraph matching problem based on graph theory. Then, the initial solution of the weighted 3D matching problem is solved based on conflict graph by maximum independent set. Finally, the close-to-optimal solution is obtained in quasipolynomial time by iterative updating independent set with m-claw searching.
To facilitate understanding of the channel allocation scheme, some notations and definitions used are introduced as follows. 30 Definition 3 (bipartite conflict graph). A bipartite conflict graph G C = (V C1 , V C2 ) is a special construct in graph theory. Its vertex set is divided into two parts V C1 and If any a 2 V C1 and b 2 V C2 are not disjoint (e.g. adjacent), there is an edge connection between them. Thus, if two sets intersect, they are called to be adjacent (neighbor).
Definition 4 (m-claw). In bipartite conflict graph G C , a mclaw denoted as C m is an induced sub-graph of G C , which consisted of a center and talons T C m . T C m is a set consisting of m independent nodes. There exists at least one common element between the center node and the node containing talons T C m . The center node connected with all independent nodes is the common adjacent node of talons. Based on the above definition, the channel allocation process is as follows: As Figure 2 pair denote the three vertices in hyperedge E H (e.g. e m, n, k denotes the combination of the mth M-UE, nth R-UE, and kth D2D), and the throughput R m, n, k represents the weight of corresponding hyperedge, for example, v(e H m, n, k ) = R m, n, k . Therefore, the channel allocation problem can be modeled as a maximum weight threeuniform hypergraph matching problem. According to graph theory, a matching is a set of non-adjacent edges such that no edge covers common vertex. And maximum matching is defined as a matching with the largest number of edges in all possible matching results. Then, a maximum weight hypergraph matching is to capture a subset of disjoint hyperedges, such that the number of hyperedges and the sum of weight is largest. To simplify the original hypergraph, the hyperedges which do not satisfy equation (27) is removed. Thus, a feasible hypergraph G F is obtained as shown in Figure 2(b), which is a sub-graph of original hypergraph G H . For example, (e H 1, 1, 2 )-entry is one of the feasible hyperedges; however, (e H 1, 1, 1 ) entry is not a feasible hyperedge. If each hyperedge in hypergraph is denoted as a vertex, such as the hyperedges (e H 1, 1, 2 )-entry in Figure 2(c) is presented as a vertex v C 1 , then the feasible hypergraph G F can be transformed as bipartite conflict graph G C . According Definition 3, let V A and V B denote the two vertices set of conflict graph. Define V A as a maximum independent set of G C , then  N (B, S), then the adjacent ver- 2 ) = f(e2,3,2 ), (e4,3,2 ), (e1,2,3 )g. As Figure 2(c) shows, by constructing bipartite conflict graph, the maximum weight hypergraph matching is converted into a maximum independent set V A with largest weight searching problem in conflict graph. Based on Definition 4, we know that for any center node v C 1 , there exist m mutually independent adjacent nodes of v C 1 . Note that, for k-uniform hypergraph, there are at most k + 1 mutually disjoint neighbors.32,33 In other words, the conflict graph corresponding to this hypergraph is k + 1-claw free, for example, any vertex in it cannot find a k + 1-claw. As shown in Figure 2(c), there are three adjacent vertexes for v C 1 , in fact, it only has two talons, as v C 5 and v C 6 occupy the same elements. Similarly, the vertex v C 4 has five adjacent vertexes, but it only possesses three talons, since v C 5 and v C 6 occupy the same elements and v C 7 and v C 8 occupy the same elements. Based on m-claw theory, if there exists a m-claw (including talons v C 1 ,v C 2 , ..., v C m ) to enhance the overall performance for a given center vertex v C a in independent set V A of conflict graph G C , then add T C m into V A and remove all the edges intersecting with them out of V A . After that, the matching result is best.
With the above knowledge, the proposed channel allocation scheme contains the following steps: 1. Construct Hypergraph based on M, N , K, and R m, n, k ; 2. Find the Feasible Hypergraph from the original Hypergraph; 3. Construct a Bipartite Conflict Graph to map the feasible hypergraph; 4. Initialize the solution of the independent set by a heuristic greedy algorithm; 5. Search the m-claw for any vertex of the identified solution; 6. Update the result if a superior performance is achieved by any neighborhood.
The pseudo-code details are shown in the following Algorithm 1 and Algorithm 2. Note that the vertex number of V A is denoted as V A j j.
The procedure of proposed scheme As Figure 3 illustrates, the detailed communication procedure of the proposed scheme is given. The channel state information (CSI) of all H-CRAN uplinks and   D2D links by eNB and RRHs with fronthaul links, which is assumed following quasistatic block fading distribution, is perfectly acquired BBU pool node. Then, the BBU pool further performs baseband signal processing and resource allocation, which consists of three main stages. In Stage 1, the optimal transmit power is found by relaxing the constraints in equation (18), and correspondingly, the throughput matrix is obtained for each MUE, RUE, and D2D pair. In Stage 2, a weighted hypergraph model is constructed for sub-channel allocation problem. In Stage 3, the total achievable throughput system is obtained with the proposed algorithm. Finally, the calculation result is fed back to all the users through eNB and RRHs (Figure 4).

Numerical results
In this section, we present the simulation results and their corresponding analysis. We consider three-tier H-CRAN with one eNB, some RRUs, and some D2D pairs. As shown in Figure 5, the macrocell is assumed with a radius of 1000 m. The RRUs with 200-300 m coverage distance located in the edge of the macrocell. The D2D pairs are randomly distributed in the coverage of the macrocell, and D2D communication range is within 20-60 m. Both M-UEs and R-UEs are correspondingly distributed within the coverage of eNB and RRUs, respectively. Specific simulation parameters and values are listed in Table 3.
To evaluate the performance, we compared the performance of our proposed algorithm with the iterative Kuhn and Munkres (IKM) algorithm in Liu et al., 34 which transforms the 3D matching problem into three iterative two-dimensional matching problems, and the greedy 3D matching algorithm (GDM) in Wang et al., 35 which arrange all the hyperedges in descending order on the basis of their weights and then obtain the sub-optimal solution in a greedy manner with low complexity. In the following simulation results, greedy iterative hypergraph matching (GIHM) is used to represent our proposed algorithm. They all are used to optimize the three-dimensional matching problem. We also evaluate the performance of the three schemes in terms of system throughput and the admitted rate with different parameter settings. All the results averaged over 200 random trials. Figure 6 shows the total throughput at different D2D radius and RRH coverage radius. It is observed that the throughput gain decreases when the distance of the D2D pair or RRH-RUEs become farther apart. The reason is that the D2D link would be weaker with increasing separation distance of the D2D pair or the RRH-RUEs under constant transmit power. Moreover, it is seen that performance GDM is worst. The reason is that the GDM is based on the heuristic descending algorithm and searches an independent set of hyperedges with the highest weight. It ignores the effect of intersecting hyperedges on performance. The IKM outperforms GDM approaches, because it converts the 3D match into a 2D match and then solves it based on iterative algorithm, which can only guarantee the optimality of the 2D match while cannot guarantee the most positive of the three-dimensional match. The proposed scheme performs the best among other

Algorithm 2. A maximum independent set selection based on Greedy Algorithm
schemes, and this is because the scheme is modeled with hypergragh matching and solved by conflict graph based on greedy algorithm and multiple iterations. The near-optimal throughput is achieved as shown in Figure 2. When D2D distance is 20 m, the proposed algorithm can improve throughput performance by 6:25% and 30% as compared to other schemes. Figure 7 is presented to show that the proposed algorithm increases the throughput of the H-CRAN combined with D2D with different number of D2D pairs and R-UEs. Obviously, the total throughput is enhanced as the number of D2D pairs and R-UEs increase. As it was observed, the turning point of the number of D2D pairs is 18. The throughput growth rate is significant when the number of D2D in the system is smaller. However, the throughput curve rises slowly when the number of D2D in the system is larger than it. Since the min (M, N , K) in the three-partite hypergraph G H has limited the maximum number of non-intersecting hyperedges in independent sets, the total throughout grows rapidly when the vertex number in one dimension is less than that of the other two  dimensions; otherwise, the throughout growth rate is small. Moreover, it is easy to see that the proposed algorithms provide a better performance for D2D communication comparing with IKM and GDM algorithms. When D2D pairs' number is 18, the total throughput of our proposed algorithm has a gain about 6% and 34% over that of IKM and GDM algorithms, respectively.
In Figure 8, the impact of the minimum SINR requirement on the users' throughput performance is further evaluated. It is shown that the throughput increases slightly, when the minimum SINR of RUEs is less than 7. The reason is that the increase in the minimum SINR requirement of RUEs will bring larger transmit power of RUEs, which will inevitably increase throughput. When the minimum SINR of RUEs is less than 13 and more than 7, the throughput decreases significantly. This is because in order to meet the minimum SINR requirements of RUEs, the increasing power will bring too much interference to other users, which will reduce the user admitted rate. It is not difficult to see in Figure 11. In addition, the larger minimum QoS requirements of D2D pairs, the smaller total throughput, regardless of the matching algorithm are used. The reason is the similarity of increasing minimum QoS requirement of R-UEs. Once again, the performance of our proposed algorithms is better than the IKM and GDM algorithms. Figure 9 investigates the performance of three schemes in terms of the number of admitted underlaying users with different D2D radius and RRH coverage radius. The result shows that the number of admitted    underlaying users of all the schemes decreases as the D2D pairs and RRH-RUEs distance increase. However, with the increase in D2D distance, the decline rate is gradually slow, and the performance gap between algorithms also decreases. The reason is that a larger distance of D2D pairs or RRH-RUEs represents a weaker link, thus it is more likely that D2D or RRH-RUEs transmission can be accomplished without losing other constraint. When D2D distance is 20 m, the performance of the proposed algorithm increases about 3% and 21% as compared to other schemes.
From Figure 10, the performance of the schemes in terms of the number of admitted underlaying users is investigated for different number of D2D pairs and R-UEs. It is easy to find that the more D2D pairs or RUEs are there, the higher number of admitted underlaying users is, which is consistent with that shown in Figure 7, further verifying that underlaying communication can improve the system throughput. Besides, the number of admitted rate increases distinctly with the enhancement of the number of D2D pairs in small region. Since in small region, the number of D2D pairs' constraint becomes the major factor that narrows the feasible region of the number of admitted underlaying users. While the number of D2D pairs increases gradually, the raise in the admitted rate tends to be gentle owing to the fact that the number of RUEs constrains the admitted rate increment.
From Figure 11, the performance of the schemes in terms of the access number of underlaying users is investigated for different minimum SINR requirements of D2D pairs and R-UEs. It is easy to observe that when g N min is not sufficiently large, the admitted rate decreases slowly with the increase in g N min because the admitted rate of the underlaying users is less affected by low QoS requirement. However, the feasible region of the problem shrinks significantly for large g N min , which causes the fast decline of the admitted rate in system. Besides, it also shows that the admitted rate declines sharply as the minimum SINR of D2D pairs g K min is large, which demonstrates that both g N min and g K min parameters affect admitted rate at the same time. The reason is that it is more difficult for the D2D pairs and RUEs to find suitable channel which satisfy the QoS constrains of MUEs, RUEs, and D2D pairs simultaneously. When g K min is 5, the total throughput of our proposed algorithm surpasses that of IKM and GDM algorithm with amount about 6:3% and 41:5%, respectively.

Conclusion
In this article, a joint optimization problem to maximize the system capacity of hybrid scenario of the H-   CRAN and D2D communication is presented. By stepby-step scheme, the original NP hard problem is decomposed into two sub-problems, for example, the power allocation and the channel allocation. In the first problem, the 3D optimal power programming problem is solved by the sub-channel allocation based on Geometric Vertex Search. Then, with bipartite conflict graph and the m-claw searching, the optimal result of the channel allocation problem based on the 3D hypergraph is achieved. Simulation shows that system throughput and access rate are affected in terms of D2D number, D2D distance, and RUE SINR requirements. The throughput of the proposed scheme has about 6:2% gain over other scheme as well as and the access rate of the proposed scheme has at least 3% gain over other scheme. Therefore, the proposed scheme is beneficial.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.