Научная статья на тему 'Быстрое умножение элементов в конечных двупорождённых группах периода пять'

Быстрое умножение элементов в конечных двупорождённых группах периода пять Текст научной статьи по специальности «Математика»

CC BY
150
44
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ПЕРИОДИЧЕСКАЯ ГРУППА / СОБИРАТЕЛЬНЫЙ ПРОЦЕСС / ПОЛИНОМЫ ХОЛЛА / HALL''S POLYNOMIALS / PERIODICAL GROUPS / GATHERING PROCESS

Аннотация научной статьи по математике, автор научной работы — Кузнецов Александр Алексеевич, Кузнецова Александра Сергеевна

Представлен алгоритм, основанный на полиномах Холла, для быстрого умножения элементов в конечных двупорождённых группах периода пять.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Fast multiplication in finite two-generated groups of exponent five

An algorithm for fast multiplication of elements in finite two-generated groups of exponent five based on Hall's polynomials is described.

Текст научной работы на тему «Быстрое умножение элементов в конечных двупорождённых группах периода пять»

2013 Вычислительные методы в дискретной математике №1(19)

УДК 512.54

БЫСТРОЕ УМНОЖЕНИЕ ЭЛЕМЕНТОВ В КОНЕЧНЫХ ДВУПОРОЖДЁННЫХ ГРУППАХ ПЕРИОДА ПЯТЬ1

А. А. Кузнецов, А. С. Кузнецова Сибирский государственный аэрокосмический университет, г. Красноярск, Россия

E-mail: [email protected]

Представлен алгоритм, основанный на полиномах Холла, для быстрого умножения элементов в конечных двупорождённых группах периода пять.

Ключевые слова: периодическая группа, собирательный процесс, полиномы Холла.

Введение

Пусть В0(2,5,k) —максимальная конечная двупорождённая бернсайдова группа периода 5 ступени нильпотентности k. В данном классе групп наибольшей является группа B0(2, 5,12), порядок которой равен 534 [1]. Для каждой из B0(2, 5, k) известны PC-представления (Power Commutator presentation [2, 3]), которые несложно получить, используя систему компьютерной алгебры GAP. Для исследования различных свойств B0(2, 5,k), например таких, как функция роста или диаметр Кэли группы, часто требуется вычислять произведения её элементов. Пусть n^1 ... „П" и n^1 ... пП" — два произвольных элемента в группе В0(2, 5,k), записанных в коммутаторном виде. Тогда их произведение равно

„Ж1 „Ж" _ ау1 Пуп _ „z1 nzn

Cv 1 • • • Ci 1 • • • Ci 1 • • • *

Основой для нахождения коэффициентов zi является собирательный процесс [2, 3], который реализован в системах компьютерной алгебры GAP и MAGMA. Однако это не единственный способ для вычисления произведений элементов группы. В работе Ф. Холла [4] показано, что zi представляют собой полиномиальные функции (в нашем случае над полем Z5), зависящие от переменных xi,... , xi, y1,... , yi, которые называют сейчас полиномами Холла. Согласно [4],

Zi _ Xi + yi + pi(x 1, . . .,Xi-i,yi,... ,yi-1).

Ч. Симс в [5] анонсировал, что вычислил указанные полиномиальные функции, однако они так и не были нигде опубликованы. Кроме того, авторам не удалось найти компьютерные программы, в которых реализован метод полиномов Холла для вычисления произведений элементов группы. Для восполнения указанного пробела в настоящей работе приведено подробное описание алгоритма для вычисления полиномов Холла в группах В0(2,5, k). Полиномы zi найдены при помощи машинных вычислений в системе компьютерной алгебры MATLAB, после чего полученные формулы реализованы в виде отдельной программы на языке Си. Сравнение скорости вычисления около 104 произведений случайно выбранных элементов в каждой из групп В0(2, 5, k)

1Работа выполнена при финансовой поддержке гранта РФФИ (код проекта № 10-01-00509-а).

показало, что метод полиномов Холла имеет значительное преимущество перед собирательным процессом. В среднем использование полиномов позволяет получить произведение двух элементов группы на порядок (т. е. в 10 раз) быстрее собирательного процесса.

