The concept of optimal planning of a linearly oriented segment of the 5G network

In the article, the extreme problem of finding the optimal placement plan of 5G base stations at certain points within a linear area of finite length is set. A fundamental feature of the author’s formulation of the extreme problem is that it takes into account not only the points of potential placement of base stations but also the possibility of selecting instances of stations to be placed at a specific point from a defined excess set, as well as the aspect of inseparable interaction of placed 5G base stations within the framework of SON. The formulation of this extreme problem is brought to the form of a specific combinatorial model. The article proposes an adapted branch-and-bounds method, which allows the process of synthesis of the architecture of a linearly oriented segment of a 5G network to select the best options for the placement of base stations for further evaluation of the received placement plans in the metric of defined performance indicators. As the final stage of the synthesis of the optimal plan of a linearly oriented wireless network segment based on the sequence of the best placements, it is proposed to expand the parametric space of the design task due to the specific technical parameters characteristic of the 5G platform. The article presents a numerical example of solving an instance of the corresponding extremal problem. It is shown that the presented mathematical apparatus allows for the formation of a set of optimal placements taking into account the size of the non-coverage of the target area. To calculate this characteristic parameter, both exact and two approximate approaches are formalized. The results of the experiment showed that for high-dimensional problems, the approximate approach allows for reducing the computational complexity of implementing the adapted branch-and-bounds method by more than six times, with a slight loss of accuracy of the optimal solution. The structure of the article includes Section 1 (introduction and state-of-the-art), Section 2 (statement of the research, proposed models and methods devoted to the research topic), Section 3 (numerical experiment and analysis of results), and Section 4 (conclusions and further research).


