Правильные скобочные последовательности

Определения

Определение:

Скобочная последовательность — класс комбинаторных объектов, представляющих собой последовательность скобочных символов.

Примеры скобочных последовательностей

• [math](())))([/math]
• [math])()()))()(()())[/math]

Определение:

Правильная скобочная последовательность — частный случай скобочной последовательности, определяющийся следующими образами: "" (пустая строка) есть правильная скобочная последовательность; пусть [math]S[/math] — правильная скобочная последовательность, тогда [math](S)[/math] есть правильная скобочная последовательность; пусть [math]S1[/math], [math]S2[/math] — правильные скобочные последовательности, тогда [math]S1S2[/math] есть правильная скобочная последовательность;

Примеры правильных скобочный последовательностей

• [math]((()()()()))[/math]
• [math](())(()())[/math]

Алгоритм проверки правильности скобочной последовательности

Пусть нам дана скобочная последовательность, записанная в строку [math]s[/math]. Возьмем переменную [math]counter[/math], [math]counter = 0[/math]. Будем последовательно перебирать все символы этой строки. Если мы встречаем открывающуюся скобку, то увеличиваем [math]counter[/math] на [math]1[/math], закрывающую — уменьшаем на [math]1[/math]. Если на протяжении всего перебора [math]counter[/math] было неотрицательным и после завершения осталось нулем, то скобочная последовательность правильна.

Псевдокод

 boolean check(s: string)
   counter = 0
   for i = 1 to length(s)
     if s[i] == '('
       counter++
     else
       counter-- 
     if counter < 0
       return false 
   if counter == 0
     return true
   else
     return false

Надо отметить, что скобочные последовательности могут состоять не только из одного типа скобок. При этом недопустимо такое расположение, когда один тип скобок закрывает другой:

Примеры скобочных последовательностей с несколькими типами скобок

• [math]()[()]\{()()[]\}[/math] — верно
• [math][(]\{\})[/math] — неверно

В этом случае для проверки надо будет использовать стек.

Лексикографический порядок порядок правильных скобочных последовательностей

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

Примеры лексикографического порядка для [math]n[/math] и [math]k[/math], где [math]n[/math] — число открывающихся скобок, а [math]k[/math] — число видов скобок

[math]n = 3[/math]		[math]k = 1[/math]
[math]((()))[/math]	[math](()())[/math]	[math](())()[/math]	[math]()(())[/math]	[math]()()()[/math]

[math]n = 2[/math]		[math]k = 2[/math]
[math]()[][/math]	[math]([])[/math]	[math][()][/math]	[math][]()[/math]

Алгоритм генерации лексикографического порядка будет предложен ниже.

Алгоритмы генерации

Рекурсивный алгоритм получения лексикографического порядка

Пусть нам известно число [math]n[/math]. Надо вывести все правильные скобочные последовательности в лексикографическом порядке с [math]n[/math] открывающимися скобками:

Для запуска алгоритма необходимо сделать вызов gen(n, 0, 0, "").

• ans - строка, в которой мы считаем ответ
• counter_open - количество открывающих скобок в данный момент
• counter_open - количество закрывающих скобок в данный момент

 function gen (n: int, counter_open: int, counter_close: int, ans: string)
   if counter_open + counter_close == 2 [math]\cdot[/math] n
     print ans
     return
   if counter_open < n
     gen(n, counter_open + 1, counter_close, ans + "(")
   if counter_open > counter_close
     gen(n, counter_open, counter_close + 1, ans + ")")

Генерация следующей скобочной последовательности

Пусть нам известна строка [math]s[/math], представляющая собой правильную скобочную последовательность. Нам необходимо вывести следующую скобочную последовательность, а если ее нет, то вывести "No solution". Чтобы получить следующую скобочную последовательность надо найти последнюю открывающуюся скобку, которую можно заменить, заменить ее на закрывающуюся, а оставшиеся в конце скобки (если они есть) заменить на минимально возможную последовательность скобок:

 string next(s: string)
   counter_close = 0
   counter_open = 0
   for i = length(s) downto 1
       if s[i] == '('
         counter_open++ 
         if counter_close > counter_open
           break
       else 
         counter_close++
   delete(s, length(s) - counter_open - counter_close + 1, counter_close + l)
   if s == ""
     return "No Solution"
   else
    s = s +')'
    for j = 1 to counter_open
      s = s + '('
    for j = 1 to counter_close - 1
      s = s + ')'
    return s;

