English The Algo Lab · 论文导读

平面上的 Tusnády 问题

给 nn 个点染红蓝两色,要让每个矩形里红蓝尽量一样多。最好能做到多好?下界几十年来停在 log⁡n\log n,上界在 2017 年降到 log⁡3/2n\log^{3/2} n,中间还差一个 log⁡n\sqrt{\log n}。这份导读从一个染色游戏讲起,一步一步推到这个缺口是怎么补上的。


依据 · Tight Bounds for Tusnády's Problem in the Plane作者 · Zhewei Wei稿件 · r65 版arXiv · 2610.08130Lean · wzskytop/tusnady-lean
结论
Δ2(n)=Θ(log⁡3/2n)\Delta_2(n)=\Theta(\log^{3/2} n)
下界从 log⁡n\log n 抬到 log⁡3/2n\log^{3/2} n,与已知上界吻合
最坏情形
随机撒点
对任意 A>0A>0,以至少 1−e−An1-e^{-An} 的概率,任何染色都不够好
证明的主角
局部振荡
只用到阶梯函数的最大最小值、条件平均和几个初等不等式
目录
  1. I问题
    1. 1一个染色游戏
    2. 2从一维到二维
    3. 3记号与增长速度
  2. II之前为什么不行
    1. 4上下界的历史
    2. 5上界从哪来
    3. 6下界与转移原理
    4. 7其他下界工具
  3. III换一个问法
    1. 8一次对付一种染色
    2. 93/2 是怎么猜的
    3. 10逐位揭示
  4. IV局部振荡
    1. 11阶段与矩形
    2. 12端点平均与振荡
    3. 13先平均,再比较
    4. 14取期望:基准和增益
    5. 15一位数字的影响
    6. 16把 h 层接起来
    7. 17点不会太挤
    8. 18定理 1.1 的证明
  5. V回顾
    1. 193/2 的账
    2. 20为什么这次可以
    3. 21开放问题
    4. 22依赖图与记号表
导言

这份导读回答三个问题

论文的主结论只有一行:平面上轴平行矩形的组合差异度是 Θ(log⁡3/2n)\Theta(\log^{3/2} n)。围绕它有三个问题,全文按这个顺序展开。

之前为什么不行

旧的下界方法都停在 log⁡n\log n。第二部分逐个说明它们是怎么工作的,天花板在哪里。

这次为什么可以

论文换了问法:点随机撒;先固定一种染色,证明它“染得好”的概率小到 e−cne^{-cn},再把 2n2^n 种染色一并排除。第三部分讲这个转换。

有什么新技术

一个叫“局部振荡”的量。在好层上,横向细分一级、纵向合并一级,会带来正的条件期望增益。固定已经揭示的信息,只改变一个点的下一位横坐标数字,新势函数只变一点点。第四部分给出完整推导。

需要的预备知识

全文只假设大二水平:求和与不等式,指数和对数,概率、期望与独立性。条件期望、Cauchy–Schwarz 不等式这样的工具,第一次出现时都有一个标着“工具”的方框把它讲清楚,并给出证明。标着“推导”的折叠块是完整证明,默认展开;标着“选读”的块可以跳过,不影响主线。每节标题上方标了难度:入门、进阶、核心推导。

全景:五句话的证明

第 8 节在单位正方形里随机撒 nn 个点。固定一种染色 χ\chi,记 Fχ(x,y)F_\chi(x,y) 为锚定矩形 [0,x]×[0,y][0,x]\times[0,y] 里的颜色和。只要“FχF_\chi 处处都小”的概率比 2−n2^{-n} 小得多,就能用并集界一次排除全部 2n2^n 种染色。
第 11–12 节把正方形按 bb 进制切开。第 jj 阶段,横向区间宽 b−jb^{-j},纵向的块高 bj−hb^{j-h}。对每个横向区间取两个端点处的平均 aI(y)a_I(y),量它在每个块里的最大值减最小值,再对所有区间和块取平均,得到 ZjZ_j。有 Zj≥0Z_j\ge0 和 Zj≤2max⁡∣Fχ∣Z_j\le2\max|F_\chi|。
第 13–15 节横向细分一级、纵向合并一级,叫过一层,一共 hh 层。在每个“比较矩形”里,细分之后的振荡之和平均不降,而且多出一截,至少是 14mR\frac14\sqrt{m_R},mRm_R 是比较矩形的内部区域里的点数。另一方面,改动一个点的下一位数字,ZZ 至多变 1M\frac1M,其中 MM 是每层比较矩形的个数。
第 16–17 节把 ZZ 的总涨幅拆成各层的“条件增益”和一个“波动项”。只要点摊得开,一层的条件增益至少是 n/M\sqrt{n/M} 的常数倍,nM\frac nM 是一个比较矩形里的平均点数;hh 层上点都摊得开,这件事以极高的概率成立,而且与染色无关。波动项是大的负数的概率,由有界差分不等式控制,小到乘上 2n2^n 之后仍然很小。
第 18 节取层数 hh 为 log⁡n\log n 的量级、nM\frac nM 为 hh 的量级。于是除了概率至多 e−Ane^{-An} 的例外,每种染色都有 max⁡∣Fχ∣>h64b3/2n/M\max|F_\chi|>\frac h{64b^{3/2}}\sqrt{n/M},这是 h3/2h^{3/2} 也就是 log⁡3/2n\log^{3/2}n 的量级。

这些记号后面都会定义,第一遍只需要看个大意。

这份导读依据作者提供的论文 r65 版撰写。第四部分逐节对应论文第 2 节的证明:式 (1)–(18)、引理 2.1 和 2.3–2.6、Fact 2.2、命题 2.7,以及附录的 Fact A.1,编号都与 r65 版一致。论文的证明写得很紧凑,这里补上了省略的步骤,并加了例子和交互图;第 20 节有论文的四步概括和各节的对应。导读里补充的推导、例子和交互演示不能代替形式化验证,也不能代替独立的审阅。所有结论以论文为准;标明“背景”的内容来自相关文献,不在论文正文里。

Part I · 第 1–3 节

问题

先把问题本身弄清楚:染色游戏怎么玩,一维为什么容易,二维难在哪里,以及“log⁡3/2n\log^{3/2} n”这样的记号到底在说什么。

§ 1 · 入门一个染色游戏

平面上有 nn 个点。你的任务是把每个点染成红色或蓝色。染完之后,对手可以任选一个边与坐标轴平行的矩形,数一数矩形里红点比蓝点多几个,或者少几个。这个差的绝对值就是你的损失。

你想让损失尽量小,而且颜色必须在对手挑矩形之前定好。所以一种染色的好坏,要看它在最坏的那个矩形上的表现。

交互 · 01染色游乐场

点一下某个点可以翻转它的颜色。方框标出当前染色下最坏的矩形。试着手动把差异度降到最优值。

多试几组点会发现两件事。第一,点数很少时,最优值通常只有 1 或 2。第二,随手染的颜色往往离最优差得不少。整篇论文关心的问题是:点数 nn 很大时,这个最优值最坏会有多大。

为了方便计算,以后把红色记作 +1+1,蓝色记作 −1-1。矩形里“红点数减蓝点数”就是矩形内所有点的颜色之和。

§ 2 · 入门从一维到二维

一维:交替染色就够了

先看简单的情形。如果所有点都在一条直线上,“矩形”就退化成区间。把点从左到右排好,红蓝交替地染。

这个区间里:4 红 3 蓝,和为 +1 实心 = +1(红) 空心 = −1(蓝)
一个区间里的点在排序中是连续的一段。连续一段交替的 ±1\pm1 相加,结果只能是 −1-1、00 或 +1+1。
一维的答案热身

直线上任意 n≥1n\ge1 个点,按顺序交替染色后,每个区间里的颜色之和的绝对值不超过 11。所以一维情形的差异度恒为 11,与 nn 无关。

二维:只按一个坐标交替不够

二维的点有两个顺序:按 xx 坐标排一次,按 yy 坐标排一次。如果只按 xx 坐标交替染色,那么每个竖直的条带都很平衡,但水平条带可能很糟。上面的游乐场里有一个“特制点集”按钮:它把 xx 顺序中排奇数位的点都放在下半部分,偶数位的都放在上半部分。这时按 xx 交替染色,下半个正方形里就全是红点。

更麻烦的是,矩形同时对 xx 和 yy 提要求。落在一个矩形里的点,在 xx 顺序里不一定连续,在 yy 顺序里也不一定连续。一维那个“连续一段”的论证完全用不上。

Tusnády 当年问的就是:二维能不能也做到某个与 nn 无关的常数。Beck 在 1981 年给出了否定的回答,同时证明差异度至多是 log⁡4n\log^4 n 量级。从那以后,问题变成了确定它的准确增长速度。

只看“锚定矩形”就够了

先把点集平移到第一象限,让所有坐标都是正的。左下角固定在原点的矩形 A(x,y)=[0,x]×[0,y]A(x,y)=[0,x]\times[0,y] 叫锚定矩形。它只有两个参数,比一般矩形好处理得多。而且只看锚定矩形最多损失 4 倍。

R 减去 A(x1, y2) 减去 A(x2, y1) 加回 A(x1, y1) x1 x2 y1 y2 0 外框是 A(x2, y2)
左下角那块被减了两次,所以要加回一次。这就是容斥原理。

用式子写出来:对矩形 R=(x1,x2]×(y1,y2]R=(x_1,x_2]\times(y_1,y_2] 和任意染色 χ\chi,

χ(P∩R)=χ(P∩A(x2,y2))−χ(P∩A(x1,y2))−χ(P∩A(x2,y1))+χ(P∩A(x1,y1)).\begin{aligned}\chi(P\cap R)={}&\chi\bigl(P\cap A(x_2,y_2)\bigr)-\chi\bigl(P\cap A(x_1,y_2)\bigr)\\&-\chi\bigl(P\cap A(x_2,y_1)\bigr)+\chi\bigl(P\cap A(x_1,y_1)\bigr).\end{aligned}

右边四项每一项的绝对值都不超过“锚定矩形上的最大差”,所以一般矩形上的差至多是它的 4 倍。反过来,锚定矩形本身就是矩形。因此两种定义下的差异度只差一个不超过 4 的常数因子,研究增长速度时可以只看锚定矩形。论文的下界甚至直接对锚定矩形成立。

§ 3 · 入门记号与增长速度

三层的 min 和 max

设 PP 是平面上的有限点集。一个染色是函数 χ:P→{+1,−1}\chi:P\to\{+1,-1\}。对平面上的集合 RR,记

χ(P∩R)=∑p∈P∩Rχ(p).\chi(P\cap R)=\sum_{p\in P\cap R}\chi(p).

点集 PP 的差异度,和 nn 个点的最坏差异度,定义为

disc⁡(P)=min⁡χ max⁡R ∣χ(P∩R)∣,Δ2(n)=max⁡∣P∣=ndisc⁡(P),\disc(P)=\min_{\chi}\ \max_{R}\ \bigl|\chi(P\cap R)\bigr|,\qquad \Delta_2(n)=\max_{|P|=n}\disc(P),

其中 RR 取遍所有轴平行矩形。下标 2 表示平面。这三层可以读成一场三步的博弈:

max对手先选 nn 个点的位置,想让你难办。
min你看到点之后选染色 χ\chi,想让最坏的矩形也平衡。
max对手看到染色之后选矩形 RR,挑差最大的那个。
易错

证明上界和证明下界是两种完全不同的任务,后面会反复用到这个区别。

上界 Δ2(n)≤K\Delta_2(n)\le K:对每一个 nn 点集,都要找出一种染色,使所有矩形的差不超过 KK。

下界 Δ2(n)>K\Delta_2(n)> K:要找出一个 nn 点集,使每一种染色都有某个矩形的差超过 KK。染色有 2n2^n 种,要一个不漏地全部排除。

渐近记号

我们关心的是 n→∞n\to\infty 时的增长速度,不关心常数倍。三个记号的意思是:

  • f(n)=O(g(n))f(n)=O(g(n)):存在常数 CC,使 nn 足够大时 f(n)≤C g(n)f(n)\le C\,g(n)。读作“至多是 gg 的量级”。
  • f(n)=Ω(g(n))f(n)=\Omega(g(n)):存在常数 c>0c>0,使 nn 足够大时 f(n)≥c g(n)f(n)\ge c\,g(n)。读作“至少是 gg 的量级”。
  • f(n)=Θ(g(n))f(n)=\Theta(g(n)):以上两条同时成立,量级完全确定。

另外,log⁡3/2n\log^{3/2}n 是 (log⁡n)3/2(\log n)^{3/2} 的简写。对数取什么底不影响量级,因为换底只差常数倍。需要具体数值的地方,会写明用的是以 2 为底的 log⁡2\log_2 还是自然对数 ln⁡\ln。

交互 · 02log⁡n\log n 和 log⁡3/2n\log^{3/2}n 差多少

拖动滑块改变点数 nn。这里忽略所有常数因子,只比较函数本身的大小。

从数值上看,两者只差 log⁡2n\sqrt{\log_2 n} 倍:一百万个点时大约是 4.5 倍。和随机染色能保证的 nlog⁡n\sqrt{n\log n} 相比,两者都小得可怜。但数学上要问的是真实的增长速度到底是哪一个,这不能靠“差不多”来回答。

本文的结论论文 Corollary 1.2,由 Theorem 1.1 和 Nikolov 的上界推出

存在常数 c>0c>0,使得对每个 n≥1n\ge1 有 Δ2(n)≥clog⁡3/2n\Delta_2(n)\ge c\log^{3/2}n。结合已知的上界,Δ2(n)=Θ(log⁡3/2n)\Delta_2(n)=\Theta(\log^{3/2}n)。

论文实际证明的是一个更强的说法:随机撒的点集几乎总是这么难染。准确的表述在第 8 节。

Part II · 第 4–7 节

之前知道什么,为什么卡住

上界一路从 log⁡4n\log^4 n 降到 log⁡3/2n\log^{3/2}n,下界却一直停在 log⁡n\log n。这一部分把两边的方法都过一遍,重点看下界的几种工具各自停在哪里,为什么。

§ 4 · 入门上下界的历史

这一部分是背景材料。论文只在引言和第 3 节里各用一段话概括了这些工作;这里的推导取自标准文献,个别地方是这份导读补充的,会在出现的地方注明。

上界和下界是用两套完全不同的工具分别推进的。把已知结果画在一张图上,纵轴是 log⁡n\log n 的指数。

交互 · 03Δ2(n)\Delta_2(n) 的上下界

把鼠标移到圆点上(或用 Tab 键聚焦)可以看到出处。两条线之间的区域是当时还不知道的范围。

表格视图(论文 Table 1)
界出处年份备注
O(log⁡4n)O(\log^4 n)Beck1981
O(log⁡3.5+εn)O(\log^{3.5+\varepsilon} n)Beck1989对每个固定的 ε>0\varepsilon>0
O(log⁡3n)O(\log^3 n)Bohus1990
O(log⁡2.5n)O(\log^{2.5} n)Srinivasan1997
O(log⁡2n)O(\log^2 n)Bansal–Garg2017构造性的,能给出算法
O(log⁡3/2n)O(\log^{3/2} n)Nikolov2017
Ω(log⁡n)\Omega(\log n)Schmidt1972经转移原理得到
Ω(log⁡n)\Omega(\log n)Matoušek–Nikolov2015分解范数;dd 维是 Ω(log⁡d−1n)\Omega(\log^{d-1}n)
Ω(log⁡3/2n)\Omega(\log^{3/2} n)本文随机点即可达到

图里能读出三件事。

  • 上界那条线一级一级往下降,最后在 2017 年到达 log⁡3/2n\log^{3/2}n。
  • 下界那条线一直是平的。Schmidt 1972 年的定理说的是连续差异度,要经过第 6 节的转移原理才变成染色问题的下界。2015 年换了一种全新的方法,得到的仍然是 log⁡n\log n。
  • 于是剩下一个 log⁡n\sqrt{\log n} 的缺口。哪一边是对的,事先并不清楚:论文提到,Bansal 在 2022 年国际数学家大会的报告里猜想下界 log⁡n\log n 才是对的。论文的结论是另一边:上界 log⁡3/2n\log^{3/2}n 是对的。Bansal 和 Nikolov 2026 年的综述也讨论了平面上剩下的这个缺口。

缺口两头都有道理。log⁡n\log n 是另一个经典问题(连续差异度,见第 6 节)的正确答案。log⁡3/2n\log^{3/2}n 则是最好的染色方法做到的。到底是染色方法还能改进,还是下界的证明不够强,一直没有定论。

§ 5 · 进阶上界从哪来

本文证明的是下界。但上界的想法里有一个结构,下界的证明里会再遇到:把问题拆成 log⁡n\log n 层尺度,每一层付一份代价。所以值得先看一眼。

最朴素的办法:随机染色

每个点独立地抛硬币决定颜色。这样能做到什么程度?需要三个概率工具,其中 Markov 不等式和并集界后面还会再用。

工具 · Markov 不等式

设随机变量 X≥0X\ge0,a>0a>0。则 Pr⁡[X≥a]≤E[X]/a\Pr[X\ge a]\le \EE[X]/a。

证明只有一行:E[X]≥E[X⋅1{X≥a}]≥a⋅Pr⁡[X≥a]\EE[X]\ge \EE\bigl[X\cdot\one_{\{X\ge a\}}\bigr]\ge a\cdot\Pr[X\ge a]。意思是,一个非负的量如果平均值很小,它就不太可能很大。

工具 · 并集界

对任意事件 E1,…,ENE_1,\dots,E_N,有 Pr⁡[E1∪⋯∪EN]≤∑iPr⁡[Ei]\Pr[E_1\cup\dots\cup E_N]\le\sum_i\Pr[E_i]。“至少发生一件坏事”的概率,不超过各件坏事概率之和。证明:并集的示性变量不超过各事件示性变量之和,两边取期望。

工具 · ±1\pm1 之和的尾概率

设 S=ϵ1+⋯+ϵsS=\eps_1+\dots+\eps_s,各 ϵi\eps_i 独立,等概率取 ±1\pm1。则对任意 a>0a>0,Pr⁡[ ∣S∣≥a ]≤2e−a2/(2s)\Pr[\,|S|\ge a\,]\le 2e^{-a^2/(2s)}。

推导尾概率的证明
  1. 对任意 t>0t>0,由独立性,E[etS]=∏iE[etϵi]=(cosh⁡t)s\EE[e^{tS}]=\prod_i\EE[e^{t\eps_i}]=(\cosh t)^s,其中 cosh⁡t=12(et+e−t)\cosh t=\tfrac12(e^t+e^{-t})。
  2. 逐项比较幂级数:cosh⁡t=∑k≥0t2k(2k)!\cosh t=\sum_{k\ge0}\dfrac{t^{2k}}{(2k)!},et2/2=∑k≥0t2k2k k!e^{t^2/2}=\sum_{k\ge0}\dfrac{t^{2k}}{2^k\,k!}。因为 (2k)!=k!⋅(k+1)(k+2)⋯(2k)≥k!⋅2k(2k)!=k!\cdot(k+1)(k+2)\cdots(2k)\ge k!\cdot2^k,所以 cosh⁡t≤et2/2\cosh t\le e^{t^2/2}。
  3. 对非负随机变量 etSe^{tS} 用 Markov 不等式:Pr⁡[S≥a]=Pr⁡[etS≥eta]≤e−ta E[etS]≤est2/2−ta\Pr[S\ge a]=\Pr[e^{tS}\ge e^{ta}]\le e^{-ta}\,\EE[e^{tS}]\le e^{st^2/2-ta}。
  4. 取 t=a/st=a/s,指数变成 −a2/(2s)-a^2/(2s)。由对称性 Pr⁡[S≤−a]\Pr[S\le -a] 有同样的界,两者相加得结论。

现在固定一个 nn 点集和一个矩形 RR。RR 里有 s≤ns\le n 个点,随机染色下 χ(P∩R)\chi(P\cap R) 正是 ss 个独立 ±1\pm1 之和。取 a=10 nln⁡na=\sqrt{10\,n\ln n},

Pr⁡[ ∣χ(P∩R)∣≥a ]≤2e−a2/(2n)=2n−5.\Pr\bigl[\,|\chi(P\cap R)|\ge a\,\bigr]\le 2e^{-a^2/(2n)}=2n^{-5}.

矩形虽然有无穷多个,但它们截出的点子集是有限的:子集由“xx 顺序中的哪一段”和“yy 顺序中的哪一段”决定,至多 (n(n+1)2)2≤n4\bigl(\tfrac{n(n+1)}2\bigr)^2\le n^4 种。由并集界,存在某个矩形超标的概率至多 n4⋅2n−5=2/nn^4\cdot2n^{-5}=2/n,当 n≥3n\ge3 时小于 1。所以一定有一种染色让所有矩形都不超标:

Δ2(n)≤10 nln⁡n(n≥3).\Delta_2(n)\le\sqrt{10\,n\ln n}\qquad(n\ge3).

这个界离真相差得很远,真相是 log⁡n\log n 的幂次。但这段推导展示了一种基本手法:用概率证明存在性。第三部分会把同样的手法反过来用。

把前缀拆成二进区间

要做到 log⁡n\log n 的幂次,需要利用矩形的结构。把 xx 方向分成 m=2hm=2^h 列。长度为 22 的幂、位置对齐的区间叫二进区间:第 0 层是整个 [1,m][1,m],第 1 层是左右两半,一直到第 hh 层的单独一列。

交互 · 04前缀的二进分解

拖动滑块改变前缀 [1,v][1,v]。高亮的块是它的分解;每一层至多用到一块。

有两个事实,后面都会用到:

  • 每一列恰好属于 h+1h+1 个二进区间,每层一个。
  • 每个前缀 [1,v][1,v] 都能写成至多 h+1h+1 个互不相交的二进区间的并,每层至多一个。这其实就是 vv 的二进制展开。

于是一个锚定矩形(宽度是前缀 [1,v][1,v],高度到某处为止)可以拆成至多 h+1h+1 个“二进条带”。如果一种染色能让每个二进条带内部的误差都不超过 EE,那么锚定矩形的误差就不超过 (h+1)E(h+1)E。每个条带单独看是个一维问题,容易。困难在于同一种染色要同时照顾全部 h+1h+1 层。

历史上的上界大多沿着“分层”的思路走。早期的证明用的是部分染色的想法:每一轮只给一部分点定色,同时控制一族标准矩形上的误差,剩下的点留到下一轮。每一轮、每一层尺度都有损失,所以得到的是 log⁡n\log n 的高次幂。后来的改进主要是在减少这些损失。

目前最好的上界

Nikolov 在 2017 年证明了 Δ2(n)=O(log⁡3/2n)\Delta_2(n)=O(\log^{3/2}n)。证明依赖 Banaszczyk 的一个深刻定理,定理本身远超大二水平,这里把它当作黑箱,只看它怎么和二进分解搭配。

选读从 Banaszczyk 定理到 O(log⁡3/2n)O(\log^{3/2}n) 的骨架
Banaszczyk 定理(前缀和版本)背景 · 2012

