(Uniform) determinantal representations

Similar documents
On the complexity of the permanent in various computational models

MA 265 FINAL EXAM Fall 2012

Review problems for MA 54, Fall 2004.

Math 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination

The definition of a vector space (V, +, )

NOTES ON HYPERBOLICITY CONES

Permanent vs. Determinant

Announcements Wednesday, November 01

Linear Algebra Practice Problems

Lecture 2: Stable Polynomials

MATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL

Linear Algebra Practice Problems

Algebra Workshops 10 and 11

Geometric Complexity Theory and Matrix Multiplication (Tutorial)

Math 250B Final Exam Review Session Spring 2015 SOLUTIONS

MATH 369 Linear Algebra

2018 Fall 2210Q Section 013 Midterm Exam II Solution

Oeding (Auburn) tensors of rank 5 December 15, / 24

Math Linear Algebra Final Exam Review Sheet

Final Review Sheet. B = (1, 1 + 3x, 1 + x 2 ) then 2 + 3x + 6x 2

Plethysm and Kronecker Products

Final EXAM Preparation Sheet

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?

1 Counting spanning trees: A determinantal formula

Algebra 2 (2006) Correlation of the ALEKS Course Algebra 2 to the California Content Standards for Algebra 2

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #2 Solutions

Tensor decomposition and tensor rank

Announcements Wednesday, November 01

Discrete Math, Second Problem Set (June 24)

Equality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.

235 Final exam review questions

ARCS IN FINITE PROJECTIVE SPACES. Basic objects and definitions

Formalizing Randomized Matching Algorithms

Functions and Equations

Math 323 Exam 2 Sample Problems Solution Guide October 31, 2013

6 Orthogonal groups. O 2m 1 q. q 2i 1 q 2i. 1 i 1. 1 q 2i 2. O 2m q. q m m 1. 1 q 2i 1 i 1. 1 q 2i. i 1. 2 q 1 q i 1 q i 1. m 1.

Linear Algebra M1 - FIB. Contents: 5. Matrices, systems of linear equations and determinants 6. Vector space 7. Linear maps 8.

A monomial or sum of monomials

Linear Algebra 1 Exam 2 Solutions 7/14/3

DIAGONALIZATION. In order to see the implications of this definition, let us consider the following example Example 1. Consider the matrix

2.3. VECTOR SPACES 25

Calculating determinants for larger matrices

MATH 315 Linear Algebra Homework #1 Assigned: August 20, 2018

MAT 242 CHAPTER 4: SUBSPACES OF R n

Permanent versus determinant: not via saturations

Lecture 7: Schwartz-Zippel Lemma, Perfect Matching. 1.1 Polynomial Identity Testing and Schwartz-Zippel Lemma

Problem Set (T) If A is an m n matrix, B is an n p matrix and D is a p s matrix, then show

Final Examination 201-NYC-05 - Linear Algebra I December 8 th, and b = 4. Find the value(s) of a for which the equation Ax = b

SYMMETRIC DETERMINANTAL REPRESENTATIONS IN CHARACTERISTIC 2

Conceptual Questions for Review

Math 601 Solutions to Homework 3

Dimension. Eigenvalue and eigenvector

Math 113 Winter 2013 Prof. Church Midterm Solutions

MH1200 Final 2014/2015

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #1. July 11, 2013 Solutions

Geometric Modeling Summer Semester 2010 Mathematical Tools (1)

Linear Algebra. Min Yan

18.06 Professor Johnson Quiz 1 October 3, 2007

Eigenvectors. Prop-Defn

Math Exam 2, October 14, 2008

Topics in Graph Theory

1 Last time: determinants

MATH 320: PRACTICE PROBLEMS FOR THE FINAL AND SOLUTIONS

Math 396. An application of Gram-Schmidt to prove connectedness

(b) If a multiple of one row of A is added to another row to produce B then det(b) =det(a).

Multiplication of Polynomials

Remark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.

