精读笔记

Problem Setting

这篇论文实际解决的是 MPCC 正则化方法中的“表述退化”问题，而不是单纯提出另一个 relaxation。MPCC 的难点在于 complementarity 约束让标准 NLP 约束资格天然失效，导致 KKT 理论、乘子稳定性和二阶结构都不能直接套用。已有正则化路线通常通过放松 complementarity 来得到可由 NLP solver 处理的近似问题，但核心矛盾是：放松后的可行集是否既足够接近原 MPCC，又保留一个稳定的临界点结构。

Scholtes 正则化数值强，但理论收敛到 S-stationarity 通常要额外 SONC；Kanzow-Schwartz 正则化理论上更强，在 MPCC-LICQ 和 ND2 下可无 SONC 得到 S-stationarity，但数值表现并不相称。本文的判断是：KS 的问题不在可行集，而在函数描述。它把一个本质上带 disjunctive corner 的集合硬写成 NLP inequality，导致在关键区域梯度消失、LICQ 失败、乘子可能灾难性放大。本文真正要解决的是：如何保留 KS 可行集的收敛优势，同时换一种不会破坏约束资格和拓扑结构的描述。

Motivation

作者的核心观察非常直接：KS(t) 的可行集和某个 OR/max-type 约束描述的可行集完全相同，但 KS 的 φ/Φ 函数描述在 re-entering corner 和 H12 点处制造了不必要的退化。也就是说，理论上有优势的不是 KS 的光滑函数 Φ，而是它定义的几何集合；数值上拖后腿的恰恰是这个函数表达。

已有路线缺的是一种“几何正确”的正则化表达。Scholtes 通过 product ≤ t 形成平滑膨胀区域，但会改变临界点 index 的组织方式；KS 给出更紧的正则化集合，但把 disjunctive geometry 错当作 NLP geometry。本文因此选择不改集合、只改语言：把同一个集合显式写成 t-F1 ≥ 0 or t-F2 ≥ 0，或等价 max{t-F1,t-F2} ≥ 0。这个方向的关键不是更平滑，而是更忠实于集合的局部拓扑。

Core Idea

核心思想是：正则化的本质不是让 MPCC 变成一个“看起来光滑”的 NLP，而是让近似问题的约束表示与其真实局部几何一致。D(t) 与 KS(t) 有相同可行集，但 D(t) 用 disjunctive constraints 表达：对每个 complementarity pair，要求至少一个分量不超过 t，同时保持非负性。这使正则化集合仍然排除双正的大块区域，但避免了 KS 在 F1=F2=t 处由 Φ 梯度为零造成的退化。

这个改变引入的 inductive bias 是“把 biactive/branch 结构显式保留下来”。在 D(t) 中，活跃约束被分解为 H1、H2、H12、N1、N2 等结构化集合，乘子符号和活跃分支可以直接对应到 MPCC 的 a01、a10、a00 及其乘子符号。这样做的理论后果是 MPCC-LICQ 能局部继承为 DISJ-LICQ，C/M/S stationarity 的极限可以被控制，C-index 也能被追踪。和 prior 的本质区别不是可行域更好，而是同一可行域的法向锥表示更正确。

Method

第一，D(t) 用 max/OR 约束替代 KS 的 Φ≤0。它解决的是 KS 函数描述导致的梯度消失和 LICQ 失效问题。核心变化是：同一可行集的边界法向不再通过一个退化的复合函数给出，而是通过 disjunctive branch 的活动约束给出。

第二，作者建立 MPCC-LICQ 到 DISJ-LICQ 的继承。这个步骤不是技术装饰，而是整个方法成立的入口：没有约束资格继承，乘子唯一性、stationarity 极限、非退化临界点追踪都不稳。其核心机制是，靠近 MPCC 可行点且 t 足够小时，D(t) 的活跃约束索引只能落在 MPCC 对应的 active complementarity 集合中，因此线性无关性可由 MPCC-LICQ 传递。

第三，证明 D(t) 的 C/M-stationary points 收敛到 MPCC 的 C/M-stationary points，并在 MPCC-ND2 下得到 S-stationarity。这里关键是重命名 disjunctive 乘子为 MPCC 乘子，并利用活跃集合互斥关系控制乘子符号。相比 Scholtes，它避免了依赖 SONC 才到 S-stationarity 的问题。

第四，追踪 C-index。作者不仅关心收敛到什么 stationarity，还关心临界点拓扑类型是否保持。结论是，在 MPCC-ND4 下，quadratic index 和 biactive index 都对应；若 ND4 失败，quadratic index 可出现至多由 vanishing non-biactive multipliers 数量控制的 shift。这个部分解释了为什么 D(t) 的临界点结构比 Scholtes 更稳定。

Key Insight / Why It Works

最重要的 insight 是：MPCC 正则化的成败很大程度上取决于“法向结构是否被正确编码”，而不是可行集是否更平滑。KS 的可行集并不差，甚至理论上比 Scholtes 更接近 MPCC 的强 stationarity；问题是它的 NLP 描述在关键角点把法向信息压塌了。D(t) 通过 disjunctive 表述恢复了多个分支的法向，使乘子不需要通过接近零的 Φ 梯度或差分项间接表达，从而避免了 KS 中 ε-stationarity 的 unbounded multiplier 和 catastrophic cancellation。

真正有效的部分是 functional description 的改变带来的 constraint qualification inheritance 和 multiplier/sign structure alignment。C-index 分析是更深层的理论贡献：它说明 D(t) 不是仅仅“数值上好一点”，而是在拓扑类型层面以更自然的方式逼近 MPCC。特别是对非退化 S-stationary/minimizer，D(t) 的逼近点通常避开 H12 re-entering corner，局部表现为普通 KKT 点；这很好地解释了为什么它可能比 KS 数值稳定。