Introduction and state-of-the-art
The strategy for automating the control processes of wireless communication networks is implemented by 3GPP (the main developer of the Technical Specifications for 5G+ generation mobile communications) in the form of a technological concept for the use of Self-Organizing Networks (SON) [1][2][3][4][5], covering all stages of the life cycle of the corresponding communication networks.As defined in the 3GPP technical specifications, SON is a network used by a dedicated control and operational system to automate the management of a 5G network with minimal human intervention.The implementation of SON algorithms in mobile communication networks of the 5G generation is primarily focused on achieving such goals as [6][7][8]: • reduction of capital and operating costs of operators (CAPEX/OPEX) during deployment, operation and optimization of the network; • increasing spectral efficiency, productivity and capacity, improving radio coverage, and increasing the efficiency of resource use of Next Generation Radio Access Networks (NGRAN).
In the context of its intended purpose, SON is considered a driver technology for the development of 5G networks.However, malfunctioning SON network algorithms (for example, after updating cell coverage areas and when connecting new base stations) can pose a threat to the normal functioning of the wireless network and have a significant negative impact, reducing the quality of services of connected users.Considering that SON algorithms are not standardized by 3GPP technical specifications, research into issues related to achieving the abovementioned intended purpose of SON is relevant.
Let's focus our attention on designing the architecture of a 5G network segment with SON, which is deployed over a linear area of finite length.We chose this topology not by chance because the linear organization of space is typical both for urban infrastructure and for intercity communications (highways, railways).When designing a target linearly oriented segment of a 5G network, it is necessary to solve several interrelated problems, namely [9,10]: • conduct an analysis and preliminary study of the network deployment location, • select technical means and protocols, • determine the architecture of base station placement and evaluate the performance characteristics of the 5G network segment using mathematical modelling.
An effective way to increase the technical and economic indicators of such a project is to optimize the architecture of the created 5G network segment, namely, solving the problem of selecting the optimal set of stations from a given redundant set and determining their locations along the target linearly oriented aria.
The challenge of determining the optimal locations for 5G base stations can be conceptualized as a network design problem involving relays (NDPR), a subject that has garnered significant attention in recent research [11][12][13][14][15][16][17][18][19].Originating from the work of Cabral et al. in [15], this problem aims to minimize costs while ensuring that each source node can reach the sink node through a single route within limited transmission distances.Konak's extension [14] introduces a novel problem formulation incorporating a set covering constraints.This formulation forms the basis for a genetic algorithm featuring a specialized crossover/mutation operator, generating feasible paths for each commodity.The relay locations on these paths are then determined by solving the corresponding set covering problem.In [15], Xiao et al. propose a hybrid approach utilizing variable neighbourhood searches.The variable neighbourhood algorithm identifies routes for each source node, and an implicit enumeration algorithm determines optimal relay locations for a given set of routes.Li et al. [12] put forward an iterative meta-heuristic algorithm based on tabu search, addressing the NDPR problem in two steps.However, these studies primarily focus on transmission distance constraints and do not fully account for the operational dynamics of wireless sensor networks.In contrast, Bagaa et al. [18] and Li et al. [19] establish models for optimal relay node placement, specifically targeting the minimization of energy consumption in wireless signal transmission.Various works in the literature [20,21] delve into relay node location problems considering signal coverage constraints.Naveen et al. [21] define signal coverage as a location within the communication range of a node eventually connected, possibly via multiple hops, to a sink node.Existing models studying signal coverage often rely on the ideal 0-1 disk channel model [22] and the latest Signal-to-Noise Ratio (SNR) model [16].
Nevertheless, these investigations predominantly concentrate on ascertaining whether all designated areas can be covered, with limited attention given to service quality and timely reliability [23].In the realm of wireless networks, the paramount concern is whether user needs can be satisfied within specified timeframes.Hence, it becomes crucial to incorporate timely reliability considerations when modelling the location of base stations.Notably, Nigam et al. [11] and Naveen et al. [21] have delved into the optimal placement of relay nodes in wireless sensor networks, taking delay constraints into account.Their objective is to formulate a multihop wireless network that meets transmission path delays, ultimately minimizing the number of required relay base stations.However, these investigations uniformly employ the number of hops as the sole criterion for measuring network transmission delay and reliability.This approach overlooks the intricacies of the actual data transmission processand user requirements.Consequently, the results obtained from these studies may not directly offer actionable insights for communication operators.
One of the major challenges faced by current heuristic-based models in addressing the optimization of 5G base station locations is the accurate simulation of Line-of-Sight (LOS) propagation and the coverage of 5G signals in urban settings.Many existing location optimization models for cellular network planning were originally designed for 2G/3G/4G networks [24].In urban areas, the penetration loss of 2G/3G/4G signals is considerably lower than that of 5G signals [25].Consequently, the location optimization models developed for earlier generations of cellular networks did not need to account for the simulation of LOS propagation of signals.
Recent research endeavours have employed Geographic Information Systems (GIS) to facilitate the simulation of Line-of-Sight (LOS) propagation and coverage of 5G signals, offering support for optimization and enabling intuitive visual analyses in the planning of 5G networks.For instance, [26] utilized GIS alongside a computational geometry approach to optimize the locations of below-rooftop wall-mounted 5G base stations in urban outdoor areas, with GIS providing crucial support [27].However, it is worth noting that the computational geometry approach may not be suitable for global optimization.Conversely, recent studies have introduced spatially implicit heuristic algorithm-based models to aid in 5G cellular network planning, particularly in specialized 5G application scenarios [28].For example, the deployment of 5G base stations on unmanned aerial vehicles (UAVs) has been explored as an effective means to mitigate the penetration loss of 5G signals in urban areas.Consequently, these optimization models for UAV base station deployment do not involve the simulation of LOS propagation of 5G signals.
The papers [29][30][31][32] discuss the strategic planning and implementation of wireless heterogeneous networks for smart metering, employing a cross-layer approach.Our methodology integrates considerations from the network layer, addressing routing and flow demands at each link within the network, while simultaneously accounting for limitations at the physical layer, specifically the capacity constraints of short-range technology in a multi-hop configuration.The proposed model utilizes a column generation approach to address the capacitated multicommodity flow problem, incorporating wireless link capacities, coverage parameters, and associated costs.Those studies incorporate multi-hop routing of packets within a mesh network composed of smart meters and concentrators connected to a cellular network through base stations.Each link's traffic is represented in a multigraph with occupation percentages, and a scalable routing tree is constructed on a georeferenced map to simulate a real-world deployment.The outcomes illustrate the performance of the proposed model concerning concentrator traffic load, network coverage, and energy consumption reduction.The findings underscore that the integration of multi-hop short-range technology leads to a cost-effective infrastructure by diminishing the number of smart meters requiring a direct connection to cellular technology.
Unlike the previously mentioned studies, our research focuses on a specific optimization problem formulated in the context of organizing a wireless network based on the 5G platform.The specificity is that the mentioned network is deployed along a linear segment, for example, a road, pipeline, or railway.A special feature of the 5G platform is, among other things, the ability to implement several virtual network segments on the base station hardware using Network Slicing technology.In the model proposed in Section 2.1, it is assumed that the frequency resource of the base stations (the cost of which is known) is divided into two virtual segments, one of which is used to serve subscribers, and the second to ensure the network connectivity (for communication with other base stations).The first segment is characterized by higher throughput and a smaller coverage radius, and the second-vice versa.For the deployed network, parameters such as total cost and total delay are taken into account.Within the target linear segment, the coordinates of the points at which the base station can be located are known.The special same type of base stations (gateways) are located at the ends of the target segment.The optimization problem is that it is necessary to place base stations on the target segment in such a way as to maximize the size of the telecommunications coverage of the segment while fulfilling the requirement for each station to communicate with gateways through a system of hosted stations, as well as fulfilling restrictions on the amount of end-to-end delay and the total cost of deployed stations.Thus, the wireless network structure obtained as a result of solving such an optimization problem will be reliable (due to constant communication between base stations and gateways) and efficient (due to maximizing coverage focused on serving subscribers).In addition, the gateway concept allows the integration of multiple line segments into a single network.This allows us to limit ourselves to the analysis of dozens of points of potential location of base stations in the process of searching for the optimal network configuration.In turn, this makes it possible to solve such an optimization problem using exact methods, among which the most suitable is the branch and bound method (subject to its adaptation to the above-formulated specifics of setting the optimization problem).In addition, the proposed mathematical apparatus takes into account that base stations (not gateways) can be of different types (this allows us to take into account the possibility of upgrading already functioning network architectures).Naturally, the last clarification further complicates the final optimization problem.All of the above stated in the complex is a new research problem, as a result of which it is possible to design linearly organized network segments based on 5G base stations that meet the conditions of high reliability, high quality of service and rational use of funds.
Characteristic of the mentioned studies is the variety of methods used to solve the problems of optimal design of wireless network architectures, which does not eliminate the general drawback: the authors deliberately narrow the range of characteristics of the deployed base stations taken into account (most often, to two parameters-coverage semidiameter and cost).The authors also consider the characteristics of the placed base stations to be unified, ignoring that the actual values of the corresponding characteristics depend on the location of the device.At the same time, the aspect of network connectivity (the continuity of interaction between located base stations within the SON) is not considered.The motive of the authors is clear-they strive to move from multi-parameter optimization problems to single-parameter ones.This makes it possible to increase the probability of the existence of a feasible solutions region and simplifies the solution process as such, but makes the resulting models extremelyabstracted from the real object.
Note that the final stage of designing a wireless network segment is estimating its performance characteristics to verify the compliance of the values of these assessments with the specified requirements [9,10].This estimation most often uses simulation modelling, which is a computationally intensive process.If the network characteristics do not meet the specified requirements, we have to return to the previous stages and repeat the iterative design process with new input data to obtain a new architecture for the wireless network segment.
Next, we will present the main elements that determine the orientation, structure and novelty of the presented research.
The object of the research is the process of designing the architecture of a linearly oriented segment of the 5G network in the parametric space of controlled qualitative parameters with constraints.
The research subject includes the theory and methods of mathematical modelling, functional analysis, combinatorics and operations research.
The purpose of the research is formulated as follows: to create a computationally efficient, holistic concept for finding the optimal 5G base station placement plan at specified points within a linear area of finite length.
The tasks of the research will be defined as follows: • present the research object in the form of an extremal problem on a finite set in the form of a specialized combinatorial model, • to solve the formulated problem, adopt the branch-and-bounds method, which consists of taking into account the specifics of the object of research, presented in the form of a corresponding mathematical model, • to formulate variants of the procedure for finding a sequence of optimal placement plans of base 5G stations, taking into account the technological characteristics of the latter ones, which differ in terms of computational complexity, • implement a numerical experiment, the results of which will justify the effectiveness of the created mathematical apparatus.Now let's formulate the main contribution of the research.In the article, the extreme problem of finding the optimal placement plan of 5G base stations at certain points within a linear area of finite length is set.A fundamental feature of the author's formulation of the extreme problem is that it takes into account not only the points of potential placement of base stations but also the possibility of selecting instances of stations to be placed at a specific point from a defined excess set, as well as the aspect of inseparable interaction of placed 5G base stations within the framework of SON.The formulation of this extreme problem is brought to the form of a specific combinatorial model.The article proposes an adapted branch-and-bounds method, which allows the process of synthesis of the architecture of a linearly oriented segment of a 5G network to select the best options for the placement of base stations for further evaluation of the received plans in the metric of defined performance indicators.As the final stage of the synthesis of the optimal plan of a linearly oriented wireless network segment based on the sequence of the best placements, it is proposed to expand the parametric space of the design task due to the specific technical parameters characteristic of the 5G platform.

