Lattice polygons. P : lattice polygon in R 2 (vertices Z 2, no self-intersections)

Similar documents
Enumerating integer points in polytopes: applications to number theory. Matthias Beck San Francisco State University math.sfsu.

Ehrhart Positivity. Federico Castillo. December 15, University of California, Davis. Joint work with Fu Liu. Ehrhart Positivity

Computing the continuous discretely: The magic quest for a volume

Combinatorial Reciprocity Theorems

10 Years BADGeometry: Progress and Open Problems in Ehrhart Theory

The partial-fractions method for counting solutions to integral linear systems

Higher Spin Alternating Sign Matrices

h -polynomials of dilated lattice polytopes

and Other Combinatorial Reciprocity Instances

Combinatorial Reciprocity Theorems

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices

Nonlinear Discrete Optimization

Generalized Ehrhart polynomials

A Geometric Approach to Graph Isomorphism

Real-rooted h -polynomials

Classification of Ehrhart polynomials of integral simplices

High-dimensional permutations

MAXIMAL PERIODS OF (EHRHART) QUASI-POLYNOMIALS

Hodge theory for combinatorial geometries

A New Approximation Algorithm for the Asymmetric TSP with Triangle Inequality By Markus Bläser

3. Vector spaces 3.1 Linear dependence and independence 3.2 Basis and dimension. 5. Extreme points and basic feasible solutions

Tropical decomposition of symmetric tensors

A Questionable Distance-Regular Graph

Acyclic Digraphs arising from Complete Intersections

Ehrhart Polynomials of Zonotopes

An Introduction to Transversal Matroids

arxiv: v1 [math.co] 10 Aug 2016

Characterizations of Strongly Regular Graphs: Part II: Bose-Mesner algebras of graphs. Sung Y. Song Iowa State University

THE MIRROR SYMMETRY CONJECTURE FOR NETWORKS ON SURFACES

Dependencies in discrete tomography. Rob Tijdeman, Milan, Italy, 20 April

The Arithmetic of Graph Polynomials. Maryam Farahmand. A dissertation submitted in partial satisfaction of the. requirements for the degree of

Tropicalizations of Positive Parts of Cluster Algebras The conjectures of Fock and Goncharov David Speyer

A Finite Calculus Approach to Ehrhart Polynomials. Kevin Woods, Oberlin College (joint work with Steven Sam, MIT)

arxiv:math/ v2 [math.co] 9 Oct 2008

CS675: Convex and Combinatorial Optimization Fall 2016 Combinatorial Problems as Linear and Convex Programs. Instructor: Shaddin Dughmi

Auerbach bases and minimal volume sufficient enlargements

Linear equations in linear algebra

Linear algebra and applications to graphs Part 1

Gordon solutions. n=1 1. p n

arxiv: v3 [math.co] 1 Oct 2018

High-dimensional permutations and discrepancy

Chordal structure in computer algebra: Permanents

A 2 G 2 A 1 A 1. (3) A double edge pointing from α i to α j if α i, α j are not perpendicular and α i 2 = 2 α j 2

Degrees and Trees. Garth Isaak Lehigh University. 47th SEICCGTC at FAU, March Acknowledgements to: Kathleen Ryan,...

Mathematics 530. Practice Problems. n + 1 }

Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5

APPENDIX: MATHEMATICAL INDUCTION AND OTHER FORMS OF PROOF

Sampling Contingency Tables

Math 40510, Algebraic Geometry

Reconstruction and Higher Dimensional Geometry

A thesis presented to the faculty of San Francisco State University In partial fulfilment of The Requirements for The Degree

Ehrhart polynome: how to compute the highest degree coefficients and the knapsack problem.

LMI MODELLING 4. CONVEX LMI MODELLING. Didier HENRION. LAAS-CNRS Toulouse, FR Czech Tech Univ Prague, CZ. Universidad de Valladolid, SP March 2009

