graph structure in polynomial systems chordal networks
play

Graph structure in polynomial systems: chordal networks Pablo A. - PowerPoint PPT Presentation

Graph structure in polynomial systems: chordal networks Pablo A. Parrilo Laboratory for Information and Decision Systems Electrical Engineering and Computer Science Massachusetts Institute of Technology Based on joint work with Diego Cifuentes


  1. Graph structure in polynomial systems: chordal networks Pablo A. Parrilo Laboratory for Information and Decision Systems Electrical Engineering and Computer Science Massachusetts Institute of Technology Based on joint work with Diego Cifuentes (MIT) Lund University - LCCC - June 2017 Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 1 / 26

  2. Background: structured polynomial systems Many application domains require the solution of large-scale systems of polynomial equations. Among others: robotics, power systems, chemical en- gineering, cryptography, etc. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 2 / 26

  3. Polynomial systems and graphs A polynomial system defined by m equations in n variables: f i ( x 0 , . . . , x n − 1 ) = 0 , i = 1 , . . . , m Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 3 / 26

  4. Polynomial systems and graphs A polynomial system defined by m equations in n variables: f i ( x 0 , . . . , x n − 1 ) = 0 , i = 1 , . . . , m Construct a graph G (“primal graph”) with n nodes: Nodes are variables { x 0 , . . . , x n − 1 } . For each equation, add a clique connecting the variables appearing in that equation Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 3 / 26

  5. Polynomial systems and graphs A polynomial system defined by m equations in n variables: f i ( x 0 , . . . , x n − 1 ) = 0 , i = 1 , . . . , m Construct a graph G (“primal graph”) with n nodes: Nodes are variables { x 0 , . . . , x n − 1 } . For each equation, add a clique connecting the variables appearing in that equation Example: I = � x 2 x 2 x 2 x 3 � 0 x 1 x 2 + 2 x 1 + 1 , 1 + x 2 , x 1 + x 2 , Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 3 / 26

  6. Questions “Abstracted” the polynomial system to a (hyper)graph. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 4 / 26

  7. Questions “Abstracted” the polynomial system to a (hyper)graph. Can the graph structure help solve this system? For instance, to optimize, or to compute Groebner bases? Or, perhaps we can do something better ? Preserve graph (sparsity) structure? Complexity aspects? Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 4 / 26

  8. (Hyper)Graphical modelling Pervasive idea in many areas, in particular: numerical linear algebra, graphical models, constraint satisfaction, database theory, . . . Key notions: chordality and treewidth. Many names: Arnborg, Beeri/Fagin/Maier/Yannakakis, Blair/Peyton, Bodlaender, Courcelle, Dechter, Freuder, Lauritzen/Spiegelhalter, Pearl, Robertson/Seymour, . . . Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 5 / 26

  9. (Hyper)Graphical modelling Pervasive idea in many areas, in particular: numerical linear algebra, graphical models, constraint satisfaction, database theory, . . . Key notions: chordality and treewidth. Many names: Arnborg, Beeri/Fagin/Maier/Yannakakis, Blair/Peyton, Bodlaender, Courcelle, Dechter, Freuder, Lauritzen/Spiegelhalter, Pearl, Robertson/Seymour, . . . Remarkably (AFAIK) almost no work in computational algebraic geometry exploits this structure. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 5 / 26

  10. (Hyper)Graphical modelling Pervasive idea in many areas, in particular: numerical linear algebra, graphical models, constraint satisfaction, database theory, . . . Key notions: chordality and treewidth. Many names: Arnborg, Beeri/Fagin/Maier/Yannakakis, Blair/Peyton, Bodlaender, Courcelle, Dechter, Freuder, Lauritzen/Spiegelhalter, Pearl, Robertson/Seymour, . . . Remarkably (AFAIK) almost no work in computational algebraic geometry exploits this structure. Reasonably well-known in discrete (0/1) optimization, what happens in the continuous side? (e.g., Waki et al., Lasserre, Bienstock, Vandenberghe, Lavaei, etc) Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 5 / 26

  11. Chordality Let G be a graph with vertices x 0 , . . . , x n − 1 . A vertex ordering x 0 > x 1 > · · · > x n − 1 is a perfect elimination ordering if for all ℓ , the set X ℓ := { x ℓ } ∪ { x m : x m is adjacent to x ℓ , x ℓ > x m } is such that the restriction G | X ℓ is a clique. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 6 / 26

  12. Chordality Let G be a graph with vertices x 0 , . . . , x n − 1 . A vertex ordering x 0 > x 1 > · · · > x n − 1 is a perfect elimination ordering if for all ℓ , the set X ℓ := { x ℓ } ∪ { x m : x m is adjacent to x ℓ , x ℓ > x m } is such that the restriction G | X ℓ is a clique. A graph is chordal if it has a perfect elimination ordering. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 6 / 26

  13. Chordality Let G be a graph with vertices x 0 , . . . , x n − 1 . A vertex ordering x 0 > x 1 > · · · > x n − 1 is a perfect elimination ordering if for all ℓ , the set X ℓ := { x ℓ } ∪ { x m : x m is adjacent to x ℓ , x ℓ > x m } is such that the restriction G | X ℓ is a clique. A graph is chordal if it has a perfect elimination ordering. (Equivalently, in numerical linear algebra: Cholesky factorization has no “fill-in”) Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 6 / 26

  14. Chordality, treewidth, and a meta-theorem A chordal completion of G is a chordal graph with the same vertex set as G , and which contains all edges of G . Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 7 / 26

  15. Chordality, treewidth, and a meta-theorem A chordal completion of G is a chordal graph with the same vertex set as G , and which contains all edges of G . The treewidth of a graph is the clique number (minus one) of its smallest chordal completion. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 7 / 26

  16. Chordality, treewidth, and a meta-theorem A chordal completion of G is a chordal graph with the same vertex set as G , and which contains all edges of G . The treewidth of a graph is the clique number (minus one) of its smallest chordal completion. Informally, treewidth quantitatively measures how “tree-like” a graph is. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 7 / 26

  17. Chordality, treewidth, and a meta-theorem A chordal completion of G is a chordal graph with the same vertex set as G , and which contains all edges of G . The treewidth of a graph is the clique number (minus one) of its smallest chordal completion. Informally, treewidth quantitatively measures how “tree-like” a graph is. Meta-theorem: NP-complete problems are “easy” on graphs of small treewidth. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 7 / 26

  18. (Simple) example: stable set on trees Given a graph, a stable (or independent ) set is a subset of vertices, such that no two are pairwise neighbors. STABLE SET problem: Compute a stable set of maximum cardinality. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 8 / 26

  19. (Simple) example: stable set on trees Given a graph, a stable (or independent ) set is a subset of vertices, such that no two are pairwise neighbors. STABLE SET problem: Compute a stable set of maximum cardinality. For general graphs, NP-complete. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 8 / 26

  20. (Simple) example: stable set on trees Given a graph, a stable (or independent ) set is a subset of vertices, such that no two are pairwise neighbors. STABLE SET problem: Compute a stable set of maximum cardinality. For general graphs, NP-complete. On trees, linear-time solvable! Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 8 / 26

  21. (Simple) example: stable set on trees Given a graph, a stable (or independent ) set is a subset of vertices, such that no two are pairwise neighbors. STABLE SET problem: Compute a stable set of maximum cardinality. For general graphs, NP-complete. On trees, linear-time solvable! Fix a root, and solve this recursion starting from the leaves: � � S ( i ) = max( S ( j ) , 1 + S ( j ) ) , j ∈ children ( i ) j ∈ grandchildren ( i ) S (leaf) = 1 , where S ( i ) represents the size of the largest independent set of the corresponding subtree. Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 8 / 26

  22. Bad news? (I) Recall the subset sum problem, with data A = { a 1 , . . . , a n } ⊂ Z . Is there a subset of A that adds up to 0? Letting s i be the partial sums, we can write a polynomial system: 0 = s 0 0 = ( s i − s i − 1 )( s i − s i − 1 − a i ) 0 = s n The graph associated with these equations is a path (treewidth=1) s 0 — s 1 — s 2 — · · · — s n But, subset sum is NP-complete... :( Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 9 / 26

  23. Bad news? (II) For linear equations, “good” elimination preserves graph structure (perfect!) Cifuentes, Parrilo (MIT) Graph structure in polynomial systems LCCC 2017 10 / 26

Recommend


More recommend