Skip to content

定理Theorem

Ehrhart 伸缩计数定理

Ehrhart theorem · Ehrhart polynomial · Ehrhart quasipolynomial · Ehrhart多项式 · Ehrhart准多项式

证明有理多胞形的全部非负整数伸缩格点数按剩余类呈多项式,整数顶点时恢复单一多项式;以同高整圆锥和半开基本域给生成函数证书,处理周期折叠、低维切片和负系数。

把一张图形放大两倍,面积乘四,格点数却不恰好乘四:边界和原点也有贡献。若顶点坐标是有理数,边界还会周期性地经过格点。Ehrhart定理说明,这些误差不是任意波动,而由有限组多项式精确控制。

形式陈述 ​

输入、有理分母与计数函数 ​

设 P⊂Rd为非空有理多胞形,仿射维数为 r。有理表示全部顶点坐标属于 Q;多胞形有界,故每个伸缩只含有限多个标准格点。

取正整数 q使 qP的顶点全为整数。最小这样的 q称为 P的分母,等于所有顶点坐标既约分母的最小公倍数。对整数 t≥0定义

LP(t)=|tP∩Zd|,tP={tx:x∈P}.

Ehrhart定理断言:存在 q个有理系数多项式 f0,…,fq−1,使

LP(qk+s)=fs(k)(k≥0, 0≤s<q).

各 fs的次数至多为 r,其中最大次数恰为 r。这样的函数称为准多项式;等价地,LP(t)=cr(t)tr+⋯+c0(t),各系数 cj都是整数变量上的周期函数。

q是一个可用周期,最小周期可能更小;精确地说,最小周期整除 q。若顶点全为整数,可取 q=1,于是 LP是次数 r的一个普通多项式。由于 0P={0},总有 LP(0)=1。这里结论从零开始成立,不只是“大到一定程度后成立”。

若 r=d,每个剩余类以 t为变量时的最高项系数都等于普通体积 vold(P)。低维有理多胞形可能在某些剩余类完全没有格点,不能把这一句不加修改地照搬。

一个有限的生成函数证书 ​

证明还给出普通生成函数

(1)EP(z)=∑t≥0LP(t)zt=H(z)(1−zq)r+1,H(z)=∑h=0q(r+1)−1ahzh,

其中 ah为非负整数。分子由若干半开基本域的有限格点计数得到,不需要猜测多项式。

由 (1−w)−r−1=∑j≥0(j+rr)wj抽系数,可直接恢复

(2)fs(k)=∑j=0ras+qj(k−j+rr).

这里二项式写成次数 r的多项式。若 0≤k<j≤r,其上参数在 0,…,r−1之间,乘积本身为零;因此公式连最初几项也正确。

直觉

多放一个坐标,把放大次数变成高度 ​

令 vi为顶点,把它们改成整向量

ui=(qvi,q)∈Zd+1.

它们生成的闭圆锥 C在高度 t>0的截面恰为 tP×{t},高度零只有原点。于是数所有伸缩图形,等于逐层数同一个圆锥。

若一个小锥由 r+1条线性无关整向量生成,将每个坐标拆为非负整数部分和基本区间中的剩余部分,就把无限计数缩成有限计数。每沿任一生成元走一步,高度增加 q,因此出现 r+1个因子 (1−zq)−1。

切块时,边界必须只归一块 ​

闭三角形可以共用一条边。若直接把各块的闭点数相加,这条边会重复;若把每块全部边界删掉,共用边又会漏掉。正确办法是选一个统一的微小移动方向,让共享边界上的每一点只归给移动后进入的那一块。

这产生的块通常是半开的:有些面保留,有些面删除。半开块是计数的分工方式,不等于原多胞形的内部;下一页的互反正要利用两种不同边界规则之间的关系。

例子与边界

六个剩余类,而非一个普通多项式 ​

取

P={(x,y):x,y≥0, 2x+3y≤1}.

顶点为 (0,0),(1/2,0),(0,1/3),分母六。对固定 t,加入整数松弛量 w=t−2x−3y≥0,就有

EP(z)=1(1−z)(1−z2)(1−z3).

统一分母得

H(z)=(1+z+⋯+z5)(1+z2+z4)(1+z3),EP(z)=H(z)(1−z6)3.

