CC BY
21
ПРИКЛАДНАЯ ДИСКРЕТНАЯ МАТЕМАТИКА 2016 Вычислительные методы в дискретной математике №3(33)
ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ В ДИСКРЕТНОЙ МАТЕМАТИКЕ
УДК 512.54
ОБ ОДНОМ АЛГОРИТМЕ ВЫЧИСЛЕНИЯ ФУНКЦИЙ РОСТА В КОНЕЧНЫХ ДВУПОРОЖДЁННЫХ ГРУППАХ ПЕРИОДА 51
А. А. Кузнецов
Сибирский государственный аэрокосмический университет, г. Красноярск, Россия
Пусть B0(2, 5) = (ai,a2) —наибольшая конечная двупорождённая бернсайдова группа периода 5, порядок которой равен 534. Для каждого элемента данной
ai ao
группы существует уникальное коммутаторное представление вида ai1 ■ a22 ■ ■ ... ■ aQ434, где ai G Z5, i = 1,2,..., 34. Здесь a1 и a2 — порождающие элементы Bo(2, 5); аз,..., аз4 — коммутаторы, которые вычисляются рекурсивно через a1 и а2. Определим фактор-группу группы B0(2, 5) следующего вида: Bk = = Bo (2, 5)/(afc+i ,...,a34). Очевидно, что |Bk | = 5k. Предложен новый алгоритм, при помощи которого вычислены функции роста Bk относительно порождающих множеств {a1,a2} и {а1, а-1, а2, а-1} для k = 15, 16, 17. На основе полученных данных вычислены диаметр и средний диаметр соответствующих графов Кэли.
Ключевые слова: функция роста, группа Бернсайда, граф Кэли. DOI 10.17223/20710410/33/10
AN ALGORITHM FOR COMPUTATION OF THE GROWTH FUNCTIONS IN FINITE TWO-GENERATED GROUPS OF EXPONENT 5
A. A. Kuznetsov
Siberian State Aerospace University, Krasnoyarsk, Russia
E-mail: alex_kuznetsov80@mail.ru
Let B0(2, 5) = (a1,a2) be the largest 2-generator Burnside group of exponent 5. It has the order 534. There is a power commutator representation of B0(2, 5). In this case, every element of the group can be uniquely represented as aQ1 ■ aQo ■ ... ■ aQ|4, where ai G Z5, ai G B0(2, 5), i = 1, 2,..., 34. Here, a1 and a2 are generators of B0(2, 5), commutators a3,...,a34 are recursively defined by a1 and a2. We define Bk = B0(2, 5)/(ak+1,..., a34) as a quotient of B0(2, 5). It is clearly that |Bk| = 5k. A new algorithm for computing the growth function of Bk is created. Using this algorithm, we calculated the growth functions of Bk relative to generating sets {a1, a2} and {a1, a-1, a2, a-1} for k = 15, 16, 17.
Keywords: Burnside group, the growth function.
1 Работа поддержана грантом Президента РФ (проект №МД-3952.2015.9).
Введение
Одним из важных инструментов для определения строения группы является изучение её роста относительно фиксированного порождающего множества. Пусть G = (X). Шаром Ks радиуса s группы G будем называть множество всех её элементов, которые могут быть представлены в виде групповых слов в алфавите X длины не больше s. Для каждого целого неотрицательного s можно определить функцию F(s) роста группы, которая равна числу элементов группы G относительно X, представимых в виде несократимых групповых слов длины s. Таким образом,
F(0) = |Ко| = 1, F(s) = |Ks| - |Ks-i| при s e N.
Как правило, функцию роста конечной группы представляют в виде таблицы, в которую записывают ненулевые значения F(s).
К сожалению, вычисление функции роста большой конечной группы является хотя и разрешимой, но весьма сложной проблемой. Это связано с тем, что в общем случае задача по определению минимального слова элемента группы, как показали С. Ивен и О. Голдрейх [1], является NP-трудной. Таким образом, в наихудшем случае количество элементарных операций, которые необходимо выполнить для решения указанной задачи, представляет собой экспоненциальную функцию от |X |.
Отметим также, что при вычислении функции роста группы мы параллельно выясняем характеристики соответствующего графа Кэли, например диаметр и средний диаметр [2]. Пусть F(s0) > 0, но F(s0 + 1) = 0, тогда s0 является диаметром графа Кэли группы G в алфавите порождающих X, который будем обозначать Dx(G).
__1 so
Соответственно средний диаметр DX (G) равен —— s • F(s).
|G| s=0
Как известно, благодаря своим замечательным свойствам графы Кэли нашли применение при моделировании топологий многопроцессорных вычислительных систем (МВС)—суперкомпьютеров [3], а также дата-центров [4]. И вполне возможно, что в недалёком будущем понадобятся знания о сверхбольших графах, которые будут востребованы при проектировании распределённых систем, пиковая производительность которых будет достигать 1 экзафлопс и выше.
Одной из широко применяемых топологий МВС является k-мерный гиперкуб, который задаётся группой B(k, 2) — k-порождённой бернсайдовой группой периода 2. Группа B(k, 2) имеет простую структуру и равна прямому произведению k экземпляров циклической группы порядка 2. Обобщением гиперкуба является n-мерный тор, который порождается прямым произведением n экземпляров циклических подгрупп, порядки которых могут не совпадать. В работе [5] исследованы графы Кэли групп B(k, 3), т. е. групп периода 3, а также проведён сравнительный анализ данных графов с гиперкубами. Анализ выявил, что графы B(k, 3) обладают более предпочтительными характеристиками, чем B(k, 2). Это означает, что при парном сравнении графов B(k1, 3) и B(k2, 2) с примерно одинаковым количеством вершин первые будут иметь меньшие диаметры, средние диаметры и степени. В связи с этим представляет интерес задача по исследованию графов Кэли конечных бернсайдовых групп другого периода.
Пусть B0(2, 5) = (а1,а2) —максимальная конечная двупорождённая бернсайдова группа периода 5, порядок которой равен 534 [6]. Используя систему компьютерной алгебры GAP, нетрудно получить рс-представление (Power Commutator presentation) данной группы [3, 7]. В этом случае каждый элемент g e B0(2, 5) может быть однозначно записан в следующем виде:
g = а?1 • а£2 • ... • «ОТ, а e Z5, i =1, 2,..., 34.
Здесь a1 и a2 —порождающие элементы Во (2, 5); a3,...,a34 — коммутаторы, которые вычисляются рекурсивно через a1 и a2 [6].
Обозначим через Вк фактор-группу В0 (2, 5) следующего вида:
Вк = Во(2,5)/(afc+i,... ,аз4). Очевидно, что |Bk | = 5к и для всех g Е В к выполняется
g = aa1 ■ aa2 ■ ... ■ af.
Помимо прикладного интереса, существует и другая причина повышенного внимания исследователей к функциям роста Вк и особенно к В34 = В0(2, 5). Это объясняется тем, что полученная информация может оказаться полезной при решении открытой проблемы о конечности В (2, 5) — свободной двупорождённой бернсайдовой группы периода 5. Если В(2, 5) конечна, то В34 = В(2, 5). Однако в обозримом будущем вычислить функцию роста Вз4 едва ли возможно, поскольку количество её элементов очень велико: 534 = 582076609134674072265625 w 5 ■ 1023.
К. А. Филипповым в работе [8] вычислена функция роста группы В14 в алфавите A2 = {a1, a2}. Аналогичная задача для симметричного порождающего множества A4 = = {ai , a- , a2, a2 } решена Ч. Симсом [9]. В обоих случаях решение было получено при помощи длительных компьютерных вычислений. Отметим также, что Симс в своём алгоритме использовал элегантный трюк, основанный на применении к В14 группы автоморфизмов, сохраняющих длину группового слова. Это дало возможность свести задачу к рассмотрению только 1/16 части элементов группы. Подобный трюк для несимметричного порождающего множества менее эффективен, поскольку позволяет сократить количество элементов не более чем в 2 раза.
Дальнейшее изучение функций роста групп В к для k > 14 сталкивается с серьёзными вычислительными трудностями ввиду их больших порядков.
В данной работе представлен эффективный алгоритм для вычисления функций роста указанных групп, который даёт возможность преодолеть существующий барьер. Компьютерная реализация алгоритма позволила вычислить функции роста групп Вк для k = 15, 16 и 17.
В п. 1 рассматривается базовый алгоритм для вычисления функции роста конечной группы, а также исследуется его сложность. В п. 2 алгоритм модифицируется для исследования групп В к. В п. 3 приведены результаты компьютерных вычислений по модифицированному алгоритму.
1. Описание базового алгоритма
В качестве отправной точки далее представлен общий алгоритм A-I, вычисляющий функцию роста произвольной конечной группы G, заданной фиксированным порождающим множеством X.
Лемма 1. Алгоритм A-I корректен, т.е. он за конечное число шагов вычислит функцию роста любой конечной группы G, заданной фиксированным порождающим множеством X.
Доказательство. По построению A-I выражает каждый элемент группы G в виде группового слова наименьшей длины в алфавите X. После каждого прохода от п. 2 до п. 8 множество Ks представляет собой шар радиуса s группы G относительно X. Поскольку количество элементов группы конечно, через конечное число шагов, при некотором s = s0 + 1, все g Е G будут записаны в виде слов, длины которых меньше
Алгоритм 1. А-I
Вход: X = |ж1, х2,..., хт} —порождающее множество С Выход: функция роста ^(з) группы С 1: ^ := 0, ^(0) := 1, Ко := {е}, Р := Ко 2: 8 := 8 + 1 3: К := К5_1
4: Для всех X € X и р € Р 5: $ = X ■ р (или д = р ■ х) 6: Если д € К8, то К := К и {д} 7: Р := К - Кв-1 8: ^(з) := |Р| 9: Если ^(8) > 0, то
переход в п. 2, 10: иначе
11: s0 := s — 1, переход в п. 12
__1 so
12: Dx(G) := so, Dx(G) := -— E s ■ F(s), выход
|Kso 1 s=0
s0 + 1. Получим F(s0 + 1) = 0 и |Ks01 = |G|, что означает остановку алгоритма. При этом __1 so
Dx(G) = s0 и Dx(G) = -— E s ■ F(s). ■
|Kso | s=0
Несмотря на то, что алгоритм A-I хорошо известен специалистам, автору не удалось найти публикаций, в которых были бы вычислены оценки его временной сложности. Например, в системе компьютерной алгебры GAP реализован указанный алгоритм, вычисляющий функцию роста группы подстановок. Однако в документации к системе отсутствуют сведения о его сложности. Для устранения указанного пробела воспользуемся машиной с произвольным доступом к памяти — RAM-машиной, а также асимптотическим анализом сложности [10]. Далее термины «временная сложность» и «сложность» будем считать равнозначными; кроме этого, введём обозначения:
— T(|G|) — сложность алгоритма A-I;
— O(/ (|G|)) —верхняя асимптотическая оценка сложности;
— П(/(|G|)) —нижняя асимптотическая оценка сложности;
— в(/(|G|)) —одновременно верхняя и нижняя оценки сложности.
Нам понадобятся следующие вспомогательные утверждения.
Лемма 2. В процессе работы алгоритма A-I необходимо вычислить и проверить |X| ■ |G| групповых слов.
Доказательство. При вычислении слов длины s ^ 1 следует обработать F(s — 1)|X| элементов. В итоге получим
Dx (G)+1 Dx (G)+1
E F(s — 1)|X| = |X| E F(s — 1) = |X| ■ |G|.
s=1 s=1
Лемма доказана. ■
Лемма 3. Пусть N(|G|) —суммарное количество проверок вида g Е Ks алгоритма A-I. Тогда верно следующее неравенство:
2|G|2 + (|X| - 3) |G| + 1 ^ N(|G|) ^ (jX| - 2) |G|2 + 2|G|.
Доказательство. Рассмотрим сначала те групповые слова, которые войдут в Ks. Всего таких элементов, не считая единицы, (|G| — 1). Для первого из них необходимо осуществить одну проверку (с единицей группы), для второго — две и т. д. В итоге получим 1 + 2 + ... + (|G| — 1) = |G|(|G| — 1)/2 операций.
Теперь рассмотрим те слова, которые не войдут в Ks. Согласно лемме 2, количество таких элементов равно |X| ■ |G| — (|G| — 1). В лучшем случае каждый из этих элементов будет проверен один раз и исключён, а в худшем — для каждого из них потребуется осуществить |G| сравнений.
Таким образом, общее число проверок в худшем и лучшем случаях равно |G|(|G| — 1)/2 + |X| ■ |G| — (|G| — 1) и |G|(|G| — 1)/2 + (|X| ■ |G| — (|G| — 1))|G| соответственно. После элементарных преобразований получим нужное неравенство. ■
Теорема 1. T(|G|) Е ft(|G|2) и T(|G|) Е O(|X| ■ |G|2).
Доказательство. При |G| ^ то сложностью операций всех пунктов алгоритма A-I, кроме четвертого, следует пренебречь. Поскольку количество элементарных операций, необходимое для вычисления произведения двух элементов группы, конечно и ограничено ©(1), то суммарная сложность M(|G|) п. 5 ограничена ©(|X| ■ |G|). Согласно лемме 3, получим |G|2/2 ^ N(|G|) ^ N(G) Е ft(|G|2) и N(|G|) ^ |X| ■ |G|2 ^ ^ N(|G|) Е O(|X| ■ |G|2). Так как T(|G|) = M(|G|) + N(|G|), получим T(|G|) Е ft(|G|2) и T(|G|) Е O(|X| ■ |G|2). ■
Следствие 1. Если |X| < |G|, то T(|G|) Е ©(|G|2).
Следствие 2. Если |X| ~ |G|, то T(|G|) Е O(|G|3).
Если порядок группы G достаточно большой, то наиболее критическими участками алгоритма A-I будут пп. 5 и 6, т.е. умножение элементов группы и проверка, позволяющая определить, встречался ранее такой элемент или нет.
2. Модификация алгоритма A-I для исследования групп Вк
В данном случае X = A2 или A4. Рассмотрим п. 5 алгоритма.
Существуют два способа умножения элементов в группах В к: собирательный процесс [3, 7] и метод полиномов Холла [11]. Вычислительные эксперименты показали, что второй способ позволяет значительно быстрее собирательного процесса (по крайней мере на порядок) умножать элементы в группах Вк [12].
Следующий вопрос: как эффективнее умножать порождающие x Е X на элементы группы p Е Вк — слева (1 способ) или справа (2 способ), т. е. x■ p или p■ x? Кроме этого, была предложена идея при умножении слева использовать пять вариантов полиномов в зависимости от значения степени а1 элемента p (3 способ). Для ответа на данный вопрос необходимо вычислить суммарное число операций t(k) в полиномах Холла, необходимое для выполнения одного произведения. В полиномах имеются три вида действий: сложение, умножение и остаток от деления числа на 5. На рис. 1 представлены графики t(k) для каждого из трёх способов (в предположении, что указанные действия имеют одинаковую сложность равную 1).
Рис. 1. Графики сложности t(k) для одного произведения
Рис. 1 наглядно показывает, что третий способ самый эффективный. Этот факт был учтён при реализации п. 5 алгоритма.
Теперь обратимся к п. 6 алгоритма A-I. Для его оптимизации потребуется каждому
элементу g £ Bk присвоить уникальный порядковый номер. Пусть g = а^1 произвольный элемент Bk. Определим биективное отображение ф следующего вида:
Ф
g <—> ng = «i ... .
Здесь ng представляет собой целое неотрицательное число, записанное в пятеричной системе счисления, которое возьмём в качестве порядкового номера g. Легко увидеть, что ng пробегает все значения от 0 до (5k — 1).
Модифицируем A-I следующим образом. В п. 1 добавим булев вектор V размерности 5k, все элементы которого первоначально равны 0. Для удобства индексация элементов V начинается с 0, т.е. последний элемент вектора имеет индекс (5k — 1). Заменим п. 6 алгоритма A-I на следующий: 6: Если Vng = 0, то := 1 и Ks := Ks U {g}.
Данная модификация позволяет избежать трудоёмкой процедуры поиска элемента g в множестве Ks. Будем считать, что для выполнения п. 6 алгоритма требуется не более t(g) элементарных операций.
Несмотря на то, что порядки Bk ограничены сверху, имеет смысл рассмотреть асимптотическую оценку сложности модифицированной версии алгоритма, поскольку данные группы содержат достаточно большое число элементов при k > 8.
Теорема 2. Пусть T(|Bk|) —вычислительная сложность модифицированного алгоритма A-I. Тогда верно, что T(|Bk|) £ 0(|Bk|).
Доказательство. При вычислении асимптотической оценки следует учитывать сложность только наиболее трудоёмких фрагментов алгоритма, т. е. пп. 5 и 6. Поскольку T(|G|) = (t(k) + t(g))|X| • |Bk|, |X| << |Bk|, t(k) << |Bk| и t(g) << |Bk|, получим T(G) £ e(|Bfc|). ■
Также отметим, что в п. 5 все элементы g вычисляются независимо друг от друга, поэтому этот участок алгоритма можно легко распараллелить. В этом случае сначала параллельно вычисляются все произведения g, а затем для каждого элемента получившегося массива последовательно выполняется п. 6.
,«2
ак
Легко увидеть, что перечисленные модификации алгоритма А-1 применительно к группам Вк значительно повышают его производительность.
3. Компьютерная реализация алгоритма
Модифицированный алгоритм А-1 реализован на языке С++ и апробирован на 4-процессорном компьютере (на каждом по 16 ядер) с общей памятью объёмом 256 Гбайт. В результате вычислены функции роста групп Вк для к = 15, 16 и 17 для порождающих множеств А2 = {«1,0.2} и А4 = {а,!,«-1 , а^а-1}. Затраченное на расчёты время примерно равно 1 час для к = 15, 5 часов — для к =16 и около суток — для к =17. На рис. 2-7 представлены графики функций роста указанных групп.
В таблице приведены характеристики графов Кэли групп Вк, вычисленные в настоящей работе, а также, для полноты картины, уже известные. Значения средних диаметров округлены до ближайших целых.
Отметим, что для решения задачи маршрутизации в графе Кэли, после того как будут получены минимальные слова группы в результате работы алгоритма А-1, можно будет воспользоваться алгоритмом из [4], который имеет сложность 0(^Х(£)).
Длина
Рис. 2. График функции роста группы В15 в алфавите А2
Длина
Рис. 4. График функции роста группы В17 в алфавите А
О 5 10 15 20 25 30 35
Длина
Рис. 5. График функции роста группы В15 в алфавите А4
О 5 10 15 20 25 30 35
Длина
1.5 -
о
ь
о ° 1 О- 1 -
0.5 -
0 5 10 15 20 25 30 35
Длина
Рис. 7. График функции роста группы В17 в алфавите А4
Характеристики графов Кэли групп Bk
k 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
DA2 (Bk) 8 10 13 20 20 25 30 31 32 35 40 40 45 46 50 55
Da2 (Bk) 4 6 8 12 14 16 20 22 24 27 29 31 34 36 39 42
DA4 (Bk) 4 6 9 10 13 15 19 20 20 24 26 28 30 30 33 34
Da4 (Bk) 2 4 5 7 9 10 13 14 15 17 19 21 22 23 26 27
ЛИТЕРАТУРА
1. EvenS. and Goldreich O. The Minimum Length Generator Sequence is NP-hard // J. Algorithms. 1981. No. 2. P. 311-313.
2. Кузнецов А. А., Кузнецова А. С. Параллельный алгоритм для исследования графов Кэли групп подстановок // Вестник СибГАУ. 2014. №1. C. 34-39.
3. Holt D., EickB., and O'Brien E. Handbook of Computational Group Theory. Boca Raton: Chapman & Hall/CRC Press, 2005. 514 p.
4. Camelo M, Papadimitriou D., Fabrega L., and Vila P. Efficient routing in data center with underlying Cayley graph // Proc. 5th Workshop on Complex Networks CompleNet. 2014. P. 189-197.
5. Кузнецов А. А. Графы Кэли бернсайдовых групп периода 3 // Сибирские электронные математические известия. 2015. T. 12. C. 248-254.
6. Havas G., Wall G., and Wamsley J. The two generator restricted Burnside group of exponent five // Bull. Austral. Math. Soc. 1974. No. 10. P. 459-470.
7. Sims C. Fast multiplication and growth in groups // Proc. 1998 Intern. Symp. Symbolic Algebraic Computation. 1998. P. 165-170.
8. Филиппов К.А. О диаметре Кэли одной подгруппы группы B0(2, 5) // Вестник СибГАУ. 2012. №1. С. 234-236.
9. Sims C. Computation with Finitely Presented Groups. Cambridge: Cambridge University Press, 1994. 628 p.
10. Скиена C. Алгоритмы. Руководство по разработке. СПб.: БХВ-Петербург, 2014. 720 с.
11. HallP. Nilpotent Groups: Notes of Lectures Given at the Canadian Mathematical Congress Summer Seminar, University of Alberta, 12-30 August, 1957. London: Queen Mary College, 1969.
12. Кузнецов А. А., Кузнецова А. С. Быстрое умножение элементов в конечных двупорождён-ных группах периода пять // Прикладная дискретная математика. 2013. № 1. C. 110-116.
REFERENCES
1. Even S. and Goldreich O. The Minimum Length Generator Sequence is NP-hard. J. Algorithms, 1981, no. 2, pp. 311-313.
2. Kuznetsov A. A. and Kuznetsova A. S. Parallel'nyy algoritm dlya issledovaniya grafov Keli grupp podstanovok [A parallel algorithm for study of the Cayley graphs of permutetion groups]. Vestnik SibGAU, 2014, no. 1, pp. 34-39. (in Russian)
3. HoltD., EickB., and O'Brien E. Handbook of Computational Group Theory. Boca Raton, Chapman & Hall/CRC Press, 2005. 514 p.
4. Camelo M., Papadimitriou D, Fabrega L., and Vila P. Efficient routing in data center with underlying Cayley graph. Proc. 5th Workshop on Complex Networks CompleNet, 2014, pp. 189-197.
5. Kuznetsov A. A. Grafy Keli bernsaydovykh grupp perioda 3 [The Cayley graphs of Burnside groups of exponent 3]. Siberian Electronic Math. Reports, 2015, vol.12, pp.248-254. (in Russian)
6. Havas G., Wall G., and Wamsley J. The two generator restricted Burnside group of exponent five. Bull. Austral. Math. Soc., 1974, no. 10, pp. 459-470.
7. Sims C. Fast multiplication and growth in groups. Proc. 1998 Intern. Symp. Symbolic Algebraic Computation, 1998, pp. 165-170.
8. Filippov K.A. O diametre Keli odnoy podgruppy gruppy B0(2, 5) [About Cayley diameter of one subgroup of the group B0(2, 5)]. Vestnik SibGAU, 2012, no. 1, pp. 234-236. (in Russian)
9. Sims C. Computation with Finitely Presented Groups. Cambridge, Cambridge University Press, 1994. 628 p.
10. Skiena C. Algoritmy. Rukovodstvo po razrabotke [Algorithms. Guide for the Development]. SPb., BKhV-Petersburg Publ., 2014. (in Russian)
11. Hall P. Nilpotent Groups: Notes of Lectures Given at the Canadian Mathematical Congress Summer Seminar, University of Alberta, 12-30 August, 1957. London, Queen Mary College, 1969.
12. Kuznetsov A. A. and Kuznetsova A. S. Bystroe umnozhenie elementov v konechnykh dvuporozhdennykh gruppakh perioda pyat' [Fast multiplication in finite two-generated groups of exponent five]. Prikladnaya Diskretnaya Matematika, 2013, no. 1, pp. 110-116. (in Russian)