Statement of the research
Suppose that there is a set of multiband base stations B = {b i }, i ¼ 1; n, each element of which is characterized by parameters where s i is the coverage semidiameter of the service supported by the base station b i , focused on users' tasks maintenance; Sij The maximum amount of funds that can be spent on the deployment of the network architecture B is represented by the value of the parameter V.The maximum permissible delay in the transmission of a data packet between the users of the network B is represented by the value of the parameter D.
The investigated network has a linear topology, which is deployed on a segment ρ = [p 0 , p m +1 ] of length l, where p 0 and p m+1 are the starting and ending points of the segment ρ, respectively.Within the segment ρ, a set of points is defined that are admissible for placing base stations of the network B: The point p k corresponds to the coordinate x k within the segment ρ, 8k.The coordinates of the start and end points of the segment ρ are defined as x 0 = x 0 (p 0 ) = 0 and x m+1 = x m+1 (p m+1 ) = l, respectively.
The logic of forming such a network infrastructure assumes that specialized base stations b 0 and b n+1 are located at points p 0 and p m+1 , which are not included in the set of typical base stations B: {b 0 , b n+1 } = 2 B, and act as gateways.At the same time, the coverage semidiameters of specialized base stations {b 0 , b n+1 } for connection with base stations from the set B are represented by the values of the parameters S 0i and S (n+1)i , respectively.The costs of specialized base stations {b 0 , b n+1 } are not taken into account when forming the value of the parameter V for an arbitrary implementation of the network B.
We need to place base stations from a set of B in admissible points P at the segment ρ of length l in such a way as to maximize the aggregate communication coverage of the latter.With: • each base station from the set B must constantly maintain communication with gateways {b 0 , b n+1 } (directly or through other base stations of the same type), • the cost of implementing the found base stations B placement plan should not exceed the value of V, • the maximum allowable delay of data packet transmission between users within the coverage provided by the found base stations B placement plan should not exceed the value of D.
Let's generalize the formulated requirements in the form of an extremal problem on a finite set.An acceptable placement plan of base stations from the set B is the following set of pairs increasing in value of the coordinate x k , for which the following requirements are met: 1.For each pair (p k , b i ): • on the left: either there is a pair of (p q , b j ) such that x k − x q � S ij and x k − x q � S ji , or x k − x 0 � S i0 and x k − x 0 � S 0i ; • on the right: either there is a pair of (p g , b f ) such that 2. At each point from the set P, no more than one base station from the set B can be located; 3. The inequality , where B δ is the set of deployed base stations from the set B, and D i is the average data packet transmission delay determined for the base stations from the set B δ ; 4. The total cost of base stations from the set B δ does not exceed the value of V.
Note that the fulfilment of requirement 1 guarantees that each base station from the set B δ will communicate with the gateways {b 0 , b n+1 } either directly or through other stations of this set.
Let's match each acceptable placement plan A with the corresponding coverage value y(A), which is interpreted as the total coverage for all base stations from the set B included in the placement plan A. In this context, we define the concept of non-coverage of a segment ρ as h We denote the set of all admissible placement plans A by W. Accordingly, let us present the extremal problem described above in the combinatorial form: Finally, we denote the entire set of options for placing base stations from the set B on the given set P as O. Wherein, B δ 2 O.

