『笔记』2-SAT

前置

\(SAT\) 是适定性( \(Satisfiability\) )问题的简称。一般形式为 \(k \ -\) 适定性问题,简称 \(k-SAT\) 。而当 \(k>2\) 时该问题为 \(NP\) 完全的。所以我们只研究 \(k=2\) 的情况。

定义

\(2-SAT\) ,简单的说就是给出 \(n\) 个集合,每个集合有两个元素,已知若干个 \(<a,b>\) ,表示 \(a\)\(b\) 矛盾(其中 \(a\)\(b\) 属于不同的集合)。然后从每个集合选择一个元素,判断能否一共选 \(n\) 个两两不矛盾的元素。显然可能有多种选择方案,一般题中只需要求出一种即可。

有学长者云:

\(n\)\(01\) 变量 \(x_1∼x_n\),另有 \(m\) 个变量取值需要满足的限制。

每个限制是一个 $$ 元组 \((x_{p1},x_{p2},\dots,x_{pk})\) ,满足 \(x_{p1} \oplus x_{p2} \oplus \dots \oplus x_{pk} = a\) 。其中 \(a\)\(0/1\)\(\oplus\) 是某种二元 \(bool\) 运算。

要求构造一种满足所有限制的变量的赋值方案。

\(2-SAT\) 问题是通过建立图论模型,在 \(O(n+m)\) 的时间复杂度内判断是否有解,若有解可以构造出一组合法解。

思路

值得注意的是在不同的题目中二元 \(bool\) 运算可能有差异,但是建图的基本思路大致相同。

来观摩一组 OI-Wiki 的例子:

比如邀请人来吃喜酒,夫妻二人必须去一个,然而某些人之间有矛盾(比如 \(A\) 先生与 \(B\) 女士有矛盾, \(C\) 女士不想和 \(D\) 先生在一起),那么我们要确定能否避免来人之间没有矛盾,有时需要方案。这是一类生活中常见的问题。

使用布尔方程表示上述问题。设 \(a\) 表示 \(A\) 先生去参加,那么 \(B\) 女士就不能参加( \(\lnot a\)); \(b\) 表示 C 女士参加,那么 \(\lnot b\) 也一定成立( \(D\) 先生不参加)。总结一下,即 \((a \lor b)\) (变量 \(a\)\(b\) 至少满足一个)。对这些变量关系建有向图,则有:\(\lnot a \Rightarrow b \land \lnot b \Rightarrow a\)\(a\) 不成立则 \(b\) 一定成立;同理,\(b\) 不成立则 \(a\) 一定成立)。建图之后,我们就可以使用缩点算法来求解 \(2-SAT\) 问题了。

核心

Tarjan 缩点大法

\(Tarjan\) 大法好!

主要还是考虑如何更合适地建图

再来一组例子:

假设有 \(a_1\)\(b_2\)\(a_2\)\(b_1\) 两对,已知 \(a_1\)\(b_2\) 间有矛盾,于是为了方案自洽,由于两者中必须选一个,所以我们就要拉两条有向边 \((a_1,b_1)\)\((b_2,a_2)\) 表示选了 \(a_1\) 则必须选 \(b_1\) ,选了 \(b_2\) 则必须选 \(a_2\) 才能够自洽。

然后通过这样建边再跑一遍 \(Tarjan\) 判断是否有一个集合中的两个元素在同一个强连通分量中,若有则不可能,否则输出方案。构造方案只需要把几个不矛盾的强连通分量拼起来就好了。

  • 输出方案时可以通过变量在图中的拓扑序确定该变量的取值。如果变量 \(\lnot x\) 的拓扑序在 \(x\) 之后,那么取 \(x\) 值为真。应用到 \(Tarjan\) 算法的缩点,即 \(x\) 所在强连通分量编号在 \(\lnot x\) 之前时,取 \(x\) 为真。因为 \(Tarjan\) 算法求强连通分量时使用了栈,所以 \(Tarjan\) 求得的强连通分量编号相当于反拓扑序。

  • 时间复杂度为 \(O(n+m)\)

暴力DFS

\(Tarjan\) 好? \(DFS\) 表示不服。

\(DFS\) 大法妙!

直接选取图上一个点,沿着一条路径搜下去,如果一个点被选择了,那么这条路径以后的点都将被选择。

