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

Материал из Викиконспекты
Перейти к: навигация, поиск
(Разбиение на слагаемые: исправлен внешний цикл)
(Разбиение на слагаемые: исправлена грамматика и косметические изменения)
Строка 66: Строка 66:
  
 
== Разбиение на слагаемые ==
 
== Разбиение на слагаемые ==
Рассмотрим алгоритм получения номера в лексикографическом порядке данного разбиение на слагаемые числа <tex>N</tex>. Нужно помнить о том, что разбиения, отличающиеся только порядком слагаемых, считаются одинаковыми. Из всех разбиений, получаемых перестановками слагаемых, выберем то, где слагаемые упорядочены лексикографически и будем строить его.   
+
Рассмотрим алгоритм получения номера в лексикографическом порядке данного разбиения на слагаемые числа <tex>N</tex>. Нужно помнить о том, что разбиения, отличающиеся только порядком слагаемых, считаются одинаковыми. Из всех разбиений, получаемых перестановками слагаемых, выберем то, где слагаемые упорядочены лексикографически, и будем строить его.   
  
 
*<tex>\mathtt{numOfPart}</tex> {{---}} искомый номер разбиения
 
*<tex>\mathtt{numOfPart}</tex> {{---}} искомый номер разбиения
Строка 74: Строка 74:
 
*<tex>\mathtt{d[i][j]}</tex> {{---}} количество разбиений числа <tex>i</tex> на слагаемые, где каждое слагаемое <tex>\geqslant j</tex>.   
 
*<tex>\mathtt{d[i][j]}</tex> {{---}} количество разбиений числа <tex>i</tex> на слагаемые, где каждое слагаемое <tex>\geqslant j</tex>.   
  
Пересчитывать <tex>\mathtt{d[i][j]}</tex> будем по возрастанию <tex>i</tex>, а при равенстве <tex>i</tex> по убыванию <tex>j</tex>.  
+
Пересчитывать <tex>\mathtt{d[i][j]}</tex> будем по возрастанию <tex>i</tex>, а при равенстве <tex>i</tex> {{---}} по убыванию <tex>j</tex>.  
  
Разбиение числа, в котором каждое слагаемое <tex> \geqslant j</tex> может либо содержать слагаемое <tex>j</tex>, таких разбиений <tex>\mathtt{d[i - j][j]}</tex>, либо не содержать, таких разбиений <tex>\mathtt{d[i][j + 1]}</tex>.  
+
Разбиение числа, в котором каждое слагаемое <tex> \geqslant j</tex> может либо содержать слагаемое <tex>j</tex> (таких разбиений <tex>\mathtt{d[i - j][j]}</tex>), либо не содержать (таких разбиений <tex>\mathtt{d[i][j + 1]}</tex>).  
  
 
Получаем рекуррентное соотношение для подсчёта <tex>d</tex>:
 
Получаем рекуррентное соотношение для подсчёта <tex>d</tex>:
Строка 87: Строка 87:
  
  
  '''int''' part2num(part: '''list<int>'''):
+
'''int''' part2num(part: '''list<int>'''):
    numOfPart = 0
+
  numOfPart = 0
    '''for''' i = 1 '''to''' K  <font color=green>// <tex>K</tex> {{---}} количество слагаемых в исходном разбиении </font>   
+
  '''for''' i = 1 '''to''' part.size
      '''for''' j = last '''to''' part[i] - 1 <font color=green>// перебираем все элементы, лексикографически меньше нашего, но больше или равны предыдущего</font>     
+
    '''for''' j = last '''to''' part[i] - 1           <font color=green>// перебираем все элементы, лексикографически меньшие нашего, но большие или равные предыдущего</font>     
        numOfPart += d[N - sum - j][j] <font color=green>// прибавляем количество перестановок, которые могли начинаться с <tex>j</tex></font>
+
      numOfPart += d[N - sum - j][j]       <font color=green>// прибавляем количество перестановок, которые могли начинаться с <tex>j</tex></font>
      sum += part[i] <font color=green>// увеличиваем уже поставленную сумму</font>
+
    sum += part[i]                         <font color=green>// увеличиваем уже поставленную сумму</font>
      last = part[i] <font color=green>// обновляем последний поставленный элемент </font>
+
    last = part[i]                         <font color=green>// обновляем последний поставленный элемент </font>
    '''return''' numOfPart <font color=green>// возвращаем ответ</font>
+
  '''return''' numOfPart                         <font color=green>// возвращаем ответ</font>
  
 
Стоит отметить, что количество итераций вложенного цикла не более, чем <tex>N</tex>, так как всего количество возможных слагаемых {{---}} <tex>N</tex>, и ни какое из них цикл не обработает дважды, поскольку каждый раз начинает с <tex>last</tex>, которое больше чем любое из обработанных чисел. Поэтому асимптотика алгоритма {{---}} <tex>O(N)</tex>.
 