Representation of the set O in the form of a branching tree
We will focus the material of Section 2 on the formalization of the solution of the extremal problem (1) by the branch-and-bounds method [33].In this context, it is necessary to determine the procedure for forming a branching tree (binary search tree) for the set O to which the root vertex corresponds.The basis of the procedure we are looking for is the rule of dividing the original set into subsets.Each subset will correspond to a vertex on the branching tree.The tree itself is oriented relative to the root vertex.So, let's determine the rule for dividing the original set into subsets and the procedure for moving along the branching tree.
Let's redefine the set O as W c = 0 , where the subscript c fixes the iteration number.Starting from c = 0, at each c-th iteration we will split the original subset W c into two derived subsets are fulfilled for derived subsets (look at the Fig 1).As a rule that determines the division, we introduce a binary value β ki .Let β ki = 1, if the base station b i 2 B will be located at the point p k .Accordingly, β ki = 0 if no base station from the set B is located at point p k .Next, we will associate the partitioning rule β ki = 1 with a subset of W 1 c , and we will associate the partitioning rule β ki = 0 with a subset of W 2 c .
At the stage of partitioning any initial set W c , the entire set of partitioning rules B = {β ki } can be divided into three subsets: • the subset B + , which includes all cases when β ki = 1, • the subset B − , which includes all cases when β ki = 0, • the subset B ~, which includes all cases that are not affected by the value of β ki at the current c-th iteration.
To split the original set W c , at each iteration, a variable from the subset B ~with the smallest index i among all variables with the smallest index k is selected.Thus, a free location point p k with the smallest index k is first selected and the option of placing a base station b i 2 B with the smallest index i at this point is investigated.
After splitting the next source set W c into derived subsets W 1 c and W 2 c (Fig 2 ), the latter ones on the branching tree are redefined as W c+1 and W c+2 , respectively.This movement into the depth of the branching tree is typical.At the same time, if an empty set W c or a set consisting of one element is obtained as a result of the next partition, then such a vertex of the tree is considered closed, and after its analysis, a reverse transition is made from the closed set W c to one of the previously formed subsets.
To solve the extremal problem (1), we will use the LIFO rule to move through the branching tree.According to this rule, we will implement movement in depth until we find a closed vertex.At the same time, of the two sets W 1 c and W 2 c , the set W 1 c will be the first to be examined for the possibility of closing the corresponding vertex.If the found vertex does not turn out to be closed, then movement in depth along the same branch will be implemented from it.Otherwise, a reverse transition will be made: the last formed vertex from the unconsidered ones (the unclosed vertex with the largest sequence number c) will be selected to continue the movement.The LIFO branching procedure will end when all vertices of the tree are closed.Note that the above-defined vertex selection rule for the implementation of movement in depth gives any set W c an important property.Suppose that for the studied set W c , c 2 N, a base station from the set B has already been placed at the point p q following the partitioning rule β ki , which separated the derived set W c from the original set W 0 .Then, for all points p k located to the right of the point p q : k < q, the location of base stations from the set B has already been determined (at the same time, some of these points may not contain base stations).