Online Exercises for Linear Algebra XM511

MATH 304 Linear Algebra Lecture 33: Bases of eigenvectors. Diagonalization.

OHSx XM511 Linear Algebra: Solutions to Online True/False Exercises

ADVANCED TOPICS IN ALGEBRAIC GEOMETRY

EIGENVALUES AND EIGENVECTORS 3

AN INTRODUCTION TO AFFINE TORIC VARIETIES: EMBEDDINGS AND IDEALS

PROBLEMS ON LINEAR ALGEBRA

The converse is clear, since

ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA

SYLLABUS. 1 Linear maps and matrices

Minimum Polynomials of Linear Transformations

Linear Algebra Primer

Symmetric Determinantal Representation of Weakly-Skew Circuits

Econ Slides from Lecture 8

Topic 2 Quiz 2. choice C implies B and B implies C. correct-choice C implies B, but B does not imply C

Rings. EE 387, Notes 7, Handout #10

The GCT chasm I. Ketan D. Mulmuley. The University of Chicago. The GCT chasm I p. 1

Linear and Bilinear Algebra (2WF04) Jan Draisma

Describing convex semialgebraic sets by linear matrix inequalities. Markus Schweighofer. Universität Konstanz

arxiv: v1 [math.ag] 25 Dec 2015

Unless otherwise specified, V denotes an arbitrary finite-dimensional vector space.

LINEAR ALGEBRA QUESTION BANK

Operations w/polynomials 4.0 Class:

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det

Abstract Vector Spaces

(II.B) Basis and dimension

Section 7.3: SYMMETRIC MATRICES AND ORTHOGONAL DIAGONALIZATION

Copositive matrices and periodic dynamical systems

Lecture 11: Finish Gaussian elimination and applications; intro to eigenvalues and eigenvectors (1)

Achievement Level Descriptors Mathematics Algebra 1

Ma/CS 6b Class 23: Eigenvalues in Regular Graphs

Transcription:

1 (Uniform) determinantal representations Jan Draisma Universität Bern October 216, Kolloquium Bern

Determinantal representations 2 R := C[x 1,..., x n ] and R d := {p R deg p d} Definition A determinantal representation of p R of size N is a matrix M R N N 1 with det(m) = p.

Determinantal representations 2 R := C[x 1,..., x n ] and R d := {p R deg p d} Definition A determinantal representation of p R of size N is a matrix M R N N 1 with det(m) = p. n = 1: companion matrices x 1 x 1 det...... x 1 a a 1 a n 2 a n 1 + a n x = a + a 1 x +... + a n x n

Determinantal representations 2 R := C[x 1,..., x n ] and R d := {p R deg p d} Definition A determinantal representation of p R of size N is a matrix M R N N 1 with det(m) = p. A bivariate example x 1 det y 1 a + bx + cy dx + ey f y = a + bx + cy + dx2 + exy + f y 2

Determinantal representations 2 R := C[x 1,..., x n ] and R d := {p R deg p d} Definition A determinantal representation of p R of size N is a matrix M R N N 1 with det(m) = p. A bivariate example x 1 det y 1 a + bx + cy dx + ey f y = a + bx + cy + dx2 + exy + f y 2 Determinantal representations always exist, but how small? the determinantal complexity dc(p) is the smallest N.

Why? 3

Motivation I: permanent versus determinant 4 If p has a determinantal representation M of small size N, then p can be evaluated efficiently using Gaussian elimination.

Motivation I: permanent versus determinant 4 If p has a determinantal representation M of small size N, then p can be evaluated efficiently using Gaussian elimination. Definition perm m := π S m x 1π(1) x mπ(m) is the m m permanent. Example 1 1 perm 3 1 1 1 1 1 = 3 counts perfect matchings:

Motivation I: permanent versus determinant 4 If p has a determinantal representation M of small size N, then p can be evaluated efficiently using Gaussian elimination. Definition perm m := π S m x 1π(1) x mπ(m) is the m m permanent. Example 1 1 perm 3 1 1 1 1 1 = 3 counts perfect matchings: Counting matchings in bipartite graphs is believed hard, so dc(perm m ) should be large!

