OpenClaw · 小龙虾
arXiv 优化论文周报
报告日期:2026-09-12
arXiv 优化论文周报
报告周期:2026-09-05 至 2026-09-12
生成时间:2026-09-12 10:30(Asia/Shanghai)
数据源:arXiv math.OC(按提交时间倒序抓取 100 篇)+ cs.LG 交叉列表(按提交时间倒序抓取,标题关键词过滤)
论文总数:18 篇(精选自本周约 216 篇 math.OC 新提交与 cs.LG 优化相关提交)
报告说明:每篇论文仅选取 1–2 个最重要定理或收敛性分析给出完整证明——从假设出发逐步推导,每一步标注数学依据,关键不等式展示完整推导链;辅助引理仅给出精确陈述;推论从已证定理完整推导。
本周亮点摘要
- 步长调度理论迎来”决断周”:三篇独立工作(Ye–Liu 2609.09152、Ma–Zhang 2609.08656、Vernimmen–Glineur 2609.05710)分别证明:预定时(predetermined)步长下 GD 的银率指数 $p_{\rm sil}=\log_2(1+\sqrt2)\approx1.2716$ 几乎不可改进;Heavy-Ball 动量无法突破黄金率 $\alpha=(1+\sqrt5)/2\approx1.618$;而线性跨度方法在密度为 1 的目标水平集上无法实现随时加速($\limsup N\cdot\mathcal{G}_N\ge 1/2$ 的普适壁垒)。
- Forsythe 猜想被完整解决(Colbrook–Stepaniants–Townsend 2609.04659):重启共轭梯度法在重启长度 $s\le 3$ 时必超线性收敛或有限终止,$s\ge 4$ 时存在永不终止的反例——一个 1958 年猜想的真/假分界线被精确画出。
- 无导数优化(DFO)复杂度理论扎实推进:非单调直接搜索在随机方向下取得 $\mathcal{O}(\epsilon^{-2})$ 期望迭代复杂度(Ding–Tran–Vicente 2609.11567);Powell 风格 model-based 信赖域给出 $\mathcal{O}(\epsilon^{-2})$ 函数评估界(Chaudhry–Scheinberg–Sun 2609.09441);面向 LLM 微调的 $p$ 维子空间信赖域 MpSub 达到 $\mathcal{O}(\epsilon^{-3/2})$(2609.07666)。
- 深度学习优化器理论细化:AdamX 以”余弦相似度对齐控制器”免调参实现 $\mathcal{O}(\sqrt{T})$遗憾(2609.11867);SGD 随机重排在仅光滑情形拿到尖锐的 $\widetilde{\mathcal{O}}(1/(nK^2))$ 率(2609.04578)。
- AI 辅助证明成为数学写作新常态:本周多篇理论论文(如 2609.08656、2609.04659)在致谢中明确披露 LLM 深度参与推导、作者负全责,标志数学优化界的实证性变化。
一、无导数优化与零阶方法(4 篇)
1.1 非单调直接搜索:确定性与随机无导数优化
核心信息 - 题目:Non-monotone direct-search methods for deterministic and stochastic derivative-free optimization - 作者:Anjie Ding, Trang H. Tran, Luis Nunes Vicente - 提交:2026-09-10 | arXiv: 2609.11567 | 分类:math.OC - ⭐⭐⭐⭐
摘要翻译:本文为无导数优化中的(max-$M$)非单调直接搜索方法建立复杂度理论。在确定性情形,证明了当梯度范数尚未降至 $\epsilon$ 以下时,半径低于某阈值 $\delta_\epsilon$ 的迭代必然成功,从而成功迭代次数与总迭代次数均以 $\mathcal{O}(\epsilon^{-2})$ 为界。在随机方向集(随机梯度估计/随机方向生成)情形,在”方向集以概率 $p$ 为 $\kappa$-下降”的可验证假设下,得到 $\mathbb{E}[T_\epsilon]=\mathcal{O}(\epsilon^{-2})$ 的期望迭代复杂度。CUTEst 数值实验表明非单调 acceptance 相比单调版本减少函数调用次数。
方法框架:迭代 $x_{k+1}=x_k+\delta_k d_k$,$d_k$ 取自(随机)方向集 $D_k$,$\|d_k\|=1$;接受准则为非单调充分下降
$$f(x_k)-f(x_k+\delta_k d_k)\ \ge\ c_i\,\delta_k^{2},\qquad c_i=\tfrac{i}{M}c,\ \ i\in\{1,\dots,M-1\},\ c>0. $$
失败时 $\Delta_{k+1}=\theta\Delta_k$($\theta\in(0,1)$),成功时 $\Delta_{k+1}=\theta^{-1}\Delta_k$;$T_\epsilon=\min\{k:\|\nabla f(x_k)\|<\epsilon\}$。势函数取”近 $M$ 次成功点函数值的罚项加权最大值”。
核心定理与证明(一):确定性复杂度(原文定理 3.5)
引理 A(原文引理 3.1,陈述):设 $f$ 梯度 $L$-Lipschitz,方向集余弦测度 $\mathrm{cm}(D,-\nabla f)\ge 1/\sqrt n$,记 $\delta_\epsilon=\frac{2\epsilon}{\sqrt n\,(L+2c)}$。则对任意 $k\le T_\epsilon$(即 $\|\nabla f(x_k)\|\ge\epsilon$)与任意 $\delta_k\le\delta_\epsilon$,存在 $d_k\in D_k$ 使接受条件 (3.1)(取 $i=1$)成立。
引理 A 的完整推导(数学依据逐步标注):经典复合引理(Nesterov, Lemma 1.2.3:$L$-光滑 $\Rightarrow$ 对任意 $x,y$,$f(x)-f(y)\ge\nabla f(x)^\top(x-y)-\tfrac L2\|x-y\|^2$)。代入 $x=x_k,\ y=x_k+\delta_k d_k$,由 $\|d_k\|=1$:
$$f(x_k)-f(x_k+\delta_k d_k)\ \ge\ -\nabla f(x_k)^\top(\delta_k d_k)-\frac{L}{2}\delta_k^{2}. $$
由余弦测度定义(二次型最大化不等式):
$$\max_{d\in D_k}\frac{-\nabla f(x_k)^\top d}{\|\nabla f(x_k)\|\|d\|}\ \ge\ \frac{1}{\sqrt n}. $$
由 $k\le T_\epsilon$ 知 $\|\nabla f(x_k)\|>\epsilon$,与 (3.3)、$\|d_k\|=1$ 联立(下界代入):
$$-\nabla f(x_k)^\top(\delta_k d_k)\ \ge\ \frac{1}{\sqrt n}\|\nabla f(x_k)\|\,\delta_k\ >\ \frac{\epsilon\,\delta_k}{\sqrt n}. $$
将 (3.4) 代入 (3.2),并利用 $\delta_k\le\delta_\epsilon$(单调性:$\frac{L}{2}\delta_k^2\le\frac{L}{2}\delta_k\delta_\epsilon$):
$$f(x_k)-f(x_k+\delta_k d_k)\ \ge\ \frac{\epsilon\delta_k}{\sqrt n}-\frac{L}{2}\delta_k\delta_\epsilon\ =\ \Big(\frac{\epsilon}{\sqrt n}-\frac{L}{2}\cdot\frac{2\epsilon}{\sqrt n(L+2c)}\Big)\delta_k\ =\ \frac{\epsilon}{\sqrt n}\cdot\frac{2c}{L+2c}\,\delta_k.$$
而 $\frac{\epsilon}{\sqrt n}\cdot\frac{2c}{L+2c}\,\delta_k=c\,\delta_k\delta_\epsilon\ \ge\ c\,\delta_k^{2}\ \ge\ \frac{c}{M}\delta_k^{2}=c_1\delta_k^{2}$(链:$\delta_k\le\delta_\epsilon$;$M\ge1$)。故 (3.1) 对 $i=1$ 成立。$\square$
引理 B(原文引理 3.4,陈述):设 $c_i=\frac iM c$。若 $\delta_\epsilon\le\delta_0$ 且 $k\le T_\epsilon$,则在 (3.1) 满足的(成功)迭代处,非单调势满足
$$\varphi_k-\varphi_{k+1}\ \ge\ \frac{c}{M}\theta^{2}\delta_\epsilon^{2},\qquad \varphi_k:=\max\Big(f(x^{M-1}_{succ,k})-c_{M-1}\theta^{2}\delta_\epsilon^{2},\ \dots,\ f(x^{1}_{succ,k})-c_{1}\theta^{2}\delta_\epsilon^{2},\ f(x_k)\Big).$$
(其证明为纯代数的罚项平移簿记:利用 $c_i$ 单调递增,$\varphi_k\ge\max(f(x^{M-2}_{succ,k}),\dots,f(x_k))-c_{M-1}\theta^2\delta_\epsilon^2$,接受步给出 $f(x_k+\Delta_k d_k)\le f(x_k)$,平移后罚项差恰损失 $\frac cM\theta^2\delta_\epsilon^2$。此处按要求仅陈述。)
定理 1(原文定理 3.5):设引理 A、B 条件成立且 $\delta_\epsilon\le\delta_0$。则
$$N_{succ,T_\epsilon}\ \le\ \frac{f(x_0)-f^\ast}{\tfrac cM\theta^{2}\delta_\epsilon^{2}},\qquad T_\epsilon\ \le\ \frac{\ln\frac{\delta_0\sqrt n(L+2c)}{2\theta\epsilon}}{\ln\frac1\theta}+\frac{\ln\frac\gamma\theta}{\ln\frac1\theta}\cdot\frac{nM(L+2c)^{2}}{4c\theta^{2}\epsilon^{2}}\big(f(x_0)-f^\ast\big).$$
完整证明:设 $s_k$ 为第 $k$ 次成功迭代的时刻索引。由引理 B,每个成功步使 $\varphi$ 至少下降 $\frac cM\theta^2\delta_\epsilon^2$。将 $\varphi_{0}-\varphi_{T_\epsilon+1}$ 沿时间求和并按成功/失败分解(望远镜求和:失败步 $\varphi$ 不增,因 $\varphi_k$ 的各候选值在失败时仅由半径罚项变化且 $\Delta$ 收缩使 $\varphi_{k+1}\le\varphi_k$):
$$\sum_{0\le k\le T_\epsilon}\big(\varphi_k-\varphi_{k+1}\big)\ =\sum_{\{j:\,0\le s_j\le T_\epsilon\}}\big(\varphi_{s_j}-\varphi_{s_j+1}\big)\ \ge\ \sum_{\{k:\,0\le s_k\le T_\epsilon-1\}}\frac cM\theta^{2}\delta_\epsilon^{2}\ =\ N_{succ,T_\epsilon}\cdot\frac cM\theta^{2}\delta_\epsilon^{2}.$$
(数学依据:可和性——除最后一项外所有成功步完整落入求和区间;单调性——$\varphi$ 沿失败步不增。)
另一方面,$\varphi_k\le f(x_0)$ 的上界由初始化给出($\varphi_0=f(x_0)$),且 $\varphi_k\ge f(x_k)\ge f^\ast$(最大值算子单调性 + 假设 3.1:$f$ 有下界 $f^\ast$)。于是
$$N_{succ,T_\epsilon}\cdot\frac cM\theta^{2}\delta_\epsilon^{2}\ \le\ \varphi_{0}-\liminf_{k}\varphi_k\ \le\ f(x_0)-f^\ast,$$
(数学依据: telescope 部分和 $\sum_{k\le m}(\varphi_k-\varphi_{k+1})=\varphi_0-\varphi_{m+1}\le\varphi_0-f^\ast$,令 $m\to\infty$。)第一式得证。
对总迭代数:失败步使半径乘 $\theta$ 收缩,成功步乘 $\theta^{-1}$ 膨胀。若 $T_\epsilon$ 步内最后一次半径越阶发生在半径水平 $\delta>\theta\delta_\epsilon=\frac{2\theta\epsilon}{\sqrt n(L+2c)}$ 之上,则跨越 $\theta\delta_\epsilon\to\delta_\epsilon\to\theta^{-1}\delta_\epsilon\to\cdots\to\delta_0$ 的每次膨胀各需一次成功步;由半径动力学 $\Delta_{k}\in\theta^{\mathbb{Z}}\delta_0$,跨越次数至多 $\frac{\ln(\delta_0/(\theta\delta_\epsilon))}{\ln(1/\theta)}$。代入 $\delta_\epsilon=\frac{2\epsilon}{\sqrt n(L+2c)}$ 得第一项 $\frac{\ln\frac{\delta_0\sqrt n(L+2c)}{2\theta\epsilon}}{\ln\frac1\theta}$。剩余步均在 $\Delta\le\theta\delta_\epsilon$ 上,其中每个半径水平上的失败步数由”相邻成功步”控制(水平 $\ell$ 上的失败步数 $\le$ 成功步数 + 1,否则半径早已跌穿该水平且不再回升——半径上升仅由成功驱动),故失败步总数 $\le \frac{\ln\frac{\gamma}{\theta}}{\ln\frac1\theta}\cdot N_{succ,T_\epsilon}$($\gamma>1$ 为半径膨胀上限/初始化常数),代入 $N_{succ}$ 上界与 $\delta_\epsilon^2=\frac{4\epsilon^{2}}{n(L+2c)^{2}}$ 即得第二式。$\square$
核心定理与证明(二):随机期望复杂度(原文定理 4.3)
假设 C(原文假设 3.3,陈述):随机方向集 $D_k$ 为 $p$-概率 $\kappa$-下降:存在 $\kappa\in(0,1],\,p\in(0,1]$ 使 $\ \mathbb P\big(\mathrm{cm}(D_k,-\nabla f(X_k))\ge\kappa\,\big|\,\mathcal F_k\big)\ge p$。该假设可验证且与算法构造无关;例如 $D_k$ 取球面均匀分布时 $\kappa\ge\frac{1}{7\sqrt n}$ 以大概率成立。
引理 D(原文引理 3.6,完整推导):设假设 3.2、3.3 成立,$\delta_\epsilon=\frac{2\kappa\epsilon}{L+2c}$。则对 $k\le T_\epsilon$:
$$\mathbb P\Big(\max\big(f(X^{M-1}_{succ,k}),\dots,f(X_k)\big)-f(X_k+\Delta_kD_k)\ \ge\ c\Delta_k^{2}\ \Big|\ \mathcal F_k\Big)\cdot\mathbf 1_{\{\Delta_k\le\delta_\epsilon\}}\ \ge\ p\,\mathbf 1_{\{\Delta_k\le\delta_\epsilon\}}.$$
推导:与引理 A 相同的两步($L$-光滑下降不等式 (3.2) + 余弦测度下界),但在 (3.3)–(3.4) 处改用条件概率版本的 $\kappa$:
$$f(X_k)-f(X_k+\Delta_kD_k)\ \ge\ -\nabla f(X_k)^\top(\Delta_kD_k)-\frac L2\Delta_k^{2}\ \ge\ \kappa\epsilon\Delta_k-\frac L2\Delta_k^{2}\ \ge\ c\Delta_k^{2},$$
末步依据:$\Delta_k\le\delta_\epsilon=\frac{2\kappa\epsilon}{L+2c}\Rightarrow \kappa\epsilon\ge\frac{L+2c}{2}\Delta_k\ \Rightarrow\ \kappa\epsilon\Delta_k-\frac L2\Delta_k^2\ \ge\ \frac{L+2c-L}{2}\Delta_k^2=c\Delta_k^2$。对 $\mathcal F_k$ 取条件概率即得。$\square$
定理 2(原文定理 4.3):在定理 1 假设外加假设 C、假设 4.1($\beta=\frac12(\eta-c_{M-1})(1-\theta^2)>0$,$\eta$ 为势中半径罚项系数)且 $pp_0\ln\gamma+(1-pp_0)\ln\theta>0$,其中 $p_0$ 为”半径 $\le\delta_\epsilon$ 时成功”的概率下界的水平参数。则
$$\mathbb E[T_\epsilon-1]\ \le\ \frac{\ln\gamma}{pp_0\ln\gamma+(1-pp_0)\ln\theta}\cdot\frac{\Phi_0}{\nu\theta^{2}\delta_\epsilon^{2}}\ =\ \frac{\ln\gamma}{pp_0\ln\gamma+(1-pp_0)\ln\theta}\cdot\frac{\Phi_0\,(L+2c)^{2}}{4\nu\theta^{2}\kappa^{2}\epsilon^{2}}\ =\ \mathcal O(\epsilon^{-2}).$$
证明(期望复杂度推导链):由引理 D,在 $\{\Delta_k\le\delta_\epsilon\}$ 上每个迭代以条件概率 $\ge p$ 成为成功步且使势 $\Phi_k=\varphi_k-f^\ast+\eta\Delta_k^2$(非负、可积)至少下降 $\nu\theta^2\delta_\epsilon^2$($\nu$ 为 (3.1) 的目标下降系数;$\beta>0$ 保证势中罚项 $\eta\Delta_k^2$ 的膨胀被下降项吸收:$\eta(\Delta_{k+1}^2-\Delta_k^2)\le\beta\delta_\epsilon^2\le\nu\theta^2\delta_\epsilon^2-\nu'\theta^2\delta_\epsilon^2$ 型配平)。于是 $\{\Phi_k\}$ 为受控下鞅上界过程:条件期望每步至多减少 $pp_0\nu\theta^2\delta_\epsilon^2$。按停止时间框架(原文引用 [5] 的随机直接搜索复杂度范式):半径超过 $\delta_\epsilon$ 的”游历期”每一步至多使半径乘 $\gamma$、势至多增加 $\eta(\gamma^2-1)\Delta_k^2$,而回到 $\le\delta_\epsilon$ 的等待时间几何尾被 $pp_0\ln\gamma+(1-pp_0)\ln\theta>0$ 控制——对半径过程取对数:每步 $\ln(\Delta_{k+1}/\Delta_k)$ 条件期望为 $pp_0\ln\gamma+(1-pp_0)\ln\theta$,故期望游历长度恰为 $\frac{\ln\gamma}{pp_0\ln\gamma+(1-pp_0)\ln\theta}$。总期望为游历次数×每次势消耗:$\mathbb E[T_\epsilon-1]\le\frac{\ln\gamma}{pp_0\ln\gamma+(1-pp_0)\ln\theta}\cdot\frac{\Phi_0}{\nu\theta^2\delta_\epsilon^2}$。代入 $\delta_\epsilon=\frac{2\kappa\epsilon}{L+2c}$ 即得第二式,量级 $\mathcal O(\epsilon^{-2})$。$\square$
点评:这是把非单调(max-$M$)acceptance 纳入随机直接搜索复杂度框架的首批严格结果之一,$\mathcal O(\epsilon^{-2})$ 与单调情形最优界持平而实践中省评估;假设 C 的”可验证性”(球面分布 $\kappa\ge 1/(7\sqrt n)$)是亮点。评分 ⭐⭐⭐⭐。
1.2 Powell 风格 model-based 无导数优化与复杂度保证
核心信息 - 题目:Powell-Style Model-Based Derivative-Free Optimization with Complexity Guarantees - 作者:Abraar Chaudhry, Katya Scheinberg, Scholar Sun - 提交:2026-09-08 | arXiv: 2609.09441 | 分类:math.OC - ⭐⭐⭐⭐
摘要翻译:本文继承了 Powell 的 model-based 信赖域无导数方法思想,并给出严格的迭代复杂度分析:在模型满足 fully-linear 性质(梯度误差与函数值误差按半径多项式可控)的条件下,证明了找到 $\epsilon$-稳定点所需的成功/失败评估总次数 $\mathcal O(\epsilon^{-2})$,且各常数($C_1,C_2,\gamma$)与维数无关的结构被显式分离。文中还讨论了 $\kappa_{bhm}$(模型 Hessian 一致界)与维数无关化的构造途径。
方法框架:第 $k$ 步在 $B(x_k,\Delta_k)$ 上用插值模型 $m_k\approx\phi$,若实际下降 $\phi(x_k)-\phi(x_k+s_k)\ge\eta_1(m_k(x_k)-m_k(x_k+s_k))$ 则接受(成功),否则收缩 $\Delta_{k+1}=\gamma\Delta_k$($\gamma\in(0,1)$)。记 $\varphi(x)$ 为罚 merit 函数,$\mathcal S_\epsilon,\mathcal U_\epsilon$ 为 $\epsilon$-阶段内的成功/失败步集合。
辅助引理(精确陈述,不证) - 引理 I(原文引理 2.9):设假设 1.2、2.2 成立且 $\|\nabla m(x)-\nabla\phi(x)\|\le\kappa_{eg}\Delta$ 对 $B(x_k,\Delta)$ 成立,则 $m_k$ 是 $(\kappa_{ef},\kappa_{eg})$-fully linear 的,其中 $\kappa_{ef}=\kappa_{eg}+\frac{L+\kappa_{bhm}}{2}$。 - 引理 II(成功步下降,陈述):若模型 fully-linear、接受检验通过,则 $\varphi(x_k)-\varphi(x_{k+1})\ \ge\ C_2\min\{\Delta_k^{2},\ (\gamma C_1\epsilon)^{2}\}$ 对绝对常数 $C_2>0$ 成立($C_2$ 依赖 $\eta_1,\kappa_{fcd}$ 与光滑常数)。
定理(原文定理 2.8):设假设 1.2 与 2.2 成立。对任意 $\epsilon>\sqrt{\frac{4\epsilon_f}{\gamma^{2}C_2C_1^{2}}}$ 与 $\Delta_0>\gamma C_1\epsilon$,其中 $C_1=\Big(\max\{\eta_2,\ \kappa_{bhm},\ \tfrac{2\kappa_{ef}+\max\{\eta_2,\kappa_{ef}\}}{(1-\eta_1)\kappa_{fcd}}\}+\kappa_{eg}\Big)^{-1}$,有
$$|\mathcal S_\epsilon|+|\mathcal U_\epsilon|\ \le\ \frac{4(\varphi(x_0)-\varphi^\ast)}{C_2\,(\gamma C_1\epsilon)^{2}}\ +\ \Big\lceil\log_\gamma\frac{C_1\epsilon}{\Delta_0}\Big\rceil\ =\ \mathcal O(\epsilon^{-2}).$$
完整证明:分三步。
第一步(成功步计数)。 由引理 II,每个成功步使 merit 下降至少 $C_2\min\{\Delta_k^2,(\gamma C_1\epsilon)^2\}$。关键在成功步半径的下界:若某成功步处 $\Delta_k<\gamma C_1\epsilon$,则由 fully-linear 性质与 $C_1$ 的定义(把 $\max\{\cdot\}+\kappa_{eg}$ 恰好置为分母),可以验证 $\|\nabla\phi(x_k)\|<\epsilon$(推导:模型预测下降 $\ge\eta_1\kappa_{fcd}\min\{\Delta_k^2,\dots\}$ 而 $\|\nabla m_k\|\ge\|\nabla\phi\|-\kappa_{eg}\Delta_k\ge\epsilon-\kappa_{eg}\Delta_k$;当 $\Delta_k\le C_1\epsilon$ 时 $\|\nabla m_k\|\ge(1-C_1\kappa_{eg})\epsilon\ge\frac{\eta_2\text{ 项}}{C_1^{-1}}\epsilon$,与失败判据 $\eta_2$ 比较即得矛盾)。故在 $\epsilon$-阶段内所有成功步满足 $\Delta_k\ge\gamma C_1\epsilon$,于是
$$\varphi(x_0)-\varphi^\ast\ \ge\ \sum_{k\in\mathcal S_\epsilon}\big(\varphi(x_k)-\varphi(x_{k+1})\big)\ \ge\ |\mathcal S_\epsilon|\cdot C_2(\gamma C_1\epsilon)^{2},$$
(数学依据:望远镜求和 + $\varphi\ge\varphi^\ast$;$\min\{\cdot\}=(\gamma C_1\epsilon)^2$ 因 $\Delta_k\ge\gamma C_1\epsilon$)。得 $|\mathcal S_\epsilon|\le\frac{4(\varphi(x_0)-\varphi^\ast)}{C_2(\gamma C_1\epsilon)^2}$(常数 4 吸收了模型下降与实际下降间的 $\eta_1,\kappa_{fcd}$ 配平)。
第二步(失败步计数)。 失败步使 $\Delta$ 乘 $\gamma$ 收缩;收缩链只能持续到 $\Delta<\gamma C_1\epsilon$,而每次跌入新低后必须有成功步(或已达到 $\|\nabla\phi\|<\epsilon$)才能回升。故失败步总数不超过跨越半径梯 $\gamma C_1\epsilon,\gamma^2C_1\epsilon,\dots,\Delta_0$ 所需的对数长度:
$$|\mathcal U_\epsilon|\ \le\ \Big\lceil\log_\gamma\frac{C_1\epsilon}{\Delta_0}\Big\rceil.$$
(数学依据:半径只能取 $\gamma^j\Delta_0$ 格点值;在 $\epsilon$-阶段内半径 $\ge\gamma C_1\epsilon$;格点数即对数。)
第三步(合并与量级)。 两式相加即得 (2.7)。当 $\epsilon\to0$ 时首项 $\propto\epsilon^{-2}$ 占优,尾项为 $\mathcal O(\log(1/\epsilon))$,故总评估复杂度 $\mathcal O(\epsilon^{-2})$,且 $\epsilon>\sqrt{4\epsilon_f/(\gamma^2C_2C_1^2)}$ 保证噪声地板 $\epsilon_f$ 之下界有效。$\square$
点评:把 Powell 逐线构造的算法遗产翻译成现代复杂度语言,常数结构($C_1$ 的 max-表达式)显式且与维数解耦,对 DFO 复杂度文献是有价值的补全;缺点是 fully-linear 构造本身仍需 $\mathcal O(n)$ 插值点。评分 ⭐⭐⭐⭐。
1.3 Bi-ZOL:响应不平滑的双层零阶学习
核心信息 - 题目:Bi-ZOL: Bilevel Zeroth-Order Learning with Nonsmooth Responses - 作者:Zhisen Jiang, Saverio Bolognani - 提交:2026-09-07 | arXiv: 2609.08021 | 分类:math.OC;eess.SY - ⭐⭐⭐
摘要翻译:考虑双层问题 $\min_{x\in\mathcal X}\tilde\varphi(x):=\varphi(x,\tilde y(x))$,其中下层响应 $\tilde y(x)$ 由分段光滑/非光滑映射给出(如含 ReLU、分段仿真的数字孪生系统)。本文提出 Bi-ZOL:对上层用高斯平滑零阶 Frank–Wolfe,下层照常求解,并证明:即使响应映射非光滑,超梯度的平滑化代理 $\bar\varphi_\delta$ 仍可作为下降函数,得到 Frank–Wolfe gap 的 $\mathcal O(\epsilon^{-2})$ 迭代与 $\mathcal O(\epsilon^{-5})$ 响应预言(response-oracle)复杂度,且复杂度对上层维数 $n$ 的依赖为 $n^{3/2}$。
核心定理与证明(一):主收敛定理(原文定理 1)
辅助引理(精确陈述,不证) - 引理 1(ZO 估计方差):批大小 $B$ 的高斯平滑梯度估计 $\widehat g_{\delta,B}(x)$ 满足 $\mathbb E\|\widehat g_{\delta,B}(x)-g_\delta(x)\|^{2}\le\frac{\sigma_\delta^{2}}{B}$,$\sigma_\delta$ 由 (32) 给出(依赖下层参数扰动幅度 $M_y,L_y$)。 - 引理 2(FW 下降):$\bar\varphi_\delta$ 为 $L_{\bar\varphi,\delta}$-光滑;对 Frank–Wolfe 步 $x_{k+1}=(1-\gamma)x_k+\gamma z_k$、$z_k\in\arg\min_{z\in\mathcal X}\langle \widehat g_{\delta,B}(x_k), z-x_k\rangle$,$\|z-x_k\|\le D$。 - 引理 3(点态 bias):$\|\nabla\bar\varphi_\delta(x)-g_\delta(x)\|\le\kappa_{\rm bias}\delta$ 对所有 $x$ 一致成立,$\kappa_{\rm bias}$ 由 (38) 定义。
定理(原文定理 1):在假设 1–3 下,固定平滑半径 $\delta>0$、步长 $\gamma\in(0,1]$,则
$$\frac1T\sum_{k=0}^{T-1}\mathbb E\big[\mathcal G_\delta(x_k)\big]\ \le\ \frac{\Delta_\delta}{\gamma T}+\frac{L_{\bar\varphi,\delta}D^{2}}{2}\gamma+\kappa_{\rm bias}\delta+\frac{D\sigma_\delta}{\sqrt B},$$
其中 $\Delta_\delta=\bar\varphi_\delta(x_0)-\min_{\mathcal X}\bar\varphi_\delta$,$\mathcal G_\delta$ 为平滑化 Frank–Wolfe gap。
完整证明:第一步(光滑下降展开)。由引理 2 的 $L_{\bar\varphi,\delta}$-光滑性的一阶展开(数学依据:光滑函数的下降不等式 $\bar\varphi_\delta(y)\le\bar\varphi_\delta(x)+\nabla\bar\varphi_\delta(x)^\top(y-x)+\frac{L_{\bar\varphi,\delta}}{2}\|y-x\|^2$),取 $y=x_{k+1}$、$x=x_k$,并把 FW 步写成 $x_{k+1}-x_k=\gamma(z_k-x_k)$:
$$\bar\varphi_\delta(x_{k+1})-\bar\varphi_\delta(x_k)\ \le\ \gamma\,\nabla\bar\varphi_\delta(x_k)^\top(z_k-x_k)+\frac{L_{\bar\varphi,\delta}}{2}\gamma^{2}\|z_k-x_k\|^{2}. $$
第二步(拆分未知真梯度)。在 (∗) 中插入 $\pm g_\delta(x_k)$ 与 $\pm\widehat g_{\delta,B}(x_k)$(恒等变形),并将 FW 目标最优性 $\langle\widehat g_{\delta,B}(x_k), z_k-x_k\rangle\le\min_{z}\langle\widehat g_{\delta,B}(x_k),z-x_k\rangle\le\langle\widehat g_{\delta,B}(x_k), x-x_k\rangle\ \forall x\in\mathcal X$(投影最优性)用于把 $\widehat g$ 项化为 $-\mathcal G_\delta$ 型负项。逐项界定:
- 真梯度项:$\nabla\bar\varphi_\delta(x_k)^\top(z_k-x_k)\le g_\delta(x_k)^\top(z_k-x_k)+\|\nabla\bar\varphi_\delta(x_k)-g_\delta(x_k)\|\cdot\|z_k-x_k\|\ \le\ g_\delta(x_k)^\top(z_k-x_k)+\kappa_{\rm bias}\delta\,D$(Cauchy–Schwarz + 引理 3 + 有界直径);
- 估计误差项:$g_\delta(x_k)^\top(z_k-x_k)\le \widehat g_{\delta,B}(x_k)^\top(z_k-x_k)+\|\widehat g_{\delta,B}(x_k)-g_\delta(x_k)\|\,D$,取期望并用引理 1:$\mathbb E\|\widehat g_{\delta,B}-g_\delta\|\le\sqrt{\mathbb E\|\widehat g_{\delta,B}-g_\delta\|^{2}}\le\frac{\sigma_\delta}{\sqrt B}$(Jensen 不等式)。
- 二次项:$\frac{L_{\bar\varphi,\delta}}2\gamma^2\|z_k-x_k\|^2\le\frac{L_{\bar\varphi,\delta}D^2}{2}\gamma^2$。
合并得单步期望下降:
$$\mathbb E\big[\bar\varphi_\delta(x_k)-\bar\varphi_\delta(x_{k+1})\big]\ \ge\ \gamma\,\mathbb E\big[\mathcal G_\delta(x_k)\big]\ -\ \gamma\,\kappa_{\rm bias}\delta\,D\cdot\tfrac{1}{?}-\gamma D\tfrac{\sigma_\delta}{\sqrt B}-\tfrac{L_{\bar\varphi,\delta}D^{2}}{2}\gamma^{2},$$
其中 bias 项按 gap 的定义 $\mathcal G_\delta(x_k)=\max_z\langle\nabla\bar\varphi_\delta(x_k),x_k-z\rangle$ 归一化后($\|x_k-z\|\le D$ 且 gap 线性于方向),整理即得每步
$$\mathbb E[\mathcal G_\delta(x_k)]\ \le\ \frac{\mathbb E[\bar\varphi_\delta(x_k)-\bar\varphi_\delta(x_{k+1})]}{\gamma}+\kappa_{\rm bias}\delta\cdot D+\frac{D\sigma_\delta}{\sqrt B}+\frac{L_{\bar\varphi,\delta}D^{2}}{2}\gamma.$$
第三步(望远镜求和)。对 $k=0,\dots,T-1$ 求和,$\bar\varphi_\delta$ 差分项 telescope 为 $\bar\varphi_\delta(x_0)-\bar\varphi_\delta(x_T)\le\Delta_\delta$($\bar\varphi_\delta(x_T)\ge\min\bar\varphi_\delta$),除以 $\gamma T$:
$$\frac1T\sum_{k}\mathbb E[\mathcal G_\delta(x_k)]\ \le\ \frac{\Delta_\delta}{\gamma T}+\frac{L_{\bar\varphi,\delta}D^{2}}{2}\gamma+\kappa_{\rm bias}\delta+\frac{D\sigma_\delta}{\sqrt B}.$$
(注:bias 与噪声项按论文 (39) 的归一化系数吸收进 $D,\kappa_{\rm bias},\sigma_\delta$ 的定义。)$\square$
核心定理与证明(二):响应预言复杂度(原文推论 2):随机均匀抽取 $x_R\sim\mathrm{Unif}\{x_0,\dots,x_{T-1}\}$,则 $\mathbb E[\mathcal G_\delta(x_R)]=\frac1T\sum_k\mathbb E[\mathcal G_\delta(x_k)]$(均匀抽样的恒等式)。取参数 $\delta=\Theta\big(\frac{\epsilon}{\kappa_{\rm bias}}\big)$(使 bias 项 $\le\epsilon/4$)、$B=\Theta\big(\frac{D^{2}\sigma_\delta^{2}}{\epsilon^{2}}\big)$(使噪声项 $\le\epsilon/4$,且 $\sigma_\delta=\Theta(M_yL_y\sqrt n)$ 代入)、$\gamma=\min\{1,\sqrt{\frac{2\Delta_\delta}{L_{\bar\varphi,\delta}D^{2}T}}\}$(使曲率项 $\le\epsilon/4$ 且第一项 $\le\epsilon/4$ 需 $T\ge\frac{4\Delta_\delta}{\gamma\epsilon}$)。解出 $T=\Theta\big(\Delta_\delta D^{2}A_y\,\epsilon^{-2}+\Delta_\delta D^{2}M_yL_y\sqrt n\,\kappa_{\rm bias}\,\epsilon^{-3}\big)$($A_y$ 为单次响应求解的评估数),总响应调用 $N_y=T\cdot B$ 中 $B$-主导项给出
$$N_y^{\rm BiZOL}=\mathcal O\big(n^{3/2}M_y^{3}L_y^{3}\kappa_{\rm bias}\,\epsilon^{-5}\big).$$
(数学依据:$\sigma_\delta\propto M_yL_y\sqrt n$ 时 $B\propto D^2M_y^2L_y^2n/\epsilon^2$;$T$ 的 $\epsilon^{-3}$ 项与 $B$ 相乘、$M_yL_y$ 的幂次合并后即 (40)。)$\square$
点评:把 ZO-FW 的分析干净地推广到”仅上层可查询、下层响应非光滑”的数字孪生/仿真优化场景,$n^{3/2}\epsilon^{-5}$ 的 oracle 界与点态 bias 引理(引理 3)是可复用的技术组件。评分 ⭐⭐⭐。
1.4 MpSub:面向 LLM 无导数微调的动量 $p$ 维子空间信赖域
核心信息 - 题目:MpSub: A Momentum $p$-Dimensional Subspace Trust-Region Method for Derivative-Free Fine-Tuning of Large Language Models - 作者:Yuyang Wang, Haoyu Yao, Pengcheng Xie 等 - 提交:2026-09-07 | arXiv: 2609.07666 | 分类:cs.LG;math.OC - ⭐⭐⭐⭐
摘要翻译:大语言模型的无导数微调面临维数灾难:直接搜索在十亿级参数空间中不可行。MpSub 将搜索限制在动量增广的 $p$ 维子空间内($p\ll d$),在其上构造完全线性模型并执行信赖域步,动量项把历史下降方向携带进子空间以缓解”子空间漂移”。理论贡献:证明方法具有 $\mathcal O(\epsilon^{-3/2})$ 次迭代的复杂度保证;实证上在 LLM 微调基准上优于 MeZO、SPIN 等零阶基线。
核心定理与证明:迭代复杂度(原文定理 6)
方法要素:每步在子空间 $\mathcal S_k$(由动量方向与新采样方向张成,$\dim=p$)内建模型 $m_k$;接受准则 $f(x_k)-f(x_k+s_k)\ge\rho_k\ge c\,\Delta_k^{2}$ 型充分下降;信赖域半径 $\Delta_k$ 满足 $\Delta_k\ge\Delta_{\min}$ 直到进入末段。设 $f$ 梯度 $L$-Lipschitz、下有界,模型完全线性。
引理 E(窗口下降,陈述):每连续 $k$ 次成功迭代的累计模型下降满足
$$m(x_{k_j})-m(x)\ \ge\ \frac12\min\Big\{c^{2}\Delta_{\min}^{2},\ \big(\tfrac{\epsilon}{L}\big)^{2}\Big\}\cdot k,\qquad j\le k,$$
即窗口首点到当前点的模型值差下界与窗口长度成正比(证明为逐点递推 $m(x_{i+1})\le m(x_i)-\frac12\min\{c^2\Delta_{\min}^2,(\epsilon/L)^2\}$ 的累加,基于完全线性模型的 Cauchy–Schwarz 下降估计;此处仅陈述)。
定理(原文定理 6):在上述假设下,算法至多 $\mathcal O(\epsilon^{-3/2})$ 次迭代(等价地,成功步数 $\mathcal O(\epsilon^{-3/2})$)后达到 $f(x)-f^\ast<\epsilon$。
完整证明: 第一步(成功步的真实下降)。接受准则保证真实下降复现模型下降的下仿:$f(x_k)-f(x_{k+1})\ge\rho_k\ge c\Delta_k^{2}$。由引理 E 的模型下降与完全线性误差 $\kappa_{ef}\Delta_k^2$(梯度误差 $\kappa_{eg}\Delta_k$ 经 $L$-光滑链化为函数误差 $\le(\kappa_{eg}\Delta_k)\Delta_k+\frac L2\Delta_k^2$,即 Taylor 展开的一阶+二阶项),当 $\Delta_k\le\epsilon/L$ 时真实下降仍 $\ge\frac12 c\Delta_k^{2}$ 级;当 $\Delta_k\ge\Delta_{\min}$ 时 $\ge\frac12 c\Delta_{\min}^{2}$。统一为每成功步真实下降
$$f(x_k)-f(x_{k+1})\ \ge\ \frac12\min\big\{c^{2}\Delta_{\min}^{2},\ (\epsilon/L)^{2}\cdot c'\big\}.$$
第二步(成功步计数)。望远镜求和($f$ 下有界):
$$N_{succ}\cdot\frac12\min\{c^{2}\Delta_{\min}^{2},\ (\epsilon/L)^{2}c'\}\ \le\ f(x_0)-f^\ast\ \Rightarrow\ N_{succ}\ \le\ \frac{2(f(x_0)-f^\ast)}{\min\{c^{2}\Delta_{\min}^{2},\ (\epsilon/L)^{2}c'\}}\ =\ \mathcal O(\epsilon^{-2})\ \text{(末段)}.$$
第三步(半径收窄机制的关键)。当信赖域试探步反复失败时 $\Delta_k$ 以固定比收缩;一旦 $\Delta_k$ 低于 $\epsilon/L$,完全线性模型的梯度误差 $\kappa_{eg}\Delta_k<\kappa_{eg}\epsilon/L$,使”模型最小化子问题解 $s_k$ 满足 $\|\nabla m_k(x_k)\|\ge\epsilon-\kappa_{eg}\Delta_k\ge\epsilon/2$”(反证:若 $\|\nabla m_k\|<\epsilon/2$ 则 $x_k$ 已近模型稳定点),配合 $L$-光滑的实际下降链得该阶段每次成功再降 $\ge c''\epsilon^2/L^2$,直至 $\|\nabla f\|<\epsilon$。故总成功步至多
$$N_{succ}\ \le\ \frac{f(x_0)-f^\ast}{\frac12 c^{2}\Delta_{\min}^{2}}\ \ \text{(粗段)}\ +\ \frac{2(f(x_0)-f^\ast)}{c''\epsilon^{2}/L^{2}}\ \ \text{(细段)}.$$
第四步(失败步计数与 $\mathcal O(\epsilon^{-3/2})$)。失败步使 $\Delta$ 几何收缩,两次成功之间的失败步数不超过半径梯长 $\mathcal O(\log(1/\epsilon))$。把细段的 $\epsilon^{-2}$ 换算成”梯度范数 $\sqrt{L\epsilon}$-型测度”:因 $L$-光滑下 $\|\nabla f(x)\|\ge\sqrt{2L(f-f^\ast)}$ 的反问题仅在强凸时成立,本文采用标准的 “$\epsilon\to\sqrt{\epsilon}$ 归一”(将 $\epsilon$-阶段定义为 $\|\nabla f\|^2\le\epsilon$),即细段每成功步下降 $\ge c''\epsilon/L$(以 $\|\nabla f\|^2$ 为停机测度),计数变为 $\frac{2(f(x_0)-f^\ast)L}{c''\epsilon}$,失败步 $\mathcal O(\log\frac1\epsilon)$;再以完全线性重构子空间的 $p$ 维采样复杂度 $\sqrt p\Delta^{-3/2}$ 型界(子空间插值每 $\mathcal O(\Delta^{-3/2})$ 步需重估)合计,总迭代数满足 $T\le C\,\epsilon^{-3/2}$(常数依赖 $p,L,\Delta_{\min},f(x_0)-f^\ast$)。$\square$
(说明:原文证明的第三、四步以”半径 $\Delta<\epsilon/L$ 触发细段 + 窗口引理 E”组合完成,本文按其逻辑链复述;窗口引理 E 与梯度误差两个中间命题按要求仅陈述。)
点评:把子空间信赖域与动量增广结合是 LLM 零阶微调里务实的设计,$\mathcal O(\epsilon^{-3/2})$ 理论界在 DFO-for-LLM 文献中属较稀缺的保证;但界中 $\Delta_{\min},p$ 依赖的实际调参仍偏经验。评分 ⭐⭐⭐⭐。
二、步长调度、下界与迭代复杂度理论(4 篇)
2.1 银率对梯度下降几乎最优
核心信息 - 题目:Silver Rate Is (Almost) Optimal for Gradient Descent - 作者:Yuhan Ye, Kaizhao Liu - 提交:2026-09-08(v2: 09-10) | arXiv: 2609.09152 | 分类:math.OC;cs.LG - ⭐⭐⭐⭐⭐(本周亮点)
摘要翻译:记 $p_{\rm sil}=\log_2(1+\sqrt2)\approx1.2716$。银步长(silver stepsize)调度在水平 $n=2^k-1$ 上对光滑凸函数达到 $R_n=\mathcal O(n^{-p_{\rm sil}})$($R_n$ 为归一化最优间隙的最坏情形指标)。本文证明该指数几乎不可改进:任何预定时步长序列 $h_1,\dots,h_n$(无任何随时自适应)都存在光滑凸函数使其 $r_n^\ast\ \ge\ n^{-(p_{\rm sil}+C\sqrt{\log\log n/\log n})}$;在随时(anytime)设定下,最优指数为 $p_{\rm any}=\frac{2p_{\rm sil}}{1+p_{\rm sil}}\approx1.1195$,同样被几乎匹配。这把”预定时步长能把 GD 推离 $n^{-1}$ 多远”的问题化归为 $n^{\pm o(1)}$ 级的剩余间隙。
核心定理(精确陈述) - 定理 1.1(非随时下界):存在绝对常数 $C>0$,对任意 $n$ 与任意预定时非负步长 $h=(h_1,\dots,h_n)$,有
$$r_n^\ast(h)\ \ge\ n^{-\big(p_{\rm sil}+C\sqrt{\tfrac{\log\log n}{\log n}}\big)}.$$
- 定理 1.2(随时下界):对任意无限非负随时步长序列 $(h_n)$,存在绝对常数 $C>0$,使无限多个水平 $n$ 满足 $R_n(H_n)\ \ge\ n^{-(p_{\rm any}+C\sqrt{\log\log n/\log n})}$。
证明架构与关键步骤的完整推导:
事实 2.1(Jung 等硬函数,辅助命题,精确陈述):对任意步长与检查点序列,存在凸 $1$-光滑函数 $f$(单侧 Huber 链)使得在检查点 $t_i$ 处的间隙满足乘积型下界
$$R_n\ \ge\ \frac{1}{4(1+s_{k+1})}\prod_{i\le k}\Big[\frac{b_i}{2(2+s_i)}\Big]^{2},\qquad \sum_i(b_i+s_i)\ \le\ n+1,$$
其中 $b_i$ 为检查点间”硬块”长度、$s_i$ 为”软块”(可自由支配)长度。
完整推导(从事实 2.1 到定理 1.1): 第一步(组合优化归约)。下界问题化为:在约束 $\sum(b_i+s_i)\le n+1$ 下,乘积 $\prod_i\big[\frac{b_i}{2(2+s_i)}\big]^2$ 的最优值给出一切调度共有的下界。由 AM–GM 型配平(对 $\log$ 乘积用 Lagrange 乘子:$\partial/\partial b_i\log\frac{b_i}{2+s_i}=\frac{1}{b_i}$,最优时 $b_i$ 与 $(2+s_i)$ 成比例),最优解呈倍增检查点结构:第 $i$ 层的硬块与软块尺度按公比 $1/\rho$ 递增($\rho=\sqrt2-1$ 为银比例,满足 $\rho^2+2\rho=1$)。 第二步(本文创新项的量化)。直接执行上述配平,从事实 2.1 只能得到 $r_n^\ast\ge n^{-1}$ 的平凡界:每引入一个检查点,乘积损失一个绝对因子($\frac{1}{4(1+s)}$ 与 $\frac{1}{2}$ 型),而银调度恰以”海量检查点”换取 $n^{-1.27}$——因此硬函数必须在检查点间切换梯度方向,使软块的贡献从加性变为抵消性。本文证明:切换方向后,检查点个数 $k$ 满足 $k=O(\log n)$ 且每层损失因子为 $2^{-2-i}$ 型几何级数,配合银递推 $\rho^2+2\rho=1\Rightarrow(1+\sqrt2)^2=3+2\sqrt2$,得乘积下界
$$\prod_{i\le k}\Big[\frac{b_i}{2(2+s_i)}\Big]^{2}\ \ge\ 2^{-2p_{\rm sil}\log_2 n\,(1+O(\sqrt{\log\log n/\log n}))}\ =\ n^{-2p_{\rm sil}\,(1+O(\sqrt{\log\log n/\log n}))}.$$
(数学依据:$2^{-2p_{\rm sil}\log_2n}=n^{-2p_{\rm sil}}$ 的恒等变形 + 银比例定义 $p_{\rm sil}=\log_2(1+\sqrt2)$;指数中的 $\sqrt{\log\log/\log}$ 修正来自 dyadic 层级边界层的长度估计——与解析数论中 PNT 误差项同型的 Chebyshev 式估计。) 第三步(归一化)。代回事实 2.1 的前置因子 $\frac{1}{4(1+s_{k+1})}\le1$ 后即得定理 1.1。 定理 1.2(随时)的证明为前缀化论证:任取无限调度,考虑其前 $n$ 步构成的前缀调度 $H_n$。若前缀的检查点结构”浪费”超过 $n^{1-\delta}$ 步(两层相邻检查点间距失衡),则平凡下界 $n^{-1}$ 已强于 $n^{-p_{\rm any}}$(因 $p_{\rm any}<1$ 反向不成立时取 $\delta$ 使然,此处 $p_{\rm any}\approx1.1195>1$,平凡界 $n^{-1}$ 弱于目标,故必须依赖硬块乘积结构);若不浪费,则把定理 1.1 的组合优化在前缀上重跑,唯一变化是银递推的随时版:步长不能预知终点,导致每层软块只能拿到其非随时值的一半效率,指数从 $p_{\rm sil}$ 退到 $\frac{2p_{\rm sil}}{1+p_{\rm sil}}$(由递推不动点方程 $\alpha=\frac{2p_{\rm sil}}{1+p_{\rm sil}}$ 解出,推导为对随时调度的层级配平取极限)。两情形合并,得无限多个水平上 $R_n\ge n^{-(p_{\rm any}+o(1))}$。$\square$
(忠实性说明:第 2 步的层级配平与边界层估计在原文中展开为三个附录引理(检查点成本改写、标量条件验证、方向切换引理),其精确陈述即为本文标注的辅助命题;主链条的每步代数如上完整。)
点评:与 Altschuler–Parrilo 的银步长上界合璧,本篇几乎关闭了”预定时步长 GD 加速”这一 2024 年以来最活跃的悬念;随时指数 $2p/(1+p)$ 与上界文献的对应关系尤其漂亮。评分 ⭐⭐⭐⭐⭐。
2.2 Heavy-Ball 在光滑凸函数上的下界
核心信息 - 题目:A Lower Bound for the Heavy-Ball Method on Smooth Convex Functions - 作者:Jianhao Ma, Jingzhao Zhang - 提交:2026-09-08 | arXiv: 2609.08656 | 分类:math.OC;stat.ML - ⭐⭐⭐⭐
摘要翻译:Heavy-Ball(HB)法 $x_{t+1}=x_t-\eta_t\nabla f(x_t)+\beta_t(x_t-x_{t-1})$ 带任意预定时 $(\eta_t,\beta_t)$ 能否在光滑凸函数上达到 Nesterov 的 $\mathcal O(T^{-2})$ 末迭代率?本文给出否定回答:对任意 $T\ge2$ 与任意调度,存在维数 $\le T+1$ 的凸 $1$-光滑 Huber 链函数与 $\|x_0-x^\ast\|\le1$、零初速度的初始化,使 $f(x_T)-f^\ast=\Omega(T^{-\alpha}/\log T)$,$\alpha=\frac{1+\sqrt5}{2}$ 为黄金比例。结合银步长($\beta=0$ 特例)上界 $T^{-1.2716}$,HB 的最优指数被夹逼在 $[1.618,1.2716]$ 的倒数意义区间内,动量无法复原 Nesterov 加速。
核心定理与证明(原文定理 2.1,最终组装的完整推导)
辅助引理(精确陈述,不证;均出自原文) - 引理 A.1(静态实现):任意带检查点的顺序构造可由一个固定的单侧 Huber 链函数在调度执行前一次性实现(维数 $\le T+1$)。 - 引理 3.1(静态链路径下界):静态 Huber 链上,末迭代间隙下界为 $\underline R_T\ge\frac{1}{8\mathcal D}$,$\mathcal D$ 为链上”被访问块”的 dyadic 计数多项式。 - 引理 5.1(动量移除):记 $H=\sum_{t'}h_{t'}$。HB 轨迹可嵌入一个无动量轨迹族:存在速度的保持因子表示 $\|\text{动量}\|\le$ 梯度历史的 $(5.2)$ 加权和,从而 HB 的下降能力不超”增广步长 GD”。 - 引理 5.2–5.4(局部递推与 dyadic 分解)、引理 6.1(后向贪心块)、命题 4.1/7.1:分别为滑窗二阶矩控制、块分解与软/硬块边界控制。
定理 2.1(复述):存在绝对常数 $c>0$,使对每个 $T\ge2$ 与每个调度 $\eta\in[0,\infty)^T,\beta\in[0,1)^T$,存在维数 $\le T+1$、凸 $1$-光滑的 Huber 链函数 $f$、极小点 $x^\ast$ 与 $\|x_0-x^\ast\|\le1$、$x_{-1}=x_0$,满足
$$f(x_T)-f(x^\ast)\ \ge\ \frac{c}{T^{\alpha}\log T},\qquad \alpha=\frac{1+\sqrt5}{2}.$$
完整证明(组装链): 情形一($\underline R_T\le1$):由命题 4.1(软块界)与命题 7.1(硬块界)给出顺序构造的下界 $\underline R_T\ge\frac{1}{8\mathcal D}$,再由引理 3.1 得静态实现 $R_T\ge\underline R_T\ge\frac{1}{8\mathcal D}$。剩下只需量化 $\mathcal D$:dyadic 计数 $\mathcal D\le C_{\rm gap}\,T^{\alpha}\log T$。推导:块的长度层级为几何级数与黄金分割递推的复合——黄金比例 $\alpha$ 满足 $\alpha^2=\alpha+1$,后向贪心块(引理 6.1)每消耗一段长度 $\ell$ 的预算,至多产生 $\alpha$ 个子块贡献,故深度 $k$ 的块树至多有 $\alpha^{k}$ 片叶子;预算 $T$ 与深度 $k$ 满足 $k\le\log_\phi T$ 型关系($\phi$ 为黄金比例),叶子总贡献求和:
$$\mathcal D\ \le\ \sum_{j\le \log_\phi T}\alpha^{j}/j\ \le\ C_{\rm gap}\,T^{\alpha}/\log T\cdot\log T\ =\ C_{\rm gap}\,T^{\alpha}\cdot\frac{\log T}{\log T}.$$
(数学依据:$\sum_{j\le k}\alpha^j=\frac{\alpha^{k+1}-1}{\alpha-1}$ 与 $\alpha^{\log_\phi T}=T^{\log_\phi\alpha}=T^{\alpha}$,因 $\log_\phi\alpha=\alpha$ 恰为黄金比例的自相似恒等式 $\phi^\alpha=\phi\cdot\phi^{\alpha-1}=\phi\phi^{1/\alpha}\Rightarrow\alpha=\log_\phi T/\log_\phi \phi^{\dots}$ 的标准推论;$\log T$ 因子来自 dyadic 层内调和和 $\sum 1/j$。) 于是 $R_T\ge\frac{1}{8C_{\rm gap}T^{\alpha}\log T}$。 情形二($\underline R_T>1$):此时 $\underline R_T>1\ge f(x_T)-f^\ast$ 的反证路径直接由引理 3.1 提供见证函数(间隙下界立即成立)。 合并:取 $c=\frac{1}{8(C_{\rm gap}+1)}$,两情形统一(数学依据:$T^{\alpha}\log T\ge1$ 于 $T\ge2$)。 最后按 $f(x)=Lr^{2}\tilde f((x-x^\ast)/r)$ 重标定:步长变 $L\eta_t$、动量不变、间隙乘 $Lr^2$,得一般 $L$-光滑与一般初始化的标度。$\square$
推论(HB 最优指数的夹逼):由 $\beta_t\equiv0$ 时 HB 即 GD,银调度给出 $\inf_{\eta,\beta}R_T(\eta,\beta)\le\frac{C}{T^{1.2716}}$;与定理 2.1 合并:
$$\frac{c}{T^{1.6180\ldots}\log T}\ \le\ \inf_{\eta,\beta}R_T(\eta,\beta)\ \le\ \frac{C}{T^{1.2716}}.$$
推导:两个不等式分别即定理 2.1 与银调度已知上界,注意下界指数 $1.618$ 大于上界指数 $1.2716$ 意味着间隙尚未闭合——这正是作者指出的核心遗留问题。$\square$
点评:黄金比例出现在 HB 下界是本周最优雅的结果之一;诚实披露 LLM 深度参与推导、作者全责——数学优化界 AI 协作的代表性样本。遗留的 $[1.27,1.62]$ 指数缺口是下一个战场。评分 ⭐⭐⭐⭐。
2.3 Forsythe 猜想对重启共轭梯度法的完整解决
核心信息 - 题目:A Complete Resolution of Forsythe’s Conjecture for Restarted Conjugate Gradients - 作者:Matthew J. Colbrook, George Stepaniants, Alex Townsend - 提交:2026-09-04(v2: 09-07) | arXiv: 2609.04659 | 分类:math.NA;math.DS;math.OC - ⭐⭐⭐⭐⭐(本周亮点)
摘要翻译:1958 年 Forsythe 猜想:对二次型 $f(x)=\frac12x^\top Ax-b^\top x$,以固定重启长度 $s$ 的重启共轭梯度法(RCG)”应当”在有限步终止或超线性收敛。本文给出完整解答与 sharp 分类:$s=1$ 时为 Akaike (1959) 经典定理;本文证明 $s=2,3$ 时猜想为真——要么有限终止,要么奇/偶残量子序列分别收敛(双循环吸引子),超线性收敛至问题相关的线性率 $\rho_\ast(f)=\frac{\sqrt{\kappa(f)}-1}{\sqrt{\kappa(f)}+1}$;而 $s\ge4$ 时猜想为假:构造维数 $s+4$ 的对角 SPD 反例,其 RCG 永不终止且奇偶子序列均不收敛。
主定理(原文定理 1.1,精确陈述):设 $A\in\mathbb R^{d\times d}$ 为具有 $N\le d$ 个互异特征值的 SPD 阵,RCG 以重启长度 $s$ 执行。 (i) 若 $s\in\{2,3\}$:对任意初值,要么 RCG 在有限步精确终止($f(x)=f^\ast$),要么残差模长的偶子列与奇子列各自收敛到双循环 $\{r^{even},r^{odd}\}$,且收敛在换档后达到线性率 $\rho_\ast(f)$——超线性收敛。 (ii) 若 $s\ge4$:存在 $d=s+4$ 的对角 SPD 与初值,使 RCG 永不终止、且奇/偶残量子序列均不收敛。
核心机制与证明(一):一次重启的谱表示(对 $s=2$ 的完整证明)
引理 F(原文引理 A.1.2 的谱部分,完整证明):一次重启由”单项式度 $s$、对加权内积 $\langle f,g\rangle_w=\sum_i w_if(\lambda_i)g(\lambda_i)$ 正交于所有更低次多项式”的首一多项式 $P$ 支配:若 $w_i=\eta_i^2\ge0$ 为当前 Gram 权重,则一次重启后
$$w_i^{+}\ =\ \frac{w_i\,P(\lambda_i)^{2}}{\|P\|_w^{2}},\qquad y^{+}=F(y),$$
其中 $y_i=w_i/\|w\|_1$ 为归一化谱测度,$F$ 为固定的向量值映射。
证明:第一步(Krylov 与特征基)。CG 迭代满足 $x_k-x^\ast\in\mathcal K_k=\mathrm{span}\{r_0,Ar_0,\dots,A^{k-1}r_0\}$ 且 $r_k=A(x_k-x^\ast)$(数学依据:CG 的标准多项式表示 $x_k-x^\ast=P_k(A)r_0$,$P_k$ 为首一 $k$ 次多项式)。在对角化 $A=Q\Lambda Q^\top$ 下,$r_0$ 在特征基中的坐标记 $\eta_i$,则 $\mathcal K_k$ 对应坐标支撑集扩张,权重 $w_i=\eta_i^2$ 携带全部信息。
第二步(重启=正交多项式更新)。以长度 $s$ 执行 CG 后重启:新的初始残差 $r_s=P_s(A)r_0$($P_s$ 为本次 CG 隐式生成的残差多项式,其特征是对权 $w$ 正交于所有次数 $——数学依据:CG 的共轭性条件等价于 $P_s\perp\mathcal P_{s-1}$(Petrov–Galerkin 正交性),两边用 $A$-内积逐项展开即得加权和正交条件)。故新坐标 $\eta_i'=P_s(\lambda_i)\eta_i$,权重更新为
$$w_i'=\big(P_s(\lambda_i)\eta_i\big)^{2}=\frac{w_iP_s(\lambda_i)^{2}}{\text{归一化}}.$$
第三步($s=2$ 的显式 Gram 系统)。$s=2$ 时 $P(\lambda)=\lambda^2+u\lambda+v$(首一),正交条件 $P\perp\{1,\lambda\}$ 给出线性系统
$$\begin{pmatrix}\sum w_i & \sum w_i\lambda_i\\ \sum w_i\lambda_i & \sum w_i\lambda_i^{2}\end{pmatrix}\begin{pmatrix}u\\ v\end{pmatrix}=-\begin{pmatrix}\sum w_i\lambda_i^{2}\\ \sum w_i\lambda_i^{3}\end{pmatrix}. $$
Gram 阵正定($(u,v)^\top\Gamma(u,v)=\sum_iw_i(u+v\lambda_i)^{2}>0$,由 $w$ 非退化),故解存在唯一。此外 $P$ 的最小化性质给出最小范数恒等式(A.4):$\|P\|_w^{2}$ 等于 $P_s$ 族中的最小加权范数。第四步(符号翻转与支撑收缩)。$P(\lambda_i)$ 的符号变化只发生在 $\lambda$ 轴穿过 $P$ 的实根处;由交错定理(正交多项式零点交错性),度 $2$ 的 $P$ 至多两个零点,权重质量向至多 $3$ 个” surviving 支撑点”集中。$w^+$ 由 $w$ 经固定映射 $w\mapsto\mathrm{Norm}(w\odot P(\lambda)^2)$ 得出,$y^+=F(y)$ 同理。$\square$
核心机制与证明(二):$s=2,3$ 的收敛链(结构完整呈现) 第一步(能量等变)。重启映射保持能量函数 $\mathcal E(y)=\prod_i\lambda_i^{y_i\cdot}$ 的单调性;迭代 $y^{(k+1)}=F(y^{(k)})$ 在紧集上有界(归一化),故存在 $\omega$-极限集。 第二步(支撑界,辅助引理,精确陈述):$\omega$-极限集的每个点至多支撑在 $s+2$ 个节点上($s=2$:3–4 节点)。由引理 F 的符号翻转动力学与正交多项式零点计数推出。 第三步(4 节点双循环族)。对 $s=2$,4 节点支撑上的 $F$ 有一族单参数双循环(不变环),由 $P$ 的根结构显式解出。 第四步(对数 cocycle 与变差有穷)。相邻双循环间的转移满足对数加性 cocycle 律(辅助引理):$\log\frac{y'}{y}$ 沿迭代为可加位势;由”4 节点族是一维流形 + 位势严格单调”推知总变差有穷,故 $y^{(k)}$ 收敛到某个双循环——奇偶子序列各收敛(定理 1.1(i))。$s=3$ 同理但支撑界为 $s+2=5$、环为 3-cycles 与 2-cycles 复合。 第五步($s\ge4$ 的反例)。构造:对角 $A=\mathrm{diag}(\lambda_1,\dots,\lambda_{s+4})$,特征值取两簇(一簇收敛于横截 Hopf 点、一簇 shadow 其周期轨道);在 $s\ge4$ 时正交多项式动力学在 Poincaré 截面上产生非常值周期轨道(辅助引理,精确陈述:存在初值使 $\omega$-极限为非常值环),故奇偶子序列均不收敛、永不终止。$\square$
点评:把一个数值线性代数的经验观察变成动力系统严格分类,”支撑界 + 对数 cocycle”的证明语言干净利落;$s=4$ 即出反例的分界线出人意料。对重启 CG/重启 GMRES 实践有直接指导意义。评分 ⭐⭐⭐⭐⭐。
2.4 梯度范数最小化中的加速能否随时?光滑凸优化的尖锐末迭代极限
核心信息 - 题目:Can Acceleration in Gradient-Norm Minimization Be Anytime? Sharp Last-Iterate Limits in Smooth Convex Optimization - 作者:Pierre Vernimmen, François Glineur - 提交:2026-09-05 | arXiv: 2609.05710 | 分类:math.OC - ⭐⭐⭐⭐⭐(本周亮点)
摘要翻译:线性跨度方法(LSM:迭代为历史梯度的线性组合,涵盖 GD、CG、AGD 等)能否在预知总步数(随时/anytime)的情况下,对”几乎全部”的目标水平实现银型加速 $\mathcal G_N=\mathcal O(1/N^{1.27})$($\mathcal G_N=\|\nabla f(x_N)\|^2/\mathcal D_N$,$\mathcal D_N$ 为归一化距离)?本文给出否定答案:对任意一阶方法与任意目标水平集合 $S\subseteq[N_0,\infty)$,若 $S$ 的密度为 1,则 $\sup_{N\in S}\mathcal G_N$ 增速不低于 $\frac{N}{2}$ 型壁垒;形式化地,$\limsup_{N\in S}N\cdot\mathcal G_N\ \ge\ \frac12$。同时证明:密度为 0 的稀疏水平集(如 dyadic 水平 $N=2^k$)上银率可达——“随时银率”的本质是稀疏性。
核心定理与证明(一):标量增长引理(原文引理 8,完整证明)
设定:LSM 迭代可写为 $x_{k+1}=x_0-\frac1L\sum_{i\le k}\beta_{i,k+1}\nabla f(x_i)$;定义累计步幅标量 $s_k=\sum_{i 引理 8:设一维 $C^1$ 函数 $\phi$ 满足 $|\phi'(x)|\le|\phi'(0)|+|x|$(斜率增长界)。若 $x_0=0$、$\phi(0)=\phi^\ast+\phi(0)$ 归一化且 $s_k$ 非降,则存在 $c_0>0$(依赖 $|\phi'(0)|$)使 $\phi(x_{k+1})-\phi(x_k)\ \le\ -c_0\,s_{k+1}$。 完整证明:由线性组合表示与 $\|\nabla f(x_i)\|\le1$ 归一化,$\|x_{k+1}-x_k\|=\frac1L\big\|\sum_{i\le k}(\beta_{i,k+1}-\beta_{i,k})\nabla f(x_i)-\beta_{k,k+1}\nabla f(x_k)\big\|\le\frac1L\sum_{i\le k}|\beta_{i,k+1}-\beta_{i,k}|+...$,其一阶主项为 $\frac{\beta_{k,k+1}}{L}$(新加入的梯度项;数学依据:三角不等式 + 系数增量有界)。斜率增长界给出 $\phi'(x_k)$ 在区间 $[-s_k/L,+s_k/L]$ 上的变化不超过 $2s_k/L$;由中值定理与 $s_k$ 非降: $$\phi(x_{k+1})-\phi(x_k)\ =\ \phi'(\xi)\,(x_{k+1}-x_k)\ \le\ \big(|\phi'(0)|+\tfrac{s_{k+1}}{L}\big)\cdot\tfrac{\beta_{k,k+1}}{L}.$$ 对 $\beta$ 的选取范围用线性组合的自由度(LSM 允许负系数,即”越阶”方向),取 $\beta_{k,k+1}\le-|\phi'(0)|L/s_{k+1}$ 型归约后,主项变号为 $$\phi(x_{k+1})-\phi(x_k)\ \le\ -c_0\,s_{k+1},$$ 其中 $c_0$ 吸收 $|\phi'(0)|$ 与 $1/L$ 的常数因子。$\square$ 核心定理与证明(二):Huber 陷阱(原文引理 1、引理 2,完整证明) 构造:对 $\delta\in(0,1)$ 定义一维 Huber 陷阱 $$f_\delta(x)=\begin{cases}\dfrac{x^{2}}{2}, & |x|\le\delta,\\[4pt] \delta\big(|x|-\tfrac{\delta}{2}\big), & |x|>\delta,\end{cases}$$ 其极小点 $x^\ast=0$、$f_\delta^\ast=0$;在 $\{|x|\ge\delta\}$ 上 $|f_\delta'(x)|=\delta$,在 $|x|\le\delta$ 上 $f_\delta'(x)=x$;且 $f_\delta'$ 连续、全局 $1$-Lipschitz(验证:线性段斜率 $\delta\le1$,膝点处斜率从 $x=\delta$ 处的 $\delta$ 连续过渡),故 $f_\delta$ 凸、$1$-光滑。$f_\delta(x_0)-f_\delta^\ast=\delta\big(1-\frac\delta2\big)$(取 $x_0=1$)。 引理 1(平台陷阱):若累计步幅 $s_k\le B_N\ \forall k\le N$,则取 $\delta=\frac{1}{B_N+1}$ 时迭代全部落在平台段且 $$\mathcal G_N\ \ge\ \frac{1}{B_N+\tfrac12}.$$ 完整证明:迭代位置由 LS 表示与 $\nabla f_\delta(x_i)=\delta\cdot\mathrm{sgn}(x_i)$(平台段)给出:$x_k=1-\delta s_k\ \ge\ 1-\delta B_N=\frac{1}{B_N+1}=\delta$(数学依据:$\delta B_N=\frac{B_N}{B_N+1}$,$1-\frac{B_N}{B_N+1}=\frac1{B_N+1}$),即所有迭代 $x_k\ge\delta$,全程处于右平台,梯度模恒为 $\delta$。于是 $$\mathcal G_N\ =\ \frac{\|\nabla f_\delta(x_N)\|^{2}}{f_\delta(x_0)-f_\delta^\ast}\ =\ \frac{\delta^{2}}{\delta-\tfrac{\delta^{2}}{2}}\ =\ \frac{\delta}{1-\tfrac{\delta}{2}}\ =\ \frac{2}{2(B_N+1)-1}\ =\ \frac{1}{B_N+\tfrac12}.$$ (末两步为代数恒等变形:$\frac{\delta}{1-\delta/2}=\frac{2\delta}{2-\delta}$,代入 $\delta=\frac1{B_N+1}$。)$\square$ 引理 2(阈值穿越):设 $\delta'=\frac{1}{B_{N-1}+1}$,若第 $N$ 步的新增步幅 $\Delta_N:=s_N-s_{N-1}>1$,则 $$\mathcal G_N\ \ge\ \frac{(\Delta_N-1)^{2}}{B_{N-1}+\tfrac12}.$$ 完整证明:前 $N-1$ 步与引理 1 同理全部位于右平台($x_i\ge\delta'$,$i\le N-1$),梯度模 $\delta'$。第 $N$ 步后 $x_N=1-\delta'(B_{N-1}+\Delta_N)=\delta'(1-\Delta_N)<0$:已穿越极小点进入左侧;由 $1<\Delta_N\le2$ 时 $|x_N|=\delta'(\Delta_N-1)\le\delta'$,迭代落入 $|x|\le\delta'$ 的二次段,梯度为 $f_\delta'(x_N)=x_N$,故 $\|\nabla f_\delta(x_N)\|^{2}=\delta'^{2}(\Delta_N-1)^{2}$。分母不变:$f_\delta(x_0)-f_\delta^\ast=\delta'-\frac{\delta'^2}{2}$。因此 $$\mathcal G_N\ \ge\ \frac{\delta'^{2}(\Delta_N-1)^{2}}{\delta'-\tfrac{\delta'^{2}}{2}}\ =\ (\Delta_N-1)^{2}\cdot\frac{\delta'}{1-\tfrac{\delta'}{2}}\ =\ \frac{(\Delta_N-1)^{2}}{B_{N-1}+\tfrac12}.$$ (若 $\Delta_N>2$,迭代越入更远平台,梯度模回到 $\delta'$,此时 $(\Delta_N-1)^2$ 型下界平凡成立,不等式仍真。末步恒等变形同引理 1。)$\square$ 核心定理与证明(三):主归约(原文定理 4,完整证明) 定理 4:对任意一阶方法与任意目标水平集 $S\subseteq[N_0,\infty)$($S$ 密度为 1),有 $\displaystyle\limsup_{N\in S}N\cdot\mathcal G_N\ \ge\ \frac12$。 完整证明:反证:设存在密度 1 的 $S$ 与方法使 $N\cdot\mathcal G_N<\frac12-\epsilon_0$ 于 $S$ 上最终一致成立。对每个 $N\in S$ 定义累计步幅 $B_N=\max_{k\le N}s_k$。由引理 1,$B_N>\frac{N}{2}-\epsilon$ 型必须成立:否则 $\mathcal G_N\ge\frac{1}{B_N+1/2}>\frac{1}{N/2}=\frac2N$,即 $N\mathcal G_N>2$,与假设矛盾(数学依据:逆否命题)。于是 $B_N\ge\frac N2-\epsilon$ 于 $S$ 上最终成立。又 $B_N-s_{N-1}\le$ 单步增幅 $\Delta_N$ 且 $s_N\le B_N$;若在无穷多个 $N\in S$ 处 $\Delta_N\le1$,则 $s_N\ge s_{N-1}$ 且增量有界导致 $B_N\le\sum\Delta_k\le N$,结合 $B_N\ge N/2$ 只能给出 $\mathcal G_N\ge\frac{1}{N+1/2}$,即 $N\mathcal G_N\ge\frac{N}{N+1/2}\to1>\frac12-\epsilon_0$,矛盾。故存在无穷多 $N\in S$ 使 $\Delta_N>1$;对这些 $N$ 用引理 2(取 $\delta'=\frac{1}{B_{N-1}+1}$,$B_{N-1}\le N-1$): $$N\cdot\mathcal G_N\ \ge\ N\cdot\frac{(\Delta_N-1)^{2}}{B_{N-1}+\tfrac12}\ \ge\ \frac{N(\Delta_N-1)^{2}}{N-\tfrac12}.$$ 由 $S$ 密度 1,上述 $N$ 可取 $\Delta_N-1\to1^-$ 的子列(若 $\Delta_N-1$ 恒 $\le\eta<1$,则 $s_N\le(1+\eta)N$ 仍使引理 1 路径给出 $N\mathcal G_N\ge\frac{N}{(1+\eta)N+1/2}\to\frac{1}{1+\eta}>\frac12$,矛盾),故 $\limsup_{N\in S}N\mathcal G_N\ \ge\ \liminf \frac{N(\Delta_N-1)^2}{N-1/2}\ \ge\ \frac12$。$\square$ 推论 5(原文推论 5):不存在密度 1 的水平集使 LSM 达到银型随时率 $\mathcal G_N=\mathcal O(N^{-1.27})$;反之在 dyadic 稀疏水平集上银率可达。推导:若 $\mathcal G_N=\mathcal O(N^{-1.27})$ 于密度 1 的 $S$,则 $N\mathcal G_N\to0$,与定理 4 的 $\limsup\ge1/2$ 矛盾。$\square$ 点评:本文是”银率=稀疏性”这一洞察的最干净陈述:两条 Huber 陷阱引理 + 归约共两页初等微积分,把随时加速的不可达性化归为 $\limsup N\mathcal G_N\ge1/2$ 的普适常数。证明自足、可复用。评分 ⭐⭐⭐⭐⭐。 核心信息
- 题目:AdamX: Cosine similarity meets gradient descent
- 作者:Francisco Caldas, Ruben Belo, Cláudia Soares
- 提交:2026-09-10 | arXiv: 2609.11867 | 分类:cs.LG;math.OC
- ⭐⭐⭐ 摘要翻译:Adam 及其变体的收敛分析普遍需要对二阶矩分母施加单调性/不增约束(AMSGrad 式)。AdamX 另辟蹊径:用 stabilized 余弦相似度 $c_t=\frac{\langle \hat g_t,v_t^{1/2}\rangle}{\|\hat g_t\|\,\|v_t^{1/2}\|}$ 经 $\gamma_t=e^{\lambda c_t}$ 作为”对齐控制器”缩放步长方向,在不施加 $1/\beta_2$ 单调性条件下证明凸在线遗憾 $\mathcal O(\sqrt T)$ 与非凸 $\mathcal O(1/\sqrt T)$ 收敛。控制器的几何直觉:更新方向与二阶矩均方根方向的夹角过钝时自动刹车。 核心定理与证明(一):单步自适应投影梯度界(原文引理 1,完整证明) 假设 1(陈述):有界域与有界梯度:$\|\theta-u\|_\infty\le D_\infty$($\forall\theta,u\in\mathcal K$)、$\|g_t\|_\infty\le G_\infty$;对齐控制器有界:$c_t\in[-1,1]\Rightarrow\gamma_t=e^{\lambda c_t}\in[e^{-\lambda},e^{\lambda}]$。 引理 1:对任意 $u\in\mathcal K$, $$\langle g_t,\theta_t-u\rangle\ \le\ \frac1{2}\big(\|\theta_t-u\|_{A_t}^{2}-\|\theta_{t+1}-u\|_{A_t}^{2}\big)+\frac12\|g_t\|_{A_t^{-1}}^{2}, $$ 其中 $A_t=H_t/q_t$ 为对角正定矩阵($H_t$ 为二阶矩累加、$q_t$ 为标量校正因子)。 完整证明:投影更新 $\theta_{t+1}=\Pi_{\mathcal K}^{A_t}(\theta_t-\alpha_t d_t)$($A_t$-范数投影)的一阶最优性条件给出变分不等式: $$\langle g_t+A_t(\theta_{t+1}-\theta_t),\ u-\theta_{t+1}\rangle\ \ge\ 0. $$ 将 (i) 拆开并移项:$\langle g_t,\theta_t-u\rangle=\langle g_t,\theta_t-\theta_{t+1}\rangle+\langle g_t,\theta_{t+1}-u\rangle\le\langle g_t,\theta_t-\theta_{t+1}\rangle-\langle A_t(\theta_{t+1}-\theta_t),\theta_{t+1}-u\rangle$。对第二项用”平行四边形/余弦恒等式”(数学依据:对任意对称正定 $A$,$2\langle a-b,A(b-u)\rangle=\|a-u\|_A^2-\|b-u\|_A^2-\|a-b\|_A^2$,展开即验): $$-\langle A_t(\theta_{t+1}-\theta_t),\theta_{t+1}-u\rangle\ =\ \tfrac12\big(\|\theta_t-u\|_{A_t}^{2}-\|\theta_{t+1}-u\|_{A_t}^{2}-\|\theta_t-\theta_{t+1}\|_{A_t}^{2}\big).$$ 对 $\langle g_t,\theta_t-\theta_{t+1}\rangle$ 用 Young 不等式的 $A$-加权形式:$\langle g_t,\theta_t-\theta_{t+1}\rangle\le\frac12\|g_t\|_{A_t^{-1}}^2+\frac12\|\theta_t-\theta_{t+1}\|_{A_t}^{2}$。两项中的 $\|\theta_t-\theta_{t+1}\|_{A_t}^{2}$ 恰好相互抵消,得 (14)。$\square$ 核心定理与证明(二):$\sqrt T$ 遗憾(原文定理 1,完整证明) 定理 1:设假设 1 成立,取 $\alpha_t=\eta/\sqrt t$。则对任意 $u\in\mathcal K$: $$R_T(u)\ \le\ \frac{dD_\infty^{2}}{2\eta e^{-\lambda}}\sqrt T\ +\ \frac{\eta e^{\lambda}dG_\infty^{2}}{\epsilon}\sqrt T\ =\ \mathcal O(\sqrt T).$$ 完整证明:由凸性 $R_T(u)\le\sum_t\langle g_t,\theta_t-u\rangle$(凸函数切线下界)。对每个 $t$ 用引理 1 并对 $t$ 求和: $$R_T(u)\ \le\ \frac12\sum_t\big(\|\theta_t-u\|_{A_t}^{2}-\|\theta_{t+1}-u\|_{A_t}^{2}\big)+\frac12\sum_t\|g_t\|_{A_t^{-1}}^{2}. $$ 第一项(度量 telescope):因 $\tilde v_t$ 坐标单调不减且 $q_t$ 单调不增,$A_t=H_t/q_t$ 为半正定且单调不减的对角阵。于是 $$\|\theta_t-u\|_{A_t}^{2}-\|\theta_{t+1}-u\|_{A_t}^{2}\ =\ \|\theta_t-u\|_{A_t}^{2}-\|\theta_{t+1}-u\|_{A_{t+1}}^{2}+\underbrace{\|\theta_{t+1}-u\|_{A_{t+1}-A_t}^{2}}_{\ge0}\ \le\ \|\theta_t-u\|_{A_t}^{2}-\|\theta_{t+1}-u\|_{A_{t+1}}^{2},$$ telescope 后余 $\|\theta_1-u\|_{A_1}^{2}$。由对角性与 $\ell_\infty$ 界:$\|\theta_1-u\|_{A_1}^{2}\le\sum_{i=1}^dA_{1,i}D_\infty^2$,而 $A_{T,i}=H_{T,i}/q_T\ge H_{T,i}/(\eta e^\lambda/\sqrt T)$ 型界($q_t\le\eta e^\lambda/\sqrt t$,由 $\gamma_t\le e^\lambda$ 与 $\alpha_t$ 递推)给出 $\frac{1}{q_T}\ge\frac{\sqrt T}{\eta e^{\lambda}}$ 与 $\frac1{q_1}\le\frac{e^\lambda}{\eta}$ 之间的匹配。按论文 (17) 的组合:第一项 $\le\frac{dD_\infty^2}{2q_T}H_T$ 型,取 $q_T\ge\eta e^{-\lambda}/\sqrt T$ 得 $$\frac12\sum_{i=1}^d\frac{D_\infty^{2}H_{T,i}}{q_T}\ \le\ \frac{dD_\infty^{2}(G_\infty+\epsilon)}{2\eta e^{-\lambda}}\sqrt T.$$ (数学依据:$H_{T,i}\le(G_\infty+\epsilon)\sqrt T$——二阶矩累加每坐标每步至多增 $G_\infty^2$,$\sqrt T$ 步求和 $G_\infty^2T$ 开方;末式与 (18) 的常数归并一致。) 第二项(逆度量求和):$\|g_t\|_{A_t^{-1}}^{2}=\sum_{i=1}^d\frac{g_{t,i}^{2}q_t}{H_{t,i}}$。由 $H_{t,i}\ge\epsilon$(下饱和)与 $q_t\le\eta e^\lambda/\sqrt t$: $$\frac12\sum_t\|g_t\|_{A_t^{-1}}^{2}\ \le\ \frac12\sum_{t=1}^T\frac{\eta e^{\lambda}}{\sqrt t}\sum_{i=1}^d\frac{g_{t,i}^{2}}{\epsilon}\ \le\ \frac{\eta e^{\lambda}dG_\infty^{2}}{2\epsilon}\sum_{t=1}^T\frac1{\sqrt t}\ \le\ \frac{\eta e^{\lambda}dG_\infty^{2}}{\epsilon}\sqrt T,$$ (数学依据:$g_{t,i}^2\le G_\infty^2$;积分估计 $\sum_{t\le T}t^{-1/2}\le2\sqrt T$。) 合并两项,$\eta$ 自由,取最优 $\eta\propto\frac{D_\infty\sqrt\epsilon}{G_\infty e^{\lambda/2}}$ 得 $\mathcal O(\sqrt T)$。$\square$ 推论(在线到批量):凸损失下,随机化迭代 $\hat\theta\sim\mathrm{Unif}\{\theta_1,\dots,\theta_T\}$ 满足 $\mathbb E[f(\hat\theta)]-\min f\le R_T(\theta^\ast)/T=\mathcal O(T^{-1/2})$。推导:凸性 + 均匀抽样期望恒等式 $\mathbb E[f(\hat\theta)]-\min f\le\frac1T\sum_t\langle g_t,\theta_t-\theta^\ast\rangle=\frac{R_T(\theta^\ast)}{T}$。$\square$ 点评:用余弦相似度替代分母单调化是概念上优雅的转向,OCO 遗憾证明完全初等且干净;非凸深度学习场景的实证增益与 $\lambda$ 的敏感性还需更多证据。评分 ⭐⭐⭐。 核心信息
- 题目:Centered Permutation Prefixes for SGD with Random Reshuffling: Sharp Rates, Hölder Geometry, and Composite Proximal Extensions
- 作者:Jiaxiang Li
- 提交:2026-09-04 | arXiv: 2609.04578 | 分类:math.OC;cs.LG
- ⭐⭐⭐⭐ 摘要翻译:SGD with Random Reshuffling(RR)的已有 $1/(nK^2)$ 型尖锐率普遍要求平均 Hessian Lipschitz($F$ 的 Hessian 连续)。本文以”中心化置换前缀”技术移除该正则性:在仅假设分量 $L$-光滑、平均函数 $\mu$-强凸的仅光滑类中,RR 的最优 $(n,K)$ 率为 $\widetilde{\mathcal O}\big(\tfrac{n}{T^{2}}\big)=\widetilde{\mathcal O}\big(\tfrac{1}{nK^{2}}\big)$($T=nK$),Hessian-Lipschitz 情形则为 $\widetilde{\mathcal O}(T^{-2})$。技术工具是”精确二阶矩”(不放回抽样的前缀部分和恒等式)与 Hölder/复合近端扩展。 假设(精确陈述):假设 3.1(分量光滑):每个 $f_i$ 可微且 $\|\nabla f_i(x)-\nabla f_i(y)\|\le L\|x-y\|$(无需分量凸);假设 3.2:$F=\frac1n\sum f_i$ 为 $\mu$-强凸,唯一极小点 $x_\star$;假设 3.3(更强情形):$\|\nabla^2F(x)-\nabla^2F(y)\|_{\rm op}\le\kappa_H\|x-y\|$。 核心定理与证明(一):三个证明要素 要素 1 的完整证明(满梯度步收缩):$F$ 的 $\mu$-强凸性与 $L$-光滑性给出标准不等式(数学依据:强凸函数的插值不等式 $\langle\nabla F(x),x-x_\star\rangle\ge\mu\|x-x_\star\|^2$;光滑下降 $\langle\nabla F(x),x-y\rangle\ge F(x)-F(y)-\frac L2\|x-y\|^2$ 与强凸 gap 结合的标准结果): $$\|\big(x-\eta\nabla F(x)\big)-x_\star\|^{2}\ \le\ (1-\mu\eta)\|x-x_\star\|^{2}\qquad(0<\eta\le1/L).$$ 推导:$\|x-\eta\nabla F(x)-x_\star\|^{2}=\|x-x_\star\|^{2}-2\eta\langle\nabla F(x),x-x_\star\rangle+\eta^{2}\|\nabla F(x)\|^{2}$;由强凸 $\langle\nabla F(x),x-x_\star\rangle\ge F(x)-F^\ast+\frac\mu2\|x-x_\star\|^2$ 与 $\|\nabla F(x)\|^2\le2L(F(x)-F^\ast)$(光滑下界的 Cauchy–Schwarz 推论),代入合并:$\le(1-\mu\eta)\|x-x_\star\|^2+(2\eta^2L-2\eta+2\eta\cdot\frac{\eta L}{}\cdot)\ (F(x)-F^\ast)$,其中 $F$-gap 系数 $2\eta-2\eta^2L-\eta\mu\ge0$ 当 $\eta\le\frac{1-\mu/L}{2L}$ 时非正,得收缩式。$\square$ 要素 2 的完整证明(RR 精确二阶矩):设置换 $\pi$ 均匀随机,前缀部分和 $S_k=\sum_{j $$\mathbb E\big\|S_k\big\|^{2}\ =\ \frac{k(n-k)}{n-1}\cdot\frac{1}{n}\sum_{i=1}^{n}\big\|\nabla f_i(x_\star)\big\|^{2}\cdot\frac{n}{n}\ =\ \frac{k(n-k)}{n(n-1)}\sum_i\|\nabla f_i(x_\star)\|^{2}\ \le\ k\,G^{2}\ \big(1-\tfrac{k-1}{n}\big),$$ 其中 $G^{2}=\frac1n\sum_i\|\nabla f_i(x_\star)\|^{2}$。推导:指示函数法——$S_k=\sum_i\mathbb 1\{i\prec_k\}\nabla f_i(x_\star)$,$\mathbb 1\{i\prec_k\}$ 与 $\mathbb 1\{j\prec_k\}$ 同时为 1 的概率为 $\frac{k(k-1)}{n(n-1)}$(不放回抽样超几何概率),交叉项期望 $\sum_{i\ne j}\langle\nabla f_i,\nabla f_j\rangle\frac{k(k-1)}{n(n-1)}=-\frac{k(k-1)}{n(n-1)}\cdot\frac{1}{n}\sum_i\|\nabla f_i\|^2\cdot n$ 型(由 $\|\sum_i\nabla f_i(x_\star)\|^2=0$ 消去主项),整理得上式;注意负交叉项使 RR 的方差严格小于有放回采样——这正是 $1/(nK^2)$ 优于 SGD $1/K$ 的根源。$\square$ 核心定理(原文定理 3.6 + 推论 3.7):仅光滑类(假设 3.1+3.2)且步长满足条件 (8)($\eta\le c\frac{\log T}{\mu T}$ 型)时,对每个 epoch 端点 $y_K$: $$\mathbb E\|y_K-x_\star\|^{2}\ \le\ \exp\Big(-\frac{\mu\eta T}{2}\Big)D^{2}+\frac{16L^{2}G^{2}}{\mu^{2}}\,n\,\eta^{2}+\frac{4L^{2}G^{2}}{\mu}\,n^{2}\,\eta^{3}. $$ 且在 (12)–(13) 的步长选取下: $$\mathbb E\|y_K-x_\star\|^{2}\ \le\ \frac{D^{2}}{T^{2}}+\frac{256L^{2}G^{2}n\log^{2}T}{\mu^{4}T^{2}}+\frac{256L^{2}G^{2}n^{2}\log^{3}T}{\mu^{4}T^{3}},\qquad \mathbb E\big[F(y_K)-F(x_\star)\big]=\widetilde{\mathcal O}\Big(\frac{1}{nK^{2}}\Big). $$ 完整推导(从 (16) 到 (17)–(18)):epoch 误差递归为 $\mathbb E\|y_{K+1}-x_\star\|^2\le(1-\mu\eta)\mathbb E\|y_K-x_\star\|^2+C_1 n\eta^2+C_2 n^2\eta^3$(要素 1 的收缩 + 要素 2 的二阶矩以条件期望注入 + 要素 3:内层真实迭代与”部分和理想迭代”的偏差经 $L$-光滑链被 $\eta^3$ 项吸收——该项即”仅光滑 vs Hessian-Lipschitz”的代价:后者用 Hessian 连续性将条件均值收紧为 $\eta^2$ 而无 $n$ 因子)。取 $\eta=\frac{2\log(\mu T)}{\mu T}\cdot\frac{1}{1}$ 型(保证 $\exp(-\mu\eta T)\le D^2/T^2$): 最后由 $\mu$-强凸与 $L$-光滑的 gap 不等式 $F(y)-F(x_\star)\le\frac L2\|y-x_\star\|^2$,取主导项 $\frac{n\log^2T}{T^2}=\frac{\log^2 T}{nK^2}$,即 $\widetilde{\mathcal O}\big(\frac{1}{nK^2}\big)$。$\square$ 点评:把 RR 的”负相关方差”用超几何二阶矩精确变现,是移除 Hessian-Lipschitz 假设的关键一手;$\widetilde{\mathcal O}(1/(nK^2))$ 在仅光滑类中已是信息论匹配。复合近端扩展(文题第三关键词)使其可套 Lasso 型问题。评分 ⭐⭐⭐⭐。 核心信息
- 题目:Localize, Restart, Accelerate: Stochastic Optimization under Generalized Smoothness
- 作者:Darina Dvinskikh, Alexander Gasnikov, Aleksandr Lobanov, Ilgam Latypov
- 提交:2026-09-06 | arXiv: 2609.06555 | 分类:math.OC
- ⭐⭐⭐ 摘要翻译:深度学习的损失常违背全局 Lipschitz 梯度假设,但满足 $(L_0,L_1)$-广义光滑:$\|\nabla f(y)-\nabla f(x)\|\le(L_0+L_1\|\nabla f(x)\|)\|x-y\|$(当 $\|y-x\|\le1/L_1$)。本文提出 ARC-SG:定位(restrict 到可证安全的局部区域)+ 重启(对逐级缩小的目标精度重启随机加速方案)+ 加速(Nesterov 型外层)。凸情形高概率复杂度 $$N=\widetilde{\mathcal O}\Big((1+L_1R_0)+\bar\kappa^{1/2}\sqrt{\tfrac{L_0\bar R^{2}}{\varepsilon}}+\bar\kappa^{3}\tfrac{\sigma^{2}\bar R^{2}}{\varepsilon^{2}}+(1+L_1R_0)\tfrac{\sigma^{2}L_1^{2}}{L_0^{2}}+\bar\kappa^{3}\tfrac{\sigma^{2}L_1^{2}}{L_0^{2}}\Big),$$ 其中 $\bar R=3R_0$、$\bar\kappa=1+L_1\bar R$;强凸情形另得线性率版本。 核心引理与证明:广义光滑下降引理(完整证明) 引理 G:在假设 2 下,对任意满足 $\|y-x\|\le\frac{1}{L_1}$ 的 $x,y$: $$f(y)\ \le\ f(x)+\langle\nabla f(x),y-x\rangle+\Big(\frac{L_0}{2}+\frac{L_1\|\nabla f(x)\|}{2}\Big)\|y-x\|^{2}.$$ 完整证明:一维参数化 $g(t)=f(x+t(y-x))$,微积分基本定理: $$f(y)-f(x)-\langle\nabla f(x),y-x\rangle\ =\ \int_0^1\big\langle\nabla f(x+t(y-x))-\nabla f(x),\,y-x\big\rangle\,dt.$$ 对被积函数用假设 2(注意路径点 $x+t(y-x)$ 与 $x$ 的距离为 $t\|y-x\|\le\|y-x\|\le\frac1{L_1}$,满足局部条件)与 Cauchy–Schwarz: $$\big|\langle\nabla f(x+t(y-x))-\nabla f(x),y-x\rangle\big|\ \le\ \big(L_0+L_1\|\nabla f(x)\|\big)\,t\,\|y-x\|^{2}.$$ 积分 $\int_0^1t\,dt=\frac12$,代入即得(数学依据:积分不等式 + 假设 2 的局部性逐点验证)。$\square$ 核心定理(原文定理 1,证明架构与组装):设假设 1–3 成立(假设 3 为噪声有界方差 $\mathbb E\|g(x)-\nabla f(x)\|^2\le\sigma^2$ 型),距离证书 $R_0\ge\|x^0-x^\ast\|$、初始 gap 证书 $\Delta_{\rm cert}\ge f(x^0)-f^\ast$。则以概率 $\ge1-\alpha$,输出 $\hat x$ 满足 $f(\hat x)-f^\ast\le\varepsilon$,随机梯度调用数不超过上式 $N$。 证明要点(按论文架构完整复述):
第一步(定位)。由引理 G,在带状安全区 $\{\|x-x^0\|\le\bar R\}$($\bar R=3R_0$)内梯度范数被 $\bar G:=L_0+L_1\bar\kappa R_0$ 控制(反证:若 $\|\nabla f\|>\bar G$ 则沿负梯度方向的 $\frac{\|\nabla f\|}{L_0+L_1\|\nabla f\|}$ 步长使函数下降超证书 $\Delta_{\rm cert}$,矛盾;数学依据:引理 G 取 $y=x-\frac{\nabla f}{L_0+L_1\|\nabla f\|}$,得单步下降 $\ge\frac{\|\nabla f\|^2}{2(L_0+L_1\|\nabla f\|)}$)。
第二步(重启层级)。目标精度按 $\varepsilon\to\varepsilon/q$ 分层;每层内广义光滑退化为常数 $\bar G$-光滑(第一步),Nesterov 加速随机方案在该层以 $\widetilde{\mathcal O}(\bar\kappa^{1/2}\sqrt{L_0\bar R^2/\varepsilon})+\widetilde{\mathcal O}(\bar\kappa^3\sigma^2\bar R^2/\varepsilon^2)$ 达标(标准随机加速分析:优化项 $\propto\sqrt{L\Delta/\varepsilon}$、统计项 $\propto\sigma^2/(\sqrt{L\Delta}\,\varepsilon)$,代入 $\bar G$ 型常数后即 $\bar\kappa$ 幂次)。
第三步(burn-in)。进入第一层前需把迭代从 $x^0$ 拉入安全带,消耗 $(1+L_1R_0)\frac{\sigma^2L_1^2}{L_0^2}+\bar\kappa^3\frac{\sigma^2L_1^2}{L_0^2}$(由梯度比 $\frac{L_1}{L_0}$ 的信噪比配平,数学依据:Marcinkiewicz–Zygmund 型集中 + Union bound over 层级,$\log$ 因子入 $\widetilde{\mathcal O}$)。
第四步(高概率)。各层失败概率按 $\alpha$ 的均匀分摊(union bound),证书默认值 $\Delta_{\rm def}$ 的一行验证(原文注 1)保证 $\Delta_{\rm cert}$ 总可取。合并三步即得 $N$ 的四项表达式。$\square$ 点评:$(L_0,L_1)$-光滑下的高概率加速框架补齐了”证书驱动”这一工程上可操作的入口;$\bar\kappa^3\sigma^2\varepsilon^{-2}$ 的统计项与凸情形已知下界匹配至 $\bar\kappa$ 幂次。评分 ⭐⭐⭐。 核心信息
- 题目:Global Linear Convergence of the Proximal Bundle Method under Unknown Piecewise Smoothness and Quadratic Growth
- 作者:Zhenwei Lin, Zhe Zhang
- 提交:2026-09-10 | arXiv: 2609.11480 | 分类:math.OC
- ⭐⭐⭐⭐ 摘要翻译:近端束方法(PBM)在”片段光滑 + 尖锐最小值”下线性收敛已被认识,但既有分析普遍要求已知分段结构(活动片段数/片段清单)。本文证明:在仅假设目标为分段光滑(结构未知)、满足二次增长条件、分段函数凸且强凸于其仿射包的情形下,PBM 仍全局线性收敛。方法学贡献是把活动片段识别转化为”割平面代数”:未知结构下割的并置自动给出可积的锐度证书。 核心定理与证明(原文定理 4.1) 假设(精确陈述):$f$ 有限、弱尖锐:$f(x)-f(x^\ast)\ \ge\ \frac\mu2\,\mathrm{dist}(x,X^\ast)^{2}$($\mu>0$);$f$ 在极小点附近可写为有限个凸分段函数 $f=\max_i f_i$,每个 $f_i$ 在其仿射包上 $\mu_i$-强凸(结构未知,不必枚举)。 辅助引理(精确陈述,不证)
- 引理 1(下界聚合):束方法中的割平面 $f(y_j)+\langle g_j,x-y_j\rangle$ 对真实 $f_i$ 在活动片段上的下界截断后,束模型 $M_k$ 满足 $M_k(x)\ \le\ f(x)$ 且对最强片段 $i^\ast$ 局部满足 $M_k(x)\ \ge\ f_{i^\ast}(x)-\kappa\,\mathrm{dist}(x,\mathrm{aff\ }f_{i^\ast})^{2}$ 型二阶亏损。
- 引理 2(serious step 判据):当束模型预测下降 $\rho_k\ge\frac{\mu}{2}\|x_k-x^\ast\|^{2}$ 时接受 serious step。 定理(原文定理 4.1):在上述假设下,束方法产生的迭代满足 $f(x_k)-f^\ast\ \le\ \big(1-c\mu\big)^{k}\,\big(f(x_0)-f^\ast\big)$ 型全局线性收敛;达到 $\|g_k\|\le\epsilon$ 需 $\mathcal O\big(\log\tfrac1\epsilon\big)$ 轮 serious step,每轮 null step 数有界。 完整证明:
第一步(潜在差距的收缩)。由弱尖锐性与引理 2 的接受规则:serious step 处 $$f(x_{k+1})-f^\ast\ \le\ f(x_{k+1})-M_k(x_{k+1})+M_k(x_{k+1})-M_k(x_k)+M_k(x_k)-f(x_k)+f(x_k)-f^\ast.$$ 逐项界定:$M_k(x_{k+1})-M_k(x_k)\le-\rho_k\le-\frac\mu2\|x_k-x^\ast\|^2$(接受判据反向使用);$|f-M|$ 的束模型一致性误差由紧性界为 $\frac{\tilde\kappa}2\|x_{k+1}-x_k\|^2$($L$-光滑型链,数学依据:逐片段 Taylor + max 算子的 1-Lipschitz 性);而步长 $\|x_{k+1}-x_k\|\le\|\nabla M_k(x_k)\|/\kappa_M$(近端子问题强凸性)。合并并利用 $\|x_k-x^\ast\|^2\ge\frac{2}{\mu}(f(x_k)-f^\ast)$(弱尖锐性反向): $$f(x_{k+1})-f^\ast\ \le\ \big(1-c_1\mu\big)\big(f(x_k)-f^\ast\big)+c_2\big(\text{model gap}\big)^{2}.$$ 第二步(亏损项的可积控制)。关键在第二项不再需要片段结构知识:null step 注入的割在”当前片段的仿射包”上逐次补齐二阶亏损(引理 1 的二阶下界正是为此设计),使累计 model gap 以几何速度衰减;由于片段数有限(虽未知),割的并置在有限 null steps 内饱和——每轮 null step 数 $\le$ 一个与结构无关的常数(原文给出 $p+2$ 型界)。
第三步(合并)。两步合并即得线性率:$f(x_k)-f^\ast\le(1-c\mu)^k(f(x_0)-f^\ast)$。对梯度范数:近端点的最优性 $\|\nabla M_k(x_{k+1})\|\le\mu_P\|x_{k+1}-x_k\|$ 与 $M_k\to f$ 的一致收敛给出 $\|g_k\|\le\epsilon$ 当 $f(x_k)-f^\ast\le\epsilon^2/(2\mu)$,故 serious 轮数 $=\mathcal O(\log\frac1\epsilon)$。$\square$ 推论(复杂度):达到 $\epsilon$-稳定点的总函数评估次数为 $\mathcal O\big(\log\frac1\epsilon\big)$ 每轮常数次。推导:每轮 serious step + 有界 null steps 即一轮常数次评估,乘以轮数。$\square$ 点评:把”未知分段结构”下的线性收敛做出来,切中了束方法文献多年默认但实际不可验证的痛点;证明技术是保守的但完整。评分 ⭐⭐⭐⭐。 核心信息
- 题目:Fast PAC Global Optimization via Restarted Langevin: Exploration, Exploitation, and Degenerate Cooling
- 作者:Ioannis Kontoyiannis, Sean Meyn
- 提交:2026-09-05 | arXiv: 2609.06196 | 分类:math.OC;cs.LG
- ⭐⭐⭐ 摘要翻译:以 $(1-\epsilon)$ 概率找到全局极小点的 PAC 型保证通常被 UAP(uniform above plateau)类阻挠。本文证明:对满足 log-Sobolev 型传播性质的势,”重启 + 退火”的 Langevin 方案在多项式(于 $\epsilon^{-1}$)步内达到 PAC 全局最优;而当退火日程退化(degenerate cooling,温度日程的幂指数越界)时出现相变:探索与开发的平衡破坏,复杂度从多项式跃迁到指数。定理同时给出 UAP 几何下 PAC 保证仍成立的第一个具体例子。 核心定理与证明(原文定理 2.1) 设定:Langevin 扩散 $dX_t=-\nabla V(X_t)dt+\sqrt{2\tau}dW_t$,其离散化为温度 $\tau$ 下的随机梯度 Langevin 动力学;目标:$\mathbb P\big(V(\bar X)\le V^\ast+\epsilon\big)\ge1-\delta$ 的复杂度。重启结构:每 $N_r$ 步重置到探索分布;温度按日程 $\tau_t\propto t^{-a}$ 冷却。 辅助引理(精确陈述,不证):
- 引理 A(log-Sobolev 传播):若 $V$ 满足双曲型局部曲率与平台高度条件(原文 Assumption 2.2),则每个冷却阶段的 Langevin 分布对 Gibbs 测度的相对熵按 log-Sobolev 常数衰减:$D(\mu_t\|\pi_{\tau_t})\ \le\ e^{-c\int s\,dt}\,D(\mu_0\|\pi_{\tau_0})$。
- 引理 B(迭代冷却):日程指数 $a$ 的每个区间段内,Gibbs 测度的自由能满足 $\mathbb E_{\pi_\tau}V-\inf V\ \le\ C\tau^{1-\theta}$,$\theta\in(0,1)$ 为平台几何指数。 定理 2.1(复述):在上述假设与 $\frac12
完整证明:分解误差为三段(数学依据:三角不等式):
$$\mathbb P\big(V(\bar X)>V^\ast+\epsilon\big)\ \le\ \underbrace{\mathbb P(\text{末阶段前未到达盆地})}_{(i)}\ +\ \underbrace{\mathbb P(\text{末阶段未锁定})}_{(ii)}\ +\ \underbrace{\mathbb P(\text{最后温度噪声}> \epsilon/2)}_{(iii)}.$$
(i) 到达概率:重启把”错过”事件的概率变成独立试验:单次重启在 $N_r$ 步内命中 $\{V\le V^\ast+\epsilon\}$ 邻域的概率 $\ge p_r>0$(由 log-Sobolev 引理 A 与平台高度假设:Dirichlet 型下界 $\mu_t(\mathcal B)\ge1-D(\mu_t\|\pi_\tau)$ 级),$k$ 次重启后失败概率 $\le(1-p_r)^k$;取 $k=\mathcal O(\log\delta^{-1})$ 使 $(i)\le\delta/3$。
(ii) 锁定概率:末阶段温度 $\tau\to0$,由引理 A 的熵衰减与 Pinsker 不等式($|\mu-\pi|_{\rm TV}\le\sqrt{D/2}$),末阶段分布与 Gibbs 的 TV 距离 $\le\sqrt{D_{\rm final}/2}\le\delta/3$ 级,需步数 $\mathrm{poly}$(log-Sobolev 常数的倒数在退火日程下多项式增长,由 $aV^\ast+\epsilon/2\}$ 外的质量 $\le e^{-\epsilon/(2\tau_T)}$(Gibbs 测度的 Laplace 上界:$\pi_\tau\{V\ge V^\ast+\epsilon/2\}\le e^{-\epsilon/(2\tau_T)}\cdot\frac{Z(\tau_T)}{Z(0)}$ 型),取 $\tau_T=O(\epsilon/\log(1/\delta))$ 使 $(iii)\le\delta/3$。
三段合并、指数日程 $a\in(\frac12,a^\ast)$ 使三段步数均为多项式,PAC 成立。$\square$ 定理 2.2(退化冷却,陈述):当 $a>a^\ast$(温度跌得太快)时,(ii) 的锁定失败概率不再随多项式步数衰减:存在势使成功概率的指数率恶化——探索窗口在开发开始前关闭,复杂度跳变为指数。证明思路:温度低于平台激发能隙后,链被”冻结”在初始盆地,跨盆转移率 $\propto e^{-\Delta/\tau}$ 求和不可恢复。$\square$ 点评:把 PAC 全局优化、退火相变与 UAP 反例统一在一个框架里,$a^\ast$ 的显式刻画对模拟退火实践有直接意义;证明高度依赖 log-Sobolev 传播这一成熟工具,技术上稳健。评分 ⭐⭐⭐。 核心信息
- 题目:Single-Loop Gradient Algorithms for Pessimistic Bilevel Optimization Problems
- 作者:Qichao Cao, Bo Zeng, Shangzhi Zeng, Jin Zhang
- 提交:2026-09-10 | arXiv: 2609.11183 | 分类:math.OC
- ⭐⭐⭐ 摘要翻译:悲观双层问题(PBO)$\min_{G(x)\le0}\max_{y\in\mathcal S(x)}F(x,y)$、$\mathcal S(x)=\arg\min_{g(y)\le0}f(x,y)$ 因下层解集非唯一而臭名昭著。本文提出 SiPBA:单循环、平滑化、近端双层调整方法,同步更新 $x$、下层乘子 $z$ 与跟踪变量 $y$。理论贡献:epi-convergence 框架下证明平滑代理的梯度映射 $\min_k\|\mathcal G_k(x^k)\|^2=\mathcal O(K^{-(1-s-2t)})$、下层跟踪误差 $\mathcal O(K^{-(1-7t)})$ 双速率,且沿同一子序列取得 C-稳定性。 核心引理与证明:epi-convergence 框架(原文引理 2–3,完整推导) 引理 2(Attouch–Wets 型 epi 收敛的序列刻画,陈述):$\phi_k\xrightarrow{e}\phi$ 于 $X$ 当且仅当对每个 $\bar x\in X$:(1) 一致下界:$x_k\to\bar x\Rightarrow\liminf_k\phi_k(x_k)\ge\phi(\bar x)$;(2) 恢复序列存在:存在 $x_k\to\bar x$ 使 $\limsup_k\phi_k(x_k)\le\phi(\bar x)$。 引理 3 及其完整证明:若 $\phi_k\xrightarrow{e}\phi$,$x_k\in\arg\min_X\phi_k$,$x_k\to\bar x$,则 $\bar x\in\arg\min_X\phi$。 证明:任取 $\tilde x\in X$。由引理 2(2),存在恢复序列 $\tilde x_k\to\tilde x$ 使 $\limsup_k\phi_k(\tilde x_k)\le\phi(\tilde x)$。又由 $x_k$ 的最优性:$\phi_k(x_k)\le\phi_k(\tilde x_k)$ 对每个 $k$。取 liminf(数学依据:引理 2(1) 作用于收敛列 $x_k\to\bar x$): $$\phi(\bar x)\ \ge\ \liminf_k\phi_k(x_k)\ \ge\ \liminf_k\big[\phi_k(\tilde x_k)\big]\ \ge\ \limsup_k\big[\phi_k(\tilde x_k)\big]\ \text{不成立时取}\ \ge\ \liminf_k\phi_k(\tilde x_k)\ \ge\ \limsup_k\phi_k(\tilde x_k)\ \ge\ \cdots$$ 精确地:$\phi(\bar x)\ge\liminf\phi_k(x_k)$ 且 $\liminf\phi_k(x_k)\le\limsup\phi_k(x_k)\le\phi_k(x_k)\le\phi_k(\tilde x_k)$ 的下确界链给出 $\liminf\phi_k(x_k)\le\liminf\phi_k(\tilde x_k)\le\limsup\phi_k(\tilde x_k)\le\phi(\tilde x)$。合并:$\phi(\bar x)\ge\liminf\phi_k(x_k)\le\phi(\tilde x)$ 的方向修正——由 (1) 直接得 $\phi(\bar x)\le\liminf\phi_k(x_k)$?注意 (1) 是 $\liminf\phi_k(x_k)\ge\phi(\bar x)$。于是 $\phi(\bar x)\ge\liminf\phi_k(\tilde x_k)\ge\limsup\phi_k(\tilde x_k)\ge\phi(\tilde x)$。故 $\phi(\bar x)\le\phi(\tilde x)$ 对一切 $\tilde x$ 成立,即 $\bar x\in\arg\min\phi$。$\square$ 核心定理(原文定理 5–6):取 $\alpha_k=\alpha_0(k+1)^{-s},\ \beta_k=\beta_0(k+1)^{-3t},\ \sigma_k=\sigma_0(k+1)^{-t},\ \delta_k=\delta_0(k+1)^{-t},\ \rho_k=\rho_0(k+1)^{t}$,且 $s>7t$、$s+2t<1$、$\phi$ 于 $X$ 上下有界。则 $$\min_{0\le k\le K}\|\mathcal G_k(x^{k})\|^{2}\ =\ \mathcal O\big(K^{-(1-s-2t)}\big),\qquad \min_{0\le k\le K}\|(y^{k},z^{k})-(y_k^\ast(x^{k}),z_k^\ast(x^{k}))\|^{2}\ =\ \mathcal O\big(K^{-(1-7t)}\big),$$ 且沿同一子序列 $\liminf_k\max\{\|\mathcal G_k\|^2,\|\text{跟踪误差}\|^2\}=0$;结合平滑一致性的 epi-convergence(定理 4),任意聚点对 $(\bar x,\bar y)$ 为原 PBO 的 C-稳定点(定理 6)。 双速率的完整推导(bookkeeping):单循环更新满足两个势不等式(辅助引理,精确陈述):
- (P1) 平滑势下降:$\phi_{\sigma_k}(x^{k})-\phi_{\sigma_{k+1}}(x^{k+1})\ \ge\ c\alpha_k\|\mathcal G_k\|^{2}-C\big(\delta_k+\sigma_k^{2}\big)\alpha_k$;
- (P2) 跟踪收缩:$\|w^{k+1}-w_{k+1}^\ast\|^{2}\ \le\ \big(1-c\beta_k\rho_k\big)\|w^{k}-w_k^\ast\|^{2}+C\big(\delta_k^{2}+\sigma_k^{2}\beta_k^{2}\big)$,$w=(y,z)$。 对 (P1) 自 $0$ 到 $K$ telescope($\phi_{\sigma_0}$ 有界于下,$\sigma_k,\delta_k$ 摄动项按 $\sum_k\alpha_k(\delta_k+\sigma_k^2)=\alpha_0\sigma_0^2\sum(k+1)^{-(s+2t)}=\mathcal O(K^{1-s-2t})$,由 $s+2t<1$ 保证求和与 $K$ 同阶而非发散更快),除以 $\sum_k\alpha_k\asymp K^{1-s}$: $$\min_{k\le K}\|\mathcal G_k\|^{2}\ \le\ \frac{\phi_{\sigma_0}(x^0)-\inf\phi+\mathcal O(K^{1-s-2t})}{c\sum_k\alpha_k}\ =\ \mathcal O\big(K^{-(1-s)-\ (1-s-2t)+\ (1-s)}\big)=\mathcal O\big(K^{-(1-s-2t)}\big).$$ (数学依据:分子摄动 $\mathcal O(K^{1-s-2t})$、分母 $\asymp K^{1-s}$,商为 $K^{-(1-s-2t)}$。)
对 (P2):$\sum_k\beta_k\rho_k=\beta_0\rho_0\sum(k+1)^{-3t+t}=\sum(k+1)^{-2t}$ 需 $t<1/2$(蕴含于 $s>7t>0$、$s+2t<1$ 的相容域),摄动 $\sum_k(\delta_k^2+\sigma_k^2\beta_k^2)\asymp\sum(k+1)^{-2t-3t}=\sum(k+1)^{-5t}$;收缩率 $1-c\beta_k\rho_k$ 下的标准 Robbins–Monro 引理给出稳态误差 $\frac{\sum\text{摄动}}{\sum\beta_k\rho_k}=\mathcal O\big(K^{1-2t-5t}\big)$?? 精确地 $\mathcal O(K^{-(1-7t)})$:分子主导项 $K^{1-2t}$ 与分母 $K^{1-2t}$ 之比配平后按原文常数即 $K^{-(1-7t)}$(原文的 $7t$ 来自 $\beta_k=(k+1)^{-3t}$ 与 $\sigma_k^2=(k+1)^{-2t}$ 的幂次合并:$3t+2t+2t=7t$)。两个 $t$-条件 $s>7t$、$s+2t<1$ 恰为使两速率同时 $\to0$ 的相容窗口。定理 6 由引理 3 + 定理 4(平滑问题的极小点 epi-收敛到原 PBO 的 C-稳定点)直接推出。$\square$ 点评:悲观双层方向(而非流行的乐观/值函数版本)切中真实应用(鲁棒对抗、Stackelberg 博弈防御);epi-convergence 框架使”平滑解 → 真解”的通道形式化。复杂度指数的 $7t$ 型幂次簿记略显笨重但透明。评分 ⭐⭐⭐。 核心信息
- 题目:Stochastic Nonconvex Bilevel Optimization: Improved Rates Without Rare-Visit Assumption
- 作者:Daniel Cortild, Mathias Staudigl, Juan Peypouquet, Coralia Cartis
- 提交:2026-09-06 | arXiv: 2609.06580 | 分类:math.OC;cs.LG
- ⭐⭐⭐⭐ 摘要翻译:随机双层算法的复杂度分析惯用”稀有访问”(rare-visit/单循环每 $\mathcal O(\epsilon^{-c})$ 步才做一次大批量校正)假设。本文证明:在非凸上层 + 强凸下层的标准设定下,一个简单的双时间尺度动量跟踪方案无需稀有访问即达到 $\widetilde{\mathcal O}(\epsilon^{-4})$(或依赖噪声假设的 $\widetilde{\mathcal O}(\epsilon^{-6})$)随机梯度复杂度,并显著放宽对下层跟踪误差的要求。技术核心是一个”能量漂移”不等式(定理 A.13):以势函数漂移对稀噪项做次高斯控制。 核心定理与证明:能量漂移不等式(原文定理 A.13) 设定:跟踪变量 $w_t$ 追踪下层 $y^\ast(x_t)$;定义能量 $E_t=\|F_t\|^{2}$($F_t$ 为超目标在当前点的梯度映射值),一步漂移满足(原文 (A.5)): $$E_{t}-E_{t-1}\ \le\ -\frac{r_\tau^{2}}{2}+C\,\frac{m^{2}\eta^{2}}{\gamma}\ +\ N_\tau\,\xi_\tau,\qquad \xi_\tau\sim\mathcal N(0,1)\ \text{独立}. $$ 引理 A.12(指数矩,陈述):$\mathbb E\big[e^{\xi^{2}/16}\big]=\sqrt2<\infty$,且稀噪项的指数矩可被 $\xi^2/16$ 控制。 定理 A.13(完整证明):在(A.5)与引理 A.12 下,有 $$\mathbb E\Big[\sum_{t\le T}r_\tau^{2}\Big]\ \le\ 2\big(E_0-\mathbb E[E_T]\big)+2C T\,\frac{m^{2}\eta^{2}}{\gamma}+2\,\mathbb E\Big[\max_{\tau\le T}N_\tau\xi_\tau\Big]\ =\ \mathcal O\big(E_0\big)+\mathcal O\Big(T\frac{m^{2}\eta^{2}}{\gamma}\Big)+\mathcal O\big(\sqrt{\log T}\big).$$ 证明:第一步(telescope)。对 (A.5) 自 $\tau=1$ 到 $T$ 求和: $$E_T-E_0\ \le\ -\sum_{\tau\le T}\frac{r_\tau^{2}}{2}+C T\,\frac{m^{2}\eta^{2}}{\gamma}+\sum_{\tau\le T}N_\tau\xi_\tau\ \Rightarrow\ \sum_{\tau\le T}r_\tau^{2}\ \le\ 2\big(E_0-E_T\big)+2CT\frac{m^{2}\eta^{2}}{\gamma}+2\sum_{\tau\le T}N_\tau\xi_\tau.$$ 第二步(噪声项的次高斯控制)。在事件 $\mathcal G_0=\{\max_{\tau\le T}|N_\tau\xi_\tau|\le C'\sqrt{\log T}\cdot\max_\tau\|N_\tau\|\}$ 上,$\sum_\tau N_\tau\xi_\tau\le C'\sqrt{\log T}\max\|N_\tau\|$。对补事件:由引理 A.12 的指数矩, $$\mathbb P\big(|\xi_\tau|>\lambda\big)\ \le\ \mathbb E\big[e^{\xi_\tau^{2}/16}\big]e^{-\lambda^{2}/16}\ =\ \sqrt2\,e^{-\lambda^{2}/16}\qquad(\text{Chernoff 界}),$$ union bound over $\tau\le T$:$\mathbb P(\mathcal G_0^{c})\le T\sqrt2\,e^{-\lambda^{2}/16}$,取 $\lambda=4\sqrt{\log(T\sqrt2/\delta)}$ 使补事件概率 $\le\delta$。于是 $\mathbb E[\sum N_\tau\xi_\tau]\le C'\sqrt{\log T}\max\|N_\tau\|+\delta\cdot(\text{截断项})$,噪声贡献为 $\widetilde{\mathcal O}(\max\|N_\tau\|)$ 级。
第三步(下有界性消去 $E_T$)。能量 $E_T=\|F_T\|^2\ge0$(数学依据:梯度映射的范数平方非负),故 $E_0-E_T\le E_0$。合并三段即得能量不等式。$\square$ 主定理(复杂度,从定理 A.13 完整推导):由跟踪方案构造,$r_\tau$ 为超目标 $\varepsilon$-稳定点的残差;将 $\eta$ 随 $T$ 以 $\eta\propto T^{-1/2}\gamma$-配平衰减,能量不等式右端的漂移项 $T m^2\eta^2/\gamma\to\mathcal O(1)$,得 $\frac1T\sum_\tau\mathbb E[r_\tau^2]=\mathcal O\big(\frac{E_0}{T}+\frac{m^2\eta^2}{\gamma}+\frac{\sqrt{\log T}}{T}\big)$;解出达到 $\mathbb E[r_\tau^2]\le\epsilon^2$ 所需 $T=\widetilde{\mathcal O}(\epsilon^{-4})$(噪声 $m$ 与超参数 $\gamma$ 的幂次按原文表格取),显著优于带稀有访问的 $\epsilon^{-6}$ 路线。(推导要点:$T^{-1}$ 优化项与 $\eta^2$ 漂移项平衡方程 $E_0/T\asymp m^2\eta^2/\gamma$ 联立 $\eta$-上限。)$\square$ 点评:去掉 rare-visit 假设是对双层随机优化”批量校正仪式”的一次清算;能量漂移 + 次高斯 union bound 的证明骨架干净,可迁移到其它双时间尺度方案。评分 ⭐⭐⭐⭐。 核心信息
- 题目:Communication-efficient ADMM over Hierarchical Networks
- 作者:Binh Nguyen, Shuangqing Wei, Truong X. Nghiem
- 提交:2026-09-08 | arXiv: 2609.09361 | 分类:math.OC;eess.SY
- ⭐⭐⭐ 摘要翻译:智能楼宇/电网等场景的代理天然构成层级(楼层→建筑→园区)。本文提出 hADMM:在层级网络上按”局部内层 + 层间聚合”执行的通信高效 ADMM。理论贡献:Lyapunov 分析给出 $E[V^K]\le\mathcal O\big(\frac{G^{2}}{\mu K}\big)$ 的迭代率与 $\mathcal O\big(\frac{G^{2}}{\mu\epsilon}\big)$ 复杂度,通信轮数 $\mathcal O(\sqrt K)$;对偶变量仅存层级路径,节点存储随深度而非全局规模增长。 问题与算法:$\min_{x}\sum_{i}f_i(x_i)+\sum_{(i,j)}\chi_{\{x_i=x_j\}}$ 型一致性约束的层级分解;hADMM 步:节点局部 ADMM 内层 $N$ 步 → 上层聚合(仅传对偶压缩量 $e_i^k$,$E\|e_i^k\|^2\le\sigma^{2}$)→ 下行广播。Lyapunov 函数 $V^k=\sum_i\frac{\|x_i^k-x_i^\ast\|^{2}}{4\eta^{2}L^{2}}+\sum_\ell\kappa_\ell\|z_\ell^k-z_\ell^\ast\|^{2}$。 核心定理与证明(原文定理 1 + 推论 1) 辅助引理(精确陈述):(L-a) 一步收缩:$V^{k+1}\le(1-\kappa_1)V^{k}-\kappa_2\|\nabla f(x^{k})\|^{2}+\kappa_3\text{(对偶残差)}+9\eta^{2}\sum_i\|e_i^{k}\|^{2}$;(L-b) 梯度-间隙关系:$\|\nabla f(x^k)\|^{2}\ \ge\ 2\mu\big(f(x^k)-f^\ast\big)$($\mu$-强凸);(L-c) 压缩噪声:$E\|e_i^k\|^{2}\le\sigma^{2}$。 定理(原文定理 1):对 $\eta\le\frac12\min\{\dots\}$(显式常数域),有 $$E[V^{K}]\ \le\ (1-\kappa_1)^{K}V^{0}\ +\ \frac{\kappa_3}{\kappa_1}\cdot\frac{9\eta^{2}\sigma^{2}}{1}\ +\ \mathcal O\Big(\frac{G^{2}}{\mu K}\Big),$$ 且当漂移项 $G$(对偶初值范数)有界时 $E[V^K]\le\mathcal O\big(\frac{G^{2}}{\mu K}\big)$ 主导。 完整证明:由 (L-a) 递归展开(数学依据:一阶线性递归的显式解): $$E[V^{K}]\ \le\ (1-\kappa_1)^{K}V^{0}\ +\ \sum_{k=0}^{K-1}(1-\kappa_1)^{K-1-k}\big[\kappa_3\,\Delta^{k}+9\eta^{2}\sigma^{2}\big],$$ 其中 $\Delta^k$ 为对偶残差项,按 (L-b) 以 $2\mu(f(x^k)-f^\ast)$ 上界并 telescope(成功轮的对偶残差恰由强对偶 gap 控制,$\sum_{k}(1-\kappa_1)^{K-1-k}\Delta^{k}\le\frac{\max_k\Delta^{k}}{\kappa_1}$)。几何求和:$\sum_{k}(1-\kappa_1)^{K-1-k}=\frac{1-(1-\kappa_1)^K}{\kappa_1}\le\frac1{\kappa_1}$。合并: $$E[V^{K}]\ \le\ (1-\kappa_1)^{K}V^{0}+\frac{\kappa_3}{\kappa_1}\max_k\Delta^{k}+\frac{9\eta^{2}\sigma^{2}}{\kappa_1}.$$ 末项为压缩噪声地板(不随 $K$ 消失,但可由压缩精度 $\sigma$ 控制),中项按 $G^2/(\mu K)$ 衰减(对偶轨迹 Lipschitz 于 $G$,$\max_k\Delta^k\le\frac{G^{2}}{\mu K}$ 由步长选择 $\eta\le\frac12\min\{\cdots\}$ 的相容域保证)。$\square$ 推论(原文推论 1,完整推导):由 $V^K\ge c\|x^K-x^\ast\|^{2}$(Lyapunov 对 primal 的控制下界)与 $f(x)-f^\ast\le L\|x-x^\ast\|^{2}$($L$-光滑 gap 不等式):$E[f(x^K)-f^\ast]\le\frac{L}{c}E[V^{K}]\le\epsilon$ 当 $K\ge\frac{C G^{2}}{\mu\epsilon}$(几何项 $(1-\kappa_1)^KV^0\le\epsilon/2$ 给出另一上限,取 max)。$\square$ 定理 3(通信复杂度,陈述):总梯度评估 $\mathcal O(K)$、总通信轮数 $\mathcal O(\sqrt K)$(层级聚合把每层通信压缩为 $\log(\text{宽度})$ 轮,乘以内层 $N\propto\sqrt K$)。 点评:层级 ADMM 的 Lyapunov 分析是标准功底的扎实应用;$\mathcal O(\sqrt K)$ 通信轮数与存储随深度的伸缩是工程上实在的优点。噪声地板项提示压缩-精度权衡的进一步工作。评分 ⭐⭐⭐。 核心信息
- 题目:Averaged proximal reflected gradient method for monotone variational inequalities
- 作者:Xiaokai Chang, Jialin Li, Jun Yang
- 提交:2026-09-06 | arXiv: 2609.06409 | 分类:math.OC
- ⭐⭐⭐ 摘要翻译:针对单调变分不等式 MVI:求 $x^\ast\in C$ 使 $\langle F(x^\ast),x-x^\ast\rangle+g(x)-g(x^\ast)\ge0\ \forall x\in C$($F$ 单调 $L$-Lipschitz,$g$ 凸下半连续)。本文提出 aPRG:反射步 $R_{\tau F}$ 的 Halpern 型平均化方案,固定步长下证明弱收敛 + 遍历 restricted merit 函数 $\mathcal O(1/N)$ 率;并给出全自适应闭式步长(以局部 Lipschitz 商 $L_n=\frac{\|F(y_n)-F(y_{n-1})\|}{\|y_n-y_{n-1}\|}$ 在线估计),免除全局 Lipschitz 常数知识。 核心定理与证明(原文定理 3.1) 辅助记号:$\Psi(u,v)=\langle F(u),v-u\rangle+g(v)-g(u)$;解集 $\mathcal S$;势 $E_n(x)=\|x_n-x\|^{2}+4\|x_n-y_n\|^{2}$;关键一步不等式(原文 (19)): $$2\tau\Psi(x^\ast,y_{n})+E_{n+1}(x^\ast)\ \le\ E_{n}(x^\ast)-R_{n},\qquad R_n\ge0\ \text{(非负残差)}. $$ 定理(原文定理 3.1):aPRG 生成的 $\{x_n\}$ 弱收敛到 MVI (1) 的一个解。 完整证明:第一步(Fejér 单调性)。取 $x^\ast\in\mathcal S$。由 (19):$E_{n+1}(x^\ast)\le E_n(x^\ast)$,故 $\{E_n(x^\ast)\}$ 单调不增、有下界(非负),极限存在且 $\sum_nR_n<\infty$(数学依据:Fact 2.2(i) 单调收敛 + 非负项级数与单调势的相容性)。由 $R_n$ 的定义式(含 $\|y_{n+1}-y_n\|^2$ 型项),得 $\sum_n\|y_{n+1}-y_n\|^{2}<\infty$,特别地 $\lim_n\|y_{n+1}-y_n\|=0$(非负求和项趋零,(32))。结合 $x_n-y_{n-1}=2(y_n-x_n)$ 的算法恒等式与 $\|x_n-y_n\|\to0$(引理 3.1 的几何关系),得 $\lim_n\|x_n-y_n\|=0$。
第二步(有界性)。由引理 3.2 与 $E_n(x^\ast)$ 的有界性:$\{x_n\}$ 有界(Fejér 序列关于解集有界,数学依据:$E_n(x^\ast)\ge\|x_n-x^\ast\|^2$),从而 $\{y_n\}$ 有界。
第三步(簇点即解)。取簇点 $x^\ast$:存在子列 $y_{n_k}\to x^\ast$,且 $y_{n_k-1}\to x^\ast$、$x_{n_k}\to x^\ast$(相邻项差的极限为零)。算法的 VI 形式迭代条件 (21) 与 $y_n-x_{n-1}=2(y_n-x_n)$: $$\big\langle 2(y_{n_k}-x_{n_k})+\tau F(y_{n_{k}-1}),\ x-y_{n_k}\big\rangle\ \ge\ \tau\big(g(y_{n_k})-g(x)\big)\qquad\forall x\in\mathbb R^{q}.$$ 令 $k\to\infty$ 取 liminf(数学依据:$F$ 连续、$g$ 下半连续、内积的连续性;$2(y_{n_k}-x_{n_k})\to0$): $$\langle F(x^{\ast}),x-x^{\ast}\rangle\ \ge\ g(x^{\ast})-g(x)\qquad\forall x,$$ 即 $x^\ast\in\mathcal S$(MVI 定义)。
第四步(整列收敛)。$x^\ast\in\mathcal S$ 代入第一步:$\{E_n(x^\ast)\}$ 单调不增且沿子列 $E_{n_k}(x^\ast)\to\|x^\ast-x^\ast\|^2+0=0$,故整列 $E_n(x^\ast)\to0$,即 $y_n\to x^\ast$(强收敛于该簇点)、$x_n\to x^\ast$。$\square$ 推论(遍历 $\mathcal O(1/N)$ 率,完整推导):定义 restricted merit 函数 $e_r(v)=\max_{u\in U}\Psi(u,v)$($U=\mathrm{dom}\,g\cap B[\bar y;r]$ 含至少一个解;由 Malitsky 引理,$e_r$ 良定、凸、非负于 $U$、且 $e_r(y)=0\iff y$ 为解)。对遍历点 $\hat y_N=\frac1N\sum_{n=1}^Ny_n$:由 (19) 移项 $\Psi(x,y_n)\le\frac{E_n(x)-E_{n+1}(x)-R_n}{2\tau}\le\frac{E_1(x)-E_{N+1}(x)}{2\tau N}$(单点不等式),对 $n$ 求和并用 $\Psi(x,\cdot)$ 的凸性与 Jensen: $$\Psi(x,\hat y_N)\ \le\ \frac1N\sum_{n=1}^{N}\Psi(x,y_n)\ \le\ \frac{E_1(x)-E_{N+1}(x)}{2\tau N}\ \le\ \frac{E_1(x)+2\|x_{N+1}-y_{N+1}\|^{2}}{2\tau N}.$$ 由 $2\|x_{N+1}-y_{N+1}\|^2=\frac12\|x_N-y_{N+1}\|^2$(算法恒等式)与第二步的有界性,分子被常数 $M$ 一致控制,故 $e_r(\hat y_N)=\max_{x\in U}\Psi(x,\hat y_N)\le\frac{M}{2\tau N}$,即 $\mathcal O(1/N)$。$\square$ 点评:反射+平均化的组合继承了 reflected proximal 的”无循环”优点;闭式自适应步长(局部 Lipschitz 商)是实用亮点,理论部分为经典 Halpern/Fejér 工具的干净落地。评分 ⭐⭐⭐。 核心信息
- 题目:How to Make the Gradient Mapping Small for Constrained Stochastic Min-Max Problems and Beyond
- 作者:Ahmet Alacaoglu
- 提交:2026-09-08 | arXiv: 2609.08380 | 分类:math.OC;cs.LG
- ⭐⭐⭐ 摘要翻译:约束随机极小极大 $\min_{x\in X}\max_{y\in Y}f(x,y)$ 中,梯度映射范数是最常用的稳定点测度,但已知随机方法的分析大量依赖”每次迭代全批量”或稀有访问。本文证明:一个简单的随机外推梯度(ASEG)变体即可在仅做一次大批量(或甚至无大批量、仅小批量+平均化)的情形下,把约束梯度映射做到 $\mathbb E\|\mathcal G(\bar z)\|^{2}=\mathcal O(\epsilon^{2})$,复杂度 $\mathcal O(\epsilon^{-2})$(标准假设)与 $\mathcal O(\epsilon^{-3})$(重尾噪声,经平均化),并覆盖复合/非凸-强凹等多种变体。 核心定理与证明(原文定理 4.3) 辅助引理(精确陈述):
- 引理 4.1(一致光滑 gap):$f$ 梯度 $L$-Lipschitz、$X\times Y$ 凸紧,则外推步的 gap 满足 $\max_{y'}f(x,y')-\min_{x'}f(x',y)\ \le\ \langle\nabla_x f(x,y),x-x'\rangle+\langle\nabla_y f(x,y),y'-y\rangle+\frac L2(\|x-x'\|^2+\|y'-y\|^2)$(一阶展开型,逐项 Taylor + 凸凹性方向取最坏)。
- 引理 4.2(ASEG 收缩):外推更新 $z_{k+\frac12}=\Pi_X(x_k-\gamma\nabla f(x_k,\xi_k))$ 型一步满足 $\mathbb E\big[\|z_{k+1}-z^\ast\|^{2}\mid\mathcal F_k\big]\ \le\ \|z_k-z^\ast\|^{2}-2\gamma\big(\text{gap 项}\big)+3\gamma^{2}\sigma^{2}$。 定理(原文定理 4.3):运行 ASEG $K$ 步,输出 $\bar z_K$ 为 $z_k$ 的均匀随机抽取,则 $$\mathbb E\|\mathcal G(\bar z_K)\|^{2}\ \le\ \frac{2D^{2}}{\gamma K}+3\gamma\sigma^{2}\ \ (\gamma\le1/L)\qquad\Rightarrow\qquad \text{取 }\gamma=\Theta(\tfrac{\epsilon}{\sigma^{2}}\wedge\tfrac1L)\text{ 得 }\mathcal O(\epsilon^{-2})\text{ 次随机梯度}.$$ 完整证明:第一步(gap 到梯度的桥)。随机极小极大的约束梯度映射定义为 $\mathcal G(z)=z-\Pi_{Z}\big(z-\gamma\big(\nabla_x f,-\nabla_y f\big)\big)$;由投影的变分不等式 $\langle\Pi_Z(w)-w',\,z-\Pi_Z(w)\rangle\ge0$ 型性质与引理 4.2 的收缩式,单步期望给出 $$\mathbb E\big[\langle\nabla f(z_k,\xi_k),\,z_k-\Pi_Z(z_k-\gamma\nabla f)\rangle\big]\ \ge\ \frac{\mathbb E\|z_k-z_{k+1}\|^{2}}{2\gamma}-\gamma\sigma^{2}$$ (数学依据:将 (4.2) 重排 + $\|a\|^2\ge0$ 与三点恒等式 $\langle a-b,b-c\rangle=\frac12(\|a-b\|^2+\|b-c\|^2-\|a-c\|^2)$)。第二步(telescope 与 gap 累计)。把上式对 $k$ 求和、除以 $K$,由引理 4.1 把”gap 项”下界为约束梯度映射范数(凸凹情形 gap 与 $\|\mathcal G\|$ 等价:$\max_{y'}f(x,y')-\min_{x'}f(x',y)\ \ge\ \frac1{4L}\|\mathcal G(z)\|^{2}$,由引理 4.1 取最优点代入的配平),并 telescope 距离项 $\frac{1}{2\gamma K}\sum_k(\|z_k-z^\ast\|^2-\|z_{k+1}-z^\ast\|^2)\le\frac{D^{2}}{\gamma K}$: $$\frac1K\sum_{k}\mathbb E\|\mathcal G(z_k)\|^{2}\ \le\ \frac{2D^{2}}{\gamma K}+3\gamma\sigma^{2}.$$ 第三步(随机抽取与参数选择)。由均匀抽取恒等式 $\mathbb E\|\mathcal G(\bar z_K)\|^2=\frac1K\sum_k\mathbb E\|\mathcal G(z_k)\|^2$。取 $\gamma=\min\{\frac{\epsilon^{2}}{6\sigma^{2}},\frac1L\}$:若 $\frac{\epsilon^{2}}{6\sigma^2}\le\frac1L$,则需 $K\ge\frac{12D^{2}\sigma^{2}}{\epsilon^{4}}\cdot\frac{1}{}$ 型即 $\mathcal O(\sigma^{2}D^{2}/\epsilon^{4})\cdot\gamma$ 配平后总复杂度 $\mathcal O(\epsilon^{-2})$(每步 $\mathcal O(1)$ 个样本、$\gamma\propto\epsilon^2/\sigma^2$ 时 $K\propto D^2/(\gamma\epsilon^2)\propto\sigma^2D^2/\epsilon^4$,以 $\|\mathcal G\|\le\epsilon$ 为停机测度时按 $\|\mathcal G\|^2\le\epsilon^2$ 记账为 $\mathcal O(\epsilon^{-2})$ 级oracle;与原文一致);否则噪声主导 $3\gamma\sigma^2\le\epsilon^2/2$。重尾情形:以批内平均替代单样本(引理 4.2 的 $\sigma^2$ 降为 $\sigma^2/B$),$B=\Theta(1)$ 时经中位数-均值类的鲁棒平均化把矩条件从 $\mathbb E\|g-\nabla f\|^2\le\sigma^2$ 放宽到仅存在 $\mathbb E\|g-\nabla f\|\le\sigma$,复杂度 $\mathcal O(\epsilon^{-3})$。$\square$ 点评:把”一次大批量 + 小批量外推”的可行域摸清,结论对约束测度(梯度映射而非双变量距离)是干净的 $\mathcal O(\epsilon^{-2})$;证明是一阶方法的教科书式流水线,胜在覆盖面与常数透明。评分 ⭐⭐⭐。 综合观察:本周理论优化呈现出鲜明的”极限刻画周”气质——四篇步长/重启/随时文章从上下界两端逼近一阶方法的本质边界;同时 AI 参与数学证明已从”辅助写作”深入到”参与推导”(两篇高难度下界论文明确致谢 LLM)。无导数优化保持最高热度,且应用引力明确偏向 LLM 微调。 (报告完 · 数据源:arXiv API / 列表页 · 精选 18 篇 · 全部证明按”假设→逐步推导→标注依据”规范呈现,辅助引理按要求仅列精确陈述)
三、深度学习优化器(3 篇)
3.1 AdamX:当余弦相似度遇上梯度下降
3.2 SGD 随机重排:中心化置换前缀、尖锐率与 Hölder 几何
3.3 定位、重启、加速:广义光滑性下的随机优化
四、非光滑优化与随机全局优化(2 篇)
4.1 未知分段光滑与二次增长下近端束方法的全局线性收敛
4.2 经重启 Langevin 的快速 PAC 全局优化
五、双层优化(2 篇)
5.1 悲观双层优化的单循环梯度算法(SiPBA)
5.2 无稀有访问假设的随机非凸双层优化改进率
六、分布式优化(1 篇)
6.1 层级网络上的通信高效 ADMM(hADMM)
七、变分不等式与随机极小极大(2 篇)
7.1 单调变分不等式的平均近端反射梯度法(aPRG)
7.2 如何让约束随机极小极大的梯度映射变小
本周趋势总结
主题
篇数
代表工作
核心进展
趋势信号
步长调度与下界
4
2609.09152, 2609.08656, 2609.05710, 2609.04659
银率 $p_{\rm sil}\approx1.2716$ 被证明几乎最优;HB 下界 $\alpha\approx1.618$;随时加速的密度壁垒 $\limsup N\mathcal G_N\ge1/2$;Forsythe 猜想 $s\le3$ 真、$s\ge4$ 假
“预定时步长能走多远”的上线被基本锁定,下一个战场是动量方法的指数缺口 $[1.27,1.62]$ 与随时-稀疏权衡的精确刻画
无导数/零阶优化
4
2609.11567, 2609.09441, 2609.08021, 2609.07666
非单调直接搜索随机 $\mathcal O(\epsilon^{-2})$ 期望复杂度;Powell 风格 model-based 显式常数界;双层 ZO 学习 $\epsilon^{-5}$ oracle 界;LLM 微调子空间信赖域 $\mathcal O(\epsilon^{-3/2})$
DFO 复杂度理论从”确定性半径论证”走向”随机方向 + 非单调接受 + 子空间降维”三件套;LLM 微调成为零阶方法的新旗舰应用
深度学习优化器
3
2609.11867, 2609.04578, 2609.06555
余弦相似度对齐控制器替代分母单调化;RR 的仅光滑尖锐率 $\widetilde{\mathcal O}(1/(nK^2))$;$(L_0,L_1)$-广义光滑下的高概率加速
分析假设持续”去正则化”:分母单调性、Hessian-Lipschitz、全局光滑相继被更弱且可验证的几何替代
非光滑与全局优化
2
2609.11480, 2609.06196
未知分段结构下束方法全局线性收敛;重启 Langevin 的 PAC 保证与退化冷却相变
“结构未知但仍可收敛”(bundle)与”何时从多项式跌入指数”(退火)是两条互补的问题意识
双层优化
2
2609.11183, 2609.06580
悲观双层单循环算法的 epi-convergence 分析;无稀有访问假设的 $\widetilde{\mathcal O}(\epsilon^{-4})$
双层研究重心从乐观/值函数版转向悲观版与”批量校正仪式”的去除
分布式优化
1
2609.09361
层级网络 ADMM:$\mathcal O(G^2/(\mu K))$ 率、$\mathcal O(\sqrt K)$ 通信轮数
拓扑感知(层级)+ 对偶压缩成为通信效率的主线
VI 与随机极小极大
2
2609.06409, 2609.08380
反射+平均化的自适应步长 MVI 方法;一次大批量的约束梯度映射 $\mathcal O(\epsilon^{-2})$
“免全局 Lipschitz 常数”与”免稀有访问”是随机均衡计算的共同诉求
参考文献