Adaptation of the branch-and-bounds method to the research object
The adaptation of the branch-and-bounds method for solving the extremal problem (1) taking into account the procedure of movement along the branching tree described in Section 2.2 consists in the formalization of the approach to the study of the vertices of the tree for the possibility of their closure.
Let's examine the vertex that corresponds to the set W c .According to the branch-andbounds method, the vertex can be closed in the following cases: Case 1.The set W c is empty, i.e. it is proved that there is no admissible placement plan A in the set W c with the specified branching rule β ki ; Case 2. It is proved that the set W c contains an admissible placement plan A, which corresponds to a smaller value of the objective function (1) than the best version of the placement plan Â among those already found.The value of the non-coverage function h Â À � will be considered extreme, and the placement plan Â≔A will be considered optimal.We consider , we take the length l of the segment ρ; Case 3. The optimum Â of the extremal problem (1) on the set W c is found.
Let's examine these cases in more detail.
Case 1. Checking whether the current derived set W c is empty is based on requirements 1-4 defined at the end of Section 2.1.
Consider the conditions under which requirement 1 for the set W c , c 2 N, is fulfilled.Let the set W c be formed as a result of splitting the original set using rule β qg = 1.We check whether each of the semidiameters S gz and S zg is greater than the distance l q − l d , where z is the index of the base station from the set B located at the point closest to the left p d .If the point closest to the left is p 0 , then the check is made for semidiameters S q0 and S 0q .If the described checks do not pass, then the set W c is incompatible, the vertex corresponding to it is considered closed, and the reverse transition is performed according to the rule for the formation of derived subsets defined at the beginning of Section 2.2.
If the set W c is formed as a result of partitioning the original set W 0 using the rule β qg = 0 and p d is the point with the largest index among the points from the set P, in which the base stations from the set B are already placed (if there are no base stations placed, then p d = p 0 ), then it should be checked whether among the undistributed base stations from the set B (excluding b g ) there is such a base station b j for which the distance between the points p q and p d does not exceed S jz and S zj .If the described check is not passed, then the set W c is incompatible, the vertex corresponding to it is considered closed and the reverse transition is carried out according to the rule for the formation of derived subsets defined at the beginning in Section 2.2.
Requirement 2 is fulfilled as a result of the placement of the next base station from the set B. Requirements 3 and 4 are fulfilled as a result of adding the values and delays characteristic of the base stations from the set B located on the segment ρ.
Case 2. Let's estimate the amount of non-coverage of the segment ρ of length l for the derived set W c obtained from the original set when applying the partitioning rule β qg = 1.For this, we introduce the expression where γ is the partial non-coverage function, which is defined for any two points p q and p d from the set P at q > d, in which the base stations b t and b g are located, provided that there are no other base stations between these points.For an arbitrary location plan A, the non-coverage function h(A) is defined as the sum of the values of functions ( 4) determined for all base stations located according to the placement plan A, including gateways {b 0 , b n+1 }.We define the lower bound of the non-coverage function h(A) for an arbitrary placement plan A, which corresponds to the set W c : holds, then the set W c cannot contain a placement plan A better than the already found Â.Therefore, the corresponding set W c on the search tree is closed and the reverse transition is performed according to the rule for the formation of derived subsets defined at the beginning in Section 2.2.
In the context of ( 5), the estimate of the non-covering function for the derived set W c obtained from the original set when applying the partitioning rule β qg = 1 is represented by the expression where the term r 1 (W c ) is calculated as the sum of the values of functions (4) to the left of the point p q and the value of the coverage semidiameter of the base station b g located at this point.The term r 2 (W c ) is calculated as the sum of the values of functions (4) to the right of the point p q , that is, on some fragment π to the end of the segment ρ (that is, to the point p m+1 ).We will determine the estimate r 2 (W c ) by relaxing the conditions that characterize the possible location of the base stations from the set B on the fragment π.For this, we need, taking into account the requirements 1-4 defined at the end of Section 2.1, to define a subset B π 2 B that contains not yet placed base stations and ensures minimal non-coverage on the fragment π.We will determine the subset Bπ as a result of solving the corresponding Boolean programming problem.Let's introduce Boolean variables u i ¼ 18b i 2 B m ; 08else: ( .We formulate the desired Boolean optimization problem in the parametric space of controlled variables u i 2 {0, 1} as follows: where |π| is the length of the fragment π and n is the number of unoccupied points for placing base stations from the set B on the fragment π.Note that if constraint (8) or constraint ( 9) is neglected, the optimization problem (7) converges to the integer Knapsack Problem, for the solution of which there is an efficient pseudopolynomial method [34].
The effectiveness of applying assessment in the branch-and-bounds method is determined by its accuracy and the timing of its application.There is no universal solution method for the integer linear programming problem (7).However, based on the above observation, there is a potential opportunity to formalize two less precise estimates, which can be calculated relatively easily.At the same time, the more realistic option is to ignore constraint (9), because, in reality, the number of points for placing base stations significantly exceeds their available number.Next, we will call problem (7) free of constraint ( 9) reduced.Also, the removal of the restriction on integers inherent in the basic formulation of the problem (7) appears to be potentially promising from the point of view of simplifying the solution process.We will call this version of problem (7) common, because in this interpretation problem (7) becomes an ordinary problem of linear programming for the solution of which there is a universal simplex method.Both the reduced and the usual variation of the formulation of the optimization problem (7) make it possible to find the desired estimates of the term r 2 (W c ) from the expression (6) as a result of the solution.
Let's finish the analysis of Case 2 by considering the variant when the derivative set W c obtained from the original set by applying the partitioning rule β qg = 0. Obviously, for this variant, the estimate (6) will be equal to the estimate determined for the original set.
Case 3.This case can be implemented only for the set W c for which there is a single variant of the placement plan A. The non-coverage function h(A) for such a placement plan is defined as the sum of the values of all partial non-coverage functions (4) calculated for the points of the segment ρ in which the base stations are located (including gateways {b 0 , b n+1 }).
If the requirements 1-4 defined at the end of Section 2.1 are met for the found placement plan A and inequality is satisfied, then the value of the function h(A) becomes extreme:h Â À � ≔h A ð Þ, and the placement plan A becomes a new optimal one: Â≔A.Next, the reverse transition is implemented according to the rule for the formation of derived subsets presented at the beginning of Section 2.2.If the inequality ( 10) is not satisfied, then the values of h Â À � , Â are not updated, but the reverse transition is simply implemented according to the rule for the formation of derived subsets presented at the beginning of Section 2.2.
The operation of the branch-and-bound method ends when all vertices of the search tree are closed.The result of the method is a tuple consisting of

