## Green Broadcast Transmission in Cellular Networks: A Game Theoretic Approach

Cengis Hasan , Jean-Marie Gorce † † and Eitan Altman ‡

† Inria, University of Lyon, INSA-Lyon, 6 Avenue des Arts 69621 Villeurbanne Cedex, France ‡ Inria, 2004 Route des Lucioles, 06902 Sophia-Antipolis Cedex, France { cengis.hasan, eitan.altman } @inria.fr, { jean-marie.gorce } @insa-lyon.fr

Abstract -This paper addresses the mobile assignment problem in a multi-cell broadcast transmission seeking minimal total power consumption by considering both transmission and operational powers. While the large scale nature of the problem entails to find distributed solutions, game theory appears to be a natural tool. We propose a novel distributed algorithm based on group formation games, called the hedonic decision algorithm . This formalism is constructive: a new class of group formation games is introduced where the utility of players within a group is separable and symmetric being a generalized version of parity-affiliation games. The proposed hedonic decision algorithm is also suitable for any set-covering problem. To evaluate the performance of our algorithm, we propose other approaches to which our algorithm is compared. We first develop a centralized recursive algorithm called the hold minimum being able to find the optimal assignments. However, because of the NP-hard complexity of the mobile assignment problem, we propose a centralized polynomial-time heuristic algorithm called the column control producing near-optimal solutions when the operational power costs of base stations are taken into account. Starting from this efficient centralized approach, a distributed column control algorithm is also proposed and compared to the hedonic decision algorithm . We also implement the nearest base station algorithm which is very simple and intuitive and efficiently manage fastmoving users served by macro BSs. Extensive simulation results are provided and highlight the relative performance of these algorithms. The simulated scenarios are done according to Poisson point processes for both mobiles and base stations.

We consider broadcasting under a green-aware objective aiming at reducing the energy consumption which is an important issue in wireless environments[12]. Broadcasting may bring a strong improvement in wireless channels since a common resource (in frequency and/or time) may be used for all destinations. The transmission cost for a base station (BS) to reach all nodes in a multicast group is assumed to be proportional to the power needed to reach the worst mobile among the group, where the worst refers to the mobile receiving the weaker signal which relies on its distance and on additional shadowing effects. We thus consider the situation where there is one common information that every mobile m ∈ M is interested to receive, and which can be obtained from any one of n BSs. The objective is then to achieve a mobiles assignment which minimizes the total power consumption .

Index Terms -broadcast transmission, green networking, combinatorial optimization, game theory

## I. INTRODUCTION

Broadcast scenarios have been widely studied for video or audio broadcasting. More recently, the Multimedia Broadcast/Multicast Service (MBMS) [1] became a requirement of the Long-term Evolution (LTE) specifications to support the delivery of broadcast/multicast data in LTE systems. Broadcast and multicast downlink transmissions make no significant difference at the physical layer. Basically, broadcast services are available to all users without the need of subscribing to a particular service. Therefore, multicasting can thus be seen as 'broadcast via subscription', with the possibility of charging for the subscription [1]. MBMS is intended to be used for some content, such as streaming transmission of a sport or cultural event, but broadcasting may also be of interest to transmit some signalling such as a beacon for time synchronization or for power control purposes.

This work was supported by the ECOSCells project. It was done within the INRIA &amp; Alcatel Lucent Bell-Labs joint research laboratory.

As mentioned above, this problem is relevant not only for streaming data transmission but also for signalling. A corollary question is indeed how many base stations should be kept active to ensure the signalling to a group of mobiles even when they are not transmitting. In low activity periods, such as during the night, it could be relevant to keep active only few BS ensuring the coverage of few active mobiles. Indeed, energy consumption and electromagnetic pollution are main societal and economical challenges that developed countries have to handle. The evolution of cellular networks toward smaller cells offering theoretically higher capacity could in turn lead to an unacceptable increase of the energy expenditure of wireless systems. When decreasing the cell size, the energy consumed for data transmission becomes lower compared to the operational power costs (e.g. power amplifiers, cooler, etc.) of a typical BS. Switching off a BS may then bring significant improvements in energy efficiency. Therefore, we take into account the switching on/off operation in the problem formulation. The overarching problem studied in the sequel is then finding energy-efficient broadcast transmission techniques to reduce spurious energy using distributed schemes.

The mobile assignment problem (MAP) in the context of broadcast transmission that we study in this work is actually a special case of the simple plant location problem (SPLP) [2]. SPLP lies within clustering problems . In the MAP, basically, the objective is to assign the points to at most k clusters so that the sum of all distances between points in the same cluster ( k -clustering) is minimized. In [3], the typical cost for a BS-mobile pair is assumed to be only a function of distance between the BS and the mobile, formulated as

∑ N j ∈C max i ∈ N j d α ij + P j 0 . Here, N j is a cluster of mobiles assigned to BS j , C is the set of clusters, d ij is the distance between mobile i and BS j , α is the path loss exponent and P j 0 is the operational power cost loaded to BS j . This formulation is modified in order to consider the effect of shadowing leading to the following total cost ∑ N j ∈C max i ∈ N j d α ij / Ψ ij + P j 0 where Ψ ij denotes the shadowing effect between mobile i and BS j . Thus, this modification turns the MAP into the SPLP.

While finding the global minimum of the MAP may be identified as an NP-hard problem from SPLP literature, the large scale nature of the cellular network further requires to solve it in a decentralized manner. Thus, game theory appears as a natural tool to cope with both features: distributed decision and NP-hardness. We address this problem by considering the mobiles as players being able to make strategic decisions and the BSs as the strategy identifiers: each mobile has to choose the best BS to be served.

## A. Related Work

Computational geometric approaches to the MAP can be found in [3], [4], [5], [6]. In [4], the authors examined the 1-dimensional version of the MAP, where the effects of shadowing and operational power cost are not taken into account. Polynomial time solutions via dynamic programming are proposed. In [5], authors suggested approximation algorithms (and an algebraic intractability result) for selecting an optimal line on which to place BSs to cover mobiles, and a proof of NP-hardness for any path loss exponent α &gt; 1 .

The papers [7], [8], [9], [10] focused on source-initiated broadcasting of data in static all-wireless networks. Data are distributed from a source node to each node in a network. The main objective is to construct a minimum-energy broadcast tree rooted at the source node. Multi-hop routing is not the scope of our paper.

In [11], the combined problem of (i) deciding what subset of the mobiles would be assigned to each BS, and then (ii) sharing the BSs' cost of multicast among the mobiles is studied. The subset that is wished to assign to a given BS is said to be its target set of mobiles. This problem can be conceived as a coalitional pricing game played by mobiles which is called the association game of mobiles .

## B. Our Contribution

We propose algorithmic solutions for mobile assignment in the context of broadcasting, in order to minimize the overall energy consumption related to transmission and operational powers. In this context, switching off some fraction of BSs is considered to be a way of decreasing dramatically the total energy consumption. Note however that heterogeneous networks include macro and small-cells with or without coordination. It is reasonable to assume that small-cells are subject to switching off operation while macro-cells are always turned on. They can indeed serve moving mobiles in order to decrease the number of hand-offs. Further, since the small-cells are deployed intensively, their transmission power is lower than those of macro-cells, while their circuit power dominates.

Comparing the transmission power about few milliWatts with the operational power costs which may approach tens of Watts, turning off a fraction of BSs is appealing for reducing the total energy footprint of the network. The efforts for turning off some BSs will be concentrated on small-cells serving fixed mobiles[13].

The referred literature mostly concentrates on the geometric aspects of the MAP where basically, the coverage area of a BS is assumed to be a disc which issues from omnidirectional antenna pattern. However, the effect of shadowing , special designed antenna patterns as well as the operational power costs may impact the BS-mobile assignments. In this paper, we take into account these effects by introducing a power cost matrix containing all BS-mobile pairing power costs.

Furthermore, several papers working on coverage optimization deal with static optimization and planning from a centralized point of view. But the dynamic switching-off process associated to the large-scale nature of the network induces to find distributed solutions. To this end, we deal with this problem through a group formation game formulation. Subsequently, we introduce a new algorithm based on group formation games, called hedonic decision algorithm . This formalism is constructive: a new class of group formation games is introduced where the utility of players within a group is separable and symmetric. This is a generalization of parity affiliation games and this hedonic decision algorithm is in fact applicable for any set covering problem.

To prove the efficiency of this approach, we then derive four other methods allowing to solve the initial problem. First, we propose a recursive algorithm called the hold minimum algorithm which solves the considered problem optimally. However, the hold minimum algorithm operates in a centralized way since it requires the whole knowledge for each BSmobile pairing power cost. We then adapt an approach from the SPLP literature to the MAP: a centralized polynomial-time heuristic algorithm is proposed called the the column control which produces optimal assignments when taking into account the operational power cost. This algorithm is also extended to a distributed approach, where each mobile gathers the local information from the BSs located in its range. On the other hand, the nearest base station algorithm , a distributed greedy algorithm which runs in polynomial-time is also evaluated. This algorithm is not efficient if the operational power cost is large, but is very efficient for the fast-moving users served by macro BSs.

The rest of the paper is organized as follows. In section II, the MAP is formulated mathematically as a clustering problem and different formulations are then proposed. In section III, the game framework and the hedonic decision algorithm are proposed. In section IV, we derive other algorithmic solutions for the MAP and their complexity is anlayzed in section V. Finally, we present simulation results in section VI and we expose some conclusions in section VII.

## II. THE GENERIC MAP PROBLEM

We consider the coverage problem in the case of broadcast transmission in cellular networks. We assume that each BS

Fig. 1. Broadcast transmission in cellular networks.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-5c45d6ba-786e-4bc8-b0be-929e1ffb1871.jpg)

transmits simultaneously to the mobiles. The distance between the mobile i and BS j is represented by d ij . The power needed to receive the transmission is given by P r . We consider basic signal propagation model capturing path loss as well as shadowing effect formulated as

$$P _ { i j } = P _ { r } \frac { d _ { i j } ^ { \alpha } } { \Psi _ { i j } }, \quad \quad \ \ ( 1 ) \quad \text{.}$$

