组合与概率

不等式证明中的图论方法:团数、谱与路径结构的应用

正文文字 6609 字公式 973 处

整理了昨天比较多人想看的“不等式证明中的图论方法”,联赛之前会优先整理代数板块的专题(手上的资料和做过的题目都相对更多一些).

本系列文章均由我提供完整思路、所有例题,并交给GPT整理排版而成.


本文讨论一种常见的建模视角:当一个多元不等式中的项由“哪两个指标发生关系”决定时,可以把指标看成顶点,把相应关系看成边.由此,原来的代数问题往往可转化为团数、独立集、路径、哈密顿结构或矩阵谱的估计.图论在此提供了一种整理“哪些项出现、哪些项不能同时出现”等离散关系的有效语言.

除非另有说明,本文中的图均指有限简单图;𝜔(𝐺)表示图 𝐺的团数,即最大完全子图的顶点数;𝐴(𝐺)与 𝐿(𝐺)分别表示邻接矩阵和图拉普拉斯矩阵.下标按模 𝑛理解时,会在相应题目中另行说明.

从约束形式选择工具

不同归一化条件通常对应不同工具:

  1. 𝑥_(𝑣)≥0且 ∑_(𝑣)𝑥_(𝑣)固定:优先考虑边权和、团数、着色和 Motzkin–Straus 定理;

  2. ∑_(𝑣)𝑥2𝑣固定:优先考虑邻接矩阵、图拉普拉斯与 Rayleigh 商;

  3. 出现局部比值且有明确的“向前候选”:优先考虑有向路径与乘积链;

  4. 条件涉及任意排列中的相邻项:优先考虑哈密顿路、哈密顿圈及其平均.

关键在于使图的定义准确反映原式:先确定“顶点对应什么、边表示什么”,再确定所需控制的是团、路径、谱还是平均出现次数.

一、团数与边权和:Motzkin–Straus 定理

“

方法要点 若目标式为若干非负乘积 𝑥_(𝑢)𝑥_(𝑣)之和,且 ∑_(𝑣)𝑥_(𝑣)固定,可将出现的指标对作为边.问题随即转为顶点加权图的边权和极值;团数、着色以及等权情形下的边数都可作为控制参数.

例 1.1|Motzkin–Straus 定理(加权形式)[1]

设 𝐺=(𝑉,𝐸)的团数为 𝜔=𝜔(𝐺).非负实数 𝑥_(𝑣) (𝑣∈𝑉)满足∑_(𝑣∈𝑉)𝑥_(𝑣)=1.证明:

∑_(𝑢𝑣∈𝐸)𝑥_(𝑢)𝑥_(𝑣)≤(𝜔−1)/(2𝜔),

并说明等号可以达到.

分析 在一组使目标函数达到最大值的权重中,若有两个不相邻的正权顶点,可将其中一个顶点的全部权重合并到另一个顶点,并选择使目标函数不减的合并方向.反复进行这一操作后,正权顶点组成一个团,问题化为固定和条件下完全图的二次型估计.

证明 记

𝐹(𝑥)=∑_(𝑢𝑣∈𝐸)𝑥_(𝑢)𝑥_(𝑣).

由于可行集闭且有界,而 𝐹为连续函数,故最大值存在.在所有达到最大值的权重中,取正权顶点数最少的一组.若正权顶点 𝑢,𝑣不相邻,令 𝐴_(𝑢)=∑_(𝑢𝑤∈𝐸)𝑥_(𝑤), 𝐴_(𝑣)=∑_(𝑣𝑤∈𝐸)𝑥_(𝑤).由于 𝑢𝑣∉𝐸,把 𝑥_(𝑣)全部并入 𝑥_(𝑢)时不会额外产生 𝑥_(𝑢)𝑥_(𝑣)项.不妨设 𝐴_(𝑢)≥𝐴_(𝑣).将 𝑥_(𝑢)改为 𝑥_(𝑢)+𝑥_(𝑣)、将 𝑥_(𝑣)改为 0,其余权重不变,则 𝐹(𝑥^(′))−𝐹(𝑥)=𝑥_(𝑣)(𝐴_(𝑢)−𝐴_(𝑣))≥0.由于 𝐹(𝑥)已是最大值,必有 𝐹(𝑥^(′))=𝐹(𝑥);但 𝑥^(′)的正权顶点数更少,矛盾.因此正权顶点组成一个 𝑟阶团,其中 𝑟≤𝜔.于是

𝐹 =∑_(1≤𝑖<𝑗≤𝑟)𝑥_(𝑖)𝑥_(𝑗); =(1)/(2)(1−∑_(𝑖=1)^(𝑟)𝑥2𝑖); ≤(1)/(2)(1−(1)/(𝑟)); ≤(𝜔−1)/(2𝜔).;

在一个最大团的 𝜔个顶点上取 𝑥_(𝑣)=1/𝜔,其余顶点取 𝑥_(𝑣)=0,等号成立.

“

注 当 𝐺无三角形时 𝜔(𝐺)≤2.若总权重 𝑊=∑_(𝑣∈𝑉)𝑥_(𝑣)>0,将各权重除以 𝑊后应用上式,再乘回 𝑊^(2);𝑊=0时结论成立.因此对任意非负权重都有

∑_(𝑢𝑣∈𝐸)𝑥_(𝑢)𝑥_(𝑣)≤(1)/(4)(∑_(𝑣∈𝑉)𝑥_(𝑣))^(2).

特别地,𝑚≥4时圈 𝐶_(𝑚)无三角形,若约定 𝑥_(𝑚+1)=𝑥_(1),则

∑_(𝑖=1)^(𝑚)𝑥_(𝑖)𝑥_(𝑖+1)≤(1)/(4)(𝑥_(1)+⋯+𝑥_(𝑚))^(2).

下文将这一直接推论简称为 1/4引理.

补充:Turán 定理与等权极值

Motzkin–Straus 定理给出非负顶点权的边权和上界.取所有顶点权相等时,可得到禁完全图条件下的边数估计.记 𝑇_(𝑟)(𝑁)为将 𝑁个顶点尽量均匀地分成 𝑟部(允许空部)、同部内无边而异部间全连边所得的完全 𝑟部图,并记 𝑡_(𝑟)(𝑁)=|𝐸(𝑇_(𝑟)(𝑁))|.若𝑁=𝑞𝑟+𝑠,0≤𝑠<𝑟,则 𝑇_(𝑟)(𝑁)有 𝑠个大小为 𝑞+1的部和 𝑟−𝑠个大小为 𝑞的部,从而

