Abstract 1 Introduction 2 Main results References

Integer Points in Dilates of Polytopes

Shubhangi Saraf ORCID University of Toronto, Canada    Narmada Varadarajan ORCID University of Toronto, Canada
Abstract

In this paper we study how the number of integer points in a polytope grows as we dilate the polytope. We prove new and essentially tight bounds on this quantity by specifically studying dilates of the Hadamard polytope.

The motivation for studying this question comes from the question of understanding the maximal number of monomials in a factor of a multivariate polynomial of s monomials. A recent result by Bhargava, Saraf and Volkovich [4] showed that if f is an n-variate polynomial, where each variable has degree d, and f has s monomials, then any factor of f has at most sO(d2logn) monomials. The key technical ingredient of their proof was to show that any polytope with s vertices, where each vertex lies in {0,..,d}n can have at most sO(d2logn) integer points. The precise dependence on d of the number of integer points was left open. We show that this bound, and in particular the dependence on d is essentially tight by studying dilates of the Hadamard polytope and proving new lower bounds on its number of integer points.

Keywords and phrases:
Computational geometry, Complexity theory, Integer polytopes
Copyright and License:
[Uncaptioned image] © Shubhangi Saraf and Narmada Varadarajan; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation
Related Version:
Full Version: https://arxiv.org/pdf/2510.16481
Acknowledgements:
We thank the reviewers for their helpful suggestions.
Funding:
Research partially supported by the McLean Award and an NSERC Discovery Grant.
Editors:
Hee-Kap Ahn, Michael Hoffmann, and Amir Nayyeri

1 Introduction

Integer polytopes play a starring role in several areas of mathematics, and have recently found very interesting applications in theoretical computer science. In particular, the problem of counting integer points in polytopes lies at the intersection of combinatorics and convex geometry, with applications in algebraic geometry [7, 8], number theory [3, 6], optimization [2, 1], cryptography [10], and algebraic complexity theory [4], just to name a few. Perhaps the most celebrated result on this topic is the Ehrhart polynomial: given an integer polytope Pn, the number of integer points in its dilate, |dPn|, is a polynomial in d [9]. When P is an n-dimensional polytope, the leading term of its Ehrhart polynomial is vol(P)dn, which is a good approximation to the number of integer points for large values of d. However, the coefficients of the Ehrhart polynomial are notoriously difficult to compute (see, for example, [11] for examples of polytopes with negative Ehrhart coefficients), so we have little information about integer points when the dilation factor d is much smaller than the dimension n.

Motivation from algebraic complexity theory

Our interest in estimating |dPn| when d is smaller than n comes from the problem of sparse polynomial factorization in complexity theory. A multivariate polynomial f𝔽[X1,,Xn] where each variable Xi appears with degree d is sparse if it only has poly(n) monomials. Polynomial factorization is a fundamental question in computational algebra, and sparse polynomial factorization is a special case that has been studied for several decades now [4, 5, 13]. The first randomized algorithm for sparse polynomial factorization exhibited a dependence on the sparsity of the factors of the polynomial [13], which raised the natural question: do factors of sparse polynomials have to be sparse?

We know that this is not true in general: if g is a factor of some polynomial f, there is no polynomial upper bound for the number of monomials of g in terms of the number of monomials of f. In characteristic 0, we have the counterexample f=i=1n(XidYid) and g=i=1n(Xid1+Xid2Yi++Yid1). Here, f has 2n monomials, but g is a factor of f with dn monomials. In characteristic p, this gap can be more stark. Consider f=i=1nXip and g=(i=1nXi)p/2, so f has n monomials, but g is a factor of f with at least np/2 monomials. In general, we do not have know better than the trivial upper bound of dn for the number of monomials in a factor of a sparse polynomial. However, note that in the setting where we assume that the individual degree d is an absolute constant and n can be arbitrarily large, these two counterexamples no longer exhibit significant sparsity gaps.

The current best-known bound for sparse factorization in this setting is due to [4], who show that if f has s monomials and g is a factor of f, then g has at most sO(d2logn) monomials. The key technique of their proof is to reduce the problem of counting monomials in f and g into one of counting integer points in their Newton polytopes. In this paper we further study the problem of counting integer points in polytopes, since a better upper bound there would directly imply a better bound for the factoring problem. Unfortunately, our results indicate that there is no better upper bound for polytopes, so any attempt to improve the bound for sparsity will require a new approach.

Connections to studying dilates of polytopes

