arxiv:math/ v1 [math.ag] 14 Feb 2007

Transcription

1 arxiv:math/ v [math.ag] 4 Feb 2007 SIMPLE EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA JEAN B. LASSERRE AND EDUARDO S. ZERON Abstract. Given z C n and A Z m n, we consider the problem of evaluating the counting function h(y;z) := P {z x x Z n ;Ax=y,x 0}. We provide an explicit expression for h(y; z) as well as an algorithm with possibly numerous but very simple calculations. In addition, we exhibit finitely many fixed convex cones of R n explicitly and exclusively defined by A such that for any y Z m, the sum h(y;z) can be obtained by a simple formula involving the evaluation of P z x over the integral points of those cones only. At last, we also provide an alternative (and different) formula from a decomposition of the generating function into simpler rational fractions, easy to invert. Keywords: Computational geometry; lattice polytopes.. introduction Consider the (not necessarily compact) polyhedron (.) Ω(y) := {x R n Ax = y; x 0}, with y Z m and A Z m n of maximal rank for n m; besides, given z C n, let h : Z m C be the counting function (.2) y h(y; z) := x Ω(y) Z n z x (where z x stands for k z k x k). The complex vector z C n may be chosen close enough to zero in order to ensure that h(y;z) is well defined even when Ω(y) is not compact. If Ω(y) is compact, then y h(y;z) provides us with the exact number of points in the set Ω(y) Z n by either evaluating h(y,), or even rounding h(y;z) up to the nearest integer when all the entries of z are close enough to one. Computation of h has attracted a lot of attention in recent years, from both theoretical and practical computation viewpoints. Barvinok and Pommersheim [4], Brion and Vergne [8], have provided nice exact (theoretical) formulas for h(y; z); see also Szenes and Vergne [5]. For instance, Barvinok considers z h(y;z) as the generating function (evaluated at z := e c C n ) of the indicator function x I Ω(y) Z n(x) of the set Ω(y) Z n and provides a decomposition into a sum of simpler generating functions associated with supporting cones (themselves having a signed decomposition into unimodular cones). We call this a primal approach because y is fixed, and one works

2 2 JEAN B. LASSERRE AND EDUARDO S. ZERON in the primal space R n in which Ω(y) is defined. Remarkably, Barvinok s counting algorithm which is implemented in the software LattE (see De Loera et al. [0]) runs in time polynomial in the problem size when the dimension n is fixed. The software developed by Verdoolaege [7] extends the LattE software to handle parametric polytopes. On the other hand, Brion and Vergne [8] consider the generating function H : C m C of y h(y;z), that is, (.3) w H(w) := y Z m h(y;z)w y = n k= z k w A k. They provide a generalized residue formula, and so obtain h(y; z) in closed form by inversion. We call this latter approach dual because z is fixed, and one works in the space C m of variables w associated with the m constraints Ax = y. As a resultof both primal and dualapproaches, h(y;z) is finally expressed as a weighted sum over the vertices of Ω(y). Similarly, Beck [5], and Beck, Diaz and Robins[6] provided a complete analysis based on residue techniques for the case of a tetrahedron (m = ). Despite its theoretical interest, Brion and Vergne s formula is not directly tractable because it contains many products with complex coefficients (roots of unity) which makes the formula difficult to evaluate numerically. However, in some cases, this formula can be exploited to yield an efficient algorithm as e.g. in [2] for flow polytopes, in [7] for transportation polytopes, and more generally when the matrix A is totally unimodular as in [9]. Finally, in [2, 3], we have provided two algorithms based on Cauchy residue techniques to invert H in (.3), and an alternative algebraic technique based on partial fraction expansion of H. A nice feature of the latter technique of [3] is to avoid computing residues. Contribution: Our contribution is twofold as it is concerned with both primal and dual approaches. On the primal side, we provide an explicit expression of h(y; z) and an algorithm which involves only elementary operations. It uses Brion s identity along with an explicit description of the supporting cones at the vertices of Ω(y). It also has a simple equivalent formulation as a (finite) group problem. Finally, we exhibit finitely many fixed convex cones of R n, explicitly and exclusively defined from A, such that for any y Z m, the sum h(y;z) is obtained by a simple formula which evaluates z x over the integral points of those cones only. On the dual side, we analyze the counting function h, via its generating function H in (.3). Inverting H is difficult in general, except if an appropriate expansion of H into simple fractions is available, as in e.g. [3]. In their landmark paper [8], Brion and Vergne provided a generalized residue formula which yields the generic expansion (.4) H(w) = J A g G Qg, [z k w A k] δ k, ρ g k q [z k w A k ] /q. k