where P ij and α denote transmitted power from BS j to mobile i and path loss exponent, respectively. The random variable Ψ is used to model slow fading effects and commonly follows a log-normal distribution.

The required transmission power depends on the mobile having the worst signal level from the BS (Figure 1). At this power level, all mobiles are guaranteed to receive a sufficient power. We also consider the operational power cost denoted as P j 0 which captures the energy expenditure of a typical BS j for operational costs (power amplifiers, cooler, etc.). So, the total power cost (transmission power + operational power cost) of a typical transmission between BS j and mobile i is denoted as

$$p _ { i j } = P _ { i j } + P _ { 0 } ^ { j }. \quad \ \ \ ( 2 ) \, \begin{array} { c } ( \text{see} \\ l o c a \end{array}$$

Let M = (1 , . . . , m ) and N = (1 , . . . , n ) be the sets of mobiles and BSs, respectively. Representing the power cost matrix P = ( p ij ) ∈ glyph[Rfractur] m n × , we assume p ij ∈ [0 , ∞ ) where if P ij &gt; P max , then p ij = ∞ ( P max denotes a maximal power, for instance, in WiFi, it is 100 mW).

## A. The MAP as a Clustering Problem

Clustering is a rich branch of combinatorial problems which have been extensively studied in many fields including database systems, image processing, data mining, molecular biology, etc. [3]. Consider the set of mobiles M ; a cluster is any non-empty subset of M and a clustering is a partition of M . Many different clustering problems can be defined. The mostly studied problems are defined through their objective which is to assign the points to at most k clusters so that either:

- · k-centre : the maximum distance from any point to its cluster centre is minimized,
- · k-median : the sum of distances from each point to its closest cluster centre is minimized,
- · k-clustering : the sum of all distances between points in the same cluster is minimized.

The problem of clustering a set of points into a specific number of clusters so as to minimize the sum of cluster sizes is referred to as min-size k -clustering problem. In [3], the typical cost for a BS-mobile pair is assumed to be only a function of the BS-mobile distance and leads to the formulation: ∑ N j ∈C max i ∈ N j d α ij + P j 0 , where N j is a cluster of mobiles assigned to BS j , C is the set of clusters, d ij is the distance between mobile i and BS j , α is the path loss exponent, P j 0 is the operational power cost loaded to BS j . Since in this paper, we add the effect of shadowing to such a cost, the cost function becomes ∑ N j ∈C max i ∈ N j d α ij / Ψ ij + P j 0 . Note that the shadowing effect breaks the monotonicity with d ij (see Figure 1) and then shifts the problem to a simple plant location problem (SPLP) [2] formulated as follows:

- · Let have n potential facility locations. A facility can be opened in any location j ; opening a facility location has a non-negative cost corresponding to P j 0 in the MAP. Each open facility can provide an unlimited amount of commodity corresponding to unlimited number of mobiles served by a BS in the MAP;
- · there are m customers that require a service. The goal is to determine a subset of the set of potential facility locations, at which to open facilities and an assignment of all clients to these facilities so as to minimize the overall total cost.

However, the SPLP formulation presents a very high complexity, and we rather propose to turn out the problem as a set covering problem, starting from a binary integer formulation.

## B. Binary Integer Formulation

Note that there are at most 2 m -1 possible subsets of M . Each subset can be associated to any BS and the number of combinations is given by η = n (2 m -1) . We note S the collection of total possibilities. The index set of S is denoted by L = (1 , . . . , η ) . Let W = ( w kl ) ∈ glyph[Rfractur] ( m n + ) × η be a 0-1 matrix with w kl = 1 if the node k (i.e. mobile or BS) belongs to the set ( S l ; j ) . Let q = ( q l ) ∈ glyph[Rfractur] η be an η -dimensional vector. The value q l represents the optimal power of ( S l ; j ) ∈ S by which we denote a pair which consists of a set of mobiles S l assigned to BS j . Clearly

$$q _ { l } = \max _ { i \in S _ { l } } p _ { i j }. \quad \quad \quad ( 3 ) \quad \text{succeq} \quad \text{succeq}.$$

The formulation of this problem is given by

$$( P ) \quad & p = \min \sum _ { l \in L } q _ { l } x _ { l } & \quad \text{$\alpha$} \quad \text{$\alpha$} \quad \text{$\alpha$} \quad \text{$\text{of assc}$} \\ & \quad \text{s.t.} \quad \sum _ { l \in L } w _ { k l } x _ { l } = 1, \quad k \in M \cup N, \quad ( 4 ) \quad \text{turns t} \\ & \quad x _ { l } \in \{ 0, 1 \}, \quad l \in L \\ \cdot \quad \cdot \quad \sigma \quad \cdot \quad \cdot \quad \cdot \quad \sigma \quad \sigma \cdot \quad \text{and} \ ( 7$$

where the term ∑ l ∈ L w kl x l = 1 imposes that only one BS is associated to a mobile. By this way, we do not let a mobile to be assigned to several BSs. It follows that the optimal clustering is denoted as C ∗ ⊂ S such that C ∗ = { ( S ∗ ; j ∗ ) ∈ S ∀ , j ∈ N } where ( S ∗ ; j ∗ ) is the optimal pairing.

Thanks to this formulation, we can show that the MAP may be derived as a set partitioning problem. In the MAP, the set M is associated with another set N . Therefore, the collection S contains those subsets of M each of which is associated with every element of the set N . Consider the following example.

Example 2.1: Let us have a power cost matrix given by

$$P = \begin{bmatrix} 3 & 6 \\ 5 & 1 \end{bmatrix}.$$

The collection of total possibilities:

$$\mathcal { S } = \{ ( 1 ; 1 ), ( 2 ; 1 ), ( 1, 2 ; 1 ), ( 1 ; 2 ), ( 2 ; 2 ), ( 1, 2 ; 2 ) \}. \quad ( 6 ) \ \text{Then, } \varrho$$

Recall that ( S l ; j ) denotes the cluster of mobiles S l assigned to BS j . The optimal values for each possibility is given by q = ( q l ) = (3 5 5 6 1 6) , , , , , . Then, we define the following matrix:

$$\mathbf w = \begin{bmatrix} 1 & 0 & 1 & 1 & 0 & 1 \\ 0 & 1 & 1 & 0 & 1 & 1 \\ 1 & 1 & 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & 1 & 1 \end{bmatrix}. \quad \quad ( 7 ) & \text{with} \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The optimal total power is thus calculated by the following binary integer program:

$$\text{may merge program.} & & \text{For ex} \\ p = & \min ( 3 x _ { 1 } + 5 x _ { 2 } + 5 x _ { 3 } + 6 x _ { 4 } + x _ { 5 } + 6 x _ { 6 } ) & & ( 1, 2, 3 \\ s. t. & x _ { 1 } + x _ { 3 } + x _ { 4 } + x _ { 6 } = 1, & & ( 2, 3 ; 1 \\ & x _ { 2 } + x _ { 3 } + x _ { 5 } + x _ { 6 } = 1, & & ( 8 ) & & ( 1 0, 2 7 \\ & x _ { 1 } + x _ { 2 } + x _ { 3 } = 1, & & ( 1, 3 ; 1 \\ & x _ { 4 } + x _ { 5 } + x _ { 6 } = 1, & & \text{Note} \\ & x _ { l } \in \{ 0, 1 \}, \quad l \in ( 1, \dots, 6 ). & & \text{is the ef}$$

The values x 1 = 1 and x 5 = 1 result in the optimal total power of the example scenario, i.e., p = 3 + 1 = 4 with the optimal clustering C ∗ = (1;1) (2; 2) { , } .

However, set partitioning problems are well known to be NP-hard [14]. Consequently, the MAP being a special set partitioning problem is also NP-hard.In the next section, we propose to reduce this complexity.

## C. The MAP as a Set Covering Problem

The previous formulation stated that a unique BS is allowed to serve a mobile. However in terms of pure coverage considerations, the optimal solution may feature some mobiles to be covered by several BS, no matter to which BS the mobile eventually associate with. We then relax the condition of associating only one BS to a cluster of mobiles in ( P ) such that it is now possible to have a cluster of mobiles covered by more than one BS: ∑ l ∈ L w kl x l ≥ 1 . Thus, this arrangement turns the MAP into so called set covering problem.

Consider a set of mobiles S assigned to BS j noted as ( S j ; ) . When a group of mobiles T ⊂ S deviates to another BS k then the cost due to S becomes additive. Let the cost of ( S j ; ) and ( T k ; ) be q S = max i ∈ ( S j ; ) p ij and q T = max i ∈ ( T k ; ) p ik , respectively. We denote the total cost before deviation of T as p and after deviation of T as p ′ , respectively 1 , which can be given by

$$p = p _ { r } + q _ { S },$$

$$p ^ { \prime } = p ^ { \prime } _ { r } + q _ { S \vee T } + q _ { T }, \text{ \quad \ \ } ( 1 0 )$$

where p r and p ′ r are the remaining costs before and after deviation, respectively. There is always a potential (probability) increasing the total cost when a deviation occurs, i.e. p ′ ≥ p . For better observation, let us consider the following power cost matrix:

$$\mathbf P = \begin{bmatrix} 1 0 & 1 5 & 2 5 \\ 2 7 & 2 0 & 3 3 \\ 3 2 & 3 1 & 3 0 \end{bmatrix}.$$

Let ( S j ; ) = (1 2 3; 1) , , and ( T k ; ) = (1 2; 3) , , respectively. Then, q S = max(10 27 32) = 32 , , , q T = max(25 33) = 33 , , and q S T \ = max(32) = 32 resulting in the following total costs:

$$p = 3 2$$

$$p ^ { \prime } = 3 2 + 3 3 = 6 5,$$

where p r = p ′ = 0

r . Utilizing this property, we delete from the collection S all those assignments ( S T j \ ; ) whenever the cost of ( S j ; ) is equal to the cost of ( S T j \ ; ) such that T ⊂ S .

For example, let us consider the last example where M = (1 , 2 3) , and N = (1 2 3) , , . For j = 1 , all possible assignments are (1; 1) , (2; 1) , (3; 1) , (1 2; 1) , , (1 3; 1) , , (2 , 3; 1) , (1 2 3; 1) , , corresponding to the cost vector q = (10 27 32 27 32 32 32) , , , , , , . Note that the cost of (3; 1) , (1 , 3; 1) , (2 , 3; 1) , (1 2 3; 1) , , are equal each other. Therefore,

1 Note that these costs are not optimal. We only would like to show what is the effect of deviation of a group of mobiles from theirs current BS.

we remove (3; 1) , (1 3; 1) , , (2 3; 1) , from the collection. The reduced collection of assignments becomes as following:

$$\mathcal { S } ^ { \prime } = \{ ( 1 ; 1 ), & ( 1, 2 ; 1 ), ( 1, 2, 3 ; 1 ), ( 1 ; 2 ), ( 1, 2 ; 2 ), & \stackrel { \circ } { \langle M, N ^ { n } } \\ & ( 1, 2, 3 ; 2 ), ( 1 ; 3 ), ( 1, 3 ; 3 ), ( 1, 2, 3 ; 3 ) \}. & ( 1 4 ) \quad \text{plays}$$

Note that this property is an extension of the geometric 'coverage range'. If the shadowing was not considered, this approach would reduced to defining the maximal coverage range. But with shadowing the mobile requiring the most power is not always the further one. Thus, the binary integer program of finding the solution of the problem is given by

$$p = & \min ( 1 0 x _ { 1 } + 2 7 x _ { 2 } + 3 2 x _ { 3 } + 1 5 x _ { 4 } + 2 0 x _ { 5 } + 3 1 x _ { 6 } + \int _ { \mathcal { W } } ^ { p \omega } \\ & 2 5 x _ { 7 } + 3 0 x _ { 8 } + 3 3 x _ { 9 } ) & \text{see 1}$$

$$( 1 5 )$$

$$& \varepsilon \omega _ { l } + \varepsilon \omega _ { 8 } + \varepsilon \omega _ { 9 l } & \quad \text{se} \, + \\ \text{s.t. } & x _ { 1 } + x _ { 2 } + x _ { 3 } + x _ { 4 } + x _ { 5 } + x _ { 6 } + x _ { 7 } + x _ { 8 } + x _ { 9 } \geq 1, & \quad \text{cond} \\ & x _ { 2 } + x _ { 3 } + x _ { 5 } + x _ { 6 } + x _ { 9 } \geq 1, & \quad \text{cond} \\ & x _ { 3 } + x _ { 6 } + x _ { 8 } + x _ { 9 } \geq 1, \\ & x _ { l } \in \{ 0, 1 \}, \quad l \in ( 1, \dots, 9 ). & \quad ( 1 5 )$$

The solution of this problem is found to be x 6 = 1 and x l = 0 , ∀ l ∈ (1 , . . . , 9) which fits to the optimal one.

By such an elimination, the size of the collection of assignments reduces from n (2 m -1) to nm . This prove that the set-covering formulation is much more efficient than the set-partitioning one.

## D. Brute-force Search Solution

Enumerating all possible solutions and choosing the one which produces the lowest cost is known as brute-force search or generate and test.

We represent by A = ( a ij ) ∈ glyph[Rfractur] m n × the assignment matrix where a ij ∈ (0 , 1) . If mobile i is assigned to BS j , then a ij = 1 , otherwise a ij = 0 . Notice that each row of the assignment matrix includes only unique '1' which means that a mobile is served by only one BS, i.e. ∑ j a ij = 1 . This is not in contradiction with our former remark about the possibility of having a mobile covered by several BS. Here now, we decide to associate a mobile to a BS. So if a mobile is covered by several BSs, the serving BS can be anyone of this covering set. Denoting the collection of the assignment matrices A , actually, we formalize the problem as following:

$$p = \min _ { \mathbf A \in \mathcal { A } } \left ( \sum _ { i \in M } \max _ { j \in N } \mathbf A \otimes \mathbf P \right ), \quad \quad ( 1 6 ) \, \stackrel { \text{alizat} } { T h } \,,$$

where ⊗ is the element-wise product. Note that the total number of possibilities of assignment matrices can be calculated as |A| = n m .

## III. DECENTRALIZED SOLUTION: THE HEDONIC DECISION (HD) ALGORITHM

We now turn to the study of decentralized methods for solving the MAP. Our approach is based on group formation game (see [19], [18], [20]).

A.

The Game Model

A group formation game is represented by a triple G = 〈 M,N m , ( u x x ) ∈ M 〉 where M = { 1 2 , , . . . , m } is the set of players (i.e. the mobiles), N is the set of strategies (i.e. the BSs) shared by all the players and u x : N m →glyph[Rfractur] is the utility function of player x ∈ M . Each player x ∈ M chooses exactly one element from the n alternatives in N . The choices of players are represented by σ = { s , s 1 2 . . . , s m } ⊆ N m which is called the strategy-tuple ( s x shows the strategy chosen by player x ). A partition of players according to strategy-tuple σ is denoted as Π( σ ) = { ( G σ j ) j ∈ N } where G σ j is the group of players choosing the strategy j .

We assume the two following conditions and we will see later how we can define the utilities to achieve these conditions:

- 1) Separability : The utility of any player in any group is said to be separable if its utility can always be split as a sum according to:

$$u _ { x } ( s _ { x }, \sigma _ { - x } ) = \sum _ { y \in G _ { s _ { x } } ^ { \sigma } } v _ { x } ( y ; s _ { x } ), \quad \ \ ( 1 7 )$$

where v x : { M N ; } → glyph[Rfractur] may be interpreted as the gain of player x from player y if x chooses strategy s x . Note that v x ( x s ; x ) is the utility of player x when it is the only player choosing strategy s x . Thus, the separability property states that utility transfers among a group of users sharing the same strategy is done such that the utility granted to one user is a sum of utilities granted individually by each partner in the group.

- 2) Symmetricity : The utility is said to be symmetric if the individual gain of x from y is equal to the gain of y from x when they both share the same strategy s :

$$v _ { x } ( y ; s ) = v _ { y } ( x ; s ), \ \forall x, y \in M. \quad ( 1 8 )$$

Therefore, this symmetric utility can be referred to as: v x ( y s ; ) = v y ( x s ; ) = v x, y ( ; s ) , ∀ x, y ∈ M v x,y s . ( ; ) is called the symmetric bipartite utility of player x and y while the common strategy is s .

Actually, the game defined above is a straightforward generalization of party affiliation games [23].

Theorem 3.1: G is a potential game.

Proof: A non-cooperative game is a potential game [21] whenever there exists a function Φ

$$\Phi ( s _ { x }, \sigma _ { - x } ) - \Phi ( s ^ { \prime } _ { x }, \sigma _ { - x } ) = u _ { x } ( s _ { x }, \sigma _ { - x } ) - u _ { x } ( s ^ { \prime } _ { x }, \sigma _ { - x } ), \ ( 1 9 )$$

meaning that when player x switches from strategy s x to s ′ x the difference of its utility can be given by the difference of a function Φ . This function is called a potential function .

The strategy-tuple that maximizes the potential function is a Nash equilibrium in the game.

Let us choose as following the potential function Φ :

glyph[negationslash]

glyph[negationslash]

$$\text{Le us conose as following one potential function } \Psi \colon \quad & \text{$b$. $sepa$} \\ \Phi ( \sigma ) & = \sum _ { j \in N } \left [ \sum _ { a \in G _ { j } ^ { \sigma } } v ( a, a ; j ) + \frac { 1 } { 2 } \sum _ { a \in G _ { j } ^ { \sigma } } \sum _ { b \in G _ { j } ^ { \sigma } ; b \neq a } v ( a, b ; j ) \right ] \quad & \text{Recal} \\ & = \sum _ { j \in N } \sum _ { a \in G _ { j } ^ { \sigma } } v ( a, a ; j ) + \frac { 1 } { 2 } \sum _ { j \in N } \sum _ { a \in G _ { j } ^ { \sigma } } \sum _ { b \in G _ { j } ^ { \sigma } ; b \neq a } v ( a, b ; j ). \quad & \text{groug C} \\ & = \text{motor} \\ \text{Aditalv here we take the sum of cinale lower initialec and}$$

Actually, here we take the sum of single player utilities and the half of total symmetric bipartite utilities. Let us rewrite the potential function as following:

glyph[negationslash]

$$\text{ne potential function as following:} \\ \Phi ( s _ { x }, \sigma _ { - x } ) = v ( x, x ; s _ { x } ) + \sum _ { y \in G _ { z } ^ { g } _ { x } } v ( x, y ; s _ { x } ) \\ + \underbrace { \sum _ { j \in N } \sum _ { a \in G _ { j } ^ { g } \vee x } v ( a, a ; j ) + \frac { 1 } { 2 } \sum _ { j \in N } \sum _ { a \in G _ { j } ^ { g } \vee x } \sum _ { b \in G _ { j } ^ { g } \vee x ; b \neq a } v ( a, b ; j ) } _ { I } \,. \quad \text{while} \quad \text{this, as:} \\ \text{When never } \tau \text{ such these from } \tau \text{ to } \tau ^ { \prime } \text{ the other lower} \quad.$$

When player x switches from σ x to σ ′ x (the other players do not change their strategies), then the strategy-tuple is transformed from σ to σ ′ , and the potential becomes glyph[negationslash]

$$& \text{transformed from $\sigma$ to $\sigma$}, and the potential becomes \\ & \Phi ( s ^ { \prime } _ { x }, \sigma _ { - x } ) = v ( x, x ; s ^ { \prime } _ { x } ) + \sum _ { y \in G ^ { \sigma } _ { s ^ { \prime } _ { x } } } v ( x, y ; s ^ { \prime } _ { x } ) } \quad \text{where $\theta$ is} \\ & + \underbrace { \sum _ { j \in N } \sum _ { a \in G ^ { \sigma } _ { j } ^ { \prime } \vee x } v ( a, a ; j ) + \frac { 1 } { 2 } \sum _ { j \in N } \sum _ { a \in G ^ { \sigma } _ { j } ^ { \prime } \vee x } \sum _ { b \in G ^ { \sigma } _ { j } ^ { \prime } \vee x \colon b \neq a } v ( a, b ; j ) } _ { I ^ { \prime } }. \text{system. $T$} \\ & \text{$given by $t$} \\ & \text{$u_{x}(j,$} \\ & \text{Note that $I=I^{\prime}$ since the total utility due to the other nlayers}$$

Note that I = I ′ since the total utility due to the other players is equal both in σ and σ ′ . Thus, the difference of potentials is given by

$$& \Phi ( s _ { x }, \sigma _ { - x } ) - \Phi ( s ^ { \prime } _ { x }, \sigma _ { - x } ) = v ( x, x ; s _ { x } ) - v ( x, x ; s ^ { \prime } _ { x } ) \\ & + \sum _ { y \in G ^ { \sigma } _ { s _ { x } } } v ( x, y ; s _ { x } ) - \sum _ { y \in G ^ { \sigma } _ { s ^ { \prime } _ { x } } } v ( x, y ; s ^ { \prime } _ { x } ). ( 2 3 ) \quad \text{symmet} \\ & \quad \text{express}$$

On the other hand, the difference of the utility of player x is calculated as

$$& u _ { x } ( s _ { x }, \sigma _ { - x } ) - u _ { x } ( s _ { x } ^ { \prime }, \sigma _ { - x } ) = v ( x, x ; s _ { x } ) - v ( x, x ; s _ { x } ^ { \prime } ) \\ & + \sum _ { y \in G _ { s _ { x } } ^ { \sigma } } v ( x, y ; s _ { x } ) - \sum _ { y \in G _ { s _ { x } ^ { \prime } } ^ { \sigma } } v ( x, y ; s _ { x } ^ { \prime } ). \quad & ( 2 4 ) \quad \text{Now, le} \\ & \dots \cdots.$$

By this result we conclude that

$$\Phi ( s _ { x }, \sigma _ { - x } ) - \Phi ( s _ { x } ^ { \prime }, \sigma _ { - x } ) = u _ { x } ( s _ { x }, \sigma _ { - x } ) - u _ { x } ( s _ { x } ^ { \prime }, \sigma _ { - x } ) \ ( 2 5 )$$

which proves that G is a potential game. Thus, G admits a Nash equilibrium in pure strategies σ ∗ which results in the partition Π( σ ∗ ) and which maximizes the potential function Φ .

Corollary 3.1: The proof 3.1 is constructive: any group formation game possessing separable and symmetric utility gain of players within a group always converges to a pure Nash equilibrium.

.

.

## B. Separable and Symmetric Gain Allocation

Recall that the required power for serving the group of mobiles G σ j by BS j is denoted as max x ∈ G σ j p xj . We represent by U G ( σ j ) = -max x ∈ G σ j p xj ≤ 0 , as a utility arising due to group G σ j . Note that U G ( σ j ) ≤ 0 , ∀ G σ j ⊆ M, ∀ j ∈ N is a monotonically decreasing function [16].

The clustering profit due to mobile x and y in BS j is given by

$$\Delta ( x, y ; j ) = U ( x, y ; j ) - [ U ( x ; j ) + U ( y ; j ) ] \quad ( 2 6 )$$

where U x j ( ; ) = -p xj is the utility of player x when is served alone in BS j . Therefore, ∆( x, y ; j ) = p xj + p yj -max( p xj , p yj ) = min( p xj , p yj ) .

Remark 3.1: Note that the clustering profit is a useful metric for evaluating a group of mobiles. Whenever a group of mobiles are near each other, then the clustering profit is high; thus, assigning this group to only one BS is almost always efficient.

To ensure separability and symmetricity, we propose to choose the symmetric bipartite utility according to:

glyph[negationslash]

$$v ( x, y ; j ) = \begin{cases} \theta \Delta ( x, y ; j ), & \text{if $x\neq y$} \\ - p _ { x j }, & \text{if $x=y$}, \end{cases} \quad ( 2 7 )$$

where θ is called as clustering weight which is a parameter that must be adjusted according to the environment. We will show later how it can impact the convergence point of the system. Thus, the utility function of any player x ∈ M is given by u x ( j, σ -x ) = ∑ y ∈ G σ j v x, y ( ; j ) then, glyph[negationslash]

$$\overline { \vdots 2 \rangle } \quad u _ { x } ( j, \sigma _ { - x } ) = \quad & \quad \ \ ( 2 8 ) \\ \underset { \text{is} } { \overline { \text{rs} } } \quad \begin{cases} \theta \sum _ { y \in G _ { j } ^ { \sigma } } \min ( p _ { x j }, p _ { y j } ) - p _ { x j }, & \text{if $G_{j}^{\sigma} \neq x$} \\ - p _ { x j }, & \text{if $G_{j}^{\sigma} = x$}. \end{cases}$$

## C. Interpretation of Clustering Weight

Let us rewrite the potential function according to defined symmetric bipartite utility. Then, using eq. (20), we may express the corresponding potential function which is equal to

$$\text{er } x \text{ is } \quad \Phi ( \sigma ) = \sum _ { j \in N } \left [ - \sum _ { a \in G _ { j } ^ { \sigma } } p _ { a j } + \frac { \theta } { 2 } \sum _ { a \in G _ { j } ^ { \sigma } } \sum _ { b \in G _ { j } ^ { \sigma } \wedge a } \min ( p _ { a j }, p _ { b j } ) \right ]. \\ \text{$(24)} \quad \text{Now let us consider that we can order the mobiles associated}$$

Now, let us consider that we can order the mobiles associated to BS j , according to the required power. The ordered set of mobiles related to BS j is represented by ˜ G σ j . Observe that we can now compute the potential function as following:

$$\text{its a} _ { \substack { \mathbf i t i o n \\ \text{in the} \\ \text{action} } } \quad \Phi ( \sigma ) = \sum _ { j \in N } \left [ - \sum _ { a \in \tilde { G } _ { j } ^ { \sigma } } p _ { a j } + \theta \sum _ { i = 1 } ^ { | \tilde { G } _ { j } ^ { \sigma } | - 1 } ( | \tilde { G } _ { j } ^ { \sigma } | - i ) p _ { i j } \right ]. \ \ ( 3 0 )$$

According to this result and the description of individual costs, we can state the following properties:

- 1) If θ is very small, i.e. the dominant term in the potential function and in the symmetric bipartite utility is the individual power p aj , then with a very low θ , each

mobile will privilege an association to the nearest BS. This is exactly the case when θ = 0 .

- 2) When θ increases a mobile may decide to leave its nearest neighbour if the lost in power is compensated by the second term. Suppose that the mobile wants to associate to a BS k , where all other mobiles currently associated with experience a better channel. Then the gain to associate to this second BS for this user will be θ ∑ p ij , i.e. the sum of all powers of mobiles already associated with this BS, weighted by θ . It is clear that the mobile will be joining either a cell having already strong power terms or a huge number of users.
- 3) If θ becomes very large, we can expect that all mobiles will converge to the same BS which means that only one BS is active.

## D. Network Topology as a Result of Best-reply Dynamics

Consider the setting in which only one player decides its strategy. It is called as best-reply dynamics when a player chooses the strategy which maximizes its utility. When there is no any player which can improve its utility, then this network topology, i.e. Π( σ glyph[star] ) , corresponds to a Nash equilibrium. Note that any local maximum in the potential function is a Nash equilibrium. Therefore, the network topology obtained by bestreply dynamics accounts for a local maximum of the potential function. Total power cost related to Π( σ glyph[star] ) can be given by

$$p = \sum _ { G _ { j } \in \Pi ( \sigma ^ { * } ) } \max _ { i \in G _ { j } } p _ { i j } \quad \quad ( 3 1 ) \quad \quad A. \ \, \mathcal { C }$$

where G j is the group of mobiles associated with BS j in the case of stable strategy-tuple σ glyph[star] .

Assuming that each mobile is capable to discover those BSs that can transmit to it, we can produce a scheduler in the following way: each BS generates a random clock-time for all those mobiles that it can transmit; then each mobile selects randomly a clock-time from those BSs that it can discover. We need to produce the clock-times by such a way that the collision of the turns of mobiles is minimal. In case of a collision, the clock-times of the corresponding mobiles are regenerated by corresponding BSs.

In Algorithm 1, the pseudo-code of the HD is given. Note that this is an algorithm performed in both BS and mobile sides by an exchange of the information in a separated channel.

Corollary 3.2: In the literature, the use of game models for set covering problems is called as set covering games [24], [25], [26]. The HD algorithm is a novel approach for set covering games. This algorithm is suitable for any set covering problem and facility activation problems where the agents are allowed to make strategic decisions

## IV. EFFICIENT ALGORITHMS FOR THE MAP

In this section, we propose different algorithmic solutions to evaluate and compare the efficiency of the distributed algorithm. Centralized algorithms exploit the set covering problem formulation since the search space is the smallest. However, binary integer linear programs are known to be NP-complete and thus we introduce two algorithms based on

## Algorithm 1 The Hedonic Decision

Base Station :

Check stability

while there is no stability do

Send information to each mobile about the current parti- tion

Check stability

end while

Mobile :

while there is no stability do

Determine the preferred BS according to eq. (28) Send information to the preferred BS

end while

dynamic programming: the hold minimum algorithm and the column control algorithm . Then, we develop the distributed version of the column control algorithm. A greedy solution of the problem is introduced as the nearest base station approach. The nearest BS and the column control are known and already used in the literature for SPLP problems. We adapt these algorithms to the MAP.

Because of the large scale nature of the collection set S , we rather develop the algorithms by making all operations on the power cost matrix. This approach foster the iterative removals of elements in the collection set, and ensure a faster convergence.

## A. Optimal Solution: The Hold Minimum (HM) Algorithm

The HM algorithm solves the problem optimally . We explain the algorithm by an example. Consider the power cost matrix which is given by

$$\begin{array} { c c c } y \\ & & P = \begin{bmatrix} 9 & 3 \\ 1 & 4 \\ 2 & 8 \end{bmatrix}. \quad & & ( 3 2 ) \end{array}$$

The power cost matrix can also be shown as P = ( p , p 1 2 , . . . , p n ) , where p j = ( p 1 j , p 2 j , . . . , p mj ) T . In each step, the algorithm removes a group of values p ij of the power cost matrix. Removing p ij means that we eliminate those clusters that include the mobile i and BS j from the collection set. The algorithm compares maximum n clusterings and holds only the clustering minimizing the total cost. Thus, it terminates in a step Q where each mobile is assigned to only one BS. In step s , the power cost matrix and collection set is denoted as P [ s ] = ( p 1 [ s , p ] 2 [ s , . . . , p ] n [ s ]) and S [ s ] , respectively.

Let us now turn to the example. In the initial step s = 0 , we assume that P [0] = P and S [0] = S given by

$$\mathcal { S } [ 0 ] = \{ ( 1 ; 1 ), ( 2 ; 1 ), ( 3 ; 1 ), ( 1, 2 ; 1 ), ( 1, 3 ; 1 ), ( 2, 3 ; 1 ), \\ ( 1, 2, 3 ; 1 ), ( 1 ; 2 ), ( 2 ; 2 ), ( 3 ; 2 ), ( 1, 2 ; 2 ), ( 1, 3 ; 2 ), \\ \text{ited}$$

Recall that assigning a cluster of mobiles S l to BS j has a cost max i ∈ S l p ij . Therefore, if we find the maximum value of p j , we obtain the total cost in case of all mobiles in column j are assigned to BS j . For example, max p 1 = max(9 1 2) = 9 , , .

This means that if all mobiles are assigned to BS 1, then the total cost is 9 .

The algorithm runs as following: in step s = 1 , we find the maximum value of each column of power cost matrix, then eliminate all values in power cost matrix except minimum of the calculated maximum values. Namely, max(9 1 2) , , = 9 and max(3 4 8) = 8 , , , then 9 is eliminated by putting an ∞

$$\mathbf P [ 1 ] = \begin{bmatrix} \infty & 3 \\ 1 & 4 \\ 2 & 8 \end{bmatrix}. \quad \begin{array} { c } \text{In e} \\ ( 3 4 ) & \text{pow} \\ \text{total} \\ \dots \end{array}$$

Thus, the collection set reduces to the following

$$\mathcal { S } [ 1 ] = \{ ( 2 ; 1 ), ( 3 ; 1 ), ( 2, 3 ; 1 ), ( 1 ; 2 ), ( 2 ; 2 ), ( 3 ; 2 ), ( 1, 2 ; 2 ), \quad & \text{exactly} \\ ( 1, 3 ; 2 ), ( 2, 3 ; 2 ), ( 1, 2, 3 ; 2 ) \}. \quad & ( 3 5 ) \quad & \frac { } { \text{$\text{$\text{$\text{$\text{$al$ori$it}$}$}$} } }$$

First column contains an ∞ which means that mobile 1 must be assigned to another BS (i.e. in this example, obviously BS 2 ). In fact, this represents the recursiveness of the algorithm where we run the algorithm for a sub power cost matrix. In this simple example, the sub power cost matrix is 3 . In general case, the algorithm does the following

$$P _ { j } [ s ] = \begin{cases} \ \max p _ { j } [ s ] + \mathfrak { h } ( \mathbf P _ { j } ^ { s u b } [ s ] ), & \text{if sub power cost matrix;} \\ \ \max p _ { j } [ s ], & \text{otherwise.} \end{cases} \quad \ \frac { \cdot } { \cdot }$$

(36)

where P j [ s ] represents the total cost if we assign all mobiles to BS j except those that can not be assigned to, and the optimal cost occurring due to sub power cost matrix P sub j [ s ] in step s . Here, h : glyph[Rfractur] m n × → {glyph[Rfractur] , A} is the function which gives the optimal value and assignments obtained by running HM algorithm.

For s = 2 , we calculate P 1 [2] = max(1 2)+ , h (3) = 2+3 = 5 , where P sub 1 [2] = (3) . On the other hand, we do not need to calculate P 2 [2] since it is kept in the memory. Therefore, P 2 [2] = P 2 [1] . Then, the algorithm holds minimum value of min( P 1 [2] , P 2 [2]) = P 1 [2] = 5 meaning that we remove 8 resulting in the following

$$\mathbf P [ 2 ] = \begin{bmatrix} \infty & 3 \\ 1 & 4 \\ 2 & \infty \end{bmatrix}, \quad \quad \ ( 3 7 ) \quad B. \ C \\ L _ { \Xi }$$

and

$$\mathcal { S } [ 2 ] = \{ ( 2 ; 1 ), ( 3 ; 1 ), ( 2, 3 ; 1 ), ( 1 ; 2 ), ( 2 ; 2 ), ( 1, 2 ; 2 ) \}. \ ( 3 8 ) \ \ P _ { 0 } \gg i$$

Then, for s = 3 , P 2 [3] = max(3 4) + , h (2) = 4 + 2 = 6 , where P sub 2 [3] = (2) , and P 1 [3] = P 1 [2] . We remove 4 , since min( P 1 [3] , P 2 [3]) = P 1 [3] . This gives the following matrix, collection set, assignments, and optimal total power,

$$\mathsf P [ 3 ] = \begin{bmatrix} \infty & 3 \\ 1 & \infty \\ 2 & \infty \end{bmatrix}, \quad \quad \quad ( 3 9 )$$

$$\mathcal { S } [ 3 ] = \{ ( 2 ; 1 ), ( 3 ; 1 ), ( 2, 3 ; 1 ), ( 1 ; 2 ) \} = \{ ( 2, 3 ; 1 ), ( 1 ; 2 ) \},$$

$$\mathfrak { h } ( \mathbf P ) \coloneqq \left \{ p = 5, \mathbf A = \begin{bmatrix} 0 & 1 \\ 1 & 0 \\ 1 & 0 \end{bmatrix} \right \}, \quad \begin{pmatrix} \mathfrak { m o u } \\ \mathfrak { w } ( 4 1 ) & | N _ { 3 } | \mathfrak { w } \\ \mathfrak { m o b i l } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathdag \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfra \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfro \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrum \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrig { \mathfrak { w } } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \math frak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w } \mathfrak { w                                                                                                                                                                                                        
                                                                                                                                                                                                        }                                                                                                                                                                                                       	                                                                                                                                                                                                        _                                                                                                                                                                                                       .                                                                                                                                                                                                      .                                                                                                                                                                                                        "                                                                                                                                                                                                        :                                                                                                                                                                                                        -                                                                                                                                                                                                        <                                                                                                                                                                                                        +$$

respectively. The pseudo-code of this algorithm is given in Algorithm 2.

Theorem 4.1: The HM algorithm terminates in finite step. At this step, the total power cost is minimum.

Proof: The power cost matrix is transformed in each step by

$$\mathbf P [ 0 ] \to \mathbf P [ 1 ] \to \dots \to \mathbf P [ Q ]. \text{ \quad \ \ } ( 4 2 )$$

In each step, at least n -1 values are removed from the power cost matrix. Removed values are those that increase the total cost. The value that is held is the minimum one in the corresponding step. So, in the terminal step s = Q , it arrives to such a cost that is the lowest. Each mobile is assigned to exactly one BS in the terminal step.

## Algorithm 2 The Hold Minimum

$$\text{st} \quad \text{function $(p,A) = \mathfrak { h } ( \mathbf P )$$

while each row includes several ' 1 's in A do

k

←

1

for j = indices of columns of P not including all ∞ do V k ← maximum of column j except ∞

if any row in column

j

includes

∞

then sub

P

←

rows including

$$\begin{smallmatrix} \mu & \cdot & \cdot & \cdot &$$

end if

k

←

k

end for

(

i

min

, j min

Hold

V

+1

, V

min causing

min using

)

←

(

min

i

min

∞

of

V

, j min

P

ij

≥

V

min in

P

Put

0

in all indices causing

## end while

p

←

∑

max

A

⊗

P

end

The complexity of this algorithm is assessed in Section ?? .

B. Greedy Solution: The Column Control (CC) Algorithm

Let us denote the CC algorithm by c : glyph[Rfractur] m n × → {glyph[Rfractur] , A} . This algorithm exploits the assumption that the operational power cost P 0 is fairly higher than the transmitted power, i.e. P 0 glyph[greatermuch] P ij . The aim here is to assign many mobiles to only one BS. Recall that cluster N j (set of mobiles that can be assigned to BS j ) has the transmission cost max i ∈ N j P ij and the operational cost P j 0 . For better understanding, consider the following power cost matrix in which the operational power cost is assumed to be 12 W :

$$\begin{array} {$$

In such a scenario, we can assign | N 1 | = 3 mobiles to BS 1 with cost 12 50 . , | N 2 | = 4 mobiles to BS 2 with cost 12 45 . , | N 3 | = 4 mobiles to BS 3 with cost 12 43 . , and | N 4 | = 3 mobiles to BS 3 with cost 12 35 . .

P

)

and put ij

P

∞

in

≥

V

min in

A

all indices

The logic behind the CC algorithm is the following:

- 1) Find how many mobiles can be assigned to each BS
- 2) Choose the BS to which it can be assigned the most mobiles
- 3) If there are multiple BSs in the state of 2), then choose the BS which can serve the mobiles with minimal cost

