Труды Карельского научного центра РАН №7. 2020. С. 110-119 DOI: 10.17076/mat1171
УДК 519.115:519.2
КОМБИНАТОРНЫЙ АНАЛИЗ ПОДСТАНОВОК С ФИКСИРОВАННЫМ ЧИСЛОМ ЦИКЛОВ
Н. Ю. Энатская
Московский институт электроники и математики,
Национальный исследовательский университет «Высшая школа экономики», Россия
Изучается схема n-размерных подстановок с к циклами по направлениям перечислительной комбинаторики. Это перечисление исходов схемы с определенной дисциплиной их нумерации, нахождение их числа, установление взаимнооднозначного соответствия видов исходов с их номерами, называемое задачей нумерации в прямой и обратной постановках, и моделирование исходов схемы.
Ключевые слова: подстановки; порционные добавления; циклы.
N. Yu. Enatskaya. COMBINATORIAL ANALYSIS OF THE PERMUTATIONS WITH A FIXED NUMBER OF CYCLES
The scheme of n-dimensional permutations with к cycles is studied in different areas of enumerative combinatorics. This includes direct enumeration of the outcomes of the scheme with a given discipline of their numbering, finding their numbers, establishing the one-to-one correspondence of the types of outcomes with their numbers, termed the direct and inverse numbering problem, and modelling of the outcomes of the scheme.
Keywords: the permutations; portion; the cycles.
Введение
Подстановка заданного размера задается перестановкой элементов того же размера своей нижней строки, пространство элементарных исходов для которой состоит из всех равновероятных перестановок ее элементов. Все исходы нашей (исследуемой) схемы совпадают с исходами схемы перестановок и составляют ее часть, поэтому вероятность появления исходов нашей схемы среди всех исходов схемы подстановок определяется отношением соответствующих численностей их исходов. Пространство элементарных исходов нашей схемы составляют совокупности составов ее непустых циклов в заданном числе. Вероятностное распределение ее исходов будем определять.
Число подстановок размера п с к циклами известно, например, из [2] и выражается через числа Стирлинга 1-го рода Б(п, к) как \Б(п,к)\. Здесь рассматриваются вопросы прямого перечисления подстановок заданного размера с фиксированным числом циклов, определение их числа, перечисление размеров циклов в них, нахождение вероятностного распределения исходов схемы и чисел циклов в случайных подстановках размера п, моделирование возможных исходов схемы с фиксированным числом циклов.
Определения и обозначения в схеме подстановок
Подстановка размера п задает последовательность взаимно-однозначных отображений всех п элементов множества {1,2,...,п} на се-
®
бя и записывается в виде двустрочного соответствия элементов отображения (12 ... п) и (¿1 ¿2 ... ¿п) в матричной форме.
Здесь в верхней строке элементы упорядочены от 1 до п, а нижняя строка задает последовательные отображения элементов верхней строки в элементы, стоящие под ними в нижней строке, т. е. отображения: 1 в ¿1 (1 ^ ¿1), 2 в ¿2 и т. д., где {¿1 ,¿2, ...,гп} е {1, 2,..., п}, ¿к = ¿т при к = т. Таким образом, подстановки можно задавать перестановками элементов ее нижней строки.
Цикл подстановки образует группа элементов, возвращающихся при данном отображении к исходному элементу.
Число участвующих в цикле элементов называется его размером, принимающим значения от 1 до п. (Сумма размеров всех циклов подстановки равна ее размеру п.)
Если размер цикла п, т. е. он единственный, то подстановку называют одноцикловой.
Цикловая структура подстановки задается вектором а = (од, а2,..., ап), где а - число циклов размера ¿, % = 1, п; ХП=1 а% = к - число циклов подстановки, а^п=1 = п.
Цикловой группой будем называть номера элементов одного цикла в порядке их отображений, начиная с наименьшего номера элемента в цикле. Иногда бывает удобно представлять подстановку не в традиционной двустрочной форме, а в виде перечисления ее цикловых групп в порядке роста минимальных номеров в них и их размеров, а каждую — в порядке отображений элементов в ней. Например, при п = 4, к =2 такая запись подстановки (1,3),(2,4) будет соответствовать перестановке ее нижней строки (3,4,1,2).
В [5] были решены задачи нахождения чис-ленностей подстановок по их заданной цикловой структуре, явного перечисления таких подстановок и моделирования их возможных значений.
Этапы перечисления подстановок фиксированного размера с заданным числом циклов
Отображение элементов подстановки задается перестановкой ее нижней строки (i1,i2,...,in) по ПРАВИЛУ: каждый элемент ¿1 перестановки в подстановке отображается в свой элемент перестановки с номером ¿1. Отсюда получаем АЛГОРИТМ** вычисления числа циклов подстановки по ее нижней строке-перестановке:
по перестановке выписываем ее элементы начиная с 1 и первого элемента перестановки (слева направо) и далее в порядке отображений до повт,ора 1 без ее включения - получаем
первый цикл, далее начинаем второй цикл с минимального элемента, не вошедшего в первый цикл (слева направо), и далее в порядке отображений до его повтора без включения -получаем второй цикл, и т. д. до последнего оставшегося элемента, а для подсчета числа циклов при каждом повторе элемента в цикле добавляем в счетчик единицу.
Приведем поясняющий пример. Пусть п = 8 и перестановка нижней строки есть последовательность: (6,7,5,1,2,8,3,4). Тогда по АЛГОРИТМУ** и по ПРАВИЛУ из нее получаем подстановку с двумя циклами (1,6,8,4)(2,7,3,5).
Комбинаторный анализ схемы построим на процедуре явного перечисления подстановок размера п с фиксированным числом циклов к, представленных своими цикловыми группами, соответствующее перечисление которых и будем здесь производить. Эта процедура будет состоять из двух этапов: деления всех элементов перестановки всеми способами на к непустых частей (циклов) без учета их порядка и перечисления для каждого результата такого деления всех допустимых исходов перестановок элементов внутри каждой части деления, приводящих к одноцикловой подстановке из элементов этой части.
Имея в виду математическую эквивалентность нашего деления совокупности элементов на части с размещением п различимых частиц по к неразличимым ячейкам без пустых ячеек, при рассмотрении решения задачи первого этапа будем ссылаться на соответствующие результаты в терминах размещения частиц по ячейкам. Поэтому решение задачи первого этапа может использовать результаты [3], описывающие перечисление всех исходов схемы (без ограничения на отсутствие пустых ячеек), с отбраковкой исходов, не отвечающих данному ограничению. Это привело бы к анализу большого числа лишних исходов и к сложности дальнейших исследований схемы. В связи с этим предлагается провести прямое перечисление исходов с данным ограничением.
В [8] приведена процедура перечисления всех исходов схемы перестановок в определенном порядке с решением задачи нумерации для них, и известно утверждение, что число одноцикловых подстановок размера Ь, а значит, и представляющих их перестановок нижних строк, равно (Ь — 1)!.
Поэтому для решения задачи второго этапа кроме отбраковки исходов перестановок элементов каждой части деления, приводящих к > 1 циклу в них, с использованием утверждения о числе одноцикловых подстановок фик-
сированного размера можно предложить прямое перечисление требуемых перестановок ее элементов по [3], начинающихся с минимального.
В завершение с каждым результатом первого этапа по конкретной цикловой структуре соответствующей подстановки, определяемой делением ее элементов на части (цикловые группы), нужно сопоставить все исходы второго этапа в виде требуемых перестановок элементов по каждой части деления, проводимого по результатам [4] и [6].
Построение процедуры перечисления первого этапа требует отдельного рассмотрения.
Определения, обозначения и связи основной и вспомогательных схем
Для краткости под различимыми элементами или частицами будем понимать их номера.
Схема А — основная схема подстановок размера п с ровно к циклами; вспомогательные схемы: схема Б — схема деления п различимых элементов на к непустых частей без учета порядка этих частей; схема Б. — схема Б с ]-й фиксацией к минимальных элементов в частях деления, где в качестве этих минимальных элементов варьируются все их наборы по схеме сочетаний; схема Р — схема всех разных допустимых перестановок элементов в частях деления по схеме Б, записывающих одноцик-ловые подстановки из них; схема С — схема сочетаний.
Схемы Б и Р являются вспомогательными для схемы А, схемы {Б.} — для схемы Б, а схема С — для схем {Б.}.
Схема Р является схемой одновременных действий схем перестановок и подробно исследована в [6] и [8], а схема С — в [7]. Для анализа схемы А остается изучить схемы V и Б..
Анализ схем будем производить на основе прямого перечисления ее исходов в терминах указанных выше размещений частиц по ячейкам. Для перечисления всех исходов схемы Б будем объединять все исходы всех схем Б., т. е. переборов всех размещений частиц с изученной в схеме сочетаний [7] с фиксацией их минимальных номеров в ячейках с последующим размещением остальных частиц при каждой такой фиксации. Начнем с рассмотрения вспомогательной схемы Б..
3,...,Скп-\)1 добавляя к нему частицу с номером 1, и принудительно разложим по всем к ячейкам по одной частице, считая их номера минимальными при размещении по ячейкам остальных частиц и расположив ячейки в порядке роста минимальных номеров выбранных частиц. Для таких выборов по схеме сочетаний в [4] решены задачи построения процедуры перечисления исходов и задача их нумерации. Далее размещаем по ячейкам всеми способами остальные (п — к) частиц таким образом, чтобы не изменить выбранные минимальные номера частиц в ячейках. Состав каждой ячейки будем заключать в круглые скобки и перечислять номера входящих в них элементов в возрастающем порядке (а весь исход схемы будем заключать в фигурные скобки).
Пусть Т(а) = (п\,...,пк) - возрастающая последовательность минимальных номеров частиц в к ячейках при их ]-й фиксации, где всегда п\ = 1. Это соответствует несущественности порядка частей (циклов) деления элементов подстановки. Будем размещать остальные (п — к) частиц по к ячейкам так, чтобы выбранные минимальные в них не изменились. Для этого будем добавлять частицы порционно по одной в последовательно поеди-нично растущие группы первых ячеек равновероятно по схеме размещения с повторением для каждой частицы по этим ячейкам. Тогда сначала в первую ячейку должны быть добавлены все частицы с номерами < п2 в количе-И)
стве = п2 — п\ — 1 = п2 — 2; потом в первые две ячейки - все частицы из остальных с но-
(а)
мерами < п3 в количестве —2 = п3 — п2 — 1; потом в первые три ячейки - все частицы из остальных с номерами < п4 в количестве
ш
а)
3
1. Комбинаторный анализ схемы Б. дает с (п — к):
= п4 — п3 — 1 и т. д., потом в первые (к — 1) ячейки - все частицы из остальных с номерами (а)
< пк в количестве \ = пк — пк-\ — 1; тогда, наконец, во все ячейки - все остальные части-
(з)
цы с номерами до п в количестве —к = п — пк.
Таким образом, т. к. к частиц было размещено предварительно, будем порционно в указанных выше количествах размещать суммарное число (п — к) частиц соответственно по одной, двум и т. д. к ячейкам. Проверяем, что сумма порционных добавлений частиц совпа-
1.1. Перечисление и число исходов схемы
Для выполнения условия отсутствия пустых ячеек случайным образом по схеме С сочетаний с СПк—\ исходами возьмем ^-й исход выбора нумерованных частиц (] = 2,
п2 — 1 — 1+пз —п2 — 1+п4 — пз — 1 + . . .+пк — пк-\
— 1+ п — пк = п — к.
Будем представлять все порционные добавления в виде к-компонентных векторов соответ-
112
ствующих добавлений частиц в к ячеек, а итоговые векторы добавлений во все ячейки при ]-й фиксации частиц с минимальными номерами в них получаются покомпонентным сложением векторов порционных добавлений в ячейки (циклы).
Теорема 1. Число N0. исходов схемы Б? вычисляется по формуле
МВз =П
„ш
(1)
i=1
,,(j)
Доказательство формулы (1) следует из принципа умножения для общего числа исходов в к схемах размещения частиц (по поеди-нично растущим от начальной ячейки группам ячеек) в количествах найденных выше размеров, заданных вектором —Я, последовательных порционных добавлений в них по схеме размещений с повторением, так что в 1-ю
(?)
ячейку добавляется —1 частиц числом спосо-
Л- 1т(3) (?)
бов 1Ш1 ; в первые две - —2 частиц числом способов 2ад23) и т. д.; во все ячейки - ча-
стиц числом способов к'Шк , откуда и следует формула (1).
Определим вид исхода схемы Б? при фиксированном наборе Т(?) минимальных номеров частиц в ячейках как перечень наборов возрастающих номеров частиц в ячейках, перечисленных в порядке роста минимальных номеров в них; составы ячеек будем заключать в круглые скобки, а весь исход — в фигурные скобки.
Для перечисления всех исходов схемы методом графов будем строить процесс последовательного поединичного равновероятного добавления частиц с растущими номерами в ячейки с начальным заполнением по одной частице с минимальными выбранными номерами в них. Исходы заполнения ячеек на каждом шаге будем нумеровать в порядке добавления частицы в ячейки с растущими минимальными номерами.
Поясним процесс перечисления исходов схемы Б? на примере (рис. 1).
Пример 1. Пусть п = 7, к = 3, Т(?) = (1,3, 5).
Визуально имеем N0- = 18 исходов, что совпадает с результатом по (1): N0. =1 • 21 • 32 = 18.
Рис. 1. Граф перечисления исходов схемы Dj примера 1
Fig. 1. The enumeration graph of outcomes of scheme Dj of example 1
1.2. Задача нумерации (ЗН) в схеме Б?
Пространство элементарных исходов схемы Б? состоит из всех к наборов составов из п различимых элементов с заданной ]-й фиксацией минимальных элементов этих составов.
Схема Б? представляет собой схему (п — к) последовательных действий, изученную в [6] (с решенной задачей нумерации при известных числах их исходов), состоящих в порционном размещении по одной (п — к) частиц (элементов) с размерами порций ,■?,...,■?, определенными в п. 1.1, по поединично растущим группам первых ячеек, расположенных по возрастанию минимальных номеров частиц в них так, что каждая частица размещается по ячейкам своей группы равновероятно по схеме размещений с повторением. Тогда число способов совершения этих
(з) (з) (з)
действий есть 1™1 2™2 • • • к'Шк , где числа исходов этих действий V = (у1,у2,...,уп-к) есть (1, 2,...,к), повторенные соответственно
1Ш1 , 2Ш2 ,...,кШк раз, а итоговым исходом схемы Б? является объединение исходов действий, т. е. суммарных порционных составов номеров добавленных в ячейки частиц, и j-х фиксированных минимальных номеров частиц в к ячейках.
Для иллюстрации приведем числовой пример.
Пример 2. Пусть в условиях примера 1 для схемы Б? решается ЗН, т. е. пусть п =
7, к = 3, Т(?) = (1, 3, 5). Тогда = 1,—2? =
(?)
1,—з = 2, откуда следует, что при числе действий (п — к) = 4 V = (1,2,3,3). Введем обозначения: N , Я* - соответствующие номер и вид итогового исхода схемы Б?; N(г), К(г) -
113
соответствующие номер и вид исхода схемы Dj на i-й итерации - добавления i-й частицы (i = 1,n — k).
Прямая ЗН. Пусть дан N* = 5. Требуется найти R*.
Ш!аги решения по формулам (2) и (3) из [6] в обозначениях нашей схемы:
1) N* = N(4) = 5; N(3) = [(5 + 3 — 1)/3] = 2; N(2) = [(2 + 3 — 1)/3] = 1; N(1) = [(2 + 2 — 1)/2] = 1 - номера исходов схемы по действиям;
2) j4 = N(4) mod v4 + 0 • v4 = 5 mod 3 +
0 • 3 = 2; аналогично j3 = 2 mod 3 + 0 • 3 = 2; j2 = 2 mod 2 + 1 • 2 = 2; ji = 1 mod 1 +
1 • 1 = 1 - номера предшествующих исходов в пучках предыдущих итераций, т. е. полученный набор значений {j(t)} = (1,1,2,2) задает номера частей деления (циклов) добавления растущих остальных (n — k) =4 элементов, откуда получаем искомый вид исхода R* = (124), (367)(5), совпадающий с результатом по графу на рис. 1.
Обратная ЗН. Пусть дан R* (124)(367)(5). Требуется найти N*.
Ш!аги решения по формуле (4) из [8] в обозначениях нашей схемы:
1) по R* = R(4) = (124)(367)(5) и T(j = (1, 3, 5) с учетом добавления остальных элементов в порядке роста имеем R(1 = (12), (3)(5), R(2) = (124)(3)(5), R(3) = (124)(36)(5), R(4) = (124)(367)(5), откуда соответственно получаем номера частей (циклов) добавления элементов в порядке их роста ji = 1, j2 = 1, j3 = 2, j4 = 2, которые являются номерами предшествующих искомому исходов итераций в пучках графа перечисления исходов схемы;
2) по (4) из [8] N* = N(4) = (1 — 1)1 + (1 — 1)2 + (2 — 1)3 + 2 = 5, что совпадает с результатом по графу на рис. 1.
2. Комбинаторный анализ схемы D
2.1. Перечисление и число исходов схемы
Теорема 2. Число Nd исходов схемы D есть
nd = Е ND
{(c )}
(2)
где сумма берется по перечислению исходов схемы С, т. е. по всем С,—1 наборам Т(а = (1,п2, ...,пк) е (1, 2,.. .,п).
Доказательство. Схема Б — это схема независимых ОД (см. [6]) схем {Б.}, ] = 1,2,..., С,¡к,
а ее исходы — объединение исходов всех схем Б. в порядке роста ], т. е. перечисление всех возможных составов к циклов исходной схемы А, откуда и следует (2).
2.2. Задача нумерации в схеме Б
Схема Б — это схема С,—\ независимых одновременных действий (см. [6]) итоговых порционных добавлений частиц в ячейки по схемам Б., для которых в [6] задача нумерации решена при числах способов совершения этих
действий, равных {N0.}, где ^ = 1,С,—\.
3. Комбинаторный анализ схемы А
Исходы схемы А (в фигурных скобках) представляют все возможные наборы (в круглых скобках) к цикловых групп подстановки размера п, где порядок номеров в них имеет значение. Они отличаются (в терминологии подстановок) от исходов схемы Б, где номера элементов цикловых групп (частиц в ячейках) перечислялись только в возрастающем порядке, всеми порядками их номеров, задающими возможные последовательные отображения внутри цикловых групп.
3.1. Перечисление и число исходов схемы
Для перечисления исходов схемы А по утверждению из введения нужно в перечислении исходов схемы Б по схеме одновременных действий (см. [6]) произвести в каждой круглой скобке, описывающей состав цикловой группы, все перестановки последовательных отображений, начинающиеся с минимального номера.
Число NA подстановок размера п с к циклами известно, и есть NA = |£(п, к)|, где Б(п,к) - числа Стирлинга 1-го рода. Покажем, как число NA можно вычислять по перечислению исходов схемы А; предложим для этого два способа.
Первый способ (по прямому перечислению исходов схемы Б)
Теорема 3. Число NA исходов схемы А вычисляется по формуле
к
Na = Е П dm
{(D)} m=1
(3)
где суммирование проводится по перечислению исходов схемы Б по размерам суммарных порционных добавлений элементов d = (й\,...,йк) в цикловые группы, где
(Ект=\ ^ = п — к).
114
Доказательство. Все исходы схемы А делятся на N0 групп одинаковых составов элементов в к ненулевых частях (циклах), каждая из которых представляет собой одноцик-ловые подстановки из элементов своего состава в известном количестве факториала суммарного размера своего порционного добавления. Таким образом, число N4 вычисляется как сумма N0 слагаемых (по перечислению разных по составу делений совокупности п различимых элементов на к ненулевых частей) произведений вариантов перестановок элементов каждой части этого деления, приводящих к одноцикловым подстановкам, т. е. произведений факториалов размеров суммарных порционных добавлений в них, что и доказывает формулу (3) с введением обозначений
й = й1,...,йк, (^'т=1 йт = п — к) для этих размеров. (Они непосредственно вычисляются из видов исходов схемы N0 при их перечислении.)
Второй способ (по перечислению исходов схемы С) алгебраическим методом
Пусть по перечислению исходов схемы С (см. [7]) задан j-й набор Т(?). Тогда при этой j-й фиксации на ¿-м шаге перечисления исходов схемы Б? (при размещении частиц по ¿ ячейкам) перечисляющей производящей функцией всех размеров порционных добавлений й1,...,йг в ячейки-части 1, 2,...^ есть полином
Рг? = (а1 + 0,2 +---+ аг)Гг
Е
-А
41+-----+?*=т,
•••йг!
а?*
(4)
в разложении по степеням переменных а1 ,...,аг, i = 1,2,...,к которого все наборы показателей их степеней соответственно совпадают с векторами размеров порционных добавлений в первые ¿ циклов общего суммарного размера добавлений гг = пг+1 — пг — 1, i = 1,2,...,к — 1, г к = п — пк номеров элементов, а коэффициенты при а?1 ... а?* дают кратность получения векторов размеров порционных добавлений. Для получения при данном j всех итоговых векторов из й нужно по всем i = 1, 2,.. .,к покомпонентно сложить все векторы размеров порционных добавлений во все циклы, перебирая их по схеме последовательных действий (см. [6]). Далее по всей совокупности этих итоговых векторов при всех j = 1,2,...,СП определяется число NA по формуле (3).
Перечислим порядок действий при вычислении числа N4.
Первый способ
1. Фиксируем минимальные номера элементов подстановки из п номеров в к циклах Т по схеме сочетаний С числом способов СП—{, т. к. всегда п1 = 1, и перечисляем все исходы этих фиксаций по [4].
2. Для каждой j-й фиксации (у = 1,2,..., СП—\) исхода схемы С перечисляем все исходы схемы Б? методом графов (см. рис. 1 для примера 1).
3. Объединяя результаты п. 2, получаем все исходы схемы Б, вычисляем {N0.}.
4. Для всех исходов п. 3 непосредственно по перечислению вычисляем все итоговые векторы размеров порционных добавлений в к циклах — получаем все {с?} по исходам схемы Б.
5. По (3) вычисляем значение N4.
Второй способ
1. Фиксируем минимальные номера элементов подстановки из п номеров в к циклах Т по схеме сочетаний С числом способов СП—\, т. к. всегда п1 = 1, и перечисляем все исходы этих фиксаций по [7].
2. Для каждой j-й фиксации (у = 1,2,..., С-1) исхода схемы С и всех i = 1,2,...,к рассматриваемых первых циклов вычисляем значения гг и выписываем вид (4) для Рг?.
3. Из вида Рг? находим все векторы размеров порционных добавлений в первые ¿ циклов.
4. Для получения всех итоговых векторов й при каждом j складываем покомпонентно (компоненты соответствуют номерам переменных в полиноме Рг?) все со всеми при разных ¿, полученных в п. 3).
5. й в схеме Б получаем объединением результатов п. 4) по всем j.
6. По й из п. 5) по (3) вычисляем N4.
Приведем числовой пример и будем проводить расчеты по приведенному порядку действий двух способов.
Пример 3. Пусть п = 4, к = 2. Требуется найти число подстановок длины п с ровно к циклами, т. е. N4.
1. СП—1 = С1 = 3 — число вариантов фиксаций минимальных номеров в двух циклах: п1 = 1,п2 = 2 при j = 1; п1 = 1,п2 = 3 при j = 2 и п1 = 1,п2 = 4 при j = 3.
г
1
Решаем первым способом.
2. По (1) ^ = 12—\—124—2 = 4; ND2 = 13—\24—3 = 2; N^3 = 14—\—\24—4 = 1. Приведем графы перечисления исходов схем {В.} (рис. 2).
При з = 3 Т = (1,4), Рз\ = а\ = а\, откуда имеем один вектор размеров порционного добавления (2);
Р32 = (а\ + а2)Г2 = (а\ + а2)0 = 1, откуда имеем один вектор размеров порционного добавления (0;0).
4. Вычисляем итоговые векторы размеров порционных добавлений:
Рис. 2. Граф перечисления исходов схем Dj примера 3
Fig. 2. The enumeration graph of outcomes of schemes Dj of example 3
3. По рис. 2 (по процедуре порционного добавления частиц в ячейки из п. 1.1) имеем итоговые векторы размеров добавления при j = 1 — (2,0), (1,1),(1,1),(0,2); при j = 2 - (2,0),(1,1); при j = 3 - (2,0).
4. По (3) NA = 2 !0 ! + 1 ! 1 ! + 1 ! 1 ! + 0 ! 2 ! + 2 ! 0 ! + 1 !1! + 2 !0 ! = 2+1 + 1 + 2 + 2+1 + 2 = 11.
Решаем вторым способом.
2 - 3. При j = 1 T(1) = (1, 2), Pii = a1 = a0, откуда имеем один вектор размеров порционного добавления (0);
P12 = (ai + a2)r2 = (ai + a2)2 = a1 + 2aia2 + a22, откуда имеем 4 вектора размеров порционного добавления (2;0),(1;1),(1;1),(0;2).
При j = 2 T^ = (1,3), P21 = a]1 = a1, откуда имеем один вектор размеров порционного добавления (1);
P22 = (a1 + a2)r2 = (a1 + a2)1 = a + a2, откуда имеем 2 вектора размеров порционного добавления (1;0),(0;1).
при з = 1 это векторы (2;0),(1;1),(1;1),(0;2);
при 3 = 2 это векторы (2;0),(1;1);
при 3 = 3 это вектор (2;0). 5. По (3) NA = 2!0! + 1!1! + 1!1! + 0!2! + 2!0! + 1111 + 2101 = 11.
Проверим результат непосредственным подсчетом чисел циклов в подстановках размера 4 по их явному перечислению по п. 1.1.
Представим полный перечень исходов схемы перестановок размера п = 4: (4321), (3421),(3241),(3214),(4231),(2431),(2341),(2314), (4213),(2413),(2143),(2134),(4312),(3412),(3142), (3124),(4132),(1432),(1342),(1324),(4123),(1423), (1243),(1234). Считая их нижними строками подстановок, перечислим соответствующие им непосредственно вычисленные числа циклов: 2,1,2,3,3,2,1,2,2,1,2,3,1,2,1,2,2,3,2,3,1,2,3,4. Отсюда видим, что число подстановок размера 4 с двумя циклами равно 11, что совпадает с результатом теоретического расчета.
Замечание 1. Вероятностное распределение числа Z циклов в подстановке размера п определяется по формуле Р^ = к) = \8(п,к)\/п\, где к = 1,2,...,п.
Выпишем по примеру 2 для п = 4 вероятностное распределение числа Z циклов, определяя вероятности каждого возможного размера циклов (от 1 до 4) делением их числа на 4! = 24. Р^ = 1) = 6/24, Р^ = 2) = 11/24, Р^ = 3) = 6/24, Р^ = 1) = 1/24.
3.2. Задача нумерации в схеме А
Схема А — это обобщенная схема двух последовательных действий (см. [6]) с числом исходов первого действия, равным N3, и числами исходов второго действия, равными числу всех перестановок во всех цикловых группах с первыми местами элементов с минимальными номерами, т. е. Пкгп=\ dш\ по каждому исходу схемы П. Таким образом, задача нумерации в схеме А теоретически принципиально решена на основании ее результатов для всех используемых здесь вспомогательных схем, исследованных в [4] и [6], и указанных связей этих схем со схемой А.
Для пояснения порядка действий при решении задачи нумерации в схеме А введем обозначения: и я(и) — соответственно номер и вид исхода с этим номером в схеме и (под и понимается любая из используемых здесь схем), ^ — знак следования. Тогда в этих обозначениях приведем решения прямой и обратной задач нумерации через известные аналогичные результаты задач нумерации всех используемых схем при заранее найденных размерах пучков графа перечисления исходов схемы А.
116!
Прямая задача: по N А найти НА).
Решение: по [4] и [6] N(А) ^ N(0) ^ Я(0); N(А) ^ N(р) Я(р); заменяя порядки номеров в цикловых группах в Я(° на полученные в Н(р), находим На.
Обратная задача: по НА найти N А).
Решение: упорядочиванием по возрастанию номеров каждой цикловой группы в Я(А получаем Я(А) ^ Я(0) ^ N(0); из Я(А) ^ N (Р). Теперь по [4] из
Однако в связи с двухэтапностью и с много-действенностью процедуры прямого перечисления перестановок заданного размера с заданным числом циклов явный результат задачи нумерации, следующей из этого перечисления, становится достаточно громоздким даже с учетом полученных ранее результатов решения для всех ее этапов. Поэтому здесь еще предлагается решать задачу нумерации численно, табличным путем, т. е. путем явного перечисления ее исходов с их табличной нумерацией всех исходов в порядке их получения при перечислении для заданных параметров схемы (размере подстановки и числе ее циклов) и хранить в памяти эту таблицу для дальнейших исследований.
3.3. Вероятностные распределения исходов схемы и числа исходов с к циклами в подстановке степени п
Все исходы схемы А представляют собой часть п! равновероятных исходов всех п-мерных подстановок, поэтому вероятность исходов схемы А среди всех исходов подстановки равна ^/п!.
С использованием числа N4 = (п, к) \ или формулы (3) получаем вероятностное распределение числа 2 циклов в п-размерной подстановке W = Р(2 = к) = N4/п!, к = 1~й.
Вместо исходов схемы п-мерных подстановок будем рассматривать задающую ее схему перестановок ее нижней строки, изученную в [8]. Будем удалять недопустимые состояния полного графа перечисления исходов этой схемы перестановок в нашей схеме - схеме А, которыми будем считать состояния всех траекторий, приводящих к подстановке с отличными от заданного числа к циклами. Вероятностное распределение исходов схемы А будем находить пропорциональным пересчетом вероятностей переходов по итерациям в каждом пучке среди оставшихся после описанного удаления траекторий, т. е. урезания недопустимых состояний, что означает пропорциональное перераспределение вероятностей итерационных
переходов среди всех остальных состояний в пучке, содержащем это состояние.
Распределение вероятностей исходов схемы А может быть получено алгоритмически в виде групп равновероятных исходов по каждой траектории одинаковых наборов вероятностей переходов по итерациям следующими шагами:
1) по методу графов (см. [9]) перечисляем все перестановки нижних строк подстановок;
2) для всех результатов п. 1) по АЛГОРИТМУ** (см. Введение) вычисляем числа циклов в соответствующих подстановках;
3) по результату п. 2) в графе п. 1) убираем все лишние траектории, ведущие к перестановкам, порождающим подстановки с числами циклов = к с пересчетом по итерациям равных вероятностей переходов в каждом оставшемся пучке среди оставшихся состояний;
4) вычисляем вероятности всех итоговых исходов с к циклами по графу п. 3) перемножением вероятностей переходов в них по итерациям, т. е. по их траекториям.
Теперь приведем пояснительный пример вычисления вероятностного распределения исходов схемы.
Пример 4. Пусть п = 4, к = 2.
1) Запишем все исходы схемы перестановок в порядке их получения по МГ по итерациям через стрелочку, начиная с i = 1:
(1) ^ (2,1), (1, 2) ^ (3, 2,1)), (2, 3,1), (2,1, 3), (3,1, 2), (1, 3, 2), (1, 2, 3) ^ (4, 3, 2,1), (3, 4, 2,1), (3, 2, 4,1), (3, 2,1, 4), (4, 2, 3,1), (2, 4, 3,1), (2, 3, 4,1), (2, 3,1, 4), (4, 2,1, 3), (2, 4,1, 3), (2,1, 4, 3), (2,1, 3, 4), (4, 3,1, 2), (3, 4,1, 2), (3,1, 4, 2), (3,1, 2, 4), (4,1, 3, 2), (1, 4, 3, 2), (1, 3, 4, 2), (1, 3, 2, 4), (4,1, 2, 3), (1, 4, 2, 3), (1, 2, 4, 3), (1, 2, 3, 4);
2) числа циклов в подстановках с нижними строками результатов п.1) соответственно: 2,1,2,2,3,2,1,2,2,1,2,3,1,2,1,2,2,3,2,3,1,2,3,4, из них 12 исходов с к = 2;
3) убираем все лишние траектории в графе п.1), получаем граф:
(1) ^ (2,1), (1, 2) ^ (3, 2,1)), (2, 3,1), (2,1, 3), (3,1, 2), (1, 3, 2), (1, 2, 3) ^ (4, 3, 2,1), (3, 2, 4,1), (3, 2,1, 4), (2, 4, 3,1), (2, 3,1, 4), (4, 2,1, 3), (2,1, 4, 3), (3, 4,1, 2), (3,1, 2, 4), (4,1, 3, 2), (1, 3, 4, 2), (1, 4, 2, 3). Здесь все исходы 3-й итерации, кроме первого (3,2,1) и последнего (1,2,3), порождают по 2 итоговых исхода нашей схемы с вероятностями переходов в них по 1/2, первый исход - 3 с вероятностью 1/3, а последний - 1 с вероятностью перехода, равной 1.
Отсюда из того, что вероятности переходов до 3-й итерации (включительно) не изменились и равны по итерациям, начиная с первой
P(4, 3, 2,1) = (1/2)(1/3)(1/3) P(2, 3,1, 4) = P(3, 4,1, 2) = P(1, 3, 4, 2) =
1: (1/2), (1/3), находим вероятности всех 12 исходов:
P(3, 2, 4,1) = P(3, 2,1, 4) = = 1/18; P(2, 4, 3,1) = P(4, 2,1, 3) = P(2,1, 4, 3) = P(3,1, 2, 4) = P(4,1, 3, 2) = (1/2)(1/3)(1/2) = 1/12; P(1, 4, 2, 3) = (1/2)(1/3) = 1/6.
Проверка на распределение: 3/18 + 8/12 + 1/6 = 1.
3.4. Моделирование подстановок размера n с ровно k циклами
Рассмотрим АЛГОРИТМ моделирования исходов нашей схемы через исходы общей схемы перестановки по следующим шагам:
1) перечисляем с нумерацией все n! равновероятных исходов схемы перестановок размера по алгоритму n, совпадающему с алгоритмом Кнута, приведенным, например, в [1] с доказанной там сложностью O(n!);
2) вычисляем числа циклов подстановок, заданных всеми перестановками своих нижних строк, и фиксируем среди них номера исходов нашей схемы (с заданным числом циклов) и их общее число L и находим p = L/n! (АЛГОРИТМ** вычисления числа циклов в подстановке, приведенный во Введении, имеет сложность O(n)) - получаем вектор Ck номеров исходов перечисленных нижних строк подстановок с числом циклов k;
3) разыгрываем номер исхода схемы п. 1);
4) проверяем его совпадение с номерами Ck исходов нашей схемы (п. 2)) на «ДА» и «НЕТ» с вероятностями соответственно W и (1 — W);
5) при «НЕТ» переходим к п. 3);
6) при «ДА» с вероятностью W по результату прямой ЗН для нашей схемы получаем смоделированный исход схемы.
Реализация последовательных действий 1-го и 2-го шагов составляет сложность O(n n!); реализация последовательных действий 3-го и 4-го шагов составляет сложность O(n!); число моделирований исходов схемы подстановок до получения M исходов нашей схемы имеет отрицательное биномиальное распределение с параметрами W, M как числа Z испытаний Бернулли до M-го успеха включительно (успех - получение исхода нашей схемы среди всех исходов схемы подстановок) с вероятностью каждого успеха W; P(Z = t) = CtM_-1Wt(1 — W)t-M и средним, равным M/W. Поэтому оценка сложности в среднем для моделирования M исходов нашей схемы получается порядка O(n n!M/W), откуда для моделирования одного исхода она имеет порядок O(n n!/W).
3.5. Приближенное оценивание числа исходов схемы
Моделируем по 1.2.3 или [3] NA исходов схемы перестановок размера п с п! равновероятными исходами. Считая каждую смоделированную перестановку нижней строкой подстановки того же размера, а исходом моделирования - такое получение подстановки, для каждого исхода вычисляем число циклов и определяем среди смоделированных количество М исходов с к циклами. Тогда искомое число исходов нашей схемы подстановок с к циклами приближенно определяется методом пропорций по формуле
N
Mn!
= N .
Исследуем качество полученной оценки N для N исходов нашей схемы, где M/Na - наблюденная частота успеха опыта - появления ее исхода среди Na исходов схемы перестановок. Число M можно представить в виде M = Xx + X2 + ... + Xna , где при i = 1,Na {Xi} -случайные величины, имеющие распределение Бернулли с вероятностью успеха p = N/n!. Тогда по уточненной по неравенству Чебышева теореме Бернулли выполняются соотношения
Y = P
MN
Na
n!
< е
P
Mn!
Na
N
< e*n!e
= P(|NN — N| < е) > 1 — P(1 — P) > 1 — (n!)2
NA(e*)2
4NAe2'
где для оценки N числа N исходов нашей схемы £ - ее точность, а 7 ^ 1 — 4мл£2 - оценка ее надежности 7 с этой точностью. Для получения нетривиальной оценки для 7 потребуем, чтобы NA > (п!)2/4£2.
Более точную оценку надежности 7 с заданной точностью £ оценки N числа N исходов нашей схемы можно получить по следствию из теоремы Муавра - Лапласа из соотношений
Y = P
MN
Na
n!
^ е*
P
Mn!
Na
N
< e*n!e
P(|N — N\ < е) и 2Ф (eVNa/p(1 — p)) ^ 2Ф(2еyNA/n!),
118,
где Ф(х) = (1^л/2П) /0 ехр 02/2 йх - функция Лапласа.
Литература
1. Кнут Д. Искусство программирования на ЭВМ. Т. 1. М.: Мир, 1976. 736 с.
2. Риордан Дж. Введение в комбинаторный анализ. М.: Изд-во иностр. лит., 1963. 287 с.
3. Энатская Н. Ю. Комбинаторное представление схемы размещения различимых частиц по неразличимым ячейкам // Дискретная математика. 2017. Т. 29, вып. 1. С. 126-135. ао1: 10.4213^т1410
4. Энатская Н. Ю. Комбинаторный анализ обобщенной схемы последовательных действий // Промышленные АСУ и контроллеры. 2016. № 4. С. 25-27.
5. Энатская Н. Ю. Комбинаторный анализ случайных подстановок заданных цикловых
структур // Промышленные АСУ и контроллеры. 2017. № 11. С. 29-34.
6. Энатская Н. Ю. Комбинаторный анализ схем одновременных и последовательных действий // Промышленные АСУ и контроллеры. 2016. № 2. С. 35-41.
7. Энатская Н. Ю. Комбинаторный анализ схемы сочетаний // Промышленные АСУ и контроллеры. 2015. № 8. С. 33-38.
8. Энатская Н. Ю, Колчин А. В. Комбинаторный анализ схемы перестановок // Труды КарНЦ РАН. 2014. № 4. С. 80-86.
9. Энатская Н. Ю, Хакимуллин Е. Р. Метод графов для решения задач перечислительной комбинаторики // Приборы и системы. Управление, контроль, диагностика. 2014. № 8. С. 15-21.
Поступила в редакцию 10.01.2020
References
1. Knuth D. Iskusstvo programmirovanya na EVM [The art of computer programming]. Moscow: Mir, 1976.
2. Riordan D. Vvedenie v kombinatornyi analiz [Introduction to combinatorial analysis]. Moscow: Izd-vo inostr. lit., 1963. 287 p.
3. Enatskaya N. Yu. Kombinatornoe predstav-lenie skhemy razmeshcheniya razlichimykh chastits po nerazlichimym yacheikam [Combinatorial representations for the scheme of allocations of distinguishable particles into indistinguishable cells]. Diskretnaya matem. [Discrete Math. Appl.]. 2017. Vol. 29, no. 1. P. 126-135. doi: 10.4213/dm1410
4. Enatskaya N. Yu. Kombinatornyi analiz obobshchennoii skhemy posledovatel'nykh deistvii [Combinatorial analysis of the generalized scheme of the sequential actions]. Promyshlennye ASU i kontrollery [Industrial ACS and Controllers]. 2016. No. 4. P. 25-27.
5. Enatskaya N. Yu. Kombinatornyi analiz sluchainykh podstanovok zadannykh tsiklovykh struktur [Combinatorial analysis of random permutations of given structures of cycle].
Promyshlennye ASU i kontrollery [Industrial ACS and Controllers]. 2017. No. 11. P. 29-34.
6. Enatskaya N. Yu. Kombinatornyi analiz skhem odnovremennykh i posledovatel'nykh deistvii [Combinatorial analysis of the schemes of simultaneous and sequential actions]. Promyshlennye ASU i kontrollery [Industrial ACS and Controllers]. 2016. No. 2. P. 35-41.
7. Enatskaya N. Yu. Kombinatornyi analiz skhemy sochetanii [Combinatorial analysis of combination scheme]. Promyshlennye ASU i kontrollery [Industrial ACS and Controllers]. 2015. No. 8. P. 33-38.
8. Enatskaya N. Yu., Kolchin A. V. Kombina-tornyi analiz skhemy perestanovok [Combinatorial analysis of a permutation scheme]. Trudy KarNTs RAN [Trans. KarRC RAS]. 2014. No. 4. P. 80-86.
9. Enatskaya N. Yu., Khakimullin E. R. Metod grafov dlya resheniya zadach perechislitel'noi kombinatoriki [Graphs method for solving enumerative combinatorics]. Pribory i sistemy. Upravlenie, kontrol', diagnostika [Instruments and Systems. Management, Monitoring, and Diagnostics]. 2014. No. 8. P. 15-21.
Received January 10, 2020
СВЕДЕНИЯ ОБ АВТОРБ:
Энатская Наталия Юрьевна
доцент Департамента прикладной математики, к. ф.-м. н.
Национальный исследовательский университет
«Высшая школа экономики», Московский институт
электроники и математики
ул. Таллинская, 34, Москва, Россия, 123458
эл. почта: nat1943@mail.ru
тел.: +79037411345
CONTRIBUTOR:
Enatskaya, Natalia
National Research University Higher School of Economics,
Moscow Institute of Electronics and Mathematics 34 Tallinskaya St., 123458 Moscow, Russia e-mail: nat1943@mail.ru tel.: +79037411345