Geometric complexity theory 5 Conjecture dc(perm m ) grows faster with m than any polynomial. [Valiant, 7s]

Geometric complexity theory 5 Conjecture dc(perm m ) grows faster with m than any polynomial. [Valiant, 7s] Best known bounds [Mignon-Ressayre 4, Grenet 12] m 2 2 dc(perm m ) 2m 1 [Alper-Bogart-Velasco 15: = 7 for m = 3]

Geometric complexity theory 5 Conjecture dc(perm m ) grows faster with m than any polynomial. [Valiant, 7s] Best known bounds [Mignon-Ressayre 4, Grenet 12] m 2 2 dc(perm m ) 2m 1 [Alper-Bogart-Velasco 15: = 7 for m = 3] Proof sketch of lower bound If ψ : C m m C N N affine-linear with det N (ψ(a)) = perm m (A), J := m + 1 1 1 1 1 1... 1 1 1 J perm m = ψ(j) det N =

Geometric complexity theory 5 Conjecture dc(perm m ) grows faster with m than any polynomial. [Valiant, 7s] Best known bounds [Mignon-Ressayre 4, Grenet 12] m 2 2 dc(perm m ) 2m 1 [Alper-Bogart-Velasco 15: = 7 for m = 3] Proof sketch of lower bound If ψ : C m m C N N affine-linear with det N (ψ(a)) = perm m (A), J := m + 1 1 1 1 1 1... ψ(j) 1 1 1 perm m = det N = q 1 (X) := quadratic part of perm m (J + X), form of rank m 2 q 2 (Y) := quadratic part of det N (ψ(j) + Y), form of rank 2N J

Geometric complexity theory 5 Conjecture dc(perm m ) grows faster with m than any polynomial. [Valiant, 7s] Best known bounds [Mignon-Ressayre 4, Grenet 12] m 2 2 dc(perm m ) 2m 1 [Alper-Bogart-Velasco 15: = 7 for m = 3] Proof sketch of lower bound If ψ : C m m C N N affine-linear with det N (ψ(a)) = perm m (A), J := m + 1 1 1 1 1 1... ψ(j) 1 1 1 perm m = det N = q 1 (X) := quadratic part of perm m (J + X), form of rank m 2 q 2 (Y) := quadratic part of det N (ψ(j) + Y), form of rank 2N Now q 1 (X) = q 2 (L(X)) where L linear part of ψ, so m 2 2N. J

Grenet s 2 m 1 construction 6 = 123 x 11 x 12 x 13 1 2 3 x 22 x 23 12 x 21 x 23 x 21 x 22 13 x 33 x 32 x 31 23 = 123 1 2 3 12 13 23 = 123 1 2 3 12 13 23 x 11 x 12 x 13 1 x 22 x 23 1 x 21 x 23 1 x 21 x 22 x 33 1 x 32 1 x 31 1 123 = x i j labels an arrow from an (i 1)-set to an i-set by adding j.

Geometric complexity theory 7 Theorem [Landsberg-Ressayre, 15] Grenet s representation is optimal among representations that preserve left multiplication with permutation and diagonal matrices.

Geometric complexity theory 7 Theorem [Landsberg-Ressayre, 15] Grenet s representation is optimal among representations that preserve left multiplication with permutation and diagonal matrices. GCT Programme [Mulmuley-Sohoni, 1-] Compare orbit closures X 1, X 2 of l N m perm m and det N inside the space of degree-n polynomials in N 2 variables under G = GL N 2; try to show that X 1 X 2 by showing that multiplicities of certain G-representations are higher in C[X 1 ] than in C[X 2 ] unless N is super-polynomial in m.