Стоит отметить, что количество итераций вложенного цикла не более, чем <tex>N</tex>, так как всего количество возможных слагаемых {{---}} <tex>N</tex>, и ни какое из них цикл не обработает дважды, поскольку каждый раз начинает с <tex>last</tex>, которое больше чем любое из обработанных чисел. Поэтому асимптотика алгоритма {{---}} <tex>O(N)</tex>.

Версия 01:30, 26 декабря 2017

Описание алгоритма

Номер данного комбинаторного объекта равен количеству меньших в лексикографическом порядке комбинаторных объектов (нумерацию ведём с [math]0[/math]). Все объекты меньшие данного можно разбить на непересекающиеся группы по длине совпадающего префикса. Тогда количество меньших объектов можно представить как сумму количеств объектов у которых префикс длины [math]i[/math] совпадает, а [math]i+1[/math] элемент лексикографически меньше [math](i+1)[/math]-го в данном объекте ([math]i = 0..n-1[/math]). Следующий алгоритм вычисляет эту сумму:

  • [math]\mathtt{numOfObject}[/math] — искомый номер комбинаторного объекта,
  • [math]\mathtt{a[1..n]}[/math] — данный комбинаторный обьект, состоящий из числовых представлений лексикографически упорядоченных элементов множества [math]A[/math],
  • [math]\mathtt{d[i][j]}[/math] — количество комбинаторных объектов с префиксом от [math]1[/math] до [math]i-1[/math] равным данному и с [math]i[/math]-м элементом равным [math]j[/math],
int object2num(a: list<A>):
  numOfObject = 0                          
  for i = 1 to n                   // перебираем элементы комбинаторного объекта
    for j = 1 to a[i] - 1          // перебираем элементы, в лексикографическом порядке меньшие рассматриваемого
      if элемент [math]j[/math] можно поставить на [math]i[/math]-e место
        numOfObject += d[i][j]
  return numOfObject

Сложность алгоритма — [math]O(nk) [/math], где [math]k[/math] — количество различных элементов, которые могут находиться в данном комбинаторном объекте. Например, для битового вектора [math]k=2,[/math] поскольку возможны только [math]0[/math] и [math]1[/math]. Количества комбинаторных объектов с заданными префиксами считаются известными, и их подсчет в сложности не учитывается. Приведем примеры способов получения номеров некоторых из комбинаторных объектов по данному объекту.

Битовые вектора

Рассмотрим алгоритм получения номера [math]i[/math] в лексикографическом порядке данного битового вектора размера [math]n[/math]. Всего существует [math]2^n[/math] битовых векторов длины [math]n[/math]. На каждой позиции может стоять один из двух элементов независимо от того, какие элементы находятся в префиксе, поэтому поиск элементов меньше рассматриваемого можно упростить до проверки элемента на равенство [math]1[/math]:

  • [math]\mathtt{bitvector[1..n]}[/math] — данный вектор,
  • [math]\mathtt{numOfBitvector}[/math] — искомый номер вектора,
int bitvector2num(bitvector: list<int>):
  numOfBitvector = 0
  for i = 1 to n                                        
    if bitvector[i] == 1  
      numOfBitvector += [math]2^{n-i}[/math]
  return numOfBitvector

Асимптотика алгоритма — [math]O(n) [/math].

Перестановки

Рассмотрим алгоритм получения номера в лексикографическом порядке по данной перестановке размера [math]n[/math],

  • [math]\mathtt{a[1..n]}[/math] — данная перестановка,
  • [math]\mathtt{P[1..n]}[/math] — количество перестановок данного размера,
  • [math]\mathtt{was[1..n]}[/math] — использовали ли мы уже эту цифру в перестановке,
int permutation2num(a: list<int>):
  numOfPermutation = 0
  for i = 1 to n                        // [math]n[/math] — количество элементов в перестановке 
    for j = 1  to a[i] - 1              // перебираем элементы, лексикографически меньшие нашего, которые могут стоять на [math]i[/math]-м месте  
      if was[j] == false                // если элемент [math]j[/math] ранее не был использован 
        numOfPermutation += P[n - i]    // все перестановки с префиксом длиной [math]i-1[/math] равным нашему, и [math]i[/math]-й элемент у которых   
                                           меньше нашего в лексикографическом порядке, идут раньше данной перестановки                
    was[a[i]] = true                    // [math]i[/math]-й элемент использован            
  return numOfPermutation

Асимптотика алгоритма — [math]O(n ^ 2) [/math] и [math]O(n) [/math] для предподсчёта.