如果出现一个集合中的两者都被选择了,那么此即为矛盾情况。

例题

P4782 【模板】2-SAT 问题

真·模板题

思路

  • 这是一道模板题。

  • 显而易见每个变量 \(x_i\) 都可以被分开存储,即拆分成 \(i\)\(i+n\) ,分别表示 \(x_i=1\)\(x_i=0\) ,则这两个事件是互斥的。

  • 对于限制 \(x_i\) 的每个命题 \(a\)\(b\) ,一定有一个为真,则可以写成

\[\lnot a \Rightarrow b \\ \lnot b \Rightarrow a \]

那么由此可连边建图 \((\lnot a,b)\)\((\lnot b,a)\)

  • 在该图中,若节点 \(i\)\(i+n\) 在同一个强连通分量中,即允许互相到达,则它们分别代表的互斥事件会同时发生,说明存在矛盾,即不存在一组合法的方案。

  • 否则则有解。

  • 构造合法解:

    • 对原图缩点得到一张 \(DAG\)

    • 对于变量 \(x_i\),考察节点 \(i\)\(i+n\) 所在强连通分量的拓扑关系。若两分量不连通,则 \(xi\) 可取任意一个值。否则只能取属于拓扑序较大的分量的值。因为若取拓扑序较小的值,可以根据逻辑关系推出取另一个值也是同时发生的。

      by @LuckyBlock

  • \(Tarjan\) 算法赋给强连通分量的编号顺序与拓扑序是相反的,上文已有说明。

CODE

/*

Name: P4782 【模板】2-SAT 问题

By Frather_

*/
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <stack>
using namespace std;
/*=========================================快读*/
int read()
{
    int x = 0, f = 1;
    char c = getchar();
    while (c < '0' || c > '9')
    {
        if (c == '-')
            f = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9')
    {
        x = (x << 3) + (x << 1) + (c ^ 48);
        c = getchar();
    }
    return x * f;
}
/*=====================================定义变量*/
int n, m;

const int _ = 5000050;

struct edge
{
    int to;
    int nxt;
} e[_];
int cnt, head[_];

int dfn[_], low[_], num;
int bel[_], b_num;

stack<int> s;
/*===================================自定义函数*/
void add(int from, int to)
{
    e[++cnt].to = to;
    e[cnt].nxt = head[from];
    head[from] = cnt;
}

void Tarjan(int u_)
{
    dfn[u_] = low[u_] = ++num;
    s.push(u_);
    for (int i = head[u_]; i; i = e[i].nxt)
    {
        int v_ = e[i].to;
        if (!dfn[v_])
        {
            Tarjan(v_);
            low[u_] = min(low[u_], low[v_]);
        }
        else if (!bel[v_])
            low[u_] = min(low[u_], dfn[v_]);
    }
    if (dfn[u_] == low[u_])
    {
        b_num++;
        while (true)
        {
            int t = s.top();
            s.pop();
            bel[t] = b_num;
            if (t == u_)
                break;
        }
    }
    return;
}
/*=======================================主函数*/
int main()
{
    n = read();
    m = read();
    for (int i = 1; i <= m; i++)
    {
        int x = read();
        int a = read();
        int y = read();
        int b = read();
        if (a && b)
        {
            add(x + n, y);
            add(y + n, x);
        }
        if (!a && b)
        {
            add(x, y);
            add(y + n, x + n);
        }
        if (!a && !b)
        {
            add(x, y + n);
            add(y, x + n);
        }
        if (a && !b)
        {
            add(x + n, y + n);
            add(y, x);
        }
    }

    for (int i = 1; i <= n * 2; i++)
    {
        if (!dfn[i])
            Tarjan(i);
        if (i <= n && bel[i] == bel[i + n])
        {
            printf("IMPOSSIBLE\n");
            return 0;
        }
    }
    printf("POSSIBLE\n");
    for (int i = 1; i <= n; i++)
        printf("%d ", bel[i] < bel[i + n]);
    return 0;
}

写在最后

  • 最近被教练抓颓抓得好苦啊\kk

  • 再次向已逝去的学长致敬(((不是

  • 鸣谢:

    • 《算法竞赛进阶指南》

    • OI-Wiki

    • @LuckyBlock

posted @ 2021-03-07 22:02  Frather  阅读(108)  评论(4编辑  收藏  举报