精读笔记

Problem Setting

题目：Cesàro means of firmly nonexpansive iterates need not converge strongly（arXiv preprint / 2026）。

论文实际处理的是 Baillon nonlinear mean ergodic theorem 的强收敛边界：对 Hilbert 空间中的非扩张映射，Cesàro means 已知弱收敛到某个固定点；问题是如果进一步假设映射 firmly nonexpansive，是否仍可能不强收敛。这个问题重要是因为 firmly nonexpansive 是 resolvent、projection、proximal map 的共同抽象，很多优化算法的平均化策略都默认“弱收敛 + averaging”能在实践或附加条件下接近强收敛。

真正困难点是同时满足三件事：轨道必须弱趋于唯一固定点；轨道范数不能趋于 0；更强的是 Cesàro means 也不能被平均操作压到 0。单纯构造弱不强收敛序列并不难，难的是让这个序列恰好是某个 firmly nonexpansive map 的 orbit。firm nonexpansiveness 对相邻点施加了强几何约束，例如需要类似 ||x_{n+1}||^2 ≤ <x_n,x_{n+1}> 的关系，这会迫使范数单调下降；因此反例必须在“范数下降”与“下降总量有限”之间做精确平衡。

以前方法卡在两个地方：Genel–Lindenstrauss 给出了 firmly nonexpansive 迭代弱不强收敛，但没有解决 Cesàro means 的强收敛失败；von Neumann 线性均值遍历定理又给出强收敛，容易让人误以为 nonlinear firm nonexpansiveness 也足够。本文的关键矛盾就是：线性结构中的平均化会产生强收敛，但非线性 firmly nonexpansive 只保留局部单调性，不保留足够的谱/正交分解结构。

Motivation

作者的核心观察是：若要让 firmly nonexpansive orbit 不强收敛，必须让点 x_n 在 Hilbert 空间中“逃向无穷多方向”，从而弱收敛到 0；但相邻点又必须高度相关，使 firm nonexpansive 所需的不等式成立，并且范数下降有限。换句话说，需要一条无限维中的弱消失曲线，沿曲线缓慢移动。

这解释了为什么 Gaussian kernel 自然出现。作者不是从已有反例的几何构型出发，而是反向设计一个 orbit：先设 x_n=rho_n u(t_n)，其中 u(t) 是单位向量曲线，t_n→∞ 保证弱趋于 0；再让相邻内积控制 rho_n 的衰减。若 <u(s),u(t)>=Phi(|s-t|)，则 rho_{n+1}=rho_n Phi(d_n)。为了 t_n 发散而 rho_n 不塌缩，需要 sum d_n=∞ 但 sum -log Phi(d_n)<∞。取 Gaussian Phi(d)=exp(-d^2) 后，问题变成选择 d_n∈l2\l1。

关键缺口不是“还没有一个例子”，而是还没有一个可以直接解释 averaging 为什么失败的机制。本文的构造把失败归因到 mesh 的局部聚团：在 harmonic mesh 下，后半段 index 对应的 t 区间宽度仍然有界，因此一大批 x_k 仍高度正相关，Cesàro 平均不会抵消。

Core Idea

核心思想是把非线性固定点迭代的反例转化为 RKHS/feature-curve 几何设计。选择一条单位向量曲线 u(t)，其内积是 Gaussian kernel exp(-(s-t)^2)。这条曲线有两个同时需要的性质：t→∞ 时 u(t) 弱趋于 0；但当 |s-t| 小时，u(s),u(t) 高度对齐。然后用一个慢速发散 mesh t_n 采样这条曲线，并用 rho_n 修正半径，使相邻点刚好满足 firm nonexpansiveness 的等号版本。

本质区别在于，作者不是从一个已有算子分析其 orbit，而是先构造一个“合法 orbit”，再通过 firmly nonexpansive extension theorem 反推出全空间算子。这个建模方式很关键：它把 operator construction 的难题拆成两个可控子问题——离散轨道上的单调性验证，以及全空间扩张。Gaussian kernel 给出的是一种 inductive bias：远处近似正交以产生弱消失，近处强相关以维持 firm nonexpansive 几何。

对 Cesàro means 的失败，核心不是范数本身有正下界，而是平均中的大块向量没有足够 cancellation。harmonic mesh 的后半段 t_n 只覆盖长度 O(delta) 的区间，因此至少一半项彼此内积有统一正下界；平均范数自然也有正下界。block variant 则说明通过调节 mesh 的局部密度，可以让平均有时被局部聚团支撑，有时被长区间扩散稀释，从而产生 liminf=0、limsup>0。

Method