Algebraic Methods in Combinatorics

The problematic art of counting

AN INTRODUCTION TO TORIC SURFACES

Ch. 5 Determinants. Ring Determinant functions Existence, Uniqueness and Properties

PARKING FUNCTIONS. Richard P. Stanley Department of Mathematics M.I.T Cambridge, MA

ON THE RANK OF A TROPICAL MATRIX

21-301, Spring 2019 Homework 4 Solutions

ACO Comprehensive Exam October 18 and 19, Analysis of Algorithms

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

Algebraic Methods in Combinatorics

Analogs of Hodge Riemann relations

Chapter 3. Vector spaces

Preliminaries and Complexity Theory

Outline. Some Reflection Group Numerology. Root Systems and Reflection Groups. Example: Symmetries of a triangle. Paul Renteln

Minkowski Length of Lattice Polytopes

Integer Programming Methods LNMB

Algebraic and Geometric ideas in the theory of Discrete Optimization

3.7 Cutting plane methods

Some notes on Coxeter groups

3. Linear Programming and Polyhedral Combinatorics

Math 530 Lecture Notes. Xi Chen

Matrix compositions. Emanuele Munarini. Dipartimento di Matematica Politecnico di Milano

ACYCLIC DIGRAPHS GIVING RISE TO COMPLETE INTERSECTIONS

Alternative Characterization of Ergodicity for Doubly Stochastic Chains

Toric Fiber Products

The Ehrhart polynomial of the Birkhoff polytope 1

Pigeonhole Principle and Ramsey Theory

Convex Geometry. Otto-von-Guericke Universität Magdeburg. Applications of the Brascamp-Lieb and Barthe inequalities. Exercise 12.

LEVEL MATRICES G. SEELINGER, P. SISSOKHO, L. SPENCE, AND C. VANDEN EYNDEN

1 Integer Decomposition Property

H. Schubert, Kalkül der abzählenden Geometrie, 1879.

Refined Inertia of Matrix Patterns

0 Sets and Induction. Sets

Ma/CS 6b Class 20: Spectral Graph Theory

Integer Linear Programs

Albert W. Marshall. Ingram Olkin Barry. C. Arnold. Inequalities: Theory. of Majorization and Its Applications. Second Edition.

PRIMARY DECOMPOSITION FOR THE INTERSECTION AXIOM

CSL361 Problem set 4: Basic linear algebra

Chapter 1. Preliminaries

Definition For a set F, a polynomial over F with variable x is of the form

Six Little Squares and How Their Numbers Grow

CS675: Convex and Combinatorial Optimization Fall 2014 Combinatorial Problems as Linear Programs. Instructor: Shaddin Dughmi

Toric Varieties in Coding Theory. Math in the Mountains Tutorial

Bichain graphs: geometric model and universal graphs

A Survey of Parking Functions

Discrete Optimization 23

Algebraic Models in Different Fields

Transcription:

Lattice polygons P : lattice polygon in R 2 (vertices Z 2, no self-intersections)

A, I, B A = area of P I = # interior points of P (= 4) B = #boundary points of P (= 10)

Pick s theorem Georg Alexander Pick (1859 1942) A = 2I + B 2 2 = 2 4 + 10 2 2 = 9

Higher dimensions? Pick s theorem (seemingly) fails in higher dimensions. Example. Let T 1 and T 2 be the tetrahedra with vertices v(t 1 ) = {(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1)} v(t 2 ) = {(0, 0, 0), (1, 1, 0), (1, 0, 1), (0, 1, 1)}.

Two tetrahedra

Two tetrahedra Then I(T 1 ) = I(T 2 ) = 0 B(T 1 ) = B(T 2 ) = 4 vol(t 1 ) = 1/6, vol(t 2 ) = 1/3.

Dilation Let P be a convex polytope (convex hull of a finite set of points) in R d. For n 1, let np = {nα : α P}.

