20071202 efficientalgorithms kulikov_lecture11

с/к “Эффективные алгоритмы”
Лекция 11: Задача выполнимости

А. Куликов

Computer Science клуб при ПОМИ
http://logic.pdmi.ras.ru/∼infclub/

А. Куликов (CS клуб при ПОМИ) 11. Задача выполнимости 1 / 30

План лекции

1 Определение задачи




2 Сведения
Японские кроссворды
Eternity
Максимальный разрез




2 Сведения
Eternity

3 Что мы узнали за сегодня




2 Сведения
Eternity



Формула в КНФ
Определение


Пропозициональной или Булевой (propositional, Boolean)
переменной называется переменная, принимающая значения true
(1) и false (0).


(1) и false (0).
Литералом (literal) называется Булева переменная x или ее
отрицание ¬x.


(1) и false (0).
Клозом (clause) называется дизъюнкция конечного множества
литералов, не содержащего одновременно переменной и ее
отрицания.


(1) и false (0).
отрицания.
k-клозом (k-clause) называется клоз, содержащий ровно k
литералов.


(1) и false (0).
отрицания.
k-клозом (k-clause) называется клоз, содержащий ровно k
литералов.
Формулой в конъюнктивной нормальной форме (КНФ) (formula in
conjunctive normal form, CNF) называется конъюнкция конечного
множества клозов.


Меры сложности формул



n(F ), N(F ) — кол-во различных переменных в F .



m(F ), K (F ) — кол-во клозов в F .



m(F ), K (F ) — кол-во клозов в F .
l (F ), L(F ) — кол-во литералов в (длина) F .


Задача выполнимости




Задача пропозициональной выполнимости (Boolean satisfiability
problem, SAT): определить, выполнима ли данная формула в
КНФ, то есть существует ли набор Булевых значений переменным
формулы, выполняющий формулу. Такой набор называют
выполняющим (satisfying assignment), а формулу, для которой
такой набор существует, — выполнимой (satisfiable).



Задача максимальной выполнимости (maximum satisfiability
problem, SAT): по данной формуле определить, какое
максимальное количество ее клозов может быть выполнено.



Задача максимальной выполнимости (maximum satisfiability
problem, SAT): по данной формуле определить, какое
максимальное количество ее клозов может быть выполнено.
k-SAT, MAX-k-SAT — частные случаи соответствующих задач,
когда все клозы входной формулы содержат не более k литералов.


Пример

Пример


Пример

Пример
F1 = (x ∨ y ∨ ¬z) ∧ (¬x) ∧ (¬y ∨ z)


Пример

Пример
F1 = (x ∨ y ∨ ¬z) ∧ (¬x) ∧ (¬y ∨ z)
n(F1 ) = 3, m(F1 ) = 3, l (F1 ) = 6


Пример

Пример
F1 = (x ∨ y ∨ ¬z) ∧ (¬x) ∧ (¬y ∨ z)
n(F1 ) = 3, m(F1 ) = 3, l (F1 ) = 6
F1 выполнима: x = 0, y = 1, z = 1


Пример

Пример
F1 = (x ∨ y ∨ ¬z) ∧ (¬x) ∧ (¬y ∨ z)
n(F1 ) = 3, m(F1 ) = 3, l (F1 ) = 6
F1 выполнима: x = 0, y = 1, z = 1
F2 = (x ∨ y ) ∧ (x ∨ ¬y ) ∧ (¬x ∨ y ) ∧ (¬x ∨ ¬y )


Пример

Пример
F1 = (x ∨ y ∨ ¬z) ∧ (¬x) ∧ (¬y ∨ z)
n(F1 ) = 3, m(F1 ) = 3, l (F1 ) = 6
F1 выполнима: x = 0, y = 1, z = 1
F2 = (x ∨ y ) ∧ (x ∨ ¬y ) ∧ (¬x ∨ y ) ∧ (¬x ∨ ¬y )
n(F2 ) = 2, m(F2 ) = 4, l (F1 ) = 8


Пример

Пример
F1 = (x ∨ y ∨ ¬z) ∧ (¬x) ∧ (¬y ∨ z)
n(F1 ) = 3, m(F1 ) = 3, l (F1 ) = 6
F1 выполнима: x = 0, y = 1, z = 1
F2 = (x ∨ y ) ∧ (x ∨ ¬y ) ∧ (¬x ∨ y ) ∧ (¬x ∨ ¬y )
n(F2 ) = 2, m(F2 ) = 4, l (F1 ) = 8
F2 невыполнима, но три клоза можно выполнить


NP-трудность



