=Paper=
{{Paper
|id=Vol-1623/paperco16
|storemode=property
|title=Separation Problem for k-parashuties
|pdfUrl=https://ceur-ws.org/Vol-1623/paperco16.pdf
|volume=Vol-1623
|authors=Inna Urazova, Ruslan Simanchev
|dblpUrl=https://dblp.org/rec/conf/door/UrazovaS16
}}
==Separation Problem for k-parashuties==
Separation Problem for k-parashuties Inna Urazova1 and Ruslan Simanchev1,2 1 Omsk State University, 55a Mira avenue, 644077 Omsk, Russia urazovainn@mail.ru 2 Omsk Scientific Center of SB RAS, 15 Marksa avenue, 644024 Omsk, Russia osiman@rambler.ru Abstract. This article continues the work [16] in which polyhedral setting of graph approximation problem is provided, support inequalities to polytope are built. In this work NP-hard of the separation problem of k-parachutes relative to M-graph polytope is proved. Keywords: polytope, facet inequality, separation problem. 1 Introduction Let Kn = (V, E) be a complete unoriented graph without loops and multipe edges. Spanning subgraph H ⊂ Kn is called M-graph if each of its connected components is a clique or single-vertex graph. We denote a set of all M-graphs in Kn through µ(V ). Let G ⊂ Kn be some a priori set spanning subgraph. Approximation problem of G consists in finding M-graph H minimizing the functional ρG (H) = |EG ∪ EH| − |EG ∩ EH| (1) on set µ(V ). Here EG and EH are sets of edges of G and H, respectively. In other words, it is necessary to find such a set of pairwise non-overlapping cliques at V which is as far as possible less (in terms of edges) different from G. Harary was the first to state graph approximation problem in [10] in 1955. In the 1960s-1970s in a number of works non-trivial classes of graphs on which the problem is polynomially solvable [7, 18] . In 1986 Krivanek and Moravek considered graph ap- proximation problem as a particular case of the tree clusterization and proved that it was NP-hard [11] . Systematic studies of graph approximation problem began in the last decade when the problem was re-discovered under different names (correla- tion clustering; Cluster editing) by different groups of authors [3, 4, 13]. In particular, the NP-hardness of its various options was established [1, 3, 13] and suggested first ap- proximate algorithms with the guaranteed accuracy evaluation [3, 5, 8]. The best of the currently known graph approximation algorithms finds a guaranteed solution by factor of max. 2.5 worse than the optimum one [2, 17]. 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 110 I. Urazova, R. Simanchev In this work which continues the work [16] we consider polyheldral setting of the graph approximation problem. The polyhedral approach to the solution of extreme com- binatorial problems consists in the correlation of the problem with a special polytope set as a convex hull of incidence vectors of the admissible solutions and, consequently, the use of convex analysis and integer programming means. In this route a special role belongs to the task of polytope description as a set of solutions to the system of linear equations and inequalities. If a full linear description of the polytope is available, the extreme combinatorial problem reduces to the linear programming problem (possibly - with exponential number of constraints) which quite often enables obtaining effective algorithms for the solution thereof. One of the most widely known examples of such a situation is problem of maximum weighted matching [6, 12]. As a rule, for NP-hard problems a full liner description of the polytope is not known (see, e.g., [9]). Never- theless, the availability of a partial description, i.e., classes of valid, support or facet inequalities enables building evaluations of the target function optimum value, develop special cutting plane procedures. To state the results of this article, we will introduce the following designations and notions. For any graph D ⊂ Kn via V D and ED we will denote a set of its vertices and edges, respectively. For an edge e ∈ E we will also use uv notation, where u and v are vertices from V incident to edge e. For D ⊆ Kn and u ∈ V via δD (u) we will denote the set of edges of the D incident to vertex u. If D = Kn , then in this notation we will omit index D. Each set of edges R ⊂ E induces in Kn some subgraph T in which ET = R and V T is a set of vertices incident to edges from R. If no ambiguity arises, a graph induced by a set of edges R will be denoted via R. For subgraphs D, F ⊆ Kn we assume D ∪ F = ED ∪ EF, D ∩ F = ED ∩ EF. We will correlate graph Kn with a Euclidean space RE of dimensions n 2−n by 2 correlating each edge with a coordinate axis in RE . This space may be considered as a space of column vectors with the components∑ indexed by elements from E. If x ∈ RE and R ⊆ E, then we denote linear form xe via x(R). The incidence vector e∈R of the arbitrary graph D ⊆ Kn without isolated vertices is called vector xD ∈ RE with components xD e = 1 at e ∈ ED and xe = 0 at e ∈ D / ED. The latter rule, obviously, sets a one-to-one correspondence between the set of all subgraphs without isolated vertices of Kn and the set of vertices of the unit cube in RE . A set P ⊂ RE is called a polytope if it is a convex hull of a finite number of points. Linear inequality aT x ≤ a0 (a, x ∈ RE , a ̸= 0, a0 ∈ R, ”T ” transposition sign) is called support to P if it is fulfilled for any point from P and there exists at least one point from P converting it into an equality. Any inequality support to P generates the set {x ∈ P |aT x = a0 } which is called the face of the polytope P . Maximum by the inclusion faces of the polytope are called facets. It is clear that a facet is the face and only the face the dimensionality of which is by 1 smaller than the dimensionality of the polytope itself. A support inequality generating a facet will, consequently, be referred to as a facet inequality. Separation Problem for k-parashuties 111 Facet inequalities are present (with the equivalency accuracy) in any linear descrip- tion of the polytope [14, 12]. Besides, they proved to be good as cutting planes for the solution of high dimensionality problems (see, e.g., [15, 9]). During the development of cutting plane procedures using support inequalities the separation problem goes to the foreground. It consists in the following. What is given is a class of inequalities L support to polytope P and point x̄ ∈ RE . It is required to find an inequality in the class L strictly separating the point x̄ from polytope P , or prove that in L such an inequality does not exist. Examples of extreme combinato- rial problems and inequality classes for which the separation problem is polynomially solvable, for example, travelling salesman problem and subtour elimination inequalities [9], the problem of building a multi-processor schedule and inequality class induced by the paths in the precedence graph [15] etc. are known. In this work NP-hardness of the separation problem of k-parachute class inequalities relative to M-graph polytope is proved. 2 M -graph polytope and k-parachutes We will refer to the set Pn = conv{xH ∈ RE |H ∈ µ(V )} as an M-graph polytope. Here ”conv” means ”convex hull”. In [16] it was prove that (0, 1)-vector x ∈ RE is the vector of incidences of M -graph if and only if it satisfies the system: −xuv + xuw + xvw ≤ 1 xuv − xuw + xvw ≤ 1, (2) xuv + xuw − xvw ≤ 1 where u, v, w ∈ V are all kinds of sets of three pairwise different vertices, xuv ≥ 0 for all uv ∈ E. (3) Besides, in the same work it was demonstrated that the target function (1) in terms of RE may be written as ρG (x) = |EG| + x(E Ḡ) − x(EG). (4) Therefore, graph approximation problem may be stated as a problem of integer linear programming (2)–(4). Let U = {u1 , u2 , . . . , uk } and W = {v1 , v2 , . . . , vp } be nonempty subsets of set V , U ∩ W = ∅, k ≥ 1, p ≥ 2 . We will designate star in Kn with the centre in vertex ui and arms ui vj , j = 1, 2, . . . , p through Ti , i = 1, 2, . . . , k. We will denote the clique ∪ k on the set of vertices W through Kp . Let us assume that T = Ti . Graph T ∪ Kp is i=1 called a k-parachute. We will correlate inequality k2 + k x(ET ) − x(EKp ) ≤ . 2 112 I. Urazova, R. Simanchev with this graph. In [16] it was proved that this inequality induced by k-parachute T ∪Kp is support inequality to polytope Pn if and only if p ≥ k , and facet inequality if and only if k = 1. 3 Separation Problem Let L be a class of support inequalities to polytope Pn induced by k-parachutes. Whereas polytope Pn completely lies in the unit cube of RE , we will state the separa- tion problem stronger than it was done in the Introduction. Problem A. Let x̄ ∈ RE and 0 ≤ x̄ ≤ 1. Does L class contain such an inequality for 2 which x̄(ET ) − x̄(EKp ) > k 2+k ? We need auxiliary fact, which is easily proved by induction. Lemma 1. If any t edges with t < n are withdrawn from the complete n-vertex graph, the remaining graph will contain clique of the order n − t. Theorem 1. The separation problem for the inequalities induced by k-parachutes rel- ative to the polytope of the graph approximation problem is NP-hard. Proof. Let us prove that this problem is NP-hard already at k = 1, i.e., at U = {u}. In this case Problem A may be re-stated as follows: Problem B. Let x̄ ∈ RE and 0 ≤ x̄ ≤ 1. Is there among inequalities of the type ∑p ∑ ∑ p−1 p xuvj − xvi vj ≤ 1 such an inequality which is violated by point x̄? j=1 i=1 j=2,j>i It is clear that in this problem u vertex may be fixed. Now at preset point x̄ ∈ RE , 0 ≤ x̄ ≤ 1 and vertex u ∈ V we will define vector c ∈ RE : { x̄e , e ∈ δ(u); ce = −x̄e , e ∈ E \ δ(u). We will denote this edge weighted graph via Kn (c, u). It is not difficult to see that Problem B is equivalent to the following problem. Problem C. Does this edge weighted graph Kn (c, u) contain such clique K that c(EK) > 1? Let us reduce CLIQUE problem stated as follows to Problem C: what is given is graph G and natural number s < |V G|. Does graph G contain clique of the order larger than s? So, let what is given be graph G, |V G| = m and a natural number s < m. Let us assume that n = m + 1 and build a complete graph Kn on a set of vertices V G ∪ {u} where u is an added vertex. Let us assign the following weights to Kn edges: 1 s , e ∈ δ(u); ce = 0, e ∈ EG; 1 − s , e ∈ E Ḡ, where Ḡ is the complement of graph G. Separation Problem for k-parashuties 113 Let us note that if in graph G there is clique K of the order larger than s, then for the clique K ′ ⊂ Kn (c, u) on the set of vertices V K ∪ {u} we have c(EK ′ ) > 1. Let us suppose that in graph Kn (c, u) there is clique K̄ ⊂ Kn (c, u) such that c(E K̄) > 1. Because only edges from δ(u) have positive weights, then necessarily u ∈ V K̄. It means that clique K̄ may be presented as K̄ = δK̄ (u) ∪ K where K is clique in G ∪ Ḡ. Clique K , in it turn may be represented as K = (K ∩ G) ∪ (K ∩ Ḡ). Therefore, K̄ = δK̄ (u) ∪ (K ∩ G) ∪ (K ∩ Ḡ). Because edge sets of these three graphs do not overlap pairwise, then c(E K̄) = c(δK̄ (u)) + c(E(K ∩ G)) + c(E(K ∩ Ḡ)) = 1 1 1 = |V K| − |E(K ∩ Ḡ)| = (|V K| − |E(K ∩ Ḡ)|) > 1. s s s Hence, particularly, it follows that |V K| > |E(K ∩ Ḡ)|. Let us note that the with- drawal of edges E(K ∩ Ḡ) from clique K gives exactly K ∩ G graph which, first, lies fully in G, and second, by virtue of the lemma, contains clique of the order |V K| − |E(K ∩ Ḡ)| > s. The theorem is proved. 4 Conclusion This article continues the work [16] in which facet inequalities to the graph approxi- mation problem polytope is built. The next step in this direction is the inclusion of the received facet inequalities into polyhedral type algorithms. The effectiveness of this in- clusion depends essentially on the computational difficulties of the separation problem of inequalities. In this work N P -hardness of the separation problem of k-parachutes relative to M -graph polytope is proved. The most natural direction for further research is the development of heuristics to separation problem. In conclusion, we will announce preliminary results of our numerical experiment. The aim of the experiment is the evaluation of the use of 1-parachutes in cutting plane and branch-bound algorithms to solve the graph approximation problem. Below are brief results which are as follows. We solve two integer linear programm. The first problem have function (4) as an objective function and have polyhedron (2)–(3) as a relaxation set. In the second problem, we change the polyhedron, adding to the con- straints (2)–(3) all 1-parachute inequalities. This, of course, imposes serious restrictions on the dimension of the problem. To solve this problem we used IBM ILOG CPLEX Optimization Studio package and Intel(R) Celeron(R) CPU N2830 2.16GHz. Solution time was limited to 2 hours. More than a hundred problems with different objective functions of the form (4) for n ≤ 25 were solved. Table 1 shows the average solving time for the first type and second type problems (columns SR and PR, respectively). As the table shows, the addition of 1-parachutes greatly reduces the time to solve the problem. The authors are grateful to A.V. Kononov for useful advice received when working on an article. 114 I. Urazova, R. Simanchev Table 1. n N umber of variables SR(mean time, sec) PR(mean time, sec) 10 45 0.16 0.23 15 105 0.66 1.99 20 190 46.91 35.35 25 300 6566.36 1464.68 References 1. Ageev A.A., Il’ev V.P., Kononov A.V., Talevnin A.S.: Computational complexity of the graph approximation problem. J. of Applied and Industrial Mathematics. vol. 1, Issue 1, 1–8 (2007) 2. Ailon N., Charikar M., Newman A.: Aggregating inconsistent information: Ranking and clustering. J. ACM. 55, 5, 1–27 (2008) 3. Bansal N., Blum A., Chawla S.: Correlation clustering. Machine Learning. 56, 89–113 (2004) 4. Ben-Dor A., Shamir R., Yakhimi Z.: Clustering gene expression patterns. J. Comput. Biol. 6, 3–4, 281–297 (1999) 5. Charikar M., Guruswami V., Wirth A.: Clustering with qualitative information. J. Comput. Syst. Sci. 71, 3, 360–383 (2005) 6. Edmonds J.: Maximum matching and a polyhedron with 0,1 - vertices. J. of Research National Bureau of Standards. Section B, 69, 125–130 (1965) 7. Fridman G. Sh.: A Problem of Graph Approximation (Russian). Upravlyaemye Sistemy. 8, 73–75 (1971) 8. Giotis I., Guruswami V.: Correlation clustering with a fixed number of clusters. Theory of Computing. 2, 1, 249–266 (2006) 9. Grotschel M., Holland O.: Solution of large-scale symmetric travelling salesman problems. Mathematical Programming. 51, 2, 141–202 (1991) 10. Harary, F.: On the notion of balance of a signed graph. Michigan Mathematical Journal. 2, 143–146 (1955) 11. Křivánek M., Morávek J.: NP-hard problems in hierarchical-tree clustering. Acta infor- matica. 23, 311–323 (1986) 12. Schrijver, A.: Combinatorial Optimization. Polyhedra and Efficiency. Springer, Heidelberg (2004) 13. Shamir R., Sharan R., Tsur D.: Cluster graph modification problems. J. Discrete Applied Mathematics. 144, 1–2, 173–182 (2004) 14. Simanchev R.Yu.: Convex Polytope and Facet inequalities (Russian). Omsk State Univer- sity, Omsk (1999) 15. Simanchev R.Yu., Urazova I.V.: Scheduling unit-time jobs on parallel processors polytope (Russian). Diskretnyi Analiz i Issledovanie Operatsii. 18(11), 85–97 (2011) 16. Simanchev, R.Yu., Urazova, I.V.: On the Polytope Faces of the Graph Approximation Problem. J. of Applied and Industrial Mathematics. vol. 9, 2, 283–291 (2015) 17. van Zuylen A., Williamson D.P.: Deterministic Pivoting Algorithms for Constrained Rank- ing and Clustering Problems. Mathematics of Operations Research. 34, 3, 594–620 (2009) 18. Zahn C. T.: Approximating symmetric relations by equivalence relations. J. of the Society for Industrial and Applied Mathematics. 12, 4, 840–847 (1964)