Dilation Let P be a convex polytope (convex hull of a finite set of points) in R d. For n 1, let np = {nα : α P}. P 3P

i(p, n) Let i(p, n) = #(np Z d ) = #{α P : nα Z d }, the number of lattice points in np.

ī(p, n) Similarly let P = interior of P = P P ī(p, n) = #(np Z d ) = #{α P : nα Z d }, the number of lattice points in the interior of np.

An example P 3P i(p,n) = (n + 1) 2 ī(p,n) = (n 1) 2 = i(p, n)

Reeve s theorem lattice polytope: polytope with integer vertices Theorem (Reeve, 1957). Let P be a three-dimensional lattice polytope. Then the volume V (P) is a certain (explicit) function of i(p, 1), ī(p, 1), and i(p, 2).

The Ehrhart-Macdonald theorem Theorem (Ehrhart 1962, Macdonald 1963). Let P = lattice polytope in R N, dim P = d. Then i(p,n) is a polynomial (the Ehrhart polynomial of P) in n of degree d.

The Ehrhart-Macdonald theorem Theorem (Ehrhart 1962, Macdonald 1963). Let P = lattice polytope in R N, dim P = d. Then i(p,n) is a polynomial (the Ehrhart polynomial of P) in n of degree d. Note. Eugène Ehrhart (1906 2000): taught at lycées in France, received Ph.D. in 1966.

Reciprocity, volume Moreover, i(p, 0) = 1 ī(p,n) = ( 1) d i(p, n), n > 0 (reciprocity).

Reciprocity, volume Moreover, If d = N then i(p, 0) = 1 ī(p,n) = ( 1) d i(p, n), n > 0 (reciprocity). i(p,n) = V (P)n d + lower order terms, where V (P) is the volume of P.

Reciprocity, volume Moreover, If d = N then i(p, 0) = 1 ī(p,n) = ( 1) d i(p, n), n > 0 (reciprocity). i(p,n) = V (P)n d + lower order terms, where V (P) is the volume of P. ( relative volume for d < N)

Generalized Pick s theorem Corollary (generalized Pick s theorem). Let P R d and dim P = d. Knowing any d of i(p,n) or ī(p,n) for n > 0 determines V (P).

Generalized Pick s theorem Corollary (generalized Pick s theorem). Let P R d and dim P = d. Knowing any d of i(p,n) or ī(p,n) for n > 0 determines V (P). Proof. Together with i(p, 0) = 1, this data determines d + 1 values of the polynomial i(p,n) of degree d. This uniquely determines i(p, n) and hence its leading coefficient V (P).

Reeve s theorem redux Example. When d = 3, V (P) is determined by i(p, 1) = #(P Z 3 ) i(p, 2) = #(2P Z 3 ) ī(p, 1) = #(P Z 3 ), which gives Reeve s theorem.

The Birkhoff polytope Example (magic squares). Let B M R M M be the Birkhoff polytope of all M M doubly-stochastic matrices A = (a ij ), i.e., i a ij 0 a ij = 1 (column sums 1) a ij = 1 (row sums 1). j

Integer points in B M Note. B = (b ij ) nb M Z M M if and only if b ij N = {0, 1, 2,... } b ij = n b ij = n. i j

Integer points in B M Note. B = (b ij ) nb M Z M M if and only if b ij N = {0, 1, 2,... } b ij = n b ij = n. i 2 1 0 4 3 1 1 2 1 3 2 1 1 2 4 0 j (M = 4, n = 7)

H M (n) H M (n) := #{M M N-matrices, line sums n} = i(b M,n).

H M (n) H M (n) := #{M M N-matrices, line sums n} = i(b M,n). H 1 (n) = 1 H 2 (n) = n + 1

H M (n) H M (n) := #{M M N-matrices, line sums n} = i(b M,n). H 1 (n) = 1 H 2 (n) = n + 1 [ a n a n a a ], 0 a n.