3 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA 3 Here, J A whenever A is invertible, q is the smallest common multiple of all deta 0, ρ q = e 2πi/q is the q-root of unity, δ k, {0,/q}, and Q g, C. The finite group G has q m elements. The coefficients Q g, are difficult to evaluate. Our contribution is to expand H in (.3) in the form (.5) H(w) = [ z J A j j w A j ] R 2 (;z) u Z n m µ z η[,u ] w Aη[,u ], where: Z µ = {0,,...,µ }, µ = deta, each η[,u ] Z n and: (.6) z R 2 (;z) := [ ( ) ] z k z A A µ k. k/ Identity (.5) is a nontrivial simplification of the residue formula (.4) becausetheη[,u ] s aregiven explicitly. Andsothecoefficients oftherational fraction (.5) in w are very simple to evaluate with no root of unity involved (it can also be done symbolically); however this task can be tedious as for each J A one has deta n m terms η[,u ] to determine. But once determined, (.5) is easy to invert and provides h(y;z) for any y Z m. 2. Brion s decomposition 2.. Notation and definitions. The notation C, R and Z stand for the usual sets of complex, real and integer numbers, respectively. Moreover, the set of natural numbers {0,,2,...} is denoted by N, and for every natural number µ N, the finite set {0,,...,µ } of cardinality µ is denoted by Z µ. The notation B stands for the transpose of a matrix (or vector) B R s t ; and the kth column of the matrix B is denoted by B k := (B,k,...,B s,k ). When y = 0, the cone Ω(0) in (.) is convex, and its dual cone is given by, (2.) Ω(0) := {b R n b x 0 for every x Ω(0)}. Notice that Ω(0) R n if Ω(0) = {0}, which is the case if Ω(y) is compact. Definition. Let A Z m n be of maximal rank. An ordered set = {,..., m } of natural numbers is said to be a basis if it has cardinality = m, the sequence of inequalities < 2 < < m n holds, and the square [m m] submatrix : (2.2) A := [A A 2 A m ] is invertible. We denote the set of all bases by J A. Definition 2. Given a maximal rank matrix A Z m n, and any basis J A, the complementary matrices A Z m n and A Z m (n m) stand for [A k ] k and [A k ] k/, respectively. Similarly, given z C n, the complementary vectors z C m and z C n m stand for (z k ) k and (z k ) k/, respectively.