Geometric complexity theory 7 Theorem [Landsberg-Ressayre, 15] Grenet s representation is optimal among representations that preserve left multiplication with permutation and diagonal matrices. GCT Programme [Mulmuley-Sohoni, 1-] Compare orbit closures X 1, X 2 of l N m perm m and det N inside the space of degree-n polynomials in N 2 variables under G = GL N 2; try to show that X 1 X 2 by showing that multiplicities of certain G-representations are higher in C[X 1 ] than in C[X 2 ] unless N is super-polynomial in m. Theorem [Bürgisser-Ikenmeyer-Panova, 16] This approach does not work if higher than is restricted to 1 > (so-called occurrence obstructions).

Motivation II: Solving systems of equations 8 In numerics, solving a univariate equation p(x) = is often done by finding the eigenvalues of the companion matrix of p.

Motivation II: Solving systems of equations 8 In numerics, solving a univariate equation p(x) = is often done by finding the eigenvalues of the companion matrix of p. Proposal [Plestenjak-Hochstenbach, 16] To solve p(x, y) = q(x, y) = write p = det(a + xa 1 + ya 2 ) and q = det(b + xb 1 x + yb 2 ) and solve the two-parameter eigenvalue problem (A + xa 1 + ya 2 )u = and (B + xb 1 + yb 2 )v =.

Motivation II: Solving systems of equations 8 In numerics, solving a univariate equation p(x) = is often done by finding the eigenvalues of the companion matrix of p. Proposal [Plestenjak-Hochstenbach, 16] To solve p(x, y) = q(x, y) = write p = det(a + xa 1 + ya 2 ) and q = det(b + xb 1 x + yb 2 ) and solve the two-parameter eigenvalue problem (A + xa 1 + ya 2 )u = and (B + xb 1 + yb 2 )v =. translates to a joint pair of generalised eigenvalue problems: ( 1 x )w = and ( 2 y )w = where w = u v and = A 1 B 2 A 2 B 1, 1 = A 2 B A B 2, 2 = A B 1 A 1 B

Motivation II: Solving systems of equations 8 In numerics, solving a univariate equation p(x) = is often done by finding the eigenvalues of the companion matrix of p. Proposal [Plestenjak-Hochstenbach, 16] To solve p(x, y) = q(x, y) = write p = det(a + xa 1 + ya 2 ) and q = det(b + xb 1 x + yb 2 ) and solve the two-parameter eigenvalue problem (A + xa 1 + ya 2 )u = and (B + xb 1 + yb 2 )v =. translates to a joint pair of generalised eigenvalue problems: ( 1 x )w = and ( 2 y )w = where w = u v and = A 1 B 2 A 2 B 1, 1 = A 2 B A B 2, 2 = A B 1 A 1 B If the sizes are N, then i have size N 2, and solving takes (N 2 ) 3... (plane curves have det rep of size = deg, but harder to compute).

Determinantal complexity of general polynomials 9 Theorem [Boralevi-v Doornmalen-D-Hochstenbach-Plestenjak, 16] For n fixed, there exist C 1, C 2 such that a sufficiently general p R d has dc(p) C 1 d n/2 and any p R d has dc(p) C 2 d n/2.

Determinantal complexity of general polynomials 9 Theorem [Boralevi-v Doornmalen-D-Hochstenbach-Plestenjak, 16] For n fixed, there exist C 1, C 2 such that a sufficiently general p R d has dc(p) C 1 d n/2 and any p R d has dc(p) C 2 d n/2. For the upper bound, the determinantal representation can be chosen to depend bi-affine-linearly on x 1,..., x n and on the coefficients of p; these are uniform determinantal representations.

Determinantal complexity of general polynomials 9 Theorem [Boralevi-v Doornmalen-D-Hochstenbach-Plestenjak, 16] For n fixed, there exist C 1, C 2 such that a sufficiently general p R d has dc(p) C 1 d n/2 and any p R d has dc(p) C 2 d n/2. For the upper bound, the determinantal representation can be chosen to depend bi-affine-linearly on x 1,..., x n and on the coefficients of p; these are uniform determinantal representations. Proof of lower bound If sufficiently general p R d have dc(p) N, then the map det : R N N 1 R N contains R d in the closure of its image. Comparing ( ) n + d dimensions, find N 2 (n + 1) dim C R d =. n