При описании алгоритма ограничимся случаями к ^ 4. Для к > 4 полиномы вычисляются аналогично, однако их вывод занимает значительно больше места, что делает невозможным проверить доказательство «вручную». Полиномы для наибольшего случая к =12 содержат также полиномы для всех других к и занимают 42 страницы текста, вследствие чего в работе не приводятся. Их можно легко получить, связавшись с авторами по электронной почте.

Из описания алгоритма, приведённого далее, следует

Теорема 1. Пусть а^1 ... аПп и а^1 ... аП —два произвольных элемента в группе Во(2, 5, к), записанных в коммутаторном виде, где к Е N и к ^ 4. Тогда их произведение равно аХ1 ... аПп ■ ау1 ... аПп = а^1 ... аП, где Хг Е Ъ5 — полиномы Холла, задаваемые формулами (1), (2) при к =1, (1)-(3) при к = 2, (1)-(5) при к = 3 и (1)-(8) при к = 4:

XI = XI + уь (1)

х2 = х2 + У2; (2)

хз = хз + уз + Х2У1; (3)

Х4 = Х4 + у4 + + хз У1; (4)

х5 = х5 + У5 + х2У1У2 + ^^ У1 + Х3У2; (5)

Хб = Хб + Уб + Ж2^ *) + *) + Ж4У1; (6)

х7 = х7 + У7 + **(*) У2 + (*) + Х3У1У2 + Х4У2 + Х5Уъ (7)

Х8 = х8 + У8 + Х2У^^ У1У2 + У1 + хз ^+ х5У2. (8)

Описание алгоритма

Вычислим полиномы Холла для группы В0(2, 5, 4). Работая с данной группой, автоматически получим полиномы для групп В0(2, 5,1), В0(2,5, 2) и В0(2, 5, 3), поэтому нет необходимости отдельно рассматривать эти случаи. При к < 4 коммутаторы, вес которых больше к, не рассматриваются, поскольку они равны единице группы по определению.

Запишем РС-представление группы В0(2, 5,4) согласно [1].

Коммутаторы веса 1:

а1 , а2 — образующие группы.

Коммутаторы веса 2: Коммутаторы веса 3:

Коммутаторы веса 4:

аз = [а2, а1].

а4 = [аз,а1] = [а2,а1,а1], а5 = [аз, а2] = [а2,а1,а2].

аб = [а4,а1] = [а2, а1,а1, а1],

07 — [05,01] — [02,01,02,01],

08 — [05,02] — [02,01,02,02].

Список определяющих соотношений Я для базисных коммутаторов:

05 — 1 (1 ^ г ^ 8), [02,01] — 03, [03, 01] — 04, [03, 02] — 05, [04, 01] — 06,

[04, 02] — 07, [04, 03] — 1, [05, 01] — 07, [05, 02] — 08, [05, 03] — 1, [05, 04] — 1,

[06.01] — 1, [06,02] — 1, [06,03] — 1, [06,04] — 1, [06,05] — 1, [07,01] — 1,

[07.02] — 1, [07,03] — 1, [07,04] — 1, [07,05] — 1, [07,06] — 1, [08,01] — 1,

[08.02] — 1, [08,03] — 1, [08,04] — 1, [08,05] — 1, [08,06] — 1, [08,07] — 1.

Таким образом,

Во(2, 5,4) — (01,02,03,04,05,06,07,08 | Я}.

Каждый элемент группы выражается единственным образом в виде нормального коммутаторного слова:

Уд е В0(2, 5, 4) д — 0Х10Х20Х30Х40^®0Х60Х70Х8, х € ^5.

Иногда будем писать д — (х1,... , х8).

Для того чтобы определить функции ¿¿, сначала необходимо вычислить произведения вида 0у 0Х для всех 1 ^ г < 7 ^ 8, х, у — 1, 2, 3, 4. Для пары (7, г) требуется по 16 значениям произведений (у, х) найти интерполяционный полином для каждого из 8 коммутаторов.

Начнем с первого случая 0^ 0Х:

о2о1 = (1,1,1, 0,0, 0, 0, 0), о2о1 — (2,1, 2,1, 0, 0, 0, 0),

о2о1 = (3,1, 3, 3,0,1, 0, 0), о2о1 = (4,1, 4,1, 0, 4, 0, 0),

о2о1 — (1, 2, 2, 0,1, 0, 0, 0), 2 2 о2о: = (2, 2, 4, 2, 2, 0,1, 0),

