Шаг 17.
Сложность алгоритма.
Возведение целого числа без знака в положительную целую степень

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

    Вычисление у = хn для вещественных x и n не представляет никакого труда: y = en*ln(x). Но для целых чисел этот способ неприменим, поскольку вычисления экспоненты и логарифма на цифровой вычислительной машине принципиально являются приближенными и не будут давать целочисленный результат. Не поможет здесь и округление, так как крайне трудно достоверно вычислить величину ошибки при приближенных расчетах, следовательно, не известно ни в какую сторону округлять, ни является ли истинным результатом ближайшее целое или ошибка составляет несколько единиц.

    Особенно важна эта проблема для длинных целых x и больших n, что встречается в некоторых критических приложениях, например, в шифровании информации.

    Правильный результат будет получен, если реализуем возведение в степень через целочисленное умножение. Решение "в лоб" - перемножение n экземпляров x, на что требуется (n-1) операция умножения.

    Более остроумный метод состоит в том, чтобы на каждом шаге вычислений активно использовать достижения предыдущих шагов. Пусть, например, требуется вычислить x1024. Вместо 1023 умножений можно построить следующую цепочку: х = х1, х2 = (x1)2, х4 = (x2)2, x8 = (x4)2, х16, х32, х64, х128, х256, х512, х1024. На каждом шаге получается очередная степень x, которая на следующем шаге умножается сама на себя, т. е. предыдущее значение возводится в квадрат. Результат, таким образом, достигается за десять шагов.

    Этот метод хорош, если n равно степени двойки. Его сложность равна log n умножений. Для других показателей он требует видоизменения.

    Пусть, например, требуется вычислить х1023. Очевидным подходом является достижение за девять умножений значения х512 и перемножение еще за девять операций всех ранее найденных степеней: 1023 = 1 + 2 + 4 + 8 + ... + 256 + 512. Интересно, что для меньшей степени (1023 < 1024) требуется на восемь умножений больше!

    Для произвольного n можно использовать так называемый бинарный метод. Представим n его двоичной записью. Читая запись слева направо, начиная с первой значащей цифры (единицы), будем производить следующие действия: первую единицу игнорируем; переменной y присваиваем начальное значение, равное x; далее в цикле пока не закончились цифры - умножаем y на себя, а если очередная цифра равна единице, то дополнительно умножаем y на x. По окончании цикла в y находится значение xn.

    Эти действия можно описать следующей процедурой.

procedure BinaryMethod (x, n: integer; var у: integer); 
var 
  В: integer; 
begin
  у:= x;
  В:= NextBit(n);
  while В = 0 do 
        В:= NextBit(n); {В B - первая единица.} 
  В:= NextBit(n);
  while not Eon(n) do {Пока не закончилась запись числа n.} 
    begin
      у:= у*у; 
      if В = 1 then у:= у*х;
    end;
end;

    Сложность процедуры зависит от количества единиц λ(n) в записи числа n и может быть определена формулой:

Т(n) = log n + λ(n)-1.

    Минимальная сложность получается для степеней двойки, n = 2m-1, λ(n)=1. Максимальная сложность - для чисел n = 2m-1:

λ(n) = log2(n+l).

    Что общего между рассмотренными методами оптимизации алгоритмов?

    Стандартные алгоритмы решения задач перемножения матриц, длинных чисел, возведения в степень сводили задачу к длинной последовательности однородных подзадач малого размера. Эта последовательность записывалась в виде итеративного алгоритма, складывающего простым (стандартным) образом общее решение из решений независимых малых подзадач. Малые подзадачи отличаются по своей постановке от исходной задачи.

    Оптимизированные алгоритмы делят задачу на малое число подзадач относительно большого размера (но, разумеется, меньшего, чем исходная задача). Подзадачи являются уменьшенными копиями исходной задачи, поэтому алгоритм удобно сформулировать в рекурсивной форме. Однако не всегда этот подход приводит к желаемому ускорению. Часто требуется некий нестандартный, усложненный метод соединения подзадач, позволяющий обойтись меньшим их числом. При этом подзадачи оказываются зависимыми.

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




Предыдущий шаг Содержание