More examples H 3 (n) = ( ) n + 2 4 + ( ) n + 3 4 + ( ) n + 4 4 M 0 (MacMahon) H M (0) = 1 H M (1) = M! (permutation matrices) H M (2) xm M! 2 = e x/2 1 x (Anand-Dumir-Gupta, 1966)

Anand-Dumir-Gupta conjecture Theorem (Birkhoff-von Neumann). The vertices of B M consist of the M! M M permutation matrices. Hence B M is a lattice polytope.

Anand-Dumir-Gupta conjecture Theorem (Birkhoff-von Neumann). The vertices of B M consist of the M! M M permutation matrices. Hence B M is a lattice polytope. Corollary (Anand-Dumir-Gupta conjecture). H M (n) is a polynomial in n (of degree (M 1) 2 ).

H 4 (n) Example. H 4 (n) = 1 11340 ( 11n 9 + 198n 8 + 1596n 7 +7560n 6 + 23289n 5 + 48762n 5 + 70234n 4 + 68220n 2 +40950n + 11340).

H 4 (n) Example. H 4 (n) = 1 11340 ( 11n 9 + 198n 8 + 1596n 7 +7560n 6 + 23289n 5 + 48762n 5 + 70234n 4 + 68220n 2 +40950n + 11340). Open: positive coefficients

Positive magic squares Reciprocity ±H M ( n) = #{M M matrices B of positive integers, line sum n}. But every such B can be obtained from an M M matrix A of nonnegative integers by adding 1 to each entry.

Reciprocity for magic squares Corollary. H M ( 1) = H M ( 2) = = H M ( M + 1) = 0 H M ( M n) = ( 1) M 1 H M (n) (greatly reduces computation)

Reciprocity for magic squares Corollary. H M ( 1) = H M ( 2) = = H M ( M + 1) = 0 H M ( M n) = ( 1) M 1 H M (n) (greatly reduces computation) Applications e.g. to statistics (contingency tables) by Diaconis, et al.

Zeros (roots) of H 9 (n) Zeros of H_9(n) 3 2 1 8 6 4 2 0 1 2 3

Explicit calculation For what polytopes P can i(p,n) be explicitly calculated or related to other interesting mathematics? Main topic of subsequent talks.

Zonotopes Let v 1,...,v k R d. The zonotope Z(v 1,...,v k ) generated by v 1,...,v k : Z(v 1,..., v k ) = {λ 1 v 1 + + λ k v k : 0 λ i 1}

An example Example. v 1 = (4, 0), v 2 = (3, 1), v 3 = (1, 2) (1,2) (3,1) (0,0) (4,0)