Если эта функция после выполнения выводит [math]true[/math], тогда надо напечатать полученную строку [math]s[/math], если [math]false[/math], то следует вывести "No solution".

Получение лексикографического порядка

Пусть нам известно число [math]n[/math]. Надо вывести все правильные скобочные последовательности в лексикографическом порядке с [math]n[/math] открывающимися скобками:

 function order(n: int)
   s = ""
   for j = 1 to n
     s = s + '('
   for j = 1 to n
     s = s + ')'
   print s
   t = next(s)
   while t != false
     print s
     t = next(s)
   return

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

Получение номера последовательности

Пусть [math]n[/math] — количество пар скобок в последовательности. Требуется по заданной правильной скобочной последовательности найти её номер в списке лексикографически упорядоченных правильных скобочных последовательностей.

Научимся считать вспомогательную динамику [math]d[i][j][/math], где [math]i[/math] — длина скобочной последовательности (она "полуправильная": всякой закрывающей скобке соответствует парная открывающая, но не все открытые скобки закрыты), [math]j[/math] — баланс (т.е. разность между количеством открывающих и закрывающих скобок), [math]d[i][j][/math] — количество таких последовательностей. При подсчёте этой динамики мы считаем, что скобки бывают только одного типа.

Считать эту динамику можно следующим образом. Пусть [math]d[i][j][/math] — величина, которую мы хотим посчитать. Если [math]i = 0[/math], то ответ понятен сразу: [math]d[0][0] = 1[/math], все остальные [math]d[0][j] = 0[/math]. Пусть теперь [math]i \gt 0[/math], тогда переберём, чему мог быть равен последний символ этой последовательности. Если он был равен '(', то до этого символа мы находились в состоянии [math](i-1,j-1)[/math]. Если он был равен ')', то предыдущим было состояние [math](i-1,j+1)[/math]. Таким образом, получаем формулу:

[math]d[i][j] = d[i-1][j-1] + d[i-1][j+1][/math]

(считается, что все значения [math]d[i][j][/math] при отрицательном [math]j[/math] равны нулю). Таким образом, эту динамику мы можем посчитать за [math]O(n^2)[/math].

Перейдём теперь к решению самой задачи. Сначала пусть допустимы только скобки одного типа:

 int getNumber(s: string)
    num = 0
    depth = 0
    for i = 0 to 2 * n - 1
      if s[i] == '('
        depth++
      else
        num += d[2 * n - i - 1][depth + 1]
        depth--
    return num

Пусть теперь разрешены скобки [math]k[/math] типов. Тогда при рассмотрении текущего символа [math]s[i][/math] до пересчёта [math]\rm depth[/math] мы должны перебирать все скобки, которые меньше текущего символа в установленном ранее порядке, пробовать ставить эту скобку в текущую позицию (получая тем самым новый баланс [math]\rm ndepth = \rm depth \pm 1[/math]), и прибавлять к ответу количество соответствующих "хвостов" - завершений (которые имеют длину [math]2n - i - 1[/math], баланс [math]\rm ndepth[/math] и [math]k[/math] типов скобок). Утверждается, что формула для этого количества имеет вид:

[math]d[2n - i - 1][ndepth] \cdot k^{(2n - i - 1 - ndepth) / 2}[/math]

Эта формула выводится из следующих соображений. Сначала мы "забываем" про то, что скобки бывают нескольких типов, и просто берём ответ из [math]d[2n - i - 1][{\rm ndepth}][/math]. Теперь посчитаем, как изменится ответ из-за наличия [math]k[/math] типов скобок. У нас имеется [math]2n - i - 1[/math] неопределённых позиций, из которых [math]\rm ndepth[/math] являются скобками, закрывающими какие-то из открытых ранее, — значит, тип таких скобок мы варьировать не можем. А вот все остальные скобки (а их будет [math](2n - i - 1 - {\rm ndepth}) / 2[/math] пар) могут быть любого из [math]k[/math] типов, поэтому ответ умножается на эту степень числа [math]k[/math].

Сложность данного алгоритма [math]O(n^2 + n \cdot k)[/math].

Получение k-й последовательности

Пусть [math]n[/math] — количество пар скобок в последовательности. В данной задаче по заданному [math]k[/math] требуется найти [math]k[/math]-ую правильную скобочную последовательность в списке лексикографически упорядоченных последовательностей.