1. Gaussian feature curve：解决的是如何在无限维空间中获得“弱趋于 0 但局部连续相关”的单位向量族。作者给出两种 realization：基于正交基的指数展开，以及 L2(R) 中平移 Gaussian。机制上只需要 kernel identity <u(s),u(t)>=exp(-(s-t)^2)。这不是装饰性构造，而是后续所有内积估计的坐标系。

2. l2\l1 mesh：解决的是“走到无穷远”和“总衰减有限”的冲突。d_n>0、sum d_n=∞ 让 t_n→∞；sum d_n^2<∞ 让 rho_n=exp(-sum_{k<n}d_k^2) 收敛到正数。附加条件 d_{n+1}≤d_n/(1+64d_n^2) 是技术性但重要的单调性控制，用于证明离散 shift map 的 firm nonexpansive 条件。

3. Orbit construction：设 x_n=rho_n u(t_n)。由 Gaussian kernel 得到 <x_{n+1},x_n>=||x_{n+1}||^2，即 firm nonexpansiveness 所需的相邻等号关系。更复杂的 Lemma 6.3 验证任意 m<n 的二点 firm nonexpansive monotonicity 条件，这是能扩张为全空间 firmly nonexpansive map 的关键。

4. Extension step：在离散集 {x_n}∪{0} 上定义 T_0 x_n=x_{n+1}, T_0 0=0。证明 T_0 firmly nonexpansive 后，用 refined Kirszbraun–Valentine theorem 扩张到 T:X→X，同时控制 range 落在闭凸包 C 中。这个步骤把“构造 orbit”升级为“存在真实算子”。

5. 固定点唯一性：仅有扩张不够，因为可能引入额外 fixed points。作者用 cone K 的 acute property 和一个单调内积引理排除非零 fixed point。这里的机制是：若 y 是 fixed point，则 firm nonexpansiveness 强迫 <x_n,y> 单调上升，但 x_n weakly→0 又迫使极限为 0，结合 y∈K 与 K acute 得 y=0。

6. Cesàro analysis：harmonic case 用后半段 index 的局部 t-直径有界证明平均范数正下界；block case 用交替的高密局部段和长宽块段制造 limsup 正、liminf 零。这部分不是 operator theory 的核心，而是 mesh 设计对 averaging 行为的精细控制。

Key Insight / Why It Works

最核心贡献是一个非常干净的几何解耦：firmly nonexpansive orbit 的构造可以先在离散轨道上完成，再通过扩张定理变成全空间算子。这样避免了直接写出复杂 nonlinear operator。真正有效的原因不是某个技术不等式，而是 Gaussian kernel + l2\l1 mesh 同时满足三种相互冲突的需求：弱消失、半径不塌缩、相邻 firm geometry。

Gaussian kernel 是关键选择。若 kernel 衰减太快，rho_n 会塌缩；若衰减太慢，远处向量相关性不够可控。exp(-d^2) 与 d_n∈l2\l1 正好匹配：路径长度无穷但能量有限。这和很多优化反例中的“infinite travel with finite energy”机制本质相同，只是这里通过 kernelized Hilbert geometry 表达得更透明。

Cesàro 失败的真正机制是局部相关块的质量占比不消失。对 harmonic mesh，前 n 项的后半部分落在有限 t-宽度内，因此这些向量几乎朝同一方向簇集；平均不是在大量近正交方向上均匀扩散，而是保留一个宏观 coherent block。这个 insight 很值得迁移：平均法是否强收敛，不只取决于单点弱收敛，还取决于窗口内相关结构的衰减速度。

哪些部分可能只是辅助？常数 64、1/8、各种指数不等式主要是保证 monotonicity 验证闭合的工程化数学常数，不是概念贡献。block construction 的具体 w_k=k+1、Q_k 递推也更像 proof engineering；真正思想是通过 mesh density 控制 Cesàro window 的 coherence/dispersion。

这不是 scaling，也不是数据覆盖，也不是 benchmark trick；它是一个纯理论反例。但如果用机器学习语言类比，它相当于构造了一个 feature trajectory，其中 representation 在弱拓扑下逃逸，但局部 minibatch/averaging window 保持高 alignment。所谓 averaging 失效，本质是 representation alignment 没有跨足够尺度 decorrelate。

Relation To Prior Work

最接近的是 Genel–Lindenstrauss 1975 的 firmly nonexpansive 弱不强收敛反例。本文与其本质差异在于：Genel–Lindenstrauss 主要回答 iterates 本身是否强收敛；本文进一步打击 Cesàro means，并且给出更可计算的 kernel/mesh 构造。它不是简单换一个反例，而是把反例机制从几何存在性改写成可调的相关结构设计。

与 von Neumann 线性均值遍历定理的关系是负向分界：线性非扩张算子的强 Cesàro 收敛依赖线性谱分解/闭包 range 分解；firm nonexpansiveness 虽强于 nonexpansiveness，但在非线性情形中不提供这种全局代数结构。本文说明 Baillon 定理不能按线性直觉强化。