Applying these rules to the last example, it turns out that BS 3 can cover the most mobiles | N 3 | = 4 which are N 3 = (1 , 2 3 4) , , with the lowest cost 12 43 . . Then, the algorithm assigns only a cluster of mobiles. In the following step, the CC algorithm performs the same rules to the remained mobiles which produces the sub power cost matrix P sub . In the last example, it is given by P sub = [ ∞ ∞ ∞ 12 29 . ] . This results in the assignment of mobile 5 to BS 4, because mobile 5 can not be assigned to other BSs. This shows the recursiveness of the CC algorithm. A pseudo-code is given in Algorithm 3.

Formally, in step s , we denote as following the set of mobiles assigned to BS j : R s [ ] = { Assigned mobiles in step s } . Thus, the collection set is reduced as following in step s :

$$\mathcal { S } [ s ] & = \{ \mathcal { S } [ s -$$

Considering the last example, in step s = 1 , R [1] = (1 , 2 3 4) , , . By assigning these mobiles to BS 3, the CC algorithm removes those assignments which include mobiles (1 2 3 4) , , , except ( S k ; ) = { (1 , 2 3 4) 3 , , , } . Thus, the total power and assignments are given by

$$\mathfrak { c } ( \mathbf P ) \coloneqq \left \{ p$$

Moreover, for any power cost matrix, we conclude that in a final step Q , the CC algorithm converges to the case where each mobile is assigned to one BS.

## Algorithm 3 The Column Control

function ( p, A ) = c ( P )

v

′

=

Find how many

v

min

∞

←

min

v

′

v

′′

=

Find which row of

v

′

is equal to

v

min for l = columns of P determined by | v ′′ | do

V

l

←

maximum value of column

l

of

P

end for

Find V min = min V l and the corresponding column l min if | v ′ | == 0 then

Put 1s to the column l min in A

else

A ′ ← Put 1s to the column l min in A

Find sub power cost matrix P sub which is composed by those mobiles that are not assigned

A

sub

←

Find assignments by running

c

(

P

sub

)

## end if

$$\mathbf A \leftarrow \text{Combine } \mathbf A ^ { s u b$$

$$p \leftarrow \sum \max \, \mathbf A \otimes \mathbf P$$

end has each column of

P

## C. Distributed Column Control (DCC) Algorithm

Assume that each BS broadcasts its own power vector and identities of mobiles that it can serve. Recall that we denote the power vector of BS j as p j . The power vector of BS 1 given in the power cost matrix of eq. (43) is p 1 = (12 50 12 30 12 20 . , . , . , ∞ ∞ , ) T . From a practical point of view, the BS broadcasts only the identity and associated costs of the mobiles which are associated with it. The infinity values are not broadcast. Therefore, the BS broadcasts the power vector as p 1 = (12 50 12 30 12 20) . , . , . T , and an identity vector represented as h j = ( h jk ) ∈ N | p j | . For example, h 1 = (1 2 3) , , T .

Moreover, each mobile receives power vectors from all BSs that can transmit to it. Then, each mobile generates the power cost matrix from received power vectors. For example, mobile 1 receives from BS 1, BS 2, and BS 3 the power vectors

$$p _ { 1 } = ( 1 2. 5 0, 1 2. 3 0, 1 2. 2 0 ) ^ { T },$$

$$p _ { 2 } = ( 1 2. 4 0, 1 2. 3 0, 1 2. 4 5, 1 2. 4 3 ) ^ { T },$$

and

$$p _ { 3 } = ( 1 2. 3 2, 1 2. 4 3, 1 2. 1 5, 1 2. 2 5 ) ^ { T }$$

with identity vectors h 1 = (1 2 3) , , T , h 2 = (1 2 3 4) , , , T , and h 3 = (1 2 3 4) , , , T , respectively. Mobile 1 decides that the power cost matrix is as following:

$$\begin{bmatrix} & & \hookrightarrow & & \\ 1 2. 5 0 & 1 2. 4 0 & 1 2. 3 2 \\ 1 2. 3 0 & 1 2. 3 0 & 1 2. 4 3 \\ 1 2. 2 0 & 1 2. 4 5 & 1 2. 1 5 \\ \infty & 1 2. 4 3 & 1 2. 2 5 \end{bmatrix}.$$

Note that mobile 1 realizes from h 2 and h 3 that BS 1 can not transmit to mobile 4. Therefore, it puts an infinite cost corresponding to mobile 4 and BS 1 in the power cost matrix.

By this rule each mobile determines its own power cost matrix. Thus, each mobile finds the assignments according to CC algorithm, and selects the BS from which it will receive data. For example, mobile 1 obtains the following assignment matrix by running CC algorithm

$$\begin{bmatrix} 0 & 0 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 1 \end{bmatrix}.$$

It turns out that mobile 1 chooses BS 3 for reception the broadcast data.

Remark 4.1: Through the advantage of the decentralization, the DCC algorithm makes possible the following: if the mobiles do not send any assignment information to a BS, then the corresponding BS is considered to be switched off . On the other hand, the CC algorithm is centralized, therefore, the network determines according to the assignments which BS is switched off .

## D. Greedy Solution: The Nearest BS (NBS) Algorithm

To solve the problem heuristically, the easiest way is to assign each mobile to the nearest BS. By 'nearness', we do not mean a geographical measure, instead, it is the lowest

power cost that the corresponding mobile needs from the corresponding BS. So, the mobile selects the BS transmitting with the lowest power. We assume that a mobile is capable to know the power costs corresponding to those BSs that can transmit to it.

Clearly, mobile i knows the vector p i = ( p i 1 , p i 2 , . . . , p in i ) where n i denotes the number of BSs that mobile i can be served. Then, mobile i only calculates the minimal value of p i , and chooses the corresponding BS,

$$a _ { i, j } = 1 \colon \ j = \arg \min _ { j } p _ { i },$$

where a i,j ∈ A is the BS that mobile i selects. Actually, this corresponds to remove from collection set S all assignments related to mobile i and the BSs being out BS j .

The NBS algorithm is very efficient and quick, and in most cases, it gives optimal assignments when the operational power cost P 0 is neglected (See Section VI).

## E. Greedy Set-Cover Algorithm (greedy-SC)

In [17], a greedy heuristic for the set covering problem is proposed by Chvatal. We utilize this algorithm in simulation results (section VI) to compare it with the proposed algorithms.

## V. TIME COMPLEXITY ANALYSIS

In this section, we estimate and compare the time complexity of the proposed algorithms.

Theorem 5.1: The time complexity of the HM algorithm is O m n ( 3 ) + ( n -1) ∑ m s =1 T m ( -s, n -1) + mO n ( ) .

Proof:

Note that at most n -1 operations are required to find the maximum value of a column of the power cost matrix in one step, and a single operation for finding the minimum of a vector having n values. In each step, there is at least one element in a column that is removed. At most m steps are needed until a convergence to the optimal case. Selection algorithms (finding maximum or minimum) have a time complexity in O n ( ) [27]. We can express as following the time complexity including the recursive property of the algorithm:

We also simulate the last recurrence resulting in the following graphical characteristics:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-211edcc6-beb2-4269-b624-c2e93b50583c.jpg)

m

In the axis denoting T m ( ) , a log-scale is used. The figure illustrates that the complexity grows exponentially.

Theorem 5.2: The time complexity of the CC algorithm is O mn ( ) .

Proof: This result is known and a complete proof is available [28].

Theorem 5.3: The time complexity of the NBS algorithm is O n ( ) .

Proof: It is straightforward because the only operation that is performed in the NBS algorithm is to find the minimum value of a vector of power cost matrix having a dimension at most n . Therefore, the time complexity is O n ( ) .

Theorem 5.4: The time complexity of the HD algorithm is PLS-complete with ρ m,n mO n ( ) ( ) .

Proof: We denote by ρ m,n ( ) the rounds needed for converging to the Nash equilibrium. Note that this is a random variable. In each step, only one mobile performs a maximization operation which causes totally mO n ( ) complexity. This is repeated during ρ m,n ( ) rounds. So, the time required can be given by T m,n ( ) = ρ m,n mO n ( ) ( ) . Actually, the complexity of finding the Nash equilibria in a potential game is known to be PLS-complete, refer to [29].

Furthermore, the computational results achieved in next section highlight that the number of rounds grows logarithmically with respect to the number of mobiles and BSs (Figure 8).

## VI. SIMULATION RESULTS

## A. Downlink System Model

For small-cells, the cellular network model consists of BSs arranged according to an homogeneous Poisson point process Φ of intensity λ b ( points/m 2 ) in the Euclidean plane [30]. For macro-cells, we use the classical honeycomb model to represent a well structured network made of large cells.

$$T ( m, n ) = \sum _ { s = 1 } ^ { m } \left [ ( n - 1 ) \left ( O ( m - s ) + T ( m - s, n - 1 ) \right ) + O ( n ) \right ] \text{ For small} \\ ( 4 9 ) \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{$$

(49)

where we assume that in each step s a sub-matrix occurs with complexity T m ( -s, n ) . Solving this kind of recurrence is tricky but further technical manipulations lead to following expression:

$$T ( m, n ) = O ( m ^ { 3 } n ) + ( n - 1 ) \sum _ { s = 1 } ^ { m } T ( m - s, n - 1 ) + m O ( n ). \begin{array} { c } \text{Poisson} \\ \text{in-med} \\ ( 50 ) & \text{non null} \end{array}$$

Also, we consider an independent collection of mobile users, located according to some independent homogeneous Poisson point process with intensity λ m ( points/m 2 ) . The main weakness of the Poisson model is that because of the independence of the Poisson point process, two BSs have a non null probability to be located very close to each other. This

weakness is balanced by two strengths: the natural inclusion of different cell sizes and shapes and the lack of edge effects, i.e. the network extends indefinitely in all directions [31]. The expected value of a homogeneous Poisson point process is E [Φ] = λA , where A ⊂ glyph[Rfractur] 2 denotes some area.

Moreover, the deployment scenario used to generate Figures 2, 3, 4, 6 corresponds to small-cells. We assume P r = -80 dBm , being the typical maximum received signal power of a wireless network as well as we set arbitrarily P ij = ∞ if P ij ≥ 20 dBm and we set the path loss exponent α = 3 . We also assume an equal operational power cost for all BS, P 0 = 12 W .

## B. Performance Results

We compare the proposed algorithms for different values of λ m and λ b . The average total power was calculated by Monte Carlo simulations by running the algorithms for different generated power cost matricesand taking the mean of the results.

We first start with macro-cells. As mentioned above, we use a honeycomb model with a cell radius equal to 4 km . Table I presents the optimal set-covering (SC) result as well as those obtained with the different proposed algorithms for different realization of power cost matrices and when the operational power cost is null, i.e. P 0 = 0 for all BSs. The transmission power P ij = ∞ if P ij ≥ 48 dBm . It turns out that the HD algorithm is very efficient and converges to nearly optimal assignments when operational power costs are neglected. This result indicates that switching off some BS does not decrease the total power significantly. The NBS algorithm also produces near optimal results in many examples. In these practical scenarios, the idea of switching off some BS is mostly interesting when the circuit power is dominant, which is not the case for macro-cells. Further, as fast moving users are associated with macro-cells in priority, it seems reasonable to keep all active. Therefore, the NBS algorithm is efficient and the gain achievable with any other optimal algorithm is marginal.

Table II compares the SC optimal results with all developed algorithms introduced in the paper for small-cells scenarios. The operational power cost is set to P 0 = 12 W and the BS density is increased. CC and DCC algorithms produce optimal assignments for almost all examples. However, the DCC algorithm naturally performs worse when the number of BSs increases. Moreover, the NBS algorithm exhibits worst results which highlights the interest of switching off some BSs when the operational power cost is higher than the transmission power cost.

Figure 2 illustrates the results achieved with the NBS algorithm, for different densities of BS and mobiles. This figure highlights an intuitive property of the NBS approach. When no operational cost is considered (see Figure 2), the power consumption decreases with the BS density, since the average distance between mobiles and BSs decreases accordingly. On the opposite, when an operational cost is considered, the NBS algorithm leads to an increased energy consumption since the number of active BS increases accordingly.

Let us now switch to the CC and DCC algorithms. Figure 3 focuses only on scenarios with circuit power, because they correspond to the more relevant cases. the upper curves show that the CC algorithm achieves much better results than NBS, especially when the BS density is high. This algorithm privileges solutions with larger cells. The distributed version DCC performs worst than the centralised one but still better than the NBS.

Figure 4 plots the change of the average total power with respect to the intensity of mobiles for small-cells scenario. The assumptions are as following: λ b = 1 11 . × 10 -5 points m 2 , θ = 0 002 . (in Figure 6, the optimal θ is found), and area A = 6 25 . km 2 . Note that the HD algorithm performs efficiently even though it is decentralized. For example, in case of λ m = 8 × 10 -8 , the average number of mobiles is given by 8 × 10 -8 · 6 25 . × 10 6 = 50 ; thus, the average power used per mobile is calculated as following: a) the HD algorithm: 63 50 = 1 26 / . W , b) the CC algorithm: 55 50 = 1 1 / . W , c) the greedy-SC algorithm: 50 50 = 1 / W , d) the SC algorithm: 43 83 . / 50 = 0 88 . W .