𝑡_(𝑟)(𝑁)=(𝑟−1)/(2𝑟)𝑁^(2)−(𝑠(𝑟−𝑠))/(2𝑟).

Turán 定理 [2]

若 𝑁个顶点的图 𝐺不含 𝐾_(𝑟+1),则|𝐸(𝐺)|≤𝑡_(𝑟)(𝑁),且等号由 𝑇_(𝑟)(𝑁)取得.

证明 对参数 𝑟+𝑁作强归纳.当 𝑟=1时图中不能有边;当 𝑁≤𝑟时 𝑇_(𝑟)(𝑁)=𝐾_(𝑁),结论显然成立.以下设 𝑟≥2且 𝑁>𝑟.若 𝐺不含 𝐾_(𝑟),则对参数 (𝑟−1,𝑁)应用归纳假设,得到

|𝐸(𝐺)|≤𝑡_(𝑟−1)(𝑁)≤𝑡_(𝑟)(𝑁).

将 𝑇_(𝑟−1)(𝑁)视为一个第 𝑟部为空的完全 𝑟部图.一般地,完全 𝑟部图中若两部大小相差至少 2,将较大部的一个顶点移入较小部,边数增加;反复调整后得到平衡完全 𝑟部图.因此 𝑡_(𝑟−1)(𝑁)≤𝑡_(𝑟)(𝑁). 若 𝐺含有一个 𝑟阶团 𝑆,则每个 𝑆外顶点至多与 𝑆中 𝑟−1个顶点相邻.因此,至少有一个端点在 𝑆中的边数至多为

((𝑟)/(2))+(𝑁−𝑟)(𝑟−1)=(𝑟−1)(𝑁−(𝑟)/(2)).

图 𝐺−𝑆仍不含 𝐾_(𝑟+1),且顶点数减少了 𝑟,故对参数 (𝑟,𝑁−𝑟)应用归纳假设,得到

|𝐸(𝐺−𝑆)|≤𝑡_(𝑟)(𝑁−𝑟).

又因 𝑁与 𝑁−𝑟除以 𝑟的余数相同,直接由 𝑡_(𝑟)的公式得到

𝑡_(𝑟)(𝑁)−𝑡_(𝑟)(𝑁−𝑟)=(𝑟−1)(𝑁−(𝑟)/(2)).

将以上两部分边数相加,即得 |𝐸(𝐺)|≤𝑡_(𝑟)(𝑁).

“

注 在例 1.1 中令每个顶点权均为 1/𝑁,若 𝜔(𝐺)≤𝑟,则

|𝐸(𝐺)|≤(𝑟−1)/(2𝑟)𝑁^(2).

该式与 Turán 定理具有相同的主项,Turán 定理还给出由整数分部产生的修正项.特别地,𝑟=2时得到 Mantel 定理 |𝐸(𝐺)|≤⌊𝑁^(2)/4⌋;相应的加权结论即前述 1/4引理.本文中,Turán 定理作为 Motzkin–Straus 定理在等权计数问题中的补充.


二、关系图:由整除链确定团数

“

方法要点 整除、包含、可比等关系本身具有组合结构.若两两相邻对应一条链,可由链长估计团数;若能分成若干同组内无边的集合,则可利用着色或直接合并各组权重.

例 2.1|2013 Korean Mathematical Olympiad Final Round, Problem 3

设整数 𝑛≥2,记

𝑇={(𝑖,𝑗)∣1≤𝑖<𝑗≤𝑛, 𝑖∣𝑗}.

非负实数 𝑥_(1),𝑥_(2),…,𝑥_(𝑛)满足 𝑥_(1)+⋯+𝑥_(𝑛)=1.求

∑_((𝑖,𝑗)∈𝑇)𝑥_(𝑖)𝑥_(𝑗)

的最大值.

分析 以 1,2,…,𝑛为顶点,当两数中较小者整除较大者时连边.目标式即所得关系图的边权和.团中的整数按递增顺序排列后形成严格整除链,而严格整除链中后一项至少为前一项的 2倍.因此可先确定该图的团数,再应用例 1.1.

证明 令

𝑘=⌊log_(2)⁡𝑛⌋,2^(𝑘)≤𝑛<2^(𝑘+1).

若 𝑎_(1)<𝑎_(2)<⋯<𝑎_(𝑟)构成一个团,则

𝑎_(1)∣𝑎_(2)∣⋯∣𝑎_(𝑟).

相邻两项不同,故 𝑎_(𝑡+1)≥2𝑎_(𝑡),从而

𝑎_(𝑟)≥2^(𝑟−1)𝑎_(1)≥2^(𝑟−1).

由 𝑎_(𝑟)≤𝑛<2^(𝑘+1)得 𝑟≤𝑘+1.另一方面,1,2,4,…,2^(𝑘)构成一个 𝑘+1阶团,故该图的团数为 𝑘+1.由例 1.1,

∑_((𝑖,𝑗)∈𝑇)𝑥_(𝑖)𝑥_(𝑗)≤(𝑘)/(2(𝑘+1)).

在 1,2,4,…,2^(𝑘)上各取权重 1/(𝑘+1),其余权重取 0,等号成立.因此最大值为

(𝑘)/(2(𝑘+1)),𝑘=⌊log_(2)⁡𝑛⌋.

另一种证明:分组法

将 {1,2,…,𝑛}分成 𝑘+1组:

𝐵_(𝑡)={2^(𝑡−1),2^(𝑡−1)+1,…,2^(𝑡)−1}(1≤𝑡≤𝑘),

以及

𝐵_(𝑘+1)={2^(𝑘),2^(𝑘)+1,…,𝑛}.

同一组中若 𝑎<𝑏,则 𝑏<2𝑎;若又有 𝑎∣𝑏,则因 𝑎,𝑏为不同正整数必有 𝑏≥2𝑎,矛盾.因此同一组中任意两个不同整数均不存在整除关系.

记

𝐴_(𝑡)=∑_(𝑖∈𝐵_(𝑡))𝑥_(𝑖)(1≤𝑡≤𝑘+1).

每个 (𝑖,𝑗)∈𝑇的两个端点分属不同的组,故

