Обратное преобразование Барроуза-Уиллера — различия между версиями
| Строка 359: | Строка 359: | ||
===Сложность=== | ===Сложность=== | ||
Данный алгоритм работает за <tex>0(NlogN)</tex> действий и требует <tex>0(N)</tex> памяти. Однако, если размер алфавита не очень большой, то для выяснения 1ого столбца матрицы можно использовать сортировку подсчетом, в этом случае алгоритм работает за <tex>0(N+M)</tex> действий и требует <tex>0(N+M)</tex> памяти, где <tex>M</tex> - размер алфавита. | Данный алгоритм работает за <tex>0(NlogN)</tex> действий и требует <tex>0(N)</tex> памяти. Однако, если размер алфавита не очень большой, то для выяснения 1ого столбца матрицы можно использовать сортировку подсчетом, в этом случае алгоритм работает за <tex>0(N+M)</tex> действий и требует <tex>0(N+M)</tex> памяти, где <tex>M</tex> - размер алфавита. | ||
| + | ===Псевдокод=== | ||
| + | Пусть N - количество символов во входной строке, M - количество символов в алфавите, k - номер исходной строки в матрице перестановок, s - входящая строка, count - массив для сортировки подсчетом, t - вектор обратного преобразования. | ||
| + | |||
| + | '''//считаем частоты символов''' | ||
| + | '''for i = 0..N count[i] = 0''' | ||
| + | '''for i = 0..N count[s[i]]++''' | ||
| + | '''//упорядочиваем символы, чтобы получить первый столбец исходной матрицы''' | ||
| + | '''//count[i] указывает на первую позицию символа i в первом столбце''' | ||
| + | '''sum = 0''' | ||
| + | '''for i = 0..N''' | ||
| + | ''' sum = sum + count[i]''' | ||
| + | ''' count[i] = sum - count[i]''' | ||
| + | '''//создаем вектор обратного преобразования''' | ||
| + | '''for i = 0..N''' | ||
| + | ''' t[count[s[i]]]] = i''' | ||
| + | ''' count[s[i]]++''' | ||
| + | '''//восстанавливаем исходный текст''' | ||
| + | '''j = t[x]''' | ||
| + | '''for i = 0..N''' | ||
| + | ''' print(s[j])''' | ||
| + | ''' j = t[j]''' | ||
Версия 09:30, 2 января 2011
Содержание
Наивный алгоритм
Описание
Если отсортировать символы последнего столбца матрицы перестановок - а перед началом обратного преобразования нам известен только он - то мы получим первый столбец данной матрицы,так как строки матрицы были отсортированы в лексикографическом порядке. Известно, что символы последнего столбца предшествуют символам первого, находящимся в той же строке, поскольку использовался циклический сдвиг. Таким образом, мы имеем уже список пар, состоящих из символа последнего столбца и следующего за ним символа первого столбца. После сортировки этих пар мы обретаем знание о двух первых столбцах матрицы. И так далее, пока мы не получим всю матрицу перестановок, после чего выписываем строку с заданным номером. Рассмотрим этот алгоритм на конкретном примере. Пусть дана строка "абракадабра". В результате преобразования получаем набор ("рдакраааабб", 3).
|
|
|
|
|
|
|
|
|
|
Зная номер исходной строки - 3, мы воспроизводим входные данные - "абракадабра".
Сложность
Как несложно посчитать сложность данного алгоритма , также он требует памяти.
Оптимизация
Описание
Однако, данный алгоритм можно оптимизировать. Заметим, что при каждом проявлении доселе неизвестного столбца выполнялись одни и те же действия. А именно, из строки, начинающейся с некоторого символа последнего столбца получалась строка, в которой этот символ находится на первой позиции. Теперь если определить порядок получения символов 1ого столбца из символов последнего, то полученные значения будут являться вектором обратного преобразования. Для получения исходной строки надо выписать символы из преобразованной строки в порядке, соответствующему данному вектору, начиная с номера, содержащегося в исходных данных. Рассмотрим данный алгоритм на том же примере.
|
|
Выписывая элементы исходной строки в порядке полученного вектора, начиная с 3-его номера, получаем :
| 6 | => | 10 | => | 4 | => | 8 | => | 3 | => | 7 | => | 1 | => | 5 | => | 9 | => | 0 | => | 2 |
| а | б | р | а | к | а | д | а | б | р | а |
Сложность
Данный алгоритм работает за действий и требует памяти. Однако, если размер алфавита не очень большой, то для выяснения 1ого столбца матрицы можно использовать сортировку подсчетом, в этом случае алгоритм работает за действий и требует памяти, где - размер алфавита.
Псевдокод
Пусть N - количество символов во входной строке, M - количество символов в алфавите, k - номер исходной строки в матрице перестановок, s - входящая строка, count - массив для сортировки подсчетом, t - вектор обратного преобразования.
//считаем частоты символов for i = 0..N count[i] = 0 for i = 0..N count[s[i]]++ //упорядочиваем символы, чтобы получить первый столбец исходной матрицы //count[i] указывает на первую позицию символа i в первом столбце sum = 0 for i = 0..N sum = sum + count[i] count[i] = sum - count[i] //создаем вектор обратного преобразования for i = 0..N t[count[s[i]]]] = i count[s[i]]++ //восстанавливаем исходный текст j = t[x] for i = 0..N print(s[j]) j = t[j]