Figure 7 depicts the change of the average total power with respect to the intensity of mobiles for macro-cell deployment. λ b = 80 points 3600 km 2 , A = 3600 km 2 ( 60 km × 60 km area), and θ = 0 21 . (in Figure 7, we plot the change of average total power with respect to θ , and choose the optimal value). Here, we observe that the HD algorithm produces remarkable results. Calibrating θ properly is significant, otherwise the HD algorithm may not converge to the near optimal results. On the other hand, the NBS algorithm is also efficient in the macro-cell deployment. The drawback of greedy-SC algorithm reveals here since it works with a mechanism where the larger cells are privileged.

In Figures 6 and 7, the normalized average total power is plotted with respect to θ . From the figures and our observations in experiments performed in MATLAB, it might be considered that θ is mainly affected by the area over which the algorithm runs. For example, in Figure 7, the normalized average total power has a minimum in the same value of intensity of BSs, but it moves to a higher value when the area is enlarged from 2500 km 2 to 3600 km 2 .

Figure 8 shows the change of the average number of rounds of the HD algorithm for converging to a Nash equilibrium with respect to the area. The figure implies that the average number of rounds has a logarithmic characteristic. Moreover, when the operational power costs are zero, the average number of rounds increases since smaller cells are formed; therefore, the HD algorithm needs more rounds to converge to a Nash equilibrium.