将 t=6k+s代入式(2),或直接计算

LP(t)=∑y=0⌊t/3⌋(⌊t−3y2⌋+1),

得到

(3)LP(6k+s)=3k2+(s+3)k+cs,(c0,…,c5)=(1,1,2,3,4,5).

例如 LP(6)=7,LP(7)=8,LP(8)=10。以原变量 t写时,二次项都是 t2/12,一次项都是 t/2,常数项依次为

1,512,23,34,23,512.

这个周期常数只有每六步一次取值一,故最小周期确为六;面积为 1/12也与首项吻合。

分母二,也可能完全没有周期振荡 ​

取 Q=conv{(0,0),(2,0),(0,1/2)}。它的分母为二,但

LQ(t)=∑y=0⌊t/2⌋(2t−4y+1)=(t+1)(t+2)2.

若 t=2k,和为 (k+1)(2k+1);若 t=2k+1,和为 (k+1)(2k+3),两者都化为上面同一个多项式。这称为周期折叠,说明定理只能保证最小周期整除分母。

低维有理切片会空掉整个剩余类 ​

在 R2中取竖线段

S={(1/2,y):0≤y≤1}.

它的仿射维数为一,然而

LS(t)={t+1,t 偶,0,t 奇.

奇数伸缩时,整条线仍在半整数横坐标上。这不违背次数为一的结论:它指各剩余类次数的最大值。若改成零维点 {1/2}⊂R,计数则为偶数一、奇数零,次数为零。

有理性、整数参数与有界性各有用途 ​

对无理区间 [0,2],计数为 ⌊t2⌋+1,不可能有有限准周期。否则在 t=qk这一列上,它因线性增长只能是一阶多项式;连续整数 k处的差是整数,斜率必须为整数,但增长极限给斜率 q2,矛盾。

即使顶点为整数,也不能把多项式在任意实数 t处的值当成点数。例如 P=[0,1]的公式为 t+1,在 t=1/2处给 3/2,真实格点数为一。无界多面体则可能有无穷多格点,不能使用这里的有限计数生成函数。空集的计数恒零,应单独处理,不满足本页非空对象的常数项一。

推论与应用

第一步:只用原顶点做三角剖分 ​

在 P的 r维仿射包内工作。若它已经是单纯形,无须分割。否则给每个顶点 vi选一个额外高度 hi,考察 (vi,hi)的凸包。对每个由 r+1个仿射无关原顶点决定的仿射图像,要求其余抬起的顶点不落在同一图像上。这只排除有限个关于 hi的真线性方程,因而可以选择满足要求的有理高度。

在每条竖直线上取抬起凸包的最低点,得到原 P上单值的分片仿射下表面。每个最高维下表面片恰有 r+1个顶点,否则违反高度选择;其投影因而是一个 r维单纯形。所有投影覆盖 P,两片的交来自下表面的公共面,不会在内部交叠。这给出不增加新顶点的三角剖分。

将每个单纯形同原点作锥,得到 C的有限单纯锥分割。每个锥的 r+1个生成元都选为 (qvi,q),所以高度统一为 q。

第二步:用同一个方向分配所有共享边界 ​

在 C的相对内部选 w,使它不落在任何小锥的facet所张成的超平面内。因为要避开的超平面有限,这样的 w存在。对任意 x∈C,充分小的 ε>0使 x+εw落入唯一一个小锥的相对内部;把 x分给该块。这是对每个点取足够小的扰动,不是声称一个固定步长对所有点同时有效。

在某块的基 u0,…,ur中写

x=∑λiui,w=∑βiui.

所有 βi≠0。该块对 x的精确条件是

(4)λi≥0(βi>0),λi>0(βi<0).

这些半开块两两不交,合起来恰为闭锥 C。因为 w的高度为正而每个 ui高度为 q,至少一个 βi为正;因此每块至少有一个非严格坐标。

第三步:把半开块分成有限基本域与非负整数步 ​

对式(4)中非严格坐标取基本区间 [0,1),对严格坐标取 (0,1]。分别用下取整、上取整减一,便把每个系数唯一写成

λi=ni+ρi,ni∈Z≥0,

其中 ρi位于相应基本区间。若 x为格点,则

p=x−∑niui=∑ρiui

也为格点。基本平行多面体有界,所以其中这样的 p有限;反过来,每个基本域格点加上任意非负整数步,都得到该半开块的唯一格点。

基本域点的整数高度满足 0≤h<q(r+1):各系数至多一,且至少一个区间不含一。于是该块的生成函数为

∑pzh(p)(1−zq)r+1.

对不交的块相加,得到式(1)及非负整数分子,再由式(2)证明全部非负整数参数上的准多项式性。

次数、体积与插值能证明到哪里 ​

qP含一个由整数顶点 v0,…,vr构成的 r维单纯形。对所有 ni≥0、∑i=1rni≤k,点

kv0+∑i=1rni(vi−v0)

互异且位于 kqP,故 LP(qk)≥(k+rr)。另一方面,选择在仿射包上单射的 r个环境坐标,伸缩后的每个坐标范围为 O(k),所以点数为 O(kr)。结合已证次数上界,最大次数恰为 r。

全维时,以整数点为角点的单位立方体比较 tP体积与格点数。两者的差只来自与边界相交的立方体;它们落在有限个facet的固定厚度邻域内,数量为 O(td−1)。因此

LP(t)=vold(P)td+O(td−1),

并迫使所有剩余类的最高项系数相同。二维整数多边形还可直接调用Pick公式,得到面积、半边界和常数三项。

次数和周期已被证明后,每个剩余类取 r+1个值,就唯一确定其多项式。先看几个数据再拟合,却不能独自证明这些次数与周期承诺。分母可能很大,三角剖分和基本域格点也可能很多;本页的有限构造不声称在任意维数下具有多项式位复杂度。

最小周期整除 q也可直接检查。将剩余类公式改写为原变量 t的幂次式后,各周期系数由计数函数唯一决定:比较两个周期的公倍数剩余类,两个多项式在无限多点相同,其系数便相同。若这些系数同时以 p和 q为周期,Bézout整数线性组合说明它们也以 gcd(p,q)为周期;取 p为最小周期,便得 p∣q。

分子的非负性,不等于普通系数都非负 ​

对 m≥1,Reeve四面体 Rm=conv{(0,0,0),(1,0,0),(0,1,0),(1,1,m)}的四个同高生成元为 (vi,1)。其半开基本域中的整数点为原点,以及

(1,1,j,2),1≤j<m.

具体地,顶端系数为 j/m,两个底面方向系数均为 1−j/m,原点方向系数为 j/m;高度总和为二。这些点列尽了基本域:第三坐标确定 j,前两坐标和高度再唯一确定其余系数。

所以

ERm(z)=1+(m−1)z2(1−z)4,LRm(t)=(t+33)+(m−1)(t+13)=m6t3+t2+12−m6t+1.

当 m>12时一次项系数为负,但全部非负整数 t的真实点数仍然非负。生成函数分子、二项式基系数与普通幂次系数是不同的表示。

Hilbert函数也可能因加权分次保留周期,但一般Hilbert多项式只保证最终相等;本页的小参数精确性来自基本域高度的严格上界。内部互反将进一步解释同一准多项式在负整数处的含义。完整复算见格点计数综合任务。

参考资料
  • Matthias Beck、Sinai Robins,Computing the Continuous Discretely,第二版更新稿,2020-06-19,§3.1 Theorem3.1,pp.59–62;§3.4 Theorem3.8,pp.68–70;§3.8 Theorem3.23,pp.80–81:三角剖分、整数多项式及有理准多项式。
  • Matthias Beck、Raman Sanyal,Combinatorial Reciprocity Theorems,2018-10-04稿,Proposition4.5.1,pp.120–122;Theorems4.7.2、4.8.1,pp.129、131–133;§5.3 Lemma5.3.4及Corollary5.3.5,pp.166–167:周期系数、半开基本域与边界分配。
  • Matthias Beck、Frank Sottile,Irrational proofs for three theorems of Stanley,2005,Theorem2及§3,pp.2、5–6:有理Ehrhart级数分子的非负性。本文用同一一般方向的半开块描述,完整写出式(2)的小参数检查,不依赖数值拟合。
关系图谱14 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系