=Paper=
{{Paper
|id=Vol-1623/paperco17
|storemode=property
|title=About Local Optimum of the Weber Problem on Line with Forbidden Gaps
|pdfUrl=https://ceur-ws.org/Vol-1623/paperco17.pdf
|volume=Vol-1623
|authors=Gennady Zabudsky, Natalia Veremchuk
|dblpUrl=https://dblp.org/rec/conf/door/ZabudskyV16
}}
==About Local Optimum of the Weber Problem on Line with Forbidden Gaps==
About Local Optimum of the Weber Problem on Line with Forbidden Gaps Gennady Zabudsky⋆ and Natalia Veremchuk Sobolev Institute of Mathematics Siberian Branch of RAS, Omsk, Russia zabudsky@ofim.oscsbras.ru,n-veremchuk@rambler.ru Abstract. The location problem of interconnected facilities on the line with forbidden gaps is considered. The located facilities are connected among them- selves and with gaps. Location in forbidden gaps is not allowed. It is need to minimize the total cost of connections between facilities and between facilities and gaps. It is known that the initial continuous problem is reduced to series discrete subproblems of smaller dimension. In this paper the definition of the local optimum of the problem is introduced. In order that to obtain the local op- timum it is necessary to solve some subproblems. The variants of lower bounds of the goal function of the subproblems are proposed. The bounds can be used in the branch and bounds algorithm for solving the subproblems. Keywords: forbidden gaps, local optimum, lower bounds, Weber problem 1 Introduction The analysis and decision of location problems are intensively developing directions in operations research [2, 4, 5, 17, 18]. Problems of such class have important applications and arise in various fields of activity: at designing of master plans of the enterprisers, arrangement of processing equipment in shops, definition of the locations of the service stations, etc. The various statements of such problems are defined by size of facilities, area in which they are located, various restrictions and types of criteria and so on. An important subclass of the location problems of the interconnected facilities is the Weber problem with criterion of minimum of total cost of connections [8, 12, 17, 18]. The Weber problem is an adequate model of many practical applications. For example, in automation of design of the complex systems can be the problem of placement of structural elements in the area with implementation of certain requirements [7, 16]. At designing of the plan of the petrochemical enterprise the facilities are the equipment [7, 16]. The equipment can be connected among themselves by various communications, for example, systems of the pipelines. Regularity of location often is required to ensure the straightforward paths and convenient maintenance along so-called “red” lines [15]. ⋆ This work was supported by grant 16–01–00740 from the Russian Foundation for Basic Research. Copyright ⃝ c by the paper’s authors. Copying permitted for private and academic purposes. In: A. Kononov et al. (eds.): DOOR 2016, Vladivostok, Russia, published at http://ceur-ws.org 116 G. Zabudsky and N. Veremchuk The Weber problem with restrictions on facility location and/or traveling has been considered widely in recent years [4, 6]. Depending on the type of restrictions, such problems are divided into three categories. The category considers planar facility lo- cation problems in which restrictions come from existence of forbidden regions (gaps). Forbidden gaps refer to prohibition of facility placement but allowance of free travel- ing. Moreover, some forbidden gaps may contain the equipment, for example, during modernization of the enterprise. In the second category, restrictions are imposed by barriers. Barriers are defined as regions where neither locating a facility on not passing through is allowed. The last category considers restricted planar location problems with congested regions. Congested regions are bounded areas in the plane that forbid facility location however passing through their interior is possible at some extra traveling cost. A broad overview on facility location problems with forbidden gaps is provided in [6]. If the facilities are commensurable to the area of location then they often are ap- proximated by the rectangles; otherwise, we may consider them as material points. Various approaches to the problem of optimal location of the rectangles are developed [1, 8–10, 13, 15]. In case the rectangles are unconnected, the two–dimensional problem of packing of the rectangles into the strip of the minimal length is considered in [9]. The problem is formulated as a mixed integer nonlinear programming problem. For solving this problem the probabilistic tabu search algorithm is developed. In [15] the location problem of the rectangles on the parallel lines is considered. For constructing a set of Pareto–optimal solutions, integer optimization and dynamic programming are applied. In [8, 10], the algorithms of local optimization on the plane and dynamic programming on the line without forbidden gaps are described. For the problem in which the rect- angles can be rotated and it is possible to create routes for establishing connections, the heuristic algorithm is proposed [1]. An important subclass of problems of optimal location of rectangles is the VLSI chip design problems. The wide review of practical applications and various approaches for solving such problems is provided in [13]. In this paper, we study the Weber problem on the line with forbidden gaps. The located facilities are segments, for example, projection of the rectangles onto the line. It is known that the initial continuous problem is reduced to series discrete subproblems of smaller dimension. All subproblems have identical structure. The local optimum of the problem is obtained by solving some subproblems. The variants of lower bounds of the goal function of the subproblems are proposed. The bounds can be used in algorithms for solving of the subproblems, for example, in the branch and bounds algorithm. 2 Problem Formulation and Basic Properties The Weber problem on the line with forbidden gaps is formulated as follows. Consider a straight–line segment of length LS containing some fixed rectilinear areas (forbidden gaps). It is necessary to locate on this line facilities which are the segments, whose centers are connected with each other and with the forbidden gaps. Location in forbid- den gaps is not allowed. The problem is to locate facilities on the segment of length LS outside the forbidden gaps so that they do not intersect each other and, moreover, the total cost of connections between the facilities and between facilities and gaps is About Local Optimum of the Weber Problem 117 minimized [17]. Without loss of generality, we can assume that the left border of the segment of length LS is the origin. Let Xi and Fj denote the facilities and gaps with the centers at xi and bj and lengths of li and pj , where i ∈ I = {1, . . . , n} and j ∈ J = {1, . . . , m}; while wij ≥ 0, uik ≥ 0 are the specific costs of connections between Xi and Fj , Xi and Xk for i, k ∈ I, j ∈ J, and i < k. The target is to locate the facilities X1 , . . . , Xn on the segment outside gaps F1 , . . . , Fm and so that they do not intersect with each other and the total cost of the connections between the facilities and between facilities and gaps is minimized. The mathematical model under consideration is as follows: ∑ n ∑ m ∑ n−1 ∑ n G(x) = wij |xi − bj | + uik |xi − xk | → min, (1) i=1 j=1 i=1 k=i+1 li + pj |xi − bj | ≥ , i ∈ I, j ∈ J, (2) 2 li + lk |xi − xk | ≥ , i, k ∈ I, i < k, (3) 2 li li ≤ xi ≤ LS − , i ∈ I. (4) 2 2 The first component in (1) defines the total cost of connections between facilities and gaps; the second part defines the total cost of connections between the facilities them- selves; and (2) and (3) are the conditions of disjointness. The feasible area B is disconnected and consists of the set of r disjoint ∪ segments (blocks) Bk of length Lk that contain the facilities Xi , i ∈ I, B = k=1,r Bk . The problem (1)–(4) is NP-hard; a feasible solution to the problem can be found by con- struction of a one–dimensional packing into containers [3]. In this case, the facilities with lengths of li , i ∈ I, are packed into containers with sizes Lk , k = 1, r. Moreover, if there are no gaps then (1)–(4) is an optimal linear ordering problem [11] which is NP-hard for an arbitrary set of connections between the facilities. In [17] the algorithm for obtaining an approximate decision in two phases is of- fered. In the first phase, we find an feasible partition of the facilities into blocks, and in the second phase, the facilities in the blocks are interchanged for the purpose of mini- mization of total cost of connections. A computational experiment with this algorithm and solving of the problem using IBM ILOG CPLEX and a mixed-integer linear pro- gramming model is considered. This work is continuation of researches of the problem (1)–(4). Given a feasible location, a remainder in the block Bk is a segment between two adjacent facilities that do not have a common border or between the border of Bk and an adjacent block. A block without facilities is called a block with remainder. Two elements (facilities, gaps, remainders) are called glued if they have a common border. Let JL (Bk ) and JR (Bk ) denote the set of gaps to the left and to the right of the block Bk , let IL (A) and IR (A) denote the set of facilities to the left and to the right of 118 G. Zabudsky and N. Veremchuk the block (gap) A correspondingly. Given an facility Xi in Bk , let Lwi and Rwi denote the total cost of connections for Xi defined as ∑ ∑ ∑ ∑ Lwi = wij + uik , Rwi = wij + uik . j∈JL (Bk ) k∈IL (Bk ) j∈JR (Bk ) k∈IR (Bk ) Let x = (x1 , . . . , xn ) be a feasible solution of (1)–(4) that uniquely defines the partition of X1∪, . . . , Xn into blocks. Let Ik (x) denote the set of the facility’s numbers in Bk , I = k=1,r Ik (x), and let Hk (x) denote the set of remainders in Bk for x. Let nk denote the size of set Ik (x), then |Hk (x)| ≤ nk + 1. Note that x can be represented as x = (x1 , . . . , xr ), where xk is the coordinates of the centers of the facilities that are located in Bk and have numbers from Ik (x). Proposition 1. (see [17]) Given a feasible solution x of (1)–(4), we can find a feasible solution x′ such that |Hk (x′ )| ≤ 1, k = 1, . . . , r, G(x′ ) ≤ G(x). Thus in each block Bk it is possible to consider no more than one remainder, length of which we denote by ∆k . Let LBk and RBk denote the coordinates of the left and right borders of Bk . Then, for a fixed partition of facilities into blocks, the goal function G(x) can be presented as ∑ r G(x) = Gk (xk ) + C, (5) k=1 where ∑ ∑ ∑ ( ∑ Gk (xk ) = ust |xs − xt | + |xs − LBk | wsj + s∈Ik (x) t∈Ik (x),t>s s∈Ik (x) j∈JL (Bk ) ∑ ) ∑ ( ∑ ∑ ) + usi + |xt − RBk | wtj + uti , i∈IL (Bk ) t∈Ik (x) j∈JR (Bk ) i∈IR (Bk ) C is some constant. The first component of Gk (xk ) is the sum of costs of connections between the facilities in Bk , the second is the total cost of connections between the facilities from Bk and LBk , and the third component is the total cost of connections between the facilities from Bk and RBk . In the case when the area of location on the segment is bounded on the left and on the right by gaps, we have ∑ ( ∑ 1 ∑ ∑ m−1 1 C = p1 wk1 + pj wkj + wkj + 2 2 j=2 k∈IR (F1 ) k∈IL (Fj ) k∈IR (Fj ) ∑ ∑ ∑ ∑ ∑ ∑ ) +2 wkt + 2 wks + 2 ukt + k∈IL (Fj ) t∈J,t>j k∈IR (Fj ) s∈J,ss s∈Ik (x) ∑ + wtR |xt − RBk | → min, (6) t∈Ik (x) li + lk |xi − xk | ≥ , i, k ∈ Ik (x), i < k, (7) 2 li li LBk + ≤ xi ≤ RBk − , i ∈ Ik (x). (8) 2 2 It is need to find coordinates xk of the centers facilities in Bk , so that the total cost of connections between located facilities among themselves and with FL and FR is minimized. ∩ As Ik (x) Il (x) = ∅ for all k, l = 1, . . . , r, then to find a local optimum for some fixed partition of the facilities into blocks, it is sufficiently to find the minimum in r independent subproblems (6)–(8). These subproblems have identical structure. Proposition 2. For search of the local minimum of the problem (1)–(4) it is suffi- ciently to solve r subproblems (6)–(8). The proof of proposition 2 follows from representation of the goal function (5). Note that the process of search of the local minimum can be parallelized as solving r independent subproblems (6)–(8). The properties stated above allowed to reduce the initial continuous problem to series of discrete subproblems of smaller dimension. For finding of the local optimum of the problem it is need to solve r subproblems. 120 G. Zabudsky and N. Veremchuk Solving of the problem (1)–(4) is sequential performance of two phases. In the first phase, we find feasible partition of the facilities into blocks by means of the algorithm [17]. In the second phase, for obtained partition the series of the subproblems of smaller dimension are solved. Thus we fined the local optimum of the problem (1)–(4). The local optimum with minimal value of the goal function is the exact decision of the problem (1)–(4). For the search of approximate decision of the problem (1)–(4) the stopping criteria of the algorithm may be running time, number of iterations, finding of the exact decision or given estimate of accuracy. Note that the number of possible partitions Xi , i ∈ I into the blocks B1 , . . . , Br is at most rn . The remainder in Bk can be considered as the additional located facility Xnk +1 , for which lnk +1 = ∆k and Lwnk +1 = Rwnk +1 = 0. Note that if ust = 0, ∀s, t ∈ Ik (x), s < t, then in the second phase the exact decision for the subproblem (6)–(8) is found by the polynomial algorithm [17]. In this case the graph, whose tops are the left and the right borders of the block Bk and the facilities in Bk , has a series–parallel type with one top on each chain between the source and the sink. For such structure of the graph the effective algorithm is offered in [14]. Generally, if ∃s, t ∈ Ik (x) : ust > 0, then for solving the subproblem (6)–(8) for small value nk , it is possible to use nk ! permutations of the facilities into block. In general case it is possible to apply, for example, the branch and bounds algorithm. In the algorithm calculation of the lower bounds of the goal function is important. 3 Lower Bounds Let some located facilities in Bk are glued to the left border of Bk , some located facilities are glued to the right border of Bk , and some facilities aren’t located yet. Denote ∪ by N Fl , N Fr the sets of the located facility’s numbers in Bk , then Ik (x)\{N Fl N Fr } is the set of the unlocated facility’s numbers. Without loss of generality, we shall consider that the facilities in the set N Fl have numbers from 1 to s, and in the set N Fr the facility’s numbers are from t + 1 to nk . Denote by D the set of the admissible locations of the facilities in Bk when N Fl and N Fr are defined (partial location). Let ξ(D) be the lower bound of function Gk (xk ) for the set D. Then ξ(D) can be presented as the sum of three components. The first component ξ1 (D) is the total cost of connections between the located facilities themselves and with FL and FR . The second component ξ2 (D) is lower bound of total cost of connections the unlocated facilities with FL , FR and with the located facilities in Bk . The third component ξ3 (D) is lower bound of total cost of connections between the unlocated facilities themselves. Thus ξ(D) can be presented as ξ(D) = ξ1 (D) + ξ2 (D) + ξ3 (D). The coordinates of the centers of the located facilities are known and ξ1 (D) can be calculated as follows ∑ ∑ ∑ ξ1 (D) = ust |xs − xt | + |xs − LBk |· ∪ ∪ ∪ s∈{N Fl N Fr } t∈{N Fl N Fr },t>s s∈{N Fl N Fr } ( ∑ ∑ ) ∑ ( ∑ ∑ ) · wsj + usi + |xt − RBk | wtj + uti . ∪ j∈JL (Bk ) i∈IL (Bk ) t∈{N Fl N Fr } j∈JR (Bk ) i∈IR (Bk ) About Local Optimum of the Weber Problem 121 Further two variants of calculation ξ2 (D) are∪offered. The first variant. For each i ∈ Ik (x)\{N Fl N Fr } the total cost of connections between the located facilities themselves and with FL and FR is determined as follows ∑ ∑ SL(i) = Lwi + uik , SR(i) = Rwi + uik . k∈N Fl k∈N Fr Further location of the unlocated facilities in Bk is defined in two ways. In the first way, the facilities are ordered by not increase of the relations SL(i)/li . The facilities consistently are pasted together in that order with the most left located facility in Bk . Let, for simplicity of designations, the glued unlocated facilities have numbers from s + 1 to t. In the second way, the unlocated facilities are ordered by not increase of the relations SR(i)/li . The facilities consistently are pasted together in that order with the most right located facility in Bk . Let the glued unlocated facilities have numbers from t to s + 1. Then ξ2 (D) = ξ2L (D) + ξ2R (D), where ξ2L (D) and ξ2R (D) are the lower bounds of total cost of the connections unlo- cated facilities with FL , FR respectively and with the located facilities in Bk . In both cases we receive some permutation π of facilities in the block Bk . Consider the permutation π and two distinct facilities Xi and Xj in Bk . The distance between Xi and Xj with respect to this permutation, assumed to be taken their centers, is equal to the half-length of Xi , plus the lengths lk of all facilities which are between Xi and Xj in π, plus the half-length of Xj . The half-length of Xi plus the half-length of Xj is the constant. Thus, at the calculation of the distance between the facilities Xi and Xj it is sufficient to consider only the lengths lk of all facilities which are between Xi and Xj in π. Then ξ2L (D) and ξ2R (D) can be calculated as follows: t ( ∑ ∑ q−1 ∑ s ∑ q−1 ) ξ2L (D) = Lwq lg + uqi lk , q=s+1 g=1 i=1 k=i+1 t ( ∑ ∑ nk ∑ nk ∑ i−1 ) ξ2R (D) = Rwq lg + uqi lk . q=s+1 g=q+1 i=t+1 k=q+1 The proof that ξ2L (D) and ξ2R (D) are the lower bounds of the total cost of connections unlocated facilities with FL , FR and with the located facilities in Bk is similar to the proof in [10]. ∪ The second variant. ∪ ∪The set Ik (x)\{N Fl N Fr } can be presented as union of non– crossing sets NL NC NR , where by NL , NC , NR are denoted the sets of facility’s numbers for which SL(i) > SR(i), SL(i) = SR(i), SL(i) < SR(i) respectively. Further facilities with numbers from NL are ordered by not increase of the relations (SL(i) − SR(i))/li . The facilities consistently are pasted together in that order with the most left located facility in Bk . The facilities with numbers from NR are ordered by not increase of the relations (SR(i) − SL(i))/li , and they consistently are pasted together in that order with the most right located facility. The facilities with numbers 122 G. Zabudsky and N. Veremchuk from NC are located between sets of the facilities with numbers from NL and from NR in any order. ∪ Thus, for∪each i ∈ Ik (x)\{N Fl N Fr } coordinate of the center is determined. Let Ik (x)\{N Fl N Fr } = {s + 1, . . . , t} as well as earlier. Determine the following value Z as t ( ∑ ∑ q−1 ∑ s ∑ q−1 ∑ nk ∑ nk ∑ j−1 ) Z= Lwq lg + uqi lk + Rwq lh + uqj lv . (9) q=s+1 g=1 i=1 k=i+1 h=q+1 j=t+1 v=q+1 Proposition 3. The value Z is the lower bound of the total cost of connections the unlocated facilities with FL , FR and with the located facilities in Bk . Proof. Let NL ̸= ∅ and NR = ∅ and NC = ∅. Without loss of generality, we consider the case when N Fl = ∅ and N Fr = ∅. Otherwise it is possible to redefine the cost of connections the unlocated facilities with FL , FR and with the located facilities. Then the set of the unlocated facility’s numbers is Ik (x) and Lwi > Rwi for each i ∈ Ik (x). Let (Lw1 − Rw1 )/l1 > . . . > (Lwnk − Rwnk )/lnk , then the facilities are located in order X1 , . . . , Xnk according to the second variant of the calculation ξ2 (D). Denote by Π ∗ = (1, . . . , nk ) the permutation of numbers facilities corresponding to order X1 , . . . , Xnk and denote by Z(Π ∗ ) the value Z for Π ∗ . On the formula (9) we receive Z(Π ∗ ) = Lw2 l1 + . . . + Lwt (l1 + . . . + lt−1 ) + Lwt+1 (l1 + . . . + lt ) + + . . . + Lwnk (l1 + . . . + lnk −1 ) + Rw1 (l2 + . . . + lnk ) + . . . + +Rwt (lt+1 + . . . + lnk ) + Rwt+1 (lt+2 + . . . + lnk ) + . . . + Rwnk −1 lnk . To prove that Z(Π ∗ ) is the lower bound for Z(Π), where Π is any permutation of numbers facilities from Ik (x). Assume, that Z(Π ∗ ) is not minimum. Then there is other permutation Π, such that Z(Π ∗ ) > Z(Π). Any permutation can be received by series of transposition of adjacent numbers, then it is sufficiently to consider the proof for change the position of two adjacent facilities. Denote by Π = (1, . . . t − 1, t + 1, t, t + 2, . . . , nk ), then on the formula (9) we receive Z(Π) = Lw2 l1 + . . . + Lwt+1 (l1 + . . . + lt−1 ) + Lwt (l1 + . . . + lt−1 + lt+1 ) + + . . . + Lwnk (l1 + . . . + lnk −1 ) + Rw1 (l2 + . . . + lnk ) + . . . + +Rwt+1 (lt + lt+2 + . . . + lnk ) + Rwt (lt+2 + . . . + lnk ) + . . . + Rwnk −1 lnk . By assumption Z(Π ∗ ) − Z(Π) > 0, then Z(Π ∗ ) − Z(Π) = Lw2 l1 + . . . + Lwt (l1 + . . . + lt−1 ) + Lwt+1 (l1 + . . . + lt ) + + . . . + Lwnk (l1 + . . . + lnk −1 ) + Rw1 (l2 + . . . + lnk ) + . . . + +Rwt (lt+1 + . . . + lnk ) + Rwt+1 (lt+2 + . . . + lnk ) + . . . + Rwnk −1 lnk − −Lw2 l1 − . . . − Lwt+1 (l1 + . . . + lt−1 ) − Lwt (l1 + . . . + lt−1 + lt+1 ) − . . . − −Lwnk (l1 + . . . + lnk −1 ) − Rw1 (l2 + . . . + lnk ) − . . . − About Local Optimum of the Weber Problem 123 −Rwt+1 (lt + lt+2 + . . . + lnk ) − Rwt (lt+2 + . . . + lnk ) − . . . − Rwnk −1 lnk = = (Lwt+1 − Rwt+1 )lt − (Lwt − Rwt )lt+1 > 0. Hence Lwt+1 − Rwt+1 Lwt − Rwt > . lt+1 lt It is the contradiction with inequalities (Lw1 − Rw1 )/l1 > . . . > (Lwt − Rwt )/lt > (Lwt+1 − Rwt+1 )/lt+1 > . . . > (Lwnk − Rwnk )/lnk . The proof for the subset NR ̸= ∅ is similarly. Note that facilities with numbers from NC can be located in arbitrary order. Proposition is proved. Therefore, in the capacity of ξ2 (D) it is possible to use the value Z. Following way may ∪ be use for the calculation of ξ3 (D). Let for the facilities Xs , Xt , Xq , s, t, q ∈ Ik (x)\{N Fl N Fr } take place inequalities ust > 0, usq > 0, utq > 0. Consid- ering any three of the interconnected facilities, value of ξ3 (D) can be calculated, for example, as follows ∑ ξ3 (D) = min{ust lq ; usq lt ; utq ls }. ∪ s,t,q∈Ik (x)\{N Fl N Fr } 4 Conclusion The NP–hard location problem of interconnected facilities on the line with forbidden gaps is considered. It is need to minimize the total cost of connections between the facilities and between the facilities and gaps. It is known that the initial continuous problem is reduced to series discrete subproblems of smaller dimension. For finding of the local optimum of the problem it is need to solve some subproblems. Two variants of lower bounds of the goal function of the subproblems are proposed. The bounds can be used in the branch and bounds algorithm for solving the subproblems. References 1. Erzin, A., I., Cho, J., D.: Concurrent Placement and Routing in the Design of Integrated Circuits. Avtomat. i Telemekh. No. 12, 177–190 (2003) [Automat. Remote Control 64(12), 1988–1999 (2003)] 2. Farahani, R., Z., Hekmatfar, M.: Facility location: Concepts, models, algorithms and case studies. Heidelberg: Physica-Verlag, (2009) 3. Garey, M., R., Johnson ,D., S.: Computers and Intractability: A Guide to the Theory of NP–Completeness. (Freeman, San Francisco, 1979; Mir, Moscow, 1982) 4. Klamroth, K.: Single-Facility Location Problems with Barriers. Springer Series in Oper- ations Research, (2002) 5. Kochetov, Yu., A., Panin, A., A., Plyasunov, A., V.: Comparison of metaheuristics for the bilevel facility location and mill pricing problem. Diskret. Anal. Issled. Oper. 22(3), 36-54 (2015) DOI: 10.17377/daio.2015.22.465 6. Nickel, S., Puerto, J.: Location theory. A unified approach. Berlin: Springer, (2005) 124 G. Zabudsky and N. Veremchuk 7. Legkih, S., A., Nagornay, Z., E., Zabudsky, G., G.: Automation of design of plans of production sites and shops of the sewing enterprises (in Russian). Nat. Tech. Scien. No. 4, pp. 261–266 (2005) 8. Panyukov, A., V.: The problem of locating rectangular plants with minimal cost for the connecting network. Diskret. Anal. Issled. Oper. Ser. 2, 8(1), 70–87 (2001) 9. Rudnev, A., S.: Probabilistic tabu search algorithm for the packing circles and rectangles into the strip. Diskret. Anal. Issled. Oper. 16(4), 61–86 (2009) 10. Simmons, D., M.: One-dimensional space allocation: an ordering algorithm. Oper. Res. 17(5), 812–826 (1969) 11. Suresh, G., Sahu, S.: Multiobjective Facility Layout Using Simulated Annealing. Int. J. Prod. Econom. 32(2), 239–254 (1993) 12. Trubin, V., A.: Effective algorithm for Weber problem with a rectangular metrics. Kiber- net. No. 6, 67–70 (1978) [Cybernetics 14(6), 874–878 (1978)] 13. Wai–Kai, Chen: The VLSI Handbook. CRC Press. (2000) 14. Zabudsky, G., G.: On the problem of the linear ordering of vertices of parallel-sequential graphs. Diskret. Anal. Issled. Oper. 7(1), 61-64 (2000) 15. Zabudskii, G., G., Amzin, I., V.: Algorithms of compact location for technological equip- ment on parallel lines (in Russian). Sib. Zh. Ind. Mat. 16(3), 86–94 (2013) 16. Zabudsky, G., G., Legkih, S., A.: Mathematical model of optimization of flexible modules of processing equipment (in Russian). Appl. Math. Inf. Technol. Omsk: Publishing House of OMGTU, pp.20–28 (2005) 17. Zabudskii, G., G., Veremchuk, N., S.: An Algorithm for Finding an Approxi- mate Solution to the Weber Problem on a Line with Forbidden Gaps. Diskret. Anal. Issled. Oper. 23(1), 82-96 (2016) [J. Appl. Ind. Math. 10(1), 136–144 (2016) DOI: 10.1134/S1990478916010154] 18. Zabudsky, G., G., Veremchuk, N., S.: Solving Weber problem on plane with minimax criterion and forbidden gaps (in Russian). IIGU Ser. Matematika, Vol. 9, pp. 10–25 (2014)