## VII. CONCLUSION AND FUTURE WORKS

This paper addressed the MAP problem in the context of broadcast transmission. We introduced a novel decentralized solution based on group formation games, which we named the hedonic decision (HD) algorithm. This formalism is constructive: a new class of group formation games is introduced where the utility of players within a group is separable and symmetric being a generalization of party affiliation games. We proposed a centralized optimal recursive algorithm (the HM) as well as a centralized polynomial-time heuristic algorithm (the CC).

TABLE I A = 2500 km 2 , λ b = 6 points 2500 km 2 , λ m = 1 point 25 km 2 , P 0 = 0 W θ , = 0 11 . .

|   Ex. i |   m |   n |          |      HM |      SC |      CC |     DCC |     NBS |      HD |
|---------|-----|-----|----------|---------|---------|---------|---------|---------|---------|
|       1 |  54 |   6 | 10 - 7 × | 294.036 | 294.036 | 294.036 | 295.57  | 294.036 | 294.036 |
|       2 |  43 |   6 | 10 - 7 × | 299.816 | 299.816 | 323.619 | 352.166 | 299.816 | 299.816 |
|       3 |  46 |   6 | 10 - 7 × | 250.483 | 250.483 | 271.483 | 271.483 | 250.483 | 250.483 |
|       4 |  41 |   6 | 10 - 7 × | 270.442 | 270.442 | 302.815 | 287.968 | 284.574 | 283.074 |
|       5 |  51 |   6 | 10 - 7 × | 307.723 | 307.723 | 361.767 | 361.767 | 320.266 | 307.723 |
|       6 |  45 |   6 | 10 - 7 × | 278.521 | 278.521 | 317.509 | 317.509 | 278.521 | 278.521 |
|       7 |  41 |   6 | 10 - 7 × | 305.124 | 305.124 | 345.31  | 345.31  | 306.892 | 312.211 |
|       8 |  34 |   6 | 10 - 7 × | 221.568 | 221.568 | 236.717 | 236.717 | 256.074 | 221.568 |
|       9 |  58 |   6 | 10 - 7 × | 360.356 | 360.356 | 363.188 | 363.188 | 360.356 | 360.356 |
|      10 |  52 |   6 | 10 - 7 × | 310.072 | 310.072 | 310.072 | 310.072 | 332.074 | 332.074 |
|      11 |  44 |   6 | 10 - 7 × | 283.924 | 283.924 | 313.812 | 339.655 | 283.924 | 283.924 |
|      12 |  49 |   6 | 10 - 7 × | 312.972 | 312.972 | 312.972 | 312.972 | 325.58  | 337.994 |
|      13 |  53 |   6 | 10 - 7 × | 229.934 | 229.934 | 255.312 | 255.312 | 229.934 | 232.653 |
|      14 |  60 |   6 | 10 - 7 × | 308.6   | 308.6   | 308.6   | 308.6   | 320.206 | 320.206 |
|      15 |  57 |   6 | 10 - 7 × | 329.267 | 329.267 | 342.846 | 354.065 | 329.267 | 335.953 |