∑_((𝑖,𝑗)∈𝑇)𝑥_(𝑖)𝑥_(𝑗) ≤∑_(1≤𝑠<𝑡≤𝑘+1)𝐴_(𝑠)𝐴_(𝑡); =(1)/(2)⎡⎣(∑_(𝑡=1)^(𝑘+1)𝐴_(𝑡))^(2)−∑_(𝑡=1)^(𝑘+1)𝐴2𝑡⎤⎦; ≤(1)/(2)(1−(1)/(𝑘+1)); =(𝑘)/(2(𝑘+1)).;

最后一个不等式由柯西不等式 ∑𝐴2𝑡≥1/(𝑘+1)得到.取 𝑥_(1)=𝑥_(2)=𝑥_(4)=⋯=𝑥_(2^(𝑘))=1/(𝑘+1),其余 𝑥_(𝑖)=0,等号成立.

“

注 上述分组给出了整除图的一个显式 (𝑘+1)-着色;初等证明实际上只使用“同组内无边”这一性质,并不需要知道整张图的全部边.一般地,若一个图的顶点可分成 𝑟个独立集,且各组顶点权之和为 𝐴_(1),…,𝐴_(𝑟),则

∑_(𝑢𝑣∈𝐸)𝑥_(𝑢)𝑥_(𝑣)≤∑_(𝑠<𝑡)𝐴_(𝑠)𝐴_(𝑡)≤(𝑟−1)/(2𝑟)(∑_(𝑣)𝑥_(𝑣))^(2).

这一估计是完全 𝑟部图情形的直接加权形式.


三、局部选边:将 max项化为边权和

“

方法要点 局部 max或 min只在有限候选中取值时,可在变量取定后保留使极值取得的候选边.所得图虽依赖于变量,但选边规则固定,因此可统一控制所有可能选边图中的三角形、完全图等局部构型.

例 3.1|数之谜原创代数 542

给定整数 𝑛≥2.非负实数 𝑥_(1),𝑥_(2),…,𝑥_(𝑛)满足 ∑_(𝑛)^(𝑖=1)𝑥_(𝑖)=1.下标按模 𝑛理解.求

𝑆_(𝑛)=∑_(𝑖=1)^(𝑛)𝑥_(𝑖)max{𝑥_(𝑖+1),𝑥_(𝑖+2)}

的最大值.

分析 对每个 𝑖,在 𝑖+1,𝑖+2中选取使最大值取得的顶点,并将 𝑖与该顶点连边.𝑛≥5时所得无向边互不重复,因而 𝑆_(𝑛)可写成一张简单图的边权和.接下来不必精确描述所有可能的选边图,只需控制它们的最大团即可.𝑛=2,3,4时会出现重边或短周期现象,必须分别计算.

证明 当 𝑛=2时,𝑆_(2)=(𝑥_(1)+𝑥_(2))max{𝑥_(1),𝑥_(2)}≤1,故最大值为 1.

当 𝑛=3时,不妨记三个数为 𝑋≥𝑌≥𝑍.则

𝑆_(3)=2𝑋𝑌+𝑋𝑍=𝑋(1−𝑋+𝑌).

若 𝑋≤(1)/(2),由 𝑌≤𝑋得 𝑆_(3)≤𝑋≤(1)/(2);若 𝑋≥(1)/(2),由 𝑌≤1−𝑋得 𝑆_(3)≤2𝑋(1−𝑋)≤(1)/(2).所以最大值为 1/2.

当 𝑛=4时,令 𝐴=𝑥_(1)+𝑥_(3), 𝐵=𝑥_(2)+𝑥_(4),于是 𝐴+𝐵=1.由 max{𝑢,𝑣}≤𝑢+𝑣,

𝑆_(4) ≤∑_(𝑖=1)^(4)𝑥_(𝑖)(𝑥_(𝑖+1)+𝑥_(𝑖+2)); =𝐴𝐵+2(𝑥_(1)𝑥_(3)+𝑥_(2)𝑥_(4)); ≤𝐴𝐵+(𝐴^(2)+𝐵^(2))/(2); =(1)/(2).;

取 (𝑥_(1),𝑥_(2),𝑥_(3),𝑥_(4))=((1)/(2),0,(1)/(2),0)时等号成立.

以下设 𝑛≥5.对每个 𝑖选 𝑗_(𝑖)∈{𝑖+1,𝑖+2},使 𝑥_(𝑗_(𝑖))=max{𝑥_(𝑖+1),𝑥_(𝑖+2)},并连无向边 {𝑖,𝑗_(𝑖)}.若同一条无向边被两个起点重复选中,则只能由这条边的两个端点相向选取,因此两个前进步长之和等于 𝑛;而步长都属于 {1,2},故 𝑛≥5时不可能.因此所得图 𝐺有 𝑛条不同的边,且 𝑆_(𝑛)=∑_(𝑢𝑣∈𝐸(𝐺))𝑥_(𝑢)𝑥_(𝑣).当 𝑛=5时,𝐺只有 5条边,而 𝐾_(4)有 6条边,故 𝐺不含 𝐾_(4),从而 𝜔(𝐺)≤3.记 𝐶2𝑛为圈 𝐶_(𝑛)的平方,即循环距离不超过 2的两个不同顶点之间均连边.当 𝑛=6时,𝐺⊂𝐶26,而 𝐶26的补图恰由三条互不相交的对径边组成,所以任一团至多从每对对径点中取一个顶点,故 𝜔(𝐶26)=3.于是由例 1.1,𝑆_(𝑛)≤(1)/(3) (𝑛=5,6).取 𝑥_(1)=𝑥_(3)=𝑥_(5)=1/3,其余为 0,直接代回可得等号.

若 𝑛≥7,考察三角形的任一顶点 𝑖,另外两点均须位于 𝑖±1,𝑖±2中;两点之间的循环距离还须不超过 2.因此 𝐶2𝑛中的三角形只能由三个循环连续顶点 {𝑖,𝑖+1,𝑖+2}构成.对这样的三元组,边 {𝑖,𝑖+1}与 {𝑖,𝑖+2}都只能由起点 𝑖选出:从另一端反向选取所需的步长分别为 𝑛−1与 𝑛−2,均不属于 {1,2}.因此若三边同时出现,顶点 𝑖就必须同时选择两个候选点,与“每个起点只选一条边”的构造矛盾.因此 𝐺无三角形.由 1/4引理,𝑆_(𝑛)≤(1)/(4).取 𝑥_(1)=𝑥_(2)=1/2,其余为 0,等号成立.综上,