Determination of the optimal architecture and quality characteristics of the elements of the linearly oriented segment of the 5G network
Let us recall the formulation of the extremal problem (1).It is necessary to find such an admissible placement plan A* of base stations from the set B, for which h(A*) = min{h(A), A 2 W}.
Let us determine for this problem the sequence of admissible placement plans O = A 1 , A 2 , . .., A k of the set W for the given k.We will form the sequence O so that each placement plan in its composition is no better than the previous one and no worse than the next one: Starting with A 1 , we will sequentially evaluate each subsequent placement plan from the sequence O.As soon as we find a placement plan that satisfies all the requirements formulated in Section 2.1, this plan will represent the optimal architecture of the target linear segment of the 5G network.Indeed, for all previous plans, these requirements are not fulfilled, and all subsequent placement plans after the one found in the sequence will not exceed the last one in terms of the value of the criterion h(A).
Let us pay attention to the issue of forming the sequence O using the interpretation of the branch-and-bound method described in Section 2.3.Let us reduce inequality (10) . By sequentially arranging all the optimal placement plans found during the method's operation, we will obtain a sequence in which each placement plan will be no worse than the previous one and no better than the next one.To obtain the desired sequence, it is enough to write down the sequence found as a result of the branch-and-bounds method in reverse order.
The advantage of the proposed approach to obtaining the sequence O is simplicity.The disadvantage is that only placement plans no worse than A 1 will be included in the resulting sequence.The found sequence O may not contain a placement plan satisfying all optimality criteria.To eliminate this shortcoming, we can resort to the expansion of the set O. Let's say that as a result of solving the extremal problem (1), we want to get not only the optimal solution but also all solutions that are worse than the optimal one by a value no greater than Δ = lε > 0: where ε is the percentage deviation.In practice, the described procedure is implemented by replacing constraint (9) with inequality (11) in the optimization problem (7).By varying the value of Δ, we can get a close to optimal sequence of O, devoid of the above-mentioned drawback.Now let's formalize the calculation of such qualitative characteristics of the obtained optimal architecture of the linearly oriented segment of the 5G network as connectivity semidiameters S, coverage semidiameters s and transfer delay between network endpoints (gateways) D p2p .
To determine the first two indicators, the primary thing is to establish the energy potential of the network architecture.We characterize this indicator by the equality where P tran is the transmitter power; L tran , L rec are signal losses on the transmitter and receiver sides, respectively; G tran , G rec are the main characteristics of the transmitter and the receiver, respectively; L link is a signal loss in the communication channel (air); S rec is the sensitivity of the receiver; SFM is the Signal Fading Margin.The value of the indicator L link is determined using the Friis formula [33]: where F 5G is the central frequency of the 5G communication system; S t−r is the distance between the transmitter and the receiver; K(F, S) is a constant, the value of which depends on the units of measurement of the characteristics F and S (for 5G: Generalizing the expressions (12), and (13), we will express the indicator S t−r .After simplification, we get: Expression ( 14) is the basis for calculating semidiameters s i , S ij for base stations from the set B. Now let's formulate an expression for calculating the transfer delay between the endpoints (gateways) of the network.An acceptable level of abstraction will be the simulation of the investigated segment of the 5G network by a queuing system with cross-traffic and M/M/ 1-type nodes.According to Burke's theorem [34], the flow at the exit of the node, which means at the entrance of each subsequent phase, is Poisson.The intensity at the output of each phase is equal to the sum of the intensities of all incoming flows (as noted in Section 2.1, such an intensity is equal to μ).We assume that the real bandwidth of a typical base station from the set B is half of that declared in the specification.
Taking into account the described initial postulates, we characterize the service intensity for the b i -th base station by the expression η i = 0.5t i /r, where t i is the throughput of the b i -th base station and r is the average length of the data packet.For 8b i 2 B, the load factor is defined as where n is is the number of incoming flows.Condition t i < 1 is a necessary and sufficient condition for the existence of a stationary mode of operation of the described queuing system.With the specified parameters t i , according to Little's law, the average transfer delay for the b i -th base station of the linearly oriented segment of the 5G network is equal to Therefore, the transfer delay between endpoints (gateways) of the studied network is described by the expression

Results
We will demonstrate the functionality of the mathematical apparatus proposed in Section 2 on an example.Let a linear area ρ with a length of l = 2.  12)), and the superscript {S, s} means the orientation of the corresponding parameter, namely: • the index S means directing the power to ensure the connectivity of base stations in the studied network segment (maintenance of communication between base stations), • the index s means directing the power to cover the area around the base stations from the set B (supporting communication between users); • parameter L t/r = L train = L rec characterizes signal loss; parameter t characterizes a throughput of the base station; parameter v characterizes the cost of placing base station b i 2 B at point FThe maximum amount of funds that can be spent on the deployment of the network architecture B is represented by the value of the parameter V = 12 [c.u.]The maximum permissible delay in the transmission of a data packet between network users is represented by the value of the parameter D = 1.5 [ms].We assume that the intensity of data packet arrival is μ = 100 [1/s], and the average length of a data packet is r = 1500 [bytes].
According to expression (14) (taking into account expressions (8), ( 12), ( 13)), we calculate the semidiameters s i , S ij , i ¼ 0; 6, j ¼ 0; 6, i 6 ¼ j, for all base stations that should be placed on the area ρ (5 base stations from the set B: b 1 , b 2 , . ..,b 5 , and two gateways b 0 , b 6 ).When calculating semidiameters s i , the allowance for signal fading is SFM = 14, and when calculating semidiameters S ij , we took SFM = 20.For both series of calculations, the centre frequency is F 5G = 3.4 [GHz].The results of the calculations are visualized in the form of diagrams in Figs 3 and 4.
Now we know all the initial data for the application of the adapted branch-and-bound method presented in Section 2. • the value of the non-coverage function h Â À � , • the cost of implementation of the placement plan V Â À � , • the delay D Â À � , • the number of the vertex of the search tree N Â À � , on the best placement plan Â found on the N i -th iteration of the implementation of the adapted branch-and-bounds method.
From Fig 5, we see that the optimal value of the objective function h(A*) = 0 is reached at the N i = 9-th iteration in the 213-th vertex of the search tree, subject to the constraints of V (A*) � V = 12, D(A*) � D = 1.5.The optimal value of the objective function h(A*) = 0 found on the N i = 9-th iteration means that if the placement plan A* is implemented in the area ρ, there will be no fragments not covered by 5G communication.
Corresponding to the one presented in