о2о1 = (3, 2,1,1, 3, 2, 3, 0), о2о1 = (4, 2, 3, 2, 4, 3,1, 0),

о2о} = (1, 3, 3, 0, 3, 0, 0,1), о2о1 = (2, 3,1, 3,1, 0, 3, 2),

о22о'1 = (3, 3, 4, 4,4, 3, 4, 3), о^о^ = (4, 3, 2, 3, 2, 2, 3, 4),

о2о} = (1, 4, 4, 0,1, 0, 0, 4), 42 о2о: = (2, 4, 3, 4, 2, 0,1, 3),

о2о1 = (3, 4, 2, 2, 3, 4, 3, 2), 44 о2о: = (4, 4, 1, 4, 4, 1, 1, 1).

Запишем

о| о?

пх оУ о/з1’2)(х’У) 0/i1,2)(x>y)

О^ П2 П3 П4

... о:

/81,2)(х>у) '8 ,

4 4

где /г(1,2) (х, у) — Е Е водхру9 —некоторые полиномы над полем Z5. Для их определе-р=19=1

ния выполним интерполяцию для каждого коммутатора г — 3, 4,..., 8.

Для нахождения /г(1’2)(х,у) требуется решить систему линейных уравнений над указанным полем:

Е Е хру9 — ХГ, х,у — 1, 2, 3, 4, (9)

р=19=1

где — значение г-го коммутатора для пары (у, х). Данная система имеет 16 пере-

менных и состоит из 16 уравнений.

Покажем на примере 8-го коммутатора нахождение /8^2)(ж, у). Для краткости вместо в8 будем писать врд• Подставив в систему (9) все значения г|Х, получим

1 1 1

2 4 2

3 4 3

4 1 4 224 433 1 3 2 321 334 1 2 3 422 231 441 312 213 1 4 4

111

342

243

414

243

111

414

342

342

414

111

243

414

243

342

111

111 1 3 4 1 2 4 141 243 222 232 213 342 323 333 312 414 4 3 1 4 2 1 444

111

213

312 414 143 244 341

442 142 241 344

443 1 1 4 212

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

313 4 1 1

11

41

41 11 13 43

43

13 12

42 42 12

14

44 44 14

11

31

21

41

11

31

21

41

11

31

21

41

11

31

21

41

/ви\ /ОХ

в21 0

в12 О

в31 0

в22 0

в13 О

в41 0

в32 0

в23 1

в14 2

в42 3

в33 4

в24 4

в43 3

в34 2

^в44 ) 1

Ранг матрицы системы равен 16, поэтому она имеет единственное решение ви — 2, в12 — 2, в13 — 1, все остальные коэффициенты равны нулю. Следовательно,

/81,2) (^, у) — 2ху + 2ху2 + ху3.

Аналогичным образом вычисляются все полиномы /г(1,2)(х, у). Перечислим их:

/1(1,2)(х,у) — x,

/21,2)(х,у) — У,

/31,2)(х,у) — ^

С(1,2) ,

Л ’ Чж,у) = 2жу + 3ж2у = ^2 )у,

/51’2)(ж,у) = 2жу + 3жу2 = ж

/61’2) (ж У) = 2жУ + 2ж2у + ж3у = У

/71,2) (ж, у) = 4жу + жу2 + ж2у + 4ж2у2 =

/81’2) (ж, у) = 2жу + 2жу2 + жу3 = ж ^

Таким образом,

ж2 у2

< = ( ж,у,жуЛ ж)у,42

Вторая пара ау аХ:

а3а} = (1, 0,1,1,0, 0, 0, 0), а^а^ — (2, 0, 1, 2, 0, 1, 0, 0

а3а} = (3, 0,1, 3,0, 3, 0, 0), а^а} — (4, 0, 1, 4, 0, 1, 0, 0

а|а} = (1, 0, 2, 2,0, 0, 0, 0), 2 2 а3а} — (2, 0, 2, 4, 0, 2, 0, 0

а3а} — (3, 0, 2,1,0,1, 0, 0), а3а} — (4, 0, 2, 3, 0, 2, 0, 0

а^а} = (1, 0, 3, 3,0, 0, 0, 0), а3а} — (2, 0, 3, 1, 0, 3, 0, 0

а3а} — (3, 0, 3, 4,0, 4, 0, 0), а3а} — (4, 0, 3, 2, 0, 3, 0, 0

а^а} = (1, 0, 4, 4,0, 0, 0, 0), 4 2 а3а} — (2, 0, 4, 3, 0, 4, 0, 0

а^а} = (3, 0, 4, 2,0, 2, 0, 0), 4 4 а3а} — (4, 0, 4, 1, 0, 4, 0, 0

а| аХ = ( х, 0,у, ху, 0, 2ху + 3х2у, 0,0).

Третья пара ау аХ:

а3 а2 — (0, 1, 1, 0,1, 0, 0, 0) а3а2 2 0 1 2 (0 0, 0,1),

а3 а32 — (0, 3, 1, 0, 3, 0, 0, 3) а3а2 — (0 ,4 , 1 ,0 ,4 , 0, 0,1),

а со ю а2 — (0, 1, 2, 0, 2, 0, 0, 0) а со ю а ю ю — (0 ,2 ,2 ,0 ,4 , 0, 0, 2),

а со ю а32 — (0, 3, 2, 0,1, 0, 0,1) 24 а3а2 — (0, 4, 2, 0, 3, 0, 0, 2),

а со со а2 — (0, 1, 3, 0, 3, 0, 0, 0) а со со а ю ю — (0, 2, 3, 0,1, 0, 0, 3),

а со со а32 — (0, 3, 3, 4) 0, 0, 4, 0, а3а2 2, 0, 3, 4, (0, 0,0, 3),

а3 а2 = (0, 1, 4, 0, 4, 0, 0, 0) 22 а 43 а 3, 0, 4, 2, (0, 0,0, 4),

а3 а32 — (0, 3, 4, 0, 2, 0, 0, 2) 44 а 3а 2 1 0, 4, 4, (0, 0, 0, 4);

У а3 а Х 2— = (0,х^ 0 ху , 0 , 0 , 2ху + 3х2у).

Четвёртая пара а4 аХ:

а (1, 0, 0, 1,0,1, 0, 0), а4а} — (2, 0, 0, 1, 0, 2, 0, 0),

а4 а: — (3, 0, 0, 1,0, 3, 0, 0), а4а} — (4, 0, 0, 1, 0, 4, 0, 0),

а4 а: — (1, 0, 0, 2, 0, 2, 0, 0), 2 2 а4а} — (2, 0, 0, 2, 0, 4, 0, 0),

а4а} — (3, 0, 0, 2,0,1, 0, 0), 24 а4а} — (4, 0, 0, 2, 0, 3, 0, 0),

а4а} — (1, 0, 0, 3, 0, 3, 0, 0), а4а} — (2, 0, 0, 3, 0, 1,0, 0),

а4а} — (3, 0, 0, 3, 0, 4, 0, 0), а4а} — (4, 0, 0, 3, 0, 2,0, 0),

а4а} — (1, 0, 0, 4, 0, 4, 0, 0), 42 а4а} — (2, 0, 0, 4, 0, 3,0, 0),

а4а} — (3, 0, 0, 4, 0, 2, 0, 0), 44 а4а} — (4, 0, 0, 4, 0, 1,0, 0);

а 4 аХ — (х, 0, 0, у, 0, ху, 0, 0)

Пятая пара а4 af •