与 Baillon nonlinear mean ergodic theorem 的关系更直接：Baillon 给弱收敛，本文证明强收敛一般不可能，即使 T firmly nonexpansive 且固定点唯一。这个补全了 theorem 的 sharpness。

与 Hundal、Güler、proximal point counterexamples 的关系需要谨慎。本文的 T 是 firmly nonexpansive，因此可视为某个 maximally monotone operator 的 resolvent；但文中明确指出对应 A 非 paramonotone，不是 subdifferential，T 不能是 proximal mapping。因此它没有覆盖 Güler 那类 proximal 算法负结果，反而说明 firm nonexpansive resolvent 类内部仍有结构层级：一般 resolvent 足以失败，但 proximal/subdifferential 情形是否同样 Cesàro 失败仍是更强问题。

看似新的部分中，Kirszbraun–Valentine extension 是已有工具，Gaussian RKHS realization 也是标准 kernel 技术；实质创新在于把这些工具组合成一个可调 orbit，并用 mesh 的局部/全局尺度分离控制 Cesàro 行为。

Dataset / Evaluation

这是纯数学论文，没有 dataset、实验或 benchmark。evaluation 的形式是构造性证明，claim 与证据高度一致：作者要证明“need not converge strongly”，只需给出一个合法反例；Theorem 7.1、Example 7.4 和 Section 8 分别覆盖 iterates 非强收敛、Cesàro means 统一远离 0、Cesàro norm 振荡三种强度。

如果把 evaluation 理解为理论覆盖范围，它验证的是一般 infinite-dimensional Hilbert space 中 firmly nonexpansive mappings 的最坏情形，而不是实际优化算法的典型行为。由于构造依赖无限维和非 proximal resolvent，它不能直接说明常见 finite-dimensional proximal algorithms 或带 compactness/regularity 条件的方法会失败。

文中没有 empirical claim，因此不存在 benchmark leakage、数据重叠或 hidden supervision 问题。需要注意的是，标签中出现“学习理论、模仿学习”与正文内容并不匹配；本文本质属于 nonlinear functional analysis / monotone operator theory / optimization convergence counterexample，而不是机器学习实验论文。

Limitation

第一，构造依赖无限维。有限维中弱收敛等价于强收敛，因此整个反例机制不存在。这意味着负结果针对的是无限维优化、函数空间、PDE/variational setting 或抽象 Hilbert 空间算法，而不是有限维数值优化的直接反例。

第二，T 的全空间形式通过 extension theorem 获得，虽然 orbit 显式，但 operator 本身不够 constructive。对于想要实现一个具体算法或 proximal step 的读者，这个反例的可操作性有限。它证明存在性与机制，但不是一个自然产生的 optimization model。

第三，T 不是 proximal mapping。文中通过非 paramonotonicity 说明对应 maximally monotone operator 不是 subdifferential，因此不能把结论直接迁移到 convex minimization 的 proximal point algorithm。若目标是理解实际 convex optimization 中 averaging 是否失败，这篇只给出下界式警告，而不是最终答案。

第四，关键技术条件 d_{n+1}≤d_n/(1+64d_n^2) 和若干常数设计有 proof-engineering 色彩。它们保证 firm nonexpansiveness 的全局二点不等式，但未必是 sharp。文中未充分说明这些 mesh 条件的必要性；更一般的 kernel/mesh characterization 仍缺失。

第五，block oscillation 证明 Cesàro norm 可 liminf=0、limsup>0，但这是通过刻意设计 mesh density 实现的。它说明 Baillon 结论不能强化为强收敛，但不说明这种振荡在自然算法中常见。这里的增益来源非常清楚：不是新算法能力，而是 adversarial geometry。

Takeaway

• 1. Baillon nonlinear mean ergodic theorem 在一般非线性 Hilbert setting 下基本已经 sharp：firm nonexpansiveness 和唯一 fixed point 仍不足以推出 Cesàro strong convergence。
• 2. 判断 averaged iterates 是否强收敛，不能只看单点轨道弱收敛和范数单调；必须看 Cesàro window 内的相关结构。
• 局部 coherent blocks 足以抵消 averaging 的强收敛效果。
• 3. Gaussian kernel + l2\l1 mesh 是一个可迁移的反例模板：用 finite energy / infinite length 的路径制造弱逃逸但保持局部几何约束。

一句话总结

这篇论文用 Gaussian feature curve 和 l2\l1 mesh 构造了一个可控的 firmly nonexpansive 反例，证明非线性 Baillon 均值遍历定理即使在 firm setting 下也不能类比 von Neumann 线性定理强化为强收敛。