Сочетания

Рассмотрим алгоритм получения номера в лексикографическом порядке данного сочетания из [math]n[/math] по [math]k[/math]. Как известно, количество сочетаний из [math]n[/math] по [math]k[/math] обозначается как [math]\binom{n}{k}[/math]. Тогда число сочетаний, в которых на позиции [math]1[/math] стоит значение [math]val_1[/math], равно [math]$$\sum\limits^{val_1-1}_{i=1} {\binom{n-i}{k-1}}$$[/math]; число сочетаний, в которых на позиции [math]2[/math] стоит значение [math]val_2[/math], равно [math]$$\sum\limits^{val_2-1}_{i=val_1+1} {\binom{n-i}{k-2}}$$[/math]. Аналогично продолжаем по следующим позициям:

  • [math]\mathtt{numOfChoose}[/math] — искомый номер сочетания,
  • [math]\mathtt{C[n][k]}[/math] — количество сочетаний из [math]n[/math] по [math]k[/math], [math]\mathtt{C[n][0] = 1}[/math],
  • [math]\mathtt{choose[1..K]}[/math] — данное сочетание, состоящее из [math]K[/math] чисел от [math]1[/math] до [math]N[/math], из технических соображений припишем ноль в начало сочетания: [math]\mathtt{choose[0] = 0}[/math],
int choose2num(choose: list<int>):
  numOfChoose = 0
  for i = 1 to K                                       
    for  j = choose[i - 1] + 1 to choose[i] - 1
      numOfChoose += C[N - j][K - i]
  return numOfChoose

Асимптотика алгоритма — [math]O(K \cdot N) [/math] и [math]O(K \cdot N) [/math] для предподсчёта.

Разбиение на слагаемые

Рассмотрим алгоритм получения номера в лексикографическом порядке данного разбиения на слагаемые числа [math]N[/math]. Нужно помнить о том, что разбиения, отличающиеся только порядком слагаемых, считаются одинаковыми. Из всех разбиений, получаемых перестановками слагаемых, выберем то, где слагаемые упорядочены лексикографически, и будем строить его.

  • [math]\mathtt{numOfPart}[/math] — искомый номер разбиения
  • [math]\mathtt{last}[/math] — последнее поставленное число в разбиении.
  • [math]\mathtt{sum}[/math] — сумма, которую мы уже поставили.
  • [math]\mathtt{part[1 \ldots N]}[/math] — данное разбиение
  • [math]\mathtt{d[i][j]}[/math] — количество разбиений числа [math]i[/math] на слагаемые, где каждое слагаемое [math]\geqslant j[/math].

Пересчитывать [math]\mathtt{d[i][j]}[/math] будем по возрастанию [math]i[/math], а при равенстве [math]i[/math] — по убыванию [math]j[/math].

Разбиение числа, в котором каждое слагаемое [math] \geqslant j[/math] может либо содержать слагаемое [math]j[/math] (таких разбиений [math]\mathtt{d[i - j][j]}[/math]), либо не содержать (таких разбиений [math]\mathtt{d[i][j + 1]}[/math]).

Получаем рекуррентное соотношение для подсчёта [math]d[/math]:

[math]d[i][j] = \left \{\begin{array}{ll} 1, & i = j, \\ 0, & i \lt j \\ d[i][j + 1] + d[i - j][j], & i \gt j \end{array} \right. [/math]


int part2num(part: list<int>):
  numOfPart = 0
  for i = 1 to part.size
    for j = last to part[i] - 1            // перебираем все элементы, лексикографически меньшие нашего, но большие или равные предыдущего    
      numOfPart += d[N - sum - j][j]       // прибавляем количество перестановок, которые могли начинаться с [math]j[/math]
    sum += part[i]                         // увеличиваем уже поставленную сумму
    last = part[i]                         // обновляем последний поставленный элемент 
  return numOfPart                         // возвращаем ответ

Стоит отметить, что количество итераций вложенного цикла не более, чем [math]N[/math], так как всего количество возможных слагаемых — [math]N[/math], и ни какое из них цикл не обработает дважды, поскольку каждый раз начинает с [math]last[/math], которое больше чем любое из обработанных чисел. Поэтому асимптотика алгоритма — [math]O(N)[/math].

Асимптотика алгоритма — [math] O (N)[/math] и [math]O(N^2)[/math] на предподсчёт.

См. также

Источники информации

  • Программирование в алгоритмах / С. М. Окулов. — М.: БИНОМ. Лаборатория знаний, 2002. стр.31
  • Дискретная математика. Теория и практика решения задач по информатике / С. М. Окулов. — М.: БИНОМ. Лаборатория знаний, 2008.