Small regularizers can preserve linear programming solutions exactly. This paper provides the first average-case analysis of exact regularization: with a standard Gaussian cost vector and fixed constraint set, bounds are established for the probability that exact regularization succeeds as a function of regularization strength. Failure is characterized via the Gaussian measure of inner cones, controlled by novel two-sided bounds on the measure of shifted cones. Results reveal dimension-dependent scaling laws and connect exact regularization of linear programs to their polyhedral geometry via the normal fan and the Gaussian (solid-angle) measure of its cones. Computable bounds are obtained in several canonical settings, including regularized optimal transport. Numerical experiments corroborate the predicted scalings and thresholds.
論文ID : 2510.13083タイトル : Average-case thresholds for exact regularization of linear programs著者 : Michael P. Friedlander, Sharvaj Kubal, Yaniv Plan, Matthew S. Scott分類 : math.OC cs.NA math.NA math.PR発表日 : 2025年10月15日論文リンク : https://arxiv.org/abs/2510.13083 小さな正則化項は線形計画法の解を正確に保存することができます。本論文は正確な正則化の初めての平均ケース分析を提供します。標準ガウス成本ベクトルと固定制約集合の下で、正則化強度に関する正確な正則化成功確率の界を確立します。失敗は内錐のガウス測度によって特徴付けられ、シフトされた錐測度の新しい両側界によって制御されます。結果は次元依存のスケーリング則を明らかにし、法錐扇とその錐のガウス(立体角)測度を通じて線形計画法の正確な正則化をその多面体幾何と結びつけます。いくつかの典型的な設定で計算可能な界が得られ、正則化最適輸送を含みます。数値実験は予測されたスケーリングと閾値を確認します。
本論文が研究する中核的な問題は正確な正則化 現象です。線形計画法問題
(P0) min { − g ⋅ x ∣ A x ≤ b } \text{(P0)} \quad \min \{-g \cdot x | Ax \leq b\} (P0) min { − g ⋅ x ∣ A x ≤ b }
およびその正則化版
(P ε ) min { − g ⋅ x + ε ψ ( x ) ∣ A x ≤ b } \text{(P}_\varepsilon\text{)} \quad \min \{-g \cdot x + \varepsilon\psi(x) | Ax \leq b\} (P ε ) min { − g ⋅ x + ε ψ ( x ) ∣ A x ≤ b }
において、正則化パラメータε \varepsilon ε が十分に小さいとき、正則化問題の解は依然として元の問題の解であり、すなわちSol ( P ε ) ⊆ Sol ( P 0 ) \text{Sol}(\text{P}_\varepsilon) \subseteq \text{Sol}(\text{P}_0) Sol ( P ε ) ⊆ Sol ( P 0 ) が成立します。
計算上の課題 :線形計画法は非一意解と極端なデータ摂動感度を持つ可能性があり、正則化はこれらの問題を緩和できます理論的空白 :既存の決定論的分析は正確な正則化閾値ε ˉ \bar{\varepsilon} ε ˉ を決定するために元の問題を先に解く必要があります実用的需要 :最適輸送などの応用では、二次正則化はエントロピー正則化よりも疎性をよく保存します幾何学的洞察 :確率分析を通じて正確な正則化と多面体幾何の深い関連性を明らかにしますMangasarianとMeyerの決定論的理論はP0と選択問題P_ψを同時に解く必要があります González-SanzとNutzの界は平均ケースではあまりに緩く(n倍改善されました) 次元依存のスケーリング則分析が不足しています シフトされた錐のガウス測度界 :定理1.3を確立し、錐V+wのガウス測度の両側界を提供します幾何学的特徴付け :正確な正則化確率がすべての頂点における内錐ガウス測度の合計に等しいことを証明します(定理1.5)内錐界のスイート :メンバーシップ条件(定理2.1)と表現ベクトル(定理2.5)を通じて包括的な界を開発します辺界分析 :面構造を通じて辺界の両側界を提供します(補題3.2、定理3.3)具体的応用 :∞-球制約と正則化最適輸送に対して最適な界を提供します(定数まで)多面体可行領域Q = { x ∈ R n ∣ A x ≤ b } Q = \{x \in \mathbb{R}^n | Ax \leq b\} Q = { x ∈ R n ∣ A x ≤ b } と正則化関数ψ \psi ψ が与えられたとき、成本ベクトルg ∼ N ( 0 , I n ) g \sim N(0, I_n) g ∼ N ( 0 , I n ) のとき、正確な正則化事象ER ( ε ) \text{ER}(\varepsilon) ER ( ε ) の確率を分析します。
x ∈ Q x \in Q x ∈ Q に対して、法錐は以下のように定義されます:
N ( x ) : = { v ∈ R n ∣ v ⋅ ( y − x ) ≤ 0 for all y ∈ Q } N(x) := \{v \in \mathbb{R}^n | v \cdot (y-x) \leq 0 \text{ for all } y \in Q\} N ( x ) := { v ∈ R n ∣ v ⋅ ( y − x ) ≤ 0 for all y ∈ Q }
主要性質(命題1.1):
N ( z ) N(z) N ( z ) がn次元であることとz ∈ Vert ( Q ) z \in \text{Vert}(Q) z ∈ Vert ( Q ) は同値です頂点における法錐はほぼ全空間を分割します:⋃ z ∈ Vert ( Q ) N ( z ) = R n \bigcup_{z \in \text{Vert}(Q)} N(z) = \mathbb{R}^n ⋃ z ∈ Vert ( Q ) N ( z ) = R n P0の最適性:g ∈ N ( z ) g \in N(z) g ∈ N ( z ) P_εの最適性:g ∈ N ( z ) + ∂ ( ε ψ ) ( z ) g \in N(z) + \partial(\varepsilon\psi)(z) g ∈ N ( z ) + ∂ ( ε ψ ) ( z ) 定理1.3 (中核的技術結果):錐V ⊆ R d V \subseteq \mathbb{R}^d V ⊆ R d とベクトルw ∈ R d w \in \mathbb{R}^d w ∈ R d に対して、
γ ( V + w ) ∈ [ γ ( V ) exp ( − 1 2 ∥ w ∥ 2 − dist ( w , − V ∗ ) d ) , γ ( V ) exp ( − 1 2 ∥ Π V ∗ w ∥ 2 + dist ( w , V ∗ ) d ) ] \gamma(V + w) \in \left[\gamma(V) \exp\left(-\frac{1}{2}\|w\|^2 - \text{dist}(w, -V^*)\sqrt{d}\right), \gamma(V) \exp\left(-\frac{1}{2}\|\Pi_{V^*}w\|^2 + \text{dist}(w, V^*)\sqrt{d}\right)\right] γ ( V + w ) ∈ [ γ ( V ) exp ( − 2 1 ∥ w ∥ 2 − dist ( w , − V ∗ ) d ) , γ ( V ) exp ( − 2 1 ∥ Π V ∗ w ∥ 2 + dist ( w , V ∗ ) d ) ]
証明技法:
上界:Hölder不等式と弱双対性を使用し、分散パラメータκ \kappa κ を最適化 下界:Jensen不等式とMoreau分解を組み合わせ ∂ ( ε ψ ) ( z ) ∩ N ( z ) ≠ ∅ \partial(\varepsilon\psi)(z) \cap N(z) \neq \emptyset ∂ ( ε ψ ) ( z ) ∩ N ( z ) = ∅ のとき、内錐はN ( z ) + w N(z) + w N ( z ) + w に簡略化され、定理1.3を直接適用します。
メンバーシップ条件を満たさない場合、表現ベクトルw ~ = S T − 1 ( S T w ) + \tilde{w} = S_T^{-1}(S_T w)_+ w ~ = S T − 1 ( S T w ) + を構築し、以下を満たします:
M ( T , w ) = M ( T , w ~ ) M(T, w) = M(T, \tilde{w}) M ( T , w ) = M ( T , w ~ )
ここでM ( T , w ) = T ∖ ( T + w ) M(T, w) = T \setminus (T + w) M ( T , w ) = T ∖ ( T + w ) は辺界です。
補題3.1 :辺界は以下のように分解できます
γ ( M ( T , w ) ) ≤ ∑ i = 1 ℓ γ ( P i ) \gamma(M(T, w)) \leq \sum_{i=1}^\ell \gamma(P^i) γ ( M ( T , w )) ≤ ∑ i = 1 ℓ γ ( P i )
ここでP i = ⋃ t ∈ [ 0 , 1 ) [ T i + t w ] P^i = \bigcup_{t \in [0,1)} [T^i + tw] P i = ⋃ t ∈ [ 0 , 1 ) [ T i + tw ] (s i ⋅ w > 0 s_i \cdot w > 0 s i ⋅ w > 0 のとき)
Birkhoff多面体 (二重確率行列)上の二次正則化∞-球 上の二次および線形正則化次元範囲:n ∈ [ 10 3 , 10 4 ] n \in [10^3, 10^4] n ∈ [ 1 0 3 , 1 0 4 ] 各( n , ε ) (n, \varepsilon) ( n , ε ) ペアに対して20回の独立試行 正確な正則化失敗確率P ( ER c ( ε ) ) P(\text{ER}^c(\varepsilon)) P ( ER c ( ε )) 期待正則化閾値E [ ε ˉ ] E[\bar{\varepsilon}] E [ ε ˉ ] 理論界と経験的観察の比較 理論的予測:
失敗確率界:P ( ER c ( ε ) ) ≤ 1 2 ε 2 n + ε n 3 / 4 P(\text{ER}^c(\varepsilon)) \leq \frac{1}{2}\varepsilon^2\sqrt{n} + \varepsilon n^{3/4} P ( ER c ( ε )) ≤ 2 1 ε 2 n + ε n 3/4 期待閾値下界:E [ ε ˉ ] ≥ 1 − e − 4 n 2 n 3 / 4 E[\bar{\varepsilon}] \geq \frac{1-e^{-4n}}{2n^{3/4}} E [ ε ˉ ] ≥ 2 n 3/4 1 − e − 4 n 実験検証:
水平曲線ε n 3 / 4 = 2 \varepsilon n^{3/4} = 2 ε n 3/4 = 2 での経験的失敗確率は約0.5で、理論界0.875と一致 スケーリング関係は対数スケールで直線を示し、次元依存性を確認 二次正則化 :
理論界:P ( ER c ( ε ) ) ≤ ε n + 1 2 ε 2 n P(\text{ER}^c(\varepsilon)) \leq \varepsilon n + \frac{1}{2}\varepsilon^2 n P ( ER c ( ε )) ≤ ε n + 2 1 ε 2 n ε < n − 1 \varepsilon < n^{-1} ε < n − 1 のとき、第一項が支配的で、界は定数因子2 π e \sqrt{2\pi e} 2 π e 内で最適線形正則化 :
辺界手法はより厳密な界を与えます:P ( ER c ( ε ) ) ≤ 1 2 π ε ∥ p ∥ 1 exp ( ε n − 1 ∥ p ∥ ) P(\text{ER}^c(\varepsilon)) \leq \frac{1}{\sqrt{2\pi}}\varepsilon\|p\|_1 \exp(\varepsilon\sqrt{n-1}\|p\|) P ( ER c ( ε )) ≤ 2 π 1 ε ∥ p ∥ 1 exp ( ε n − 1 ∥ p ∥ ) ほぼ疎なベクトルp p p (比率∥ p ∥ 1 / ( n ∥ p ∥ ) \|p\|_1/(\sqrt{n}\|p\|) ∥ p ∥ 1 / ( n ∥ p ∥ ) で定量化)に対してより効果的 命題4.1 は∞-球上の下界を証明します:
P ( ER c ( ε ) ) ≥ 1 − exp ( − 2 π e n ε ) P(\text{ER}^c(\varepsilon)) \geq 1 - \exp\left(-\sqrt{\frac{2}{\pi e}}n\varepsilon\right) P ( ER c ( ε )) ≥ 1 − exp ( − π e 2 n ε ) ε ≤ ( π e ) / 2 ⋅ 1 2 n \varepsilon \leq \sqrt{(\pi e)/2} \cdot \frac{1}{2n} ε ≤ ( π e ) /2 ⋅ 2 n 1 に対して、P ( ER c ( ε ) ) ≥ 1 2 π e n ε P(\text{ER}^c(\varepsilon)) \geq \frac{1}{\sqrt{2\pi e}}n\varepsilon P ( ER c ( ε )) ≥ 2 π e 1 n ε
Mangasarian and Meyer 1979 :基本条件の確立 Friedlander and Tseng 2008 :一般凸計画法の特徴付け 双対選択問題P ψ \text{P}_\psi P ψ を通じた閾値ε ˉ \bar{\varepsilon} ε ˉ の特徴付け Cuturi 2013 :エントロピー正則化のSinkhorアルゴリズム González-Sanz and Nutz 2024 :二次正則化の正確性 本論文は彼らの界E [ ε ˉ ] ≥ c / ( 4 n 2 ) E[\bar{\varepsilon}] \geq c/(4n^2) E [ ε ˉ ] ≥ c / ( 4 n 2 ) をE [ ε ˉ ] ≥ 1 / ( 4 n ) E[\bar{\varepsilon}] \geq 1/(4n) E [ ε ˉ ] ≥ 1/ ( 4 n ) に改善 次元スケーリング則 :正確な正則化閾値は明確な次元依存性を持ちます。例えば∝ n − 3 / 4 \propto n^{-3/4} ∝ n − 3/4 (Birkhoff多面体)および∝ n − 1 \propto n^{-1} ∝ n − 1 (∞-球)幾何学的関連性 :法錐扇のガウス測度を通じて正確な正則化と多面体幾何の深い関連性を確立ソフト相転移 :明確な相転移閾値が存在し、小さいε \varepsilon ε では高確率で成功、大きいε \varepsilon ε では高確率で失敗多面体制限 :分析は多面体可行領域に限定されますガウス仮定 :成本ベクトルはガウス分布である必要があります計算複雑性 :いくつかの界はすべての頂点の法錐を計算する必要があります解距離界 :正確な正則化が失敗するとき、dist ( x ε , Sol ( P 0 ) ) \text{dist}(x_\varepsilon, \text{Sol}(\text{P}_0)) dist ( x ε , Sol ( P 0 )) の確率界支持単調性 :二次制約正則化LPにおける支持単調性の確率の定量化一般錐計画法 :二次および半定値制約への拡張理論的革新 :正確な正則化の初めての平均ケース分析で、重要な理論的空白を埋めます技術的深さ :シフトされた錐ガウス測度の両側界は中核的技術貢献です実用的価値 :正則化最適輸送などの応用に計算可能な界を提供幾何学的洞察 :正確な正則化と多面体幾何の深い関連性を明らかに実験検証 :数値実験は理論的予測を十分に検証適用範囲 :多面体制約とガウス成本ベクトルに限定定数最適化 :いくつかの界の定数は十分に厳密でない可能性計算複雑性 :実際の応用ではすべての頂点法錐の計算が困難な場合がある理論的貢献 :確率的線形計画法理論に新しいツールを提供応用価値 :最適輸送と疎最適化に実用的な指導を提供方法論 :シフトされた錐分析技術は他の問題に一般化可能正則化効果の理解が必要な線形計画法応用 最適輸送における正則化パラメータの選択 高次元線形計画法のロバスト性分析 疎最適化における正確な回復保証 本論文は36篇の関連文献を引用し、主に以下を含みます:
古典的凸最適化理論(Rockafellar, Hiriart-Urruty & Lemaréchal) 正確正則化理論(Mangasarian & Meyer, Friedlander & Tseng) 最適輸送(Cuturi, González-Sanz & Nutz) 高次元確率(Vershynin, Ledoux & Talagrand)