Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2022. Т. 26, № 3. С. 573-591 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1947
EDN: AWBYGO
УДК 519.115
Вероятностные модели для анализа обратных экстремальных задач комбинаторики
Н. Ю. Энатская
Национальный исследовательский университет «Высшая школа экономики», Московский институт электроники и математики им. А. Н. Тихонова, Россия, 123458, Москва, ул. Таллинская, 34.
Аннотация
В обратной экстремальной задаче для комбинаторной схемы при заданном значении целевой функции вида определенного экстремального значения ее характеристики строится вероятностная модель, обеспечивающая получения этого значения в ее исходах. Рассматривается два типа таких характеристик, относящихся к каждому или совокупности исходов схемы.
Доасимптотический анализ такой модели проводится авторским перечислительным методом. Его основу составляет построение итерационного случайного процесса с итерациями последовательных этапов нумерованного бесповторного перечисления и формирования исходов схемы. Итерационное развитие процесса представляется вероятностным графом.
Исследование исходов схемы по модели в перечислительном методе проводится по следующим направлениям: визуального нумерованного представления исходов схемы, нахождения их числа, установления взаимно-однозначного соответствия между видами и номерами исходов схемы, получения их (управляемого случайным процессом перечисления исходов схемы) вероятностного распределения и их моделирования с этим распределением.
Наряду с непосредственным исследованием схем по указанным направлениям предлагаются алгоритмы получения результатов для них путем их частичного пересчета из результатов аналогичного анализа более общих, ранее изученных схем с меньшими ограничениями на значения рассматриваемых характеристик.
Ключевые слова: обратная экстремальная задача, экстремальное значение характеристики, доасимптотический анализ схемы.
Математическое моделирование, численные методы и комплексы программ Научная статья
© Коллектив авторов, 2022 © СамГТУ, 2022 (составление, дизайн, макет)
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Энатская Н. Ю. Вероятностные модели для анализа обратных экстремальных задач комбинаторики // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2022. Т. 26, № 3. С. 573-591. EDN: AWBYGO. DOI: 10.14498/vsgtu1947. Сведения об авторе
Наталия Юрьевна Энатская А https://orcid.org/0000-0003-1241-7543 кандидат физико-математических наук, доцент; доцент; деп. прикладной математики; e-mail: nat1943@mail.ru
Получение: 12 августа 2022 г. / Исправление: 25 августа 2022 г. / Принятие: 31 августа 2022 г. / Публикация онлайн: 20 сентября 2022 г.
1. Введение. В литературе известны экстремальные задачи комбинаторики, состоящие в отыскании среди конечного множества альтернатив одной, отвечающей заданному значению целевой функции1 (см., например, [1,2]). В работах [3-5] проводился асимптотический анализ вероятностного поведения максимальных уровней заполнения ячеек в схемах равновероятного размещения различимых частиц по различимым ячейкам с разными условиями.
Здесь (при конечных значениях параметров) изучается процедура построения модели случайного процесса перечисления исходов комбинаторных схем, приводящих к заранее заданным экстремальным (минимальным или максимальным) значениям их определенных характеристик. Такие задачи можно рассматривать как обратные экстремальные задачи комбинаторики с заданными значениями этих характеристик, являющихся целевыми функциями.
2. Постановка задачи. В доасимптотических условиях значений параметров комбинаторной схемы при заданном экстремальном значении его определенной характеристики нужно построить вероятностную модель бесповторного нумерованного перечисления ее исходов и провести их анализ по всем направлениям ПМ.
Основные определения, обозначения и понятия:
Е — изучаемая схема с заданным экстремальным значением ее характеристики;
5 — начальная схема для схемы Е, соответствующая более общей схеме (схему 5 при пересчете из нее результатов для схемы Е считаем изученной);
итерация — шаг перехода в графе перечисления исходов схемы к следующему этапу (шагу) их перечисления;
ПП — процедура прямого перечисления всех исходов схемы с учетом ее специфики;
ПМ — перечислительный метод исследования моделей по приведенным в аннотации направлениям; он состоит в построении вероятностных моделей на формирующих исходы итерационных случайных процессах по-единичного добавления этапов перечисления исходов ранее изученных схем с управляемыми (через вероятности итерационных переходов) вероятностными распределениями их исходов;
ЗН — задача нумерации по установлению взаимно-однозначного соответствия между всеми номерами и видами исходов схемы;
ПЗН —прямая задача нумерации, состоящая в нахождении вида исхода схемы по его номеру;
ОЗН —обратная задача нумерации, состоящая в нахождении номера исхода схемы по его виду;
траектория исхода т — последовательность исходов итераций, ведущая к нему в графе от начала их перечисления;
хКлассический пример — известная задача коммивояжера: нахождения кратчайшего кольцевого маршрута возвращения в начальный город с посещением заданных городов по разу.
пучок в графе процесса — совокупность переходов из каждого исхода каждой итерации;
размер пучка на каждой итерации — число переходов из каждого исхода итерации;
пучковая структура графа — перечисление размеров пучков всех итераций;
УМ —универсальный (единый) метод моделирования исходов схемы, состоящий в разыгрывании одним случайным числом его номера с сформированным вероятностным распределением исходов схемы и получении его смоделированного вида по результату решения ПЗН.
3. Вспомогательные результаты. Приведем краткие достаточные для дальнейшего чтения сведения основных используемых далее опубликованных результатов.
3.1. Перечислительный метод (ПМ). ПМ разработан для проведения доасимптотического анализа комбинаторных схем, т.е. при конечных значениях их параметров [6]. Его суть состоит в организации получения качественной информации об исходах схемы и переводе ее в количественную — результатов ее анализа. Эта качественная информация представляет собой исходы итерационного случайного процесса реализации комбинаторной схемы в результате последовательного поединичного (итерационного) добавления элементов схемы до заданного параметром значения или этапов перечисления составляющих схему более простых ранее изученных схем. Инструментами такого перевода качественной информации о видах всех исходов схемы являются метод графов (МГ), задача нумерации (ЗН), задание итерационных переходных вероятностей (управляющих распределением итоговых исходов схемы) и универсальное моделирование (УМ) исходов с этим их вероятностным распределением. Одним из важных приемов ПМ для получения новых результатов доасимптотического анализа комбинаторных схем является операция по перечислению (ОПП). Ее введение дает возможность вычисления характеристик схем, зависящих от видов предсостояний исходов при их формировании итерационным случайным процессом как функций от них. Целью применения ПМ является изучение схем Е по указанным направлениям.
3.2. Обобщенная схема последовательных действий (ОПД) Под действиями схемы понимаются реализации исходов итераций при последовательном добавлении этапов перечисления всех промежуточных и итоговых исходов схемы.
Схема ОПД с приводимыми здесь результатами анализа из [7] возникает, когда каждому следующему действию (итерации) подвергаются исходы предыдущего действия и числа исходов каждого следующем действия могут быть неодинаковыми, т.е. зависят не только от характера действия, но и от вида предыдущего исхода. Результатом этого являются возможно разные размеры пучков внутри итераций в графе перечисления исходов схемы при переходе от исходов предыдущего действия к последующему.
При анализе схемы ОПД считаем известными результаты всех направлений исследования ПМ составляющих ее комбинаторных схем действий.
В схеме проводится к последовательных действий, г-тое из которых совершается N « способами, % = 1 ,к. Число пучков каждой итерации равно числу исходов предыдущей итерации. Тогда число исходов этих к действий складывается из N(й-1) пучков размерами (¿(г) = ,й^,..., ), т.е. общее
число N = Nисходов схемы получается из рекуррентного соотношения при г = к и N = N(0) = 1
N
N« = £ 1=\
Вид исхода после совершения г действий будет формироваться из принятых видов исходов последовательных действий с обозначениями: г — номер действия, а ]г — номер исхода в результате его совершения.
Задача нумерации решается для нашей схемы при решенной ЗН для каждого из к действий при известной пучковой структуре графа перечисления исходов нашей схемы.
Прямая и обратная задачи нумерации решены следующими теоремами
(к)
Теорема 1. Совершается к действий и задан номер исхода Щ . Тогда его вид, определяемый номерами исходов траектории т в содержащих их пучках 1^} от первой до к-той итераций, вычисляется по рекуррентной формуле для
N ГХ)-1
зг = тг) - £ г = 1Л
1=1
где все пучковые структуры действий й(г) заданы и
м(к-1) = 5 + тахI : ^ ^ = Ак < м(к)^ ,
где 5 = 0 при = и 5 = 1 при < .
По решенной ЗН для всех действий находим по {¿г} виды их исходов, из которых получаем искомый вид исхода К* .
Теорема 2. Совершается к действий и задан вид исхода = {]1,... , }. Тогда его номер Ы*^ определяется по рекуррентной формуле при г = к, г = 1к,
N ГХ)-1
^ = £ + зг, 1=1
начиная с г = 1 при Ж*(1) = з1.
4. Два типа характеристик в схемах Е. Будем различать два типа характеристик Т — типа А и типа В:
1) характеристика типа А для исхода схемы — определяемая и сравниваемая внутри исхода схемы и принимающая в нем ^ 1 значений;
2) характеристика типа В для совокупности исходов схемы — однозначно определяемая каждым исходом схемы и используемая для ее сравнения в совокупности исходов.
Процедура прямого перечисления всех исходов схемы в этих случаях проводится по-разному.
В первом случае перечисляются все исходы, каждый из которых включает в себя часть, фиксирующую заданное экстремальное значение характеристики с добавлением всех вариантов остальной его части. Тогда все перечисляемые таким образом исходы схемы будут иметь заданное экстремальное значение характеристики.
Во втором случае перечисляются все исходы схемы со значениями характеристики Т до ее заданного экстремального значения t включительно, т.е. при задании ее максимального значения t — все исходы со значениями Т ^ t, а при задании ее минимального значения t — все исходы со значениями Т ^ t. Тогда данная характеристика совокупности всех перечисленных исходов схемы будет иметь заданное экстремальное значение.
5. Примеры процедур перечисления исходов схем и определения их численностей в схемах Е. Приведем примеры описания конкретных схем Е с разными указанными типами экстремальных характеристик их исходов.
Схемы типа А с фиксацией экстремального значения характеристики исходов первого типа — типа А:
АЛ. Схема перестановок с фиксированной максимальной (минимальной) длиной серии поединично растущих (убывающих) элементов. А.2. Схема перестановок с фиксированной максимальной (минимальной) длиной серии растущих (убывающих) элементов. A.3. Схема перестановок с фиксированной максимальной (минимальной) разностью между ее соседними элементами. АА. Схема подстановок с фиксированным максимальным (минимальным) размером цикла.2
A.5. Схемы размещения частиц по ячейкам с фиксированным значением максимального (минимального) уровня заполнения ячеек.
Схемы типа В с фиксацией экстремального значения характеристики исходов второго типа — типа В:
B.1. Схема перестановок с фиксированным максимальным (минимальным)
» 3
числом инверсий.3
В.2. Схема подстановок с фиксированным максимальным (минимальным) числом циклов.
B.3. Схема перестановок с фиксированным максимальным (минимальным)
элементом на фиксированном месте. В.4. Схема перестановок с фиксированной максимальной (минимальной) длиной серии поединично растущих (убывающих) элементов начиная с данного.
В.5. Схема перестановок с фиксированной максимальной длиной серии без
инверсий начиная с данного элемента. В.6. Схема сочетаний с фиксированным максимальным (минимальным) раз-
4
махом.4
2Цикл подстановки — перестановка подмножества элементов подстановки в порядке их последовательных отображений с замыканием, т.е. последний элемент перестановки отображается в первый. Размер цикла — число входящих в него элементов.
3Число инверсий элемента перестановки равно числу больших ее элементов, стоящих правее, а число инверсий перестановки есть сумма инверсий всех ее элементов.
В.7. Схема сочетаний с фиксированным максимальным (минимальным) элементом.
В.8. Схема сочетаний с фиксированным максимальным (минимальным) числом элементов из данного диапазона значений элементов исхода схемы.
Подобных примеров схем Е можно привести много, а анализ каждого требует отдельного рассмотрения, начальным этапом которого является построение процедуры прямого перебора исходов с учетом природы каждой схемы. Здесь для примера приведем идеи прямого перечисления исходов комбинаторных схем Е в некоторых наиболее простых случаях (в основном в схемах типа В). Все другие схемы Е, требующие более сложных ПП исходов схем с более глубоким погружением в их специфику, оставим для отдельных исследований.
Примеры получения формул для чисел N исходов схем из их ПП (с указанием типа схемы в скобках). Для получения аналитических результатов исследования схемы Е требуется построение процедуры бесповторного ПП ее исходов методом графов. Она отражает закономерности связей между номерами и видами исходов, выявление которых составляет результаты решения ЗН, необходимые для дальнейшего исследования схемы Е по определенным здесь направлениям. Поэтому бесповторное ПП исходов схемы дает основу для проведения аналитического изучения схемы и является первой целью рассмотрения схем Е. Однако это потребует отдельного анализа каждой конкретной схемы Е с учетом ее специфики.
Все N исходов схем типа А представляют совокупность исходов с достижением характеристикой заданного экстремального значения £ в каждом исходе.
Все N исходов схем типа В состоят из двух совокупностей — первая из исходов с достижением характеристикой заданного экстремального значения £ и вторая — на не противоречащих этому экстремальному значению Ь без его достижения.
Введем обозначения для параметров схем в следующих ниже примерах, называя нумерованные элементы их номерами:
- в схеме сочетаний производится выбор из п различимых элементов по г без возвращения и без учета их порядка, тогда число исходов схемы равно Сгп;
- в схеме перестановок производится выбор из п различимых элементов по п без возвращения и с учетом их порядка, тогда число исходов схемы равно п!.5
Пример 1 (Тип В). Схема сочетаний с параметрами п, г и с фиксированным минимальным элементом Т = Ь ^ (п — г + 1).
ПП: Для перечисления первой совокупности исходов включаем в выборку элемент £ с добавлением выборок по схеме сочетаний остальных (г — 1) элементов из (п — ¿) элементов от (£ + 1) до п числом способов Сгп—\, а для второй совокупности исходов набираем всю выборку из элементов > £ числом способов откуда N = Сгп—\ + .
4Размах в исходе схемы сочетаний — разность между его максимальным и минимальным элементами.
5Перечисления исходов схем сочетаний, перестановок и ОПД см. в [7].
Пример 2 (Тип В). Схема сочетаний с параметрами п, г и с фиксированным максимальным элементом Т = £ ^ г.
ПП: Для перечисления первой совокупности включаем в выборку элемент Ь с добавлением выборок по схеме сочетаний остальных (г — 1) элементов из (£ — 1) элемента от 1 до (£ — 1) числом способов С^-1, а для второй совокупности исходов набираем всю выборку из элементов < £ числом способов 1, откуда N = С1-1 + С[_1.
Пример 3 (Тип В). Схема сочетаний с параметрами п, г и с фиксированным минимальным числом Т = Ь элементов от 1 до в, где тах(0, в — п + г) ^ ^ £ ^ тт^, г) = к.
ПП: Включаем в совокупность исходов схемы все выборки по схеме сочетаний из элементов от 1 до в по ¿, а остальные (г — ¿) элементов — по схеме сочетаний из элементов от (в + 1) до п числом способов С*Далее, варьируя значение г (число вошедших в выборку элементов от 1 до в) от £ до к, получаем перечисление всей искомой совокупности с общим числом исходов
N = тИ=ь с*
Пример 4 (Тип В). Схема сочетаний с параметрами п, г и с фиксированным максимальным числом Т = Ь элементов от 1 до в, где к = тах(0, 8—п+г) ^ ^ £ ^ тт(з, г).
ПП: Включаем в совокупность исходов схемы все выборки по схеме сочетаний из элементов от 1 до 8 по ¿, а остальные (г — ¿) элементов — по схеме сочетаний из элементов от (в + 1) до п числом способов ^ = С*. Далее, варьируя значение г (число вошедших в выборку элементов от 1 до в) от к до ¿, получаем перечисление всей искомой совокупности с общим числом исходов N = С]
Пример 5 (Тип А). Размещение г неразличимых частиц по п различимым ячейкам с фиксированным минимальным уровнем Т = ¿, г ^ Ш заполнения ячеек.
ПП: По схеме сочетаний Сгп способами выбираем г ячеек из п. По этим г ячейкам размещаем (г — пЬ) частиц без пустых ячеек по схеме сочетаний с повторением числом способов С^Г1 с последующим добавлением остальных пЬ частиц по Ь частиц в каждую из п ячеек при возможных значениях г в диапазоне к = тах(0, п — г + пЪ) ^ г ^ тт(п — 1, г — Ы) = К. Варьируя все эти значения г из [к, К] и суммируя соответствующие им числа исходов, получаем число исходов схемы N = ^=к СгпС^Г.1 с достижением значения Т = £ в каждом исходе.
Пример 6 (Тип В). Схема перестановок элементов с номерами от 1 до п с фиксированной минимальной длиной Т = Ь серии, начиная с 1, подряд идущих элементов.
ПП: Переставляем данную серию длиной £ элементов как один номер с остальными (п — ¿) элементами по схеме перестановок числом способов N = (п — Ь + 1)! и перечисляем их по схеме перестановок [4]. Тогда все перечисленные исходы схемы будут содержать серию заданного вида или длины > ¿, если сразу после нее стоит элемент (£ + 1), или содержать серию длины, равной ¿, в противном случае, т.е. минимальная длина серии в перечисленных исходах схемы будет равна ¿.
Пример 7 (Тип В). Схема перестановок элементов с номерами от 1 до п с фиксированной максимальной длиной Т = Ь серии, начиная с 1, подряд идущих номеров.
ПП: В совокупности исходов нашей схемы длины г серий данного вида ^ ¿. Переставляем данную серию из г элементов как один элемент с остальными (п — г — 1) элементами без элемента (г + 1) по схеме перестановок числом способов (п — г)! и перечисляем их по схеме перестановок. Далее вставляем элемент (г + 1) в перестановке до, между и после всех элементов полученных перестановок, кроме места сразу после последнего элемента серии. Число таких мест для него есть (п — г), откуда получаем число исходов (п — г)!(п — г). Теперь, варьируя длину г серии от 1 до заданного Ь и объединяя все такие исходы, получаем все исходы нашей схемы численностью N = ^¿=1^ — г)!(п — г). Тогда максимальная длина серии полученной совокупности исходов окажется равной ¿.
6. Решение ЗН в схемах Е. Решение ЗН в схемах Е проводится как непосредственным выявлением закономерностей связи при бесповторном ПП между видами и номерами ее исходов для каждой конкретной схемы отдельно с учетом ее специфики, так и пересчетом ранее полученных для нее результатов решения ЗН для начальной схемы 5. Представим эти два подхода к решению ЗН на конкретных примерах.
6.1. Первый подход к решению ЗН. Рассмотрим по одному примеру решения ЗН для схем Е с заданными экстремальными значениями характеристик двух введенных типов через результаты решения ЗН составляющих их схем.
Пример 8 (Тип А). Дана схема подстановок размера п = 5, максимальный размер цикла Т = 2. Провести прямое перечисление исходов схемы и решить ЗН.
Построим процедуру перечисления исходов схемы типа А. Это значит, что исходы схемы будут включать в себя все подстановки с размерами циклов ^ 2 с достижением размера 2 в каждом исходе.
Исходы схемы различаются наборами допустимых размеров циклов, их составами и порядками отображений в них.
Процедура прямого перечисления всех исходов схемы проводится тремя итерациями (этапами):
1) набором составов допустимых размеров циклов;
2) перечислением составов элементов в этих циклах;
3) перебором порядков отображений в циклах.
Очевидно, что здесь возможны два варианта наборов первой итерации размеров циклов (2, 2,1) и (2,1,1,1), в каждом из которых порядок перечисления составов циклов совпадающих размеров циклов несущественен, а при Т = 2 порядок отображений на 3-й итерации у них единственный.6 Поэтому процедура бесповторного ПП исходов схемы при каждом наборе их размеров сводится к делению всеми способами всех 5-ти элементов между циклами по этим наборам в заранее определенном порядке. Будем их перечислять, например, в порядке убывания их размеров, а среди циклов совпадающих размеров
6 Число перестановок в порядке отображений из т элементов есть (т — 1)! с фиксированным начальным элементом, при гп = 2 имеем (т — 1)! = 1! = 1.
(для бесповторности перечисления составов циклов) — в порядке роста минимальных элементов в них начиная с минимального в каждом составе цикла. Поэтому состав каждого цикла должен начинаться с минимального элемента подстановки из не вошедших в предыдущие циклы. Представим графом на рис. 1 процесс формирования и перечисления исходов схемы начиная с первой итерации при начальных исходах схемы совокупностей составов исходов циклов (по схемам сочетаний) с их размерами (2, 2,1) и (2,1,1,1) с подряд идущей общей нумерацией. Численности исходов схемы с каждым составом размеров циклов вычисляются по схеме перестановки с повторением (представляемой через схемы сочетаний [7]) без учета порядка совпадающих по размеру циклов. Они вычисляются делением числа исходов каждой схемы перестановок с повторением [7] на произведение факториалов чисел повторов размеров циклов. Это и будет учитывать неразличимость их порядков. Тогда получаем соответственно числа исходов из данных начальных составов размеров циклов = 5!/(2! ■ 2! ■ 1! ■ 2!) = 15, И2 = 5!/(2! ■ 1! ■ 1! ■ 1! ■ 3!) = 10 с общим числом исходов N = Ы1 + Ы2 = 25. Здесь Ы1 = 15 и Ы2 = 10 — размеры пучков предытоговой итерации, объединяющих наборы всех циклов в один этап перечисления исходов по каждому набору их размеров.
Считая решенной ЗН для этапов, составляющих нашу схему, и имея полный перечень всех исходов схемы и их пучковую структуру, будем решать ЗН, сверяя результаты с рис. 1.
Прямая ЗН. Номер исхода схемы ЫЕ = 18. Требуется найти его вид В*Е. По рис. 1 В*Е = (1, 4)(2)(3)(5). Найдем его из условий схемы ОПД по теореме 1:
1) на первой итерации найдем номер к1 пучка, содержащего искомый исход схемы к1 = 2, т.к. 15 < 18;
2) найдем номер к2 искомого исхода во втором пучке к2 = 18 — 15 = 3;
3) из результата решения прямой ЗН для второго этапа перечисления исходов схемы с составом размеров циклов (2,1,1,1) по номеру исхода к2 =3 находим Д(2) = В*Е = (1, 4)(2)(3)(5) (в соответствии с принятым порядком нумерации исходов при их перечислении), что совпадает с результатом по рис. 1.
Обратная ЗН. Вид исхода схемы В*Е = (1, 4)(2)(3)(5). Требуется найти его номер ЫЕ. По рис. 1 ЫЕ = 18. Найдем его из условий схемы ОПД по теореме 2:
1) на первой итерации по составу размеров циклов найдем номер к1 пучка, содержащего искомый исход схемы, по размерам его пучков — у нас к1 = 2;
2) по обратной ЗН во втором пучке для второй составляющей схемы с составом размеров циклов (2,1,1,1) находим к2 его номер в пучке к2 = 3;
3) искомый номер КЕ по теореме 2 находим как сумму размеров предшествующего пучка и к2, т.е. = 15 + 3 = 18, что совпадает с результатом по рис. 1.
Пример 9 (Тип В). Дана схема сочетаний из п = 9 по г = 3 с максимальным размахом V = 4. Исходы схемы сочетаний перечисляются в возрастаю-
(1,2)
(1,3)
(1,4)
(1,5)-
(2.3)-
(2.4)-
(2.5)-
3.4
3.5 4,5
2.4
2.5 4,5 2,3 2,5 3,5
2.3 2,5
3.4
4.5 3,5 3,4
(5)
(4)
(3)
(5)
(4) (2)
(5)
(3) (2)
(4) (3) (2) (1) (1) (1)
(1,2) — —(з) — —(4) — —(5) 16
(1,3) — —(2) — —(4) — —(5) 17
(1,4) — —(2) — —(3) — —(5) 18
(1,5) — —(2) — —(з) — —(4) 19
(2,3) — —(1) — —(4) — —(5) 20
(2,4) — —(1) — —(3) — —(5) 21
(2,5) — —(1) — —(з) — —(4) 22
(3,4) — —(1) — —(2) — —(5) 23
(3,5) — —(1) — —(2) — —(4) 24
(4,5) — —(1) — —(2) — —(3) 25
numbers of outcomes 1 2
3
4
5
6
7
8
9
10 11 12
13
14
15
Рис. 1. Граф перечисления исходов схемы примера 8 [Figure 1. Outcomes enumeration graph for the scheme in Example 8]
щем порядке номеров ее элементов.7 Провести прямое перечисление исходов схемы и решить ЗН.
Построим процедуру перечисления исходов схемы типа В. Это значит, что исходы схемы будут включать все исходы схемы сочетаний с размахами выборок ^ 4.8
ПП исходов схемы проводится 2-мя итерациями:
7Размах исхода схемы сочетаний равен разности номеров последнего и первого элементов.
8Номера элементов называем их значениями. 582
- первая перечисляет допустимые поединично растущие значения первых элементов и для каждого из них — поединично растущие допустимые значения третьих элементов;
- вторая производит варьирование возможных значений вторых элементов в диапазоне: больше первого и меньше третьего их значений.
Получаем следующие размеры пучков первой итерации: =6, N2 = 6, N3 = 6, N4 = 6, N5 = 6, N6 = 3, N7 = 1, составляющих исходы итоговой итерации — всего 34 исхода.
Представим это графом на рис. 2 в качестве исходов схемы 7-ми пучков начальной итерации в указанном выше порядке подряд с общей нумерацией.
Теперь по результатам ЗН для ОПД (п. 3.2) и имея полный перечень всех исходов схемы и их пучковую структуру, будем решать ЗН, сверяя результаты с рис. 2.
ПЗН. Номер исхода схемы ЫЕ = 14. Требуется найти его вид В*Е. По рис. 2 В*Е = (3, 4, 6). Определим его из условий схемы ОПД по теореме 1:
1) найдем номер к1 пучка искомого исхода схемы к1 = 3, т.к. 6 + 6 = 12 < < 14 < 6 + 6 + 6 = 18;
2) его номер к2 в пучке найдем из к2 = 14 — (6 + 6) = 2, вычитая из
= 14 суммарный размер двух предшествующих пучков итерации;
3) по результату решения ПЗН в схеме третьего пучка находим его вид второго исхода, являющийся итоговым для изучаемой схемы КЕ = = (3, 4, 6), что совпадает с результатом по рис. 2.
ОЗН. Вид исхода схемы В*Е = (3, 4, 6). Требуется найти его номер ИЕ. По рис. 2 ЫЕ = 14. Найдем его из условий схемы по теореме 2:
1) по первому элементу вида В*Е = (3, 4, 6) находим номер к1 =3 пучка искомого исхода схемы;
2) его номер к2 = 2 в пучке находим из результата обратной ЗН в схеме третьего пучка;
3) вычисляем искомый номер КЕ = (6+6)+2 = 14, суммируя номер к2 = 2 исхода в пучке с суммарным размером предшествующих пучков этой итерации, что совпадает с результатом по рис. 2.
6.2. Второй подход к решению ЗН. Исходы схемы Е могут стоять на любых местах в перечислении исходов схемы 5.
Прямая ЗН. Дан ЫЕ, требуется найти В*Е.
Алгоритм 1.
1. Проверяем исходы схемы 5 на выполнение в них условия схемы Е и считаем число А таких проверок, просматривая их с начала подряд до МЕ-го выполнения включительно.
2. По результату ПЗН в схеме 5 находим в ней вид исхода с номером
= А, совпадающий с искомым видом К*Е в схеме 5 после его записи в форме исходов в схеме Е.
Обратная ЗН. Дан В*Е, требуется найти ИЕ.
Алгоритм 2.
1. Записываем все исходы схемы 5 в форме исходов схемы Е.
2. При просмотре исходов схемы 5 подряд с начала их перечисления сравниваем их с данным видом исхода КЕ, считая число А таких сравнений до совпадения К*3 = КЕ (включительно).
(2)
(3)
(4)
(5)
of out-
comes
1 2 3) 1
1 2 4) 2
1 3 4) 3
1 2 5) 4
1 3 5) 5
1 4 5) 6
2 3 4) 7
2 3 5) 8
2 4 5) 9
2 3 6) 10
2 4 6) 11
2 5 6) 12
3 4 5) 13
3 4 6) 14
3 5 6) 15
3 4 7) 16
3 5 7) 17
3 6 7) 18
4 5 6) 19
4 5 7) 20
4 6 7) 21
4 5 8) 22
4 6 8) 23
4 7 8) 24
5 6 7) 25
5 6 8) 26
5 7 8) 27
5 6 9) 28
5 7 9) 29
5 8 9) 30
6 7 8) 31
6 7 9) 32
6 8 9) 33
7 8 9) 34
numbers
Рис. 2. Граф перечисления исходов схемы примера 9 [Figure 2. Outcomes enumeration graph for the scheme in Example 9]
3. Первые А исходов схемы 5 проверяем на выполнение в них условия схемы Е и считаем число таких выполнений совпадающим с искомым номером Ы*.
7. Вероятностное распределение исходов схем Е. Вероятностное распределение исходов схемы Е пересчитывается из вероятностного распределением схемы 5 по вероятностям траекторий графа перечисления исходов схемы 5, ведущих к исходам схемы Е после удаления из начального графа лишних для Е траекторий и пропорционального пересчета по пучкам каждой итерации вероятностей итерационных переходов. Эта процедура описана в перечислительном методе в [6,7].
8. Моделирование исходов схем Е. При решенной ЗН для схем Е с учетом специфики каждой при определенном распределении вероятностей всех исходов схемы, по которому и результатам ПЗН может быть проведено универсальное моделирование (УМ) исходов схемы. Оно состоит в том, что для моделирования каждого исхода его номер (в перечислении) разыгрывается случайным числом с заданным вероятностным распределением, а по нему из результата решения ПЗН определяется его смоделированный вид. Эта процедура описана в перечислительном методе в [6,7].
9. Примеры областей применения схем Е с примером анализа
по ПМ. Примеры областей применения схем Е:
1) в математических исследованиях, например, а) деления целого числа на определенное количество целых слагаемых с их заданным минимальным значением; б) деления целого числа на целые слагаемые в заданном минимальном количестве;
2) при составлении расписаний разного рода мероприятий с заданными экстремальными значениями их характеристик, например, недельного расписания определенного числа аудиторных часов занятий преподавателя с заданной максимальной нормой присутственных часов на работе в течение рабочего дня института;
3) в криминалистике при проведении следственных действий по свидетельским показаниям, связанным с экстремальными значениями характеристик числовых последовательностей, например, для построения полного перечисления номеров автомобилей с заданной минимальной цифрой.
Пример анализа схемы 1 (тип А) по ПМ. Рассматривается схема деления целого числа г на п слагаемых с каждым слагаемым ^ к при г ^ гп. В наглядных терминах размещения частиц по ячейкам она соответствует схеме размещения г неразличимых частиц по п неразличимым ячейкам с минимальным уровнем заполнения ячеек к при г ^ гп. Вид исхода схемы представляет собой перечисление уровней заполнения ячеек г? = (щ,... ..., цп) (это учитывает неразличимость частиц) в заранее заданном порядке, например, в порядке неубывания этих уровней (это учитывает неразличимость ячеек).
Алгоритм 3 перечисления исходов схемы.
1. Во все ячейки кладем по к неразличимых частиц, а остальные г — пк неразличимые частицы размещаем всеми N * = способами
по схеме сочетаний с повторением по различимым ячейкам со второй ячейки до п-й (этим на первом месте оставляем минимальный уро-
вень заполнения ячеек, равный к). В результате получаем места V = = (ь1,... ,Ьг,... ,ьп-2) расстановки (п — 2) внутренних границ ячеек между лежащими в ряд г* = г — пк частицами, т.е. среди г — пк + п — 2 общего числа мест границ и частиц.
2. Переводим результат п. 1 в термины ?у по формулам
щ = к, щ = к + ^г+1 — Ví — 1, г = 2,п.
3. Упорядочиваем результаты п. 2 по неубыванию уровней заполнения ячеек с отбраковкой повторов — получаем бесповторное перечисление исходов нашей схемы.
Теорема 3. Число N исходов схемы определяется по формуле
N = £ I(щ),
( (Д})
где суммирование (ОПП) производится по перечислению маркировок составов уровней заполнения ячеек в исходах схемы по их совпадениям:
Ц = (^1,...,^з), « ^ п,
а I= 1 при Е > 0 и I(2) = 0 при 2 = 0.
Доказательство. Искомое N есть число разных вышеописанных составов уровней заполнения ячеек — результатов п. 3 алгоритма 3 упорядоченных наборов уровней заполнения ячеек в исходах схемы сочетаний с повторением с минимальным уровнем заполнения ячеек, равным к. Оно совпадает с числом ненулевых частот при их маркировке. Тогда число N исходов схемы записывается доказываемой формулой. □
Остальные направления исследования схемы проводятся по общим вышеприведенным методикам.
Разберем численный пример расчета данной схемы по всем направлениям ПМ.
Пример 10. Пусть в нашей схеме п = 3, г = 10, к = 2.
Проведем ПП исходов схемы по шагам алгоритма 3:
1) г — пк = 4, N * = С|+2-1 = С5 =5 — это исходы мест 1, 2, 3, 4, 5 одной внутренней границы двух ячеек;
2) приводящие соответственно к исходам схемы в форме наборов уровней заполнения ячеек:
(2, 2, 2) + (0, 0, 4) = (2, 2, 6); (2, 2, 2) + (0,1, 3) = (2, 3, 5);
(2, 2, 2) + (0, 2, 2) = (2, 4, 4); (2, 2, 2) + (0, 3,1) = (2, 5, 3);
(2, 2, 2) + (0, 4, 0) = (2, 6, 2);
3) после требуемого упорядочивания с отбраковкой повторов получаем N = 3 исхода схемы: (2, 2, 6), (2, 3, 5), (2, 4, 4).
Вычислим значение N по теореме 3. Для этого промаркируем упорядоченные исходы схемы по их совпадениям — получим частоты (2, 2,1) соответствующих повторяющихся исходов схемы (2, 2, 6), (2, 3, 5) и (2, 4, 4), откуда по теореме 3 получаем N = I(2) +1(2) +1(1) = 3, что совпало с полученным выше результатом.
Решение ЗН. Применим прием решения ЗН для данной схемы Е из результатов ее решения в более общей схеме 5 с тем же минимальным уровнем заполнения ячеек без обязательного его достижения в каждом исходе и с учетом взаимных порядков ячеек. Перечисление ее исходов определяется в дополнительной схеме И уровнями заполнения трех различимых ячеек при размещении по ним четырех неразличимых частиц по схеме сочетаний с повторением с числом исходов С2 = 15. Перечисление этих исходов в порядке перечисления исходов схемы сочетаний (см. [7]) означает расстановку всеми способами двух внутренних границ ячеек между лежащими в ряд четырьмя неразличимыми частицами, т.е. выбор двух мест среди шести мест для четырех частиц и двух границ. В терминах этих мест записываем исходы:
(1, 2), (1, 3), (1, 4), (1, 5), (1, 6),
(2, 3), (2, 4), (2, 5), (2, 6), (3, 4),
(3, 5), (3, 6), (4, 5), (4, 6), (5, 6).
Переводим их в термины уровней заполнения ячеек по формулам п. 2 алгоритма 3 — получаем в порядке перечисления исходов схемы сочетаний видов
(0, 0, 4), (0,1, 3), (0, 2, 2), (0, 3,1), (0, 4, 0), (1, 0, 3), (1,1, 2), (1, 2,1), (1, 3, 0), (2, 0, 2), (2,1,1), (2, 2, 0), (3, 0,1), (3,1, 0), (4, 0, 0)
в этой дополнительной схеме И. Итоговые исходы схемы 5 в виде уровней заполнения трех ячеек в порядке ячеек находим сложением по трем ячейкам уровней в исходах схемы И с предварительными уровнями по две частицы в каждой из трех ячеек. Приходим к перечислению исходов в схеме 5:
(2, 2, 6), (2, 3, 5), (2, 4, 4), (2, 5, 3), (2, 6, 2), (3, 2, 5), (3, 3, 4), (3, 4, 3), (3, 5, 2), (4, 2, 4), (4, 3, 3), (4, 4, 2), (5, 2, 3), (5, 3, 2), (6, 2, 2).
Прямая ЗН. Дан Ы** = 3, требуется найти В*Е.
По алгоритму 1:
1) А = Ж* = М**, = 3;
2) отсюда по ПЗН для исходов схемы 5 получаем вид третьего исхода схемы 5: К*3 = В*Е = (2, 4, 4), что совпадает с полученным выше по ПП результатом.
Обратная ЗН. Дан В*Е = (2, 4, 4), требуется найти .
По алгоритму 2:
1) получаем исходы схемы 5 в порядке неубывания уровней заполнения ячеек:
(2, 2, 6), (2, 3, 5), (2, 4, 4), (2, 3, 5), (2, 2, 6), (2, 3, 5), (3, 3, 4), (3, 3, 4), (2, 3, 5), (2, 4, 4), (3, 3, 4), (2, 4, 4), (2, 3, 5), (2, 3, 5), (2, 2, 6);
2) А = 3;
3) в первых трех исходах выполняется условие схемы Е, откуда получаем = 3, что совпадает с полученным выше по ПП результатом.
Вероятностное распределение исходов схемы. Определим вероятности исходов нашей схемы Е при равновероятных исходах схемы 5. Исходы схемы Е отличаются от исходов схемы 5 обязательным достижением в каждом исходе минимального значения 2 уровня заполнения ячеек и несущественностью порядка перечисления уровней их заполнения. Поэтому из перечисления исходов схемы Б с вероятностями исходов по 1/С| = 1/15, однозначно соответствующих исходам схемы И, удаляем исходы, полученные из них без пустых ячеек. Это исходы
(3,3, 4), (3,4, 3), (4,3, 3)
без достижения заданного минимального уровня заполнения ячеек, равного двум. В результате пропорционального пересчета по [6] получаем вероятности оставшихся исходов схемы 5 по 1/12. Для учета несущественности порядка перечисления вероятностей составов уровней заполнения ячеек в исходах схемы Е для их вычисления кратно (числу разных порядков каждого состава этих уровней в схеме 5) увеличиваем их вероятности:
Р(2, 2, 6) = 3 ■ (1/12) = 0.25; Р(2, 3, 5) = 3! ■ (1/12) = 0.5;
Р(2, 4, 4) = 3 ■ (1/12) = 0.25. Проверка на распределение: Р(2, 2, 6) + Р(2, 3, 5) + Р(2, 4, 4) = 1.
Моделирование исхода схемы. Разыгрываем номер X исхода нашей схемы Е в порядке их перечисления с полученным вероятностным распределением по ряду распределения:
Р (X = 1) = 0.25; Р (X = 2) = 0.5; Р (X = 3) = 0.25.
Для этого датчиком случайных чисел генерируем значение случайной величины с равномерным распределением на отрезке [0,1]. Пусть это число 0.78 € (0.75,1]). По результату решения ПЗН это значит, что смоделирован третий исход нашей схемы К = (2, 4, 4).
10. Заключение. Введен новый класс обратных экстремальных задач комбинаторики на комбинаторных схемах Е с фиксированными экстремальными значениями характеристик двух типов, относящихся к исходам и сравниваемых внутри каждого исхода, или — к совокупности исходов и сравниваемых между исходами.
Анализ схем Е проведен новым перечислительным методом (ПМ) (см. [6]), состоящим в организации получения качественной информации об их исходах и ее переводе в количественные результаты по направлениям: явного прямого бесповторного нумерованного перечисления исходов схем, вычисления их числа, решения для них ЗН в прямой и обратной постановках, установления их вероятностного распределения и построения процедуры их моделирования.
В схемах Е предложено два подхода к решению ЗН: прямым перечислением их исходов по схеме ОПД (см. [6]) и пересчетом результатов решения ЗН для исходов соответствующих начальных схем 5.
Моделирование исходов схемы из Е предложено проводить в ПМ универсальным моделированием (УМ) с вероятностным распределением ее исходов, найденным из установленного случайным процессом соответствующего веро-ятностого распределения исходов в ее начальной схеме 5 в ПМ, по результату решения прямой ЗН для исходов изучаемой схемы Е.
Обсуждены примеры областей использования схем Е. По всем направлениям ПМ проведен анализ конкретного примера решения математической задачи.
Конкурирующие интересы. Конкурирующих интересов не имею. Авторский вклад и ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.
Финансирование. Исследование выполнялось без финансирования. Библиографический список
1. Афанасьев М. Ю., Суворов Б. П. Исследование операций в экономике. Модели, задачи, решения. М.: ИНФРА-М, 2003. 202 с.
2. Баранов В. И., Стечкин Б. С. Экстремальные комбинаторные задачи и их приложения. М.: Физматлит, 2004. 240 с.. ЕБЫ: ЕХОТЦО.
3. Колчин В. Ф. О предельном поведении крайних членов вариационного ряда в полиномиальной схеме // Теория вероятн. и ее примен., 1969. Т. 14, №3. С. 476-487.
4. Хакимуллин Е. Р. О предельном поведении максимального заполнения в равновероятной схеме размещения частиц комплектами // Матем. заметки, 1981. Т. 30, №2. С. 277289.
5. Викторова И. И. Об асимптотическом поведении максимума в равновероятной полиномиальной схеме// Матем. заметки, 1969. Т. 5, №3. С. 305-316.
6. Энатская Н. Ю. Вероятностные модели комбинаторных схем // Вестн. ЮУрГУ. Сер. Матем. моделирование и программирование, 2020. Т. 13, №3. С. 103-111. ЕБЫ: 0КМ1иС. 001: https://doi.org/10.14529/mmp200312.
7. Энатская Н. Ю. Анализ комбинаторных схем в доасимптотической области изменения параметров// Труды Карельского научного центра РАН, 2018. №7. С. 117-133. ЕБЫ: хадгга. 001: https://doi.org/10.17076/mat750.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2022, vol. 26, no. 3, pp. 573-591 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1947
MSC: 60F15
Probabilistic models for the analysis of inverse extremal problems in combinatorics
N. Yu. Enatskaya
National Research University "Higher School of Economics", Moscow Institute of Electronics and Mathematics, 34, Tallinskay st, Moscow, 123458, Russian Federation.
Abstract
In an inverse extremal problem for a combinatorial scheme with a given value of the objective function of the form of a certain extreme value of its characteristic, a probabilistic model is developed that ensures that this value is obtained in its outcomes. Two types of such characteristics are considered, relating to each of the schemes or to a set of outcomes.
The pre-asymptotic analysis of such a model is carried out by the author's enumerative method. It is based on the construction of an iterative random process with iterations of successive stages of a numbered non-repetitive enumeration and the formation of outcomes of the scheme. The iterative development of the process is represented by a probabilistic graph.
The study of the outcomes of the scheme according to the model in the enumerative method is carried out in the following areas: visual numbering of the outcomes of the scheme, finding their number, establishing a one-to-one correspondence between the types and numbers of outcomes of the scheme, obtaining their probabilistic distribution (controlled by a random process of listing the outcomes of the scheme), and modeling them with this distribution.
Along with the direct study of circuits in these areas, algorithms are proposed to obtain results for them by partially recalculating them from the results of a similar analysis of more general, previously studied circuits without restrictions or with less restrictions on the values of the characteristics under consideration.
Keywords: inverse extremal problem, extremal value of a characteristic, pre-asymptotic analysis of a circuit.
Received: 12th August, 2022 / Revised: 25th August, 2022 / Accepted: 31st August, 2022 / First online: 20th September, 2022
Mathematical Modeling, Numerical Methods and Software Complexes Research Article
© Authors, 2022
© Samara State Technical University, 2022 (Compilation, Design, and Layout) 3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
EnatskayaN. Yu. Probabilistic models for the analysis of inverse extremal problems in combinatorics, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2022, vol. 26, no. 3, pp. 573-591. EDN: AWBYGO. DOI: 10.14498/vsgtu1947 (In Russian). Authors' Details:
Nataliya Yu. Enatskaya https://orcid.org/0000-0003-1241-7543 Cand. Phys. & Math. Sci.; Associate Professor; Dept. of Applied Mathematics; e-mail: nat1943@mail.ru
Competing interests. I have no competing interests.
Author's Responsibilities. I take full responsibility for submit the final manuscript
to print. I approved the final version of the manuscript.
Funding. This research received no external funding.
References
1. Afanas'ev M. Yu., Suvorov B. P. Issledovanie operatsii v ekonomike. Modeli, zadachi, resh-eniia [Operations Research in Economics. Models, Tasks, Solutions]. Moscow, INFRA-M, 2003, 202 pp. (In Russian)
2. Baranov V. I., Stechkin B. S. Ekstremal'nye kombinatornye zadachi i ikh prilozheniia [Extremal Combinatorial Problems and Their Applications]. Moscow, Fizmatlit, 2004, 240 pp. (In Russian). EDN: RXGTQD.
3. Kolchin V. F. On the limiting behavior of extreme order statistics in a polynomial scheme, Theory Probab. Appl., 1969, vol.14, no. 3, pp. 458-469. DOI: https://doi.org/ 10.1137/1114058.
4. Khakimullin E.R. Asymptotic behavior of the maximum occupancy in an equiprobable scheme of allocation of particles by complexes, Math. Notes, 1981, vol. 30, no. 2, pp. 626633. DOI: https://doi.org/10.1007/BF01708846.
5. Viktorova I. I. Asymptotic behavior of maximum of an equiprobable polynomial scheme, Math. Notes, 1969, vol.5, no. 3, pp. 184-191. EDN: ZZYIEX. DOI: https://doi.org/10.1007/ BF01388624.
6. Enatskaya N. Yu. Probabilistic models of combinatorial schemes, Bulletin of the South Ural State University, Ser. Mathematical Modelling, Programming and Computer Software, 2020, vol.13, no. 3, pp. 103-111 (In Russian). EDN: DKMIUC. DOI: https://doi.org/ 10.14529/mmp200312.
7. Enatskaya N. Yu. Analysis of combinatorial schemes in the pre-asymptotic region of parameter change, Proceedings of the Karelian Research Centre of the Russian Academy of Sciences, 2018, no. 7, pp. 117-133 (In Russian). EDN: XRQZZR. DOI: https://doi.org/ 10.17076/mat750.