Discussion
The results of solving the problem of optimal planning of a linearly oriented segment of the urban 5G network with the main output parameters l = 2.3, |P| = 5, |B| = 6 presented in the previous section testify to the functionality of the mathematical apparatus proposed in Section 2. Let's pay attention to those shown in Fig 5 graphs.The linear nature of the dependencies D = f(N i ) and V = f(N i ) indicates that the work of the method began with a near-optimal superposition of the values of the controlled parameters.The close-to-linear nature of the dependences of h = f(N i ) and N = f(N i ) is a sign of the computational efficiency of the adaptation of the We will investigate the dependence of such qualitative indicators as the solution time T [s], non-coverage h [m] and solution accuracy Acc [%], for the three options for calculating the term r 2 (W c ) from expression (6) described in Section 2.3, which determines the evaluation of the not-covering function h(A) when solving the extremal problem (1) by the adapted branchand-bounds method.For the definition r 2 (W c ), we proposed three variants of the formulation of the optimization problem (7): • base (BOP, optimization problem (7) with constraints (8), ( 9): integer Linear Programming Problem), • reduced (ROP, optimization problem (7) with constraint (8): integer Knapsack Problem), • common (COP, optimization problem (7) with constraints (8), ( 9): ordinary Linear Programming Problem).From those presented in Figs 6-8 results, it can be concluded that to estimate the term r 2 (W c ) from the expression (6) with a large dimension of the optimization problem (1), the statement COP should be used.At the same time, we will play slightly lower in accuracy (see Fig 3), but we will significantly (more than 6 times for large pairs (|P|, |B μ |)) win in the speed of obtaining a solution.