TABLE II A = 4 km 2 , λ b = 6 points 4 km 2 , λ m = 18 points 4 km 2 , P 0 = 12 W θ , = 0 003 . .

|   Ex. i |   m |   n |     HM |     SC |     CC |    DCC |    NBS |     HD |
|---------|-----|-----|--------|--------|--------|--------|--------|--------|
|       1 |  14 |   4 | 24.153 | 24.153 | 24.153 | 24.153 | 48.183 | 24.153 |
|       2 |  14 |   7 | 24.14  | 24.14  | 36.226 | 36.226 | 72.139 | 24.152 |
|       3 |  23 |   7 | 36.138 | 36.138 | 48.177 | 48.275 | 84.042 | 48.185 |
|       4 |  16 |   7 | 36.195 | 36.195 | 36.195 | 36.257 | 72.208 | 48.204 |
|       5 |  19 |   4 | 24.109 | 24.109 | 24.109 | 36.174 | 48.126 | 36.184 |
|       6 |   7 |   3 | 12.07  | 12.07  | 12.07  | 12.07  | 24.073 | 12.07  |
|       7 |  15 |   8 | 36.178 | 36.178 | 36.178 | 36.222 | 84.241 | 36.222 |
|       8 |  15 |   7 | 24.11  | 24.11  | 24.11  | 24.144 | 72.101 | 24.154 |
|       9 |  18 |   9 | 36.099 | 36.099 | 36.148 | 36.253 | 84.042 | 36.148 |
|      10 |  24 |   6 | 36.111 | 36.111 | 36.191 | 36.191 | 60.134 | 36.191 |
|      11 |  21 |   3 | 24.149 | 24.149 | 24.149 | 24.149 | 36.162 | 24.149 |
|      12 |  18 |  10 | 36.111 | 36.111 | 36.111 | 48.141 | 84.053 | 48.152 |
|      13 |  23 |   5 | 36.105 | 36.105 | 36.105 | 48.143 | 60.218 | 36.105 |
|      14 |  13 |   4 | 24.082 | 24.082 | 24.104 | 24.104 | 48.123 | 24.11  |
|      15 |  17 |   8 | 36.19  | 36.19  | 36.19  | 36.217 | 60.217 | 36.19  |

