Амортизационный анализ

Основные определения

Определение:

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

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

Определение:

Средняя амортизационная стоимость операций — величина [math]a[/math], находящаяся по формуле: [math]a = \frac{\sum\limits^{n}_{i = 1} {t_i}}{n}[/math], где [math]t_1,t_2, ... t_n[/math] - время выполнения операций [math]1,2, ... , n,[/math] совершённых над структурой данных.

Амортизационный анализ использует следующие методы:

1. Метод усреднения (метод группового анализа).

2. Метод потенциалов.

3. Метод предоплаты (метод бухгалтерского учета).

Метод усреднения

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

Пример

Рассмотрим стек с операцией [math]multipop(a)[/math] — извлечение из стека [math]a[/math] элементов. В худшем случае она работает за [math]O(n)[/math] времени, если удаляются все элементы массива. Однако прежде чем удалить элемент, его нужно добавить в стек. Итак, если в стеке было не более [math]n[/math] элементов, то в худшем случае с каждым из них могли быть произведены 2 операции - добавление в стек и извлечение из него. Например, если было [math]n[/math] операций [math]push[/math] - добавление в стек, стоимость каждой [math]O(1)[/math], и одна операция [math]multipop(n)[/math], то суммарное время всех операций — [math]O(n)[/math], всего операций [math]n + 1[/math], а значит, амортизационная стоимость операции — [math]O(1)[/math].

Математическое обоснование:

Пусть [math]n[/math] — количество операций, [math]m[/math] — количество элементов, задействованных в этих операциях. Очевидно [math]m \leq n[/math] Тогда:

[math]a = \frac{\sum\limits^{n}_{i=1} {t_i}}{n} = \frac{\sum\limits^{n}_{i=1} \sum\limits^{m}_{j=1} {t_{ij}}}{n} = \frac{\sum\limits^{m}_{j=1} \sum\limits^{n}_{i=1} {t_{ij}}}{n},[/math] где [math]{t_{ij}}[/math] — стоимость [math]i[/math]-ой операции над [math]j[/math]-ым элементом. Величина [math]\sum\limits^{n}_{i=1} {t_{ij}}[/math] не превосходит 2, т. к. над элементом можно совершить только 2 операции, стоимость которых равна 1 — добавление и удаление. Тогда:

[math]a \leq \frac{2m}{n} \leq 2,[/math] так как [math]m \leq n[/math].

Таким образом, cредняя амортизационная стоимость операций [math]a = O(1)[/math].

Метод потенциалов

Теорема (О методе потенциалов):

Введём для каждого состояния структуры данных величину [math]\Phi[/math] — потенциал. Изначально потенциал равен [math]\Phi_0[/math], а после выполнения [math]i[/math]-ой операции — [math]\Phi_i[/math]. Стоимость [math]i[/math]-ой операции обозначим [math]a_i = t_i + \Phi_i - \Phi_{i-1}[/math]. Пусть [math]n[/math] — количество операций, [math]m[/math] — размер структуры данных. Тогда средняя амортизационная стоимость операций [math]a = O(f(n, m)),[/math] если выполнены два условия: 1) Для любого [math]i: \enskip a_i = O(f(n, m))[/math] 2) Для любого [math]i: \enskip \Phi_i = O(n \relax f(n, m))[/math]

Доказательство:

[math]\triangleright[/math]

[math]a = \frac{\sum\limits^{n}_{i = 1} {t_i}}{n} = \frac{\sum\limits^{n}_{i = 1} {a_i} + \sum\limits^{n - 1}_{i = 0} {\Phi_i} - \sum\limits^{n}_{i = 1} {\Phi_i} }{n} = \frac{n \relax O(f(n, m)) + \Phi_0 - \Phi_n}{n} = O(f(n, m))[/math]

[math]\triangleleft[/math]

Пример

В качестве примера вновь рассмотрим стек с операцией [math]multipop(a)[/math]. Пусть потенциал — это количество элементов в стеке. Тогда:

1.1) [math]a_{push} = 1 + 1 = 2,[/math] т. к. время выполнения операции [math]push[/math] — 1, и изменение потенциала — тоже 1.

1.2) [math]a_{pop} = 1 - 1 = 0,[/math] т. к. время выполнения операции [math]pop[/math] — 1, а изменение потенциала — -1.

1.3) [math]a_{multipop} = k - k = 0,[/math] т. к. время выполнения операции [math]multipop(k)[/math] — k, а изменение потенциала — -k.

2) Для любого [math]i: \enskip \Phi_i = O(n),[/math] так как элементов в стеке не может быть больше [math]n[/math]

Таким образом, [math]f(n, m) = 1[/math], а значит, cредняя амортизационная стоимость операций [math]a = O(1)[/math].

Метод предоплаты

Представим, что использование определенного количества времени равносильно использованию определенного количества монет (плата за выполнение каждой операции). В методе предоплаты каждому типу операций присваивается своя учётная стоимость. Эта стоимость может быть больше фактической, в таком случае лишние монеты используются как резерв для выполнения других операций в будущем, а может быть меньше, тогда гарантируется, что текущего накопленного резерва достаточно для выполнения операции. Для доказательства оценки средней амортизационной стоимости [math]O(f(n, m))[/math] нужно построить учётные стоимости так, что для каждой операции она будет составлять [math]O(f(n, m))[/math]. Тогда для последовательности из [math]n[/math] операций суммарно будет затрачено [math]n \cdot O(f(n, m))[/math] монет, следовательно, cредняя амортизационная стоимость операций будет [math]a = \frac{\sum\limits^{n}_{i = 1} {t_i}}{n} = \frac{n \cdot O(f(n, m))}{n}[/math] [math]= O(f(n, m))[/math].

Пример

Опять же рассмотрим стек с операцией [math]multipop(a)[/math]. При выполнении операции [math]push[/math] будем использовать две монеты — одну для самой операции, а вторую — в качестве резерва. Тогда для операций [math]pop[/math] и [math]multipop[/math] учётную стоимость можно принять равной нулю и использовать для удаления элемента монету, оставшуюся после операции [math]push[/math].

Таким образом, для каждой операции требуется [math]O(1)[/math] монет, а значит, cредняя амортизационная стоимость операций [math]a = O(1)[/math].

Литература

Томас Кормен. Алгоритмы. Построение и анализ. - Санкт-Петербург, 2005. стр. 483-491.