Spaces connected to 1 1 Definition Given a nonzero subspace V R write V d := V R d. V is connected to 1 if V d+1 R 1 V d for all d.

Spaces connected to 1 1 Definition Given a nonzero subspace V R write V d := V R d. V is connected to 1 if V d+1 R 1 V d for all d. Example For n = 2, V spanned by these monomials is connected to 1: 1 x y x 2 xy y 2 x 3 x 2 y xy 2 y 3 x 4 x 3 y x 2 y 2 xy 3 y 4

Spaces connected to 1 1 Definition Given a nonzero subspace V R write V d := V R d. V is connected to 1 if V d+1 R 1 V d for all d. Example For n = 2, V spanned by these monomials is connected to 1: 1 x y x 2 xy y 2 x 3 x 2 y xy 2 y 3 x 4 x 3 y x 2 y 2 xy 3 y 4 Lemma V connected to 1, with basis 1 = f 1, f 2,..., f m of ascending degrees, write f i = i 1 j=1 l i j f j with l i j R 1. Then V = the span of the (m 1) (m 1)-subdeterminants of M(V) := l 21 1 l 31 l 32 1....... l m1 l m2 l m,m 1 1

First construction 11 Proposition Let V R be connected to 1, of dimension m, and such that R 1 V R d. Then there is a uniform determinantal representation of size m for the polynomials in R d.

First construction 11 Proposition Let V R be connected to 1, of dimension m, and such that R 1 V R d. Then there is a uniform determinantal representation of size m for the polynomials in R d. Example x 1 det y 1 a + bx + cy dx + ey f y = a + bx + cy + dx 2 + exy + f y 2 M(V) V : 1 x y x 2 xy y 2

First construction 11 Proposition Let V R be connected to 1, of dimension m, and such that R 1 V R d. Then there is a uniform determinantal representation of size m for the polynomials in R d. Example x 1 det y 1 a + bx + cy dx + ey f y = a + bx + cy + dx 2 + exy + f y 2 M(V) V : 1 x y x 2 xy y 2 Theorem [Hochstenbach-Plestenjak 16] For n = 2 there exist uniform det representations of size d2 4.

Analysis of first construction 12 V connected to 1 and R 1 V R d imply dim V 1 n ( n + d n )

Analysis of first construction 12 V connected to 1 and R 1 V R d imply dim V 1 n ( n + d n ) Proposition For fixed n, uniform determinantal representation of size dn n n!.

Analysis of first construction 12 V connected to 1 and R 1 V R d imply dim V 1 n ( n + d n ) Proposition For fixed n, uniform determinantal representation of size dn n n!. Construction uses the lattice of type A n 1 with generating matrix 2 1. 1 2........ 1 1 2 (David Madore, YouTube)

Analysis of first construction 12 V connected to 1 and R 1 V R d imply dim V 1 n ( n + d n ) Proposition For fixed n, uniform determinantal representation of size dn n n!. Construction uses the lattice of type A n 1 with generating matrix 2 1. 1 2........ 1 1 2 But the exponent of d is n rather than n/2. (David Madore, YouTube)

Second construction: divide and conquer! 13 Proposition Suppose V 1, V 2 R connected to 1 such that R 1 V 1 V 2 R d. Then there is a uniform det representation of degree-d polynomials of size 1 + dim V 1 + dim V 2.

Second construction: divide and conquer! 13 Proposition Suppose V 1, V 2 R connected to 1 such that R 1 V 1 V 2 R d. Then there is a uniform det representation of degree-d polynomials of size 1 + dim V 1 + dim V 2. M(V 1 ) Example x 1 x 1 det c c 1 c 2 y = i+ j 2 c i j x i y j c 1 c 11 1 y M(V c 2 1 2 ) T