4 4 JEAN B. LASSERRE AND EDUARDO S. ZERON For each basis J A with associated matrix A Z m m, introduce the indicator function δ : Z m N defined by : { if A (2.3) y δ (y) := y Zm, 0 otherwise. Notice that δ is a multi-periodic function with periods A and µ := deta, meaning that δ (y+a q) = δ (y+µ q) = δ (y) for all y,q Z m. Finally, given a triplet (z,x,u) C n Z n R n, introduce the notation : (2.4) z x := z x zx 2 2 z s xn, z := max{ z, z 2,..., z n }, ln z := (ln(z ),ln(z 2 ),...,ln(z n )). Notice that z x = z x z x, for all bases J A and all z C n, x Z n Brion s decomposition. Let Ω(y) be the convex polyhedron in (.) with y Z m,a Z m n being of maximal rank, and let h : Z m C be the counting function in (.2), with z <. Obviously h(y; z) = 0 whenever the equation Ax = y has no solution x N n. The main idea is to decompose the function h following Brion s ideas. Given any convex rational polyhedron P R n, let [P] : R n {0,} be its characteristic function, and f[p] : C C its associated rational function, such that (2.5) z f[p,z] := x P Z n z x, holds whenever the sum converges absolutely. For every vertex V of P, define Co(P,V) R n to be the supporting cone of P at V. Then, Brion s formula yields the decomposition : (2.6) [P] = [Co(P,V)], vertices V modulo the group generated by the characteristic functions of convex polyhedra which contain affine lines. And so, (2.7) f[p,z] = f[co(p,v),z]. vertices V The above summation is formal because in general there is no z C n for which the series {z x x P Z n } and {z x x Co(P,V) Z n } converge absolutely for all vertices V. The notation E stands for the sum of all elements of a countable set E C. It is a complex number whenever the resulting series converges absolutely; otherwise it stands for a formal series.

5 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA 5 Example: LetP := [0,] RsothatCo(P,{0}) = [0,+ )andco(p,{}) = (,]. Simple enumeration yields f[p,z] = z 0 +z = +z, but one also has: f[p,z] = f[co(p,{0}),z]+f[(p,{}),z] = /( z)+z 2 /(z) = +z. 3. Computing h(y; z): A primal approach Let C(J A ) := {Ax x N n } R m be the cone generated by the columns of A, and for any basis J A, let C() R m be the cone generated by the columns A k with k. As A has maximal rank, C(J A ) is the union of all C(), J A. With any y C(J A ) associate the intersection of all cones C() that contain y. This defines a subdivision of C(J A ) into polyhedral cones. The interiors of the maximal subdivisions are called chambers. In each chamber γ, the polyhedron Ω(y) is simple, i.e. A y > 0 for all J A such that A y 0. For any chamber γ, define, (3.) B(J A,γ) := { J A γ C()}. The intersection of all C() with B(J A,γ) is the closure γ of γ. Back to our original problem, and setting P := Ω(y), the rational function f[p,z] is equal to h(y;z) in (.2) whenever z <. We next provide an explicit description of the rational function f[co(p, V), z] for every vertex V of P. Let δ be the function defined in (2.3), and let Z µ := {0,,...,µ } with µ := deta. A vector V R n is a vertex of P = Ω(y) if and only if there exists a basis J A such that : (3.2) V = A y 0 and V = 0, where V and V are given in Definition 2. Moreover, the supporting cone of P at the vertex V is described by : (3.3) Co(Ω(y),V) := {x R n Ax = y; x k 0 if V k = 0}. Let us now define the larger set (3.4) C(Ω(y),) := {x R n A x +A x = y; x 0}, so that Co(Ω(y),V) is a subcone of C(Ω(y),) for all bases J A and vertex V of Ω(y) which satisfy V = 0 (recall (3.2)). Besides, when V = 0 and y γ for some chamber γ, then C(Ω(y),) and Co(Ω(y),V) are identical because Ω(y) is a simple polytope, and so A y > 0 for all J A. Recall that A Z m n and A Z m (n m) stand for [A k ] k and [A k ] k/, respectively. Similarly, given a vector x Z n, the vectors x and x stand for (x k ) k and (x k ) k/ respectively. The following result is from [8, p. 88].

6 6 JEAN B. LASSERRE AND EDUARDO S. ZERON Proposition. Let y R m and let Ω(y) be as in (.), and let y γ for some chamber γ. Then, (3.5) [Ω(y)] = [C(Ω(y),)], B(J A,γ) modulo the group generated by the characteristic functions of convex polyhedra which contain affine lines. Proof. Using notation of [8, p. 87], define the linear mapping p : R n R m with p(x) = Ax, so that the polyhedra P (y) and Ω(y) are identical. Moreover, for every basis B(J A,γ), v : R m R n is the linear mapping: y [v (y)] = A y and [v (y)] = 0, y R m. Finally, for every ˆx R n with ˆx 0, ρ (ˆx) := ˆx v (Aˆx) satisfies, [ρ (ˆx)] = A A ˆx and [ρ (ˆx)] = ˆx. Therefore, the cone [v (y)+ρ (C)] in [8] is the set of points x R m such that x 0 and x = A (y A x ); and so this cone is just [C(Ω(y),)] in (3.4). Therefore a direct application of (3.2.) in [8, p. 88] yields (3.5). Theorem 2. Let y Z m, z C n with z <, and let y γ for some chamber γ. Recall the set of bases B(J A,γ) defined in (3.). With P := Ω(y), the rational function h defined in (.2) can be written: (3.6) h(y;z) = f[c(ω(y),),z] = R (y,;z) R 2 (;z), B(J A,γ) B(J A,γ) (3.7) (3.8) with and z R (y,;z) := z A y z R 2 (;z) := k/ u Z n m µ [ ( z k z A δ (y A u)z u, z A A u A k ) µ ]. The pair {R,R 2 } is well defined whenever z C n satisfies z k 0 and z k z A A k for every basis J A which does not contain the index k. Proof. By a direct application of Brion s theorem to the sum (3.5), the associated rational functions f[ω(y), z] and f[c(ω(y), ), z] satisfy: (3.9) h(y,z) = f[ω(y),z] = f[c(ω(y),),z]. B(J A,γ) Therefore, inordertoshow(3.6), oneonlyneedstoprovethat therational function R (y,;z) R 2 (;z) is equal to f[c(ω(y),),z], i.e., (3.0) R (y,;z) R 2 (;z) = {z x x C(Ω(y),) Z n },

7 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA 7 on the domain D = {z C n > z k z A A k for each k }. Notice that R 2 (;z) = k/ ( z k z A A k ) µ = = [ ] z µv k k = z µv, k/ v k N z A A k z µa A v v N n m on D. On the other hand, according to (3.4), the integer vector x Z n lies inside the cone C(P,V ) if and only if : x = A (y A x ), δ (y A x ) = and x = u+µ v, with u Z n m µ and v N n m. But from the definition (3.7) of R (y,;z) and z x = z x zx (3.) R (y,;z) R 2 (;z) = z A y u Z n m µ v N n m = za (y A x ), δ (y A u)z x z A A x = {z x x C(Ω(y),) Z n } = f[c(ω(y)),z], which is exactly (3.0). Notice that x = u + µ v, and so δ (y A u) = δ (y A x ) because of the definition (2.3) of δ. Finally, using (3.) in (3.9) yields that (3.6) holds whenever z < and R (y,;z) and R 2 (;z) are all well defined. Notice that R 2 is constant with respect to y, and from the definition (2.3) of δ, R is quasiperiodic with periods A and µ, meaning that (3.2) R (y +A q,;z) = R (y,;z)z q R (y +µ q,;z) = R (y,;z) ( z A and q) µ hold for all y,q Z m. Obviously, the more expensive part in calculating R 2 ( ) in (3.8) is to compute the determinant µ = deta. On the other hand, computing R ( ) in (3.7) may become quite expensive when µ is large, as one must evaluate µ n m terms, the cardinality of Zµ n m. However, as detailed below, a more careful analysis of (3.7) yields some simplifications. 3.. Simplifications via group theory. From the proof of Theorem 2, the closed forms (3.7) (3.8) for R ( ) and R 2 ( ) are deduced from (3.), i.e., R (y,;z) R 2 (;z) = z A y x Z n m δ (y A x )z x z A A x after setting x = u + µ v and recalling that δ (y) is a periodic function, i.e., δ (y +µ q) = δ (y) for all y,q Z m. However, we have not used yet,,

8 8 JEAN B. LASSERRE AND EDUARDO S. ZERON that δ (y +A q) = δ (y) as well. For every J A, consider the lattice : (3.3) Λ := j A j Z Z m, generated by the columns A j, j. The following quotient group G := Z m /Λ = Z m/ (3.4) A j Z j = {Ec[0,],Ec[2,],...,Ec[µ,]} iscommutative, withµ = deta elements(or, equivalenceclasses)ec[j,], and so, G is isomorphic to a finite Cartesian product of cyclic groups Z ηk, i.e., G = Zη Z η2 Z ηs. Obviously, µ = η η 2 η s, and so, G is isomorphic to the cyclic group Z µ whenever µ is a prime number. Actually, G = {0} whenever µ =. Notice that the Cartesian product Z η Z ηs can be seen as the integer space Z s modulo the vector η := (η,η 2,,η s ) N s. Hence, foreveryfinitecommutativegroupg, thereexistapositiveinteger s, a vector η N s with positive entries, and a group isomorphism, (3.5) g : G Z s mod η, where g (ξ) mod η means evaluating [g (ξ)] k mod [η ] k, for all indices k s. For every y Z m, there exists a unique equivalence class Ec[j y,] which contains y, and so we can define the following group epimorphism, (3.6) ĥ : Z m Z s mod η, y ĥ(y) := g (Ec[j y,]). On the other hand, the unit element of G is the equivalence class Ec[0,] which contains the origin, that is, Ec[0,] = {A q q Z m }. Hence, ĥ(y) = 0 if and only if there exists q Z m such that y = A q. We can then redefine the function δ as follows, { if ĥ (3.7) y δ (y) := (y) = 0, 0 otherwise, One also needs the following additional notation; given any matrix B Z m t, (3.8) ĥ (B) := [ĥ(b ) ĥ(b 2 ) ĥ(b t )] Z s t. And so, from (3.7), ĥ(y A u) ĥ(y) ĥ(a )u mod η. Finally, using (3.7) in (3.7), one obtains a simplified version of R ( ) in the form: (3.9) R (y,;z) = { z A y z u } u Zµ n m ;. ĥ (y) ĥ(a )u mod η z A A u

9 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA 9 Next, with q Z m fixed, ν q A q Zm for some integer ν q, if and only if ν q ĥ (q) = 0 mod η. If we set ν q = µ, then µ A q Zm, and µ ĥ (q) = 0 mod η, because G has µ = deta elements. Nevertheless, µ may not be the smallest positive integer with that property. So, given J A and k /, define ν k, to be order of ĥ(a k ). That is, ν k, is the smallest positive integer such that ν k, ĥ (A k ) = 0 mod η, or equivalently : (3.20) ν k, A A k Z m. Obviously ν k, µ. Moreover, µ is a multiple of ν k, for it is the order of an element in G. For example, the group Z 2 modulo η = ( 2 7) has 4 elements; and the elements b = ( ) 0, b2 = ( 0) and b3 = ( ) have respective orders : 2, 7 and 4. Notice that, 2b 7b 2 4b 3 0 mod η. But, 2b 3 2b 2 0 and 7b 3 b 0 mod η. The important observation is that δ (y ν k, A k q) = δ (y) for all q Z m and k /, which follows from (3.20) and (2.3). Thus, following step by step the proof of Theorem 2, we obtain: Corollary 3. Let y Z m, z C n with z <, and let y γ for some chamber γ. Recall the set of bases B(J A,γ) defined in (3.). With B(J A,γ), let R and R 2 be as in Theorem 2. Then (3.2) R (y,;z) R 2 (y;z) (3.22) where : R 2(;z) := k/ = R (y,;z) R 2 (y;z), [ ( z k z A A k ) νk, ], (3.23) R (y,;z) := z A y = { z A y z u z A A u u U δ (y A u )z u z A A u = u U ; ĥ (y) ĥ(a )u mod η with U := {u N n m 0 u k ν k, ; k / }. One can also obtain (3.2) by noticing that: R (y,;z) R (y,;z) = R 2(;z) R2 (;z) = ( +β ν k, + +β µ ν ) k,, k/ where β k, = z k z A A k, and µ is a multiple of ν k, Simplifications via finite number of generators. Decompose Z m into µ := deta disjoint equivalent classes, where y,ξ Z m are equivalent if and only if δ (y ξ) =. For every basis J A, let G be the },

10 0 JEAN B. LASSERRE AND EDUARDO S. ZERON quotient group defined in (3.4), that is, G = Z m/ A j Z = {Ec[0,],...,Ec[µ,]}. j Notice that y,ξ Z n belong to Ec[j,] if and only if A (y ξ) Z n, and that Z m is equal to the disjoint union of all classes Ec[j,]. Next, pick up a minimal representative element of every class, i.e., fix (3.24) ξ[j,] Ec[j,] such that A y A ξ[j,] 0, for every y Ec[j,] with A y 0. The minimal representative elements ξ[j,] in (3.24) can be computed as follows: Let d Ec[j,], arbitrary, and let d Z m be such that his k-entry d k is the smallest integer greater than or equal to the k-entry of A d. The vector ξ[j,] defined by d + A d satisfies (3.24). Notice that d + A d 0. Besides, let d,y Ec[j,] with A y 0. There exists q Z m such that y = d+a q. Hence q A d; in addition, q d follows from the above definition of d, and so A y d +A d 0. Therefore, the vector ξ[j,] := d+a d satisfies (3.24). In particular, if Ec[0,] is the class which contains the origin of Z m, then ξ[0,] = 0. Notice that for every integer vector y Z m, there exists a unique ξ[j,] such that : y = ξ[j,] + A q, for q Z m. Moreover, the extra condition A y 0 holds if and only if: (3.25) y = ξ[j,] + A q with q N m. We obtain a compact form of h(y;z) when y Z m γ, for some chamber γ. Theorem 4. Let h and ξ[j,] be as in (.2) and (3.24), respectively. Let y Z m γ, for some chamber γ. Recall the set of bases B(J A,γ) defined in (3.). For every basis B(,γ) there is a unique index 0 j[y,] < µ such that y is contained in the equivalence class Ec[j[y,],] defined in (3.4), and so: (3.26) h(y;z) = R (ξ[j[y,],],;z) z A y, R 2 (;z) B(,γ) where A y Z m is such that his k-entry is the largest integer less than or equal to the k-entry of A y. Proof. Recall that if y Z m γ h(y;z) = B(,γ) R (y,;z) R 2 (y;z) Next, recalling the definition (3.) of B(J A,γ), A y 0 for every basis B(J A,γ) with y γ. Recall that there is a unique index j[y,] < µ

11 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA such that y = ξ[j[y,],] +A q with q N m ; see (3.25) and the comment just before. To obtain the vector q N m, recall that the minimal representative element ξ[j[y,],] in (3.24) is the sum y +A y where y Z m is such that his k-entry yk is the smallest integer greater than or equal to A y, for we only need to fix d = y in the paragraph that follows (3.24). In particular, A y = y, and ξ[j[y,],] = y A A y, which when used in (3.7) and (3.2), yields, And so (3.6) implies (3.26). R (y,;z) = R (ξ[j[y,],],;z)z A y. Theorem 4 explicitly shows that it suffices to compute R (v,;z) for finitely many values v = ξ[j,], with B(,γ) and 0 j < µ, in order to calculate h(y;z) for arbitrary values y Z m γ, via (3.26). In other words, in theclosure γ of a chamber γ, one only needs to consider finitely many fixed convex cones C(Ω(ξ[j,]),) R n, where B(,γ) and 0 j < µ, and compute their associated rational function (3.26). The counting function h(y; z) is then obtained as follows. Input: y Z m γ, z C n. Output ρ = h(y;z). Set ρ := 0. For every B(,γ) : Compute ξ[j[y,],] := y A A y Z m. Read the value R (ξ[j[y,],],;z)/r 2 (;z), and update ρ by: ρ := ρ+ R (ξ[j[y,],],;z) R 2 (;z) z A y. For the whole space Z m it suffices to consider all chambers γ and all cones C(Ω(ξ[j,]),) R n, where B(,γ) and 0 j < µ. Finally, in view of (3.7)-(3.8), the above algorithm can be symbolic, i.e., z C m can be treated symbolically, and ρ becomes a rational fraction of z. 4. Generating function An appropriate tool for computing the exact value of h(y;z) in (.2) is the formal generating function H : C m C, (4.) s H(s) := n h(y;z)s y = z y Z m k s A, k k= where s y is defined in (2.4) and the sum is understood as a formal power series, so that we need not consider conditions for convergence. This generatingfunctionwasalreadyconsideredinbrionandvergne[8]withλ = ln s. Following notation of [8, p. 805], let 0 ˆx R n be a regular vector, i.e., no entry [A Aˆx] j vanishes for any basis J A or index j m.

12 2 JEAN B. LASSERRE AND EDUARDO S. ZERON Define : (4.2) ε j, := (4.3) { if [A Aˆx] j > 0, if [A Aˆx] j < 0. Next, for every basis J A, index j and vector u Z n m, fix : θ[j,,u ] Z : the smallest integer greater than or equal to ε j, [A A u ] j. Define also the vector η[,u ] Z n by : u j if j ; (4.4) η[,u ] j = θ[j,,u ] if j, ε j, = ; θ[j,,u ] if j, ε j, =. The following expansion can be deduced from [8]. Theorem 5. Let 0 ˆx R n be regular and consider the vectors η[,u ] Z n defined in (4.4) for J A and u Z n m. The following expansion holds: n z k s A = [ ] (4.5) k z j s A j k= R 2 (;z) J A u Z n m µ j z η[,u ] s Aη[,u ], where Z µ = {0,,...,µ }, µ = deta and: (4.6) z R 2 (;z) := [ ( z k z A k/ A k ) µ ]. Proof. From Brion and Vergne s identity [8, p. 83], n (4.7) e w = [ ε j, ]F(C ˆx +ρ (C),L), k j= J A j where F(C ˆx + ρ (C),L) is the formal power series l el added over all elements l in the intersection of the cone C ˆx +ρ (C) with the integer lattice L = Z[w,...,w n ]. Moreover, the coefficients ε j, are defined in (4.2) and the cone Cˆx is defined by the following formula [8, p. 805], { } (4.8) C ˆx = ε j, x j w j x R m, x 0. j Finally, given the real vector space W = R[w,...,w n ], every ρ : W W is a linear mapping defined by its action on each basis element w k of W, (4.9) ρ (w k ) := w k j [A A k ] j w j.

13 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA 3 Hence, ρ (w j ) = 0 for every j, and the cone ρ (C) is given by { (4.0) ρ (C) = k w k k x [A x R A x ] j w n m }, j ; x 0 j see [8, p.805]. Thus, every element in theintersection of the cone C ˆx +ρ (C) with the lattice Z[w,...,w n ] must be of the form : (4.) k w k + k x ε j, ξ j w j, with x N n m, j ξ Z m and ξ j ε j, [A A x ] j. On the other hand, for every basis, define µ = deta and : (4.2) x = u +µ v, with u Z n m µ and v N n m. Moreover, as in (4.3), fix θ[j,,u ] Z to be the smallest integer greater than or equal to ε j, [A A u ] j. Thus, we can rewrite (4.) so that the intersection of the cone C ˆx +ρ (C) with the lattice Z[w,...,w n ] must be of the form : [ uk w k +v k µ ρ(w k ) ] + [ ] (4.3) ε j, w j θ[j,,u ] j +q j, k j with u Z n m µ, v N n m and q N m. Wecandeduce(4.3)from(4.)byrecallingthedefinition(4.9)ofρ (w k ) and letting : ξ j := θ[j,,u ] + q j ε j, µ [A A v ] j. Since F(C ˆx +ρ (C),L) is the formal power series l el with summation over all elements l in (4.3), one obtains (4.4) F(C ˆx +ρ (C),L) = [ u Z n m µ j e ε j,θ[j,,u ]w j e ε j,w j ][ k ] e u kw k e µρ(w. k) With η[,u ] Z n being as in (4.4), using (4.4) into (4.7) yields the expansion n e w = [ ] (4.5) k e w j j= [ n J A u Z n m µ j=e η[,u ]jwj][ k j e µρ(w k) Finally, we defined w k := ln(z k ) + ln s A k for every index k n, where the vectors s,z C n have all their entries different from zero and ln s is the [ n] matrix defined in (2.4). So e w k = z k s A k. Moreover, ].

14 4 JEAN B. LASSERRE AND EDUARDO S. ZERON recalling the definition (4.9) of ρ (w k ), the following identities hold for all k n, (4.6) ρ (w k ) = ln(z k ) j ln(z j )[A A k ] j. Notice j A j[a A k] j = A k. A direct application of (4.6) and the identities e w k = z k s A k yields (4.5), i.e.: n k= z k s A = k J A u Z n m µ with R 2 (;z) = k/ z η[,u ] s Aη[,u ] R 2 (;z) [ ( z k z A A direct expansion of (4.5) yields the following: A k j z j s A j, ) µ ]. Theorem 6. Let 0 ˆx R n be regular, and let h and η be as in (.2) and (4.4), respectively. Let J A be the set of bases associated with A. Then for every pair of (y,z) Z m C n with z < : (4.7) h(y;z) = J A,A y 0 z A y R 2 (;z) u Z n m µ z u z A A u { ( ) if A y Aη[,u] N m, 0 otherwise, where: Z µ = {0,,...,µ }, µ = deta, (4.8) (4.9) 0 [ A Aη[,u] ] for each j j and R 2 (;z) := [ ( ) ] z k z A A µ k. k/ Proof. Recall the expansion of H(s) as a formal power series : (4.20) H(s) = n h(y;z)s y = z y Z m k s A. k k= We also have a similar formal power series for the product : z j j s A = [ ] z q j j j sa jq j = j q j = Combining the latter with (4.5) yields: (4.2) H(s) = J A q N m u Z n m µ z q saq. q N m z η[,u ] z q s Aη[,u ] s Aq. R 2 (;z)

15 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA 5 Notice that (4.20) and (4.2) are identical. Hence, if we want to obtain the exact value of h(y;z) from (4.2), we only have to sum up all terms with exponent Aη[,u ] + A q equal to y. That is, recalling that each A is invertible, (4.22) h(y;z) = J A u Z n m µ z η[,u ] z q R 2 (;z) { ( if A y Aη[,u ] ) N m, 0 otherwise. On the other hand, setting q := A (y Aη[,u ]) and recalling the definition (4.4), zη[,u ] (4.23) z η[,u ] z q z u = za y z A Aη[,u ] = za y z u z A A u We finally prove that the vector A Aη[,u ] is bounded, so that (4.8) holds, (4.24) 0 [ A Aη[,u ] ] j From (4.3) and (4.4), the following identity (4.25) for each j. [ A Aη[,u ] ] j = [ A A u ] j + { θ[j,,u ] if ε j, = ; θ[j,,u ] if ε j, = holds for every index j. Next, suppose that [A A u ] j = R+r where R Z and 0 r < are the respective integer and fractional parts. We can obtain θ[j,,u ] in (4.3) as follows: If ε j, = then θ[j,,u ] = R is the smallest integer greater than or equal to R r, and so (4.25) yields [ A Aη[,u ] ] j = r. If ε j, = and 0 < r <, then θ[j,,u ] = R+ and : [ A Aη[,u ] ] j = r. At last, if ε j, = and r = 0, then θ[j,,u ] = R and so, [ A Aη[,u ] ] j =. Therefore (4.24) holds because 0 r <. Finally, a direct application of (4.23) into (4.22) yields the following version of (4.7):. h(y;z) = J A,A y 0 R 2 (;z) u Z n m µ z A y z u z A A u { ( if A y Aη[,u ] ) N m, 0 otherwise.

16 6 JEAN B. LASSERRE AND EDUARDO S. ZERON Notice that the first sum is calculated only on those bases J A which satisfy A y 0. This follows from (4.22) combined with (4.24). Observe that (4.7) is different from (3.6) or (3.2) because in (3.6) and (3.2), the summation is over bases in the subset B(,γ) {J A ; A y 0}. References [] Baldoni-Silva, W., Vergne, M.: Residues formulae for volumes and Ehrhart polynomials of convex polytopes. arxiv:math.co/ v, 200. [2] Baldoni-Silva, W., De Loera, J.A., Vergne, M.: Counting integer flows in networks. Found. Comput. Math. 4 (2004), [3] Barvinok, A.I.: Computing the volume, counting integral points and exponentials sums. Discr. Comp. Geom. 0 (993), [4] Barvinok, A.I., Pommersheim, J.E.: An algorithmic theory of lattice points in polyhedral. in: New Perspectives in Algebraic Combinatorics, MSRI Publications 38 (999), [5] Beck, M.: Counting Lattice Points by means of the Residue Theorem. Ramanujan Journal 4 (2000), [6] Beck, M., Diaz, R., Robins, S.: The Frobenius problem, rational polytopes, and Fourier-Dedekind sums. J. Numb. Theor. 96 (2002), 2. [7] Beck, M., Pixton, D.: The Ehrhart polynomial of the Birkhoff polytope. Discr. Comp. Math. 30 (2003), [8] Brion, M., Vergne, M.: Residue formulae, vector partition functions and lattice points in rational polytopes. J. Amer. Math. Soc. 0 (997), [9] Cochet, C.: Réduction des graphes de Goretsky-Kottwitz-MacPherson; nombres de Kostka et coefficients de Littlewodd-Richardson. Thèse de Doctorat: Mathématiques, Université Paris 7, Paris, Décembre [0] De Loera, J.A., R. Hemmecke, R., Tauzer, J., Yoshida, R.: Effective lattice point counting in rational convex polytopes. J. of Symb. Comp., to appear. [] Pukhlikov, A.V., Khovanskii, A.G.: A Riemann-Roch theorem for integrals and sums of quasipolynomials over virtual polytopes. St. Petersburg Math. J. 4 (993), [2] Lasserre, J.B., E.S. Zeron, E.S.: On counting integral points in a convex rational polytope. Math. Oper. Res. 28 (2003), [3] Lasserre, J.B., Zeron, E.S.: An alternative algorithm for counting lattice points in a convex polytope. Math. Oper. Res. 30 (2005), [4] Schrijver, A.: Theory of Linear and Integer Programming. John Wiley & Sons, Chichester, 986. [5] Szenes, A., Vergne, M.: Residue formulae for vector partitions and Euler-MacLaurin sums. Adv. in Appl. Math.30 (2003), [6] Szenes, A.: Residue theorem for rational trigonometric sums and Verlinde s formula. Duke Math. J. 8 (2003), [7] Verdoolaege, S., Beyls, K., Bruynooghe, M., Seghir, R., Loechner, V.: Analytical Computation of Ehrhart Polynomials and its Applications for Embedded Systems. Technical report # 376, Computer Science Department, KUL University, Leuwen, Belgium.

17 EXPLICIT FORMULA FOR COUNTING LATTICE POINTS OF POLYHEDRA 7 LAAS-CNRS and Institute of Mathematics, LAAS 7 Avenue du Colonel Roche, 3077 Toulouse Cédex 4, France. address: lasserre@laas.fr Depto. Matemáticas, CIVESTAV-IPN, Apdo. Postal 4740, Mexico D.F , México. address: eszeron@math.cinvestav.mx