max𝑆_(𝑛)=⎧⎩⎨⎪⎪1, 𝑛=2,;(1)/(2), 𝑛=3,4,;(1)/(3), 𝑛=5,6,;(1)/(4), 𝑛≥7.;

四、谱方法:邻接矩阵、图拉普拉斯与 DFT

“

方法要点 平方和约束下的二次型适合写成矩阵形式:相邻乘积对应邻接矩阵,差分平方对应图拉普拉斯;若系数具有循环平移不变性,则可用 DFT 分离各频率.

设图 𝐺的邻接矩阵为 𝐴,顶点上的实数写成列向量 𝑥.由于每条无向边在 𝑥^(𝖳)𝐴𝑥中出现两次,

2∑_(𝑢𝑣∈𝐸(𝐺))𝑥_(𝑢)𝑥_(𝑣)=𝑥^(𝖳)𝐴𝑥.

记 𝐴的最大特征值为 𝜆_(max)(𝐴).由实对称矩阵最大特征值的 Rayleigh 商表征,

∑_(𝑢𝑣∈𝐸(𝐺))𝑥_(𝑢)𝑥_(𝑣)≤(𝜆_(max)(𝐴))/(2)∑_(𝑣∈𝑉(𝐺))𝑥2𝑣.

以下主要使用最大特征值 𝜆_(max)(𝐴).若以邻接矩阵的谱半径 𝜌(𝐴)表述,则由 Perron--Frobenius 定理有 𝜌(𝐴)=𝜆_(max)(𝐴);后文不再单独使用谱半径.若只需较弱而简便的估计,由 2𝑥_(𝑢)𝑥_(𝑣)≤𝑥2𝑢+𝑥2𝑣还可得

∑_(𝑢𝑣∈𝐸(𝐺))𝑥_(𝑢)𝑥_(𝑣)≤(1)/(2)∑_(𝑣)𝑑(𝑣)𝑥2𝑣≤(Δ(𝐺))/(2)∑_(𝑣)𝑥2𝑣,

其中 𝑑(𝑣)和 Δ(𝐺)分别表示顶点 𝑣的度数和图的最大度数.

例 4.1|路图上的精确极值

设 𝑛≥2,实数 𝑥_(1),𝑥_(2),…,𝑥_(𝑛)满足

𝑥21+𝑥22+⋯+𝑥2𝑛=1.

求

𝑆=𝑥_(1)𝑥_(2)+𝑥_(2)𝑥_(3)+⋯+𝑥_(𝑛−1)𝑥_(𝑛)

的最大值.

分析 目标式是路图 𝑃_(𝑛)的边权和,而约束是平方和.这正是 Rayleigh 商最适合的情形:只需求出 𝑃_(𝑛)邻接矩阵的最大特征值,并找出相应特征向量,即可同时得到最优常数和等号情形.

证明 设 𝐴为 𝑃_(𝑛)的邻接矩阵.若 𝑣=(𝑣_(1),…,𝑣_(𝑛))^(𝖳)是特征向量,并约定 𝑣_(0)=𝑣_(𝑛+1)=0,则

𝑣_(𝑖−1)+𝑣_(𝑖+1)=𝜆𝑣_(𝑖)(1≤𝑖≤𝑛).

取 𝑣_(𝑖)=sin⁡(𝑖𝜃),由恒等式

sin⁡((𝑖−1)𝜃)+sin⁡((𝑖+1)𝜃)=2cos⁡𝜃sin⁡(𝑖𝜃)

得到 𝜆=2cos⁡𝜃.边界条件 𝑣_(𝑛+1)=0给出

𝜃=(𝑘𝜋)/(𝑛+1)(𝑘=1,2,…,𝑛).

由此得到 𝑛个互异的特征值

2cos⁡(𝑘𝜋)/(𝑛+1)(𝑘=1,2,…,𝑛),

故它们恰为 𝐴的全部特征值.因此

𝜆_(max)(𝐴)=2cos⁡(𝜋)/(𝑛+1).

因此

𝑆≤cos⁡(𝜋)/(𝑛+1).

取

𝑥_(𝑖)=√((2)/(𝑛+1))sin⁡(𝑖𝜋)/(𝑛+1)(1≤𝑖≤𝑛),

利用

∑_(𝑖=1)^(𝑛)sin^(2)⁡(𝑖𝜋)/(𝑛+1)=(𝑛+1)/(2),

可知这组数满足平方和为 1,且对应最大特征值.故

max𝑆=cos⁡(𝜋)/(𝑛+1).

图拉普拉斯与差分平方

设 𝐷为度数对角矩阵,图拉普拉斯矩阵定义为 𝐿=𝐷−𝐴.直接展开有

𝑥^(𝖳)𝐿𝑥=∑_(𝑢𝑣∈𝐸(𝐺))(𝑥_(𝑢)−𝑥_(𝑣))^(2).

若 𝐺连通,将 𝐿的特征值按从小到大排列为

0=𝜆_(1)<𝜆_(2)≤⋯≤𝜆_(|𝑉|).

𝜆_(1)=0对应常向量 𝟏.当 ∑_(𝑣)𝑥_(𝑣)=0时,𝑥⊥𝟏,故 Rayleigh 商给出

∑_(𝑢𝑣∈𝐸(𝐺))(𝑥_(𝑢)−𝑥_(𝑣))^(2)≥𝜆_(2)∑_(𝑣)𝑥2𝑣.

𝜆_(2)常称为图的代数连通度.邻接矩阵用于估计相邻顶点乘积和,图拉普拉斯用于估计相邻顶点差分平方.

循环图与离散傅里叶变换

圈图 𝐶_(𝑛)的邻接矩阵和图拉普拉斯矩阵都是循环矩阵,可用离散傅里叶变换(DFT)同时对角化.令 i^(2)=−1,并设

𝜁=𝑒^(2𝜋i/𝑛),𝑥^(ˆ)_(𝑘)=(1)/(√(𝑛))∑_(𝑗=1)^(𝑛)𝑥_(𝑗)𝜁^(−𝑗𝑘)(𝑘=0,1,…,𝑛−1).

当 𝑥_(𝑗)均为实数时,𝑥^(ˆ)_(𝑛−𝑘)=𝑥^(ˆ)_(𝑘)^(–)(下标按模 𝑛理解).帕塞瓦尔(Parseval)等式为