The key technique to proving the upper bound of sO(d2logn) in [4] is to study integer points in polytopes. In particular if f=g×h, then one can show that the Newton polytope of f is the Minkowski sum of the Newton polytopes of g and h. The high-level idea in the proof is to show the following: If g has too many monomials, then the Newton polytope of g must necessarily have many vertices. This can then be used to show that the Newton polytope of f must also have many vertices, hence f must have many monomials. Their main technical contribution is to show that any polytope with s vertices, where each vertex lies in [d,+d]n, can have at most sO(d2logn) integer points. The precise dependence on d of the number of integer points was left open.

The main goal of our paper is to understand if this bound for integer points in polytopes is tight. Saptharishi observed that the logn in the exponent of the bound is necessary (see [4] for a proof of this) by counting integer points in Hadamard polytopes, a class of simplices with vertices in {1,1}n. Specifically, he shows that the number of integer points in the Hadamard polytope is at least nΩ(logn). In this paper, we show that the dependence of d2 in the exponent is also essentially tight by studying dilates of the Hadamard polytope.

An astute reader may recognize the Hadamard polytope as the extremal solution to Hadamard’s maximal determinant problem: the parallelepiped spanned by its vertices maximizes volume among all parallelepipeds with vertices in {1,1}n. We expect that the Hadamard polytope should also maximize the number of integer points – a discrete analog of volume – among all simplices with vertices in {1,1}n. This is the heuristic that motivates us to study the Hadamard polytope as the “worst-case simplex” for our problem.

In the case when s=poly(n) (as with factors of sparse polynomials), we can reduce our problem about counting integer points in polytopes in [d,d]n to counting integer points in dilates of polytopes as follows. First, by Carathéodory’s theorem, each vertex of P is a convex combination of at most n+1 points of {d,d}n. Replacing the vertices of P with these points, we obtain a polytope Pd containing P that still has poly(n) vertices, but they now lie in {d,d}n. So, P0=1dPd is an integer polytope with vertices in {1,1}n and |dP0n||Pn|.

Our contribution in this paper is to provide lower bounds for the number of integer points in the Hadamard polytope. Our first main result is an elegant algebraic characterization of the integer points in the Hadamard polytope by constructing an explicit bijection with affine subspaces of 𝔽2logn. This underpins part of our second main result: a lower bound for the number of integer points in small dilates. To lower bound the points in larger dilates, we use the probabilistic method. These bounds show that both the dependences on d2 and on logn in the bound |dPn||V(P)|O(d2logn) are necessary.

2 Main results

Let n=2m for some m. The n×n Hadamard matrix H has its rows and columns indexed by elements of the vector space 𝔽2m:

H(a,b) =(1)a,b,

where a,b is the standard dot product on 𝔽2m.

Definition 1.

The Hadamard polytope PHadn is the convex hull of the column vectors of H.

The columns of H are linearly independent so they form the vertices of PHad; the Hadamard polytope is an (n1)–simplex with vertices in {1,1}n. The lower bound |PHadn|nΩ(logn) was first observed by Ramprasad Saptharishi and appears in the paper [4]. Our contribution is a matching upper bound on |PHadn| through a complete characterization of all integer points in the Hadamard polytope using the vector space structure of F2m. This characterization will also be useful when we study integer points in dilates of PHad.

First, some notation: let {ha:a𝔽2m} be the set of column vectors of the Hadamard matrix H, i.e. the vertices of the Hadamard polytope PHad. We can also identify the coordinates of n with vectors in 𝔽2m, with the convention that the first coordinate in n corresponds to the zero vector in 𝔽2m. For every a𝔽2m and vn, let v(a) denote the ath coordinate of v in the standard basis, and

supp(v)={a𝔽2m:v(a)0}

We also know that every point vPHad can be expressed as a convex combination v=a𝔽2mtaha. Since the vectors {ha:a𝔽2m} are linearly independent, this expression is unique, and we can define

T(v) ={a𝔽2m:ta0}.

This can also be thought of as the support of v in the basis {ha:a𝔽2m}.

Statements of main results

Theorem 2.

For every vPHadn, supp(v) is a subspace of 𝔽2m. Further, T(v)=supp(v)+b for some b𝔽2m, and

v =1|T(v)|cT(v)hc

is actually a uniform linear combination of vertices.

Notice that this theorem is a combinatorial correspondence between affine subspaces of 𝔽2m and integer points in PHad. That is, every affine subspace S𝔽2m, let vS be the uniform convex combination of the vertices of PHad corresponding to the points of S. In other words, vS is the barycenter of the corresponding face of PHad. Then, the map SvS is a bijection between affine subspaces of 𝔽2m and integer points in PHad.

Corollary 3.

The number of integer points in PHad is |PHadn|=nΘ(logn).

Theorem 4.

