Поиск
Тема: светлая

Делить и властвовать: повышение эффективности алгоритмов. Часть 2

Прокомментировать Просмотры: 26

Ссылка на первую часть

Мастер-теорема

На примере из прошлой части, попробуем сформулировать и обобщить принцип «Разделяй и властвуй». Мы беремся за проблему, размера n, делим эту проблему на подзадачи размером n/b. Количество таких подзадач обозначим числом a. И еще имеется задача скомпоновать результаты выполнения этих a задач размером n/b в итоговый результат для задачи размера n, который будем считать задачей полиномиальной сложности степени c, O(nc) . Если задача компоновки будет не полиномиальной, то все изложение резко усложнится. Поэтому, давайте позволим задаче компоновки быть полиномиальной, тем более в это попадает очень большое количество алгоритмов.

В примере выше, произведения двух чисел, мы разделили задачу длины n на три подзадачи длины n/2 и компоновали результат линейными операциями O(n). Значит, для этого случая: a=3, b=2, c=1

Выше мы уже выводили рекуррентную формулу для рекурсивного алгоритма метода «Разделяй и властвуй», которую просто обобщим:

T(n)=aT([n/b])+O(n^c)\quad (1)

Мастер-теорема. Пусть имеются рекуррентные соотношения:.

T(n)=\begin{equation*} \begin{cases} O(n^c),\quadесли\quad c>log_b(a); \\ O(n^clog_b(n)),\quadесли\quad c=log_b(a); \\ O(n^{log_b(a)}),\quadесли\quad c<log_b(a); \end{cases} \end{equation*}

Собственно, в прошлой части мы уже решали эту задачу, просто напомню. Строим рекурсивное дерево, где из узла-задачи длины n на уровни ниже расходятся a задач длины n/b. Высота дерева logb(n). Давайте для упрощения считать, что n это целая степень b, тогда и не надо мучиться с округлением и высота дерева целой будет. Это все таки не математическая публикация, а научно-популярная, поэтому давайте сделаем это упрощение.

Рекурсивное дерево "Разделяй и властвуй"
Рекурсивное дерево «Разделяй и властвуй»

На самом нижнем уровне рекурсивного дерева у нас nlog(a) задач размером O(1)

На уровнях выше выполняются задачи по компоновке результатов предыдущего уровня. Всего на уровне k будет ak таких задач, каждая со сложностью O((n/bk)c)

Таким образом, для оценки всего алгоритма нам нужно посчитать сумму сложности по всем уровням рекурсии. Мы уже проделывали подобные вычисления в прошлой части. Это будет выражение с суммой геометрической прогрессии:

T(n)=O(n^c\cdot\sum_{i=0}^{log_bn}(\frac{a}{b^c})^i)\quad (2)

Отсюда, после суммирования, получится три варианта для итогового выражения: когда геометрическая прогрессия убывает; когда геометрическая прогрессия растет; когда прогрессия отсутствует, то есть множитель прогрессии 1

В первом случае, имеем убывающую геометрическую прогрессию. Сумма прогрессии не превышает суммы бесконечной геометрической прогрессии с этим множителем, а та не больше константы, зависящей только от a,b,c. Константу можно вынести и останется только O(nc)

\sum_{i=0}^{log_bn}(\frac{a}{b^c})^i<\sum_{i=0}^{\infty}(\frac{a}{b^c})^i=\frac{1}{1-\frac{a}{b^c}}=const(a,b,c)

А значит выражение (2) вырождается в

T(n)=O(n^c)

Второй случай, когда геометрическая прогрессия растет, уже подробно рассматривался в первой части на примере алгоритма Карацубы. Имеем сумму некоторого конечного числа элементов геометрической прогрессии, которые можем вычислить по формуле для суммы геометрической прогрессии. В итоге, обобщив ответ из части 1, мы получим формулу

T(n)=O(n^{log_b(a)})

Наконец, последний случай, когда прогрессия отсутствует. В этом случае в формуле (2) будут единицы под суммой, а сама сумма вырождается в выражение количества элементов суммы. В итоге получаем выражение

T(n)=O(n^c\cdot log_bn)

Применения

В комментариях к прошлой части были озвучены ряд алгоритмов, иллюстрирующие концепцию «Разделяй и властвуй»
Самый простой случай, это Бинарный поиск. Известен еще с древних времен. Найти в отсортированном списке значение. Алгоритм общеизвестен. Но теперь рассмотрим его с точки зрения концепции «Разделяй и властвуй»
Делим список [0,…,n] на два равных списка с индексами: [0,…, n/2-1], [n/2,….,n] Мы сравниваем искомое значение с максимальным значением из первого списка и определяем, лежит ли искомое значение в первом списке или во втором. Далее рекурсивно обходим список, уменьшенный вполовину, а на второй список не обращаем внимание.

То есть, мы делим проблему длины n, на a=1 одну проблему длины n/2, b=2. И на каждом уровне выполняем одну операцию сравнения за константное время O(1), c=0. Рекурсивное соотношение «Разделяй и властвуй» (1) для бинарного поиска выглядит так:

T(n)=T(n/2)+O(1)

Т.е. a=1, b=2, c=0 и согласно мастер-теореме, это второй случай, где

c=log_b(a)=0

И по мастер-теореме мы получаем известную формулу сложности алгоритма бинарного поиска T(n)=O(log(n))

Сортировка слиянием алгоритм эффективной сортировки. Список значений с индексам [0,…,n], который нужно отсортировать делим на два равных списка с индексами: [0,…, n/2-1], [n/2,….,n]. Рекурсивно сортируем каждый из меньших списков. Далее собираем итоговый список сливая два отсортированных списка.

Для расчета сложности алгоритма, сначала надо понять процедуру, которой можно слить два отсортированных списка в один отсортированный и какова ее сложность.Процедура выглядит так. Берем первый наименьший элемент из списка 1 и сравниваем его с первым наименьшим элементом из списка 2. Кто из них самый наименьший — поступает в итоговый слитый список и убирается из списка 1 или списка 2. Далее рекурсивно повторяем алгоритм для двух списков, один из которых на этом шаге стал меньше на один элемент.

На каждый рекурсивный вызов процедуры слияния происходит одно сравнение, т.е. действие сложности O(1). И таких действий, из описания рекурсии, будет n, где n сумма размеров первоого и второго списка. Таким образом, процедура слияния будет иметь сложность O(n)

Теперь, определившись с алгоритмом, запишем рекуррентное соотношение (1) для мастер-теоремы

T(n)=2T(n/2)+O(n)

Получается a=2, b=2, c=1 и согласно мастер теореме, это второй случай, где

c=log_b(a)=1

И итоговое выражения для алгоритма сортировки слиянием

T(n)=O(n\cdot log(n))

Умножение двух матриц

Если у нас есть две квадратные матрицы A и B размерностью n, то их произведением Z=AxB называют матрицу элементы которой равны

Z_{i,j}=\sum_{i=1}^nA_{ik}B_{kj}

В данном случае, в отличие от прошлых задач, будем считать что элементы матрицы имеют длину, ограниченную байтом, словом и т.п. В этом случае умножение двух элементов считается константным временем.

Тогда умножение двух матриц по наивному алгоритму с использование формулы выше сводится к вычислению n2 элементов Z, каждое из которых является суммированием n произведений. Т.е. итоговая сложность умножения получается O(n3). Казалось бы улучшить тут нечего. Но на помощь приходит «Разделяй и властвуй»

Если размерность матрицы n четная, то такие матрицы можно умножать блоками размерности n/2 по формуле

X=\left[ \begin{array}{cc} A & B \\ C & D \end{array} \right];Y=\left[ \begin{array}{cc} E & F \\ G & H \end{array} \right];XY=\left[ \begin{array}{cc} A & B \\ C & D \end{array} \right]\left[ \begin{array}{cc} E & F \\ G & H \end{array} \right]=\left[ \begin{array}{cc} AE+BG & AF+BH \\ CE+DG & CF+DH \end{array} \right]\quad(3)

Если воспользоваться таким разбиением для «Разделяй и властвуй» то умножение двух матриц размерности n сводится к 8 произведениям матриц размерности n/2. После чего выполняется сложение этих матриц, что требует O(n2) операций для того, чтобы вычислить каждый из n2 элементов итоговой матрицы

Соотношение (1) в этом случае будет

T(n)=8T(n/2)+O(n^2)

Данное соотношение по мастер-теореме дает итоговую сложность O(n3), то есть никакого выигрыша в сравнении с наивным алгоритмом. Но вы уже догадываетесь, что как и при скалярном произведении, можно не делать все 8 произведений, а обойтись меньшим числом. И это действительно так, правда это не так тривиально, как в скалярном случае, но можно легко проверить, что соотношения верны. Достаточно выполнить 7 произведений матриц размерности n/2

P_1=AF-H; P_2=(A+B)H; P_3=(C+D)EP_4=D(G-E); P_5=(A+D)(E+H); P_6=(B-D)(G+H)P_7=(A-C)(E+F)

Тогда произведения для формулы (3) можно получить используя только сложения и вычитания матриц P1…P7

AE+BG=P_4+P_5+P_6-P_2;\quad CE+DG=P_3+P_4AF+BH=P_1+P2;\quad CF+DH=P_1+P_5-P_3-P_7

Теперь у нас 7 произведений вместо 8. Рекуррентное соотношение будет

T(n)=7T(n/2)+O(n^2)

И итоговая сложность алгоритма умножения матриц будет

T(n)=O(n^{log_27})\approx O(n^{2.81})

В третьей части посмотрим, как применяется «Разделяй и властвуй» к важным составляющим квантовых алгоритмов

 

Источник

Похожие статьи

Comments - Cackle

Поиск по играм, новостям и статьям…

Введите не менее двух символов

Введите не менее двух символов