∑_(𝑗=1)^(𝑛)𝑥2𝑗=∑_(𝑘=0)^(𝑛−1)|𝑥^(ˆ)_(𝑘)|^(2).

循环移位 𝑥_(𝑗)↦𝑥_(𝑗+1)在第 𝑘个频率上对应乘以 𝜁^(𝑘).因此,约定 𝑥_(𝑛+1)=𝑥_(1),有

∑_(𝑗=1)^(𝑛)(𝑥_(𝑗)−𝑥_(𝑗+1))^(2)=∑_(𝑘=0)^(𝑛−1)|1−𝜁^(𝑘)|^(2)|𝑥^(ˆ)_(𝑘)|^(2)=4∑_(𝑘=0)^(𝑛−1)sin^(2)⁡(𝑘𝜋)/(𝑛)|𝑥^(ˆ)_(𝑘)|^(2).

由 Parseval 等式计算循环移位的内积并取实部,还得到

2∑_(𝑗=1)^(𝑛)𝑥_(𝑗)𝑥_(𝑗+1)=2∑_(𝑘=0)^(𝑛−1)cos⁡(2𝑘𝜋)/(𝑛)|𝑥^(ˆ)_(𝑘)|^(2).

因此 𝐴(𝐶_(𝑛))的特征值为

2cos⁡(2𝑘𝜋)/(𝑛)(𝑘=0,1,…,𝑛−1),

而 𝐿(𝐶_(𝑛))的特征值为

4sin^(2)⁡(𝑘𝜋)/(𝑛)(𝑘=0,1,…,𝑛−1).

例 4.2|Fan–Taussky–Todd 离散 Wirtinger 不等式(周期型)[3]

设 𝑛≥3,实数 𝑥_(1),…,𝑥_(𝑛)满足 𝑥_(1)+⋯+𝑥_(𝑛)=0,并约定 𝑥_(𝑛+1)=𝑥_(1).证明

∑_(𝑖=1)^(𝑛)(𝑥_(𝑖)−𝑥_(𝑖+1))^(2)≥4sin^(2)⁡(𝜋)/(𝑛)∑_(𝑖=1)^(𝑛)𝑥2𝑖,

并确定等号情形.

分析 零和条件等价于 𝑥^(ˆ)_(0)=0.在非零频率 𝑘=1,…,𝑛−1中,4sin^(2)⁡(𝑘𝜋/𝑛)的最小值为 4sin^(2)⁡(𝜋/𝑛).

证明 由上面的 DFT 恒等式及 𝑥^(ˆ)_(0)=0,

∑_(𝑖=1)^(𝑛)(𝑥_(𝑖)−𝑥_(𝑖+1))^(2) =4∑_(𝑘=1)^(𝑛−1)sin^(2)⁡(𝑘𝜋)/(𝑛)|𝑥^(ˆ)_(𝑘)|^(2); ≥4sin^(2)⁡(𝜋)/(𝑛)∑_(𝑘=1)^(𝑛−1)|𝑥^(ˆ)_(𝑘)|^(2); =4sin^(2)⁡(𝜋)/(𝑛)∑_(𝑖=1)^(𝑛)𝑥2𝑖.;

等号成立当且仅当傅里叶系数只可能出现在 𝑘=1,𝑛−1两个频率上.对实数序列,这等价于

𝑥_(𝑗)=𝐴cos⁡(2𝜋𝑗)/(𝑛)+𝐵sin⁡(2𝜋𝑗)/(𝑛)(1≤𝑗≤𝑛),

其中 𝐴,𝐵∈ℝ;𝐴=𝐵=0时对应零序列.

“

注 例 4.2 是 Fan–Taussky–Todd 离散 Wirtinger 不等式的周期型.Ky Fan、Olga Taussky 与 John Todd 在 1955 年论文 Discrete analogs of inequalities of Wirtinger中证明了周期型以及若干边界型离散不等式.其两端为零的形式为

∑_(𝑖=0)^(𝑛)(𝑥_(𝑖+1)−𝑥_(𝑖))^(2)≥4sin^(2)⁡(𝜋)/(2(𝑛+1))∑_(𝑖=1)^(𝑛)𝑥2𝑖,𝑥_(0)=𝑥_(𝑛+1)=0,

见习题 8.4.周期型的自然基底是 DFT 的复指数函数;两端为零时相应的基底为 sin⁡(𝑖𝑘𝜋/(𝑛+1)),即离散正弦变换(DST)的基底,也可由长度 2(𝑛+1)的奇延拓与 DFT 得到.因而例 4.1、例 4.2 与习题 8.4 都可理解为相应对称矩阵在正弦基或傅里叶基下的对角化.


五、有向路径:由局部选择构造乘积链

“

方法要点 比值具有明确先后次序,且每项只有有限个向前候选时,可用有向边记录选择.候选步长的上界控制路径长度,而路径上比值相乘时中间变量相消,便于随后应用均值不等式.

例 5.1|2019 IZhO Day 1, Problem 2[5]

求最大的实数 𝐶,使任意两两不同的正实数 𝑎_(1),𝑎_(2),…,𝑎_(2019)都满足

(𝑎_(1))/(|𝑎_(2)−𝑎_(3)|)+(𝑎_(2))/(|𝑎_(3)−𝑎_(4)|)+⋯+(𝑎_(2019))/(|𝑎_(1)−𝑎_(2)|)>𝐶.

分析 因为题设中的各数两两不同,所以对正数 𝑎,𝑏,𝑐且 𝑏≠𝑐有

|𝑏−𝑐|<max{𝑏,𝑐},(𝑎)/(|𝑏−𝑐|)>(𝑎)/(max{𝑏,𝑐})=min{(𝑎)/(𝑏),(𝑎)/(𝑐)}.

因此每一项都可以严格压低为两个比值中的较小者.选择分母较大的候选项,并将这一选择记为一条向前的有向边,就得到步长为 1或 2的路径.将路径末端与起点连接后,所取比值的乘积为 1;再用 AM–GM 不等式把“路径项数”转成“比值和的下界”.

证明 设原式为 𝑆.循环平移下标后可设 𝑎_(1)=min_(1≤𝑖≤2019)𝑎_(𝑖).舍去第一项,对 𝑖=2,…,2017用上述基本估计;又因 𝑎_(1)最小且各数两两不同,有 |𝑎_(2019)−𝑎_(1)|<𝑎_(2019)、|𝑎_(1)−𝑎_(2)|<𝑎_(2).

