Мастер-теорема — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
Строка 1: Строка 1:
'''Мастер теорема''' теорема позволяющая найти асимптотическое решение (с помощью [https://ru.wikipedia.org/wiki/%C2%ABO%C2%BB_%D0%B1%D0%BE%D0%BB%D1%8C%D1%88%D0%BE%D0%B5_%D0%B8_%C2%ABo%C2%BB_%D0%BC%D0%B0%D0%BB%D0%BE%D0%B5 О - большое нотации]) рекуррентных соотношений, которые могут возникнуть во многих алгоритмах, например таких как разделяй и властвуй. Однако не все рекуррентные соотношения могут быть решены через мастер теорему, ее обобщения включаются в метод Акра-Бацци.
+
'''Мастер теорема''' (Master theorem) теорема позволяющая найти асимптотическое решение рекуррентных соотношений, которые могут возникнуть в реализации многих алгоритмов. Однако не все рекуррентные соотношения могут быть решены через мастер теорему, ее обобщения включаются в метод [http://en.wikipedia.org/wiki/Akra%E2%80%93Bazzi_method Акра-Бацци].
  
 
==Формулировка и доказательство мастер-теоремы==
 
==Формулировка и доказательство мастер-теоремы==
 
{{
 
{{
 
Теорема|statement=
 
Теорема|statement=
Пусть у нас дано соотношение вида:  
+
Пусть, при реализации алгоритма мы получили соотношение такого вида:
  
<math> T(n) = \begin{cases}
+
<tex dpi = "135"> T(n) = \begin{cases}
 
   a \; T\!\left(\frac{n}{b}\right) + n^{c}  , &      n > 1\\  
 
   a \; T\!\left(\frac{n}{b}\right) + n^{c}  , &      n > 1\\  
 
   d    , &      n = 1
 
   d    , &      n = 1
 
\end{cases}
 
\end{cases}
</math>  
+
</tex>  
, где <math>a</math> — количество подзадач, на которые мы разбили нашу задачу, <math>n</math> — размер нашей задачи, <math>n / b</math> — размер подзадачи, <math> n ^ {c} </math> — стоимость работы, проделанной рекурсивными вызовами, который включает в себя стоимость деления проблемы и стоимость слияния решения подзадач, <math>d</math> — единичная стоимость для данной задачи.
+
 
Пусть <math>a</math> — <math>\mathbb N </math> число большее 1,  <math>b</math> — <math>\mathbb R </math> число большее 1, пусть также <math>c</math> — <math>\mathbb R^{+} </math> число и <math>d</math> — <math>\mathbb R^{+} </math> , тогда возможны три случая:
+
где <math>a</math> — количество подзадач, на которые мы разбиваем нашу задачу, <math>n</math> — размер нашей задачи, <tex dpi = "145">\frac{n}{b}</tex> — размер подзадачи, <math> n ^ {c} </math> — стоимость работы, проделанной рекурсивными вызовами, который включает в себя стоимость деления проблемы и стоимость слияния решения подзадач, <math>d</math> — единичная стоимость для данной задачи.
 +
Пусть <math>a</math> — <math>\mathbb N </math> число большее 1,  <math>b</math> — <math>\mathbb R </math> число большее 1, пусть также <math>c</math> — <math>\mathbb R^{+} </math> число и <math>d</math> — <math>\mathbb R^{+} </math> , тогда решение данного рекуррентного соотношения разбивается на три возможных случая:
  
 
1. Если <math>c > \log_b a</math>, то <math>T(n) = \Theta\left( n^{c} \right)</math>
 
1. Если <math>c > \log_b a</math>, то <math>T(n) = \Theta\left( n^{c} \right)</math>
Строка 20: Строка 21:
 
3. Если <math>c < \log_b a</math>, то <math>T(n) = \Theta\left( n^{\log_b a} \right)</math>
 
3. Если <math>c < \log_b a</math>, то <math>T(n) = \Theta\left( n^{\log_b a} \right)</math>
  
|proof= Для доказательства мы установим <math>d = 1</math>, это требуется для того, чтобы наши вычисления были хорошо определены при рекурсивном спуске.
+
|proof= Для доказательства мы установим <math>d = 1</math>, это требуется для того, чтобы при рекурсивном спуске не возникало огромных вычислений.  
 
Давайте рассмотрим дерево рекурсии. Всего в нем будет <math>\log_b n</math> уровней. На каждом таком уровне, количество подзадач будет умножаться на <math>a</math>, так на уровне <math>i</math> будет <math>a^i</math> подзадач. Также известно, что каждая подзадача на уровне <math>i</math> размера <math>n / b^i</math>. Подзадача размера <math>n / b^i</math> требует <math>(n / b^i) ^ c</math> дополнительных затрат, поэтому общее количество совершенных операций на уровне <math>i</math> :
 
Давайте рассмотрим дерево рекурсии. Всего в нем будет <math>\log_b n</math> уровней. На каждом таком уровне, количество подзадач будет умножаться на <math>a</math>, так на уровне <math>i</math> будет <math>a^i</math> подзадач. Также известно, что каждая подзадача на уровне <math>i</math> размера <math>n / b^i</math>. Подзадача размера <math>n / b^i</math> требует <math>(n / b^i) ^ c</math> дополнительных затрат, поэтому общее количество совершенных операций на уровне <math>i</math> :
<math>a^i(n / b^i)^c = n^c(a^i/b^(ic)) = n^c(a/b^c)^i</math>
+
<math>a^i(n / b^i)^c = n^c(a^i/b^{ic}) = n^c(a/b^c)^i</math>
Заметим, что количество занятой памяти увеличивается, уменьшается и остается константой, если <math>(a/b^c)^i</math> увеличивается, уменьшается или остается константой соответственно.
+
Заметим, что количество операций увеличивается, уменьшается и остается константой, если <math>(a/b^c)^i</math> увеличивается, уменьшается или остается константой соответственно.
Поэтому мы должны разобрать три случая, когда <math>(a/b^c)^i</math> больше 1, равен 1 или меньше 1.
+
Поэтому мы должны разобрать три случая, когда <math>(a/b^c)^i</math> больше <math>1</math>, равен <math>1</math> или меньше <math>1</math>.
 
Рассмотрим <math>(a/b^c)^i = 1</math> <math>\Leftrightarrow</math> <math>a = b^c</math> <math>\Leftrightarrow</math> <math>\log_b a = c \log_b b</math> <math>\Leftrightarrow</math>  <math>\log_b a = c</math>.
 
Рассмотрим <math>(a/b^c)^i = 1</math> <math>\Leftrightarrow</math> <math>a = b^c</math> <math>\Leftrightarrow</math> <math>\log_b a = c \log_b b</math> <math>\Leftrightarrow</math>  <math>\log_b a = c</math>.
 
Распишем всю работу в течение рекурсивного спуска:
 
Распишем всю работу в течение рекурсивного спуска:
<tex dpi = "130"> \displaystyle\sum_{i=1}^{log_b n}n^c(\frac{a}{b^c})^i = n^c\displaystyle\sum_{i=1}^{log_b n}(\frac{a}{b^c})^i</tex>
+
<tex dpi = "130"> \displaystyle\sum_{i=1}^{\log_b n}n^c(\frac{a}{b^c})^i = n^c\displaystyle\sum_{i=1}^{\log_b n}(\frac{a}{b^c})^i</tex>
 
Откуда получаем:
 
Откуда получаем:
  
 
1. <math>\log_b a < c </math> <math>\Rightarrow</math> <math>T(n) = \Theta\left( n^{c} \right)</math> (т.к. <tex dpi = "130"> (\frac{a}{b^c})^i</tex> убывающая геометрическая прогрессия)
 
1. <math>\log_b a < c </math> <math>\Rightarrow</math> <math>T(n) = \Theta\left( n^{c} \right)</math> (т.к. <tex dpi = "130"> (\frac{a}{b^c})^i</tex> убывающая геометрическая прогрессия)
  
2. <math>\log_b a = c </math> <math>\Rightarrow</math> <tex dpi = "130"> T(n) = \displaystyle\sum_{i=1}^{log_b n}n^c(\frac{a}{b^c})^i = n^c\displaystyle\sum_{i=1}^{log_b n}(\frac{a}{b^c})^i = n^c\displaystyle\sum_{i=1}^{log_b n}1^i = n^c + n^c\log_b n\ = \Theta\left( n^{c} \log n \right) </tex>
+
2. <math>\log_b a = c </math> <math>\Rightarrow</math> <tex dpi = "125"> T(n) = \displaystyle\sum_{i=1}^{\log_b n}n^c\cdot(\frac{a}{b^c})^i = </tex> <tex dpi = "125> n^c\cdot\displaystyle\sum_{i=1}^{\log_b n}(\frac{a}{b^c})^i =                                                                                                   n^c\cdot\displaystyle\sum_{i=1}^{\log_b n}1^i = n^c + n^c\log_b n = \Theta\left( n^{c} \log n \right) </tex>
  
3. <math>\log_b a > c </math> <math>\Rightarrow</math> <tex dpi = "130">T(n) = \displaystyle\sum_{i=1}^{log_b n}n^c(\frac{a}{b^c})^i = n^c\displaystyle\sum_{i=1}^{log_b n}(\frac{a}{b^c})^i = \Theta\left( n^c(\frac{a}{b^c})^{log_b n} \right)</tex>, но  <tex dpi = "150"> n^c(\frac{a}{b^c})^{log_b n} </tex> <tex dpi = "130"> =  </tex>  <tex dpi = "150">  n^c(\frac{a^{log_b n} }{(b^c)^{log_b n}})  </tex> <tex dpi = "130"> =  </tex> <tex dpi = "150">  n^c(\frac{n^{log_b a}}{n^c})</tex> <tex dpi = "130"> =  </tex>  <tex dpi = "130">  \Theta\left( n^{\log_b a} \right) </tex>
+
3. <math>\log_b a > c </math> <math>\Rightarrow</math> <tex dpi = "125"> T(n) = \displaystyle\sum_{i=1}^{\log_b n}n^c\cdot(\frac{a}{b^c})^i = n^c\cdot\displaystyle\sum_{i=1}^{\log_b n}(\frac{a}{b^c})^i = \Theta\left( n^c\cdot(\frac{a}{b^c})^{log_b n} \right)</tex>, но  <tex dpi = "150"> n^c\cdot(\frac{a}{b^c})^{log_b n} </tex> <tex dpi = "130"> =  </tex>  <tex dpi = "150">  n^c\cdot(\frac{a^{log_b n} }{(b^c)^{log_b n}})  </tex> <tex dpi = "130"> =  </tex> <tex dpi = "150">  n^c\cdot(\frac{n^{log_b a}}{n^c})</tex> <tex dpi = "130"> =  </tex>  <tex dpi = "130">  \Theta\left( n^{\log_b a} \right) </tex>  
  
 
}}
 
}}
  
 
==Примеры==   
 
==Примеры==   
1.Пусть у нас задано такое рекуррентное соотношение:
+
 
 +
=== Примеры задач ===
 +
 
 +
1.Пусть задано такое рекуррентное соотношение:
  
 
Рассчитать для <math>x = 7</math>.
 
Рассчитать для <math>x = 7</math>.
Строка 62: Строка 66:
 
</math>
 
</math>
  
<math>f(n) = n\sqrt n > n^{3/2} = O(n^{3/2}) </math>, а также  
+
<math>f(n) = n\sqrt {n + 1} > n^{3/2} = O(n^{3/2}) </math>, а также  
<math>f(n) = n\sqrt n < n\sqrt{n + n} < n\sqrt{2n} = O(n^{3/2}) </math>
+
<math>f(n) = n\sqrt {n + 1} < n\sqrt{n + n} < n\sqrt{2n} = O(n^{3/2}) </math>
  
==Недопустимые соотношения==
+
=== Недопустимые соотношения ===
 
Рассмотрим пару ошибочно-составленных соотношений:
 
Рассмотрим пару ошибочно-составленных соотношений:
 
*<math>T(n) = 2^nT\left (\frac{n}{2}\right )+n^n</math>
 
*<math>T(n) = 2^nT\left (\frac{n}{2}\right )+n^n</math>
 
*:<math>a</math> не является константой; количество подзадач может меняться
 
*:<math>a</math> не является константой; количество подзадач может меняться
 
*<math>T(n) = 2T\left (\frac{n}{2}\right )+\frac{n}{\log n}</math>
 
*<math>T(n) = 2T\left (\frac{n}{2}\right )+\frac{n}{\log n}</math>
*:не полиномиальное различие <math>f(n)</math> и <math>n^{\log_b a}</math>
+
*:не полиномиальное различие <math>f(n)</math> и <math>n^{\log_b a}</math>, т.к. <tex dpi = "145"> \frac{f(n)}{n^{\log_b a}} < n^r , для любого r > 0</tex>  
 
*<math>T(n) = 0.5T\left (\frac{n}{2}\right )+n</math>
 
*<math>T(n) = 0.5T\left (\frac{n}{2}\right )+n</math>
 
*:<math>a</math> < 1 не может быть меньше одной подзадачи
 
*:<math>a</math> < 1 не может быть меньше одной подзадачи
Строка 92: Строка 96:
 
| По мастер-теореме <math>c = \log_b a</math>, где <math>a = 1, b = 2, c = 0</math>
 
| По мастер-теореме <math>c = \log_b a</math>, где <math>a = 1, b = 2, c = 0</math>
 
|-
 
|-
| Обход бинарного дерева
+
| Обход [http://neerc.ifmo.ru/wiki/index.php?title=%D0%94%D0%B5%D1%80%D0%B5%D0%B2%D0%BE_%D0%BF%D0%BE%D0%B8%D1%81%D0%BA%D0%B0,_%D0%BD%D0%B0%D0%B8%D0%B2%D0%BD%D0%B0%D1%8F_%D1%80%D0%B5%D0%B0%D0%BB%D0%B8%D0%B7%D0%B0%D1%86%D0%B8%D1%8F двоичного дерева]
 
| <math>T(n) = 2 T\left(\frac{n}{2}\right) + O(1)</math>
 
| <math>T(n) = 2 T\left(\frac{n}{2}\right) + O(1)</math>
 
| <math>O(n)</math>
 
| <math>O(n)</math>
| По мастер-теореме <math>c < \log_b a</math> где <math>a = 2, b = 2, c = 0</math>
+
| По мастер-теореме <math>c < \log_b a</math>, где <math>a = 2, b = 2, c = 0</math>
 
|-
 
|-
 
|  [[Сортировка слиянием]]
 
|  [[Сортировка слиянием]]
Строка 103: Строка 107:
 
|}
 
|}
  
== Cсылки ==
+
 
 +
== Источники информации ==
 
* [http://en.wikipedia.org/wiki/Master_theorem Википедия — Мастер-теорема]
 
* [http://en.wikipedia.org/wiki/Master_theorem Википедия — Мастер-теорема]
 
* [https://math.dartmouth.edu/archive/m19w03/public_html/Section5-2.pdf Dartmouth university — the master theorem]
 
* [https://math.dartmouth.edu/archive/m19w03/public_html/Section5-2.pdf Dartmouth university — the master theorem]
 
== Литература ==
 
 
*''Кормен Т., Лейзерсон Ч., Ривест Р., Штайн К.'' Алгоритмы: построение и анализ, 2-е издание. М.: Издательский дом "Вильямс", 2005. ISBN 5-8459-0857-4
 
*''Кормен Т., Лейзерсон Ч., Ривест Р., Штайн К.'' Алгоритмы: построение и анализ, 2-е издание. М.: Издательский дом "Вильямс", 2005. ISBN 5-8459-0857-4
 
==См. также==
 
 
* [[Амортизационный анализ]]
 
* [[Амортизационный анализ]]
  
 
[[Категория:Дискретная математика и алгоритмы]]
 
[[Категория:Дискретная математика и алгоритмы]]
 
[[Категория:Амортизационный анализ]]
 
[[Категория:Амортизационный анализ]]

Версия 23:52, 6 мая 2015

Мастер теорема (Master theorem) теорема позволяющая найти асимптотическое решение рекуррентных соотношений, которые могут возникнуть в реализации многих алгоритмов. Однако не все рекуррентные соотношения могут быть решены через мастер теорему, ее обобщения включаются в метод Акра-Бацци.

Формулировка и доказательство мастер-теоремы

Теорема:
Пусть, при реализации алгоритма мы получили соотношение такого вида:

[math] T(n) = \begin{cases} a \; T\!\left(\frac{n}{b}\right) + n^{c} , & n \gt 1\\ d , & n = 1 \end{cases} [/math]

где [math]a[/math] — количество подзадач, на которые мы разбиваем нашу задачу, [math]n[/math] — размер нашей задачи, [math]\frac{n}{b}[/math] — размер подзадачи, [math] n ^ {c} [/math] — стоимость работы, проделанной рекурсивными вызовами, который включает в себя стоимость деления проблемы и стоимость слияния решения подзадач, [math]d[/math] — единичная стоимость для данной задачи. Пусть [math]a[/math][math]\mathbb N [/math] число большее 1, [math]b[/math][math]\mathbb R [/math] число большее 1, пусть также [math]c[/math][math]\mathbb R^{+} [/math] число и [math]d[/math][math]\mathbb R^{+} [/math] , тогда решение данного рекуррентного соотношения разбивается на три возможных случая:

1. Если [math]c \gt \log_b a[/math], то [math]T(n) = \Theta\left( n^{c} \right)[/math]

2. Если [math]c = \log_b a[/math], то [math]T(n) = \Theta\left( n^{c} \log n \right)[/math]

3. Если [math]c \lt \log_b a[/math], то [math]T(n) = \Theta\left( n^{\log_b a} \right)[/math]
Доказательство:
[math]\triangleright[/math]

Для доказательства мы установим [math]d = 1[/math], это требуется для того, чтобы при рекурсивном спуске не возникало огромных вычислений. Давайте рассмотрим дерево рекурсии. Всего в нем будет [math]\log_b n[/math] уровней. На каждом таком уровне, количество подзадач будет умножаться на [math]a[/math], так на уровне [math]i[/math] будет [math]a^i[/math] подзадач. Также известно, что каждая подзадача на уровне [math]i[/math] размера [math]n / b^i[/math]. Подзадача размера [math]n / b^i[/math] требует [math](n / b^i) ^ c[/math] дополнительных затрат, поэтому общее количество совершенных операций на уровне [math]i[/math] : [math]a^i(n / b^i)^c = n^c(a^i/b^{ic}) = n^c(a/b^c)^i[/math] Заметим, что количество операций увеличивается, уменьшается и остается константой, если [math](a/b^c)^i[/math] увеличивается, уменьшается или остается константой соответственно. Поэтому мы должны разобрать три случая, когда [math](a/b^c)^i[/math] больше [math]1[/math], равен [math]1[/math] или меньше [math]1[/math]. Рассмотрим [math](a/b^c)^i = 1[/math] [math]\Leftrightarrow[/math] [math]a = b^c[/math] [math]\Leftrightarrow[/math] [math]\log_b a = c \log_b b[/math] [math]\Leftrightarrow[/math] [math]\log_b a = c[/math]. Распишем всю работу в течение рекурсивного спуска: [math] \displaystyle\sum_{i=1}^{\log_b n}n^c(\frac{a}{b^c})^i = n^c\displaystyle\sum_{i=1}^{\log_b n}(\frac{a}{b^c})^i[/math] Откуда получаем:

1. [math]\log_b a \lt c [/math] [math]\Rightarrow[/math] [math]T(n) = \Theta\left( n^{c} \right)[/math] (т.к. [math] (\frac{a}{b^c})^i[/math] убывающая геометрическая прогрессия)

2. [math]\log_b a = c [/math] [math]\Rightarrow[/math] [math] T(n) = \displaystyle\sum_{i=1}^{\log_b n}n^c\cdot(\frac{a}{b^c})^i = [/math] [math] n^c\cdot\displaystyle\sum_{i=1}^{\log_b n}(\frac{a}{b^c})^i = n^c\cdot\displaystyle\sum_{i=1}^{\log_b n}1^i = n^c + n^c\log_b n = \Theta\left( n^{c} \log n \right) [/math]

3. [math]\log_b a \gt c [/math] [math]\Rightarrow[/math] [math] T(n) = \displaystyle\sum_{i=1}^{\log_b n}n^c\cdot(\frac{a}{b^c})^i = n^c\cdot\displaystyle\sum_{i=1}^{\log_b n}(\frac{a}{b^c})^i = \Theta\left( n^c\cdot(\frac{a}{b^c})^{log_b n} \right)[/math], но [math] n^c\cdot(\frac{a}{b^c})^{log_b n} [/math] [math] = [/math] [math] n^c\cdot(\frac{a^{log_b n} }{(b^c)^{log_b n}}) [/math] [math] = [/math] [math] n^c\cdot(\frac{n^{log_b a}}{n^c})[/math] [math] = [/math] [math] \Theta\left( n^{\log_b a} \right) [/math]
[math]\triangleleft[/math]

Примеры

Примеры задач

1.Пусть задано такое рекуррентное соотношение:

Рассчитать для [math]x = 7[/math].

[math] t(x) = \begin{cases} 3 \; t\!\left(\frac{x}{2}\right) + x^{2} , & x \gt 2\\ 5x , & 1 \lt x \lt 2 \end{cases} [/math]

Заметим, чтобы узнать [math]t(7)[/math] , мы должны знать [math]t(7/2)[/math], чтобы узнать [math]t(7/2)[/math], мы должны узнать [math]t(7/4)[/math], [math]1 \lt 7/4 \lt 2[/math], тогда [math]t(7/4) = 35/4[/math] , [math]t(7/2) = 3*35/4 + 49/4[/math], тогда [math]t(7) = 3t(7/2) + 7^2 = 329/2[/math]


2. Задано такое соотношение:

[math]f(n) =[/math] [math]n\sqrt{n + 1}[/math]

[math] T(n) = \begin{cases} 2 \; T\!\left(\frac{n}{3}\right) + f(n) , & n \gt 1\\ d , & n = 1 \end{cases} [/math]

[math]f(n) = n\sqrt {n + 1} \gt n^{3/2} = O(n^{3/2}) [/math], а также [math]f(n) = n\sqrt {n + 1} \lt n\sqrt{n + n} \lt n\sqrt{2n} = O(n^{3/2}) [/math]

Недопустимые соотношения

Рассмотрим пару ошибочно-составленных соотношений:

  • [math]T(n) = 2^nT\left (\frac{n}{2}\right )+n^n[/math]
    [math]a[/math] не является константой; количество подзадач может меняться
  • [math]T(n) = 2T\left (\frac{n}{2}\right )+\frac{n}{\log n}[/math]
    не полиномиальное различие [math]f(n)[/math] и [math]n^{\log_b a}[/math], т.к. [math] \frac{f(n)}{n^{\log_b a}} \lt n^r , для любого r \gt 0[/math]
  • [math]T(n) = 0.5T\left (\frac{n}{2}\right )+n[/math]
    [math]a[/math] < 1 не может быть меньше одной подзадачи
  • [math]T(n) = 64T\left (\frac{n}{8}\right )-n^2\log n[/math]
    [math]f(n)[/math] не положительна
  • [math]T(n) = T\left (\frac{n}{2}\right )+n(2-\cos n)[/math]
    регулярно меняющееся [math]f(n)[/math]


Приложение к известным алгоритмам

Алгоритм Рекуррентное соотношение Время работы Комментарий
Целочисленный двоичный поиск [math]T(n) = T\left(\frac{n}{2}\right) + O(1)[/math] [math]O(\log n)[/math] По мастер-теореме [math]c = \log_b a[/math], где [math]a = 1, b = 2, c = 0[/math]
Обход двоичного дерева [math]T(n) = 2 T\left(\frac{n}{2}\right) + O(1)[/math] [math]O(n)[/math] По мастер-теореме [math]c \lt \log_b a[/math], где [math]a = 2, b = 2, c = 0[/math]
Сортировка слиянием [math]T(n) = 2 T\left(\frac{n}{2}\right) + O(n)[/math] [math]O(n \log n)[/math] По мастер-теореме [math]c = \log_b a[/math], где [math]a = 2, b = 2, c = 1[/math]


Источники информации