2-SAT (two-statisfiability) 算法 学习笔记

$2-SAT$问题指的是对于若干限制求出一组可行解的问题。

考虑对于$n$个值域为${0,1}$的变量$x_1 , x_2 ,...,x_n$ 满足若干限制:

若 $x_i = p$ 则 $x_j = q ( i,j\in[1,n],p,q \in \{0,1\})$

我们考虑对于每一个变量$x_i$开一个值域$x_{i,0} , x_{i,1}$表示第$i$个变量取值为$0/1$的点。

然后考虑每一组命题, 若 $x_i = p$ 则 $x_j = q , i,j\in[1,n],p,q \in \{0,1\}$ ,

首先把$x_{i,p}$和$x_{j,q}$两个点连一条单向边。

并且 由于这个命题非常特殊,值域大小只为2 , 其逆否命题也是限制(可以显然的反证)。

于是就把$x_{j,1-q}$和$x_{i,1-p}$ 连一条单向边。

 

然后对于上面的图跑tarjan找出scc,如果$x_{i,0}$和$x_{i,1}$在同一连通块中了,

说明下列命题成立: “若$x_i = 0$则$x_i = 1$ ” ,“若$x_i = 1$则$x_i = 0$ ” 所以此时答案无解。

 

如何构造出一组解呢,由于tarjan的dfs特性,

设$c_i$表示通过tarjan算法求出的连通块编号(这本身就是逆拓扑序的)。

$val_i = c_{x_{i,0}} > c_{x_{i,1}}$  就构造出$val_i , i\in [1,n]$一组合法解了。

P4782 【模板】2-SAT 问题

设变量$x_i \in \{0,1\}$ ,给出若干组关系:$x_i = p $ 或者 $x_j = q$ 

如果$x_i , (i\in [1,n])$有解,先输出"POSIBLE"然后输出一组合法解。

如果$x_i ,  (i\in [1,n])$无解,则输出"IMPOSIBLE"即可。

对于100%的数据$n,m\leq 10^6$

Sol: 

$x_i$为$p$ 或 $x_j$为$q$ 可以拆成$2\times 2$对逻辑关系。

  • 若$x_i$为$1-p$,则$x_j$为$q$ , 若$x_j$为$1-q$,则$x_i$为$p$
  • 若$x_j$为$1-q$,则$x_i$为$p$ , 若$x_i$为$1-p,$ 则$x_j$为$q$

实现方面,只需:$[1,N]$为原来元素的$0$域,$[N+1,2N]$ 为原来元素的$1$域.

# include <bits/stdc++.h>
using namespace std;
const int N=2e6+10,M=4e6+10;
struct rec{ int pre,to;}a[M];
stack<int>s;
bool ins[N];
int cnt,tot,n,m,head[N],c[N],dfn[N],low[N],val[N];
void adde(int u,int v)
{
    a[++tot].pre=head[u];
    a[tot].to=v;
    head[u]=tot;
}
void tarjan(int u)
{
    dfn[u]=low[u]=++dfn[0];
    s.push(u);ins[u]=1;
    for (int i=head[u];i;i=a[i].pre){
        int v=a[i].to;
        if (!dfn[v]) {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        } else if (ins[v]) low[u]=min(low[u],dfn[v]);
    }
    if (dfn[u]==low[u]) {
        cnt++; int v;
        do {
            v=s.top(); s.pop();
            ins[v]=0; c[v]=cnt;
        } while (u!=v);
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++) {
        int s,ps,t,pt; scanf("%d%d%d%d",&s,&ps,&t,&pt);
        adde(s+(1-ps)*n,t+pt*n); adde(t+(1-pt)*n,s+ps*n);
        adde(t+(1-pt)*n,s+ps*n); adde(s+(1-ps)*n,t+pt*n);
    }
    for (int i=1;i<=2*n;i++) if (!dfn[i]) tarjan(i);
    for (int i=1;i<=n;i++) if (c[i]==c[i+n]) {
        puts("IMPOSSIBLE"); return 0;
    }
    for (int i=1;i<=n;i++) val[i]=c[i]>c[n+i];
    puts("POSSIBLE");
    for (int i=1;i<=n;i++) printf("%d ",val[i]); puts("");
    return 0;
}

 

转载于:https://www.cnblogs.com/ljc20020730/p/11295659.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值