设 v1,…,vn∈Rdv_1,\dots,v_n\in\R^d 的长度都不超过 1,K⊂Rd\mathcal K\subset\R^d 是关于原点对称的凸体,并且一个标准高斯随机向量落在 K\mathcal K 里的概率至少是 1−12n1-\frac1{2n}。则存在符号 ϵi∈{±1}\eps_i\in\{\pm1\},使所有前缀和 ∑i≤uϵivi\sum_{i\le u}\eps_iv_i(u=1,…,nu=1,\dots,n)都落在 5K5\mathcal K 里。

  1. 把点按高度排成 p1,…,pnp_1,\dots,p_n,记 rir_i 为 pip_i 在 xx 方向的名次。锚定矩形上的颜色和是 A(u,v)=∑i≤u, ri≤vϵiA(u,v)=\sum_{i\le u,\ r_i\le v}\eps_i。
  2. 取 h=⌈log⁡2n⌉h=\lceil\log_2n\rceil,把名次 1,…,2h1,\dots,2^h 当作列。给每个点一个向量 viv_i,坐标用二进区间 II 编号:若 ri∈Ir_i\in I,该坐标是 1/h+11/\sqrt{h+1},否则是 0。每个点恰在 h+1h+1 个二进区间里,所以 ∥vi∥=1\norm{v_i}=1。
  3. 对每个 vv,令 bvb_v 为前缀 [1,v][1,v] 的二进分解的示性向量,它至多有 h+1h+1 个 1。分解是不相交的,所以 ⟨bv, vi⟩=1{ri≤v}/h+1\ip{b_v}{v_i}=\one_{\{r_i\le v\}}/\sqrt{h+1}。
  4. 取 K={y: ∣⟨bv, y⟩∣≤t 对 v=1,…,n}\mathcal K=\{y:\ |\ip{b_v}{y}|\le t\ \text{对}\ v=1,\dots,n\},它是关于原点对称的凸集(严格地说要限制在这些 bvb_v 张成的子空间里才有界,这不影响论证)。对高斯向量 gg,⟨bv, g⟩\ip{b_v}{g} 服从方差不超过 h+1h+1 的正态分布,所以 Pr⁡[∣⟨bv, g⟩∣>t]≤2e−t2/(2(h+1))\Pr[|\ip{b_v}{g}|>t]\le2e^{-t^2/(2(h+1))}。对这 nn 个 vv 用并集界,取 t=2(h+1)ln⁡(4n2)t=\sqrt{2(h+1)\ln(4n^2)} 就满足定理的条件。
  5. 定理给出的符号满足 ∣∑i≤uϵi⟨bv, vi⟩∣≤5t\bigl|\sum_{i\le u}\eps_i\ip{b_v}{v_i}\bigr|\le5t,也就是 ∣A(u,v)∣≤5th+1=5(h+1)2ln⁡(4n2)=O(log⁡3/2n).|A(u,v)|\le5t\sqrt{h+1}=5(h+1)\sqrt{2\ln(4n^2)}=O(\log^{3/2}n).

两个因子的来源很清楚:h+1≈log⁡nh+1\approx\log n 来自二进分解的层数,log⁡n\sqrt{\log n} 来自要求 nn 个约束同时成立。原文处理的是一般维数,用的是分解范数的语言,细节以原文为准。

上界的故事到这里停住了。它留下一个印象:有 log⁡n\log n 层,每层要多付 log⁡n\sqrt{\log n}。本文证明的是,这个代价一分也省不掉。

§ 6 · 进阶下界从哪来:连续差异度与转移原理

一个更老的问题:点放得均不均匀

不染色,只问:在单位正方形 [0,1]2[0,1]^2 里放 nn 个点,怎样才算均匀。如果点是均匀的,面积为 aa 的矩形里应该有大约 anan 个点。点集 PP 的连续差异度定义为最坏的偏差

D(P)=sup⁡R ∣ ∣P∩R∣−n⋅area(R)∣,D(n)=min⁡∣P∣=nD(P),D(P)=\sup_{R}\ \Bigl|\,|P\cap R|-n\cdot\mathrm{area}(R)\Bigr|,\qquad D(n)=\min_{|P|=n}D(P),

其中 RR 取遍 [0,1]2[0,1]^2 里的轴平行矩形。注意这里是你自己放点,所以取 min⁡\min。

交互 · 05三种点集有多均匀

这里计算的是锚定矩形上的最大偏差(它与 D(P)D(P) 至多差 4 倍)。阴影是偏差最大的那个锚定矩形。

网格和随机点的偏差都在 n\sqrt n 量级。第三种点集叫二进网,它的偏差只有 log⁡n\log n 量级。下面把这个结论推出来,因为推导里又出现了“log⁡n\log n 层,每层贡献常数”的结构。

二进网:每层误差不超过 1

设 N=2sN=2^s。把整数 i∈{0,…,N−1}i\in\{0,\dots,N-1\} 写成 ss 位二进制,再把这 ss 位倒过来读成一个二进制小数,记作 φ(i)\varphi(i)。例如 s=3s=3 时 i=6=(110)2i=6=(110)_2,φ(6)=(0.011)2=3/8\varphi(6)=(0.011)_2=3/8。点集取为

P={ (i/N, φ(i)) : i=0,1,…,N−1 }.P=\bigl\{\,(i/N,\ \varphi(i))\ :\ i=0,1,\dots,N-1\,\bigr\}.
网性质背景 · Hammersley 点集

设 p+q=sp+q=s。每个形如 [a2p,a+12p)×[b2q,b+12q)\bigl[\frac a{2^p},\frac{a+1}{2^p}\bigr)\times\bigl[\frac b{2^q},\frac{b+1}{2^q}\bigr) 的盒子(面积恰为 1/N1/N)里恰好有 PP 的一个点。

推导网性质,以及 D(P)≤4(log⁡2N+1)D(P)\le 4(\log_2N+1)
  1. 网性质。i/Ni/N 落在第一个区间里,等价于 ii 的最高 pp 位是 aa 的二进制。φ(i)\varphi(i) 落在第二个区间里,等价于 φ(i)\varphi(i) 的前 qq 位小数是 bb 的二进制,也就是 ii 的最低 qq 位被确定。p+q=sp+q=s,所以 ii 的全部 ss 位都被确定,恰有一个 ii。
  2. 拆分。考虑锚定矩形 [0,X)×[0,Y)[0,X)\times[0,Y)。由二进制展开,[0,X)[0,X) 可以拆成至多 ss 个互不相交的二进区间(长度是互不相同的 2−p2^{-p},1≤p≤s1\le p\le s),外加一段长度小于 1/N1/N 的零头。X=1X=1 时它就是 [0,1)[0,1) 这一个区间。
  3. 每个二进条带。设 II 是长度 2−p2^{-p} 的二进区间。横坐标在 II 里的点有 2s−p2^{s-p} 个。由网性质(取 q=s−pq=s-p),它们的纵坐标在每个 [b2q,b+12q)\bigl[\frac b{2^q},\frac{b+1}{2^q}\bigr) 里恰有一个。所以 I×[0,Y)I\times[0,Y) 里的点数在 ⌊Y2q⌋\lfloor Y2^q\rfloor 和 ⌈Y2q⌉\lceil Y2^q\rceil 之间,而“应有的点数”是 N⋅2−p⋅Y=Y2qN\cdot2^{-p}\cdot Y=Y2^q。误差不超过 1。
  4. 零头。横坐标都是 1/N1/N 的倍数,长度小于 1/N1/N 的区间里至多一个点,应有点数小于 1。误差不超过 1。
  5. 合计。锚定矩形的偏差至多 s+1=log⁡2N+1s+1=\log_2N+1。一般矩形是四个锚定矩形的加减,所以 D(P)≤4(log⁡2N+1)D(P)\le4(\log_2N+1)。

任意点数。对任意 n≥2n\ge2,取点 (i/n,φ(i))(i/n,\varphi(i)),其中 0≤i<n0\le i<n,φ\varphi 把 ii 的二进制数字倒序写成小数(前导零不影响结果)。锚定矩形选中的横坐标构成前 kk 个指标,且 ∣k−nX∣≤1|k-nX|\le1。把这段指标分成长度为 2 的幂、起点与长度对齐的块,至多需要 ⌊log⁡2n⌋+1\lfloor\log_2 n\rfloor+1 块。每一块的倒序数字在等长二进纵区间里各放一个点,所以计数误差至多为 1。把各块误差与 ∣k−nX∣Y≤1|k-nX|Y\le1 相加,就对任意 nn 得到 O(log⁡n)O(\log n);一般矩形再付至多四倍。

反方向是 Schmidt 在 1972 年证明的定理:无论怎么放点,D(n)≥clog⁡nD(n)\ge c\log n。它的证明这里不展开。两者合起来,连续差异度的答案完全确定:

D(n)=Θ(log⁡n).D(n)=\Theta(\log n).

只看锚定矩形的版本是 sup⁡x,y∣ ∣P∩([0,x]×[0,y])∣−nxy ∣\sup_{x,y}\bigl|\,|P\cap([0,x]\times[0,y])|-nxy\,\bigr|,它的最小值同样是 Θ(log⁡n)\Theta(\log n)。把它除以 nn,就是常说的星差异度,最优量级是 (log⁡n)/n(\log n)/n。

转移原理:好染色能造出均匀点集

连续差异度和组合差异度之间有一座桥,叫转移原理。它的证明是初等的,而且很有意思:反复染色、只留红点。

转移原理(简化形式)背景 · 标准论证的简化版,常数由本导读给出

对所有 n≥2n\ge2,有 D(n)≤2 Δ2(2n2)+5D(n)\le2\,\Delta_2(2n^2)+5。也就是说,如果任何点集都能染得很平衡,就一定存在很均匀的点集。

P02k·n 个点 P12k−1·n 个点 P22k−2·n 个点 Pkn 个点 染色,留红点染色,留红点 每一步:D(Pj+1) ≤ D(Pj) / 2 + (这一步染色的差异度)
旧的偏差每步减半,新引入的偏差就是那一步染色的差异度。
推导转移原理的证明
  1. 减半。设 PP 有 NN 个点(NN 为偶数),χ\chi 是一种染色,所有矩形上 ∣χ(P∩R)∣≤f|\chi(P\cap R)|\le f。令 P+P^+ 为红点的集合。对每个矩形 RR,红点数是 ∣P+∩R∣=12(∣P∩R∣+χ(P∩R))|P^+\cap R|=\frac12\bigl(|P\cap R|+\chi(P\cap R)\bigr),所以 ∣ ∣P+∩R∣−12∣P∩R∣ ∣≤f2.\Bigl|\,|P^+\cap R|-\tfrac12|P\cap R|\,\Bigr|\le\tfrac f2. 取 RR 为整个正方形,得 ∣∣P+∣−N2∣≤f2\bigl||P^+|-\frac N2\bigr|\le\frac f2。
  2. 修正点数。P+P^+ 的点数可能不恰好是 N/2N/2,但只差至多 f/2f/2 个。删掉或补上这么多点,得到恰有 N/2N/2 个点的 P′P'。每动一个点,任何矩形里的点数至多变 1,所以 ∣∣P′∩R∣−12∣P∩R∣∣≤f\bigl||P'\cap R|-\frac12|P\cap R|\bigr|\le f。
  3. 连续差异度的递推。由三角不等式, ∣ ∣P′∩R∣−N2 area(R)∣≤∣ ∣P′∩R∣−12∣P∩R∣ ∣+12∣ ∣P∩R∣−N area(R)∣≤f+12D(P).\begin{aligned}\Bigl|\,|P'\cap R|-\tfrac N2\,\mathrm{area}(R)\Bigr|&\le\Bigl|\,|P'\cap R|-\tfrac12|P\cap R|\,\Bigr|+\tfrac12\Bigl|\,|P\cap R|-N\,\mathrm{area}(R)\Bigr|\\&\le f+\tfrac12D(P).\end{aligned} 对 RR 取上确界:D(P′)≤12D(P)+fD(P')\le\frac12D(P)+f。
  4. Δ2\Delta_2 关于点数单调不减。在一个点集的远处再添一个点,原来的每个矩形都可以缩小一点避开新点,所以差异度不会变小。
  5. 迭代。取 k=⌈log⁡2n⌉k=\lceil\log_2n\rceil,从 N0=2knN_0=2^kn 个点的 P0P_0 出发做 kk 次减半,PjP_j 有 2k−jn2^{k-j}n 个点。因为 2k<2n2^k<2n,每个 PjP_j 的点数都小于 2n22n^2,所以每一步都能找到差异度不超过 f=Δ2(2n2)f=\Delta_2(2n^2) 的染色。把递推展开, D(Pk)≤D(P0)2k+f⋅(1+12+14+… )≤D(P0)2k+2 Δ2(2n2).D(P_k)\le\frac{D(P_0)}{2^k}+f\cdot\Bigl(1+\frac12+\frac14+\dots\Bigr)\le\frac{D(P_0)}{2^k}+2\,\Delta_2(2n^2).
  6. 起点。取 P0P_0 为一个 ⌈N0⌉×⌈N0⌉\lceil\sqrt{N_0}\rceil\times\lceil\sqrt{N_0}\rceil 的网格去掉多余的点。直接数格点可得 D(P0)≤5N0D(P_0)\le5\sqrt{N_0}。于是 D(P0)/2k≤5n/2k≤5D(P_0)/2^k\le5\sqrt{n/2^k}\le5,因为 2k≥n2^k\ge n。
  7. 结论。PkP_k 恰有 nn 个点,所以 D(n)≤D(Pk)≤2 Δ2(2n2)+5D(n)\le D(P_k)\le2\,\Delta_2(2n^2)+5。

把它和 Schmidt 定理 D(n)≥clog⁡nD(n)\ge c\log n 放在一起,得到 Δ2(2n2)≥12(clog⁡n−5)\Delta_2(2n^2)\ge\frac12(c\log n-5)。再用一次单调性,对所有 NN 都有 Δ2(N)=Ω(log⁡N)\Delta_2(N)=\Omega(\log N)。这就是经典下界的来源。

天花板

这条路线的逻辑是“组合差异度不会比连续差异度小”。可是连续差异度恰好就是 Θ(log⁡n)\Theta(\log n),二进网已经做到了。所以单靠这条路线,证不出比 log⁡n\log n 更大的下界。

转移原理是单向的:好染色能造出均匀点集,但均匀点集的存在不保证有好染色。本文的结论正说明两者确实不同:平面上组合差异度比连续差异度大一个 log⁡n\sqrt{\log n} 因子。

§ 7 · 进阶其他下界工具,和它们停下的地方

除了转移原理,还有两类方法给出过 Ω(log⁡n)\Omega(\log n)。它们的思路相近:对每一种染色,用同一套办法给出“存在坏矩形”的证据。

正交函数法

这是 Roth 在 1954 年引入的方法。把偏差写成正方形上的函数 D(x,y)D(x,y),再找一个检验函数 FF,使内积 ⟨D, F⟩=∬DF\ip{D}{F}=\iint DF 很大而 FF 本身不大。由 Cauchy–Schwarz 不等式,DD 的平方平均就不可能小。检验函数用 Haar 函数拼成:每个二进盒子上放一个取值 ±1\pm1 的小波,不同形状的盒子给出互相正交的函数。

选读Roth 定理的完整推导:D(N)≥1128log⁡2ND(N)\ge\frac1{128}\sqrt{\log_2N}

设 P⊂[0,1)2P\subset[0,1)^2 有 NN 个点,D(x,y)=#(P∩[0,x)×[0,y))−NxyD(x,y)=\#\bigl(P\cap[0,x)\times[0,y)\bigr)-Nxy。

  1. 尺度。取整数 nn 使 2N≤2n<4N2N\le2^n<4N。对每个 r=(r1,r2)r=(r_1,r_2),r1+r2=nr_1+r_2=n,把正方形切成 2n2^n 个大小为 2−r1×2−r22^{-r_1}\times2^{-r_2} 的二进盒子,每个面积 2−n2^{-n}。这样的 rr 共有 n+1n+1 个。
  2. Haar 函数。对区间 II,令 hIh_I 在左半取 +1+1,右半取 −1-1,区间外为 0。对盒子 R=I×JR=I\times J,令 hR(x,y)=hI(x)hJ(y)h_R(x,y)=h_I(x)h_J(y)。
  3. 空盒子上的内积。直接积分,∫Ix hI(x) dx=−∣I∣2/4\int_I x\,h_I(x)\,dx=-|I|^2/4,所以 ∬Nxy hR=N∣R∣2/16\iint Nxy\,h_R=N|R|^2/16。另一方面,如果 RR 里没有 PP 的点,那么每个点 pp 的计数函数 1{px<x}1{py<y}\one_{\{p_x<x\}}\one_{\{p_y<y\}} 在 RR 上至少对一个变量是常数,而 hIh_I、hJh_J 的积分都是 0,所以计数部分与 hRh_R 的内积为 0。于是空盒子满足 ⟨D, hR⟩=−N∣R∣2/16\ip{D}{h_R}=-N|R|^2/16。
  4. 一层的检验函数。令 fr=−∑hRf_r=-\sum h_R,对所有类型为 rr 的空盒子求和。因为 N≤2n−1N\le2^{n-1},空盒子至少有 2n−12^{n-1} 个,所以 ⟨D, fr⟩≥2n−1⋅N 4−n16=N2−n32>1128,∥fr∥2≤2n⋅2−n=1.\ip{D}{f_r}\ge2^{n-1}\cdot\frac{N\,4^{-n}}{16}=\frac{N2^{-n}}{32}>\frac1{128},\qquad\norm{f_r}^2\le2^n\cdot2^{-n}=1.
  5. 正交。若 r≠r′r\ne r',不妨设 r1<r1′r_1<r_1'。类型 r′r' 的盒子在 xx 方向更窄,它的 xx 区间要么与 II 不相交,要么整个落在 II 的某一半里,此时 hIh_I 在其上是常数,积分为 0。所以 ⟨fr, fr′⟩=0\ip{f_r}{f_{r'}}=0。
  6. 合起来。令 F=∑rfrF=\sum_rf_r。则 ⟨D, F⟩>n+1128\ip{D}{F}>\frac{n+1}{128},∥F∥2=∑r∥fr∥2≤n+1\norm{F}^2=\sum_r\norm{f_r}^2\le n+1。由 Cauchy–Schwarz, sup⁡∣D∣ ≥ ∥D∥2 ≥ ⟨D, F⟩∥F∥>n+1128≥log⁡2N128.\sup|D|\ \ge\ \norm{D}_2\ \ge\ \frac{\ip{D}{F}}{\norm{F}}>\frac{\sqrt{n+1}}{128}\ge\frac{\sqrt{\log_2N}}{128}.

每一层 rr 贡献一个常数,共有 n+1≈log⁡Nn+1\approx\log N 层。在平方平均的意义下,这些贡献按“平方和开根号”合并,得到 log⁡N\sqrt{\log N}。Schmidt 1972 年的定理把它提高到了 log⁡N\log N。后来 Halász 在 Roth 的框架里用更精巧的检验函数重新证明了这个结果,让各层的贡献直接相加。

把同样的想法用到染色上,Wei 和 Yi 证明了二进网的组合差异度是 Ω(log⁡n)\Omega(\log n)。这类证明的收益是按尺度一层一层累加的:共 log⁡n\log n 层,每层贡献一个常数。

在第 6 节那个二进网上,这已经是能证明的极限了,因为它本身就能染得这么好。下面的验证是这份导读补充的,不在论文里,只用到第 6 节的工具。

推导第 6 节的二进网可以染到差异度 4(log⁡2N+1)4(\log_2N+1)
  1. 染色。第 ii 个点 (i/N,φ(i))(i/N,\varphi(i)) 的颜色取为 (−1)b(i)(-1)^{b(i)},其中 b(i)b(i) 是 ii 的二进制里 1 的个数。
  2. 一个二进条带里的颜色。设 II 是长度 2−p2^{-p} 的二进区间,q=s−pq=s-p。横坐标在 II 里的点是 i=a⋅2q+ji=a\cdot2^q+j,j=0,…,2q−1j=0,\dots,2^q-1。它们按纵坐标排序,就是按 jj 的 qq 位倒序 rr 排序。倒序不改变 1 的个数,所以排在第 rr 位的点的颜色是 ±(−1)b(r)\pm(-1)^{b(r)},正负号对整个条带是同一个。
  3. 前缀和不超过 1。序列 (−1)b(r)(-1)^{b(r)},r=0,1,2,…r=0,1,2,\dots,是 +,−,−,+,−,+,+,−,…+,-,-,+,-,+,+,-,\dots。b(2r+1)=b(2r)+1b(2r+1)=b(2r)+1,所以它每两项一组,每组是 (+,−)(+,-) 或 (−,+)(-,+)。任何前缀和只能是 −1-1、00 或 11。
  4. 合计。条带 I×[0,Y)I\times[0,Y) 里的点是这个排序的一段前缀,颜色和的绝对值不超过 1。像第 6 节那样把锚定矩形拆成至多 ss 个条带加一段零头,颜色和的绝对值不超过 s+1s+1。一般矩形再乘 4。

所以,只要证明针对的是这样的点集,结论就不可能超过 log⁡n\log n。要突破它,必须换一类点集。

分解范数 γ2\gamma_2

第二类方法来自线性代数。把“哪个点在哪个矩形里”写成 0–1 矩阵 AA。把 AA 分解成两个矩阵的乘积 A=BCA=BC,定义

γ2(A)=min⁡A=BC (B 的最长行向量的长度)×(C 的最长列向量的长度).\gamma_2(A)=\min_{A=BC}\ \bigl(B\ \text{的最长行向量的长度}\bigr)\times\bigl(C\ \text{的最长列向量的长度}\bigr).

Matoušek、Nikolov 和 Talwar 证明,这个量从两边夹住了差异度(准确地说是“遗传差异度”:在所有子点集上取最坏的差异度。Tusnády 问题可以化成网格上的这个量):

γ2(A)Clog⁡M≤herdisc⁡(A)≤Clog⁡M γ2(A).\frac{\gamma_2(A)}{C\log M}\le\operatorname{herdisc}(A)\le C\sqrt{\log M}\,\gamma_2(A).

其中 MM 是矩阵行数,CC 是常数。第 5 节的二进分解给出一维前缀矩阵的一个分解:最大行范数与最大列范数都至多为 h+1\sqrt{h+1},所以它们的乘积至多为 h+1h+1。在含 nn 个点的各坐标方向等长的完整 dd 维乘积网格上,锚定盒子的关联矩阵满足 γ2=Θd(log⁡dn)\gamma_2=\Theta_d(\log^d n)。下面的下界通过遗传差异度选取网格的子集,并不声称每个点集都有这个 γ2\gamma_2 值。对任意 nn 点集,把各坐标的秩嵌入网格,则得到所需的 γ2\gamma_2 上界。因此,固定 dd 时,这两个通用不等式给出

Ω(log⁡d−1n) ≤ Δd(n) ≤ O(log⁡d+1/2n).\Omega(\log^{d-1}n)\ \le\ \Delta_d(n)\ \le\ O(\log^{d+1/2}n).
天花板

左右两个不等式各带一个对数因子的损失,而且已经知道这两个损失对一般的矩阵都无法改进。对平面乘积网格,γ2≈log⁡2n\gamma_2\approx\log^2n,下界一侧除以 log⁡n\log n 之后正好又是 log⁡n\log n。真实值 log⁡3/2n\log^{3/2}n 严格落在两个界之间,这把通用的尺子量不出来。

小结:三种方法,同一个数

方法给出的下界停在这里的原因
连续差异度 + 转移原理Ω(log⁡n)\Omega(\log n)连续差异度本身就是 Θ(log⁡n)\Theta(\log n)
正交函数法Ω(log⁡n)\Omega(\log n)已有的论证每层只挣到常数;第 6 节的二进网本身就能染到 O(log⁡n)O(\log n)
分解范数 γ2\gamma_2Ω(log⁡d−1n)\Omega(\log^{d-1}n),平面即 log⁡n\log n通用不等式自带对数因子的损失
本文:随机点集,逐个染色估计概率Ω(log⁡3/2n)\Omega(\log^{3/2}n)与上界吻合