а4 а2 = = (0, 1, 0,1,0, 0, 1, 0) а4а2 — (0, 2, 0, 1, 0, 0, 2, 0

а4 а2 = = (0, 3, 0,1,0, 0, 3, 0) а4а2 — (0, 4, 0, 1, 0, 0, 4, 0

a24 а2 = = (0, 1, 0, 2,0, 0, 2, 0) 22 а 24 а 22 — (0, 2, 0, 2, 0, 0, 4, 0

а4 а32 — = (0, 3, 0, 2,0, 0, 1, 0) 24 а 24 а 42 — (0, 4, 0, 2, 0, 0, 3, 0

а4 а2 = = (0, 1, 0, 3, 0, 0, 3, 0) а 34 а 22 — (0, 2, 0, 3, 0, 0, 1, 0

а4 а32 — (0, 3, 0, 3, 0, 0, 4, 0) а4а2 — (0, 4, 0, 3, 0, 0, 2, 0

а4 а2 = (0, 1, 0, 4,0, 0, 4, 0) 42 а4а2 — (0, 2, 0, 4, 0, 0, 3, 0

а4 а32 — (0, 3, 0, 4,0, 0, 2, 0) 44 а4а2 — (0, 4, 0, 4, 0, 0, 1, 0

а4 af — (0, x 0, у, 0, 0, xy, 0)

X У

aX аЛ

Шестая пара аУ af. Коммутаторы а4 и а3 перестановочны, поэтому a4 af — а3 а4

.X h

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Теперь рассмотрим a| af:

а^ а^ — (1, 0, 0, 0,1, 0, 1, 0) а^а і — (2, 0, 0, 0, 1, 0, 2, 0),

а5а} — (3, 0, 0, 0,1, 0, 3, 0) а^а і — (4, 0, 0, 0, 1, 0,4, 0),

а^ а^ — (1, 0, 0, 0, 2, 0, 2, 0) 22 а^а^ — (2, 0, 0, 0, 2, 0,4, 0),

а^а! — (3, 0, 0, 0, 2, 0, 1, 0) 24 а5а: — (4, 0, 0, 0, 2, 0,3, 0),

а^а} — (1, 0, 0, 0, 3, 0, 3, 0) а^а} — (2, 0, 0, 0, 3, 0,1, 0),

33 а^а і — (3, 0, 0, 0, 3, 0, 4, 0) а^а^ — (4, 0, 0, 0, 3, 0,2, 0),

а^а } — (1, 0, 0, 0, 4, 0, 4, 0) 42 а5а: — (2, 0, 0, 0, 4, 0, 3, 0),

а^ а і — (3, 0, 0, 0, 4, 0, 2, 0) 44 а^а^ — (4, 0, 0, 0, 4, 0,1, 0);

а5 af — (x, 0, 0, 0,у, 0, xy, 0)

И наконец, последняя неперестановочная пара a| af:

ю а а2 — (0, 1, 0, 0,1, 0, 0,1) 22 а ю а 2, (0, 0, 0, ,2 0, 0, 1

а5 а32 — (0, 3, 0, 0,1, 0, 0, 3) а5а2 — (0, 4, 0, 0, 1, 0, 0, 4

а5 а2 — (0, 1, 0, 0, 2, 0, 0, 2) 22 а5а2 — (0, 2, 0, 0, 2, 0, 0, 4)

а5 а32 — (0, 3, 0, 0, 2, 0, 0,1) 24 а5а42 — (0, 4, 0, 0, 2, 0, 0, 3)

а5 а2= (0, 1, 0, 0, 3, 0, 0, 3) а5а2— (0, 2, 0, 0, 3, 0, 0, 1)

а5 а32 — (0, 3, 0, 0, 3, 0, 0, 4) а5а2 = (0, 4, 0, 0, 3, 0, 0, 2)

а5 а2 — (0, 1, 0, 0, 4, 0, 0, 4) 42 а5а2 — (0, 2, 0, 0, 4, 0, 0, 3)

а5 а32 — (0, 3, 0, 0, 4, 0, 0, 2) 44 а5а2 — (0, 4, 0, 0, 4, 0, 0, 1)

а5 af — (0, x, 0, 0,у, 0, 0, xy)

Все оставшиеся пары ауУ а? перестановочны, т. е. ауУ а? = а? ау.

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

ух x У /}+і(х>у) /)+,2)(х>у)

aj aj — aj aj a.^ а^+Т2

/8i,j)(x>y)

Пользуясь (10), можем вычислить произведение а^1 ... а^8 ■ ау1 После выполнения данной процедуры найдём все ^ (1)-(8). Теорема доказана.

.. aj8 — а}1

(10)

.aj8.

8

ЛИТЕРАТУРА

1. 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.

2. Sims C. Computation with Finitely Presented Groups. Cambridge: Cambridge University Press, 1994. 628 p.

3. HoltD., EickB., and O’Brien E. Handbook of computational group theory. Boca Raton: Chapman & Hall/CRC Press, 2005. 514 p.

4. 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.

5. Sims C. Fast multiplication and growth in groups // Proc. Intern. Symp. Symbolic and Algebraic Computation. New York, NY, USA, 1998. P. 165-170.

i Надоели баннеры? Вы всегда можете отключить рекламу.