CC BY
36
Доклады БГУИР
2010 № 5 (51)
УДК004.056.55; 621.391.26
КЛАССИФИКАЦИЯ ТОЧЕЧНЫХ ОБРАЗОВ И КЛАССИЧЕСКАЯ ПРОБЛЕМА
РАЗБИЕНИЯ ЧИСЕЛ
В.А. ЛИПНИЦКИЙ*, В.К. КОНОПЕЛЬКО, Н.В. СПИЧЕКОВА
*Военная академия Республики Беларусь Минск-57,220057, Беларусь
Белорусский государственный университет информатики и радиоэлектроники П.Бровки,6 Минск, 220013, Беларусь
Поступила в редакцию 4 октября 2010
Изучается действие симметрической группы на строках и столбцах квадратных матриц порядка п, п > 2, с элементами 0 и 1, содержащих в точности п единиц. Устанавливается тесная связь проблемы классификации возникающих при этом действии орбит с проблемой разбиения чисел, классической в комбинаторике и теории чисел. С помощью формулы количества разбиений числа установлены границы для количества орбит. Сформулированы необходимые условия принадлежности матриц одной орбите.
Ключевые слова: распознавание образов, разбиение чисел, формула Эйлера, действие группы на множестве, матрица, симметрическая группа, орбита.
Постановка задачи
Данная работа является прямым продолжением исследований, начатых в [1]. Развивается общий подход к изучению действия группы = Sn х Sn - квадрата симметрической группы - на множестве Рп всех квадратных (0,1)— матриц порядка п , содержащих в точности п единиц. Неподвижные точки относительно действия отсутствуют - каждая орбита содержит более одного элемента. В то же время нет и полных — орбит - содержащих по |Б2| = п! 2 элементов. Исключение составляют лишь две полные орбиты - по одной из множеств Р3 и Р4 (см. [1]). Внутреннее строение - орбит характеризуется расположением единиц по строкам и столбцам их представителей. Описание же этого расположения впрямую взаимосвязано с классической проблемой разбиения чисел, хорошо известной в теории чисел и комбинаторике [2-4]. Данная статья посвящена установлению названной взаимосвязи.
Разбиения чисел
В дальнейшем будем использовать, взятое из монографии [2], следующее
Определение 1. Разбиением натурального числа n называется последовательность положительных чисел, сумма которых равна n. Слагаемые данной суммы называются частями разбиения. Разбиения, отличающиеся только порядком частей, считаются равными.
Важную роль в теории разбиений играет функция разбиений p(n). Она определяется как количество всех не равных друг другу разбиений числа n . Если n отрицательно, то p(n) =0. Последовательность, состоящая из нулей, считается разбиением нуля. Поэтому р(0) = 1. Значения р(п) растут очень быстро с ростом п . Так, р( 1) = 1, р(2) = 2,
/>(3) = 3, р(А) = 5, />(5) = 7, />(6) = 11, />( 10) = 42, />(20) = 627, />(50) = 204226, />(200) = 3972999029388 [2].
Имеет место следующая рекуррентная формула, называемая формулой Эйлера, которая позволяет последовательно находить значения функции p(n):
р(п) = р(п-1) + р(п - 2)-р(п-5) -р(п -7) +... + (-1)
9+1
ff Т 2 Л ( о 2 , М
Р
3q -q
п-
V V 2 ,
30'+« п--
2 /У
Если и невелико, то р(п) легко вычислить, просто выписав все разбиения. Вывод этой рекуррентной формулы можно найти в [2].
Разбиения числа и векторы распределений единиц
Всякой матрице А еРп соответствуют два вектора = (\|.л2.....\л ) и сА - (с1,с2,...,сп) распределений единиц по строкам и столбцам, здесь число 5 равно весу / - й строки данной матрицы, то есть количеству единиц в этой строке; число с ■ равно весу ] - ого столбца матри-
п п
цы А . Так как ^ Л'( = ^ с, = и, то таким образом каждой матрице А е Рп соответствует пара
(5А, сА) разбиений натурального числа п , если допускать в определении 1 и нулевые части..
При таком допущении, у матриц А и В, принадлежащих одной — орбите, пары разбиений (5А, сА ) и (^в, св ) одинаковы. Ведь у пар (5А, сА ) и (вв, св ) лишь компоненты !л и !в, а также компоненты сА и св могут отличаться друг от друга перестановкой координат.
По векторам !л и сА можно соответственно рассчитать числа 5ХгА и со1А, равные числу ненулевых координат в 5Л и сА .
Через JA обозначаем - орбиту, порождённую матрицей А е Рп .
Предложение 1. Среди матриц орбиты ^ существует такая матрица К, у которой пара ($к, ск ) распределения единиц по строкам и столбцам имеет обратный лексикографический порядок координат своих компонент: > я2 > ... > ; сх >с2 > ... > сп.
Доказательство. Пара (як, сг ) получается из пары (1А, сА ) подходящей перестановкой координат у каждого из векторов , и с ,. Пусть перестановка / е Л'л такая, что /'(л7,) = \ . а перестановка g&Sn такая, что g(cA) - ск . Переставим строки матрицы А в соответствии с требованиями подстановки /, а затем переставим у полученной матрицы столбцы в соответствии с требованиями подстановки g. В результате получим требуемую матрицу К Предложение 1 полностью доказано.
Замечание. Очевидно, для всех матриц данной орбиты JА пара (^, ск) с обратным лексикографическим распределением единиц определена однозначно, хотя соответствующих матриц К может быть и несколько. Чтобы подтвердить принадлежность пары (як, сг ) именно к орбите ^ будем впредь её обозначать через (5К, сК ) .
Пример 1. Орбиту 3Е, порожденную единичной матрицей Л е Рп . характеризует пара с) векторов —ск — (1, 1,..., 1). В качестве матрицы К здесь можно взять любую из п\ матриц -орбиты 3Е.
Спектр распределения единиц в — орбитах
Определение 2. Для каждой ^ - орбиты .! , пара распределений из пред-
ложения 1 называется спектром распределения единиц по строкам и столбцам в этой орбите. Из предложения 1 и замечания к нему вытекает
Предложение 2. Спектр , с К) является однозначной характеристикой каждой - орбиты 3А е Рп.
Возникает обратная задача: пусть (5, с) - (.V,, 5и; с,, с2,..., сп) - пара разбиений числа п , при этом > я2 > ... > 5И; с1 > с2 > ... > сп; существует ли матрица А е Рп со спектром распределения единиц (5^ = (-V, с) ?
Ответ почти очевиден - не всегда. Например, не существует матрицы в Рп при п > 2, у
которой бы все п единиц были бы расположены в первой строке и в первом столбце одновременно. Только такая матрица могла бы иметь в качестве спектра распределения единиц пару (5, с) = (п, 0,... ,0; п, 0,... ,0) разбиений числа п .
Определение 3. Пара (I, с) = (51з 52с1з с2,...,сп), ^ > > ... >
сх > с2 > ... > сп, разбиений числа п называется реализуемой, если найдётся матрица А е Рп со спектром распределения единиц (5^ ,с£) = (I, с).
Весовое содержание единиц и весовой портрет ф, 1 — матрицы
Введем еще две характеристики матриц многообразия Ри .
Определение 4. Пусть одна из единиц матрицы А е Рп расположена на пересечении 1 — й строки и / — го столбца этой матриц. Пусть 5г. — вес /' - й строки матрицы А, а с ■ - вес 7 - го столбца этой матриц. Пара чисел (я, с7 ) называется весовым содержанием
названной единицы. Совокупность пар (я, с7 ) всех единиц матрицы А называется ее весовым портретом.
Непосредственно доказывается следующее
Предложение 3. Если матрица А е Рп содержит единицу с весовым содержанием (я, с), то каждая матрица орбиты содержит единицу с таким же весовым содержанием. У всех матриц — орбиты .1, одинаковый весовой портрет. Предложение 3 позволяет ввести следующее
Определение 5. Весовым портретом - орбиты ./, называется весовой портрет любой матрицы этой орбиты.
Необходимые признаки принадлежности матриц одной — орбите
Приведём небольшой и далеко не полный список общих свойств матриц, принадлежащих одной — орбите. Такие свойства являются необходимыми признаками принадлежности
двух матриц одной орбите. Если для одной из них хотя бы один из необходимых признаков выполняется, а для другой - не выполняется, то эти матрицы должны принадлежать различным £2 — орбитам.
Признак 1. Если матрицы А и В принадлежат одной - орбите, то векторы 5А и (векторы сА и св) отличаются друг от друга лишь перестановкой координат.
Признак 2. Если матрицы А и В принадлежат одной Бп - орбите, то соответствующие им параметры 81гЛ и ъ1гв (со1А и со1в) равны между собой.
Признак 3. Если матрицы А и В принадлежат одной - орбите, то их весовые портреты одинаковы.
Множество Р является частью множества всевозможных квадратных матриц порядка п с коэффициентами 0 и 1. Эти коэффициенты можно интерпретировать как элементы поля Галуа из двух элементов ОЕ{2) = 0,1 . Все поля Галуа данного порядка изоморфны друг другу, в частности, ОЕ(2) изоморфно полю классов вычетов 2 / 22 . Этот изоморфизм и определяет однозначно арифметику поля ОЕ(2) . Таким образом, Р можно считать частью линейного пространства Мп (ОЕ(2)) всех квадратных матриц порядка п с коэффициентами из поля ОЕ(2) . Тогда строки любой матрицы А е 1>п образуют систему векторов пространства строк с п координатами над полем ОЕ(2) . Ранг этой системы векторов в линейной алгебре носит
название ранга матрицы А .
Из свойств ранга матрицы непосредственно вытекает
Признак 4. Если матрицы А и В принадлежат одной Б2п - орбите, то ранги этих матриц равны между собой.
Данное условие позволяет ввести следующее
Определение 6. Рангом Б2 - орбиты <А> называется ранг любой из матриц этой орбиты и обозначается через г(< А >).
Введенные характеристики - спектр распределения единиц, ранг, весовой портрет -удобно использовать для распознавания Б2 - орбит.
Пример 2. Векторы распределений единиц по строкам и столбцам матриц
А =
'111100^ "011110"
100000 100000
000010 000000 и В = 100000 000000
000000 000000
чоооооо; чооооооу
одинаковы: яА =яв = 4, 1, 1, 0, 0, 0 ; сА=св = 2,1,1,1,1,0
Но галц(А) = Зф гап^(Н) = 2, поэтому матрицы А и В принадлежат различным Л',^ - орбитам, орбиты < А > и < В > различны.
Пример 3. Возьмем следующие три матрицы из множества Р6:
'111000^ '111000^ '101100"
110000 100100 110000
000100 000000 ; В = 010000 000000 ; С = 010000 000000
000000 000000 000000
чооооооу чооооооу чооооооу
А =
У всех трех матриц ранг равен трем, одинаковы у них и спектры распределения единиц по строкам и столбцам: - 1В = 1С = (3,2,1,0,0,0) ; сА-св-сс- (2,2,1,1,0,0).
В матрице А имеется единица с весовым содержанием (1, 1), у матриц В и С такой единицы нет. Следовательно, — орбита ./ , Ф-,//( и .] , Ф .]г . Матрица В содержит единицу с весовым содержанием (2, 1), у матрицы С такой единицы нет. Следовательно, ./,., Ф .1 г . Все три орбиты , , попарно различны.
Нижняя и верхняя оценки количества Sn — орбит
Теорема 1 .Для каждого значения /7, п > 2, мощность |РИ множества Рп / S^ всех S2 — орбит находится в границах
I ч! 1
p(n)< \Pn / S2„\<-C"n2.
I I n
Доказательство. Выше отмечено, что не каждая пара разбиений числа n реализуется в качестве спектра распределения единиц той или иной матрицы множества Ри . Однако для любого разбиения числа n , взятого в качестве вектора s , найдётся такое иное разбиение числа n , взяв которое в качестве вектора c , мы получим реализуемую пару (s, c) . Отсюда сразу же
следует неравенство: p(n)< pn /S^ |.
Докажем, что каждая Л'^ — орбита содержит не менее п матриц; тогда будет установлена вышеуказанная верхняя граница количества Л'^ — орбит. В любой Л'^ — орбите найдётся
матрица А с ненулевой первой строкой. Если в матрице А все п единиц расположены в первой строке, то орбита < А > имеет мощность п (см. в примере 1 [1]). Пусть в первой строке матрицы А расположено т единиц, где 1 < г < [///2] ([£] — целая часть числа к). В орбите < А > найдется матрица G, в которой первая строка содержит т единиц, расположенных в первых г столбцах. Тогда в силу предложения 2.5 [5] S2n —орбита должна содержать не менее п матриц.
Пусть в первой строке матрицы А расположено // единиц, где п > jU >[п/ 2]. Тогда в матрице А найдётся другая строка, содержащая 1 <т <п — /и единиц, следовательно т <[п/ 2]. Повторяя проведенные рассуждения, убеждаемся, что и в этом случае |< А >| > п . Теорема 1 полностью доказана.
Заключение
Проблема классификации точечных образов сведена к описанию орбит на множестве всех квадратных (0,1) — матриц с п единицами. Орбиты возникают под действием симметрической группы на строках и столбцах названных матриц. Установлена тесная взаимосвязь распределения единиц по строкам и столбцам с проблемой разбиения чисел. Известная рекуррентная формула Эйлера числа разбиений натуральных чисел послужила основой для установления нижней и верхней оценок количества исследуемых орбит. Введенные характеристики орбит -спектр распределения единиц, ранг, весовой портрет - удобно использовать для распознавания орбит.
POINT PATTERN CLASSIFICATION AND CLASSICAL NUMBER PARTITION
PROBLEM
V.A. LIPNITSKI, V.K. KONOPELKO, N.V. SPICHEKOVA
Abstract
The action of the symmetric group on the rows and columns of square matrices of order n, n > 2, with elements 0 and 1, containing exactly n units is studied. A close relationship between the problem of classification of orbits arising in this action with the number partition problem, classical in the theory of combinations and number theory is established. Limits for the number of orbits are found
by means of the formula for the number of partitions. The necessary conditions for matrices belonging to one orbit are formulated.
Литература
1. Конопелько В.К., Липницкий В.А., Смолякова О.Г., Спичекова Н.В. // Докл. БГУИР. 2010. № 5. С. 40-47.
2. Эндрюс Г. Теория разбиений. М., 1982.
3. Холл М. Комбинаторика. М.,1970.
4. Hardy G.H., Wright E.M. An Introduction to the Theory of Numbers. Berlin, 1973.
5. Липницкий В.А., Конопелько В.К. Норменное декодирование помехоустойчивых кодов и алгебраические уравнения. Минск, 2007.