这里回顾的具体论证有两个特点。第一,这些论证没有给出本文对随机点集的统一下界:转移原理和分解范数路线给出的坏点集藏在某个网格的子集里;前面介绍的正交函数法在组合差异度上的应用针对二进网,但该方法本身并不限于这类点集。第二,这里介绍的几种按尺度分层计算收益的论证,每层挣到的都是常数,log⁡n\log n 层合计 log⁡n\log n。本文两点都换掉了:点是随机的,每层尺度挣到的是 log⁡n\sqrt{\log n} 的量级。证明的骨架仍然是一层一层地挣,这一点和正交函数法一脉相承。不同的是每一层挣多少,以及怎样保证对全部 2n2^n 种染色同时成立。

Part III · 第 8–10 节

换一个问法

不去构造坏点集,而是随机撒点,再证明每一种染色都几乎不可能染好。这一部分讲这个转换、它对概率估计提出的要求,以及后面一直要用的一个小工具:把随机数一位一位地揭示出来。

§ 8 · 入门随机撒点,一次对付一种染色

第 3 节强调过,下界要排除全部 2n2^n 种染色。上面介绍的经典论证没有给出本文这种对随机点集、同时适用于全部染色的下界。论文换了一种安排:

点集不构造。在单位正方形 [0,1]2[0,1]^2 里独立、均匀地撒 nn 个点 P1,…,PnP_1,\dots,P_n,记 Pi=(Ui,Vi)P_i=(U_i,V_i)。
染色给第 ii 个点指定颜色 χi∈{+1,−1}\chi_i\in\{+1,-1\}。一种染色就是一个向量 χ=(χ1,…,χn)\chi=(\chi_1,\dots,\chi_n),一共 2n2^n 个。颜色是按点的编号指定的,撒点之前就定好了。
目标对每一个固定的 χ\chi,证明“这种染色在所有锚定矩形上都很平衡”的概率小得出奇。

把锚定矩形 [0,x]×[0,y][0,x]\times[0,y] 里的颜色和写成一个二元函数:

Fχ(x,y)=∑i=1nχi 1{Ui≤x, Vi≤y}.F_\chi(x,y)=\sum_{i=1}^n\chi_i\,\one_{\{U_i\le x,\ V_i\le y\}}.

这里 1{⋯ }\one_{\{\cdots\}} 是示性函数:条件成立时取 1,否则取 0。再记 ∥Fχ∥∞=max⁡0≤x,y≤1∣Fχ(x,y)∣\norm{F_\chi}_\infty=\max_{0\le x,y\le1}\abs{F_\chi(x,y)},它就是这种染色在锚定矩形上的最大差。论文的主定理是下面这句话。

定理 1.1(随机点)论文 Theorem 1.1

对每个 A>0A>0,存在常数 cA>0c_A>0,使下面的结论对每个整数 n≥1n\ge1 成立。设 P1,…,PnP_1,\dots,P_n 是 [0,1]2[0,1]^2 中独立均匀的随机点,则

Pr⁡[ 对每一个 χ∈{±1}n 都有 ∥Fχ∥∞≥cA(log⁡2n)3/2 ] ≥ 1−e−An.\Pr\Bigl[\ \text{对每一个}\ \chi\in\{\pm1\}^n\ \text{都有}\ \norm{F_\chi}_\infty\ge c_A(\log_2n)^{3/2}\ \Bigr]\ \ge\ 1-e^{-An}.

方括号里的事件说的是点集,不是某一种染色:这组点不管怎么染,都有一个锚定矩形的差至少是 cA(log⁡2n)3/2c_A(\log_2n)^{3/2}。所以染色完全可以在看到所有点之后再挑。你挑出来的那一种也是 2n2^n 个向量之一,结论对它照样成立。

从这个定理马上得到三件事。

  • 坏点集存在。一个概率为正的事件不是空的,而随机撒的点以概率 1 互不相同,所以至少有一个 nn 点集,任何染色的差都至少是 cA(log⁡2n)3/2c_A(\log_2n)^{3/2}。锚定矩形也是矩形,于是 Δ2(n)≥clog⁡3/2n\Delta_2(n)\ge c\log^{3/2}n。加上 Nikolov 的上界,Δ2(n)=Θ(log⁡3/2n)\Delta_2(n)=\Theta(\log^{3/2}n)。
  • 随机点集就是最坏的。不需要任何巧妙的构造,随手撒的点以极接近 1 的概率达到最坏的量级。
  • AA 可以任意大。想让失败概率是 e−100ne^{-100n} 也行,只是常数 cAc_A 会跟着变小。论文没有去优化 cAc_A,它非常小,所以这个界只对很大的 nn 有意义;第 18 节末尾给出了具体的数。

证明的两步

记 t=cA(log⁡2n)3/2t=c_A(\log_2n)^{3/2}。证明分成两步,第一步是全文的主体。

固定染色对每一个固定的 χ\chi,证明 Pr⁡[∥Fχ∥∞≤t]≤e−(A+ln⁡2)n\Pr\bigl[\norm{F_\chi}_\infty\le t\bigr]\le e^{-(A+\ln2)n}。
并集界“存在某个 χ\chi 染得好”是 2n2^n 个事件的并,所以它的概率至多是 2n⋅e−(A+ln⁡2)n=e−An2^n\cdot e^{-(A+\ln2)n}=e^{-An}。

第二步只用到第 5 节的并集界,不需要这 2n2^n 个事件之间有任何独立性。但它对第一步提了一个很高的要求:每一种染色“染得好”的概率都必须比 2−n2^{-n} 还小。如果只证到 0.9n0.9^n,乘上 2n2^n 之后是 1.8n1.8^n,什么也说明不了。换句话说,平均每多一个点,“染得好”的概率就得不止减半。

第四部分的证明比这两步稍微精细一点。“染得好”会被拆成两个事件:一个只和点的位置有关,对所有染色是同一个事件,排除一次就够了;另一个才需要逐个染色估计,再乘 2n2^n。但要求没有变:逐个染色估计的那一部分,概率必须比 2−n2^{-n} 小得多。

第 5 节的随机染色(上界)这里(下界)
随机的是什么染色点的位置
固定的是什么点集,和一个矩形一种染色
并集界对什么用至多 n4n^4 个矩形2n2^n 种染色
每个事件的概率要多小n−5n^{-5} 就够必须小于 2−n2^{-n}
结论存在一种好染色以极高的概率没有好染色

固定 χ\chi 之后,随机性全在点的位置上:每个点带着事先指定的颜色,独立地落在正方形里均匀的一处。另外,同一个点的横坐标 UiU_i 和纵坐标 ViV_i 也是独立的。第 10 节会利用这一点:先把全部纵坐标看作已知,只让横坐标随机。

§ 9 · 进阶该瞄准多大:32\frac32 是怎么猜出来的

动手证明之前,先估一下 tt 最大能取到多少。这一节是启发式的计算,对应论文第 3 节里“The Brownian sheet and Talagrand's proof”那一段的前半。它不是证明的一部分,但能解释 32\frac32 这个指数的来历。

固定 χ\chi,记 S=∑iχiS=\sum_i\chi_i。Fχ(x,y)F_\chi(x,y) 是 nn 个独立项之和,每一项以概率 xyxy 贡献 ±1\pm1,否则贡献 0。它的均值是 S xyS\,xy,围绕均值的起伏是 n\sqrt n 的量级。减去均值,再除以 n\sqrt n,得到的过程协方差是 min⁡(x,x′)min⁡(y,y′)−xy x′y′\min(x,x')\min(y,y')-xy\,x'y'。概率论里有一张随机曲面和它很像,叫布朗单,记作 W(x,y)W(x,y)。它是布朗运动的二维版本:W(x,y)W(x,y) 服从均值 0、方差 xyxy 的正态分布,协方差是 min⁡(x,x′)min⁡(y,y′)\min(x,x')\min(y,y'),不相交的矩形上的增量相互独立。两者差的那一项 xy x′y′xy\,x'y',相当于把布朗单在右上角钉住,换成 W(x,y)−xy W(1,1)W(x,y)-xy\,W(1,1)。下面只借布朗单来猜尺度,不区分这两个过程,也不管均值那一项。

小球估计背景 · Bass 1988,Talagrand 1994

布朗单在整个正方形上都留在 [−ε,ε][-\varepsilon,\varepsilon] 里的概率满足

−log⁡Pr⁡[∥W∥∞≤ε] ≍ ε−2log⁡3(1/ε)(ε→0),-\log\Pr\bigl[\norm{W}_\infty\le\varepsilon\bigr]\ \asymp\ \varepsilon^{-2}\log^3(1/\varepsilon)\qquad(\varepsilon\to0),

其中 ≍\asymp 表示两边相差不超过常数倍。左边的上界是 Bass 证明的,与之匹配的下界是 Talagrand 证明的。后面还会用到 ≳\gtrsim 和 ≲\lesssim:不计常数倍的“至少”和“至多”。

一维的对照更容易理解。布朗运动 B(s)B(s) 在 [0,1][0,1] 上始终留在 [−ε,ε][-\varepsilon,\varepsilon] 里的概率大约是 e−c/ε2e^{-c/\varepsilon^2}。理由是:把时间切成长度 ε2\varepsilon^2 的小段,每一小段里它走出 ±ε\pm\varepsilon 的概率是个常数,而要连续 1/ε21/\varepsilon^2 段都不出界。二维多出来的 log⁡3(1/ε)\log^3(1/\varepsilon) 没有这么简单的解释,这正是 Talagrand 定理的内容。

现在把它代进第 8 节的要求里。“∥Fχ∥∞≤t\norm{F_\chi}_\infty\le t”相当于“∥W∥∞≤ε\norm{W}_\infty\le\varepsilon”,其中 ε=t/n\varepsilon=t/\sqrt n。于是猜测

Pr⁡[∥Fχ∥∞≤t] ≈ exp⁡(−c nt2 log⁡3nt).\Pr\bigl[\norm{F_\chi}_\infty\le t\bigr]\ \approx\ \exp\Bigl(-c\,\frac{n}{t^2}\,\log^3\frac{\sqrt n}{t}\Bigr).

要让它小于 2−n2^{-n},指数必须至少是 nn 的量级:

nt2 log⁡3nt ≳ n⟺t ≲ log⁡3/2nt.\frac{n}{t^2}\,\log^3\frac{\sqrt n}{t}\ \gtrsim\ n\qquad\Longleftrightarrow\qquad t\ \lesssim\ \log^{3/2}\frac{\sqrt n}{t}.

当 tt 只是 log⁡n\log n 的幂次时,log⁡nt≈12log⁡n\log\frac{\sqrt n}t\approx\frac12\log n。所以门槛正好在 t≍log⁡3/2nt\asymp\log^{3/2}n:比它小的 tt,并集界够用;比它大的,就不够了。32\frac32 就是 3 的一半:小球估计里对数的三次方,开平方之后变成 32\frac32 次方。

问题对应的随机过程小球概率的指数解出的 tt实际答案
直线上的区间布朗运动ε−2\varepsilon^{-2}常数1(第 2 节)
平面上的矩形布朗单ε−2log⁡3(1/ε)\varepsilon^{-2}\log^3(1/\varepsilon)log⁡3/2n\log^{3/2}nΘ(log⁡3/2n)\Theta(\log^{3/2}n)(本文)

第一行是个验算:一维时指数里没有对数,解出来 tt 是常数,和第 2 节“交替染色,差恒为 1”对得上。

这只是猜测

“像布朗单”说的是分布的极限,它描述的是典型的起伏。这里要估计的概率是 e−ne^{-n} 这么小的数,弱收敛本身控制不了指数小的概率。论文的说法是:这个计算只预言了尺度;FχF_\chi 不是高斯过程,高斯极限控制不了这么小的概率。所以第四部分要直接对有限个点证明。

另外,上面的类比略去了两样东西:极限过程在右上角是钉住的;颜色不平衡时 FχF_\chi 还有均值 S xyS\,xy,全染成红色时 Fχ(1,1)=nF_\chi(1,1)=n,根本不是 n\sqrt n 的量级。猜尺度时可以先不管它们,证明里不能。

第四部分的证明完全不用高斯分布。但它最后得到的参数关系,和上面这个计算一模一样,第 19 节会回来对一下账。

§ 10 · 进阶逐位揭示:把一个随机数拆成一串骰子

证明里的随机性是一点一点放出来的。先看一个随机数怎么拆。取整数 b≥2b\ge2 作为底数。把 [0,1][0,1] 等分成 bb 段,UU 落在第几段,就把这个编号记作它的第 1 位数字 d1∈{1,…,b}d_1\in\{1,\dots,b\}。再把那一段等分成 bb 小段,UU 落在第几小段记作 d2d_2。这样一直做下去。这就是 UU 的 bb 进制小数,只是数字从 1 数起,为的是和后面“第 rr 个子区间”对上号。

01 1234 1234 1234 底数 b = 4:这个点的前三位数字是 3、2、4
每往下一行,就把上一行高亮的那一段放大成整行,再四等分。
工具 · 均匀随机数的各位数字

如果 UU 在 [0,1][0,1] 上均匀分布,那么它的各位数字 d1,d2,…d_1,d_2,\dots 相互独立,而且每一位都在 {1,…,b}\{1,\dots,b\} 上均匀分布。

证明:“前 jj 位数字依次是 r1,…,rjr_1,\dots,r_j”等价于“UU 落在某一个长度为 b−jb^{-j} 的区间里”,概率是 b−j=1b⋅1b⋯1bb^{-j}=\frac1b\cdot\frac1b\cdots\frac1b。联合概率等于各自概率的乘积,这正是独立且均匀的定义。

所以一个均匀随机数可以看成一串独立的 bb 面骰子。已经掷出前 jj 个,只知道 UU 在哪一个长 b−jb^{-j} 的区间里;下一个骰子决定它落在这个区间的哪一个子区间,而且与前面掷出的结果无关。

已知信息和条件期望

证明分阶段进行,阶段从 0 编到 hh(hh 是层数,第 11 节取定)。记 Fj\mathcal F_j 为第 jj 阶段的已知信息:全部 nn 个点的纵坐标 V1,…,VnV_1,\dots,V_n,加上每个横坐标 UiU_i 的前 jj 位数字。F0\mathcal F_0 只知道纵坐标。每往下走一个阶段,就多揭示每个点的一位数字。

因为各点相互独立,UiU_i 又与 ViV_i 独立,所以在 Fj\mathcal F_j 已知的条件下,nn 个点的下一位数字是 nn 个独立的、均匀的 bb 面骰子。从第 jj 阶段到第 j+1j+1 阶段,新增的随机性只有这些。

工具 · 条件期望

E[X∣Fj]\EE[X\mid\mathcal F_j] 的意思是:把已知信息当成常数,只对还没揭示的数字求平均。结果是一个由已知信息决定的量。后面用到它的四条性质。

线性。E[X+Y∣Fj]=E[X∣Fj]+E[Y∣Fj]\EE[X+Y\mid\mathcal F_j]=\EE[X\mid\mathcal F_j]+\EE[Y\mid\mathcal F_j],不需要独立性。

已知的量可以提出来。如果 YY 由 Fj\mathcal F_j 决定,那么 E[YX∣Fj]=Y⋅E[X∣Fj]\EE[YX\mid\mathcal F_j]=Y\cdot\EE[X\mid\mathcal F_j]。求平均时 YY 是常数。

逐层平均。先对后揭示的数字求平均,再对先揭示的求平均,等于一次对全部求平均:E[E[X∣Fj+1]∣Fj]=E[X∣Fj]\EE\bigl[\EE[X\mid\mathcal F_{j+1}]\bigm|\mathcal F_j\bigr]=\EE[X\mid\mathcal F_j],特别地 E[E[X∣Fj]]=E[X]\EE\bigl[\EE[X\mid\mathcal F_j]\bigr]=\EE[X]。

绝对值。∣E[X∣Fj]∣≤E[ ∣X∣∣Fj]\bigl|\EE[X\mid\mathcal F_j]\bigr|\le\EE\bigl[\,|X|\bigm|\mathcal F_j\bigr]。平均值的绝对值不超过绝对值的平均值,这是三角不等式。

如果只有有限种可能的结果,条件期望就是一个加权平均,这四条都可以直接验证。本文的情形正是如此:每一步只涉及 nn 个骰子,一共 bnb^n 种结果,每种的概率都是 b−nb^{-n}。只有第 13 节有一处要对连续的位置求平均,那里会说明。

Part IV · 第 11–18 节

局部振荡

这一部分是证明本身,依据论文的 r65 稿。主角是一个叫局部振荡的量 ZjZ_j。证明的骨架是三件事:过一层,ZZ 平均涨一截(引理 2.1、2.3、2.4 和 Fact 2.2);点摊得开的时候,这一截足够大,而点几乎总是摊得开(引理 2.5);实际的涨幅很难比平均值低太多(引理 2.6)。命题 2.7 把它们合起来,对全部染色用并集界。

路线是:先把正方形按阶段切成小矩形,定义 ZjZ_j(第 11、12 节)。在一个矩形里证明,细分之后振荡的期望只增不减,而且多出一截(第 13、14 节);一位数字对 ZZ 的影响很小(第 15 节)。然后把 hh 层接起来(第 16 节),说明在 M≤nM\le n 时,各层以高概率同时是好层(第 17 节),最后选参数收尾(第 18 节)。公式编号 (1)–(18)、引理 2.1、2.3–2.6、Fact 2.2、命题 2.7 和附录的 Fact A.1 都和论文一致。个别地方为了讲解调了先后,会在当处说明。

§ 11 · 入门h+1h+1 个阶段,一族等面积的矩形

从这里到第 17 节,染色 χ\chi 是固定的,把 FχF_\chi 简写成 FF。以概率 1,各点的纵坐标互不相同,也没有哪个点的坐标形如 k/brk/b^r(kk、r≥0r\ge0 是整数)。下面只在这个事件上讨论。后面所有区间的端点都是 k/brk/b^r 的形式,所以没有点落在分割线上。

取三个参数:

  • 底数 b≥3b\ge3,是整数;
  • 层数 h≥1h\ge1,是整数;
  • M=bh−1M=b^{h-1}。

它们的具体取值到第 18 节才定:bb 是一个只依赖 AA 的大常数,hh 大约是 log⁡n\log n 除以 log⁡b\log b。

阶段。对 j=0,1,…,hj=0,1,\dots,h,第 jj 阶段把正方形这样切开:

  • 横向切成 bjb^j 个长度为 b−jb^{-j} 的区间 II。前 jj 位数字相同的横坐标正好落在同一个区间里。
  • 纵向切成 bh−jb^{h-j} 个高度为 bj−hb^{j-h} 的块 BB。块含下端、不含上端。

于是第 jj 阶段共有 bj⋅bh−j=bh=bMb^j\cdot b^{h-j}=b^h=bM 个小矩形 I×BI\times B,每个的面积都是 b−j⋅bj−h=b−hb^{-j}\cdot b^{j-h}=b^{-h}。个数和面积都与 jj 无关,变的只是形状:第 0 阶段是 bhb^h 根横放的细条,第 hh 阶段是 bhb^h 根竖放的细条。

交互 · 06同一个正方形的四个阶段

为了看得清,这里取 b=4b=4、h=3h=3。拖动滑块改变阶段 jj:横向越分越细,纵向越并越粗,小矩形的个数和面积始终不变。下面的读数里写着圈出的那个点已经揭示的数字。

阶段 j=0,1,…,hj=0,1,\dots,h 一共 h+1h+1 个。从第 jj 阶段到第 j+1j+1 阶段的过渡,叫做第 jj 层,j=0,…,h−1j=0,\dots,h-1,一共 hh 层。论文把它叫做 transition j→j+1j\to j+1,把小矩形叫做 cell,把子区间叫做 child,把纵向的区间按场合叫做 vertical interval 或 block。

过一层,同时发生两件事:

  • 每个区间 II 被切成 bb 个子区间 C1,…,CbC_1,\dots,C_b,从左到右编号。一个点落在哪个子区间,由它的下一位数字决定。
  • 每 bb 个相邻的旧块 B1,…,BbB_1,\dots,B_b(从下到上)并成一个新块 JJ。
B1B2⋮⋮Bb C1C2⋯⋯Cb J II 第 j 阶段:b 个旧块第 j+1 阶段:b 个子区间 纵向合并横向细分
一次过渡发生在一个矩形 R=I×JR=I\times J 里(对应论文的 Figure 2 和 Figure 3 左)。左边的 bb 个旧矩形 I×BℓI\times B_\ell 和右边的 bb 个新矩形 Ck×JC_k\times J 铺满的是同一块地方。

这样的 R=I×JR=I\times J 叫比较矩形:II 是第 jj 阶段的区间,JJ 是第 j+1j+1 阶段的块。它一共有 bj⋅bh−j−1=Mb^j\cdot b^{h-j-1}=M 个,每个面积是 1M\frac1M。所以

nM=一个比较矩形里的平均点数.\frac nM=\text{一个比较矩形里的平均点数}.

这个比值后面会反复出现。到第 18 节会看到,它应该取成 hh 的量级。

“面积相同、形状各异的一族矩形”并不新鲜:第 7 节 Roth 的方法用的就是面积固定的二进盒子。新的地方在于,下面要在这族矩形上量一个不同的东西。

§ 12 · 进阶端点平均与振荡

端点平均

对一个横向区间 I=[s,t]I=[s,t],把 FF 在两个端点处的值取平均:

aI(y)=F(s,y)+F(t,y)2.(2)a_I(y)=\frac{F(s,y)+F(t,y)}2.\tag{2}

论文把端点记作 I−I^-、I+I^+,把这个函数记作 a(I)a^{(I)},在 II 固定的时候简写成 aa。式 (1) 是下面 ZjZ_j 的定义,论文在引言里就给出了。

它是 yy 的函数,叫做 II 的端点平均。看看每个点在里面占多大份量。在纵坐标不超过 yy 的点里,F(s,y)F(s,y) 把 II 左边的点按颜色加起来,F(t,y)F(t,y) 加的是这些点再加上 II 里面的点。所以:

左端点 s右端点 t区间 I 两个端点都数到它:系数 1 只有右端点数到:系数 ½ 都数不到:系数 0
一个点在 aIa_I 里的系数只取决于它在 II 的左边、里面还是右边。
aI(y)=∑i=1nχi wI(Ui) 1{Vi≤y},wI(u)={1,u 在 I 的左边,12,u 在 I 里,0,u 在 I 的右边.a_I(y)=\sum_{i=1}^n\chi_i\,w_I(U_i)\,\one_{\{V_i\le y\}},\qquad w_I(u)=\begin{cases}1,&u\ \text{在}\ I\ \text{的左边},\\[2pt]\tfrac12,&u\ \text{在}\ I\ \text{里},\\[2pt]0,&u\ \text{在}\ I\ \text{的右边}.\end{cases}

所以 aI(y)a_I(y) 是一个阶梯函数:yy 每往上越过一个点,它就跳 χiwI(Ui)\chi_iw_I(U_i),也就是 ±1\pm1、±12\pm\frac12 或者不动。那个 12\frac12 看起来有点怪,第 13、14 节会看到它是整个证明的支点。论文把这个系数叫做点的权重(weight)。

例

取 b=4b=4、h=1h=1。看区间 I=[14,12]I=[\frac14,\frac12],放三个点:II 的左边有一个红点,高度 0.3;II 里有一个蓝点,高度 0.55,还有一个红点,高度 0.8。它们的系数依次是 11、12\frac12、12\frac12。让 yy 从下往上走:

aI(y)={0,y<0.3,1,0.3≤y<0.55,1−12=12,0.55≤y<0.8,12+12=1,y≥0.8.a_I(y)=\begin{cases}0,&y<0.3,\\1,&0.3\le y<0.55,\\1-\frac12=\frac12,&0.55\le y<0.8,\\\frac12+\frac12=1,&y\ge0.8.\end{cases}

0.3 以下没有点;越过左边的红点,加 1;越过蓝点,减 12\frac12;越过最后一个红点,加 12\frac12。

交互 · 07端点平均 aI(y)a_I(y) 是怎么数出来的

取 b=4b=4,放 12 个点。左图把两个锚定矩形 [0,s]×[0,y][0,s]\times[0,y] 和 [0,t]×[0,y][0,t]\times[0,y] 叠在一起:II 左边的那一条被盖了两层,II 里面只盖了一层,所以系数是 1 和 12\frac12。右图把 aIa_I 画成高度的函数,纵轴和左图对齐。拖动高度 yy,看虚线每越过一个点,右边跳多少。

几件事值得动手试一下。把阶段调到 0,这时 I=[0,1]I=[0,1],左端点的矩形是空的,aI(y)=12F(1,y)a_I(y)=\frac12F(1,y),每个点的系数都是 12\frac12。把区间移到最右边,几乎所有点的系数都是 1,aIa_I 接近 F(1,y)F(1,y) 本身。翻转 II 里面一个点的颜色,aIa_I 在它上方的部分整体挪 1;翻转 II 左边的一个点,挪 2;翻转右边的点,什么也不变。

振荡

aI(y)a_I(y) 是阶梯函数,只跳有限次,所以在任何一个块上它都有最大值和最小值。对一个关于高度的函数 ff 和一个纵向的块 BB,记

osc⁡B f=max⁡Bf−min⁡Bf,mid⁡B f=max⁡Bf+min⁡Bf2,(3)\operatorname{osc}_B\,f=\max_Bf-\min_Bf,\qquad\operatorname{mid}_B\,f=\frac{\max_Bf+\min_Bf}2,\tag{3}

分别叫做 ff 在 BB 上的振荡和中点。振荡量的是 ff 在 BB 上取值的范围有多宽,中点量的是这个范围在哪里。最大值和最小值就是 mid⁡B f±12osc⁡B f\operatorname{mid}_B\,f\pm\frac12\operatorname{osc}_B\,f。块含下端、不含上端。

工具 · 振荡的三条性质

不怕加常数。整个块上同时加一个数 cc,最大值和最小值一起动,振荡不变;中点跟着加 cc。

三角不等式。osc⁡B(f+g)≤osc⁡B f+osc⁡B g\operatorname{osc}_B(f+g)\le\operatorname{osc}_B\,f+\operatorname{osc}_B\,g。理由:max⁡(f+g)≤max⁡f+max⁡g\max(f+g)\le\max f+\max g,min⁡(f+g)≥min⁡f+min⁡g\min(f+g)\ge\min f+\min g,两式相减。反复用它,对有限个函数的和同样成立。另外,乘一个常数 cc,振荡变成 ∣c∣|c| 倍:osc⁡B(cf)=∣c∣osc⁡B f\operatorname{osc}_B(cf)=|c|\operatorname{osc}_B\,f。

反向三角不等式。∣osc⁡B(f+g)−osc⁡B f∣≤osc⁡B g\bigl|\operatorname{osc}_B(f+g)-\operatorname{osc}_B\,f\bigr|\le\operatorname{osc}_B\,g。理由:把 ff 写成 (f+g)+(−g)(f+g)+(-g),再用一次三角不等式和 osc⁡B(−g)=osc⁡B g\operatorname{osc}_B(-g)=\operatorname{osc}_B\,g,得到 osc⁡B f≤osc⁡B(f+g)+osc⁡B g\operatorname{osc}_B\,f\le\operatorname{osc}_B(f+g)+\operatorname{osc}_B\,g。另一个方向就是三角不等式本身。

第一条是选振荡的原因:挪动一个高度不在 BB 里的点,BB 上的函数只是整体加了一个常数。第二条在第 13 节用,第三条在第 15 节用。

势函数 ZjZ_j

把第 jj 阶段所有小矩形上的振荡取平均:

Zj=1bM∑I×Bosc⁡B aI,(1)Z_j=\frac1{bM}\sum_{I\times B}\operatorname{osc}_B\,a_I,\tag{1}

求和遍历第 jj 阶段的 bMbM 个小矩形 I×BI\times B(论文叫 cell)。论文把 ZjZ_j 叫做势函数(potential)。在上面的例子里,II 是第 1 阶段的区间,唯一的块是整个 [0,1)[0,1),aIa_I 在它上面取到 00、12\frac12、11 三个值,振荡是 1−0=11-0=1。

ZjZ_j 的两条基本性质论文式 (9) 和它后面一段
0≤Zj≤2∥F∥∞,(9)0\le Z_j\le2\norm F_\infty,\tag{9}

并且 ZjZ_j 由 Fj\mathcal F_j 决定。

式 (9) 在论文里到 2.4 小节才出现,这里提前用了它的编号。

推导两条性质
  1. 下界。每一项都是最大值减最小值,不会是负的。
  2. 上界。由 (3),一个振荡是 aIa_I 的两个值之差;由 (2),aIa_I 的每个值是 FF 的两个值的平均,绝对值不超过 ∥F∥∞\norm F_\infty。所以每一项振荡都不超过 2∥F∥∞2\norm F_\infty,平均值也一样。
  3. 由 Fj\mathcal F_j 决定。第 jj 阶段的区间端点是 b−jb^{-j} 的整数倍。一个点在某个端点的左边还是右边,看它的前 jj 位数字就知道了。所以每个系数 wI(Ui)w_I(U_i) 都由 Fj\mathcal F_j 决定,纵坐标本来就已知。

注意 Z0Z_0 不一定是 0:第 0 阶段的块虽然很矮,里面仍然可能有点。后面只用到 Z0≥0Z_0\ge0。

式 (9) 的上界是 ZZ 和差异度之间唯一的联系,但已经够用了:

∥Fχ∥∞ ≥ 12Zh.\norm{F_\chi}_\infty\ \ge\ \tfrac12Z_h.

所以只要证明 ZhZ_h 很大,这种染色就有一个差很大的锚定矩形。而 Z0≥0Z_0\ge0,之后一共走 hh 步。接下来的计划就是证明:它每一步平均都在涨,而且实际值很难比平均值低太多。

交互 · 08ZjZ_j:全部小矩形上振荡的平均

仍然取 b=4b=4、h=3h=3,放 40 个点。每个阶段都有 bM=64bM=64 个小矩形 I×BI\times B。左图给每个小矩形上色,颜色越浓,aIa_I 在它上面的振荡越大;ZjZ_j 就是这 64 个振荡的平均。点一个小矩形,右图画出它所在区间的 aI(y)a_I(y),橙色的带子是选中的块,粗横线从块里的最小值画到最大值,长度就是振荡。

把阶段从 0 调到 3,可以看到三件事。第 0 阶段的块很矮,一个块里通常只有零到两三个点,振荡多半是 0、12\frac12 或 1,所以 Z0Z_0 很小。在这组点上,往后小矩形越来越瘦高,每个块越来越高,里面的点越来越多,颜色整体变浓,ZjZ_j 跟着变大。但不管怎么变,横条都不会越过虚线,这就是式 (9)。再换成“按 yy 顺序交替”染色:F(1,y)F(1,y) 被压得很小,可是在这个样本上 ZjZ_j 仍然涨得起来,因为窄区间上的 aIa_I 只看得见一部分点,交替的规律在那里被打乱了。

这些都是一组样本上的现象。单次实现里 ZjZ_j 完全可以下降,两个点就够:取 b=3b=3、h=1h=1,点 (45,15)(\frac45,\frac15) 染红、(910,45)(\frac9{10},\frac45) 染蓝,按定义算出 Z0=13Z_0=\frac13,Z1=16Z_1=\frac16。后面证明的是条件期望在涨,再另外控制实际值对条件期望的偏离。

为什么是端点平均,为什么是振荡

端点平均让“细分”变得干净:把 II 切成 bb 个子区间,再让每个点随机落进其中一个,子区间的端点平均再取平均,期望恰好等于原来的 aIa_I(第 14 节的式 (7))。只取右端点就没有这个性质。

振荡不怕加常数:整个块上同时加一个数,最大值和最小值一起动,差不变。所以一个点横向移动时,对每个固定的横向区间 I,只有包含该点高度的纵向块的振荡可能改变(第 15 节)。直接用 ∣aI∣|a_I| 的最大值就做不到这一点。

§ 13 · 核心推导一个矩形里的比较(上):不含概率的部分

这两节的目标

固定一个比较矩形 R=I×JR=I\times J。细分之前,RR 被切成 bb 个横条 I×BℓI\times B_\ell;细分之后,被切成 bb 个竖条 Ck×JC_k\times J。把 RR 最低和最高的两个横条去掉,剩下的部分叫做 RR 的内部区域,里面的点数记作 mRm_R:

Rin=I×(B2∪⋯∪Bb−1),mR=∣P∩Rin∣.R^{\mathrm{in}}=I\times(B_2\cup\dots\cup B_{b-1}),\qquad m_R=|P\cap R^{\mathrm{in}}|.

因为 b≥3b\ge3,内部区域不是空的。这两节要证明的是引理 2.3:竖条上的振荡之和,期望至少是横条上的振荡之和再加一笔,

E[ ∑k=1bosc⁡J aCk] ≥ ∑ℓ=1bosc⁡Bℓ aI + 14mR.\EE\Bigl[\,\sum_{k=1}^b\operatorname{osc}_J\,a_{C_k}\Bigr]\ \ge\ \sum_{\ell=1}^b\operatorname{osc}_{B_\ell}\,a_I\ +\ \frac14\sqrt{m_R}.

这一节先做不含概率的部分,得到式 (5) 和 (6)。它们对骰子的任何结果都成立。下一节再取期望。

为了少写下标,这两节把 aIa_I 简写成 aa,叫旧函数;把 aCka_{C_k} 简写成 aka_k,叫子函数。论文也是这样写的。

第一步:一个关于阶梯函数的不等式

先不管点和骰子,只看一个定义在 JJ 上、只取有限个值的函数 ff。JJ 上的振荡和 bb 个旧块上的振荡有什么关系?最容易想到的是:JJ 包含每一个旧块,所以 osc⁡J f≥osc⁡Bℓ f\operatorname{osc}_J\,f\ge\operatorname{osc}_{B_\ell}\,f,bb 个不等式加起来就是 bosc⁡J f≥∑ℓosc⁡Bℓ fb\operatorname{osc}_J\,f\ge\sum_\ell\operatorname{osc}_{B_\ell}\,f。引理 2.1 说,其实还能多出一项。

引理 2.1(纵向比较)论文 Lemma 2.1

对 JJ 上每一个只取有限个值的函数 ff,

b osc⁡J f ≥ ∑ℓ=1bosc⁡Bℓ f+2 ∣mid⁡Bbf−mid⁡B1f∣.(4)b\,\operatorname{osc}_J\,f\ \ge\ \sum_{\ell=1}^b\operatorname{osc}_{B_\ell}\,f+2\,\bigl|\operatorname{mid}_{B_b}f-\operatorname{mid}_{B_1}f\bigr|.\tag{4}

mid⁡Bbf−mid⁡B1f\operatorname{mid}_{B_b}f-\operatorname{mid}_{B_1}f 叫做最高的块相对最低的块的偏移(论文叫 offset):ff 在最高的块里取值的范围,和在最低的块里取值的范围,中心差了多少。

推导引理 2.1
  1. 首尾两个块合起来看。B1B_1 和 BbB_b 都在 JJ 里,而 ff 在一个块 BB 上的最大、最小值是 mid⁡B f±12osc⁡B f\operatorname{mid}_B\,f\pm\frac12\operatorname{osc}_B\,f。拿最高的块的最大值减最低的块的最小值,再反过来减一次: osc⁡J f ≥ max⁡Bbf−min⁡B1f=osc⁡B1f+osc⁡Bbf2+(mid⁡Bbf−mid⁡B1f),osc⁡J f ≥ max⁡B1f−min⁡Bbf=osc⁡B1f+osc⁡Bbf2−(mid⁡Bbf−mid⁡B1f).\begin{aligned}\operatorname{osc}_J\,f&\ \ge\ \max_{B_b}f-\min_{B_1}f=\frac{\operatorname{osc}_{B_1}f+\operatorname{osc}_{B_b}f}2+\bigl(\operatorname{mid}_{B_b}f-\operatorname{mid}_{B_1}f\bigr),\\[4pt]\operatorname{osc}_J\,f&\ \ge\ \max_{B_1}f-\min_{B_b}f=\frac{\operatorname{osc}_{B_1}f+\operatorname{osc}_{B_b}f}2-\bigl(\operatorname{mid}_{B_b}f-\operatorname{mid}_{B_1}f\bigr).\end{aligned}
  2. 取较大的一个。两个右边只差在偏移前面的符号上,所以 osc⁡J f ≥ osc⁡B1f+osc⁡Bbf2+∣mid⁡Bbf−mid⁡B1f∣.\operatorname{osc}_J\,f\ \ge\ \frac{\operatorname{osc}_{B_1}f+\operatorname{osc}_{B_b}f}2+\bigl|\operatorname{mid}_{B_b}f-\operatorname{mid}_{B_1}f\bigr|.
  3. 中间的块。对中间 b−2b-2 个块,各用一次 osc⁡J f≥osc⁡Bℓ f\operatorname{osc}_J\,f\ge\operatorname{osc}_{B_\ell}\,f。把这 b−2b-2 个不等式和第 2 步的两倍加起来,左边一共是 bb 个 osc⁡J f\operatorname{osc}_J\,f。
0.5 1.25 oscJ f B1 B2 B3 B4 b = 4 的示意。横轴是高度,从 J 的下端到上端;纵轴是 f 的值 粗竖线是 f 在首尾两个块里的取值范围,圆点是中点:偏移 = 1.25 − 0.5 = 0.75 J 上的振荡 2.5 ≥ B4 里的最大值 2 − B1 里的最小值 0
灰色的是最低和最高两个块。这个例子里四个块上的振荡是 1, 2, 1, 1.51,\ 2,\ 1,\ 1.5,偏移是 0.750.75,式 (4) 的右边是 5.5+2×0.75=75.5+2\times0.75=7,左边是 4×2.5=104\times2.5=10。

偏移就是后面增益的来源。给 ff 在最高的块上整体加一个常数:偏移变了,而每个块上的振荡都不变。所以式 (4) 里只有左边 JJ 上的振荡看得见它。

第二步:细分之后,对子函数取平均

现在把 II 切成 bb 个子区间 C1,…,CbC_1,\dots,C_b,每个子区间有自己的端点平均 ak=aCka_k=a_{C_k}。不去逐个研究它们,而是先把它们平均起来:

aˉ(y)=1b∑k=1bak(y).\bar a(y)=\frac1b\sum_{k=1}^ba_k(y).

由振荡的三角不等式(第 12 节),bb 个函数各自的振荡加起来,不小于它们的和的振荡;而它们的和是 b aˉb\,\bar a:

∑k=1bosc⁡J ak ≥ osc⁡J(∑k=1bak)=b osc⁡J aˉ.\sum_{k=1}^b\operatorname{osc}_J\,a_k\ \ge\ \operatorname{osc}_J\Bigl(\sum_{k=1}^ba_k\Bigr)=b\,\operatorname{osc}_J\,\bar a.

右边正好是引理 2.1 左边的形状。对 f=aˉf=\bar a 用引理 2.1:

先平均,再比较论文式 (5)
∑k=1bosc⁡J ak ≥ ∑ℓ=1bosc⁡Bℓ aˉ+2 ∣mid⁡Bbaˉ−mid⁡B1aˉ∣.(5)\sum_{k=1}^b\operatorname{osc}_J\,a_k\ \ge\ \sum_{\ell=1}^b\operatorname{osc}_{B_\ell}\,\bar a+2\,\bigl|\operatorname{mid}_{B_b}\bar a-\operatorname{mid}_{B_1}\bar a\bigr|.\tag{5}

左边是细分之后 bb 个竖条上的振荡之和,正是我们要的量。右边只剩一个函数 aˉ\bar a:它在 bb 个旧块上的振荡,加上它在首尾两个块之间的偏移的两倍。bb 个子区间各自的细节都不见了。

到这里为止没有用到任何概率。不管点落在哪个子区间,式 (5) 都成立。

第三步:一个点在 aˉ\bar a 里的权重

要比较 aˉ\bar a 和原来的 aa,先算每个点在 aˉ\bar a 里的权重。由式 (2),aˉ\bar a 数每个点时用的权重,是它在 bb 个子区间里的系数的平均;而 aa 用的是它在 II 里的系数,论文叫它的旧权重。II 左边的点,在每个 aka_k 里的系数都是 1,平均还是 1;II 右边的点都是 0。这两类点在 aˉ\bar a 里和在 aa 里一模一样。

II 里面的点旧权重是 12\frac12,新的权重要看它落在哪个子区间。设它落在 CrC_r 里:它右边的 b−rb-r 个子区间都把它算作“在左边”,各给系数 1;它自己所在的子区间给 12\frac12;它左边的子区间给 0。所以它在 aˉ\bar a 里的权重是

(b−r)+12b=12+b+1−2r2b(r=1,…,b).\frac{(b-r)+\frac12}b=\frac12+\frac{b+1-2r}{2b}\qquad(r=1,\dots,b).

这里的 rr 是子区间的编号,也就是骰子的结果。

b=4b=4 时,把每个子区间给它的系数排成一张表,行是子区间 kk,列是骰子的结果 rr;每一列的平均就是它在 aˉ\bar a 里的权重:

r=1r=1r=2r=2r=3r=3r=4r=4
a1a_1 里的系数12\frac12000
a2a_2 里的系数112\frac1200
a3a_3 里的系数1112\frac120
a4a_4 里的系数11112\frac12
aˉ\bar a 里的权重78\frac7858\frac5838\frac3818\frac18
比 12\frac12 多出的部分+38+\frac38+18+\frac18−18-\frac18−38-\frac38
交互 · 09一个点在 aˉ\bar a 里的权重

行是子区间 kk,列是骰子的结果 rr。先看一次骰子:点落进 CrC_r,表里的第 rr 列就是 bb 个子区间分别给它的系数,这一列的平均就是它在 aˉ\bar a 里的权重,写在表的下面。再连掷很多次,看权重的平均值往哪里去。

上面的式子把权重拆成了两部分。前一部分 12\frac12 正是旧权重。后一部分是多出来的偏差,可正可负。把颜色也乘进去,对落在 CrC_r 里的第 ii 个点记

Xi=χi⋅b+1−2r2b,X_i=\chi_i\cdot\frac{b+1-2r}{2b},

叫做它的增量(论文叫 jump:aˉ\bar a 在这个点的高度比 aa 多跳的那一下)。于是 aˉ\bar a 就是 aa 加上 II 里面各点的增量:

aˉ\bar a 等于 aa 加上增量论文式 (6)
aˉ(y)=a(y)+∑i: Ui∈I, Vi≤yXi.(6)\bar a(y)=a(y)+\sum_{i:\ U_i\in I,\ V_i\le y}X_i.\tag{6}

求和遍历 II 里面、高度不超过 yy 的点。用话说:aˉ\bar a 跟着旧函数走,只是在 II 里每个点的高度处多跳一下 XiX_i。这仍然是一个恒等式,没有概率。

交互 · 10平均函数 aˉ\bar a 和式 (5)

取 b=4b=4,画一个比较矩形 RR,并且只考虑 RR 里的点。先看“aa”:右图是阶梯函数 a(y)a(y),每个旧块旁边的粗线从块里的最小值画到最大值,长度就是这个块上的振荡。再看“aˉ\bar a”:实线是掷完骰子之后的 aˉ\bar a,虚线是原来的 aa;最低和最高两个块的粗线上各有一个圆点,标的是中点。图下面那条粗线从一个中点画到另一个中点,长度就是偏移的绝对值。点 C1,…,C4C_1,\dots,C_4 可以看每个子区间各自的函数和它在整个 JJ 上的振荡。

多掷几次骰子。aˉ\bar a 每次都不一样,在虚线 aa 的两侧摆动。而读数里的那串不等式每一次都成立:四个竖条上的振荡之和,不小于 4osc⁡Jaˉ4\operatorname{osc}_J\bar a,后者又不小于式 (5) 的右边。

易错

aˉ(y)\bar a(y) 本身是随机的,它依赖骰子。下一节会说它的期望是 a(y)a(y),但掷一次骰子,一个点的权重可以是 78\frac78,也可以是 18\frac18;b=4b=4 时没有哪个结果恰好给出 12\frac12(bb 是奇数时,只有正中间的子区间给出 12\frac12)。

那个 12\frac12 也不能换。如果不取端点平均,只用区间右端点处的 FF 值,点在自己所在的子区间里系数就是 1,权重变成 b−r+1b\frac{b-r+1}b,对 rr 平均是 b+12b\frac{b+1}{2b},而它在 II 里的系数是 1,两边对不上。

另一种看法:条件期望

在区间 II 里均匀地随机取一个横坐标 xx,与所有的点无关。那么 aI(y)a_I(y) 是 F(x,y)F(x,y) 在 Fj\mathcal F_j 下的条件期望:II 里面的点和 xx 都在 II 里均匀分布,谁在左边各占一半。论文第 3 节把它写成 aI(y)=E[1∣I∣∫IF(x,y) dx∣Fj]a_I(y)=\EE\bigl[\frac1{|I|}\int_IF(x,y)\,dx\bigm|\mathcal F_j\bigr],并指出这是 Talagrand 证明里“横向取平均”的条件期望版本。同样的道理,aˉ(y)\bar a(y) 是同一个量在 Fj+1\mathcal F_{j+1} 下的条件期望,只不过多知道了一位数字。(这里的平均包括各点在区间里的具体位置,不再是有限个结果的加权平均,但第 10 节的四条性质照样成立。)

§ 14 · 核心推导一个矩形里的比较(下):取期望,基准和增益

上一节得到了式 (5),它对骰子的任何结果都成立:

∑k=1bosc⁡J ak ≥ ∑ℓ=1bosc⁡Bℓ aˉ+2 ∣mid⁡Bbaˉ−mid⁡B1aˉ∣.\sum_{k=1}^b\operatorname{osc}_J\,a_k\ \ge\ \sum_{\ell=1}^b\operatorname{osc}_{B_\ell}\,\bar a+2\,\bigl|\operatorname{mid}_{B_b}\bar a-\operatorname{mid}_{B_1}\bar a\bigr|.

这一节对右边取期望,分两件事做。第一件:aˉ\bar a 在各旧块上的振荡,期望不小于 aa 的振荡,这是基准。第二件:偏移的绝对值,期望至少是 18mR\frac18\sqrt{m_R},乘以 2 就是增益。合起来就是引理 2.3。

随机的是什么

先说清楚什么是已知的,什么是随机的。各点的颜色和高度固定;每个点相对 II 的位置(在 II 的左边、在 II 里、在 II 的右边)也固定。II 里面的每个点独立地、等可能地落进 bb 个子区间之一,也就是各掷一个 bb 面的骰子。论文把这叫做局部实验(local experiment)。在局部实验里,旧函数 aa 是固定的,子函数 a1,…,aba_1,\dots,a_b 是随机的。下面的 E\EE 都是对这些骰子取的。

别的都不需要指定:II 外面的点在每个子区间的同一侧,所以子函数依赖横坐标的方式只有两样:每个点相对 II 的位置,和 II 里的点各自落进哪个子区间。第 16 节会说明,在 Fj\mathcal F_j 已知的条件下,过一层时每个比较矩形里发生的正是这个实验。

第一步:增量是对称的,所以 E aˉ=a\EE\,\bar a=a

看一个 II 里面的点的增量 Xi=χi⋅b+1−2r2bX_i=\chi_i\cdot\frac{b+1-2r}{2b},其中 rr 在 1,…,b1,\dots,b 上等可能。rr 从 1 走到 bb 时,b+1−2rb+1-2r 依次取 b−1, b−3, …, −(b−1)b-1,\ b-3,\ \dots,\ -(b-1),这组数关于 0 对称。所以不管颜色是什么,XiX_i 取每个值和取它的相反数的概率一样大:它是对称的(和 −Xi-X_i 同分布),特别地 EXi=0\EE X_i=0。不同的点用不同的骰子,各个增量相互独立。

交互 · 11一个点的增量:对称,均值为 0

条形是一个红点(χi=+1\chi_i=+1)落在各个子区间时的增量 b+1−2r2b\frac{b+1-2r}{2b}。左右对调子区间,条形上下翻转。读数先验证均值是 0,再给出后面要用的三个量:二阶矩 σ2\sigma^2、增量的最大绝对值、四阶矩。换一个 bb 看看它们怎么变。

由式 (6),aˉ(y)\bar a(y) 等于 a(y)a(y) 加上若干个增量,每个增量的期望是 0。所以

平均函数的期望是旧函数论文式 (7)
E aˉ(y)=a(y)对每个高度 y.(7)\EE\,\bar a(y)=a(y)\qquad\text{对每个高度 }y.\tag{7}

细分一步之后再取平均,期望意义下回到原来的函数。这就是第 12 节说的“端点平均让细分变得干净”。

第二步:基准

工具 · 最大值的期望不小于期望的最大值

设 ϕ(y)\phi(y) 是一个随机的函数(每个 yy 处的值都是随机变量),BB 是一个集合。对 BB 里任意一个固定的 y0y_0,总有 max⁡y∈Bϕ(y)≥ϕ(y0)\max_{y\in B}\phi(y)\ge\phi(y_0)。两边取期望:Emax⁡y∈Bϕ(y)≥E ϕ(y0)\EE\max_{y\in B}\phi(y)\ge\EE\,\phi(y_0)。这对每个 y0y_0 都成立,所以

Emax⁡y∈Bϕ(y) ≥ max⁡y0∈BE ϕ(y0).\EE\max_{y\in B}\phi(y)\ \ge\ \max_{y_0\in B}\EE\,\phi(y_0).

同样的道理,Emin⁡y∈Bϕ(y)≤min⁡y0∈BE ϕ(y0)\EE\min_{y\in B}\phi(y)\le\min_{y_0\in B}\EE\,\phi(y_0)。

直观地说:先看骰子再挑最好的高度,平均起来不会比事先挑定一个高度差。把工具用在 aˉ\bar a 和每个旧块 BℓB_\ell 上,再用式 (7):

Eosc⁡Bℓaˉ=Emax⁡Bℓaˉ−Emin⁡Bℓaˉ ≥ max⁡Bℓa−min⁡Bℓa=osc⁡Bℓa.\EE\operatorname{osc}_{B_\ell}\bar a=\EE\max_{B_\ell}\bar a-\EE\min_{B_\ell}\bar a\ \ge\ \max_{B_\ell}a-\min_{B_\ell}a=\operatorname{osc}_{B_\ell}a.

旧的振荡保得住。对 ℓ=1,…,b\ell=1,\dots,b 求和,再对式 (5) 两边取期望:

基准,加上偏移的期望论文式 (8)
E[ ∑k=1bosc⁡J ak] ≥ ∑ℓ=1bosc⁡Bℓ a+2 E ∣mid⁡Bbaˉ−mid⁡B1aˉ∣.(8)\EE\Bigl[\,\sum_{k=1}^b\operatorname{osc}_J\,a_k\Bigr]\ \ge\ \sum_{\ell=1}^b\operatorname{osc}_{B_\ell}\,a+2\,\EE\,\bigl|\operatorname{mid}_{B_b}\bar a-\operatorname{mid}_{B_1}\bar a\bigr|.\tag{8}

右边第一项是细分之前横条上的振荡之和,这一步说的是“不亏”。第二项是增益,下面三步给它一个下界。

第三步:把内部点的贡献单独拿出来

内部区域 Rin=I×(B2∪⋯∪Bb−1)R^{\mathrm{in}}=I\times(B_2\cup\dots\cup B_{b-1}) 里的点有一个特别的性质:它们都在 B1B_1 的上方、BbB_b 的下方。记

D=∑i: Pi∈RinXiD=\sum_{i:\ P_i\in R^{\mathrm{in}}}X_i

为这 mRm_R 个内部点的增量之和,再记 gg 为式 (6) 里去掉这些增量之后的函数:

g(y)=a(y)+∑i: Ui∈I, Pi∉Rin, Vi≤yXi.g(y)=a(y)+\sum_{i:\ U_i\in I,\ P_i\notin R^{\mathrm{in}},\ V_i\le y}X_i.

那么:

  • 在 B1B_1 上,内部点都还没被数到,所以 aˉ=g\bar a=g;
  • 在 BbB_b 上,内部点全都被数到了,所以 aˉ=g+D\bar a=g+D,比 gg 整体高出一个常数 DD。
B1B2B3B4B5B6B7B8 最高的块上:ā = g + D 最低的块上:ā = g 内部区域 Rin B2 到 B7。这里的点的 增量加起来是 D b = 8 的示意:这个 R 里有 14 个点,内部区域里有 10 个,mR = 10
对应论文的 Figure 3(左、中)。内部区域里的 10 个点决定 DD。上下两个灰色块里的 4 个淡色的点只进入 gg:它们的高度在 B1B_1 或 B8B_8 里,增量只影响块的一部分,不能当成常数提出来。

整个块上加一个常数 DD,最大值和最小值都跟着加 DD,中点也加 DD。所以

偏移里有一个 DD论文引理 2.3 的证明,Step 4
mid⁡Bbaˉ−mid⁡B1aˉ=θ+D,其中  θ=mid⁡Bbg−mid⁡B1g.\operatorname{mid}_{B_b}\bar a-\operatorname{mid}_{B_1}\bar a=\theta+D,\qquad\text{其中}\ \ \theta=\operatorname{mid}_{B_b}g-\operatorname{mid}_{B_1}g.

θ\theta 是 gg 自己的偏移:没有内部点的增量时,首尾两个块本来就差多少。

例

取 b=4b=4。为了简单,设一共只有两个点,都在 RR 里:一个红点在 B2B_2 里,一个蓝点在 B3B_3 里,都是内部点。细分之前,aa 在红点处从 0 升到 12\frac12,在蓝点处又降回 0。四个旧块上的振荡是 0,12,12,00,\frac12,\frac12,0,和是 1。注意 aa 在最低和最高的块上都恒等于 0:两个颜色抵消了,首尾没有偏移。

细分之后,设红点落进 CrC_r,蓝点落进 Cr′C_{r'}。没有别的点,所以 g=ag=a,θ=0\theta=0,而 D=5−2r8−5−2r′8=r′−r4D=\frac{5-2r}8-\frac{5-2r'}8=\frac{r'-r}4。aˉ\bar a 在 B1B_1 上是 0,在 B4B_4 上是常数 DD,偏移就是 DD。十六种等可能的结果里:

∣D∣|D|014\frac1412\frac1234\frac34
有几种结果4642

所以 E∣D∣=116(6⋅14+4⋅12+2⋅34)=516\EE|D|=\frac1{16}\bigl(6\cdot\frac14+4\cdot\frac12+2\cdot\frac34\bigr)=\frac5{16}。细分之前首尾恰好持平,细分之后只有四分之一的机会仍然持平。式 (8) 右边的第二项是 2⋅516=582\cdot\frac5{16}=\frac58;引理 2.3 保证的是 142≈0.35\frac14\sqrt2\approx0.35。

偏移是“一个固定的数 θ\theta,加上一个随机的和 DD”。对某一次骰子,DD 可能恰好把原来的偏移 θ\theta 抵消掉。但 DD 是若干个独立的对称增量之和,这样的和不会总挤在某个数附近,所以平均起来抵消不掉。下一步把这句话变成不等式。

第四步:对称随机变量的和离任何定点都不近

先算增量的两个矩。Xi2X_i^2 不依赖颜色,rr 在 1,…,b1,\dots,b 上等可能:

EXi2=1b∑r=1b(b+1−2r2b)2=b2−112b2=:σ2.\EE X_i^2=\frac1b\sum_{r=1}^b\Bigl(\frac{b+1-2r}{2b}\Bigr)^2=\frac{b^2-1}{12b^2}=:\sigma^2.

另外 ∣Xi∣≤b−12b|X_i|\le\frac{b-1}{2b}(r=1r=1 或 bb 时取到),所以 Xi4≤(b−12b)2Xi2X_i^4\le\bigl(\frac{b-1}{2b}\bigr)^2X_i^2,于是

EXi4≤(b−12b)2σ2≤b2−14b2 σ2=3σ4.\EE X_i^4\le\Bigl(\frac{b-1}{2b}\Bigr)^2\sigma^2\le\frac{b^2-1}{4b^2}\,\sigma^2=3\sigma^4.
推导这两个矩
  1. 二阶矩。b+1−2rb+1-2r 取 b−1, b−3, …, −(b−1)b-1,\ b-3,\ \dots,\ -(b-1),它们的平方和是 b(b2−1)3\frac{b(b^2-1)}3(比如 b=4b=4:9+1+1+9=20=4⋅1539+1+1+9=20=\frac{4\cdot15}3)。所以 EXi2=1b⋅b(b2−1)3⋅4b2=b2−112b2\EE X_i^2=\frac1b\cdot\frac{b(b^2-1)}{3\cdot4b^2}=\frac{b^2-1}{12b^2}。论文换了一种算法:Xi2=(r−b+12)2/b2X_i^2=\bigl(r-\frac{b+1}2\bigr)^2/b^2,而 rr 的方差是 (b+1)(2b+1)6−(b+1)24=b2−112\frac{(b+1)(2b+1)}6-\frac{(b+1)^2}4=\frac{b^2-1}{12}。
  2. 四阶矩。中间那个不等号是 (b−1)2≤b2−1(b-1)^2\le b^2-1,也就是 b≥1b\ge1。最后的等号是 b2−14b2=3σ2\frac{b^2-1}{4b^2}=3\sigma^2。

bb 很大时 σ2\sigma^2 接近 112\frac1{12},这是 [−12,12][-\frac12,\frac12] 上均匀分布的方差。

工具 · Cauchy–Schwarz 不等式(期望形式)

对任意随机变量 X,YX,Y,有 E∣XY∣≤E[X2]⋅E[Y2]\EE|XY|\le\sqrt{\EE[X^2]\cdot\EE[Y^2]}。

证明:对任意 s>0s>0,由 (s ∣X∣−∣Y∣/s)2≥0(\sqrt s\,|X|-|Y|/\sqrt s)^2\ge0 得 ∣XY∣≤12(sX2+Y2/s)|XY|\le\frac12\bigl(sX^2+Y^2/s\bigr)。取期望,再取 s=E[Y2]/E[X2]s=\sqrt{\EE[Y^2]/\EE[X^2]}。有一个二阶矩是 0 时结论显然。

Fact 2.2(对称随机变量的和)论文 Fact 2.2,证明在附录 A.1

设 X1,…,XmX_1,\dots,X_m 相互独立,每个都是对称的,并且 EXi2=σ2\EE X_i^2=\sigma^2,EXi4≤3σ4\EE X_i^4\le3\sigma^4,其中 σ>0\sigma>0。则对任意实数 θ\theta,

E ∣θ+X1+⋯+Xm∣ ≥ σm3.\EE\,\bigl|\theta+X_1+\dots+X_m\bigr|\ \ge\ \sigma\sqrt{\frac m3}.

对均匀的 ±1\pm1 符号,这是 Khintchine 不等式(常数差一些);四阶矩的条件对它们成立,对这里的增量上面刚验证过。

推导Fact 2.2

记 Y=X1+⋯+XmY=X_1+\dots+X_m。m=0m=0 时右边是 0,结论显然;下面设 m≥1m\ge1。

  1. 平移不帮忙。每个 XiX_i 对称,它们又独立,所以 YY 和 −Y-Y 同分布,E∣θ+Y∣=E∣θ−Y∣\EE|\theta+Y|=\EE|\theta-Y|。两者取平均,再用 ∣θ+y∣+∣θ−y∣≥2∣y∣|\theta+y|+|\theta-y|\ge2|y|,得到 E∣θ+Y∣≥E∣Y∣\EE|\theta+Y|\ge\EE|Y|。
  2. 二阶矩。EY2=∑i,jE[XiXj]=∑iEXi2=mσ2\EE Y^2=\sum_{i,j}\EE[X_iX_j]=\sum_i\EE X_i^2=m\sigma^2,因为 i≠ji\ne j 时 E[XiXj]=EXi⋅EXj=0\EE[X_iX_j]=\EE X_i\cdot\EE X_j=0。
  3. 四阶矩。展开 Y4Y^4。一项里如果有某个下标出现奇数次,由独立性和对称性,它的期望是 0。剩下两类:四个下标全相同的 Xi4X_i^4;两两配对的 Xi2Xj2X_i^2X_j^2(i<ji<j),每一对出现 (42)=6\binom42=6 次。所以 EY4=∑iEXi4+6∑i<jEXi2 EXj2≤3mσ4+3m(m−1)σ4=3m2σ4.\EE Y^4=\sum_i\EE X_i^4+6\sum_{i<j}\EE X_i^2\,\EE X_j^2\le3m\sigma^4+3m(m-1)\sigma^4=3m^2\sigma^4.
  4. 两次 Cauchy–Schwarz。把 Y2Y^2 拆成 ∣Y∣1/2⋅∣Y∣3/2|Y|^{1/2}\cdot|Y|^{3/2},再把 ∣Y∣3|Y|^3 拆成 ∣Y∣⋅Y2|Y|\cdot Y^2: EY2≤(E∣Y∣)1/2(E∣Y∣3)1/2,E∣Y∣3≤(EY2)1/2(EY4)1/2.\EE Y^2\le\bigl(\EE|Y|\bigr)^{1/2}\bigl(\EE|Y|^3\bigr)^{1/2},\qquad\EE|Y|^3\le\bigl(\EE Y^2\bigr)^{1/2}\bigl(\EE Y^4\bigr)^{1/2}. 第一式两边平方,再代入第二式:(EY2)2≤E∣Y∣⋅(EY2)1/2(EY4)1/2(\EE Y^2)^2\le\EE|Y|\cdot(\EE Y^2)^{1/2}(\EE Y^4)^{1/2},也就是 E∣Y∣≥(EY2)3/2/(EY4)1/2\EE|Y|\ge(\EE Y^2)^{3/2}/(\EE Y^4)^{1/2}。论文在这里用一次 Hölder 不等式(指数 32\frac32 和 3),结果相同。
  5. 代入。E∣Y∣≥m3/2σ33 mσ2=σm3\EE|Y|\ge\dfrac{m^{3/2}\sigma^3}{\sqrt3\,m\sigma^2}=\sigma\sqrt{\dfrac m3}。
交互 · 12Fact 2.2:增量之和离哪个点都不近

取 mm 个独立的增量,每个在 b+1−2r2b\frac{b+1-2r}{2b}(r=1,…,br=1,\dots,b)这 bb 个值上等可能。上图是它们的和 DD 的精确分布,红色竖线是 −θ-\theta 的位置,E∣θ+D∣\EE|\theta+D| 就是 DD 到这条线的平均距离。下图把这个平均距离画成 θ\theta 的函数,三条横线是证明里依次出现的三个下界。

拖动平移量:曲线在 θ=0\theta=0 处最低,这是推导的第 1 步。三条横线自上而下排列,这是第 4、5 步。把项数从 1 调到 12,三条线一起按 m\sqrt m 的速度上升:项越多,和摊得越开,离任何一个点的平均距离就越大。mm 个内部点给出 m\sqrt m 量级的增益,来源就在这里。

第五步:合起来,得到引理 2.3

把骰子分成两批:内部区域里的点的骰子,和 II 里其余的点的骰子(高度在 B1B_1 里、在 BbB_b 里、或者在 JJ 之外的点)。先把第二批看成已知的,只对第一批求平均。这时 gg 是一个固定的函数,θ\theta 是一个固定的数。而 DD 是 mRm_R 个独立的对称增量之和,它们和第二批骰子无关。

对 θ+D\theta+D 用 Fact 2.2,取 m=mRm=m_R。偏移的绝对值的条件期望至少是

σmR3=b2−16b mR ≥ 18mR.\sigma\sqrt{\frac{m_R}3}=\frac{\sqrt{b^2-1}}{6b}\,\sqrt{m_R}\ \ge\ \frac18\sqrt{m_R}.

后一个不等号是 16(b2−1)≥9b216(b^2-1)\ge9b^2,b≥2b\ge2 时成立。这个下界对任何 θ\theta 都一样,所以第二批骰子怎么掷都压不低它;再对第二批求平均(第 10 节的“逐层平均”),它仍然成立:

E ∣mid⁡Bbaˉ−mid⁡B1aˉ∣ ≥ 18mR.\EE\,\bigl|\operatorname{mid}_{B_b}\bar a-\operatorname{mid}_{B_1}\bar a\bigr|\ \ge\ \frac18\sqrt{m_R}.

代进式 (8):

引理 2.3(局部增益)论文 Lemma 2.3

在局部实验里,

E[ ∑k=1bosc⁡J ak] ≥ ∑ℓ=1bosc⁡Bℓ a + 14mR.\EE\Bigl[\,\sum_{k=1}^b\operatorname{osc}_J\,a_k\Bigr]\ \ge\ \sum_{\ell=1}^b\operatorname{osc}_{B_\ell}\,a\ +\ \frac14\sqrt{m_R}.

左边是细分之后 RR 里 bb 个竖条上的振荡之和,右边第一项是细分之前 bb 个横条上的振荡之和。把 RR 从横着切改成竖着切,振荡之和平均至少多出 14mR\frac14\sqrt{m_R}。论文说这是整个证明的核心。

这一笔有多大?一个比较矩形里平均有 nM\frac nM 个点,其中 1−2b1-\frac2b 的比例在内部区域里。所以典型情况下,每个比较矩形每细分一次,振荡之和多出 n/M\sqrt{n/M} 的量级。这就是第 5 节末尾说的那个代价:每一层多付一个平方根。

还有一点后面要用:增益的大小 mRm_R 只数点,不看颜色。染色会改变 θ\theta,也会改变各个 XiX_i 的符号;但 XiX_i 是对称的,变号之后分布不变,而 Fact 2.2 对任何 θ\theta 都成立。

交互 · 13一个比较矩形里的引理 2.3

取 b=8b=8,只画一个比较矩形 RR,并且只考虑 RR 里的点。左图里横带是 8 个旧块,竖条是 8 个子区间,灰色的两条是最低和最高的旧块。每按一次“重新掷骰子”,各点重新随机落进一个子区间。表格对比每个旧块上 aa 的振荡和这一次 aˉ\bar a 的振荡。读数给出式 (5) 右边的两部分和左边的新的和,再给出多次掷骰子的平均。

多掷几次,再按“掷 2 万次取平均”。三件事应该能看出来。aˉ\bar a 的振荡之和,长期平均不小于旧的和,这是基准。偏移的两倍的长期平均明显大于 14mR\frac14\sqrt{m_R},后者只是一个保守的下界。新的和的长期平均又比“基准加增益”大不少:振荡的三角不等式和引理 2.1 都丢掉了一些东西。把点数调大,偏移大致按点数的平方根增长。

一个矩形能给出什么,还缺什么

引理 2.3 直接用在整个正方形上,会得到什么?取 I=[0,1]I=[0,1],J=[0,1)J=[0,1),RR 就是单位正方形(去掉上边)。每个子函数的振荡是 aka_k 的两个值之差,每个值又是 FF 的两个值的平均,所以不超过 2∥F∥∞2\norm F_\infty;引理左边至多是 2b E∥F∥∞2b\,\EE\norm F_\infty。固定各点的高度,横坐标独立均匀,第一位数字做的就是局部实验,于是

E ∥F∥∞ ≥ mR8b.\EE\,\norm F_\infty\ \ge\ \frac{\sqrt{m_R}}{8b}.

典型情况下内部区域里有大约 (1−2b)n(1-\frac2b)n 个点,所以这是 n\sqrt n 的量级,比要证的 log⁡3/2n\log^{3/2}n 大得多。

易错

这不是差异的下界。上面的染色是在抽横坐标之前定好的;看过点再染色,可以做到 O(log⁡3/2n)O(\log^{3/2}n)(Nikolov)。我们要的是对全部 2n2^n 种染色同时成立的结论,所以每一种固定的染色失败的概率要远小于 2−n2^{-n}。一个矩形上比较一次,只给出一个期望,离这个要求很远。

所以后面四节不用一个大矩形,而是用 MM 个小矩形,并且在 hh 个逐渐变细的横向尺度上重复这个比较。引理 2.3 是为此设计的:它的输出(竖条上的振荡之和)和输入(横条上的振荡之和)是同一种量,上一层的输出可以直接当下一层的输入。这样做之后还缺两样东西:

点数增益取决于每个内部区域分到多少点。点会不会太挤?第 17 节,引理 2.5。
期望引理 2.3 只说期望。实际值会不会比期望低很多?第 15、16 节,引理 2.6。

§ 15 · 进阶改一位数字,影响有多大

引理 2.3 说的是平均。要把平均变成“几乎总是”,需要知道细分之后的振荡对每一个骰子有多敏感。答案是:很不敏感。这一节不含概率,只用第 12 节的两条性质:振荡不怕加常数,以及反向三角不等式 ∣osc⁡B(f+g)−osc⁡B f∣≤osc⁡B g\bigl|\operatorname{osc}_B(f+g)-\operatorname{osc}_B\,f\bigr|\le\operatorname{osc}_B\,g。论文把这件事放在引理 2.6 的证明里,是它的 Step 1。

一位数字的影响论文引理 2.6 的证明,Step 1

在第 jj 层,只改变一个点 PiP_i 的下一位数字,也就是把它在所在的旧区间 II 里挪到另一个子区间,其余的点不动。那么在第 j+1j+1 阶段的 bMbM 个小矩形里,只有 bb 个的振荡可能变,每个至多变 1。所以

Zj+1 至多变 bbM=1M.Z_{j+1}\ \text{至多变}\ \frac b{bM}=\frac1M.
推导为什么只有 bb 个小矩形受影响
  1. 哪些函数会变。设 PiP_i 从 CrC_r 挪到 Cr′C_{r'}。它的系数只在 II 的 bb 个子区间里变:在别的区间看来,它仍然在同一侧。它在 aka_k 里的系数是 11、12\frac12 或 00,挪了之后还是这三个数之一,所以 ak(y) 的变化=ek 1{Vi≤y},∣ek∣≤1,a_k(y)\ \text{的变化}=e_k\,\one_{\{V_i\le y\}},\qquad|e_k|\le1, 其中 eke_k 是 χi\chi_i 乘以新旧系数之差。这是一个只在高度 ViV_i 处跳一次的阶梯函数。
  2. 高度不在块里。看第 j+1j+1 阶段的一个小矩形 Ck×BC_k\times B,BB 是一个新块。如果 Vi∉BV_i\notin B,那么 BB 要么整个在 ViV_i 下面,要么整个在它上面,上面那个变化在 BB 上是常数。振荡不怕加常数,osc⁡B ak\operatorname{osc}_B\,a_k 不变。
  3. 高度在块里。变化在 BB 上只取 0 和 eke_k 两个值,它自己的振荡是 ∣ek∣≤1|e_k|\le1。由反向三角不等式,osc⁡B ak\operatorname{osc}_B\,a_k 的变化不超过 1。
  4. 合计。含 ViV_i 的新块只有一个,记作 JJ。受影响的小矩形就是 C1×J,…,Cb×JC_1\times J,\dots,C_b\times J 这 bb 个,它们合起来是比较矩形 I×JI\times J。Zj+1Z_{j+1} 是 bMbM 个振荡之和除以 bMbM,所以至多变 bbM\frac b{bM}。

这里看得出为什么要量振荡。固定已揭示信息,改变一个点的下一位横坐标数字,会改变它上方所有高度上的 aka_k 都变了,受影响的块可能有很多。但其中只有一个块的振荡会变,其余的块只是整体平移了一下。如果量的是 ∣ak∣|a_k| 的最大值,这些块就都得算进去。

交互 · 14挪动一个点,哪些格子的振荡会变

取 b=4b=4,画出一个区间的 4 个子区间和 4 个新块。右边的表格是这 16 个小矩形各自的振荡。点一下选中一个点,再把它挪到别的子区间:振荡发生变化的格子会被标出来。

不管怎么挪,标出来的格子都在同一行,也就是那个点所在的块;变化量的绝对值之和不超过 b=4b=4。到这里,关于“细分一步”的两件事都齐了:

平均引理 2.3:每个比较矩形里的振荡之和,平均至少多出 14mR\frac14\sqrt{m_R}。
影响改一个点的下一位数字,只有它所在的那个比较矩形里的振荡之和会变,至多变 bb;Zj+1Z_{j+1} 至多变 1M\frac1M。

§ 16 · 核心推导把 hh 层接起来:计划、条件增益和波动

前三节只看一个比较矩形、一次细分。现在回到第 12 节的 Z0,Z1,…,ZhZ_0,Z_1,\dots,Z_h,把 hh 层接起来。这一节先列出剩下的计划,然后完成其中两步:引理 2.4 说每一层平均涨多少,引理 2.6 说实际的涨幅很难比平均值低太多。

第一步:把涨幅拆成“条件增益”和“波动”

ZjZ_j 由 Fj\mathcal F_j 决定。所以每一层的实际涨幅可以写成两段:

Zj+1−Zj=E[Zj+1−Zj∣Fj]⏟条件增益:这一层平均涨多少+(Zj+1−E[Zj+1∣Fj])⏟波动:实际值减去平均值.Z_{j+1}-Z_j=\underbrace{\EE\bigl[Z_{j+1}-Z_j\bigm|\mathcal F_j\bigr]}_{\text{条件增益:这一层平均涨多少}}+\underbrace{\bigl(Z_{j+1}-\EE[Z_{j+1}\mid\mathcal F_j]\bigr)}_{\text{波动:实际值减去平均值}}.

对 j=0,…,h−1j=0,\dots,h-1 求和:

初值、条件增益和波动项论文式 (10)
Zh=Z0+∑j<hE[Zj+1−Zj∣Fj]+T,T=∑j<h(Zj+1−E[Zj+1∣Fj]).(10)Z_h=Z_0+\sum_{j<h}\EE\bigl[Z_{j+1}-Z_j\bigm|\mathcal F_j\bigr]+T,\qquad T=\sum_{j<h}\bigl(Z_{j+1}-\EE[Z_{j+1}\mid\mathcal F_j]\bigr).\tag{10}

E[Zj+1−Zj∣Fj]\EE[Z_{j+1}-Z_j\mid\mathcal F_j] 叫做第 jj 层的条件增益(conditional gain),TT 叫做波动项(fluctuation term)。由式 (9),Z0≥0Z_0\ge0,Zh≤2∥F∥∞Z_h\le2\norm F_\infty。所以只要右边后两项合起来很大,∥F∥∞\norm F_\infty 就很大。剩下的证明分三步:

(a)每一层的条件增益至少是一个数 δj\delta_j,它只依赖点的位置,不依赖染色。引理 2.4,本节。
(b)M≤nM\le n 且 bb 很大时,除了一个指数小的概率,每个 δj\delta_j 都很大。引理 2.5,第 17 节。
(c)对一种固定的染色,波动项 TT 很难是一个很大的负数。引理 2.6,本节。

命题 2.7(第 18 节)把三步合起来,再对所有染色用并集界。

第二步:每一层平均涨多少

第 jj 层有 MM 个比较矩形,面积都是 1M\frac1M。和第 13 节一样,比较矩形 RR 的内部区域是 RR 去掉最低和最高的旧小矩形,mRm_R 是里面的点数。

引理 2.4(条件增益)论文 Lemma 2.4

对 0≤j<h0\le j<h,以概率 1,

E[ Zj+1−Zj∣Fj] ≥ δj:=14bM∑RmR,(11)\EE\bigl[\,Z_{j+1}-Z_j\bigm|\mathcal F_j\bigr]\ \ge\ \delta_j:=\frac1{4bM}\sum_R\sqrt{m_R},\tag{11}

求和遍历第 jj 层的 MM 个比较矩形。

推导引理 2.4
  1. 按比较矩形分组。每个旧的小矩形是恰好一个比较矩形的横条,每个新的小矩形是恰好一个比较矩形的竖条:R=I×JR=I\times J 里有 bb 个旧的 I×BℓI\times B_\ell 和 bb 个新的 Ck×JC_k\times J。所以 bM (Zj+1−Zj)=∑R(∑kosc⁡J aCk−∑ℓosc⁡Bℓ aI).bM\,(Z_{j+1}-Z_j)=\sum_R\Bigl(\sum_k\operatorname{osc}_J\,a_{C_k}-\sum_\ell\operatorname{osc}_{B_\ell}\,a_I\Bigr).
  2. 每个比较矩形里是一次局部实验。在 Fj\mathcal F_j 已知的条件下:各点的高度已知;每个点相对每个旧区间 II 的位置已知;II 里面的点,下一位数字是独立、均匀的骰子,决定它落进哪个子区间。这正是第 14 节的局部实验,而且在所有的 RR 里同时进行。aIa_I 是旧函数,aC1,…,aCba_{C_1},\dots,a_{C_b} 是子函数。
  3. 求和。对每个 RR 用引理 2.3,括号里的量的条件期望至少是 14mR\frac14\sqrt{m_R}。全部加起来,除以 bMbM,得到 (11)。

δj\delta_j 不是负的,所以每一层平均都不降。它有多大?如果每个内部区域里都有 nM\frac nM 个左右的点,MM 个 n/M\sqrt{n/M} 加起来再除以 4bM4bM,δj\delta_j 大约是 14bn/M\frac1{4b}\sqrt{n/M}。hh 层累积起来是 hn/Mh\sqrt{n/M} 的量级。但“如果”不是保证,第 17 节处理它。

δj\delta_j 还有一个特点:mRm_R 只数点,不看颜色。所以 δj\delta_j 对 2n2^n 种染色都是同一个数。

式 (10) 中间那个和里的每一项至少是 δj≥0\delta_j\ge0。TT 里的每一项,在 Fj\mathcal F_j 已知时均值是 0,可正可负。要证 ZhZ_h 很大,只怕一件事:TT 是一个很大的负数,把增益吃掉了。下面证明这件事的概率极小。

第三步:有界差分

工具 · 有界的零均值随机变量

设 c>0c>0、EY=0\mathbb E Y=0 且 ∣Y∣≤c|Y|\le c。则对任意实数 tt,EetY≤cosh⁡(tc)≤et2c2/2\mathbb E e^{tY}\le\cosh(tc)\le e^{t^2c^2/2}。若 c=0c=0,则 Y=0Y=0 几乎处处成立,结论直接成立。

证明:函数 y↦etyy\mapsto e^{ty} 是凸的,在区间 [−c,c][-c,c] 上它位于连接两个端点的弦的下方:

etY≤c+Y2c etc+c−Y2c e−tc.e^{tY}\le\frac{c+Y}{2c}\,e^{tc}+\frac{c-Y}{2c}\,e^{-tc}.

取期望,E Y=0\EE\,Y=0,右边变成 12(etc+e−tc)=cosh⁡(tc)\frac12(e^{tc}+e^{-tc})=\cosh(tc)。第二个不等号 cosh⁡x≤ex2/2\cosh x\le e^{x^2/2} 在第 5 节证过。

条件期望的版本同样成立:如果在某些已知信息下 YY 的条件期望是 0,并且 ∣Y∣≤c|Y|\le c,那么 etYe^{tY} 的条件期望不超过 et2c2/2e^{t^2c^2/2}。证明一字不改,只是把 E\EE 都换成条件期望。

交互 · 15弦在上,曲线在下

横轴是 YY 的取值范围 [−c,c][-c,c]。曲线是 etye^{ty},直线段是连接两个端点的弦。圆圈是一个均值为 0 的分布,圆越大概率越大。这些点都在曲线上,它们的重心(菱形)横坐标是 0,纵坐标就是 E etY\EE\,e^{tY}。

圆圈都在弦的下方,所以重心也在弦的下方,而弦在 0 处的高度正是 cosh⁡(tc)\cosh(tc)。只取两个端点时,重心恰好落在弦上,这是最坏的情形。再往上一点是 et2c2/2e^{t^2c^2/2}。把 tctc 调成负数,图左右翻转,结论不变。

Fact A.1(有界差分)论文附录 A.2

设 ξ1,…,ξn\xi_1,\dots,\xi_n 相互独立,随机元素 Γ\Gamma 与向量 (ξ1,…,ξn)(\xi_1,\dots,\xi_n) 独立。设 φ\varphi 可测,f=φ(Γ,ξ1,…,ξn)f=\varphi(\Gamma,\xi_1,\dots,\xi_n) 可积,且 c>0c>0。固定 Γ\Gamma 的任意取值,只改变一个 ξi\xi_i 时,φ\varphi 的值至多变化 cc。则对任意实数 tt,以概率 1,

E[exp⁡(t (f−E[f∣Γ]))∣Γ] ≤ exp⁡(n t2c22).\EE\Bigl[\exp\bigl(t\,(f-\EE[f\mid\Gamma])\bigr)\Bigm|\Gamma\Bigr]\ \le\ \exp\Bigl(\frac{n\,t^2c^2}2\Bigr).
推导Fact A.1

Γ\Gamma 和各个 ξi\xi_i 独立,所以“给定 Γ\Gamma”就是把 φ\varphi 的第一个变量固定成一个值 γ\gamma。于是只要对每个固定的 γ\gamma,证明不带条件的 Eexp⁡(t(f−Ef))≤exp⁡(nt2c22)\EE\exp\bigl(t(f-\EE f)\bigr)\le\exp\bigl(\frac{nt^2c^2}2\bigr),其中 f=φ(γ,ξ1,…,ξn)f=\varphi(\gamma,\xi_1,\dots,\xi_n)。下面固定 γ\gamma。

  1. 一个一个地揭示。令 fi=E[f∣ξ1,…,ξi]f_i=\EE[f\mid\xi_1,\dots,\xi_i],即知道前 ii 个变量之后对 ff 的最佳猜测;f0=Eff_0=\EE f,fn=ff_n=f。记 Yi=fi−fi−1Y_i=f_i-f_{i-1},则 f−Ef=∑i=1nYif-\EE f=\sum_{i=1}^nY_i。
  2. 每一步的均值是 0。由“逐层平均”,E[fi∣ξ1,…,ξi−1]=fi−1\EE[f_i\mid\xi_1,\dots,\xi_{i-1}]=f_{i-1},所以 E[Yi∣ξ1,…,ξi−1]=0\EE[Y_i\mid\xi_1,\dots,\xi_{i-1}]=0。
  3. 每一步有界。因为变量相互独立,fif_i 是这样算出来的:把前 ii 个变量固定在揭示出的值上,对后面的变量求平均。把它看成第 ii 个变量的函数:第 ii 个变量换一个值,被平均的每一项至多变 cc,所以平均值也至多变 cc。而 fi−1f_{i-1} 是这个函数再对第 ii 个变量求平均。一个取值范围不超过 cc 的函数,离它自己的平均值不会超过 cc。所以 ∣Yi∣≤c|Y_i|\le c。
  4. 从最后一个往前剥。Y1,…,Yn−1Y_1,\dots,Y_{n-1} 由 ξ1,…,ξn−1\xi_1,\dots,\xi_{n-1} 决定,可以提出来;对 YnY_n 用上面工具的条件期望版本: E et∑i≤nYi=E[et∑i<nYi⋅E[etYn∣ξ1,…,ξn−1]]≤et2c2/2⋅E et∑i<nYi.\begin{aligned}\EE\,e^{t\sum_{i\le n}Y_i}&=\EE\Bigl[e^{t\sum_{i<n}Y_i}\cdot\EE\bigl[e^{tY_n}\bigm|\xi_1,\dots,\xi_{n-1}\bigr]\Bigr]\\[4pt]&\le e^{t^2c^2/2}\cdot\EE\,e^{t\sum_{i<n}Y_i}.\end{aligned} 重复 nn 次,得到 exp⁡(nt2c22)\exp\bigl(\frac{nt^2c^2}2\bigr)。

这是有界差分不等式(McDiarmid)的矩母函数形式,证明用的是 Azuma 的鞅论证;标准版本的常数更好,这里用不着。像 Y1,Y2,…Y_1,Y_2,\dots 这样“在此前揭示的一切都已知时均值为 0”的一串量,叫做鞅差。

交互 · 16逐个揭示:最佳猜测每一步至多动 cc

一个具体的例子:掷 10 个 8 面的骰子,ff 是“一次都没出现的点数有几种”。改掉任何一个骰子,ff 至多变 1,所以 c=1c=1。折线是 fif_i,即知道前 ii 个骰子之后对 ff 的最佳猜测:每揭示一个骰子,猜测更新一次。浅色的竖条是下一步允许落入的范围 fi−1±cf_{i-1}\pm c。

这个例子里 fif_i 有现成的公式。前 ii 个骰子出现过 ss 种点数时,剩下 8−s8-s 种里的每一种,在后面 10−i10-i 个骰子里都不出现的概率是 (78)10−i\bigl(\frac78\bigr)^{10-i},所以 fi=(8−s)(78)10−if_i=(8-s)\bigl(\frac78\bigr)^{10-i}。

图和推导对得上两处。每一步有正有负:出了新点数就往下跳,点数重复就往上挪一点,两者按概率加权平均是 0,这是第 2 步。(第一个骰子例外:它一定是新点数,这一步恰好是 0。)每一步都落在浅色竖条里,这是第 3 步。再做 30 轮,终点散在 Ef\EE f 两侧一个不大的范围里。

第四步:波动项很难是大的负数

引理 2.6(波动)论文 Lemma 2.6

固定一种染色,TT 是式 (10) 里的波动项。则对任意 λ>0\lambda>0,

Pr⁡[ T≤−λ ] ≤ exp⁡(−λ2M22hn).(13)\Pr[\,T\le-\lambda\,]\ \le\ \exp\Bigl(-\frac{\lambda^2M^2}{2hn}\Bigr).\tag{13}
推导引理 2.6

把第 jj 层的波动记作 Tj=Zj+1−E[Zj+1∣Fj]T_j=Z_{j+1}-\EE[Z_{j+1}\mid\mathcal F_j],于是 T=T0+⋯+Th−1T=T_0+\dots+T_{h-1}。

  1. 一位数字。改变一个点的下一位数字,Zj+1Z_{j+1} 至多变 1M\frac1M。这是第 15 节。
  2. 一层。Zj+1Z_{j+1} 是 nn 个下一位数字的函数,另外还依赖 Fj\mathcal F_j 里的信息(各点的高度和前 jj 位数字)。这 nn 个数字相互独立,也和 Fj\mathcal F_j 独立。用 Fact A.1:Γ\Gamma 取 Fj\mathcal F_j 里的信息,ξ1,…,ξn\xi_1,\dots,\xi_n 取下一位数字,c=1Mc=\frac1M,再把 tt 换成 −t-t。对任意 t>0t>0, E[e−tTj∣Fj] ≤ exp⁡(n t22M2).\EE\bigl[e^{-tT_j}\bigm|\mathcal F_j\bigr]\ \le\ \exp\Bigl(\frac{n\,t^2}{2M^2}\Bigr).
  3. hh 层连乘。T0,…,Th−2T_0,\dots,T_{h-2} 在 Fh−1\mathcal F_{h-1} 已知时都已经定了,可以提出来: E e−tT=E[e−t(T0+⋯+Th−2)⋅E[e−tTh−1∣Fh−1]]≤exp⁡(n t22M2)⋅E e−t(T0+⋯+Th−2).\EE\,e^{-tT}=\EE\Bigl[e^{-t(T_0+\dots+T_{h-2})}\cdot\EE\bigl[e^{-tT_{h-1}}\bigm|\mathcal F_{h-1}\bigr]\Bigr]\le\exp\Bigl(\frac{n\,t^2}{2M^2}\Bigr)\cdot\EE\,e^{-t(T_0+\dots+T_{h-2})}. 再对 Fh−2\mathcal F_{h-2}、Fh−3\mathcal F_{h-3} 一直到 F0\mathcal F_0 这样做,得到 E e−tT≤exp⁡(hn t22M2)\EE\,e^{-tT}\le\exp\bigl(\frac{hn\,t^2}{2M^2}\bigr)。
  4. Markov 不等式,再选 tt。和第 5 节一样,“T≤−λT\le-\lambda”就是“e−tT≥etλe^{-tT}\ge e^{t\lambda}”,所以 Pr⁡[T≤−λ]≤e−tλ⋅E e−tT≤exp⁡(−tλ+hn t22M2).\Pr[T\le-\lambda]\le e^{-t\lambda}\cdot\EE\,e^{-tT}\le\exp\Bigl(-t\lambda+\frac{hn\,t^2}{2M^2}\Bigr). 指数是 tt 的二次函数,在 t=λM2hnt=\frac{\lambda M^2}{hn} 处最小,最小值是 −λ2M22hn-\frac{\lambda^2M^2}{2hn}。

这一步不需要各层之间独立。给定每一层的已知信息,下一位数字是相互独立的骰子,改变其中一位,Zj+1Z_{j+1} 至多变 1M\frac1M。先对这一层用条件期望形式的指数矩界,再逐层取条件期望,hh 层各贡献一个因子 exp⁡(n t22M2)\exp\bigl(\frac{n\,t^2}{2M^2}\bigr),式 (13) 的指数里的 hnhn 和 M2M^2 就是这么累积出来的。

易错

不能把这个论证简化成“TT 是 hnhn 位数字的函数,每一位的影响至多 1M\frac1M”。改变较早的一位数字,会连带改变后面各层的 ZZ 和条件期望,对 TT 的影响可以超过 1M\frac1M。一个小例子:b=3b=3、h=2h=2、n=3n=3,所以 M=3M=3;颜色 (+1,−1,+1)(+1,-1,+1),纵坐标 (110,34,1920)(\frac1{10},\frac34,\frac{19}{20}),横坐标 (1118,118,12)(\frac{11}{18},\frac1{18},\frac12)。把第二个点的横坐标改成 1318\frac{13}{18},也就是把它的第一位数字从 1 改成 3,其余不动。直接按定义算,TT 从 1481\frac{14}{81} 变成 −2281-\frac{22}{81},变了 49>13\frac49>\frac13。

所以界是逐层用的:每一层只让当前这一位数字变动,此前的数字都算作已知。

式 (13) 是对一种固定的染色说的:ZjZ_j 依赖染色,换一种染色,TT 就是另一个随机变量。要对付全部 2n2^n 种染色,右边得乘上 2n2^n,所以 λ\lambda 要取得让指数 λ2M22hn\frac{\lambda^2M^2}{2hn} 比 nln⁡2n\ln2 大。第 18 节做这件事。

Zh − Z0 + 各层的条件增益之和 第 j 层 ≥ δj ≥ 0,式 (11) δj 只和点的位置有关,与染色无关 第 17 节:各层都是好层时,总和 ≥ hδ∗ 波动项 T 每层 n 位新数字,各使 Zj+1 至多变 1/M 式 (13):T ≤ −hδ∗/4 的概率极小 对每种染色分别用,再乘 2n 两件事都成立:Zh > hδ∗ − hδ∗/4,于是 max|F| ≥ Zh/2 > hδ∗/8
后面两节的走向。δ∗\delta_* 是第 17 节定义的一个数,量级是 n/M\sqrt{n/M}。用于下界各层条件增益的好事件只依赖点的位置,对所有染色统一处理;条件增益本身仍可依赖染色。右边的波动项逐个染色处理,再做并集界。
交互 · 17振荡逐层累积

取 b=8b=8、h=3h=3,所以 M=64M=64。随机撒 nn 个点,选一种染色方式,算出 Z0,…,Z3Z_0,\dots,Z_3。折线从 Z0Z_0 出发,把各层的 δj\delta_j 逐层加上去;折线每一段升高的量,是这一层条件增益的下界。表里后两列是第 17 节的内容:这一层有多大比例的点落在不重的内部区域里(至少一半就是好层),以及好层保证的下界 δ∗\delta_*。“多撒几十组”会重复实验,看 ZhZ_h 最小能到多少。

图里有三件事可以留意,说的都是眼前这组样本。第一,ZjZ_j 逐层上升,每一层的涨幅都远高于 δj\delta_j,表里可以直接对比。要记住 δj\delta_j 是条件期望的下界:单次实现的涨幅没有这样的保证,甚至可以是负的。第二,换成沿 Z 形曲线交替的染色,max⁡∣F∣\max|F| 小了很多,ZhZ_h 也跟着变小,但在这组样本上各层仍然在涨。第三,重复几十组,这些样本里 ZhZ_h 的最小值和平均值在同一个量级,没有哪一组掉到接近 0。这只是有限次实验的观察,说明不了尾概率;e−cne^{-cn} 这么小的概率要靠式 (13)。

§ 17 · 进阶点不会太挤:重的区域和好层

式 (11) 里的 δj=14bM∑RmR\delta_j=\frac1{4bM}\sum_R\sqrt{m_R},在每个内部区域里都有 nM\frac nM 个左右的点时,大约是 14bn/M\frac1{4b}\sqrt{n/M}。但这不是保证。如果点扎堆挤在少数几个内部区域里,或者大量落在各个比较矩形最低、最高的旧小矩形里,这个和就会很小;极端情况下所有点都不在内部区域里,δj=0\delta_j=0。这一节证明:M≤nM\le n 时,除了一个指数小的概率,每一层都能留住这个数的一个只依赖 bb 的比例。这是第 16 节计划里的 (b)。

重的区域,好层

一个比较矩形里平均有 nM\frac nM 个点。固定第 jj 层,对它的 MM 个比较矩形定义:

重内部区域里的点数 mRm_R 超过 bnM\dfrac{bn}M,即超过一个比较矩形的平均点数的 bb 倍,就说这个内部区域是重的(heavy)。
好层如果至少有 n2\dfrac n2 个点落在不重的内部区域里,就说这一层是好层(论文:the transition is good)。

好层说的是:至少一半的点摊得足够开。记 GjG_j 为“第 jj 层是好层”这个事件,G=G0∩⋯∩Gh−1G=G_0\cap\dots\cap G_{h-1} 为“hh 层都是好层”。它们只看点的位置,与染色无关。再记

好层保证的下界论文式 (12)
δ∗=18 b3/2nM.(12)\delta_*=\frac1{8\,b^{3/2}}\sqrt{\frac nM}.\tag{12}

它是“理想值” 14bn/M\frac1{4b}\sqrt{n/M} 的 12b\frac1{2\sqrt b} 倍。

引理 2.5(好层)论文 Lemma 2.5

设 0≤j<h0\le j<h。

(i) 在 GjG_j 上,δj≥δ∗\delta_j\ge\delta_*。

(ii) 如果 M≤nM\le n,那么 Pr⁡[ Gj 不发生 ]≤(48b)n/2\Pr[\,G_j\ \text{不发生}\,]\le\bigl(\frac{48}b\bigr)^{n/2},于是 Pr⁡[ G 不发生 ]≤h(48b)n/2\Pr[\,G\ \text{不发生}\,]\le h\bigl(\frac{48}b\bigr)^{n/2}。

推导(i):好层上 δj≥δ∗\delta_j\ge\delta_*
  1. 不重的区域里,平方根不比点数小太多。当 0<mR≤bn/M0<m_R\le bn/M 时,mR=mR/mR≥mR/bn/M\sqrt{m_R}=m_R/\sqrt{m_R}\ge m_R/\sqrt{bn/M}。当 mR=0m_R=0 时,直接验证最后的不等式即可:两边都是 0。
  2. 求和。只留下不重的区域。好层上它们里面的点数之和至少是 n2\frac n2: ∑RmR ≥ 1bn/M∑R 不重mR ≥ n2Mbn.\sum_R\sqrt{m_R}\ \ge\ \frac1{\sqrt{bn/M}}\sum_{R\ \text{不重}}m_R\ \ge\ \frac n2\sqrt{\frac M{bn}}.
  3. 除以 4bM4bM。δj≥14bM⋅n2Mbn=18 b3/2nM=δ∗\delta_j\ge\dfrac1{4bM}\cdot\dfrac n2\sqrt{\dfrac M{bn}}=\dfrac1{8\,b^{3/2}}\sqrt{\dfrac nM}=\delta_*。

为什么要把重的区域扔掉?平方根是凹的:同样多的点,摊在很多个区域里,平方根之和大;挤在一个区域里,平方根之和小。一个区域里的点数超过平均水平的 bb 倍之后,多出来的点对 ∑mR\sum\sqrt{m_R} 几乎没有贡献,干脆不算。

一层不是好层的概率

把每个比较矩形最低和最高的旧小矩形叫做边界带(boundary strips)。每个点要么在某个内部区域里,要么在边界带里。

推导(ii):M≤nM\le n 时,Pr⁡[ Gj 不发生 ]≤(48b)n/2\Pr[\,G_j\ \text{不发生}\,]\le\bigl(\frac{48}b\bigr)^{n/2}
  1. 不是好层意味着什么。这时落在不重的内部区域里的点不到 n2\frac n2 个。其余的点要么在重的内部区域里,要么在边界带里。所以有超过 n2\frac n2 个点落在集合 WW 里,WW 是边界带和全部重区域的并。
  2. 重的区域不多。没有哪个点同时在两个内部区域里,而每个重区域里有超过 bnM\frac{bn}M 个点,一共只有 nn 个点。所以重区域少于 n/bnM=Mbn\big/\frac{bn}M=\frac Mb 个。
  3. WW 不是固定的集合。哪些区域是重的,取决于点落在哪里。所以不能直接算“一个点落进 WW 的概率”。办法是把所有可能的情况都列出来:考虑每一组 S\mathcal S,它由少于 Mb\frac Mb 个内部区域组成,不管里面的区域重不重。记 WSW_{\mathcal S} 为边界带和 S\mathcal S 里的区域的并,ESE_{\mathcal S} 为“至少 n2\frac n2 个点落在 WSW_{\mathcal S} 里”这个事件。第 1、2 步说:如果第 jj 层不是好层,那么取 S\mathcal S 为重区域的全体,ESE_{\mathcal S} 就发生了。所以“不是好层”包含在这些 ESE_{\mathcal S} 的并里。
  4. WSW_{\mathcal S} 很小。固定 S\mathcal S。边界带是每个比较矩形的 bb 个旧小矩形里的两个,合起来面积是 2b\frac2b。S\mathcal S 里的区域少于 Mb\frac Mb 个,每个面积不到 1M\frac1M,合起来不到 1b\frac1b。所以 WSW_{\mathcal S} 是一个固定的集合,面积至多 3b\frac3b。
  5. 一组固定的点全落进去。如果 ESE_{\mathcal S} 发生,就有 ⌈n2⌉\lceil\frac n2\rceil 个下标,它们对应的点全在 WSW_{\mathcal S} 里。nn 个点独立、均匀,所以指定的一组点全落进固定的 WSW_{\mathcal S} 的概率至多是 (3b)⌈n/2⌉≤(3b)n/2\bigl(\frac3b\bigr)^{\lceil n/2\rceil}\le\bigl(\frac3b\bigr)^{n/2}(因为 3b≤1\frac3b\le1)。下标集合至多 2n2^n 个,由并集界,Pr⁡[ES]≤2n(3b)n/2\Pr[E_{\mathcal S}]\le2^n\bigl(\frac3b\bigr)^{n/2}。
  6. 对 S\mathcal S 用并集界。内部区域一共 MM 个,所以 S\mathcal S 至多有 2M2^M 种选法: Pr⁡[ Gj 不发生 ] ≤ 2M+n(3b)n/2 ≤ 4n(3b)n/2=(48b)n/2.\Pr[\,G_j\ \text{不发生}\,]\ \le\ 2^{M+n}\Bigl(\frac3b\Bigr)^{n/2}\ \le\ 4^n\Bigl(\frac3b\Bigr)^{n/2}=\Bigl(\frac{48}b\Bigr)^{n/2}. 中间用了 M≤nM\le n。再对 hh 层用一次并集界,就是 Pr⁡[ G 不发生 ]\Pr[\,G\ \text{不发生}\,] 的界。
边界带 各比较矩形最低、最高的旧小矩形 面积合计 2/b 重的内部区域 少于 M/b 个,面积合计不到 1/b 不重的内部区域 好层:至少一半的点在这里 一层的全部比较矩形(示意)
推导的第 1 到 4 步。灰色和红色两部分合起来面积至多 3b\frac3b;不是好层,就要求一半以上的点落在这两部分里。哪些区域是红的事先不知道,所以第 6 步要对所有可能的选法用并集界。
交互 · 18重的区域和好层

取 b=64b=64,一层有 M=4096M=4096 个比较矩形。左图的每个小方格代表一个内部区域,颜色越深点越多,红色是重的。右边的横条把 nn 个点分成三份:在不重的内部区域里、在重的内部区域里、在边界带里。青色那一段过半就是好层。均匀撒点时,默认的这组样本是好层;一般情形下,引理 2.5 保证的是以极高的概率是好层。用滑块把一部分点硬挤进 16 个区域,或者硬塞进边界带,看要挤到什么程度才不再是好层。

试一下三个预设。均匀撒点的这组样本里没有重的区域,97%97\% 左右的点在不重的内部区域里。挤 40% 的点进 16 个区域:这 16 个区域都变红了,但剩下的点仍然过半,还是好层,δj\delta_j 仍然不小于 δ∗\delta_*。挤进去的比例过了一半左右就不再是好层,第三个预设取的是 60%。要让随机撒的点自己变成这样,概率不超过最后一个读数,b=64b=64、n=Mn=M 时是 2−8502^{-850} 左右。

有三点值得停一下。第一,这里随机的是点的位置,而“好层”只和位置有关。所以 GG 对 2n2^n 种染色是同一个事件,只需要排除一次。第二,M≤nM\le n 这个条件用在第 6 步:S\mathcal S 的选法有 2M2^M 种,要让 (3b)n/2\bigl(\frac3b\bigr)^{n/2} 压得住,MM 不能比 nn 大。第三,各层的内部区域形状不同,各层是不是好层并不独立;证明只对 hh 层用了并集界,用不着独立性。

§ 18 · 核心推导收尾:定理 1.1

定理 1.1 的陈述在第 8 节。

三步合起来:命题 2.7

第 16 节计划里的 (a)(b)(c) 都有了。对一种染色 χ\chi,考虑两个事件:

事件一GG 不发生:有某一层不是好层。这件事只和点的位置有关,与染色无关。M≤nM\le n 时概率至多 h(48b)n/2h\bigl(\frac{48}b\bigr)^{n/2}(引理 2.5)。
事件二染色 χ\chi 的波动项 T≤−λT\le-\lambda,其中 λ=14hδ∗\lambda=\frac14h\delta_*。这件事和染色有关。概率由引理 2.6 控制。
命题 2.7(对所有染色同时成立的概率界)论文 Proposition 2.7

设 M≤nM\le n。则

Pr⁡[ 存在 χ 使 ∥Fχ∥∞≤h64 b3/2nM ] ≤ h(48b)n/2+2nexp⁡(−hM2048 b3).(14)\Pr\Bigl[\,\text{存在}\ \chi\ \text{使}\ \norm{F_\chi}_\infty\le\frac h{64\,b^{3/2}}\sqrt{\frac nM}\,\Bigr]\ \le\ h\Bigl(\frac{48}b\Bigr)^{n/2}+2^n\exp\Bigl(-\frac{hM}{2048\,b^3}\Bigr).\tag{14}
推导命题 2.7

左边的门槛就是 18hδ∗\frac18h\delta_*。先固定一种染色,并设 GG 发生,而 ∥F∥∞≤18hδ∗\norm F_\infty\le\frac18h\delta_*。

  1. 条件增益之和。在 GG 上,每一层的条件增益至少是 δj\delta_j(引理 2.4),而 δj≥δ∗\delta_j\ge\delta_*(引理 2.5 (i))。由式 (10) 和 Z0≥0Z_0\ge0, Zh ≥ h δ∗+T.Z_h\ \ge\ h\,\delta_*+T.
  2. FF 很小就逼出事件二。由式 (9),Zh≤2∥F∥∞≤14hδ∗Z_h\le2\norm F_\infty\le\frac14h\delta_*。和第 1 步合起来,T≤−34hδ∗T\le-\frac34h\delta_*,更有 T≤−λT\le-\lambda,λ=14hδ∗\lambda=\frac14h\delta_*。论文只用后者,常数没有优化。
  3. 事件二的概率。由 (12),δ∗2=n64 b3M\delta_*^2=\frac n{64\,b^3M}。在 (13) 里代入这个 λ\lambda: λ2M22hn=h δ∗2M232 n=hM2048 b3.\frac{\lambda^2M^2}{2hn}=\frac{h\,\delta_*^2M^2}{32\,n}=\frac{hM}{2048\,b^3}. 所以 Pr⁡[ G 发生且 ∥F∥∞≤18hδ∗ ]≤exp⁡(−hM2048b3)\Pr\bigl[\,G\ \text{发生且}\ \norm F_\infty\le\frac18h\delta_*\,\bigr]\le\exp\bigl(-\frac{hM}{2048b^3}\bigr)。
  4. 并集界。如果存在一种染色使 ∥Fχ∥∞≤18hδ∗\norm{F_\chi}_\infty\le\frac18h\delta_*,那么要么 GG 不发生,要么 GG 发生而这种染色的 FχF_\chi 很小。前一个事件只有一个;后一类事件每种染色有一个,共 2n2^n 个: Pr⁡[∃χ: ∥Fχ∥∞≤18hδ∗]≤Pr⁡[ G 不发生 ]+∑χPr⁡[ G 发生且 ∥Fχ∥∞≤18hδ∗].\Pr\bigl[\exists\chi:\ \norm{F_\chi}_\infty\le\tfrac18h\delta_*\bigr]\le\Pr[\,G\ \text{不发生}\,]+\sum_\chi\Pr\bigl[\,G\ \text{发生且}\ \norm{F_\chi}_\infty\le\tfrac18h\delta_*\bigr]. 代入引理 2.5 (ii) 和第 3 步。

右边的两项,只有第二项乘了 2n2^n。事件一是所有染色共用的,排除一次就够了。

选参数

现在把 bb、hh 定下来,目标是让 (14) 的右边不超过 e−Ane^{-An},同时让左边的门槛至少是 cA(log⁡2n)3/2c_A(\log_2n)^{3/2}。给定 A>0A>0,取整数 bb 大到满足两个条件(ln⁡\ln 是自然对数,论文写作 log⁡\log):

12ln⁡b48 ≥ A+2,b4096 ≥ A+2+ln⁡2.(15)\frac12\ln\frac b{48}\ \ge\ A+2,\qquad\frac b{4096}\ \ge\ A+2+\ln2.\tag{15}

这样的 bb 只依赖 AA,而且 b>48b>48,所以 64≤b464\le b^4。再取

cA=b−8(1+log⁡2b)3/2.c_A=\frac{b^{-8}}{(1+\log_2b)^{3/2}}.

对每个 n≥1n\ge1,取 hh 为满足 b−5h bh−1≤nb^{-5}h\,b^{h-1}\le n 的最大正整数,并令 M=bh−1M=b^{h-1}。h=1h=1 时左边是 b−5≤1b^{-5}\le1,而左边随 hh 增大趋于无穷,所以这样的 hh 存在。

参数之间的关系论文式 (16)、(17)、(18)
b−5hM ≤ n < 2 b−4hM,(16)b^{-5}hM\ \le\ n\ <\ 2\,b^{-4}hM,\tag{16} log⁡2n < (1+log⁡2b) h,(17)\log_2n\ <\ (1+\log_2b)\,h,\tag{17} cA(log⁡2n)3/2 < b−8h3/2.(18)c_A(\log_2n)^{3/2}\ <\ b^{-8}h^{3/2}.\tag{18}
推导式 (16)、(17) 和 (18)

(16) 左边的不等式是 hh 的定义。hh 是最大的,所以 h+1h+1 不满足条件:n<b−5(h+1)bhn<b^{-5}(h+1)b^h。用 h+1≤2hh+1\le2h 和 bh=bMb^h=bM,得到 n<2h b−4Mn<2h\,b^{-4}M,这是 (16) 的右边。对同一个不等式用 h+1≤2hh+1\le2^h 和 b−5≤1b^{-5}\le1:n<(2b)hn<(2b)^h,取以 2 为底的对数就是 (17)。把 (17) 两边取 32\frac32 次方,再乘上 cAc_A,就是 (18)。

式 (16) 说,一个比较矩形里的平均点数 nM\frac nM 在 b−5hb^{-5}h 和 2b−4h2b^{-4}h 之间,是 hh 的量级。式 (17) 说 h>log⁡2n1+log⁡2bh>\frac{\log_2n}{1+\log_2b},至少是 log⁡n\log n 的量级;反过来 hh 也不会更大,因为它的定义要求 b−5h bh−1≤nb^{-5}h\,b^{h-1}\le n。

为什么是 b−5b^{-5}?(16) 的下界让 (14) 的门槛至少是 h3/264b4\frac{h^{3/2}}{64b^4}。上界让 (14) 最后那个指数超过 bn4096\frac{bn}{4096},bb 够大时压得住 2n2^n。如果换成 b−4b^{-4},这个指数只能保证超过 n4096\frac n{4096},对付不了 2n2^n 种染色。

为什么密度 nM\frac nM 取 hh 的量级?密度再大一些,每个 δj\delta_j 会更大。但取 λ\lambda 为总增益的量级时,(13) 的指数是 hMhM 的量级(bb 固定),和密度无关;它必须超过 nln⁡2n\ln2。所以 nn 至多是 hMhM 的量级,密度至多是 hh 的量级。

定理 1.1 的证明

推导分两种情形
  1. 情形一:M>nM>n。命题 2.7 需要 M≤nM\le n,这种情形要另外处理。由 (16),b−5hM≤n<Mb^{-5}hM\le n<M,所以 h<b5h<b^5;由 (18),cA(log⁡2n)3/2<b−8⋅b15/2<1c_A(\log_2n)^{3/2}<b^{-8}\cdot b^{15/2}<1。于是只要找一个锚定矩形,在任何染色下差都是 1。以概率 1 各点的高度互不相同,最低的点只有一个,锚定矩形 [0,1]×[0,min⁡iVi][0,1]\times[0,\min_iV_i] 里只有它,颜色和是 ±1\pm1。定理在这种情形下以概率 1 成立。
  2. 情形二:M≤nM\le n。式 (14) 右边的第一项。这时 h≤2h−1≤M≤nh\le2^{h-1}\le M\le n。由 (15) 的第一个条件,(48b)n/2=exp⁡(−n2ln⁡b48)≤e−(A+2)n\bigl(\frac{48}b\bigr)^{n/2}=\exp\bigl(-\frac n2\ln\frac b{48}\bigr)\le e^{-(A+2)n},所以第一项至多是 n e−(A+2)nn\,e^{-(A+2)n}。
  3. 第二项。由 (16),hM>12nb4hM>\frac12nb^4,所以 hM2048 b3>bn4096≥(A+2+ln⁡2) n\frac{hM}{2048\,b^3}>\frac{bn}{4096}\ge(A+2+\ln2)\,n,用的是 (15) 的第二个条件。于是 2nexp⁡(−hM2048b3)≤enln⁡2⋅e−(A+2+ln⁡2)n=e−(A+2)n2^n\exp\bigl(-\frac{hM}{2048b^3}\bigr)\le e^{n\ln2}\cdot e^{-(A+2+\ln2)n}=e^{-(A+2)n}。
  4. 相加。(14) 的右边不超过 (n+1) e−(A+2)n=(n+1)e−2n⋅e−An≤e−An(n+1)\,e^{-(A+2)n}=(n+1)e^{-2n}\cdot e^{-An}\le e^{-An},因为 n+1≤e2nn+1\le e^{2n}。
  5. 换算成 log⁡n\log n。在这个例外事件之外,每一种染色都满足 ∥Fχ∥∞>h64 b3/2nM ≥ h64 b3/2hb5=h3/264 b4 ≥ b−8h3/2 > cA(log⁡2n)3/2.\norm{F_\chi}_\infty>\frac h{64\,b^{3/2}}\sqrt{\frac nM}\ \ge\ \frac h{64\,b^{3/2}}\sqrt{\frac h{b^5}}=\frac{h^{3/2}}{64\,b^4}\ \ge\ b^{-8}h^{3/2}\ >\ c_A(\log_2n)^{3/2}. 依次用到 (16) 的 nM≥hb5\frac nM\ge\frac h{b^5}、64≤b464\le b^4、以及 (18)。

第 5 步是指数 32\frac32 最后出现的地方:门槛是 h⋅n/Mh\cdot\sqrt{n/M},而 nM\frac nM 至少是 hh 的一个常数倍,所以门槛是 h3/2h^{3/2} 的量级。第 4 步解释了 (15) 里为什么是 A+2A+2 而不是 AA:多出来的 e−2ne^{-2n} 用来吃掉因子 n+1n+1,它来自第一项前面的 h≤nh\le n 和第二项。

交互 · 19指数账本:式 (14) 的两项

把式 (14) 右边出现的每个因子都写成 e( ⋅ )ne^{(\,\cdot\,)n} 的形式,只看括号里的数。两个事件各记一本账。事件一的两行代价来自引理 2.5 的证明:S\mathcal S 的 2M≤2n2^M\le2^n 种选法,和下标集合的 2n2^n 种选法。事件二要对 2n2^n 种染色用并集界,所以有一行 ln⁡2\ln2。两个小计都不超过 −(A+2)-(A+2) 时,证明的第 2 到 4 步就走得通。第一项前面的因子 hh 不是指数级的,没有记进账里。为了方便拖动,这里取 b=2Lb=2^L;论文只要求 bb 是整数。

两本账各有代价和收益。事件一的收益是 12ln⁡b3\frac12\ln\frac b3,随 bb 只按对数增长;事件二的收益是 b4096\frac b{4096},随 bb 线性增长,但起步晚。AA 小的时候,卡住的是事件二:A=0.5A=0.5 时它要求 b≥13080b\ge13080,而事件一只要 b≥7124b\ge7124。AA 大的时候反过来,卡住的是事件一,因为对数增长得慢。两个小计都能随 bb 增大而变得任意地负,这就是定理里的 AA 可以任意大的原因:想要多小的失败概率,就把底数 bb 取多大。

交互 · 20选参数:hh、MM 和两种情形

横轴是点数,刻度是 log⁡2n\log_2n 的对数。纵轴是一个比较矩形里的平均点数 nM\frac nM。带子是式 (16) 给的范围 b−5h≤nM<2b−4hb^{-5}h\le\frac nM<2b^{-4}h,圆点是当前的 nn。圆点在水平线 nM=1\frac nM=1 以下是 M>nM>n 的情形,以上是 M≤nM\le n 的情形。为了画得下,这里取 b=2Lb=2^L,LL 只从 4 到 8。

拖动点数可以看到:hh 大约是 log⁡2nlog⁡2b\frac{\log_2n}{\log_2b},随 nn 增长得很慢。nM\frac nM 本身是锯齿状的:nn 增大时它跟着涨,hh 每加 1 它就掉到原来的 1b\frac1b;但它始终夹在带子里,而带子随 hh 缓慢上移。带子越过 1,需要 hh 达到 12b4\frac12b^4 到 b5b^5 的量级。最后一个读数展示表达式 b−8h3/2b^{-8}h^{3/2},它超过 1 还要更晚。图中 b = 16,…,256 不满足式 (15),这里只演示尺度变化,不给出定理的概率保证。定理要求的 bb 上万,图的形状和这里一样,只是横轴要拉长很多,下面“关于常数”给出了具体的数。

这就证完了定理 1.1。随机点集以正概率满足结论,而且以概率 1 各点互不相同,所以存在 nn 点集使任何染色在某个锚定矩形上的差至少是 cA(log⁡2n)3/2c_A(\log_2n)^{3/2}。和 Nikolov 的上界合起来:

Δ2(n)=Θ(log⁡3/2n).\Delta_2(n)=\Theta(\log^{3/2}n).
关于常数

证明给出的常数很小。以失败概率 2−n/22^{-n/2} 为例,即 A=12ln⁡2A=\frac12\ln2:满足 (15) 的最小的整数是 b=12451b=12451,这时 cA≈2−114.6c_A\approx2^{-114.6}。下界 cA(log⁡2n)3/2c_A(\log_2n)^{3/2} 要超过 1,需要 log⁡2n\log_2n 大于 cA−2/3≈276.4c_A^{-2/3}\approx2^{76.4}。在此之前,定理给出的界还不到 1,不如“差至少是 1”这件平凡的事。

这在渐近结果里很常见:定理确定的是增长的量级,不是具体数值。论文在定理后面直接说了:没有去优化 cAc_A,它非常小,所以这个界只对很大的 nn 有意义。

Part V · 第 19–22 节

回顾

证明走完了。这一部分先把 32\frac32 这个指数的账重新对一遍,再回答开头的三个问题,然后看还有什么没解决。最后附上依赖图、记号表和文献。

§ 19 · 进阶32\frac32 的账

证明里有三个数在起作用:

hh层数,也就是从一个阶段到下一个阶段的过渡次数。由 (17),它是 log⁡n\log n 的量级。
δ∗≍n/M\delta_*\asymp\sqrt{n/M}每过一个好层,ZZ 的条件增益至少是 δ∗=18b3/2n/M\delta_*=\frac1{8b^{3/2}}\sqrt{n/M}。nM\frac nM 是一个比较矩形里的平均点数。
nM≍h\frac nM\asymp h参数的选择,式 (16)。于是每个好层涨 h\sqrt h 的量级,hh 层都是好层,合计 h⋅h=h3/2h\cdot\sqrt h=h^{3/2}。
∥Fχ∥∞ > h δ∗8=h64 b3/2nM ≍ h3/2 ≍ (log⁡n)3/2.\norm{F_\chi}_\infty\ >\ \frac{h\,\delta_*}8=\frac h{64\,b^{3/2}}\sqrt{\frac nM}\ \asymp\ h^{3/2}\ \asymp\ (\log n)^{3/2}.

hδ∗8\frac{h\delta_*}8 里的 8 是这样来的:如果 ∥Fχ∥∞\norm{F_\chi}_\infty 不超过它,那么 Zh≤2∥Fχ∥∞≤hδ∗4Z_h\le2\norm{F_\chi}_\infty\le\frac{h\delta_*}4,离条件增益之和 hδ∗h\delta_* 还差 34hδ∗\frac34h\delta_*,这个差只能由波动项来补。论文只用了其中的 14hδ∗\frac14h\delta_*,常数没有优化。

为什么 nM\frac nM 不能取得更大

每层的涨幅是 n/M\sqrt{n/M},看起来 nM\frac nM 越大越好。限制来自并集界。式 (13) 里取 λ=hδ∗4\lambda=\frac{h\delta_*}4,指数是 hM2048b3\frac{hM}{2048b^3};要压住 2n2^n 种染色,它必须比 nln⁡2n\ln2 大。所以 hMhM 至少要是 nn 的 b3b^3 的若干倍,也就是

nM ≲ hb3.\frac nM\ \lesssim\ \frac h{b^3}.

第 18 节的取法 nM<2hb4\frac nM<\frac{2h}{b^4} 在这个范围里,还留了一个因子 bb 的余量,用来把失败概率压到 e−Ane^{-An}。论文在引言里是这样说的:密度再大,每层的增益会更大;但对总增益量级的偏离,波动的估计给出的指数是 hMhM 的量级,不管密度是多少,而它必须超过 nln⁡2n\ln2。所以密度 nM\frac nM 取 hh 的量级,已经是这个论证允许的最大值。

一句话的直观

每一层有 MM 个比较矩形,hh 层共 hMhM 个。对一种固定的染色,ZhZ_h 要异常地小,大多数比较矩形里的随机和都得碰巧抵消掉。如果把这 hMhM 个矩形当成各自独立、每个碰巧的概率是个常数,合起来大约是 e−c hMe^{-c\,hM}。染色有 2n2^n 种可挑,并集界要求 2n⋅e−c hM2^n\cdot e^{-c\,hM} 很小,所以要 hM≳nhM\gtrsim n,即 nM≲h\frac nM\lesssim h。

这只是关于尺度的启发式说法:各个比较矩形、各层之间并不独立。证明用的是第 16 节逐层迭代的条件指数矩界,也就是式 (13):取 λ\lambda 为总增益的量级时,它右边的指数正好是 hMhM 的常数倍。

和连续差异度对比一下。第 6 节的二进网,每个二进条带的误差不超过 1,log⁡n\log n 层合计 log⁡n\log n。这里每一层的涨幅是 log⁡n\sqrt{\log n} 的量级,合计 log⁡3/2n\log^{3/2}n。两个问题相差的 log⁡n\sqrt{\log n},就是每层多出来的这个平方根。

和第 9 节的猜测对账

第 9 节用布朗单的小球估计猜出了尺度。现在把证明里的参数翻译成那边的语言。差异度的尺度是 t≍hn/Mt\asymp h\sqrt{n/M},对应的小球半径是 ε=t/n≍h/M\varepsilon=t/\sqrt n\asymp h/\sqrt M。

布朗单的小球估计这里的证明
半径 ε\varepsilonh/Mh/\sqrt M
log⁡(1/ε)\log(1/\varepsilon)12log⁡M−log⁡h\frac12\log M-\log h,是 hh 的量级(bb 是常数)
指数 ε−2log⁡3(1/ε)\varepsilon^{-2}\log^3(1/\varepsilon)Mh2⋅h3=hM≍n\frac M{h^2}\cdot h^3=hM\asymp n

最后一行说,在这组参数下小球概率的指数正好是 nn 的量级,也就是刚好够付 2n2^n 的并集界。两边用的是同一个“三次方”:n≍ε−2h3n\asymp\varepsilon^{-2}h^3。证明没有用到布朗单,但落在了它预言的尺度上。论文第 3 节做的是第 9 节那个计算;把证明里的参数代回去对账,是这份导读加的。

§ 20 · 入门之前为什么不行,这次为什么可以

之前为什么不行

  • 借连续差异度。转移原理只能把连续差异度的下界搬过来,而连续差异度恰好是 Θ(log⁡n)\Theta(\log n)。单靠这条路,到 log⁡n\log n 为止。
  • 正交函数法。在第 7 节介绍的经典论证里,检验函数的收益按尺度逐层相加,log⁡n\log n 层各贡献一个常数。
  • 分解范数。γ2\gamma_2 与差异度之间的两个通用不等式各带对数因子的损失,平面上的下界一侧正好落在 log⁡n\log n。

这些具体论证没有给出本文对随机点集的统一下界。这里介绍的正交函数法在组合差异度上的应用针对二进网,另外两条路线给出的坏点集藏在某个网格的子集里。二进网本身就能染到 O(log⁡n)O(\log n)(见第 7 节),在它上面不可能证出更多。

这次为什么可以

换点集点随机撒。二进网这样的点集可以让各层尺度的误差同时很小;随机点的每一位数字都是新掷的骰子,每一层都带来新的起伏。
换问法一次固定一种染色,估计它“染得好”的概率,再对 2n2^n 种染色用并集界。这要求概率小到 e−cne^{-cn},而且 cc 要足够大。
换一个量不直接盯着 FF 的最大值,改看局部振荡 ZjZ_j。在好层上,它的条件期望增益至少为 δ∗=(8b3/2)−1n/M\delta_*=(8b^{3/2})^{-1}\sqrt{n/M}。固定已揭示的信息,只改变一个点新揭示的下一位横坐标数字,Zj+1Z_{j+1} 至多变化 1/M1/M。
初等的概率Cauchy–Schwarz(论文用 Hölder)、有界差分不等式、Markov 不等式和并集界。全程只用到阶梯函数的最大最小值和条件平均。论文摘要里说,这个证明“出人意料地简单和初等”。
旧方法本文
点集二进网,或者某个网格的子集随机撒点
怎样对付全部 2n2^n 种染色对每种染色给出证据(检验函数、分解范数)对每种染色估计“染得好”的概率,再用并集界
需要的关键估计检验函数的内积、矩阵的分解范数局部振荡的平均涨幅和单点影响
各层尺度的贡献常数,合计 log⁡n\log nlog⁡n\sqrt{\log n},合计 log⁡3/2n\log^{3/2}n

和旧方法的联系也值得一提。“横向变细、纵向变粗、面积不变”的矩形族,正是 Roth 和 Halász 的正交函数法里用的那一族。论文在相关工作里还指出了和 Talagrand 证明布朗单小球估计的联系:Talagrand 的证明里有一个确定性的递推,在一个坐标方向上取平均,在另一个方向上挑代表点;论文的局部比较用的是同样的一对相反的尺度。对子函数取平均,在期望意义下保住了旧函数;内部区域里的点带来的独立对称增量提供了增益。端点平均 aIa_I 则是 Talagrand 那里“横向取平均”的条件期望版本(第 13 节末尾)。

论文的四步,和这份导读的对应

论文引言把证明概括成四步。第四部分讲的就是这四步,只是顺序稍有不同:这里把“一位数字的影响”提前到了第 15 节,把好层放在了波动之后。

论文的四步内容论文里的位置这份导读
(1) Potential逐位揭示横坐标;端点平均、振荡、ZjZ_j;Zj≤2∥F∥∞Z_j\le2\norm F_\infty式 (1)–(3)、(9)§ 10–§ 12
(2) Local gain一个比较矩形里,竖条上的振荡之和平均比横条上多 14mR\frac14\sqrt{m_R};对所有比较矩形求和引理 2.1、Fact 2.2、引理 2.3、引理 2.4§ 13、§ 14、§ 16
(3) Good transitions重的区域、好层;M≤nM\le n 时每一层都是好层,例外的概率指数小,而且与染色无关引理 2.5§ 17
(4) All colorings一位数字只让 Zj+1Z_{j+1} 变 1M\frac1M;有界差分控制波动项;对 2n2^n 种染色用并集界;选参数引理 2.6、Fact A.1、命题 2.7、2.5 小节(Parameters)§ 15、§ 16、§ 18

技术清单

技术作用位置
随机点,逐个染色的并集界把“每种染色都不好”化成一个固定染色的概率估计§ 8
逐位揭示横坐标每一层只放出 nn 个独立的骰子§ 10
端点平均 aIa_I(式 (2))细分之后,子函数的平均在期望意义下回到旧函数§ 12、§ 14
局部振荡 ZjZ_j(式 (1))不怕加常数:对每个受影响的横向区间,只有包含该点高度的块可能改变振荡§ 12、§ 15
纵向比较(引理 2.1,式 (4))JJ 上的振荡看得见首尾两个块之间的偏移,单个块上的振荡看不见§ 13
先平均,再比较(式 (5)、(6))把 bb 个子函数并成一个 aˉ\bar a,它等于旧函数加上各点的增量§ 13
对称增量之和不会扎堆(Fact 2.2)内部点的增量之和 DD 加上任何固定的偏移,绝对值的期望至少是 18mR\frac18\sqrt{m_R}§ 14
条件增益加波动项(式 (10)、引理 2.4)把 ZhZ_h 写成初值、各层条件增益之和与均值为 0 的波动项;条件增益有一个与染色无关的共同下界§ 16
有界差分(Fact A.1、引理 2.6)波动项很难是大的负数,不需要各层独立§ 15、§ 16
重的区域与好层(引理 2.5)M ≤ n 时,各层以高概率同时是好层;这个事件与染色无关§ 17

§ 21 · 入门常数、形式化验证与开放问题

显式的坏点集

证明是概率性的:它说明随机点集几乎都是坏的,却指不出一个具体的。能不能明确写出一族点集,差异度是 Ω(log⁡3/2n)\Omega(\log^{3/2}n)?论文指出,主要的困难是:把独立的随机数字换成确定的选择之后,怎样保住在许多层尺度上累积起来的振荡。

高维

对 d≥3d\ge3 维的盒子,已知 Ωd(log⁡d−1n)≤Δd(n)≤Od(log⁡d−1/2n)\Omega_d(\log^{d-1}n)\le\Delta_d(n)\le O_d(\log^{d-1/2}n),缺口同样是 log⁡n\sqrt{\log n}。(第 7 节里的 log⁡d+1/2n\log^{d+1/2}n 是分解范数方法给出的较弱的上界,Nikolov 后来把它改进到了 log⁡d−1/2n\log^{d-1/2}n。)Bansal 猜想 dd 维的答案是 O(log⁡d−1n)O(\log^{d-1}n);论文的定理说明这在 d=2d=2 时不对。

高维该是哪一边?把第 9 节的计算搬到 dd 维:布朗单的小球猜想预言指数是 ε−2log⁡2d−1(1/ε)\varepsilon^{-2}\log^{2d-1}(1/\varepsilon),和 nln⁡2n\ln2 平衡之后得到 log⁡d−1/2n\log^{d-1/2}n,与已知的上界吻合。论文同时强调,这个启发式的计算并不能把结论转移到有限点集上。证明方面,论文的论证用到了横向和纵向两个相反的分辨率:细分一个坐标时,一个点的影响被限制在一个纵向块里(第 15 节)。要推广到高维,需要在几个坐标同时细分或合并时找到相应的局部比较。

常数

第 18 节末尾算过,证明给出的 cAc_A 极小:A=12ln⁡2A=\frac12\ln2 时 b=12451b=12451,cA≈2−114.6c_A\approx2^{-114.6},结论要到天文数字般的 nn 才有实质内容。论文明确说没有优化 cAc_A。

形式化验证与 AI

论文末尾有一段“AI methodology statement”,大意是:作者从 2013 年起就在考虑这个问题;2026 年 7 月开始,作者与 GPT 和 Claude 交互式地研究它,先得到一个基于“killed suffix walk”的证明,后来简化成论文第 2 节和附录里的局部振荡证明;同样的工具加上 Xean 被用来检查证明和准备稿件;主要的下界定理有一份 Lean 形式化,按论文的说法由 GPT 完成;作者核实了全部内容的正确性和原创性。形式化的代码在 github.com/wzskytop/tusnady-lean,所用工具和核查记录以仓库当前的说明为准;论文在 arXiv:2610.08130。

形式化的范围要分清。它覆盖的是论文的新下界:随机点的定理 1.1,以及由它得到的确定性下界。引用的上界、由此合成的 Θ(log⁡3/2n)\Theta(\log^{3/2}n)、历史叙述,以及这份导读补充的解释和交互演示,都不在形式化的范围里。

定理、引理与 Lean 声明的对应关系见仓库的 论文与代码对应表。

这份导读补充了推导、例子和交互演示。它们不能代替形式化验证,也不能代替独立的审阅。

其他方向

论文在相关工作里提到两个方向。一是在线差异问题:点随机到达时,已有工作研究了区间、条带和 Tusnády 问题的在线算法。对 [0,1]d[0,1]^d 里均匀随机到达的点,已知任何在线算法在前 TT 个点里以常数概率出现 Ωd(log⁡d/2T)\Omega_d(\log^{d/2}T) 的差(Bansal–Jiang–Singla–Sinha)。定理 1.1 对任何染色都成立,在线、离线都算,所以在平面上把这个界从 log⁡T\log T 抬到了 log⁡3/2T\log^{3/2}T。二是数据结构:差异度和范围搜索、近似范围计数的下界有关,而本文的难例是随机的,或许能用来得到平均情形的下界。

§ 22 · 附录依赖图、记号表与文献

证明的依赖关系

交互 · 21谁用到了谁

箭头从被引用的结论指向用到它的结论。点击任一方框跳到对应的小节。

记号表

记号含义首次出现
P, nP,\ n点集和点数§ 1
χ\chi染色,取值 ±1\pm1;χ(P∩R)\chi(P\cap R) 是 RR 内颜色之和§ 3
disc⁡(P), Δ2(n)\disc(P),\ \Delta_2(n)点集的差异度;nn 个点的最坏差异度§ 3
D(P), D(n)D(P),\ D(n)连续差异度§ 6
Pi=(Ui,Vi)P_i=(U_i,V_i)第 ii 个随机点的横坐标和纵坐标§ 8
Fχ(x,y), ∥Fχ∥∞F_\chi(x,y),\ \norm{F_\chi}_\infty锚定矩形 [0,x]×[0,y][0,x]\times[0,y] 里的颜色和;它的最大绝对值。第 11 节起简写为 FF§ 8
A, cAA,\ c_A失败概率 e−Ane^{-An} 里的指数;定理 1.1 里的常数§ 8
W, εW,\ \varepsilon布朗单;小球的半径§ 9
b, Lb,\ L底数,整数,b≥3b\ge3;定理的证明里取得很大(式 (15))。第 18 节的交互里取成 2 的幂 b=2Lb=2^L§ 10、§ 18
Fj\mathcal F_j第 jj 阶段的已知信息:全部纵坐标和每个横坐标的前 jj 位数字§ 10
h, Mh,\ M层数:阶段编号 j=0,…,hj=0,\dots,h,过渡共 hh 次,从阶段 jj 到 j+1j+1 叫第 jj 层(论文:transition j→j+1j\to j+1);M=bh−1M=b^{h-1} 是每层比较矩形的个数§ 11
I, BI,\ B第 jj 阶段的横向区间(宽 b−jb^{-j})和纵向块(高 bj−hb^{j-h});I×BI\times B 是小矩形(论文:cell)§ 11
Ck, J, RC_k,\ J,\ RII 的第 kk 个子区间(论文:child);由 bb 个旧块 B1,…,BbB_1,\dots,B_b 并成的新块;比较矩形 I×JI\times J§ 11
aI(y), wI(u)a_I(y),\ w_I(u)端点平均,FF 在 I=[s,t]I=[s,t] 的两个端点处的值的平均,式 (2);论文写作 a(I)a^{(I)}。横坐标是 uu 的点在其中的系数(权重)wI(u)w_I(u) 是 1,12,01,\frac12,0§ 12
osc⁡B f, mid⁡B f\operatorname{osc}_B\,f,\ \operatorname{mid}_B\,f振荡:ff 在 BB 上的最大值减最小值;中点:最大值和最小值的平均。式 (3)§ 12
ZjZ_j局部振荡(论文:potential):第 jj 阶段各小矩形上 aIa_I 的振荡的平均,式 (1)§ 12
Rin, mRR^{\mathrm{in}},\ m_RRR 的内部区域 I×(B2∪⋯∪Bb−1)I\times(B_2\cup\dots\cup B_{b-1});里面的点数§ 13
a, ak, aˉa,\ a_k,\ \bar a固定 II 时的简写:旧函数 a=aIa=a_I;子函数 ak=aCka_k=a_{C_k};它们的平均 aˉ=1b∑kak\bar a=\frac1b\sum_ka_k§ 13
偏移mid⁡Bbf−mid⁡B1f\operatorname{mid}_{B_b}f-\operatorname{mid}_{B_1}f:最高的块相对最低的块(论文:offset)§ 13
r, Xir,\ X_i骰子的结果(点落进 CrC_r);这个点的增量 χib+1−2r2b\chi_i\frac{b+1-2r}{2b}(论文:jump)§ 13
D, g, θD,\ g,\ \theta内部点的增量之和;aˉ\bar a 去掉这些增量之后的函数;gg 自己的偏移§ 14
σ2\sigma^2一个增量的二阶矩 b2−112b2\frac{b^2-1}{12b^2}§ 14
条件增益E[Zj+1−Zj∣Fj]\EE[Z_{j+1}-Z_j\mid\mathcal F_j]:第 jj 层平均涨多少§ 16
T, TjT,\ T_j波动项 ∑j(Zj+1−E[Zj+1∣Fj])\sum_j\bigl(Z_{j+1}-\EE[Z_{j+1}\mid\mathcal F_j]\bigr),式 (10);其中第 jj 层的那一项§ 16
δj\delta_j第 jj 层条件增益的下界 14bM∑RmR\frac1{4bM}\sum_R\sqrt{m_R},式 (11)§ 16
Γ, ξi, φ, c, t\Gamma,\ \xi_i,\ \varphi,\ c,\ tFact A.1 里的已知信息、独立变量、函数、单个变量的影响的上界、指数里的参数§ 16
λ\lambda引理 2.6 里偏离的大小;命题 2.7 里取 14hδ∗\frac14h\delta_*§ 16、§ 18
重、好层mR>bnMm_R>\frac{bn}M 的内部区域是重的;至少 n2\frac n2 个点在不重的内部区域里的层是好层§ 17
Gj, GG_j,\ G第 jj 层是好层;hh 层都是好层§ 17
δ∗\delta_*好层保证的下界 18b3/2n/M\frac1{8b^{3/2}}\sqrt{n/M},式 (12)§ 17
边界带,S, WS, ES\mathcal S,\ W_{\mathcal S},\ E_{\mathcal S}各比较矩形最低、最高的旧小矩形;少于 Mb\frac Mb 个内部区域的一组;这一组和边界带的并;“至少一半的点在 WSW_{\mathcal S} 里”这个事件。第 17 节里不带下标的 WW 是边界带和全部重区域的并,和第 9 节的布朗单无关§ 17

文献

以下书目信息取自论文 r65 版的参考文献表;标“背景”的两条是这份导读补充的。

  • Z. Wei. Tight Bounds for Tusnády's Problem in the Plane. arXiv:2610.08130. 本导读讲解的论文;导读依据的是 r65 版稿件,2026-10-07 已确认其 TeX 源文件与 arXiv v1 逐字节相同。
  • C. Aistleitner, D. Bilyk, and A. Nikolov. Tusnády's problem, the transference principle, and non-uniform QMC sampling. MCQMC 2016(Springer,2018).
  • K. Azuma. Weighted sums of certain dependent random variables. Tôhoku Mathematical Journal 19(3), 1967.
  • W. Banaszczyk. On series of signed vectors and their rearrangements. Random Structures & Algorithms 40(3), 2012.
  • N. Bansal. Discrepancy theory and related algorithms. Proceedings of the International Congress of Mathematicians 2022, Vol. 7. EMS Press, 2023.
  • N. Bansal and S. Garg. Algorithmic discrepancy beyond partial coloring. STOC 2017.
  • N. Bansal, H. Jiang, R. Meka, S. Singla, and M. Sinha. Online discrepancy minimization for stochastic arrivals. SODA 2021.
  • N. Bansal, H. Jiang, S. Singla, and M. Sinha. Online vector balancing and geometric discrepancy. STOC 2020.
  • N. Bansal and A. Nikolov. Discrepancy theory: an algorithmic and geometric perspective. arXiv:2608.00140, 2026.
  • R. F. Bass. Probability estimates for multiparameter Brownian processes. The Annals of Probability 16, 1988.
  • J. Beck. Balanced two-colorings of finite sets in the square I. Combinatorica 1(4), 1981.
  • J. Beck. Balanced two-colorings of finite sets in the cube. Discrete Mathematics 73, 1989.
  • D. Bilyk and M. Lacey. The supremum norm of the discrepancy function: recent results and connections. MCQMC 2012(Springer,2013).
  • G. Bohus. On the discrepancy of 3 permutations. Random Structures & Algorithms 1(2), 1990.
  • G. Halász. On Roth's method in the theory of irregularities of point distributions. Recent Progress in Analytic Number Theory, Vol. 2, 1981.
  • H. Jiang, J. Kulkarni, and S. Singla. Online geometric discrepancy for stochastic arrivals with applications to envy minimization. arXiv:1910.01073, 2019.
  • K. G. Larsen. On range searching in the group model and combinatorial discrepancy. SIAM Journal on Computing 43(2), 2014.
  • J. Matoušek. Geometric Discrepancy: An Illustrated Guide. Springer, 1999.
  • J. Matoušek and A. Nikolov. Combinatorial discrepancy for boxes via the γ2\gamma_2 norm. SoCG 2015.
  • J. Matoušek, A. Nikolov, and K. Talwar. Factorization norms and hereditary discrepancy. IMRN, 2020.
  • C. McDiarmid. On the method of bounded differences. Surveys in Combinatorics, 1989. Cambridge University Press, 1989.
  • A. Nikolov. Tighter bounds for the discrepancy of boxes and polytopes. Mathematika 63(3), 2017.
  • W. M. Schmidt. Irregularities of distribution. VII. Acta Arithmetica 21, 1972.
  • A. Srinivasan. Improving the discrepancy bound for sparse matrices: better approximations for sparse lattice approximation problems. SODA 1997.
  • M. Talagrand. The small ball problem for the Brownian sheet. The Annals of Probability 22(3), 1994.
  • Z. Wei and K. Yi. The space complexity of 2-dimensional approximate range counting. SODA 2013. 期刊版:Tight space bounds for two-dimensional approximate range counting, ACM Transactions on Algorithms, 2018.
  • C. Xu. Xean. GitHub 仓库,2026,github.com/chaoxu/xean.
  • 背景:K. F. Roth. On irregularities of distribution. Mathematika 1, 1954.
  • 背景:第 6、7 节关于二进网、转移原理和 γ2\gamma_2 范数的叙述,依据上面 Matoušek 的专著以及 Matoušek–Nikolov–Talwar、Nikolov 的论文。