于是

𝑆> min{(𝑎_(2))/(𝑎_(3)),(𝑎_(2))/(𝑎_(4))}+⋯; +min{(𝑎_(2017))/(𝑎_(2018)),(𝑎_(2017))/(𝑎_(2019))}; +(𝑎_(2018))/(𝑎_(2019))+(𝑎_(2019))/(𝑎_(2))=:𝑇.;

令 𝑖_(0)=2.当 𝑖_(ℓ)≤2017时,在 𝑖_(ℓ)+1与 𝑖_(ℓ)+2中选择使 𝑎_(𝑗)较大的下标 𝑗,并令 𝑖_(ℓ+1)=𝑗.设经过 𝑠次跳跃后第一次到达 {2018,2019},即

𝑖_(0)=2<𝑖_(1)<⋯<𝑖_(𝑠),𝑖_(𝑠)∈{2018,2019}.

每次跳跃长度为 1或 2.舍去 𝑇中未被这条路径使用的正项,可得

𝑇≥ (𝑎_(2))/(𝑎_(𝑖_(1)))+(𝑎_(𝑖_(1)))/(𝑎_(𝑖_(2)))+⋯; +(𝑎_(𝑖_(𝑠−1)))/(𝑎_(𝑖_(𝑠)))+(𝑎_(2018))/(𝑎_(2019))+(𝑎_(2019))/(𝑎_(2))=:𝐴.;

若 𝑖_(𝑠)=2018,则总前进距离为 2016,而每次至多前进 2,故

𝑠≥⌈(2016)/(2)⌉=1008.

此时 𝐴中沿路径的 𝑠个比值,再加上 𝑎_(2018)/𝑎_(2019)与 𝑎_(2019)/𝑎_(2),共 𝑠+2项,乘积恰为 1.由 AM–GM 不等式,

𝐴≥𝑠+2≥1010.

若 𝑖_(𝑠)=2019,则总前进距离为 2017,故

𝑠≥⌈(2017)/(2)⌉=1009.

舍去正项 𝑎_(2018)/𝑎_(2019)后,沿路径的 𝑠个比值与 𝑎_(2019)/𝑎_(2)共 𝑠+1项,乘积为 1,因而

𝐴≥𝑠+1≥1010.

因此 𝑆>1010.

下面说明常数不能增大.令 0<𝜀<1/1009,取

𝑎_(2𝑟−1)=1+𝑟𝜀(1≤𝑟≤1009),𝑎_(2𝑟)=𝑟𝜀(1≤𝑟≤1008),

并另取 𝑎_(2018)=𝜀^(2), 𝑎_(2019)=1.这些数两两不同.逐项代入原式可得

𝑆= ∑_(𝑟=1)^(1008)(1+𝑟𝜀)/(1+𝜀)+∑_(𝑟=1)^(1007)𝑟𝜀; +(1008𝜀)/(1+1009𝜀−𝜀^(2)); +(1+1009𝜀)/(1−𝜀^(2))+𝜀+1; ⟶1010.;

故原式在题设条件下的下确界为 1010,从而题目所求最大常数为 𝐶=1010.

“

注 当局部估计产生有限个候选比值时,可将所选候选项记录为有向边.若每一步的前进距离有统一上界,则可由总距离估计路径长度;在比值型问题中,再将路径首尾连接,常可得到乘积为 1的比值链,并应用 AM–GM 不等式.


六、哈密顿平均:由排列条件估计全体边权

“

方法要点 任意排列中的相邻乘积和可视为完全图上一条哈密顿路的权和;循环排列对应哈密顿圈.对全体此类结构取平均,可由完全图的对称性计算每条边的出现频率.

例 6.1|2017 IMO Shortlist A5[4]

设整数 𝑛≥3.实数 𝑥_(1),𝑥_(2),…,𝑥_(𝑛)满足:对它们的任意排列 𝑦_(1),𝑦_(2),…,𝑦_(𝑛),都有

𝑦_(1)𝑦_(2)+𝑦_(2)𝑦_(3)+⋯+𝑦_(𝑛−1)𝑦_(𝑛)≥−1.

求最大的实数 𝐾=𝐾(𝑛),使恒有

∑_(1≤𝑖<𝑗≤𝑛)𝑥_(𝑖)𝑥_(𝑗)≥𝐾.

分析 在完全图 𝐾_(𝑛)的边 𝑖𝑗上赋权 𝑥_(𝑖)𝑥_(𝑗).每个排列对应一条哈密顿路.第一步是把“路的下界”转成“圈的下界”:若圈上有一条非负边,删去它即可.只有当所有数非零且正、负数目相等时,才可能出现每条圈边都为负的正负交替圈,因此这一情形需要单独处理.第二步再对适当的哈密顿圈取平均,利用完全图的对称性计算每条边的出现概率.

证明 记 𝑄=∑_(1≤𝑖<𝑗≤𝑛)𝑥_(𝑖)𝑥_(𝑗).先设存在零项,或正项与负项的个数不同.任一哈密顿圈上必有一条权重非负的边:若圈上没有零点且所有相邻边都为负,则顶点符号必须严格正负交替,从而正、负点数相等.删去一条非负边后得到哈密顿路,故该圈总权至少为 −1.

对所有无向哈密顿圈等概率平均.一个圈含 𝑛条边,而 𝐾_(𝑛)有 𝑛(𝑛−1)/2条边;由对称性,每条固定边在随机哈密顿圈中出现的概率为 2/(𝑛−1).因此平均圈权为 (2)/(𝑛−1)𝑄≥−1,从而 𝑄≥−(𝑛−1)/(2).

余下情形中,𝑛=2𝑚,且全部 𝑥_(𝑖)非零,其中恰有 𝑚个正数与 𝑚个负数.写成 𝑝_(1),…,𝑝_(𝑚)>0, −𝑞_(1),…,−𝑞_(𝑚)<0,并记 𝑃=∑𝑝_(𝑖), 𝑅=∑𝑞_(𝑗).考虑正负交替的哈密顿圈.设某一圈全部负边权绝对值之和为 𝑊,则圈权为 −𝑊.删去绝对值为 𝑡的一条边后,所得哈密顿路权为 −𝑊+𝑡≥−1,故 𝑊−𝑡≤1.对圈上 𝑛条边求和,得到