i(z, 1) Theorem. Let where v i Z d. Then Z = Z(v 1,...,v k ) R d, i(z, 1) = X h(x), where X ranges over all linearly independent subsets of {v 1,...,v k }, and h(x) is the gcd of all j j minors (j = #X) of the matrix whose rows are the elements of X.

An example Example. v 1 = (4, 0), v 2 = (3, 1), v 3 = (1, 2) i(z, 1) = 4 0 3 1 + 4 0 1 2 + 3 1 1 2 +gcd(4, 0) + gcd(3, 1) +gcd(1, 2) + det( ) = 4 + 8 + 5 + 4 + 1 + 1 + 1 = 24.

Decomposition of a zonotope i(z, 1) = 4 0 3 1 + 4 0 1 2 + 3 1 1 2 +gcd(4, 0) + gcd(3, 1) +gcd(1, 2) + det( ) = 4 + 8 + 5 + 4 + 1 + 1 + 1 = 24.

Decomposition of a zonotope i(z, 1) = 4 0 3 1 + 4 0 1 2 + 3 1 1 2 +gcd(4, 0) + gcd(3, 1) +gcd(1, 2) + det( ) = 4 + 8 + 5 + 4 + 1 + 1 + 1 = 24.

Ehrhart polynomial of a zonotope Theorem. Let Z = Z(v 1,...,v k ) R d where v i Z d. Then i(z,n) = X h(x)n #X, where X ranges over all linearly independent subsets of {v 1,...,v k }, and h(x) is the gcd of all j j minors (j = #X) of the matrix whose rows are the elements of X.

Ehrhart polynomial of a zonotope Theorem. Let Z = Z(v 1,...,v k ) R d where v i Z d. Then i(z,n) = X h(x)n #X, where X ranges over all linearly independent subsets of {v 1,...,v k }, and h(x) is the gcd of all j j minors (j = #X) of the matrix whose rows are the elements of X. Corollary. The coefficients of i(z, n) are nonnegative integers.

Ordered degree sequences Let G be a graph (with no loops or multiple edges) on the vertex set V (G) = {1, 2,...,n}. Let d i = degree (# incident edges) of vertex i. Define the ordered degree sequence d(g) of G by d(g) = (d 1,...,d n ).

An example Example. d(g) = (2, 4, 0, 3, 2, 1) 1 2 3 4 5 6

Number of distinct d(g) Let f(n) be the number of distinct d(g), where V (G) = {1, 2,...,n}.

Number of distinct d(g) Let f(n) be the number of distinct d(g), where V (G) = {1, 2,...,n}. Example. If n 3, all d(g) are distinct, so f(1) = 1, f(2) = 2 1 = 2, f(3) = 2 3 = 8.

f(4) For n 4 we can have G H but d(g) = d(h), e.g., 1 2 1 2 1 2 3 4 3 4 3 4

f(4) For n 4 we can have G H but d(g) = d(h), e.g., 1 2 1 2 1 2 3 4 3 4 3 4 In fact, f(4) = 54 < 2 6 = 64.

The polytope of degree sequences Let conv denote convex hull, and D n = conv{d(g) : V (G) = {1,...,n}} R n, the polytope of degree sequences (Perles, Koren).

The polytope of degree sequences Let conv denote convex hull, and D n = conv{d(g) : V (G) = {1,...,n}} R n, the polytope of degree sequences (Perles, Koren). Easy fact. Let e i be the ith unit coordinate vector in R n. E.g., if n = 5 then e 2 = (0, 1, 0, 0, 0). Then D n = Z(e i + e j : 1 i < j n).

The Erdős-Gallai theorem Theorem (Erdős-Gallai). Let α = (a 1,...,a n ) Z n. Then α = d(g) for some G if and only if α D n a 1 + a 2 + + a n is even.

A generating function Fiddling around leads to: Theorem. Let F(x) = n 0 f(n) xn n! Then = 1 + 1x + 2 x2 2! + 8x3 3! + 54x4 4! +.

A nice formula ( F(x) = 1 1 + 2 n n xn 2 n! n 1 ( 1 ) n 1 xn (n 1) n! n 1 ) 1/2 + 1 ] exp n 1 n n 2 xn n!

A bad Ehrhart polynomial Let P denote the tetrahedron with vertices (0, 0, 0), (1, 0, 0), (0, 1, 0), (1, 1, 13). Then i(p,n) = 13 6 n3 + n 2 1 6 n + 1.

A bad Ehrhart polynomial Let P denote the tetrahedron with vertices (0, 0, 0), (1, 0, 0), (0, 1, 0), (1, 1, 13). Then i(p,n) = 13 6 n3 + n 2 1 6 n + 1. Thus in general the coefficients of Ehrhart polynomials are not nice. Is there a better basis?

Diagram of the bad tetrahedron y x z

The h-vector of i(p, n) Let P be a lattice polytope of dimension d. Since i(p,n) is a polynomial of degree d taking Z Z, h i Z such that n 0 i(p,n)x n = h 0 + h 1 x + + h d x d (1 x) d+1.

The h-vector of i(p, n) Let P be a lattice polytope of dimension d. Since i(p,n) is a polynomial of degree d taking Z Z, h i Z such that n 0 Definition. Define i(p,n)x n = h 0 + h 1 x + + h d x d (1 x) d+1. h(p) = (h 0,h 1,...,h d ), the h-vector of P (as an integral polytope).

An example Example. Recall i(b 4,n) = 1 11340 (11n9 +198n 8 + 1596n 7 + 7560n 6 + 23289n 5 +48762n 5 + 70234n 4 + 68220n 2 +40950n + 11340).

An example Example. Recall i(b 4,n) = 1 11340 (11n9 Then +198n 8 + 1596n 7 + 7560n 6 + 23289n 5 +48762n 5 + 70234n 4 + 68220n 2 +40950n + 11340). h(b 4 ) = (1, 14, 87, 148, 87, 14, 1, 0, 0, 0).

Properties of h(p) h(b 4 ) = (1, 14, 87, 148, 87, 14, 1, 0, 0, 0). Always h 0 = 1. Trailing 0 s i(b 4, 1) = i(b 4, 2) = i(b 4, 3) = 0. Palindromic property i(b 4, n 4) = ±i(b 4,n).

Nonnegativity and monotonicity Theorem A (nonnegativity), McMullen, RS: h i 0

Nonnegativity and monotonicity Theorem A (nonnegativity), McMullen, RS: h i 0 Theorem B (monotonicity), RS: If P and Q are lattice polytopes and Q P, then for all i: h i (Q) h i (P). B A: take Q =.

Proofs Both theorems can be proved geometrically. There are also elegant algebraic proofs based on commutative algebra. Related to toric varieties.

Fractional lattice polytopes Example. Let S M (n) denote the number of symmetric M M matrices of nonnegative integers, every row and column sum n. Then S 3 (n) = { 1 8 (2n3 + 9n 2 + 14n + 8), n even 1 8 (2n3 + 9n 2 + 14n + 7), n odd = 1 16 (4n3 + 18n 2 + 28n + 15 + ( 1) n ).

Fractional lattice polytopes Example. Let S M (n) denote the number of symmetric M M matrices of nonnegative integers, every row and column sum n. Then S 3 (n) = { 1 8 (2n3 + 9n 2 + 14n + 8), n even 1 8 (2n3 + 9n 2 + 14n + 7), n odd = 1 16 (4n3 + 18n 2 + 28n + 15 + ( 1) n ). Why a different polynomial depending on n modulo 2?

The symmetric Birkhoff polytope T M : the polytope of all M M symmetric doubly-stochastic matrices.

The symmetric Birkhoff polytope T M : the polytope of all M M symmetric doubly-stochastic matrices. Easy fact: S M (n) = # ( nt M Z M M) = i(t M,n).

The symmetric Birkhoff polytope T M : the polytope of all M M symmetric doubly-stochastic matrices. Easy fact: S M (n) = # ( nt M Z M M) = i(t M,n). Fact: vertices of T M have the form 1 2 (P + P t ), where P is a permutation matrix.

The symmetric Birkhoff polytope T M : the polytope of all M M symmetric doubly-stochastic matrices. Easy fact: S M (n) = # ( nt M Z M M) = i(t M,n). Fact: vertices of T M have the form 1 2 (P + P t ), where P is a permutation matrix. Thus if v is a vertex of T M then 2v Z M M.

S M (n) in general Theorem. There exist polynomials P M (n) and Q M (n) for which S M (n) = P M (n) + ( 1) n Q M (n), n 0. Moreover, deg P M (n) = ( M 2 ).

S M (n) in general Theorem. There exist polynomials P M (n) and Q M (n) for which S M (n) = P M (n) + ( 1) n Q M (n), n 0. Moreover, deg P M (n) = ( M 2 ). Difficult result (W. Dahmen and C. A. Micchelli, 1988): deg Q M (n) = { ( M 1 ) 2 1, M odd ) 1, M even. ( M 2 2