Conclusions
The need for optimal planning of linearly oriented segments of urban 5G networks is urgent because the landscape of a modern metropolis is organized, for the most part, in a straight line.Separately, we note that 5G base stations are quite compact and do not require separate towers to be erected for installation.Even outside the city, straight lines (highways, railways) constantly accompany a person.Note that such an optimization task can be classified as multi-criteria in the metric of service quality indicators for wireless users.
As is well known, the disadvantage of multi-criteria optimization problems is the probable lack of an analytical solution, and most solution methods are not focused on evaluating suboptimal results, which are approximate, but partially do not correspond to the imposed initial constraints.So, in the article, the extreme problem of finding the optimal placement plan of 5G base stations at certain points within a linear section of finite length is set.A fundamental feature of the author's formulation of the extreme problem is that it takes into account not only the points of potential placement of base stations but also the possibility of selecting instances of stations to be placed at a specific point from a defined excess set, as well as the aspect of inseparable interaction of placed 5G base stations within the framework of SON.This approach is more realistic, however, it is much more difficult to obtain numerical results.This approach is more realistic, however, it is much more difficult to obtain numerical results.
To overcome these difficulties, the article proposes a method that allows selecting the best options for the placement of base stations in the process of synthesizing the architecture of a linearly oriented segment of a 5G network for further evaluation of the received plans in the metric of defined performance indicators.The mentioned procedure for obtaining a sequence of the best placements is formalized in the form of an adapted branch-and-bounds method.As the final stage of the synthesis of the optimal plan of a linearly oriented 5G network segment based on the sequence of the best placements, it is proposed to expand the parametric space of the design task due to the specific technical parameters characteristic of the 5G platform.
The article presents a numerical example of the synthesis of the optimal placement plan of 5G base stations at a set of admissible points within the target linearly oriented aria of a given length.It is shown that the presented mathematical apparatus allows for the formation of a set of optimal placements taking into account the length of the non-coverage of the target area.To calculate this characteristic parameter, both exact and two approximate approaches are formalized.The results of the experiment showed that for high-dimensional problems, the approximate approach allows for a reduction in the computational complexity of implementing the adapted branch-and-bounds method by more than six times, with a slight loss of accuracy of the optimal solution.
The authors plan to direct further research to the expansion of the parametric space of the formulated extreme problem with the nomenclature of characteristic parameters inherent in the 5G platform.This will complicate the solution process, however, it will allow us to take into account the functioning of such relevant 5G technologies as OFDM, Network Slicing, etc.
n o is the coverage semidiameter of the service supported by the base station b i , focused on the connection of the i-th base station with the j-th base station, b i , b j 2 B; t i and v i are the bandwidth and cost of the base station b i , respectively.

Fig 2 .
Fig 2. Moving through the search tree.https://doi.org/10.1371/journal.pone.0299000.g002 3 for finding the optimal placement plan A* of base stations from the set B on the aria ρ.Recall that the objective function (1) is oriented to the minimization of the non-covering function h Â À � .Fig 5 presents graphs of the dependence of indicators such as: Fig 5 of the iterative process, the dynamics of the evolution of placement plans Â !A * are visualized in Fig 6.The optimal placement plan A