(𝑛−1)𝑊≤𝑛,𝑊≤(𝑛)/(𝑛−1).

在所有正负交替的无向哈密顿圈中等概率平均.固定一对正、负顶点;由对称性,一个给定正顶点的两个邻点在 𝑚个负顶点中均匀分布,因此这对顶点相邻的概率为 2/𝑚=4/𝑛.于是

(4)/(𝑛)𝑃𝑅≤(𝑛)/(𝑛−1),𝑃𝑅≤(𝑛^(2))/(4(𝑛−1)).

异号边权之和为 −𝑃𝑅,同号边权均为正,于是 𝑄>−𝑃𝑅≥−(𝑛^(2))/(4(𝑛−1))≥−(𝑛−1)/(2)(此时 𝑛≥4).故所有情形都有 𝑄≥−(𝑛−1)/2.

最后说明常数不能增大.取 𝑥_(2)=𝑥_(3)=⋯=𝑥_(𝑛)=𝑡>0, 𝑥_(1)=−(1)/(2𝑡).若负数位于哈密顿路内部,路权为 −1+(𝑛−3)𝑡^(2)≥−1;若负数位于端点,路权为 −1/2+(𝑛−2)𝑡^(2)≥−1.因此这组数满足题设,而

𝑄=−(𝑛−1)/(2)+((𝑛−1)(𝑛−2))/(2)𝑡^(2)⟶−(𝑛−1)/(2)(𝑡→0^(+)).

所以

𝐾(𝑛)=−(𝑛−1)/(2).

七、方法归纳

把一道不等式作图论化处理时,可按以下四个步骤展开:

  1. 确定顶点.通常一个变量、一个指标或一个离散对象对应一个顶点.

  2. 确定边的含义.边可以表示“这一乘积出现”“两个对象可比较”“这一候选被选中”或“排列中彼此相邻”.

  3. 把原式改写成图上的量.常见形式包括边权和、矩阵二次型、路径上的比值和以及哈密顿路/圈的权和.

  4. 估计所需的图参数.通常无须完全刻画图的结构;若能够控制团数、独立集分组、路径长度、最大特征值或边的平均出现概率,即可满足证明需要.

常见结构与工具的对应关系

  1. 非负权重、∑𝑥_(𝑣)固定、目标是 ∑_(𝑢𝑣∈𝐸)𝑥_(𝑢)𝑥_(𝑣):用 Motzkin–Straus;若已有独立集分组,也可直接合并组内权重.

  2. 各顶点权相等、问题退化为计边:用 Turán 型极值.

  3. 整除、包含、可比等固定关系:先找团所对应的链,或找独立集所对应的分层.

  4. 局部 max/min:先按变量大小“选边”,再研究所有可能选边图都满足的禁构型.

  5. ∑𝑥2𝑣固定的二次式:用邻接矩阵或图拉普拉斯;循环对称时优先考虑 DFT,路径边界则常出现正弦基.

  6. 局部比值有有限个向前候选:构造有向路径,由步长控制项数,再利用乘积消去中间变量.

  7. 任意排列的相邻项:转成哈密顿路或哈密顿圈,并对全体结构取平均.

最后还需注意三点.第一,图上的量与原式是“相等”还是仅给出“上/下界”;第二,每次使用柯西、AM–GM、Rayleigh 商或平均法时,等号条件能否同时满足;第三,若证明中出现严格不等式,最终常数通常只是上确界或下确界,需要另给逼近构造.例 5.1 的 1010是下确界而不是实际最小值;例 4.1、4.2 则由相应特征向量直接给出等号.


八、精选习题

下面六题分别对应团数、关系图、局部选边、谱方法、有向路径与哈密顿平均.作答时可先判断原式中的“成对关系”来自固定关系、局部选择、二次型还是排列相邻,再据此确定顶点、边及相应的归一化条件,最后估计所需图参数或矩阵特征值.

习题 8.1|带状邻接图

设整数 𝑛,𝑟满足 1≤𝑟<𝑛,非负实数 𝑥_(1),…,𝑥_(𝑛)满足 𝑥_(1)+⋯+𝑥_(𝑛)=1.求

∑_(1≤𝑖<𝑗≤𝑛;𝑗−𝑖≤𝑟;)𝑥_(𝑖)𝑥_(𝑗)

的最大值.

习题 8.2|子集包含关系

设 𝑚≥1.对 [𝑚]={1,2,…,𝑚}的每个子集 𝐴,给定非负实数 𝑥_(𝐴),且 ∑_(𝐴⊆[𝑚])𝑥_(𝐴)=1.求

∑_(𝐴⊊𝐵⊆[𝑚])𝑥_(𝐴)𝑥_(𝐵)

的最大值.

习题 8.3|局部选边:步长 1或 3

设整数 𝑛≥10.非负实数 𝑥_(1),…,𝑥_(𝑛)满足 ∑_(𝑛)^(𝑖=1)𝑥_(𝑖)=1,下标按模 𝑛理解.求

𝑆=∑_(𝑖=1)^(𝑛)𝑥_(𝑖)max{𝑥_(𝑖+1),𝑥_(𝑖+3)}

的最大值.

习题 8.4|Fan–Taussky–Todd 不等式(两端为零)

设 𝑛≥1,实数 𝑥_(0),𝑥_(1),…,𝑥_(𝑛),𝑥_(𝑛+1)满足 𝑥_(0)=𝑥_(𝑛+1)=0且∑_(𝑛)^(𝑖=1)𝑥2𝑖=1.求

𝐸=∑_(𝑖=0)^(𝑛)(𝑥_(𝑖+1)−𝑥_(𝑖))^(2)

的最小值.

习题 8.5|有向路径的 𝑘候选推广

设整数 𝑘≥2, 𝑚≥1,记 𝑁=𝑘𝑚+1,并设 𝑎_(1),𝑎_(2),…,𝑎_(𝑁)为正实数.证明

∑_(𝑖=1)^(𝑁−𝑘)min_(1≤𝑟≤𝑘)(𝑎_(𝑖))/(𝑎_(𝑖+𝑟)); +∑_(𝑖=𝑁−𝑘+1)^(𝑁−1)(𝑎_(𝑖))/(𝑎_(𝑖+1))+(𝑎_(𝑁))/(𝑎_(1)); ≥𝑚+1.;

