В работе рассматривается актуальная задача повышения эффективности реализации алгебраических операций в конечных алгебраических структурах. Она имеет существенное значение в информационных технологиях, в частности в обеспечении быстродействия систем цифровой обработки сигналов и криптографических протоколов. Этой задаче посвящено значительное количество как теоретических работ, так и работ практической направленности.
Можно считать, что изучение асимптотически быстрых алгоритмов умножения началось еще до появления понятия криптографии с открытым ключом, в, 1962 году, когда А. Карацуба предложил использовать метод умножения, имеющии асимптотическую сложность [33]. Несколько позднее, Тоомом и Куком в 1963;1966 г. в работах [41] и [49] было показано, что данную асимптотическую оценку можно улучшить, до oip.2^2log2″ log2 nj. А в 1966;1971 году в работах [76, 78−80] был предложенметод умножения, использующий быстрое дискретное преобразование Фурье и показано, что он имеет асимптотическую сложность 0(nog2 nog2 log2(n)).
Изучение операции инвертирования в поле началось намного позднее и нашло свое отражение в работах [11, 12, 20, 35, 65, 84]. В работах [11, 12, 65] рассматриваются алгоритмы инвертирования, использующее следствие из малой теоремы Ферма, в [35, 72, 79, 84] используются различные вариации расширенного алгоритма Евклида, имеющего асимптотическую сложность о (п3), а в [44] был предложен более эффективный вариант расширенного алгоритма, однако известен теоретически более быстрый вариант Шенхаге-Моенка алгоритма Евклида, его асимптотическая сложность оценивается как O (M (rc)log0z)), где М (и) — сложность умножения многочленов степени п [72, 79]. В работе [20] в 2006 году алгоритмы инвертирования в бинарных полях были обобщены на поля характеристики большей или равной трех.
Актуальность темы
Понижение асимптотической сложности не всегда сопровождается понижением практической сложности, под которой мы понимаем число операций или время, затрачиваемые для выполнения операции при заданных размерностях операндов. Несмотря на сравнительно малую асимптотическую сложность изученных методов, при программной реализации в оценке сложности получается довольно большая мультипликативная константа, что существенно ограничивает целесообразность их применения на практике и делает актуальным исследование других, в том числе обладающих большей асимптотической сложностью методов, которые дают лучшие оценки производительности по времени для наиболее часто используемых в современных информационных системах размерностей операндов.
Не случайно на рубеже тысячелетий появились новые интерпретации давно известного метода сдвигов и сложений для реализации умножения [70], а также новые методы деления и инвертирования в конечных полях [81], имеющие лучшую производительность, чем асимптотически оптимальные методы в диапазоне размерностей современных криптографических протоколов.
Следует также отметить, что асимптотические оценки сложности схемной и программной реализации методов одинаковы и не учитывают той специфики последней, что производительность программы зависит от взаимного расположения операндов, а также от вспомогательных действий, связанных с организацией вычислений по программе (условных переходов, циклов). Чтобы сократить объем таких действий была предложена программная реализация метода Карацубы в виде декомпозиционной схемы, имеющей структуру ХПХХ.
П].
Для получения объективных оснований выбора того или иного метода при реализации конкретной информационной системы необходимы разработка и сравнительный анализ новых и известных методов (включая их комбинирование) программной реализации алгебраических операций в тех или иных диапазонах размерностей операндов, чему и посвящена тема настоящей диссертации. Актуальность этой тематики проявляется и в бурном развитии эллиптической криптографии, где появились новые криптографические протоколы, не имеющие аналогов в классе протоколов, основанных на свойствах мультипликативной группы [7].
Целью диссертации является повышение эффективности программной реализации арифметических операций в конечных алгебраических структурах.
При этом ставились следующие задачи:
1) систематизация и программная реализация типичных алгоритмов в конечных алгебраических структурах;
2) разработка способов сравнения и оценки методов программной реализации алгебраических операций;
3) разработка методов автоматической генерации последовательных программ умножения;
4) сравнительный анализ и обоснование выбора методов программной реализации умножения в конечных полях;
5) сравнительный анализ и обоснование выбора методов программной реализации производных операций (возведения в степень, инвертирования и деления) в конечных полях;
6) сравнительный анализ методов программной реализации операций сложения и скалярного умножения в группе точек эллиптической кривой на основе рекомендуемых в диссертации методов умножения и деления в конечных полях.
Основные полученные в диссертации результаты:
1. Обоснована целесообразность комбинирования метода умножения путем сдвигов и сложений (вспомогательного метода) и метода Карацубы с уточнением порога перехода к вспомогательному методу умножения для комбинированных методов, построенных с использованием рекурсивной, циклической и последовательной реализации метода Карацубы.
2. Разработан метод автоматического синтеза последовательных программ умножения многочленов, обеспечивающий локализацию адресов, и способ автоматической верификации результата такого синтеза.
3. Разработан метод аналитического сравнения производительности последовательных программ умножения.
4. Определены наилучшие условия применения последовательных программ умножения по методу А. Карацубы в процессе умножения многочленов высоких степеней над полями различных характеристик и показано существенно большее быстродействие последовательных программ по сравнению с рекурсивными.
5. Выявлены области рекомендуемого применения методов при умножении многочленов высоких степеней над полями различных характеристик.
6. На основе сравнительного анализа эффективности алгоритмов скалярного умножения на эллиптических кривых с использованием аффинной и проективной систем представления операндов подтверждено, что существенное преимущество проективной системы сохраняется и в условиях применения наиболее эффективных алгоритмов умножения, инвертирования и деления в конечных полях.
Научная новизна полученных результатов:
1. Уточнен порог перехода для комбинированных методов, построенных с использованием рекурсивной, циклической и последовательной реализации метода Карацубы.
2. Предложен метод автоматического подтверждения правильности синтезируемых различными способами программ умножения.
3. Предложены аналитический метод оценки сложности и метод автоматического синтеза последовательных программ умножения по методу Карацубы, обеспечивающий локализацию адресов.
4. Подтверждена эффективность применения проективных координат при вычислениях в группе точек эллиптической кривой в условиях применения наиболее эффективных методов умножения и деления в конечных полях.
Методы исследования. Поставленные задачи решались с использованием методов дискретной математики, линейной алгебры, теории вычислительной сложности алгоритмов. При разработке программного обеспечения использовались методы объектно-ориентированного программирования и техники оптимизации программного обеспечения.
Практическая значимость полученных результатов. Применение предложенного метода автоматической проверки правильности синтезированных программ позволит повысить качество разрабатываемого программного обеспечения и снизить временные затраты на его разработку и тестирование за счет формальной проверки 100% исходного кода еще до компиляции.
Полученные аналитическим и экспериментальным путями оценки и интервалы оптимальности методов умножения представляют практический интерес при программировании новых производных алгоритмов, использующих базовые операции.
Применение в алгоритмах на эллиптических кривых предложенных модификаций алгоритмов умножения в расширениях полей степени 256 и выше позволило сократить время выполнения операции умножения точки эллиптической кривой на константу на 5−15% в зависимости от типа кривой. Реализация результатов. Работа выполнялась в соответствии с проектом РФФИ 08−01−632-а и заданием по разработке программного средства «Алгебраический процессор» инновационной образовательной программы. Результаты внедрены в виде математического и методологического обеспечения лабораторной работы «Сравнительный анализ временных характеристик алгоритмов операций в конечных группах, кольцах и полях» электронного образовательного ресурса «Алгебраический процессор» [43].
Апробация полученных результатов осуществлена при их обсуждении на одиннадцатой, двенадцатой,. четырнадцатой и пятнадцатой международных научно-технических конференциях студентов и аспирантов «Радиоэлектроника, электротехника и энергетика», Москва 2005, 2006, 2008, 2009; девятой, десятой и одиннадцатой Московских Международных телекоммуникационных конференциях студентов и молодых ученых «Молодежь и наука», Москва 2006, 2007, 2009.
Полученные результаты опубликованы в десяти печатных изданиях [23−31, 43], включая 3 статьи из перечня ВАК [23−25].
Основное содержание работы. Диссертация состоит из настоящего введения, четырех глав, заключения и библиографического списка.
Во введении обоснована актуальность темы диссертации, ее научная новизна и практическая значимость, сформулирована цель ' работы, рассматриваемые в ней задачи, и приведено краткое содержание диссертации по главам.
В первой главе рассмотрены известные алгоритмы выполнения базовых арифметических операций в конечных алгебраических структурах. Приводятся базовые определения и понятия, описано современное состояние проблемы, а также идеи использования комбинированных алгоритмов и дается постановка задачи.
Во второй главе обосновывается способ автоматического построения последовательных программ умножения и для последовательных программ, построенных по комбинированному методу, приводится метод выбора оптимального параметра — размерности базы рекурсии, осуществляемой при синтезе программы. и.
Идея схемной реализации метода Карацубы была предложена и обоснована в 1999;2006 г. г. в работах [4, 5] и заключается в том, что предварительное сложение и умножение фрагментов многочленов, как и последующая сборка этих результатов, осуществляется схемой без циклов. Автоматизация построения таких схем была рассмотрена в работах [22, 42]. В обеих работах на основе аналитического исследования метода строились таблицы операций и адресов операндов, а затем из таких таблиц извлекалась схема. Однако такой метод построения схем не учитывает необходимость оптимизации по объему занимаемой памяти и по размещению в ней результатов промежуточных вычислений, что делает его применение довольно ограниченным для получения последовательных программ умножения. В нашем случае под последовательной программой понимается последовательность элементарных арифметических действий и операций присваивания. Последовательная программа характеризуется размером этой последовательности, объемом используемой памяти и степенью локализации адресаций. По этим параметрам последовательные программы можно оптимизировать.
Важной особенностью таких последовательных программ является независимость хода выполнения от исходных данных. Эта особенность позволяет путем несложных преобразований, для каждого разряда многочлена — результата умножения, получить цепочку арифметических действий над разрядами операндов, производимых в процессе выполнения последовательной программы и приводящих к получению именно этого разряда результата. Такие цепочки не зависят от способов адресации, порядка выполнения не связанных между собой операций и задают математическую суть программы. Значит, в случае эквивалентности двух программ будут порождаться одинаковые цепочки действий по формированию разрядов результата умножения. Верно и обратное: в случае равенства цепочек действий программы эквивалентны.
Одним из способов автоматического построения последовательных программ умножения является использование аналитически выведенных формул и таблиц распределения адресов операндов в элементарных операциях.
Идея второго способа автоматического построения последовательного алгоритма заключается в использовании уже разработанной и проверенной рекурсивной программы умножения. В коде такой программы каждая операция выполнения элементарного действия заменяется операцией сохранения наименования действия и адресов операндов. Затем ищется минимальный адрес операнда и принимается за нулевой адрес рабочего массива последовательной программы, соответствующая поправка вносится во все адреса операндов. При данном подходе необходимо обоснование корректности получаемой в результате программы, поскольку сам подход содержит неочевидные преобразования.
Последовательные программы, получаемые в результате реализации описанных выше подходов, принадлежат одному классу, при этом отличаются только способами адресации элементов и методами распределения памяти для храненияпромежуточных результатов. Поэтому для доказательства их эквивалентности можно использовать метод сравнения цепочек действий. Таким образом, автоматизация проверки заключается в том, что при выбранных максимально возможных степенях сомножителей, задающих область применимости программ, автоматически синтезируются две последовательных программы с использованием первого подхода без оптимизации и второго подхода с оптимизацией. Затем из первой и второй программ извлекаются цепочки действий по формированию базовых слов результата умножения. Формальными преобразованиями производится автоматическое приведение цепочек к единому виду и их сравнение. При положительном заключении о равенстве цепочек формируется заключение о корректности обеих синтезированных программ.
В получаемых описанным выше методом последовательных программ умножения, можно выделить набор типовых операций, количество которых подчинено уточненной оценке:
Здесь ль П2 — параметры алгоритма, причем п = пхп2- количество разрядов в многочленах — операндах, а п2 — параметр перехода — размерность базы рекурсии, осуществляемой при синтезе программы. Таким образом, задача получения оценки сводится к задаче определения постоянных Сь. для ограниченного набора элементарных действий. Что делается путем построения семи последовательных программ и решением системы линейных уравнений для каждого выбранного действия. В итоге были получены следующие оценки:
1) для последовательных программ, построенных по комбинированному методу, использующему метод сдвигов и сложений с таблицей.
2) для последовательных программ, построенных по комбинированному методу, использующему метод сдвигов и сложений с условными переходами.
Р = С^З1022″ 1 п2 + С231оё2Щ п2 + Същп2 + С4 З1022″ 1 + С5пх + С6п2 + С7 С1) 310^'!1+4.
1022 И,.
2).
16 зб#з1о&" «+1.
16 2 8.
3).
3) для последовательных программ, построенных по методу Карацубы, использующему таблицу умножения для многочленов степени не выше 7.
31оы.11*и + 4 ^.
243 8.
Далее с использованием формул (2) и (3) функции сложности табулируются, и по таблицам определяется оптимальный параметр перехода для последовательных программ, построенных по комбинированным алгоритмам. В итоге было получено, что:
1) для последовательных программ, построенных по комбинированному методу, использующему метод сдвигов и сложений с таблицей оптимальный параметр перехода равен 32;
2) для последовательных программ, построенных по комбинированному методу, использующему метод сдвигов и сложений с условными переходами оптимальный параметр перехода равен 256.
Для обобщения данной оценки сложности операции умножения в СР (р"), где р>2, отмечается, что при реализации умножения над полями нечетной характеристики, возникает необходимость табулирования операций сложения и вычитания, для реализации которых в ОР (2п) использовались операции «исключающего или» (ХСЖ). То есть каждая операция «исключающего или» заменяется операциями адресации по таблицам сложения и вычитания в вР (р). Дополнительно с увеличением характеристики поля, последовательная программа неявно усложняется вследствие увеличения количества бинарных разрядов, которые отводятся на один элемент поля ОР (р). Повторением вычислений констант, описанных выше, было получено, что функция сложности по количеству операций адресации в последовательной программе умножения многочленов в СР (3П) имеет вид: р = 15 * З12″ 1 п22 +107,5 * З10®-2″ 1 п2 -106 * щп2 +12 Так как в указанной последовательной программе вместо всех арифметических операций использовались адресации по таблице, то для дальнейшего увеличения характеристики поля не требуется изменение алгоритма, достаточно лишь заменить таблицы выполнения элементарных арифметических действий. Иначе, если обозначить 8(и, р) — последовательную программу умножения многочленов степени п над полем СР (рп), то при условии использования таблиц операций с 8-разрядными словами и р<251, получим 8(п, 2)=8((п/8)4жГ2(рЛ)-1,р).
Полученная таким образом функция позволяет оценить сложность последовательной программы для любого значения параметра п и р<251. А ограничение р<251 определяется физическим ограничением машинной памяти на использование таблиц выполнения элементарных арифметических действий.
В третьей главе производится сравнение методов умножения по времени выполнения для колец многочленов и приводятся идеи оптимизации алгоритма за счет увеличения степени локализации промежуточных результатов вычислений.
Непосредственная реализация рекурсивного метода Карацубы при увеличении степени умножаемых многочленов влечет за собой значительное увеличение объема оперативной памяти, необходимой для хранения промежуточных результатов, что в свою очередь существенно увеличивает время, необходимое на выполнение операции умножения. Для минимизации затрат по памяти в работе предлагается модификация метода Карацубы, которая заключается во введении дополнительных операций по перемещению результатов вычисления таким образом, чтобы на каждом шаге метода перемножаемые многочлены располагались в оперативной памяти непосредственно друг за другом. Такое расположение позволяет записывать промежуточные результаты умножения многочленов меньшего порядка на место операндов, что исключает дублирование информации и позволяет более рационально использовать память за счет незначительного увеличения количества операций в алгоритме.
Т.к. данная модификация алгоритма не предполагает внесения изменений в арифметическую суть метода, то для проверки корректности реализации операции умножения модифицированным методом используется способ сравнения цепочек арифметических действий, описанный в главе 2. Результаты экспериментов показали, что наибольший выигрыш по скорости дает последовательная реализация модифицированного алгоритма Карацубы, совмещенного с использующим таблицу методом сдвигов и сложений, с параметром перехода равным 32.
В четвертой главе проводится сравнительный анализ эффективности алгоритмов умножения точки эллиптической кривой на константу с использованием аффинной и проективной систем представления операндов и усовершенствованных алгоритмов умножения, инвертирования и деления в конечных полях. Результаты экспериментов подтверждают рекомендации по применению проективных координат: использование ускоренного расширенного алгоритма Евклида существенно сокращает время деления в диапазоне степеней, характерных для криптографического применения, что, однако, оказывается недостаточным для отказа от применения проективных координат для выполнения основных операций в группе точек эллиптической кривой, позволяющих сократить число делений за счет увеличения общего числа операций умножения. При этом применение оптимизированной операции умножения многочленов для выполнения операций над точками эллиптической кривой позволило сократить время умножения точки на константу на 5−15% в зависимости от типа кривой и порядка многочлена, образующего поле, в котором производятся вычисления.
В заключении даны выводы по результатам исследования.
4.5 Выводы.
Проведенные в главе эксперименты подтвердили теоретическое предположение о целесообразности использования проективных координат для выполнения операции умножения точки эллиптической кривой на константу как для суперсингулярных эллиптических кривых, так и для несуперсингулярных, несмотря на использование быстрого алгоритма инвертирования.
Применение в алгоритмах на эллиптических кривых предложенных модификаций алгоритмов умножения в расширениях полей степени 256 и выше позволило сократить время выполнения операции умножения точки эллиптической кривой на константу на 5−15% в зависимости от типа кривой.
Заключение
.
1. В работе проведен анализ существующих методов умножения и деления в полях многочленов, обоснована целесообразность комбинирования метода умножения путем сдвигов и сложений и метода Карацубы с уточнением порога перехода, а также преимущество использования последовательных алгоритмов умножения перед рекурсивными и циклическими.
2. Разработан метод контролируемого автоматического синтеза последовательных программ умножения многочленов и метод аналитического сравнения производительности последовательных программ умножения, основанный на оценке количества выполняемых в методе элементарных операций и не требующий построения и выполнения программы.
3. Определены наилучшие условия применения последовательных программ умножения по методу А. Карацубы в процессе умножения многочленов высоких степеней над полями различных характеристик. Выбираемым параметром является глубина рекурсии при построении последовательной программы, она определяет ту максимальную степень многочленов, возникающих в процессе рекурсии, перемножение которых осуществляется одной из известных модификаций метода сдвигов и сложений, при достижении наилучшего алгоритма умножения этого типа в целом.
4. Произведено практическое сравнение методов умножения и инвертирования в поле многочленов, выявлены области рекомендуемого применения методов при умножении многочленов высоких степеней над полями различных характеристик.
5. Проведен сравнительный анализ эффективности алгоритмов умножения точки эллиптической кривой на константу с использованием аффинной и проективной систем представления операндов и усовершенствованных алгоритмов умножения, инвертирования и деления в конечных полях.
6. Экспериментально подтверждены рекомендации по применению проективных координат в операциях умножения точки эллиптической кривой на константу.
7. В процессе выполнения работы была создана библиотека классов на языке С++, реализующая основные алгебраические функции в полях многочленов и на эллиптических кривых. В связи с использованием более низкого уровня языка программирования и неуправляемого кода разработанная библиотека показала пятикратное увеличение производительности по сравнению с существующим аналогом «Алгебраический процессор» на исходных данных, характерных для современной эллиптической криптографии.