Second construction: divide and conquer! 13 Proposition Suppose V 1, V 2 R connected to 1 such that R 1 V 1 V 2 R d. Then there is a uniform det representation of degree-d polynomials of size 1 + dim V 1 + dim V 2. M(V 1 ) Example x 1 x 1 det c c 1 c 2 y = i+ j 2 c i j x i y j c 1 c 11 1 y M(V c 2 1 2 ) T Can we find V 1, V 2, connected to 1, of dim dim R d such that (R 1 )V 1 V 2 R d?

A fractal 14 Can we find V 1, V 2, connected to 1, of dim growing like dim R d such that (R 1 )V 1 V 2 R d?

A fractal 14 Can we find V 1, V 2, connected to 1, of dim growing like dim R d such that (R 1 )V 1 V 2 R d? ( ) n/2 + d For n even, split variables V 1, V 2 of dimension. n/2

A fractal 14 Can we find V 1, V 2, connected to 1, of dim growing like dim R d such that (R 1 )V 1 V 2 R d? ( ) n/2 + d For n even, split variables V 1, V 2 of dimension. n/2 For odd n, find subsets A, A 1 (Z ) n, connected to, of dimension n 2 such that A + A 1 = Z n : - start with B := j= {, 1} 2 2 j ; - B 1 := 2B so that B + B 1 = Z ; - A i := B n i ; - connect to.

A fractal 14 Can we find V 1, V 2, connected to 1, of dim growing like dim R d such that (R 1 )V 1 V 2 R d? ( ) n/2 + d For n even, split variables V 1, V 2 of dimension. n/2 For odd n, find subsets A, A 1 (Z ) n, connected to, of dimension n 2 such that A + A 1 = Z n : - start with B := j= {, 1} 2 2 j ; - B 1 := 2B so that B + B 1 = Z ; - A i := B n i ; - connect to.

A fractal 14 Can we find V 1, V 2, connected to 1, of dim growing like dim R d such that (R 1 )V 1 V 2 R d? ( ) n/2 + d For n even, split variables V 1, V 2 of dimension. n/2 For odd n, find subsets A, A 1 (Z ) n, connected to, of dimension n 2 such that A + A 1 = Z n : - start with B := j= {, 1} 2 2 j ; - B 1 := 2B so that B + B 1 = Z ; - A i := B n i ; - connect to. Take V i spanned by the monomials with exponent vectors in A i.

Outlook 15 Theorem [Boralevi-v Doornmalen-D-Hochstenbach-Plestenjak, 16] For n fixed, there exist C 1, C 2 such that a sufficiently general p R d has dc(p) C 1 d n/2 and any p R d has dc(p) C 2 d n/2. Many questions remain: what are the best constants C 1, C 2? what about the regime where d is fixed and n runs? symmetric determinantal representations?

Outlook 15 Theorem [Boralevi-v Doornmalen-D-Hochstenbach-Plestenjak, 16] For n fixed, there exist C 1, C 2 such that a sufficiently general p R d has dc(p) C 1 d n/2 and any p R d has dc(p) C 2 d n/2. Many questions remain: what are the best constants C 1, C 2? what about the regime where d is fixed and n runs? symmetric determinantal representations? Thank you!

Outlook 15 Theorem [Boralevi-v Doornmalen-D-Hochstenbach-Plestenjak, 16] For n fixed, there exist C 1, C 2 such that a sufficiently general p R d has dc(p) C 1 d n/2 and any p R d has dc(p) C 2 d n/2. Many questions remain: what are the best constants C 1, C 2? what about the regime where d is fixed and n runs? symmetric determinantal representations? Thank you! Motivation III: hyperbolic polynomials If p = det(a + i x i A i ) with A i R N N symmetric and A positive definite, then the restriction of p to any line through has only real roots. For n = 2 the converse also holds (Helton-Vinnikov).