TABLE III λ b = 0 10 . × 10 -3 points m 2 , λ m = 1 11 . × 10 -3 points m 2 , P 0 = 12 W , θ = 0 008 . .

|                      | Number of rounds   | Number of rounds   | Number of rounds   | Number of rounds   |
|----------------------|--------------------|--------------------|--------------------|--------------------|
| Example/Area( km 2 ) | 0 . 98             | 1 . 28             | 1 . 62             | 2 . 00             |
| 1                    | 3                  | 2                  | 3                  | 3                  |
| 2                    | 3                  | 3                  | 3                  | 3                  |
| 3                    | 3                  | 3                  | 2                  | 3                  |
| 4                    | 2                  | 3                  | 3                  | 3                  |
| 5                    | 3                  | 2                  | 3                  | 3                  |
| 6                    | 3                  | 3                  | 3                  | 3                  |
| 7                    | 3                  | 3                  | 3                  | 4                  |
| 8                    | 3                  | 3                  | 3                  | 4                  |
| 9                    | 3                  | 3                  | 3                  | 3                  |
| 10                   | 3                  | 3                  | 3                  | 4                  |

The results exhibit that the HD algorithm achieves very good results if the parameter θ is well chosen. The exact value of θ is not provided and may be used as a setting parameter.

The proposed HD algorithm is efficient and may be used in many other set covering problems. For instance, indoor wireless network planning has been studied for several years and the optimal BS activation is an important problem. The proposed HD algorithm could be used for planning purposes and thus optimizing the number of BS, but could also be used to optimize dynamically the number of active BS as a function of actives users. The provider may deploy a high density of BS and then could run dynamically the HD algorithm to optimize

