ПРИЛОЖЕНИЕ Сентябрь 2014
Секция 1
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ПРИКЛАДНОЙ ДИСКРЕТНОЙ МАТЕМАТИКИ
7
УДК 519.7
О МИНИМАЛЬНЫХ ПРИМИТИВНЫХ МАТРИЦАХ
Р. И. Бар-Гнар, В. М. Фомичев
Исследуется подход к минимизации реализации преобразований, степень которых реализует полное перемешивание входных данных. Введены понятия минимальной примитивной матрицы и высоты примитивной матрицы. Получены оценки числа минимальных примитивных матриц порядка п. Построены и проанализированы алгоритмы поиска минимальных примитивных матриц и оценки близости примитивной матрицы к множеству минимальных примитивных матриц.
Ключевые слова: примитивная матрица, решётка, антицепь, вычислительная сложность алгоритма.
Введение
Для оценки свойства перемешивания входов композиции преобразований часто применяется матрично-графовый подход [1, гл. 10], основанный на определении экспонентов перемешивающих матриц (графов). Обзор основных результатов по этому направлению имеется в [2].
Аппаратная и программная реализация совершенных преобразований, реализующих полное перемешивание, затруднена необходимостью реализации большого количества связей между компонентами входных и выходных векторов. В связи с этим в работе введено и исследовано понятие минимальной примитивной матрицы, позволяющей оценить предельные возможности уменьшения количества связей между входом и выходом при сохранении свойства полного перемешивания у композиции преобразований.
1. Свойства минимальных примитивных матриц
Обозначим через М0(п) множество квадратных 0,1-матриц порядка п, Б(п) —множество подстановочных матриц, Р(п) —множество примитивных матриц.
Примитивная матрица (граф) называется минимальной, если после замены любой единицы нулём (после удаления любой дуги) получается непримитивная матрица (граф). Обозначим множество минимальных примитивных матриц через Ртш(п).
Утверждение 1. Если примитивные матрицы А, В € М0(п) сопряжены в Б(п) и примитивны, то они одновременно минимальные или неминимальные.
Рассмотрим на множестве М0(п) отношение частичного порядка: А ^ В ^ а^- ^ Ь^ для всех г,] = 1,... , п, где А = (а^),В = (Ь^). Множество М0(п) образует решётку в смысле отношения ^, где матрица, все элементы которой равны 1, является наибольшим элементом, а нулевая матрица — наименьшим. Подмножество примитивных
матриц Р(п) множества М0(п) является верхней полурешёткой. Минимальные примитивные матрицы — это все минимальные элементы верхней полурешётки Р(п).
В частично упорядоченном множестве X антицепь У назовем наибольшей, если она имеет наибольшую мощность, и максимальной, если любой элемент множества X сравним с некоторым элементом антицепи У.
Утверждение 2. Рт1п (п) есть максимальная антицепь полурешётки Р(п).
Число единиц матрицы А Е М0(п) называется весом матрицы А (обозначается ||А||). Назовём к-м слоем подрешётки К решётки М0(п) подмножество (обозначается К(к)), состоящее из матриц веса к, где к = 0,1,... , п2. Заметим, что К(к) есть антицепь решётки К при всех к.
Утверждение 3. При п > 3 выполнена оценка
3п! ^ |РШ1п(п)| ^ С^72-1, где т = п2.
Нижняя оценка получена с использованием свойства Р(га+1)(п) С Рт1п(п), где
Р (п+1)
(п) содержит три класса орграфов, порядки которых равны п!:
1) полный цикл с добавленной дугой (г — 1, г + 1), г Е {1,..., п} (рис. 1,а);
2) полный цикл с добавленной петлей в вершине г Е {1,..., п} (рис. 1,б);
3) класс графов, вид которых зависит от чётности числа п: при чётном п граф состоит из двух циклов длины 2 и / = п — 1 (рис. 1,в); при нечётном — из двух циклов длины 2 и / = п (рис. 1,г). В обоих случаях (2,/) = 1, поэтому в соответствии с критерием [3, с. 226] графы примитивные.
в г
Рис. 1. Примеры минимальных графов из Р(га+1)(п)
Улучшение нижней оценки возможно на основе обобщения случая 3 (рассмотрения некоторых классов графов из слоев Р(га+г)(п) при г > 1).
2. Высота примитивной матрицы
Расстоянием по Хэммингу между 0,1-матрицами А и В (обозначается ^А,В)) называется число различных соответствующих элементов матриц А и В:
а(А,В) = £(а„- 0 Ъг,у).
г,3
а
Высотой примитивной матрицы A (обозначается h(A)) называется расстояние по Хэммингу между A и ближайшей минимальной примитивной матрицей M Е Pmin(n):
h(a) = min d(A, M).
M €Pmin(n)
Величина h(A) является определённой мерой избыточности при построении связей между элементами входа и выхода преобразований информации.
Алгоритм оценивания h(A) (метод координатного спуска):
1) последовательно просматривая элементы матрицы A, находим единицы;
2) заменяем найденный единичный элемент в матрице A нулевым и для полученной матрицы A' выполняем:
— проверяем в A' наличие нулевых строк и столбцов; если таковые есть, возвращаемся к выполнению п. 2 для следующей единицы матрицы A;
— матрицу A' без нулевых строк и столбцов проверяем на примитивность;
— если матрица A; не примитивная, возвращаемся к матрице A, восстанавливаем замененную единицу и выполняем п. 2 для следующей единицы матрицы A;
— примитивную матрицу A; проверяем на минимальность;
— если A' минимальная, то определяем m — число единиц, заменённых нулями;
— если A' не минимальная, то переходим к выполнению п. 1 для матрицы A'.
На выходе алгоритма получим минимальную примитивную матрицу M, где
h(A) ^ d(A, M) = m.
Сложность алгоритма полиномиальная. Пусть ||A|| = k, где n < k ^ n2 для примитивной матрицы A. Тогда вычислительная сложность алгоритма в битовых операциях в худшем случае не превышает величины 0(k2n3 log n). Объём памяти требуется порядка O(n2) битов. Данный метод можно использовать для поиска минимальных матриц, близких к перемешивающим матрицам раундовых подстановок блочных шифров (DES, ГОСТ 28147-89 и др.).
ЛИТЕРАТУРА
1. Фомичев B. M. Методы дискретной математики в криптологии. М.: Диалог-МИФИ, 2010.
2. Когос К. Г., Фомичев В. М. Положительные свойства неотрицательных матриц // Прикладная дискретная математика. 2012. №4(18). С. 5-13.
3. Сачков B. Н., Тараканов B. Е. Комбинаторика неотрицательных матриц. М.: ТВП, 2000.
УДК 519.1:511.2
ЧИСЛА ЭЙЛЕРА НА МНОЖЕСТВАХ ПЕРЕСТАНОВОК И АНАЛОГИ
ТЕОРЕМЫ ВИЛЬСОНА
Л. Н. Бондаренко, М. Л. Шарапова
Определяются числа Эйлера на множествах перестановок и с их помощью доказываются аналоги теоремы Вильсона для чисел стандартных полных отображений и чисел стандартных сильных полных отображений.
Ключевые слова: перестановка, числа Эйлера, полные отображения, теорема Вильсона.