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.
معرّف الورقة : 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تاريخ النشر : 15 أكتوبر 2025رابط الورقة : 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} ε ˉ الاحتياجات العملية : في تطبيقات مثل النقل الأمثل، يحافظ التنظيم التربيعي على الندرة بشكل أفضل من التنظيم الإنتروبيالرؤى الهندسية : الكشف من خلال التحليل الاحتمالي عن الارتباط العميق بين التنظيم الدقيق والهندسة متعددة الأوجهالنظرية الحتمية لمانجاساريان وماير تتطلب حل P0 واختيار المشكلة P_ψ في نفس الوقت الحدود من غونزاليز-سانز ونوتز فضفاضة جداً في الحالة المتوسطة (تحسن بمعامل n) نقص تحليل قوانين التحجيم المعتمدة على البعد حدود مقياس غوسي المخروط المزاح : إنشاء Theorem 1.3، توفير حدود ثنائية الاتجاه لمقياس غوسي المخروط V+wالتوصيف الهندسي : إثبات أن احتمالية التنظيم الدقيق تساوي مجموع مقاييس غوسي المخروط الداخلي في جميع الرؤوس (Theorem 1.5)مجموعة حدود المخروط الداخلي : تطوير حدود شاملة من خلال شروط العضوية (Theorem 2.1) وناقلات التمثيل (Theorem 2.5)تحليل المجموعات الحدية : توفير حدود ثنائية الاتجاه للمجموعات الحدية من خلال بنية الوجوه (Lemma 3.2، Theorem 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 ψ ، تحليل احتمالية حدث التنظيم الدقيق ER ( ε ) \text{ER}(\varepsilon) ER ( ε ) عندما يكون متجه التكلفة g ∼ N ( 0 , I n ) g \sim N(0, I_n) g ∼ N ( 0 , I n ) .
بالنسبة لـ 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 }
الخصائص الرئيسية (Proposition 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 ) Theorem 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 ) ]
تقنيات الإثبات:
الحد الأعلى: استخدام عدم المساواة هولدر والثنائية الضعيفة، من خلال تحسين معامل التباين κ \kappa κ الحد الأدنى: عدم المساواة جنسن مع تحليل موريو عندما يكون ∂ ( ε ψ ) ( z ) ∩ N ( z ) ≠ ∅ \partial(\varepsilon\psi)(z) \cap N(z) \neq \emptyset ∂ ( ε ψ ) ( z ) ∩ N ( z ) = ∅ ، يتم تبسيط المخروط الداخلي إلى N ( z ) + w N(z) + w N ( z ) + w ، مع تطبيق مباشر لـ Theorem 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 ) هي المجموعة الحدية.
Lemma 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 ).
متعدد برخوف (المصفوفات العشوائية المزدوجة) مع التنظيم التربيعيكرة ∞ مع التنظيم التربيعي والخطينطاق البعد: n ∈ [ 10 3 , 10 4 ] n \in [10^3, 10^4] n ∈ [ 1 0 3 , 1 0 4 ] 20 تجربة مستقلة لكل زوج ( n , ε ) (n, \varepsilon) ( n , ε ) احتمالية فشل التنظيم الدقيق 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 ∥ ) ) Proposition 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 ε .
مانجاساريان وماير 1979 : إنشاء الشروط الأساسية فريدلاندر وتسينج 2008 : توصيف البرامج المحدبة العامة توصيف العتبة ε ˉ \bar{\varepsilon} ε ˉ من خلال مشكلة الاختيار الثنائي P ψ \text{P}_\psi P ψ كوتوري 2013 : خوارزمية سينكهورن للتنظيم الإنتروبي غونزاليز-سانز ونوتز 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 (متعدد برخوف) و ∝ n − 1 \propto n^{-1} ∝ n − 1 (كرة ∞)الارتباط الهندسي : إنشاء ارتباط عميق بين التنظيم الدقيق والهندسة متعددة الأوجه من خلال مقياس غوسي مروحة المخروط الطبيعيانتقال طور ناعم : وجود عتبة انتقال طور واضحة، احتمالية نجاح عالية عند ε صغير، احتمالية فشل عالية عند ε كبيرتقييد متعدد الأوجه : التحليل مقتصر على المناطق المجدية متعددة الأوجهافتراض غوسي : يجب أن يكون متجه التكلفة موزعاً بشكل غوسيالتعقيد الحسابي : قد تتطلب بعض الحدود حساب المخاريط الطبيعية لجميع الرؤوسحدود مسافة الحل : حدود احتمالية على 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)