Как и в предыдущем разделе, посчитаем динамику [math]d[i][j][/math] — количество правильных скобочных последовательностей длины [math]i[/math] с балансом [math]j[/math].

Пусть сначала допустимы только скобки одного типа:

 string getSequence(n: int, k: int)
   depth = 0
   s = ""
   for i = 0 to 2 * n - 1
     if d[2 * n - (i + 1)][depth + 1] \geqslant k
       s += '('
       depth++
     else
       k -= d[2 * n - (i + 1)][depth + 1]
       s += ')'
       depth--
   return s

Пусть теперь разрешён не один, а [math]k[/math] типов скобок. Тогда алгоритм решения будет отличаться от предыдущего случая только тем, что мы должны домножать значение [math]d[i + 1][\rm ndepth][/math] на величину [math]k^{(2n - i - 1 - \rm ndepth) / 2}[/math], чтобы учесть, что в этом остатке могли быть скобки различных типов, а парных скобок в этом остатке будет только [math]2n - i - 1 - \rm ndepth[/math], поскольку [math]\rm ndepth[/math] скобок являются закрывающими для открывающих скобок, находящихся вне этого остатка (а потому их типы мы варьировать не можем).

Сложность данного алгоритма [math]O(n^2 + n \cdot k)[/math].

Количество правильных скобочных последовательностей. Числа Каталана

Количество правильных скобочных последовательностей со скобками одного типа совпадает с числами Каталана.

Определение:

Числа Каталана (англ. Catalan number) — последовательность чисел, выражающих: количество не изоморфных упорядоченных бинарных деревьев с корнем и [math]n + 1[/math] листьями; количество способов соединения [math]2n[/math] точек на окружности [math]n[/math] не пересекающимися хордами; количество разбиений выпуклого [math](n + 2)[/math]-угольника на треугольники не пересекающимися диагоналями; количество способов полностью разделить скобками [math]n + 1[/math] множитель; количество корректных скобочных последовательностей, состоящих из [math]n[/math] открывающих и [math]n[/math] закрывающих скобок;

Рекурентная формула

[math]C_n = \sum_{i = 0}^{n - 1} C_i C_{n - 1 - i}[/math]

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

Самой левой открывающей скобке [math]l[/math] соответствует определённая закрывающая скобка [math]r[/math], которая разбивает формулу на две части, каждая из которых в свою очередь является правильной скобочной последовательностью. Поэтому, если мы обозначим [math]i = r - l - 1[/math], то для любого фиксированного [math]r[/math] будет ровно [math]C_i C_{n-1-i}[/math] способов. Суммируя это по всем допустимым i, мы и получаем рекуррентную зависимость на [math]C_n[/math].

Аналитическая формула

[math] C_n = \frac{1}{n+1} C_{2n}^{n} [/math]

(здесь через [math]C_n^k[/math] обозначен, как обычно, биномиальный коэффициент).

Эту формулу проще всего вывести из задачи о монотонных путях. Общее количество монотонных путей в решётке размером [math]n \times n[/math] равно [math]C_{2n}^{n}[/math]. Теперь посчитаем количество монотонных путей, пересекающих диагональ. Рассмотрим какой-либо из таких путей, и найдём первое ребро, которое стоит выше диагонали. Отразим относительно диагонали весь путь, идущий после этого ребра. В результате получим монотонный путь в решётке [math](n - 1) \times (n + 1)[/math]. Но, с другой стороны, любой монотонный путь в решётке [math](n - 1) \times (n + 1)[/math] обязательно пересекает диагональ, следовательно, он получен как раз таким способом из какого-либо (причём единственного) монотонного пути, пересекающего диагональ, в решётке [math]n \times n[/math]. Монотонных путей в решётке [math](n - 1) \times (n + 1)[/math] имеется [math]C_{2n}^{n-1}[/math]. В результате получаем формулу:

[math] C_n = C_{2n}^{n} - C_{2n}^{n-1} = \frac{1}{n+1} C_{2n}^{n}[/math]

См. также

• Комбинаторные объекты
• Лексикографический порядок
• Генерация комбинаторных объектов в лексикографическом порядке
• Получение номера по объекту
• Получение объекта по номеру
• Получение следующего объекта

Источники

• Скобочные последовательности, Материал из Википедии — свободной энциклопедии

• Правильная скобочная последовательность, Материал из Википедии — свободной энциклопедии

• MAXimal :: algo :: Правильные скобочные последовательности

• MAXimal :: algo :: Числа Каталана