Получение следующего объекта — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
м (Специализация алгоритма для генерации следующего разбиения на подмножества)
(Специализация алгоритма для генерации следующего сочетания)
 
(не показано 167 промежуточных версий 16 участников)
Строка 1: Строка 1:
==Комбинаторные объекты==
 
[http://neerc.ifmo.ru/wiki/index.php?title=%D0%9A%D0%BE%D0%BC%D0%B1%D0%B8%D0%BD%D0%B0%D1%82%D0%BE%D1%80%D0%BD%D1%8B%D0%B5_%D0%BE%D0%B1%D1%8A%D0%B5%D0%BA%D1%82%D1%8B Комбинаторные объекты - Википедия]
 
  
 
== Алгоритм ==
 
== Алгоритм ==
Строка 8: Строка 6:
  
 
Отсюда понятен алгоритм:
 
Отсюда понятен алгоритм:
* Находим суффикс минимальной длины, который можно изменить без изменения префикса текущего объекта <tex>P</tex>
+
* находим суффикс минимальной длины, который можно изменить без изменения префикса текущего объекта <tex>P</tex>,
* К оставшейся части дописываем минимальный возможный элемент (чтобы было выполнено правило <tex>P < Q</tex>)
+
* к оставшейся части дописываем минимальный возможный элемент (чтобы было выполнено правило <tex>P < Q</tex>),
* Дописываем минимальный возможный хвост
+
* дописываем минимальный возможный хвост.
 
По построению получаем, что <tex>Q</tex> {{---}} минимально возможный.
 
По построению получаем, что <tex>Q</tex> {{---}} минимально возможный.
  
 
== Специализация алгоритма для генерации следующего битового вектора ==
 
== Специализация алгоритма для генерации следующего битового вектора ==
* Находим минимальный суффикс, в котором есть 0, его можно увеличить, не меняя оставшейся части
+
* Находим минимальный суффикс, в котором есть <tex>0</tex>, его можно увеличить, не меняя оставшейся части
* Вместо 0 записываем 1  
+
* Вместо <tex>0</tex> записываем <tex>1</tex>
 
* Дописываем минимально возможный хвост из нулей
 
* Дописываем минимально возможный хвост из нулей
  '''for''' i = n '''downto''' 1
+
'''int[]''' nextVector('''int[]''' a): <font color=green>// <tex>n</tex> {{---}} длина вектора</font>
      '''if''' a[i] == 0
+
  '''while''' (n >= 0) '''and''' (a[n] != 0)
          a[i] = 1
+
      a[n] = 0
          '''for''' j = i + 1 to n
+
      n--
              a[j] = 0
+
  '''if''' n == -1
          '''break'''
+
    '''return''' ''null''
 
+
  a[n] = 1
 +
  '''return''' a
 +
Приведённый алгоритм эквивалентен прибавлению единицы к битовому вектору.
 
=== Пример работы ===
 
=== Пример работы ===
 
{| class="wikitable" border = 1
 
{| class="wikitable" border = 1
|0||1||style="background:#FFCC00"|0||1||1||исходный битовый вектор
+
|0||1||0||1||style="background:#FFCC00"|1||исходный битовый вектор
 
|-
 
|-
| || ||^|| || ||находим элемент 0 (самый правый)
+
| || || || ||^|| начинаем идти с конца
 
|-
 
|-
|0||1||style="background:#FFCC00"|1||1||1||меняем его на 1
+
|0||1||0||style="background:#FFCC00"|0||style="background:#FFCC00"|0|| пока элементы равны 1, заменяем их на 0
 
|-
 
|-
|0||1||1||style="background:#FFCC00"|0||style="background:#FFCC00"|0||меняем элементы правее на нули
+
|0||1||style="background:#FFCC00"|1||0||0|| меняем первый не удовлетворяющий условию цикла элемент на 1
 
|-
 
|-
 
|'''0'''||'''1'''||'''1'''||'''0'''||'''0'''||следующий битовый вектор
 
|'''0'''||'''1'''||'''1'''||'''0'''||'''0'''||следующий битовый вектор
Строка 38: Строка 38:
  
 
== Специализация алгоритма для генерации следующей перестановки ==
 
== Специализация алгоритма для генерации следующей перестановки ==
* Двигаясь справа налево, находим элаемент, нарушающий убывающую последовательность (в обычном порядке, слева направо, см. пример)
+
* Двигаясь справа налево, находим элемент, нарушающий убывающую последовательность (в обычном порядке, слева направо, см. пример)
 
* Меняем его с минимальным элементом, большим нашего, стоящим правее
 
* Меняем его с минимальным элементом, большим нашего, стоящим правее
 
* Перевернем правую часть
 
* Перевернем правую часть
  
  '''for''' i = n - 1 '''downto''' 1
+
'''int[]''' nextPermutation('''int[]''' a): <font color=green>// <tex>n</tex> {{---}} длина перестановки</font>
      '''if''' a[i] < a[i + 1]
+
  '''for''' i = n - 2 '''downto''' 0
          // a[j] = min {a[j] > a[i], где j > i}
+
    '''if''' a[i] < a[i + 1]
          swap(a[i], a[j])
+
      min = i + 1;
          reverse(a[i + 1]..a[n])
+
      '''for''' j = i + 1 '''to''' n - 1
          '''break'''
+
        '''if''' (a[j] < a[min]) '''and''' (a[j] > a[i])
 +
          min = j
 +
      swap(a[i], a[min])
 +
      reverse(a, i + 1, n - 1)
 +
      '''return''' a
 +
  '''return''' ''null''  
  
 
=== Пример работы ===
 
=== Пример работы ===
Строка 68: Строка 73:
 
* Меняем его с минимальным элементом, большим нашего, стоящим правее.
 
* Меняем его с минимальным элементом, большим нашего, стоящим правее.
 
* Переворачиваем правую часть.
 
* Переворачиваем правую часть.
  '''procedure''' nextMultiperm(var b:array[1..N] of integer);
+
  '''int[]''' nextMultiperm('''int[]''' b): <font color=green>// <tex>n</tex> {{---}} длина мультиперестановки</font>
  '''var''' i , j : '''integer''';
+
    i = n - 2
  '''begin'''
+
    '''while''' (i >= 0) '''and''' (b[i] >= b[i + 1])  
    i := N -   1;
+
      i--
    '''while''' (i > 0) '''and''' (b[i] >= b[i + 1]) '''do'''
+
    '''if''' i >= 0  
    dec(i);
+
       j = i + 1
    '''if''' i > 0 '''then'''
+
       '''while''' (j < n - 1) '''and''' (b[j + 1] > b[i])  
    '''begin'''
+
         j++
       j := i + 1;
+
       swap(b[i] , b[j])
       '''while''' (j < N) '''and''' (b[j + 1] > b[i]) '''do'''
+
       reverse(b, i + 1, n - 1)
         inc(j);
+
       '''return''' b
       swap(b[i] , b[j]);
+
    '''else'''
       '''for''' j := i + 1 '''to''' (N + i) '''div''' 2 '''do'''
+
       '''return''' ''null''
        swap(b[j], b[N - j + i + 1]);
 
       write(b[1..N]);
 
    '''end'''
 
    '''else'''
 
    '''begin'''
 
       write('null');
 
    '''end;'''
 
  '''end;'''
 
  
 
=== Пример работы ===
 
=== Пример работы ===
Строка 105: Строка 102:
 
== Специализация алгоритма для генерации следующего сочетания ==
 
== Специализация алгоритма для генерации следующего сочетания ==
  
* Добавим в конец массива с сочетанием N+1 – максимальный элемент.
+
* Добавим в конец массива с сочетанием <tex>N+1</tex> – максимальный элемент.
* Пойдём справа налево. Будем искать номер элемента, который отличается от предыдущего на 2.
+
* Пойдём справа налево. Будем искать номер элемента, который отличается от предыдущего на <tex>2</tex> и больше.
* Увеличим найденный элемент на 1, и допишем в конец минимально возможный хвост, если такого элемента нет – на вход было дано последнее сочетание.
+
* Увеличим найденный элемент на <tex>1</tex>, и допишем в конец минимально возможный хвост, если такого элемента нет – на вход было дано последнее сочетание.
  '''procedure''' nextChoose(var a:array[1..N+1] of integer); // n,k - параметры сочетания.
+
  '''int[]''' nextChoose('''int[]''' a, '''int''' n, '''int''' k): <font color=green>// <tex>n,k </tex> {{---}} параметры сочетания</font>
  '''var''' i,j : '''integer;'''
+
  '''for''' i = 0 '''to''' k - 1
  '''begin'''
+
    b[i] = a[i]
   a[k + 1] := n + 1;
+
   b[k] = n + 1
   i := n;
+
   i = k - 1
   '''while''' (i > 0) '''and''' ((a[i + 1] - a[i]) < 2) '''do'''
+
   '''while''' (i >= 0) '''and''' (b[i + 1] - b[i] < 2)  
     i := i - 1;
+
     i--
   '''if''' i > 0 '''then'''
+
   '''if''' i >= 0  
    '''begin'''
+
      b[i]++
      a[i] := a[i] + 1;
+
      '''for''' j = i + 1 '''to''' k - 1
       '''for''' j := i + 1 '''to''' k '''do'''
+
        b[j] = b[j - 1] + 1
      a[j] := a[j - 1] + 1;
+
       '''for''' i = 0 '''to''' k - 1
       write(a[1..n]);
+
        a[i] = b[i]
    '''end'''
+
       '''return''' a
 
   '''else'''
 
   '''else'''
    write('null');
+
    '''return''' ''null''
  '''end;'''
 
  
 
=== Пример работы ===
 
=== Пример работы ===
Строка 142: Строка 138:
  
 
== Специализация алгоритма для генерации следующего разбиения на слагаемые ==
 
== Специализация алгоритма для генерации следующего разбиения на слагаемые ==
 
+
Рассматриваемый алгоритм находит следующее [[комбинаторные объекты|разбиение на слагаемые]], при этом разбиение упорядоченно по возрастанию.
* Увеличим предпоследнее число на 1, уменьшим последнее на 1.
+
* Увеличим предпоследнее слагаемое на <tex>1</tex>, уменьшим последнее слагаемое на <tex>1</tex>.
* Если предпоследний элемент стал больше последнего, то увеличиваем предпоследний элемент на величину последнего.
+
** Если предпоследнее слагаемое стало больше последнего, то увеличиваем предпоследнее слагаемое на величину последнего.
* Если предпоследний элемент меньше последнего, то проверяем, можем ли мы разбить последний элемент на сумму предпоследних. Если да – разбиваем, пока предпоследний*2 < последнего.
+
** Если предпоследнее слагаемое умноженное на 2 меньше последнего, то разбиваем последнее слагаемое <tex>s</tex> на два слагаемых <tex>a</tex> и <tex>b</tex> таких, что <tex>a</tex> равно предпоследнему слагаемому, а <tex>b = s - a</tex>. Повторяем этот процесс, пока разбиение остается корректным, то есть предпоследнее слагаемое хотя бы в два раза меньше последнего.
  
 
<code>
 
<code>
  // Алгоритм работает только для лексикографического порядка.
+
  <font color=green>// <tex>b</tex> {{---}} список, содержащий разбиение данного числа <tex>b.size</tex>{{---}} его размер </font>
// b – массив, содержащий разбиение данного числа, length – его размер.
+
  '''list<int>''' nextPartition('''list<int>''' b):  
  '''procedure''' nextPartition(var b:array[1..length] of integer);
+
     b[b.size - 1]--
  '''var''' i : '''integer;'''
+
     b[b.size - 2]++
  '''begin''' 
+
     '''if''' b[b.size - 2] > b[b.size - 1]  
     b[length] := b[length] - 1;
+
       b[b.size - 2] += b[b.size - 1]
     b[length - 1] := b[length - 1] + 1;
+
       b.remove(b.size - 1)
     '''if''' b[length - 1] > b[length] '''then'''
 
    '''begin'''
 
       b[length - 1] := b[length - 1] + b[length];
 
       length := length - 1;
 
    '''end'''
 
 
     '''else'''
 
     '''else'''
    '''begin'''
+
       '''while''' b[b.size - 2] * 2 <= b[b.size - 1]  
      i := 0;
+
         b.add(b[b.size - 1] - b[b.size - 2])
       '''while''' b[length - 1] * 2 <= b[length + i] '''do'''
+
         b[b.size - 2] = b[b.size - 3]
      '''begin'''
+
    '''return''' b
         b[length + i + 1] := b[length + i] - b[length - 1];
 
         b[length + i] := b[length - 1];
 
        i := i + 1;
 
      '''end;'''
 
      length := length + i;
 
    '''end;'''
 
    write(b[1..length];
 
  '''end;'''
 
 
</code>
 
</code>
  
Строка 179: Строка 162:
 
|1||style="background:#FFCC00"|1||style="background:#FFCC00"|7|| || ||Прибавим 1 + 1, вычтем 7 - 1.
 
|1||style="background:#FFCC00"|1||style="background:#FFCC00"|7|| || ||Прибавим 1 + 1, вычтем 7 - 1.
 
|-
 
|-
|1||style="background:#FFCC00"|2||style="background:#FFCC00"|6|| || ||Проверяем: 2<6, значит разбиваем 6 пока оно не станет <4
+
|1||style="background:#FFCC00"|2||style="background:#FFCC00"|6|| || ||Проверяем: 2 < 6, значит разбиваем 6 пока оно не станет меньше 4
 
|-
 
|-
 
|1||2||style="background:#FFCC00"|2||style="background:#FFCC00"|4|| ||
 
|1||2||style="background:#FFCC00"|2||style="background:#FFCC00"|4|| ||
Строка 186: Строка 169:
 
|-
 
|-
 
|'''1'''||'''2'''||'''2'''||'''2'''||'''2'''||Следующее разбиение на слагаемые числа 9.
 
|'''1'''||'''2'''||'''2'''||'''2'''||'''2'''||Следующее разбиение на слагаемые числа 9.
 +
|}
 +
 +
{| class="wikitable" border = 1
 +
|1||style="background:#FFCC00"|4||style="background:#FFCC00"|5||Прибавим 4 + 1, вычтем 5 - 1.
 +
|-
 +
|1||style="background:#FFCC00"|5||style="background:#FFCC00"|4||Проверяем: 5 > 4, значит прибавим к 5 + 4.
 +
|-
 +
|1||9||style="background:#FFCC00"|4||Удалим последний элемент.
 +
|-
 +
|'''1'''||'''9'''||||Следующее разбиение на слагаемые числа 10.
 
|}
 
|}
  
Строка 191: Строка 184:
  
 
Рассмотрим множество первых n натуральных чисел:<tex>N_n = \{1, 2, ..., n\}</tex>
 
Рассмотрим множество первых n натуральных чисел:<tex>N_n = \{1, 2, ..., n\}</tex>
 
{{Определение
 
|id=def1.
 
|definition='''Разбиением на множества''' называется представление множества, как объединения одного или более попарно
 
не пересекающихся подмножеств множеств.
 
}}
 
Например, для <tex>n = 5</tex> существуют следующие разбиения:
 
 
<tex> \{1, 2, 3, 4, 5\}</tex>
 
 
<tex> \{1, 2, 3\}~ \{4, 5\}</tex>
 
 
<tex> \{1, 3, 5\}~ \{2, 4\}</tex>
 
 
<tex> \{1\}~\{2\}~\{3\}~\{4\}~\{5\}</tex>
 
 
и т. д., всего таких разбиений для <tex>n = 5</tex> существует 52.
 
 
'''Примечание:'''
 
<tex> \{1, 2, 3\}~ \{4, 5\}</tex> и <tex>\{4, 5\} ~\{1, 2, 3\}</tex> - одно и то же разбиение на подмножества.
 
  
 
Упорядочим все разбиения на множества <tex>N_n</tex> лексикографически. Для этого, во-первых, в каждом разбиении упорядочим множества лексикографически. Будем говорить, что подмножество <tex> A \subset N_n </tex> лексикографически меньше подмножества <tex> B \subset N_n </tex> , если верно одно из следующих условий:
 
Упорядочим все разбиения на множества <tex>N_n</tex> лексикографически. Для этого, во-первых, в каждом разбиении упорядочим множества лексикографически. Будем говорить, что подмножество <tex> A \subset N_n </tex> лексикографически меньше подмножества <tex> B \subset N_n </tex> , если верно одно из следующих условий:
Строка 221: Строка 194:
  
 
'''Рассмотрим алгоритм нахождения лексикографически следующего разбиения на подмножества:'''
 
'''Рассмотрим алгоритм нахождения лексикографически следующего разбиения на подмножества:'''
*Будем хранить подмножества с помощью двумерного массива, например, разбиение <tex> \{1, 2, 3\}~ \{4, 5\}</tex> будет выглядеть так:
+
*Будем хранить подмножества в списке списков, например, разбиение <tex> \{1, 2, 3\}~ \{4, 5\}</tex> будет выглядеть так:
  
 
{| class="wikitable" border = 1
 
{| class="wikitable" border = 1
Строка 230: Строка 203:
  
 
* Двигаясь снизу вверх и справа налево, будем удалять элементы, записывая их в отдельный массив. Будем повторять эту операцию, пока не сможем выполнить одно из действий, описанных ниже:
 
* Двигаясь снизу вверх и справа налево, будем удалять элементы, записывая их в отдельный массив. Будем повторять эту операцию, пока не сможем выполнить одно из действий, описанных ниже:
** Заменить рассматриваемый элемент уже удаленным. Из всех подходящих элементов выбираем минимальный. '''Важное замечание''': мы не можем заменить 1ый элемент подмножества, мы можем только удалить его.
+
** Заменить рассматриваемый элемент уже удаленным. Из всех подходящих элементов выбираем минимальный. '''Важное замечание''': мы не можем заменить первый элемент подмножества, мы можем только удалить его.
 
** Дополнить рассматриваемое подмножество уже удаленным элементом. Из всех подходящих элементов выбираем минимальный.
 
** Дополнить рассматриваемое подмножество уже удаленным элементом. Из всех подходящих элементов выбираем минимальный.
 
* Допишем лексикографически минимальный хвост подмножеств из оставшихся элементов.
 
* Допишем лексикографически минимальный хвост подмножеств из оставшихся элементов.
  
 
<code>
 
<code>
  // a - матрица, содержащая подмножества
+
  '''list<list<int>>''' nextSetPartition('''list<list<int>>''' a):
// used - массив, в котором мы храним удаленные элементы
+
  <font color=green>// <tex>a</tex> {{---}} список, содержащий подмножества</font>
fl = ''false''
+
  <font color=green>// <tex>used</tex> {{---}} список, в котором мы храним удаленные элементы</font>
'''for''' i = n - 1 '''downto''' 0
+
  used = '''list<int>'''
    '''if''' можем добавить в конец подмножества элемент из used
+
  fl = ''false''
        добавляем
+
  '''for''' i = a.size - 1 '''downto''' 0
        '''break'''
+
      '''if''' (used.size != 0) '''and''' (used[used.size - 1] > a[i][a[i].size - 1])  <font color=green>// если можем добавить в конец подмножества элемент из <tex>used</tex></font>
    '''for''' j = a[i].size - 1 '''downto''' 0
+
          a[i].add(used[used.size - 1])  <font color=green>//добавляем</font>
        '''if''' можем заменить элемент, другим элементом из массива used  
+
          used.remove(used.size - 1)
            заменяем
+
          '''break'''
            fl = ''true''
+
      '''for''' j = a[i].size - 1 '''downto''' 0
            '''break'''
+
          '''if''' (used.size != 0) '''and''' (j != 0) '''and''' (used[used.size - 1] > a[i][j])    <font color=green>//если можем заменить элемент, другим элементом из списка <tex>used</tex> </font>
        used.add(a[i][j])  // удаляем элемент и добавляем его в массив
+
            a[i][j] = used[used.size - 1]  <font color=green>//заменяем</font>
    '''if''' (fl) '''break'''
+
            fl = ''true''
// далее выведем все получившиеся подмножества
+
            '''break'''
sort(used)
+
      '''if''' fl '''break'''
'''for''' i = 0 '''to''' used.size - 1
+
      used.add(a[i][j])  <font color=green>//добавляем в <tex>used</tex> <tex>j</tex> элемент <tex>i</tex>-го подмножества</font>
    println(used[i])       // выводим лексикографически минимальных хвост
+
      a[i].remove(j)   <font color=green>//удаляем <tex>j</tex> элемент <tex>i</tex>-го подмножества</font>
 +
  <font color=green>//далее выведем все получившиеся подмножества</font>
 +
  sort(used)
 +
  '''for''' i = 0 '''to''' used.size - 1
 +
    a.add('''list<int>'''(used[i]))  <font color=green>//добавляем лексикографически минимальных хвост</font>
 +
  '''return''' a
 
</code>  
 
</code>  
  
Строка 312: Строка 290:
 
|}
 
|}
  
== Ссылки ==
+
== См.также ==
 +
* [[Получение предыдущего объекта]]
 
* [[Получение объекта по номеру]]
 
* [[Получение объекта по номеру]]
 
* [[Получение номера по объекту]]
 
* [[Получение номера по объекту]]
 +
 +
== Источники информации ==
 +
 
* [http://rain.ifmo.ru/cat/view.php/vis/combinations/permutations-2000 Визуализатор перестановок]
 
* [http://rain.ifmo.ru/cat/view.php/vis/combinations/permutations-2000 Визуализатор перестановок]
 
* [http://cppalgo.blogspot.com/2011/02/episode-2.html Пример компактного кода для перестановок (С++)]
 
* [http://cppalgo.blogspot.com/2011/02/episode-2.html Пример компактного кода для перестановок (С++)]

Текущая версия на 01:25, 16 декабря 2017

Алгоритм[править]

Определение:
Получение следующего объекта — это нахождение объекта, следующего за данным в лексикографическом порядке.

Объект [math]Q[/math] называется следующим за [math]P[/math], если [math]P \lt Q[/math] и не найдется такого [math]R[/math], что [math]P \lt R \lt Q[/math].

Отсюда понятен алгоритм:

  • находим суффикс минимальной длины, который можно изменить без изменения префикса текущего объекта [math]P[/math],
  • к оставшейся части дописываем минимальный возможный элемент (чтобы было выполнено правило [math]P \lt Q[/math]),
  • дописываем минимальный возможный хвост.

По построению получаем, что [math]Q[/math] — минимально возможный.

Специализация алгоритма для генерации следующего битового вектора[править]

  • Находим минимальный суффикс, в котором есть [math]0[/math], его можно увеличить, не меняя оставшейся части
  • Вместо [math]0[/math] записываем [math]1[/math]
  • Дописываем минимально возможный хвост из нулей
int[] nextVector(int[] a): // [math]n[/math] — длина вектора
  while (n >= 0) and (a[n] != 0)
      a[n] = 0
      n--
  if n == -1
    return null
  a[n] = 1
  return a

Приведённый алгоритм эквивалентен прибавлению единицы к битовому вектору.

Пример работы[править]

0 1 0 1 1 исходный битовый вектор
^ начинаем идти с конца
0 1 0 0 0 пока элементы равны 1, заменяем их на 0
0 1 1 0 0 меняем первый не удовлетворяющий условию цикла элемент на 1
0 1 1 0 0 следующий битовый вектор

Специализация алгоритма для генерации следующей перестановки[править]

  • Двигаясь справа налево, находим элемент, нарушающий убывающую последовательность (в обычном порядке, слева направо, см. пример)
  • Меняем его с минимальным элементом, большим нашего, стоящим правее
  • Перевернем правую часть
int[] nextPermutation(int[] a): // [math]n[/math] — длина перестановки
  for i = n - 2 downto 0
    if a[i] < a[i + 1]
      min = i + 1;
      for j = i + 1 to n - 1
        if (a[j] < a[min]) and (a[j] > a[i])
          min = j
      swap(a[i], a[min])
      reverse(a, i + 1, n - 1)
      return a
  return null 

Пример работы[править]

1 3 2 5 4 исходная перестановка
^ находим элемент, нарушающий убывающую последовательность
^ минимальный элемент больше нашего
1 3 4 5 2 меняем их местами
1 3 4 2 5 разворачивам правую часть
1 3 4 2 5 следующая перестановка

Специализация алгоритма для генерации следующей мультиперестановки[править]

  • Двигаясь справа налево, находим элемент, нарушающий убывающую последовательность (в обычном порядке, слева направо, см. пример).
  • Меняем его с минимальным элементом, большим нашего, стоящим правее.
  • Переворачиваем правую часть.
int[] nextMultiperm(int[] b):  // [math]n[/math] — длина мультиперестановки
    i = n - 2
    while (i >= 0) and (b[i] >= b[i + 1]) 
      i--
    if i >= 0 
      j = i + 1
      while (j < n - 1) and (b[j + 1] > b[i]) 
        j++
      swap(b[i] , b[j])
      reverse(b, i + 1, n - 1)
      return b
    else
      return null

Пример работы[править]

1 2 3 1 2 3 Исходная перестановка.
^ Находим элемент, нарушающий убывающую последовательность.
^ Минимальный элемент больше нашего.
1 2 3 1 3 2 Меняем их местами.
1 2 3 1 3 2 Следующая мультиперестановка.

Специализация алгоритма для генерации следующего сочетания[править]

  • Добавим в конец массива с сочетанием [math]N+1[/math] – максимальный элемент.
  • Пойдём справа налево. Будем искать номер элемента, который отличается от предыдущего на [math]2[/math] и больше.
  • Увеличим найденный элемент на [math]1[/math], и допишем в конец минимально возможный хвост, если такого элемента нет – на вход было дано последнее сочетание.
int[] nextChoose(int[] a, int n, int k): // [math]n,k [/math] — параметры сочетания
  for i = 0 to k - 1 
    b[i] = a[i]
  b[k] = n + 1
  i = k - 1
  while (i >= 0) and (b[i + 1] - b[i] < 2) 
    i--
  if i >= 0 
     b[i]++
     for j = i + 1 to k - 1 
       b[j] = b[j - 1] + 1
     for i = 0 to k - 1 
       a[i] = b[i]
     return a
  else
    return null

Пример работы[править]

1 2 5 6 7 Дописываем 7 в конец сочетания.
1 2 5 6 7
^ Находим элемент i, a[i + 1] - a[ i ] >= 2
1 3 5 6 7 Увеличиваем его на 1.
1 3 4 5 6 Дописываем минимальный хвост.
1 3 4 5 Следующее сочетание.

Специализация алгоритма для генерации следующего разбиения на слагаемые[править]

Рассматриваемый алгоритм находит следующее разбиение на слагаемые, при этом разбиение упорядоченно по возрастанию.

  • Увеличим предпоследнее слагаемое на [math]1[/math], уменьшим последнее слагаемое на [math]1[/math].
    • Если предпоследнее слагаемое стало больше последнего, то увеличиваем предпоследнее слагаемое на величину последнего.
    • Если предпоследнее слагаемое умноженное на 2 меньше последнего, то разбиваем последнее слагаемое [math]s[/math] на два слагаемых [math]a[/math] и [math]b[/math] таких, что [math]a[/math] равно предпоследнему слагаемому, а [math]b = s - a[/math]. Повторяем этот процесс, пока разбиение остается корректным, то есть предпоследнее слагаемое хотя бы в два раза меньше последнего.

// [math]b[/math] — список, содержащий разбиение данного числа [math]b.size[/math]— его размер 
list<int>  nextPartition(list<int> b): 
   b[b.size - 1]--
   b[b.size - 2]++
   if b[b.size - 2] > b[b.size - 1] 
      b[b.size - 2] += b[b.size - 1]
      b.remove(b.size - 1)
   else
     while b[b.size - 2] * 2 <= b[b.size - 1] 
       b.add(b[b.size - 1] - b[b.size - 2])
       b[b.size - 2] = b[b.size - 3]
   return b

Пример работы[править]

1 1 7 Прибавим 1 + 1, вычтем 7 - 1.
1 2 6 Проверяем: 2 < 6, значит разбиваем 6 пока оно не станет меньше 4
1 2 2 4
1 2 2 2 2
1 2 2 2 2 Следующее разбиение на слагаемые числа 9.
1 4 5 Прибавим 4 + 1, вычтем 5 - 1.
1 5 4 Проверяем: 5 > 4, значит прибавим к 5 + 4.
1 9 4 Удалим последний элемент.
1 9 Следующее разбиение на слагаемые числа 10.

Специализация алгоритма для генерации следующего разбиения на подмножества[править]

Рассмотрим множество первых n натуральных чисел:[math]N_n = \{1, 2, ..., n\}[/math]

Упорядочим все разбиения на множества [math]N_n[/math] лексикографически. Для этого, во-первых, в каждом разбиении упорядочим множества лексикографически. Будем говорить, что подмножество [math] A \subset N_n [/math] лексикографически меньше подмножества [math] B \subset N_n [/math] , если верно одно из следующих условий:

  • существует [math]i[/math] такое, что [math]i \in A[/math] , [math]i \notin B[/math], для всех [math]j \lt i: j \in A[/math] если и только если [math]j \in B[/math] , и существует [math]k \gt i[/math] такое что [math]k \in B[/math];
  • [math] A \subset B [/math] и [math]i \lt j[/math] для всех [math]i \in A[/math] и [math]j \in B[/math] \ [math] A [/math].

Разбиения упорядочены лексикографически следующим образом. Разбиение [math]N_n = A_1 \cup A_2 \cup . . . \cup A_k[/math] лексикографически меньше разбиения [math]N_n = B_1 \cup B_2 \cup . . . \cup B_l[/math] если существует такое [math]i[/math], что [math]A_1 = B_1, A_2 = B_2, . . . ,A_{i - 1} = B_{i - 1}, A_i \lt B_i[/math].


Рассмотрим алгоритм нахождения лексикографически следующего разбиения на подмножества:

  • Будем хранить подмножества в списке списков, например, разбиение [math] \{1, 2, 3\}~ \{4, 5\}[/math] будет выглядеть так:
1 2 3
4 5
  • Двигаясь снизу вверх и справа налево, будем удалять элементы, записывая их в отдельный массив. Будем повторять эту операцию, пока не сможем выполнить одно из действий, описанных ниже:
    • Заменить рассматриваемый элемент уже удаленным. Из всех подходящих элементов выбираем минимальный. Важное замечание: мы не можем заменить первый элемент подмножества, мы можем только удалить его.
    • Дополнить рассматриваемое подмножество уже удаленным элементом. Из всех подходящих элементов выбираем минимальный.
  • Допишем лексикографически минимальный хвост подмножеств из оставшихся элементов.

list<list<int>> nextSetPartition(list<list<int>> a):
 // [math]a[/math] — список, содержащий подмножества
 // [math]used[/math] — список, в котором мы храним удаленные элементы
 used = list<int>
 fl = false
 for i = a.size - 1 downto 0
     if (used.size != 0) and (used[used.size - 1] > a[i][a[i].size - 1])   // если можем добавить в конец подмножества элемент из [math]used[/math]
         a[i].add(used[used.size - 1])   //добавляем
         used.remove(used.size - 1)
         break
     for j = a[i].size - 1 downto 0
         if (used.size != 0) and (j != 0) and (used[used.size - 1] > a[i][j])    //если можем заменить элемент, другим элементом из списка [math]used[/math] 
            a[i][j] = used[used.size - 1]   //заменяем
            fl = true
            break
     if fl break
     used.add(a[i][j])   //добавляем в [math]used[/math] [math]j[/math] элемент [math]i[/math]-го подмножества 
     a[i].remove(j)   //удаляем [math]j[/math] элемент [math]i[/math]-го подмножества
 //далее выведем все получившиеся подмножества
 sort(used)
 for i = 0 to used.size - 1
    a.add(list<int>(used[i]))   //добавляем лексикографически минимальных хвост
 return a

Пример работы[править]

Рассмотрим следующее разбиение:

1 2 3
4 5

1 Шаг:

1 2 3
4 5
^ Удалили элемент 5.
used


2 Шаг:

1 2 3
4
^ Удалили элемент 4. Так как он является первым в подмножестве, то мы не можем заменить его на другой.
5 used


3 Шаг:

1 2 3 4
^ Дополнили первое подмножество элементом 4
5 used


4 Шаг:

1 2 3 4
5 Дописали лексикографически минимальный хвост
used

См.также[править]

Источники информации[править]