Первая известная NP-полная задача (Кук, 1971).



3-SAT тоже NP-полна.



2-SAT может быть решена за линейное время.



2-SAT может быть решена за линейное время.
MAX-2-SAT NP-трудна, даже если каждая переменная
встречается в формуле не более трех раз.


Важность задачи



http://www.satisfiability.org/ — The International
Conferences on Theory and Applications of Satisfiability Testing.



http://www.satcompetition.org/ — The international SAT
Competitions web page.



http://www.isa.ewi.tudelft.nl/Jsat/ — Journal on
Satisfiability, Boolean Modeling and Computation.



http://www.satlib.org/ — The Satisfiability Library.



http://www.satlive.org/ — Up-to-date links for the Satisfiability
Problem.



http://www.satlive.org/ — Up-to-date links for the Satisfiability
Problem.
http://www.qbflib.org/ — The Quantified Boolean Formulas
Satisfiability Library.




2 Сведения
Eternity



Сведения


Сведения

Многие известные задачи из NP очень просто сводятся к SAT или
MAX-SAT.


Сведения

MAX-SAT.
Сведя задачу к SAT, на практике можно воспользоваться
SAT-солвером.


Сведения

MAX-SAT.
Такой подход иногда помогает, иногда — нет.




2 Сведения
Eternity




Решение японского кроссворда (japanese puzzle) заключается в
восстановлении картинки по длинам блоков подряд идущих
закрашенных клеток в строках и столбцах.


Сведение к SAT



дано поле n × m



достаточно научиться кодировать тот факт, что в i-й строке
закрашены блоки длины a1 , . . . , at



введем такие переменные:



xij — клетка (i, j) закрашена



xij — клетка (i, j) закрашена
yijp — p-й блок строки i начинается в клетке j


Сведение к SAT (продолжение)
итак, дана строка i с блоками a1 , . . . , at


p-й блок начинается ровно в одной клетке:
(yi1p ∨ yi2p ∨ · · · ∨ yimp ) ∧ {(¬yijp ∨ ¬yikp )}j̸=k



(p + 1)-й блок начинается позже конца p-го:
{(¬yij(p+1) ∨ ¬yi(j+k)p )}k≥0



{(¬yij(p+1) ∨ ¬yi(j+k)p )}k≥0

если p-й блок начинается в клетке j, то соответствующие ap
клеток закрашены:
{(¬yijp ∨ xi(j+k) )}0≤k<ap



{(¬yij(p+1) ∨ ¬yi(j+k)p )}k≥0

если p-й блок начинается в клетке j, то соответствующие ap
клеток закрашены:
{(¬yijp ∨ xi(j+k) )}0≤k<ap

если клетка не принадлежит ни одному из блоков, то она не
закрашена:
(yi(j−a1 +1)1 ∨ · · · ∨ yij1 ∨ · · · ∨ yi(j−at +1)p ∨ · · · ∨ yijt ∨ ¬xij )

Замечание


Замечание
если какой-то индекс “вылезает”, то соответствующим
переменным просто присваиваем значение 0


Замечание
если какой-то индекс “вылезает”, то соответствующим
переменным просто присваиваем значение 0
можно сводить более эффективно


Примеры решенных кроссвордов




2 Сведения
Eternity



Игра Eternity
Нужно замостить квадрат заданным набором доминошек так, чтобы
узоры на граничащих частях доминошек совпадали.



дан квадрат n × n и n2 доминошек



перенумеруем все клетки и доминошки числами от 1 до n2



заводим два типа переменных:



xij — в i-й клетке стоит j-я доминошка



xij — в i-й клетке стоит j-я доминошка
yjk — j-я доминошка в “положении” k, 1 ≤ k ≤ 4 (например, так:
k = 1 — не повернута, k = 2 — повернута на /2 по часовой,
k = 3 — на , k = 4 — на /2 против часовой)


каждая доминошка находится ровно в одном “положении”:
(yj1 ∨ yj2 ∨ yj3 ∨ yj4 ) ∧ {(¬yjp ∨ ¬yjq )}p̸=q



каждая доминошка стоит хотя бы в одной клетке:
{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ]



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ]

в каждой клетке стоит хотя бы одна доминошка:
{(xi1 ∨ xi2 ∨ · · · ∨ xin2 )}i∈[n2 ]



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ]

в каждой клетке стоит хотя бы одна доминошка:
{(xi1 ∨ xi2 ∨ · · · ∨ xin2 )}i∈[n2 ]

доминошка не стоит в двух клетках одновременно:
{(¬xpj ∨ ¬xqj )}p̸=q



в клетке одновременно не стоят две доминошки:

{(¬xip ∨ ¬xiq )}p̸=q



в клетке одновременно не стоят две доминошки:

{(¬xip ∨ ¬xiq )}p̸=q

две доминошки должны граничить равными ребрами:

{(¬xi1 j1 ∨ ¬xi2 j2 ∨ ¬yi1 k1 ∨ ¬yi2 k2 )},

если клетки i1 и i2 — соседи и при постановке туда доминошек j1
и j2 в “положениях”, соответственно, k1 и k2 нарушается узор


Не все так просто, тем не менее



итак, мы записали игру Eternity в виде задачи выполнимости
конкретной формулы при помощи полиномиального сведения



осталось напустить на полученную формулу какой-нибудь
эффективный SAT-солвер



но в чем же тогда подвох?



но в чем же тогда подвох?
давайте примерно оценим длину полученной формулы


Оценка длины

(yj1 ∨ yj2 ∨ yj3 ∨ yj4 ) ∧ {(¬yjp ∨ ¬yjq )}p̸=q (4 + 2 · 6)n2



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ] n2 · n2



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ] n2 · n2
{(xi1 ∨ xi2 ∨ · · · ∨ xin2 )}i∈[n2 ] n2 · n2



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ] n2 · n2
{(xi1 ∨ xi2 ∨ · · · ∨ xin2 )}i∈[n2 ] n2 · n2
n2 · n
(︀ )︀
{(¬xpj ∨ ¬xqj )}p̸=q 2



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ] n2 · n2
{(xi1 ∨ xi2 ∨ · · · ∨ xin2 )}i∈[n2 ] n2 · n2
n2 · n
(︀ )︀
(︀2)︀
2· n
{(¬xip ∨ ¬xiq )}p̸=q n 2



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ] n2 · n2
{(xi1 ∨ xi2 ∨ · · · ∨ xin2 )}i∈[n2 ] n2 · n2
n2 · n
(︀ )︀
(︀2)︀
2· n
{(¬xip ∨ ¬xiq )}p̸=q n 2
{(¬xi1 j1 ∨ ¬xi2 j2 ∨ ¬yi1 k1 ∨ ¬yi2 k2 )} ∼ n2 · n2 · n2



{(x1j ∨ x2j ∨ · · · ∨ xn2 j )}j∈[n2 ] n2 · n2
{(xi1 ∨ xi2 ∨ · · · ∨ xin2 )}i∈[n2 ] n2 · n2
n2 · n
(︀ )︀
(︀2)︀
2· n
{(¬xip ∨ ¬xiq )}p̸=q n 2
{(¬xi1 j1 ∨ ¬xi2 j2 ∨ ¬yi1 k1 ∨ ¬yi2 k2 )} ∼ n2 · n2 · n2
Итого, порядок длины формулы будет n6 , что при n = 16 составляеят
166 = 224 = 16 777 216.




2 Сведения
Eternity



Задача о максимальном разрезе

Задача о максимальном разрезе (maximum cut, MAX-CUT)
заключается в нахождении такого разбиения вершин графа на две
части, при котором количество ребер, концы которых принадлежат
разным частям, максимально.




Одна из знаменитых 21-й NP-полной задачи Карпа.




Одна из знаменитых 21-й NP-полной задачи Карпа.
Остается NP-полной даже на графах степени 3.


Сведение к MAX-2-SAT



каждой вершине u графа G (V , E ) поставим в соответствие
переменную xu (xu = 1 — вершина u принадлежит первой части)



для каждого ребра (u, v ) запишем два клоза (xu ∨ xv ), (¬xu ∨ ¬xv )



видно, что набор значений переменным u, v выполняет оба клоза,
когда u и v в разных частях, и выполняет только один из них —
когда в одной



видно, что набор значений переменным u, v выполняет оба клоза,
когда u и v в разных частях, и выполняет только один из них —
когда в одной
таким образом, максимальное количество одновременно
выполнимых клозов полученной формулы равно
|E | + |MAX-CUT(G )|




2 Сведения
Eternity



Что мы узнали за сегодня?



MAX-SAT.



MAX-SAT.
Сведения:



MAX-SAT.
Сведения:
японские кроссворды, игра Eternity — к SAT



MAX-SAT.
Сведения:
японские кроссворды, игра Eternity — к SAT
задача о максимальном разрезе — к MAX-SAT


Спасибо за внимание!


20071202 efficientalgorithms kulikov_lecture11

More Related Content

Similar to 20071202 efficientalgorithms kulikov_lecture11 (20)

More from Computer Science Club (20)

20071202 efficientalgorithms kulikov_lecture11