For any d<nlogn and fixed 0<ε<1/2, the number of integer points in the dilate dPHad is lower bounded by

|dPHadn|{nΩ(dlogn),dlogn4;nΩ(log2n),logn4dlog1.5nn2εd2/logn,log1.5nd(nlogn)1/2ε;(d2ε/4)n,(nlogn)1/2+εd<nlogn.

Proofs of main results

Proof of Theorem 2.

First, we need to show that supp(v) is a subspace of 𝔽2m. We will actually show something stronger: if v(a)0 and v(b)0, then v(a+b)=v(a)v(b). Note that since v is an integer point in PHad[1,1]n, its coordinates actually lie in {1,0,1}n. So, the three cases we need to consider for our proof are: v(a)=v(b)=1; v(a)=v(b)=1; and v(a)=1, v(b)=1. We will write out the details for the case when v(a)=v(b)=1, and the other cases will follow a similar computation.

Since v is a convex combination of the vertices of PHad, we can write it as v=c𝔽2mtchc, for some nonnegative real numbers (tc)c𝔽2m that satisfy c𝔽2mtc=1. From the definition of the Hadamard matrix H, we get the following identities for the coordinates a,b𝔽2m,

v(a) =c𝔽2mtchc(a)=c𝔽2m(1)c,atc=catcc⟂̸atc,
v(b) =c𝔽2mtchc(b)=c𝔽2m(1)c,btc=cbtcc⟂̸btc.

Define

sab=ca,cbtc;sa=ca,c⟂̸btc;sb=c⟂̸a,cbtc;s0=c⟂̸a,c⟂̸btc.

This allows us to rewrite our identity for v(a) as

v(a) =sab+sasbs0.

A similar identity for v(b) and the identity c𝔽2mtc=1 gives us the system of linear equations

sab+sasbs0 =v(a)=1
sabsa+sbs0 =v(b)=1
sab+sa+sb+s0 =1

The general solution to this system has the form

sabs0 =1
sa+sb =0
sb+s0 =0.

Since each of sab,sa,sb, and s0 must be nonnegative, this yields a unique solution sab=1,sa=sb=s0=0. So,

v(a+b) =c𝔽2m(1)a+b,ctc
=ca,cbtcca,c⟂̸btcc⟂̸a,cbtc+c⟂̸a,c⟂̸btc
=sabsasb+s0
=1=v(a)v(b).

The argument for the other cases follows similarly: we set up a system of equations in the variables sab,sb,sa,s0 that has a unique solution subject to nonnegativity. Computing the value of v(a+b) for this unique solution shows that v(a+b)=v(a)v(b). This shows not only that supp(v) is a basis, but also that the entries of v are determined by their values on a basis of supp(v).

Now, let b1,,bk𝔽2m be a basis for supp(v), and choose εi{0,1} so that v(bi)=(1)εi for each i=1,,k. Consider the linear map 𝔽2msupp(v) that sends c𝔽2m to b1,cb1++bk,cbksupp(v). The kernel of this map is exactly supp(v), so its range must be all of supp(v). In particular, there is some b𝔽2m whose image under this map is ε1b1++εkbk. So, v(bi)=(1)bi,b for each i=1,,k. For any asupp(v), we know that a=iIbi for some I[k], so

v(a)=iIv(bi)=(1)iIbi,b=(1)a,b.

We claim that T(v)=supp(v)+b for this vector b.

Define the vector v𝔽2m by

v =1|supp(v)|csupp(v)+bhc.

We will show that v=v by showing that v(a)=v(a) for every coordinate a𝔽2m.

First, for each asupp(v) and csupp(v)+b, we have that a,c=a,b. So,

v(a) =1|supp(v)|csupp(v)+bhc(a)
=1|supp(v)|csupp(v)+b(1)a,c
=(1)a,b
=v(a).

Next, for each asupp(v), the 𝔽2-linear map from supp(v){0,1} that sends csupp(v) to a,c is not identically zero. So, the size of the preimage of 0 under this map is equal to the size of the preimage of 1. This property still holds when we consider the translate supp(v)+b: the dot product a,c is 0 for exactly half the vectors csupp(v)+b, and is 1 for the other half. Thus, the average of hc(a) over all csupp(v)+b is zero, so

v(a)=1|supp(v)|csupp(v)+b(1)a,c=0=v(a),

since asupp(v). This shows that v does indeed have the form we claimed, with T(v)=supp(v)+b and

v=1|T(v)|cT(v)hc

Proof of Corollary 3.

Note again that the lower bound |PHadn|=nΩ(logn) was already known by Saptharishi’s contribution to [4]; we will provide a matching upper bound of nO(logn). By Theorem 2, integer points of PHadn are in one-to-one correspondence with affine subspaces of 𝔽2m. Any k-dimensional subspace V𝔽2m has exactly 2mk corresponding affine subspaces, and there are exactly [mk]2 (the Gaussian binomial coefficient) subspaces of dimension k. So, using the upper bound [mk]224m2 (see, for example, [12] for a proof of this bound),

|PHadn|=k=0m2mk[mk]2(m+1)2m+4m2=nO(logn),

which completes the proof.

Now we are ready to prove our main theorem: lower bounding the growth of integer points in the dilates dPHad.

Proof of Theorem 4.

Case 1.

d(logn)/4.

Let us introduce some more notation for our counting argument. We will show that there are nΩ(dlogn) distinct sets of d integer points in PHad whose sums are all distinct in dPHad. For vPHadn, define the exponent of v as exp(v)=dim(supp(v)), so |T(v)|=2exp(v), and

v =12exp(v)aT(v)ha.

Let ta(v) be 2exp(v) if aT(v), and 0 otherwise; that is, ta(v) is the coefficient of the vertex ha in the convex combination for v. Look at all sets {v1,,vd} such that

  1. (i)

    T(vi) is a subspace, i.e. T(vi)=supp(vi);

  2. (ii)

    logn2d<exp(vi)<logn2+d;

  3. (iii)

    exp(vi)exp(vj) if ij.

Since dlogn4, for n large enough, the number of subspaces of 𝔽2logn with dimension k for some k(logn2d,logn2+d) is [lognk]22k(lognk)=nΩ(logn). So, the number of choices for the set of subspaces {T(v1),,T(vd)} with distinct exponents in (logn2d,logn2+d) is nΩ(dlogn). We will now show that each set {v1,,vd} satisfying these properties produces a distinct point in dPHadn, proving the lower bound |dPHadn|nΩ(dlogn).

Now, suppose {v1,,vd} and {u1,,ud} are two distinct sets satisfying the above conditions. And suppose for contradiction that ivi=iui. Expanding this in terms of the vertices of PHad, we get

i=1da𝔽2mta(vi)ha =i=1da𝔽2mta(ui)ha

Since the vertices of PHad are linearly independent, we must have the following stronger identity for all a𝔽2m:

i=1dta(vi)=i=1dta(ui).

For each fixed a, since ta(vi) is either zero or 2exp(vi) and the exponents are all distinct, we can think of each sum as the unique binary representation of some dyadic rational number.

When a=0, we get t0(vi)=2exp(vi) and t0(ui)=2exp(ui) since T(vi) and T(ui) are linear subspaces containing 0, for every i=1,,d. This gives us the identity

i=1d2exp(vi)=i=1d2exp(ui).

So, the sets {exp(vi):id} and {exp(ui):id} must be equal since the above sums represent the same binary number. This means that the only way for the sums i=1dvi and i=1dui to be equal at the 0th coordinate is if their exponents are all the same. Now we will find another coordinate where they differ.

Assume without loss of generality that exp(vi)=exp(ui) for all i=1,,d, and exp(v1)<<exp(vd). There must be some k such that T(vk)T(uk). Choose aT(vk)T(uk) (which is possible since they are different subspaces of the same dimension), and now look at the identity

i=1dta(vi)=i=1dta(ui).

Again, each sum is the unique binary representation of some dyadic rational number, but ta(vk)0 while ta(uk)=0, so they cannot represent the same number, which is the desired contradiction.

Case 2.

(logn)/4dlog1.5n .

We simply reuse the bound from case 1 (note that the lower bound is monotone in d and hence we can just use the bound in the setting d=logn4),

|dPHadn||logn4PHadn|nΩ(log2n).
Case 3.

dlog1.5n.

For this final case we will use the probabilistic method. Take the projection nn1 that deletes the first coordinate of every vector, and let P0 be the image of PHad under this map. Recall that PHadn is an (n1)-simplex whose vertices all have first coordinate equal to 1, since the first coordinate corresponds to the zero vector in 𝔽2m. So, this projection is a bijection on PHad. In particular, integer points of dP0 are in one-to-one correspondence with integer points of dPHad, so |dPHadn|=|dP0n1|.

The new polytope P0 is now a full-dimensional simplex in [1,1]n1. Since the vertices of the original Hadamard polytope PHadn were orthogonal, for each vertex vV(P0)n1, the face of P0 that does not contain v is defined by the hyperplane {xn1:x,v=1}. In particular, we can write the simplex P0 as an intersection of half-spaces:

P0 =vV(P0){xn1:x,v1},

and

dP0 =vV(P0){xn1:x,vd}.

Fix some set of c coordinates S[n] and some D[d] (with c and D to be determined later) and look at the discrete hypercube HS[D] supported on those coordinates,

HS[D]=([D,D]{0})S×{0}[n]S.

We will use the probabilistic method to show that P0 contains at least 1/2 the points of HS[D] for every choice of S. Note that the hypercubes HS[D] are disjoint by construction, so

|dP0n|S[n],|S|=c|dP0HS[D]|.

Fix a vertex vV(P0), and a hypercube HS[D]. Let X be a uniform random variable taking values in HS[D]. Since v is a ±1-valued vector, 𝔼|X,v|=0. Hoeffding’s inequality implies that

(|X,v|>d) 2ed2/2cD.

By the union bound,

(XdP0)vV(P0)(|X,v|>d)2ned2/2cD.

Assuming that 2cDd22logn, the above probability is at most 1/2, so at least 1/2 the points of HS[D] must be contained in dP0. Now, summing over the disjoint sets HS[D] for all (nc) choices of S,

|dPHadn|=|dP0n1|S[n],|S|=c|HS[D]dP0|12(nc)(2D)c.

We divide our analysis into two subcases.

Case 3a.

d(nlogn)12ε, which implies that d2lognn12ε

In this case, set D=1 and c=d24logn, so

Case 3b.

d(nlogn)12+ε.

For our probabilistic argument to work, we need 2cDd22logn, which rearranges to Dd24clogn. Set D=d2ε/4 and c=n. Now,

as desired. Now,

|dPHadn| 12(nc)(2D)c(d2ε/4)n.

References

  • [1] Divesh Aggarwal, Antoine Joux, Miklos Santha, and Karol Węgrzycki. Polynomial time algorithms for integer programming and unbounded subset sum in the total regime, 2024. Submitted July 7, 2024; revised July 11, 2024. doi:10.48550/arXiv.2407.05435.
  • [2] Iskander Aliev and Martin Henk. Feasibility of integer knapsacks. SIAM Journal on Optimization, 20(6):2978–2993, 2010. doi:10.1137/090778043.
  • [3] Matthias Beck and Sinai Robins. Computing the Continuous Discretely. Undergraduate Texts in Mathematics. Springer, New York, NY, 1 edition, 2007. doi:10.1007/978-0-387-46112-0.
  • [4] Vishwas Bhargava, Shubhangi Saraf, and Ilya Volkovich. Deterministic factorization of sparse polynomials with bounded individual degree. J. ACM, 67(2), May 2020. doi:10.1145/3365667.
  • [5] P. Bisht and I. Volkovich. On solving sparse polynomial factorization related problems. Computational Complexity, 34(7), 2025. doi:10.1007/s00037-025-00268-5.
  • [6] Sheng Chen, Nan Li, and Steven V. Sam. Generalized ehrhart polynomials. Transactions of the American Mathematical Society, 364(1):551–569, 2012. doi:10.1090/S0002-9947-2011-05494-2.
  • [7] David A. Cox. Recent developments in toric geometry, 1996. arXiv:alg-geom/9606016.
  • [8] Luis Crespo, Álvaro Pelayo, and Francisco Santos. Ewald’s conjecture and integer points in algebraic and symplectic toric geometry, 2024. arXiv:2310.10366.
  • [9] Eugène Ehrhart. Sur les polyèdres rationnels homothétiques à n dimensions. Comptes rendus de l’Académie des Sciences, 254:616–618, 1962.
  • [10] Yansong Feng, Hengyi Luo, Qiyuan Chen, Abderrahmane Nitaj, and Yanbin Pan. Computing asymptotic bounds for small roots in coppersmith’s method via sumset theory. In Yael Tauman Kalai and Seny F. Kamara, editors, Advances in Cryptology – CRYPTO 2025, pages 3–32, Cham, 2025. Springer Nature Switzerland. doi:10.1007/978-3-032-01855-7_1.
  • [11] Takayuki Hibi, Akihiro Higashitani, Akiyoshi Tsuchiya, and Koutarou Yoshida. Ehrhart polynomials with negative coefficients. Graphs and Combinatorics, 35(1):363–371, 2019. doi:10.1007/S00373-018-1990-9.
  • [12] Peter M. Neumann and Cheryl E. Praeger. Cyclic matrices over finite fields. Journal of the London Mathematical Society, Series 2, 52(2):263–284, 1995. doi:10.1112/jlms/52.2.263.
  • [13] Joachim von zur Gathen and Erich Kaltofen. Factoring sparse multivariate polynomials. Journal of Computer and System Sciences, 31(2):265–287, 1985. doi:10.1016/0022-0000(85)90044-3.