Изменения

Перейти к: навигация, поиск

Динамическое программирование по профилю

7434 байта добавлено, 23:00, 18 декабря 2012
Новая страница: «'''Динамическое программирование по профилю''' <tex>-</tex> способ оптимизации перебора количе...»
'''Динамическое программирование по профилю''' <tex>-</tex> способ оптимизации перебора количества вариантов с помощью динамического программирования, когда одно из измерений не большое.
== '''Задача о замощении домино''' ==

Найти количество способов замощении таблицы <tex>n\times m</tex> с помощью домино размерами <tex>1\times 2,2\times 1</tex>.

'''Решение:'''
В качестве состояния динамики будем использовать профиля размерами n. В этом профиле 1 будет означать, что домино лежит горизонтально и заканчивается на этом столбце, иначе 0. Таких профилей будет <tex>2^n</tex>.
Теперь проверим из какого профиля в какой можно перейти.

Из профиля i в профиль j можно перейти если выполняются условия:

* Можно положить горизонтальные домино. То есть там где в j профиле стоит 1, в i профиле должен стоять 0

* Можно доложить в оставшиеся клетки вертикальные домино. То есть оставшиеся 0 в i профиле должны образовывать четные подстроки.

Пусть <tex>d[i][j] = 1</tex> если из профиля i можно перейти в j-ый, иначе 0.

Пусть так же <tex>a[k][i]</tex> - количество способов замощения первые k-1 столбцов и заканчивавшийся на i профиле.
Тогда <tex>a[k][i]=\displaystyle \sum_{j=0}^{2^n -1} a[k-1][j]\cdot d[j][i]</tex>

Ответом будет <tex> \sum a[m][i]</tex>, где i : профиль, который может быть последним (т.е. все группы из 0 имеют четные размеры)

Для удобства можно хранить профиля в виде двоичных масок

'''Реализация'''

//n, m размеры таблицы
for i = 0..(1<<n)-1
for j = 0..(1<<n)-1
if можно перейти из i в j профиль
d[i][j] = 1;
else
d[i][j] = 0;
a[0][0] = 1; //Так как мы можем начать только с профиля где все клетки 0
for k = 1..m-1
for i = 0..(1<<n)-1
for j = 0..(1<<n)-1
a[k][i] += a[k-1][j] * d[j][i];
ans = 0;
for i = 0..(1<<n)-1
if можно закончить i профилем
ans += a[m-1][i];
return ans;

''' Оценка сложности : '''
подсчет <tex>d - 2^{2n}</tex> , и подсчет <tex>a - 2^{2n}\times m</tex> в итоге <tex>O(2^{2n}\times m)</tex>

''' Оценка памяти : '''
<tex>O(2^{2n}+2^{2n}\times m)</tex>, так же можно заметить что в массиве <tex>a</tex> для k состояния нам нужно только k-1 состояние, при такой реализации нужно будет <tex>O(2^{2n})</tex>. Еще можно не считать массив d, а просто каждый раз перепроверять можем ли мы перейти в это состояние в итоге потребуется <tex>O(2\times 2^n)</tex> памяти, но нам потребуется больше времени <tex>2^{2n}\times m\times f(i,j)</tex>, где <tex>f(i,j)</tex> время проверки возможности перехода из i в j, обычно проверка равна n и тогда время получается <tex>O(2^{2n}\times m\times n)</tex>.



== '''Задача о симпатичных узорах''' ==

Дана таблица <tex>n\times m</tex>, каждая клетка которой может быть окрашена в один из двух
цветов: белый или черный. Симпатичным узором называется такая раскраска, при
которой не существует квадрата <tex>2\times 2</tex>, в котором все клетки одного цвета. Требуется найти количество симпатичных узоров для соответствующей таблицы.

'''Решение:'''
В качестве состояния динамики будем использовать профиль размеров n. В этом профиле 1 будет означать что клетка закрашена в черный цвет, и 0 если в белый.
Из профиля i в j-ый можно перейти если выполнено условие:
* если поставить i и j профиль рядом, то не должно быть квадратов <tex>2\times 2</tex> одного цвета

Пусть <tex>d[i][j] = 1</tex> если из профиля i можно перейти в j-ый, иначе 0.

Пусть так же <tex>a[k][i]</tex> - количество способов раскрашивания первые k-1 столбцов и заканчивавшийся на i профиле.
Тогда <tex>a[k][i]=\displaystyle \sum_{j=0}^{2^n -1} a[k-1][j]\cdot d[j][i]</tex>

Ответом будет <tex> \displaystyle \sum_{j=0}^{2^n -1} a[m][i]</tex>

Для удобства можно хранить профиля в виде двоичных масок

'''Реализация:'''
//n, m размеры таблицы
for i = 0..(1<<n)-1
for j = 0..(1<<n)-1
if можно перейти из i в j профиль
d[i][j] = 1;
else
d[i][j] = 0;
for i = 0..(1<<n)-1
a[i][0] = 1; //Так как мы можем начать c любого профиля
for k = 1..m-1
for i = 0..(1<<n)-1
for j = 0..(1<<n)-1
a[k][i] += a[k-1][j] * d[j][i];
ans = 0;
for i = 0..(1<<n)-1
ans += a[m-1][i]//Так как мы можем закончить любым профилем
return ans;

''' Оценка сложности : '''
подсчет <tex>d - 2^{2n}</tex> , и подсчет <tex>a - 2^{2n}\times m</tex> в итоге <tex>O(2^{2n}\times m)</tex>

''' Оценка памяти : '''
<tex>O(2^{2n}+2^{2n}\times m)</tex>, так же можно заметить что в массиве <tex>a</tex> для k состояния нам нужно только k-1 состояние, при такой реализации нужно будет <tex>O(2^{2n})</tex>. Еще можно не считать массив d, а просто каждый раз перепроверять можем ли мы перейти в это состояние в итоге потребуется <tex>O(2\times 2^n)</tex> памяти, но нам потребуется больше времени <tex>2^{2n}\times m\times f(i,j)</tex>, где <tex>f(i,j)</tex> время проверки возможности перехода из i в j, обычно проверка равна n и тогда время получается <tex>O(2^{2n}\times m\times n)</tex>.
27
правок

Навигация