相对而言，实验中的 solver 选择和 continuation 参数更多是辅助 evidence，不是核心贡献。LogMIP、SLSQP、IPOPT 的差异说明 formulation 与 solver 的相互作用很强；因此数值增益不应被解释为“一个通用算法全面更强”，而应解释为“正确几何表述减少了 solver 面对的退化”。这不是 scaling、retrieval 或 data coverage 型增益，而是 better inductive bias / latent structure alignment：把 complementarity 的分支结构显式对齐到求解问题的约束表示中。

Relation To Prior Work

最接近的是 Scholtes regularization、Kanzow-Schwartz regularization，以及 disjunctive optimization 的临界点理论。本文与 KS 的关系最微妙：D(t) 和 KS(t) 可行集相同，stationarity 上 KKT(KS) 等价于 S-stationarity(D)，因此它不是通过构造一个全新 relaxation set 获得优势，而是通过换掉描述方式获得优势。看似小改动，实质上改变了约束资格和乘子稳定性。

相对 Scholtes，D(t) 的本质差异在于 index 传递机制。Scholtes 的 KKT quadratic index 会在极限中拆分为 MPCC 的 quadratic 和 biactive parts，因此要得到 S-stationarity 往往需要额外二阶条件；D(t) 中 biactive index 与 MPCC biactive index generically 对应，quadratic index 在 ND4 下也对应。这是实质创新，不是已有正则化的简单重组。

相对 sign-type regularization，作者强调 D(t) 不表现出同样稳定 index shift 和 bifurcation 现象。相对 continuous quadrant penalty formulation，论文指出即便可行集同样可表达，若仍作为 NLP equality/penalty 表述，会在 H12 处梯度为零，再次失去 LICQ。因此本文属于“用 disjunctive critical point theory 修正 MPCC regularization 表述”的谱系，而不是传统 smooth relaxation 谱系。

Dataset / Evaluation

实验覆盖 75 个 MacMPEC 小中规模实例加文中构造例，比较 KS、D、Scholtes，并用 SLSQP、IPOPT 以及 LogMIP 做不同求解器视角的验证。评价重点包括 normalized relative objective error、maximum constraint violation 和 CPU time。总体上实验支持两个 claim：D(t) 明显优于 KS(t)；在高精度要求下，D(t) 可与甚至优于 Scholtes。

但 evaluation 仍应看作 preliminary。MacMPEC 是标准库，但规模和现代实际应用中的大规模 MPCC、动态 MPCC、nonsmooth optimal control 仍有距离。实验没有充分隔离增益来源：是 DISJ-LICQ 改善、避免 H12 退化、solver 对 max/OR formulation 更友好，还是 continuation 策略与参数选择带来的偏置，文中未完全分解。LogMIP 部分有意义，因为它验证了 disjunctive formulation 可被专门 solver 利用；但其设置为直接取 t=t_min，与 continuation NLP solver 并不完全等价。总体实验能支持“这个 formulation 值得重视”，但还不足以证明其在大规模真实部署中系统优于 Scholtes。

Limitation

第一，理论依赖明显。MPCC-LICQ、MPCC-ND2、MPCC-ND4 是很多关键结论的支撑；虽然作者强调这些条件 generic，但实际建模问题中退化并不罕见，尤其来自 KKT 嵌套、bilevel 或 equilibrium 系统的 MPCC 往往有结构性依赖。ND4 失败时 index shift 不可避免，well-posedness 也不再干净。

第二，D(t) 没有消除 disjunctive optimization 的组合性。它只是承认这个组合结构并用合适语言表达。对于通用 NLP solver，max-type/non-smooth/disjunctive 约束如何稳定处理仍是问题；对于 LogMIP 或混合逻辑求解器，scale-up 可能遇到组合爆炸。也就是说，本文把 KS 的“错误光滑化”问题转移为“正确但可能更组合”的建模问题。

第三，数值 claim 仍有归因不清。D(t) 的优势很可能来自 formulation-level conditioning 和 multiplier stability，而非新的搜索机制。对于不同 solver、不同初始化、不同 t update 策略，优势是否稳定，文中未充分说明。特别是高精度优于 Scholtes的结论很有吸引力，但是否普遍成立仍需要更大规模和更细 ablation。

第四，ε-stationarity 的讨论主要用于批判 KS 的不稳定性，但 D(t) 自身在 inexact solve 下的鲁棒 stationarity 理论并没有同等完整展开。实际 solver 总是近似求解，如果 D(t) 的近似 stationarity 如何传递到 MPCC stationarity 没有系统理论，部署时仍存在缺口。

Takeaway

• 1. MPCC 正则化中，等价可行集不等价于等价算法；函数描述决定法向锥、乘子稳定性和临界点拓扑。
• 2. KS 的理论优势来自其正则化集合，数值劣势来自把 disjunctive geometry 误写成退化 NLP。
• D(t) 的价值在于把同一个集合重新表述为正确的 OR/max structure。
• 3. C-index 追踪是本文最值得迁移的分析工具：它把“收敛到什么 stationarity”提升为“临界点类型如何随正则化保持或改变”。

一句话总结

这篇论文在 MPCC 正则化谱系中的贡献是：证明 Kanzow-Schwartz 可行集应被视为 disjunctive object 而非 NLP object，并通过正确的 OR/max 表述同时获得约束资格继承、强 stationarity 收敛和更稳定的临界点拓扑。