Fig. 2. The NBS Algorithm: Average total power ¯ p with respect to intensity of BSs λ b for increasing values of intensity of mobiles λ m = (0 16 . , 0 25 0 44 1 00) . , . , . × 10 -4 points m 2 , A = 4 km 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-603aa663-0143-47a1-9ad9-8c840ed74bc0.jpg)

then number of active BSs.

Furthermore, it could be interesting to change the game parameters. One can define different clustering weights for each BS. That setting could provide better results. It is also possible to try different symmetric bipartite utility allocations.

## REFERENCES

- [1] S. Sesia, I. Toufik, and M. Baker, LTE The UMTS Long Term Evolution: From Theory to Practice . John Wiley &amp; Sons, Ltd., 2009.
- [2] J. Krarup, P.M. Pruzan, 'The simple plant location problem: Survey and synthesis,' European Journal of Operational Research , vol. 12, pp. 3681, 1983.
- [3] V. Bil`, I. Caragiannis, C. Kaklamanis, and P. Kanellopoulos, 'Geometric o clustering to minimize the sum of cluster sizes,' In Proc. 13th European Symp. Algorithms, vol.3669 of LNCS , pp.460-471, 2005.
- [4] N. Lev-Tov and D. Peleg, 'Polynomial time approximation schemes for base station coverage with minimum total radii,' Computer Networks , vol.47, no.4, pp.489-501, Mar. 2005.

Fig. 3. Average total power ¯ p with respect to intensity of BSs λ b for increasing values of intensity of mobiles λ m = (0 16 0 25 0 44 1 00) . , . , . , . × 10 -4 points m 2 , A = 4 km 2 , P 0 = 12 W .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-541fa63b-0c91-4738-afee-1daec262d823.jpg)

Fig. 4. Small cells: Change of the average total power ¯ p with respect to the intensity of mobiles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-26baa5f1-a584-4778-96cc-9aa91b83bf47.jpg)

p

Fig. 5. Macro cells: Average total power ¯ p with respect to the intensity of mobiles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-61bef383-0101-4bca-b3e3-343bec7e67b6.jpg)

Fig. 6. Small cells: Normalized average total power with respect to θ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-11076de0-3f3d-4482-9ed8-cddf97ad2237.jpg)

- [5] H. Alt, E. M. Arkin, H. Br¨ onnimann, J. Erickson, S. P. Fekete, C. Knauer, J. Lenchner, J. S. B. Mitchell, and K. Whittlesey, 'Minimum-cost coverage of point sets by disks,' In Proceedings of ACM Symposium on Computational Geometry (SCG '06) , pp. 449-458, New York, NY, USA, 2006.
- [6] S. Funke, S. Laue, Z. Lotker, and R. Naujoks, 'Power assignment problems in wireless communication: Covering points by disks, reaching few receivers quickly, and energy-efficient travelling salesman tours,' In Proceedings of Ad Hoc Networks , pp.1028-1035, 2011.
- [7] V. Rodoplu, T.H. Meng, 'Minimum energy mobile wireless networks,' IEEE Journal on Selected Areas in Communications , vol.17, no.8, pp.1333-1344, Aug. 1999.
- [8] J. E. Wieselthier, G. D. Nguyen, and A. Ephremides, 'On the construction of energy-efficient broadcast and multicast trees in wireless networks,'.

1

point

θ

Fig. 7. Macro cells: Normalized average total power with respect to θ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-b3766489-ece3-4d45-9456-e1d250991c78.jpg)

Fig. 8. Average number of rounds with respect to the area.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0482v1-figure-e04a7551-bbf9-4b43-a5c7-b229b7bc1906.jpg)

IEEE INFOCOM 2000 , pp.586-594, Tel Aviv, Israel, 2000.

- [9] ¨ . O E˘ gecio˘ glu, and T. Gonzalez, 'Minimum-energy broadcast in simple graphs with limited node power,' IASTED International Conference on Parallel and Distributed Computing and Systems (PDCS 2001) , pp.334338, Anaheim, CA, Aug. 2001.
- [10] M. Cagalj, J. P. Hubaux, and C. Enz, 'Minimum-energy broadcast in all-wireless networks: NP-completeness and distribution issues,' In Proceedings of ACM international conference on Mobile computing and networking (MobiCom) , New York, NY, USA, pp.172-182, 2002.
- [11] C. Hasan, E. Altman, and J.M. Gorce, 'A coalition game approach to the association problem of mobiles in broadcast transmission,' WiOpt 2011 pp.236-240, 9-13 May 2011.
- [12] A. P. Bianzino, C. Chaudet, D. Rossi, J. L. Rougier, 'A survey of green networking research,' IEEE Communications Surveys and Tutorials , vol.14, no.1, pp.3-20, 2012.
- [13] H. Claussen, 'The future of small cell networks,' IEEE COMMSOC MMTC E-Lett. [Online], pp.32-36, Sept. 2010. Available:

http://committees.comsoc.org/mmc/e-news/E-Letter-September10.pdf

- [14] E. Balas and M. W. Padberg, 'Set partitioning: A survey,' SIAM Review , vol.18, pp.710-760, 1976.
- [15] C. Lund, M. Yannakakis, 'On the hardness of approximating minimization problems,' Journal of the ACM vol.41, no.5, pp.960-981, Sept. 1994.
- [16] C. K. Singh and E. Altman, 'The wireless multicast coalition game and the non-cooperative association problem,' in 30th IEEE INFOCOM , Shanghai, China, 10-15, Apr. 2011.
- [17] Chv´ atal, V., 'A greedy heuristic for the set-covering problem,' Mathematics of Operations Research , vol. 4, no. 3, pp. 233-235, 1979.
- [18] G. Hollard, 'On the existence of a pure strategy Nash equilibrium in group formation games,' Elsevier Economics Letters , vol. 66, pp. 283287, 2000.
- [19] Konishi, H., Le Breton, M., Weber, S.,'Pure strategy Nash equilibrium in a group formation game with positive externalities,' Games and Economic Behavior vol. 21, 161-182.
- [20] I. Milchtaich, 'Stability and segregation in group formation,' Games and Economic Behavior , vol. 38, pp. 318-346, 2002.
- [21] D. Monderer and L. S. Shapley, 'Potential games,' Games and Economic Behavior vol. 14, pp. 124-143, 1996.
- [22] R. W. Rosenthal, 'A class of games possessing pure-strategy Nash equilibria,' International Journal of Game Theory vol. 2, 65-67, 1973.
- [23] J. H. Aldrich, W. T. Bianco, 'A game-theoretic model of party affiliation of candidates and office holders,' Mathematical and Computer Modelling , vol. 16, no. 8-9, pp. 103-116, Aug.-Sept., 1992.
- [24] D. Leibovic and E. Willett, 'Selfish set covering,' In Midstates Conference for Undergraduate Research in Computer Science and Mathematics , Nov. 2009.
- [25] X.-Y. Li, Z. Sun, W. Wang, 'Cost sharing and strategyproof mechanisms for set cover games,' Lecture Notes in Computer Science , vol. 3404, pp. 218-230, 2005.
- [26] X.-Y. Li, Z. Sun, W. Wang, X. Chu, S. J. Tang, P. Xu, 'Mechanism design for set cover games with selfish element agents,' Theoretical Computer Science , vol. 411, no. 1, pp. 174-187, 2010.
- [27] M. Blum, R. W. Floyd, V. Pratt, R. L. Rivest, and R. E. Tarjan, 'Time bounds for selection,' J. Comput. Syst. Sci. vol.7, no.4, pp.448-461, Aug. 1973.
- [28] T. Feder, and D. Greene, 'Optimal algorithms for approximate clustering,' Proceedings of the twentieth annual ACM symposium on Theory of computing , pp. 434-444, 1988.
- [29] A. Fabrikant, C. Papadimitriou, and K. Talwar, 'The complexity of pure Nash equilibria,' In Proceedings of the thirty-sixth annual ACM symposium on Theory of computing (STOC '04) , New York, NY, USA, 604-612, 2004.
- [30] F. Baccelli, and B. Błaszczyszyn, 'Stochastic geometry and wireless networks: volume I theory,' Foundations and Trends in Networking , vol.3, pp.249-449, 2009.
- [31] J. G. Andrews, F. Bacelli, and R. K. Ganti, 'A tractable approach to coverage and rate in cellular networks,' IEEE Trans. on Commun. , vol.59, no.11, pp.3122-3134, Nov. 2011.
- [32] J. Hoydis, M. Debbah, 'Green, cost-effective, flexible, small cell networks,' IEEE ComSoc MMTC E-Letter Special Issue on 'Multimedia Over Femto Cells' , 2010.
Cengis Hasan Jean-Marie Gorce Eitan Altman
Green Broadcast Transmission in Cellular Networks: A Game Theoretic Approach
This paper addresses the mobile assignment problem in a multi-cell broadcast
transmission seeking minimal total power consumption by considering both
transmission and operational powers. While the large scale nature of the
problem entails to find distributed solutions, game theory appears to be a
natural tool. We propose a novel distributed algorithm based on group formation
games, called \textit{the hedonic decision algorithm}. This formalism is
constructive: a new class of group formation games is introduced where the
utility of players within a group is separable and symmetric being a
generalized version of parity-affiliation games. The proposed hedonic decision
algorithm is also suitable for any set-covering problem. To evaluate the
performance of our algorithm, we propose other approaches to which our
algorithm is compared. We first develop a centralized recursive algorithm
called \textit{the hold minimum} being able to find the optimal assignments.
However, because of the NP-hard complexity of the mobile assignment problem, we
propose a centralized polynomial-time heuristic algorithm called \textit{the
column control} producing near-optimal solutions when the operational power
costs of base stations are taken into account. Starting from this efficient
centralized approach, a \textit{distributed column control algorithm} is also
proposed and compared to \textit{the hedonic decision algorithm}. We also
implement the nearest base station algorithm which is very simple and intuitive
and efficiently manage fast-moving users served by macro BSs. Extensive
simulation results are provided and highlight the relative performance of these
algorithms. The simulated scenarios are done according to Poisson point
processes for both mobiles and base stations.