Difference between revisions of "Main Page"

From MWiki
Jump to: navigation, search
(Theorem of the month)
(Centre Method)
(44 intermediate revisions by the same user not shown)
Line 1: Line 1:
 +
__NOTOC__
 
= Welcome to MWiki =
 
= Welcome to MWiki =
 
== Theorem of the month ==
 
== Theorem of the month ==
Theorem: The intex method solves every solvable LP in <math>\mathcal{O}({\vartheta}^{3})</math>.
+
The centre method solves every solvable LP in <math>\mathcal{O}(\omega{\vartheta}^{2})</math>.
  
Proof and algorithm: First, we normalise and scale <math>{b}^{T}y - {d}^{T}x \le 0, Ax \le b</math> and <math>{A}^{T}y \ge d</math>. Let the ''height'' <math>h</math> have the initial value <math>{h}_{0} := |\text{min } \{{b}_{1}, ..., {b}_{m}, {-d}_{1}, ..., {-d}_{n}\}|/r</math> for the reduction factor <math>r \in \; [&frac12;, 1[</math>. Let the
+
== Proof and algorithm ==
 +
Let <math>z := \grave{m} + n</math> and <math>d \in [0, 1]</math> the density of <math>A</math>. First, normalise and scale <math>{b}^{T}y - {c}^{T}x \le 0, Ax \le b</math> as well as <math>{A}^{T}y \ge c</math>. Let <math>P_r := \{(x, y)^T \in {}^{\omega}\mathbb{R}_{\ge 0}^{z} : {b}^{T}y - {c}^{T}x \le r \in [0, \check{r}], Ax - b \le \underline{r}_m, c - {A}^{T}y \le \underline{r}_n\}</math> have the radius <math>\check{r} := s|\min \; \{b_1, ..., b_m, -c_1, ..., -c_n\}|</math> and the scaling factor <math>s \in [1, 2]</math>. It follows <math>\underline{0}_{z} \in \partial P_{\check{r}}</math>. By the strong duality theorem, the LP min <math>\{ r \in [0, \check{r}] : (x, y)^T \in P_r\}</math> solves the LPs max <math>\{{c}^{T}x : c \in {}^{\omega}\mathbb{R}^{n}, x \in {P}_{\ge 0}\}</math> and min <math>\{{b}^{T}y : y \in {}^{\omega}\mathbb{R}_{\ge 0}^{m}, {A}^{T}y \ge c\}</math>.
  
LP min <math>\{h \in [0, {h}_{0}] : x \in {}^{\omega}\mathbb{R}_{\ge 0}^{n}, y \in {}^{\omega}\mathbb{R}_{\ge 0}^{m}, {b}^{T}y - {d}^{T}x \le h, Ax - b \le (h, ..., h)^{T} \in {}^{\omega}\mathbb{R}_{\ge 0}^{m}, d - {A}^{T}y \le (h, ..., h)^{T} \in {}^{\omega}\mathbb{R}_{\ge 0}^{n}\}</math> have for <math>\underline{v} := {v}^{T}</math> the feasible interior starting point <math>v := ({\underline{x}, \underline{y}, h)}^{T} \in {}^{\omega}\mathbb{R}_{\ge 0}^{m+n+1}</math>, e.g. <math>({\underline{0}, \underline{0}, {h}_{0})}^{T}</math>.
 
  
It identifies the mutually dual LPs <math>\{{d}^{T}x : d \in {}^{\omega}\mathbb{R}^{n}, x \in {P}_{\ge 0}\}</math> and min <math>\{{b}^{T}y : y \in {}^{\omega}\mathbb{R}_{\ge 0}^{m}, {A}^{T}y \ge d\}</math>.
+
Its solution is the geometric centre <math>g</math> of the polytope <math>P_0</math>. For <math>p_k^* := (\text{min}\,p_k + \text{max}\,p_k)/2</math> and <math>k = 1, ..., \grave{z}</math> approximate <math>g</math> by <math>p_0 := (x_0, y_0, r_0)^T</math> until <math>||\Delta p||_1</math> is sufficiently small. The solution <math>t^o(x^o, y^o, r^o)^T</math> of the two-dimensional LP min <math>\{ r \in [0, \check{r}] : t \in {}^{\omega}\mathbb{R}_{&gt; 0}, t(x_0, y_0)^T \in P_r\}</math> approximates <math>g</math> better and achieves <math>r \le \check{r}/\sqrt{\grave{z}}</math>. Repeat this for <math>t^o(x^o, y^o)^T</math> until <math>g \in P_0</math> is computed in <math>\mathcal{O}({}_z\check{r} {}_e\check{r}dmn)</math> if it exists. Numbers of length <math>\mathcal{O}({\omega})</math> can only be processed in <math>\mathcal{O}(\vartheta)</math> as is generally known.
  
Hereon, we successively interpolate all <math>{v}_{k}^{*} := (\text{max } {v}_{k} + \text{min } {v}_{k})/2</math> until all <math>|\Delta{v}_{k}|</math> are sufficiently small. After that, we extrapolate then <math>v</math> via <math>{v}^{*}</math> into the boundary of the polytope. The <math>r</math>-fold of the distance exceeding <math>{v}^{*}</math> determines the new starting point <math>v</math>.
 
  
If min<math>{}_{k} {h}_{k} t = 0</math> follows from <math>t :=</math> min<math>{}_{k} \Delta{h}_{k}</math>, we end. Then we start over until min <math>h = 0</math> or min <math>h > 0</math> is certain. Since <math>h</math> at least halves itself for each iteration step in <math>\mathcal{O}({\omega\vartheta}^{2})</math>, the strong duality theorem yields the result.<math>\square</math>
+
Solving all two-dimensional LPs <math>\text{min}_k r_k</math> by bisection methods for <math>r_k \in {}^{\omega}\mathbb{R}_{\ge 0}</math> and <math>k = 1, ..., z</math> in <math>\mathcal{O}({\vartheta}^2)</math> each time determines <math>q \in {}^{\omega}\mathbb{R}^k</math> where <math>q_k := \Delta p_k \Delta r_k/r</math> and <math>r := \text{min}_k \Delta r_k</math>. Let simplified <math>|\Delta p_1| = … = |\Delta p_{z}|</math>. Here min <math>r_z</math> for <math>p^* := p + wq</math> and <math>w \in {}^{\omega}\mathbb{R}_{\ge 0}</math> would be also to solve. If <math>\text{min}_k \Delta r_k r = 0</math> follows, stop, otherwise repeat until min <math>r = 0</math> or min <math>r &gt; 0</math> is sure. If necessary, constraints are temporarily relaxed by the same small modulus.<math>\square</math>
  
 
== Recommended readings ==
 
== Recommended readings ==
[http://www.epubli.de/shop/buch/Relil-Boris-Haase-9783844208726/11049 Relil - Religion und Lebensweg]
 
 
 
[https://en.calameo.com/books/003777977258f7b4aa332 Nonstandard Mathematics]
 
[https://en.calameo.com/books/003777977258f7b4aa332 Nonstandard Mathematics]
  
 
[[de:Hauptseite]]
 
[[de:Hauptseite]]

Revision as of 01:55, 1 January 2022

Welcome to MWiki

Theorem of the month

The centre method solves every solvable LP in [math]\displaystyle{ \mathcal{O}(\omega{\vartheta}^{2}) }[/math].

Proof and algorithm

Let [math]\displaystyle{ z := \grave{m} + n }[/math] and [math]\displaystyle{ d \in [0, 1] }[/math] the density of [math]\displaystyle{ A }[/math]. First, normalise and scale [math]\displaystyle{ {b}^{T}y - {c}^{T}x \le 0, Ax \le b }[/math] as well as [math]\displaystyle{ {A}^{T}y \ge c }[/math]. Let [math]\displaystyle{ P_r := \{(x, y)^T \in {}^{\omega}\mathbb{R}_{\ge 0}^{z} : {b}^{T}y - {c}^{T}x \le r \in [0, \check{r}], Ax - b \le \underline{r}_m, c - {A}^{T}y \le \underline{r}_n\} }[/math] have the radius [math]\displaystyle{ \check{r} := s|\min \; \{b_1, ..., b_m, -c_1, ..., -c_n\}| }[/math] and the scaling factor [math]\displaystyle{ s \in [1, 2] }[/math]. It follows [math]\displaystyle{ \underline{0}_{z} \in \partial P_{\check{r}} }[/math]. By the strong duality theorem, the LP min [math]\displaystyle{ \{ r \in [0, \check{r}] : (x, y)^T \in P_r\} }[/math] solves the LPs max [math]\displaystyle{ \{{c}^{T}x : c \in {}^{\omega}\mathbb{R}^{n}, x \in {P}_{\ge 0}\} }[/math] and min [math]\displaystyle{ \{{b}^{T}y : y \in {}^{\omega}\mathbb{R}_{\ge 0}^{m}, {A}^{T}y \ge c\} }[/math].


Its solution is the geometric centre [math]\displaystyle{ g }[/math] of the polytope [math]\displaystyle{ P_0 }[/math]. For [math]\displaystyle{ p_k^* := (\text{min}\,p_k + \text{max}\,p_k)/2 }[/math] and [math]\displaystyle{ k = 1, ..., \grave{z} }[/math] approximate [math]\displaystyle{ g }[/math] by [math]\displaystyle{ p_0 := (x_0, y_0, r_0)^T }[/math] until [math]\displaystyle{ ||\Delta p||_1 }[/math] is sufficiently small. The solution [math]\displaystyle{ t^o(x^o, y^o, r^o)^T }[/math] of the two-dimensional LP min [math]\displaystyle{ \{ r \in [0, \check{r}] : t \in {}^{\omega}\mathbb{R}_{> 0}, t(x_0, y_0)^T \in P_r\} }[/math] approximates [math]\displaystyle{ g }[/math] better and achieves [math]\displaystyle{ r \le \check{r}/\sqrt{\grave{z}} }[/math]. Repeat this for [math]\displaystyle{ t^o(x^o, y^o)^T }[/math] until [math]\displaystyle{ g \in P_0 }[/math] is computed in [math]\displaystyle{ \mathcal{O}({}_z\check{r} {}_e\check{r}dmn) }[/math] if it exists. Numbers of length [math]\displaystyle{ \mathcal{O}({\omega}) }[/math] can only be processed in [math]\displaystyle{ \mathcal{O}(\vartheta) }[/math] as is generally known.


Solving all two-dimensional LPs [math]\displaystyle{ \text{min}_k r_k }[/math] by bisection methods for [math]\displaystyle{ r_k \in {}^{\omega}\mathbb{R}_{\ge 0} }[/math] and [math]\displaystyle{ k = 1, ..., z }[/math] in [math]\displaystyle{ \mathcal{O}({\vartheta}^2) }[/math] each time determines [math]\displaystyle{ q \in {}^{\omega}\mathbb{R}^k }[/math] where [math]\displaystyle{ q_k := \Delta p_k \Delta r_k/r }[/math] and [math]\displaystyle{ r := \text{min}_k \Delta r_k }[/math]. Let simplified [math]\displaystyle{ |\Delta p_1| = … = |\Delta p_{z}| }[/math]. Here min [math]\displaystyle{ r_z }[/math] for [math]\displaystyle{ p^* := p + wq }[/math] and [math]\displaystyle{ w \in {}^{\omega}\mathbb{R}_{\ge 0} }[/math] would be also to solve. If [math]\displaystyle{ \text{min}_k \Delta r_k r = 0 }[/math] follows, stop, otherwise repeat until min [math]\displaystyle{ r = 0 }[/math] or min [math]\displaystyle{ r > 0 }[/math] is sure. If necessary, constraints are temporarily relaxed by the same small modulus.[math]\displaystyle{ \square }[/math]

Recommended readings

Nonstandard Mathematics