习题 8.6|哈密顿圈平均

设 𝑛≥3,𝑥_(1),…,𝑥_(𝑛)为实数,并且对它们的任意排列 𝑦_(1),…,𝑦_(𝑛)都有

𝑦_(1)𝑦_(2)+𝑦_(2)𝑦_(3)+⋯+𝑦_(𝑛−1)𝑦_(𝑛)+𝑦_(𝑛)𝑦_(1)≤1.

求 ∑_(1≤𝑖<𝑗≤𝑛)𝑥_(𝑖)𝑥_(𝑗)可能取得的最大值.


九、习题提示与答案

8.1

在顶点集 {1,2,…,𝑛}上,若 1≤|𝑖−𝑗|≤𝑟就连边.任一团的最大、最小下标之差不超过 𝑟,故团数至多为 𝑟+1;而任意 𝑟+1个连续顶点构成 𝐾_(𝑟+1).由例 1.1,最大值为

(𝑟)/(2(𝑟+1)).

在任意 𝑟+1个连续顶点上各取权重 1/(𝑟+1),其余权重取 0,等号成立.

8.2

将 2^([𝑚])的每个子集作为顶点,两个不同子集可比较时连边.团对应严格包含链.集合大小每次至少增加 1,故链长至多为 𝑚+1;而 ∅⊂{1}⊂{1,2}⊂⋯⊂[𝑚]达到该长度.故最大值为

(𝑚)/(2(𝑚+1)).

在上述任意一条长度为 𝑚+1的最大链上各取权重 1/(𝑚+1),其余权重取 0,即可达到等号.

8.3

对每个 𝑖,在 𝑖+1与 𝑖+3中选出权重较大的顶点并连边.𝑛≥10时同一无向边不可能由两个端点相向重复选中.若所得图有三角形,沿三边取带符号步长 𝜀_(1),𝜀_(2),𝜀_(3)∈{±1,±3},则三者之和是 𝑛的倍数;但其绝对值至多为 9<𝑛,故只能为 0;而三个奇数之和为奇数,不可能等于 0,矛盾.因此图无三角形,例 1.1 给出 𝑆≤1/4.取 𝑥_(1)=𝑥_(2)=1/2、其余为 0时等号成立,所以

max𝑆=(1)/(4).

8.4

左端是两端为零的路径差分平方和.对应的 𝑛×𝑛对称矩阵主对角线为 2,两条副对角线为 −1.取

𝑣_(𝑖)=sin⁡(𝑖𝑘𝜋)/(𝑛+1)(1≤𝑖≤𝑛),

可得特征值

𝜆_(𝑘)=2−2cos⁡(𝑘𝜋)/(𝑛+1)=4sin^(2)⁡(𝑘𝜋)/(2(𝑛+1))(1≤𝑘≤𝑛).

因此最小值为

4sin^(2)⁡(𝜋)/(2(𝑛+1)).

等号在 𝑥_(𝑖)=𝑐sin⁡(𝑖𝜋)/(𝑛+1)时取得,其中 𝑐由平方和为 1确定.

8.5

从 𝑖_(0)=1出发;当 𝑖_(ℓ)≤𝑁−𝑘时,在 𝑖_(ℓ)+1,…,𝑖_(ℓ)+𝑘中选取使 𝑎_(𝑗)最大的下标 𝑗,并令 𝑖_(ℓ+1)=𝑗.于是

min_(1≤𝑟≤𝑘)(𝑎_(𝑖_(ℓ)))/(𝑎_(𝑖_(ℓ)+𝑟))=(𝑎_(𝑖_(ℓ)))/(𝑎_(𝑖_(ℓ+1))).

设经过 𝑠次跳跃后第一次进入 {𝑁−𝑘+1,…,𝑁}.由于每次至多前进 𝑘,而

(𝑁−𝑘+1)−1=𝑘(𝑚−1)+1,

故 𝑠𝑘≥𝑘(𝑚−1)+1,从而 𝑠≥𝑚.随后沿相邻比值依次连接至 𝑁,再加入 𝑎_(𝑁)/𝑎_(1).所得闭合比值链中各项都来自原式,且中间变量全部相消,乘积为 1;其中至少包含 𝑠个跳跃比值和末项 𝑎_(𝑁)/𝑎_(1),故项数至少为 𝑚+1.对这些正数应用均值不等式,其和至少为 𝑚+1;原式其余各项均为正,因此结论成立.

8.6

对所有无向哈密顿圈等概率平均.由对称性,每条边出现的次数相同;每个圈含 𝑛条边,而 𝐾_(𝑛)共有 𝑛(𝑛−1)/2条边,因此固定边出现在随机哈密顿圈中的概率为

(𝑛)/(𝑛(𝑛−1)/2)=(2)/(𝑛−1).

题设说明每个哈密顿圈的边权和至多为 1,故其平均值也至多为 1.于是

(2)/(𝑛−1)∑_(𝑖<𝑗)𝑥_(𝑖)𝑥_(𝑗) ≤1,;∑_(𝑖<𝑗)𝑥_(𝑖)𝑥_(𝑗) ≤(𝑛−1)/(2).;

取 𝑥_(1)=⋯=𝑥_(𝑛)=1/√(𝑛)时,每个哈密顿圈含 𝑛条权为 1/𝑛的边,圈权恰为 1,故上界可以达到.

说明 以下仅列出本文直接使用的经典定理原始文献及两道国际竞赛题的官方题源,供进一步查阅.

参考文献

[1] T. S. Motzkin and E. G. Straus, Maxima for Graphs and a New Proof of a Theorem of Turán, Canadian Journal of Mathematics17 (1965), 533–540.

[2] P. Turán, On the Theory of Graphs, Colloquium Mathematicum3 (1954), 19–30.

[3] Ky Fan, Olga Taussky, and John Todd, Discrete analogs of inequalities of Wirtinger, Monatshefte für Mathematik59 (1955), 73–90.

[4] International Mathematical Olympiad, Shortlisted Problems with Solutions, 58th IMO, Rio de Janeiro, 2017, Problem A5.

[5] International Zhautykov Olympiad, Mathematics Problems, 2019, Day 1, Problem 2.


留言

正常留言提交后可直接显示,无需人工审核;仍会经过防刷与垃圾内容过滤。请不要留下手机号等私人信息。解法、疑问都可以说;公式用 LaTeX:行内 $…$,整行 $$…$$。