2-SAT问题介绍求解 + 模板题P4782
什么是2-SAT问题
sat 即 Satisfiability,意思为可满足,那么2-SAT表示一些布尔变量只能取true或者false,而某两个变量之间的值存在一定的关系(如:只要a为真,b一定为假;如果a为假,b也为假),我们需要在满足所有这样的关系的情况下,求出每个变量的赋值,如果不存在,就是无解。
举个栗子
某一天,PC,YD,HL在讨论两个问题:1)winter的温度是否低于0度,2)bamboo的高度是否大于10m,众所周知,这三位大佬很强大,所以他们对于这两个问题的判断至少有一个是正确的。此时给出三人对这两个问题的判断,求出这两个问题的正确答案;如果不存在正确答案,那么可能是某个人犯糊涂了,因此我们无法得到答案,输出无解。
我们将两个问题符号化为,a:winter的温度低于0度, b :bamboo的高度大于10m,将三人的观点及其符号化表示如下:
1)PC认为:winter的温度低于0度,bamboo的高度大于10m,a∨ b
2)YD认为:wintet的温度不低于0度,bamboo的高度大于10m,¬a∨ b
3)HL认为:winter的温度不低于0度,bamboo的高度不大于10m,¬a∨ ¬b
那么a,b的取值需要满足:(a∨ b) ∧ ( ¬a∨ b) ∧ (¬a∨ ¬b),我们如何求出满足这个式子的a,b的取值这类问题就是2-SAT问题。
将2-SAT问题转为图论问题求解
我们可以将每个变量x的两个状态分别用两个点表示,记编号为i的点表示这个变量取真,i+n的点表示这个变量取假
将两个变量之间的关系用边表示,如:只要a为真,b一定为假,将此命题符号化得到¬a∨¬b,这个式子也可以写成:a → ¬b ∧ b → ¬a,表示:a为真,则b必为假 ,和b为真,则a必为假,这样我们就发现 a 和 ¬b 可以由a推出¬b, 可以由b推出¬a,在图中,我们构建这样的边表示他们的关系:由 a 向 ¬b 建一条单向边,再由 b 向 ¬a建一条单向边,总结以下将命题转化为建图的规律:
1)¬a∨¬b ----> a → ¬b ∧ b → ¬a
2)a∨b ----> ¬a → b ∧ ¬b → a
3)¬a∨b ----> a → b ∧ ¬b → ¬a
4)a∨¬b ----> ¬a → ¬b ∧ b → a
(x→y可以视作由x向y建一条单向边)
然后,我们在这个图中求强连通分量,联想到我们图中边代表的关系:由边的起点可以推出终点,所以同一强连通分量中的点真值一致,所以我们很容易想到如果 a 和 ¬a在同一强连通分量中,这说明a和¬a真值相同,这显然是不正确的,即无解;如果不存在变量x使得x和¬x在同一强连通分量中,就说明有解。
那么有解的情况下,我们如何得到每个变量的值呢?只要x所在强连通分量的拓扑序比¬x所在强连通分量的拓扑序靠后,则x为真,否则为假,而我们用tarjan算法求强连通分量的时候,对于每个强连通分量的标记是逆拓扑序的,所以 node[x] < node[¬x] 时,x取真,否则为假,这样一来,我们对每个变量进行判断并输出对应的值即可
代码区
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#include<string>
#include<fstream>
#include<vector>
#include<stack>
#include <map>
#include <iomanip> #define bug cout << "**********" << endl
#define show(x, y) cout<<"["<<x<<","<<y<<"] "
#define LOCAL = 1;
using namespace std;
typedef long long ll;
const ll inf = 2e9 + ;
const ll mod = 1e9 + ;
const int Max = 2e6 + ;
const int Max2 = 3e2 + ; struct Egde
{
int to, next;
} edge[Max]; int n, m;
int head[Max], tot;
int dfn[Max], low[Max],time_clock;
int line[Max],now;
int node[Max],sccCnt; void init()
{
memset(head, -, sizeof(head));
tot = ;
memset(dfn,,sizeof(dfn));
time_clock = ;
now = ;
memset(node,,sizeof(node));
sccCnt = ;
} void add(int u, int v)
{
edge[tot].to = v;
edge[tot].next = head[u];
head[u] = tot++;
} void tarjan(int u)
{
dfn[u] = low[u] = ++time_clock;
line[++now] = u;
for(int i = head[u] ; i != -; i = edge[i].next)
{
int v = edge[i].to;
if(!dfn[v])
{
tarjan(v);
low[u] = min(low[u],low[v]);
}
else if(!node[v])
{
low[u] = min(low[u],dfn[v]);
}
}
if(dfn[u] == low[u])
{
sccCnt++;
while(line[now] != u)
node[line[now--]] = sccCnt;
node[line[now--]] = sccCnt;
}
} int main()
{
#ifdef LOCAL
// freopen("input.txt", "r", stdin);
// freopen("output.txt", "w", stdout);
#endif
while (scanf("%d%d", &n, &m) != EOF)
{
init();
for (int i = , a, va, b, vb; i <= m; i++)
{
scanf("%d%d%d%d", &a, &va, &b, &vb);
add(a + n * (va & ), b + n * (vb ^ ));
add(b + n * (vb & ), a + n * (va ^ )); //1~n表示1,n+1,2*n表示0
}
for(int i = ;i <= (n << ) ;i ++)
{
if(!dfn[i])
tarjan(i);
}
bool ok = true;
for(int i = ;i <= n ;i ++)
{
if(node[i] == node[i+n])
{
ok = false;
break;
}
}
if(!ok)
{
printf("IMPOSSIBLE\n");
}
else
{
printf("POSSIBLE\n");
for(int i = ;i < n ;i ++)
{
printf("%d ",node[i] < node[i+n]);
}
printf("%d\n",node[n] < node[n<<]);
}
}
return ;
}
2-SAT问题介绍求解 + 模板题P4782的更多相关文章
- poj3250(单调栈模板题)
题目链接:https://vjudge.net/problem/POJ-3250 题意:求序列中每个点右边第一个>=自身的点的下标. 思路:简单介绍单调栈,主要用来求向左/右第一个小于/大于自身 ...
- bzoj2049-洞穴勘测(动态树lct模板题)
Description 辉辉热衷于洞穴勘测.某天,他按照地图来到了一片被标记为JSZX的洞穴群地区.经过初步勘测,辉辉发现这片区域由n个洞穴(分别编号为1到n)以及若干通道组成,并且每条通道连接了恰好 ...
- [置顶] 小白学习KM算法详细总结--附上模板题hdu2255
KM算法是基于匈牙利算法求最大或最小权值的完备匹配 关于KM不知道看了多久,每次都不能完全理解,今天花了很久的时间做个总结,归纳以及结合别人的总结给出自己的理解,希望自己以后来看能一目了然,也希望对刚 ...
- POJ 1985 Cow Marathon (模板题)(树的直径)
<题目链接> 题目大意: 给定一颗树,求出树的直径. 解题分析:树的直径模板题,以下程序分别用树形DP和两次BFS来求解. 树形DP: #include <cstdio> #i ...
- HDU 4825 Xor Sum (模板题)【01字典树】
<题目链接> 题目大意: 给定n个数,进行m次查找,每次查找输出n个数中与给定数异或结果最大的数. 解题分析: 01字典树模板题,01字典树在求解异或问题上十分高效.利用给定数据的二进制数 ...
- POJ 3264 Balanced Lineup(模板题)【RMQ】
<题目链接> 题目大意: 给定一段序列,进行q次询问,输出每次询问区间的最大值与最小值之差. 解题分析: RMQ模板题,用ST表求解,ST表用了倍增的原理. #include <cs ...
- SPOJ RPLN (模板题)(ST算法)【RMQ】
<题目链接> 题目大意:给你一段序列,进行q次区间查询,每次都输出询问区间内的最小值. 解题分析: RMQ模板题,下面用在线算法——ST算法求解.不懂ST算法的可以看这篇博客 >& ...
- POJ 1330 Nearest Common Ancestors (模板题)【LCA】
<题目链接> 题目大意: 给出一棵树,问任意两个点的最近公共祖先的编号. 解题分析:LCA模板题,下面用的是树上倍增求解. #include <iostream> #inclu ...
- POJ2417 Discrete Logging【BSGS】(模板题)
<题目链接> 题目大意: P是素数,然后分别给你P,B,N三个数,然你求出满足这个式子的L的最小值 : BL== N (mod P). 解题分析: 这题是bsgs算法的模板题. #incl ...
随机推荐
- Spring——JdbcTemplate
一.JdbcTemplate介绍: 为了使 JDBC 更加易于使用,Spring 在 JDBCAPI 上定义了一个抽象层, 以此建立一个JDBC存取框架,Spring Boot Spring Data ...
- mock的那点事
前言: Mock在GitHub上有12.9K的star可以看出,它在技术团队中是挺受欢迎的.这项技术被应用在不同领域的项目中. 适用场景: 下面我结合我们技术团队,列举最适合引入我们Mock服务的场景 ...
- [题解] [AtCoder2134] Zigzag MST
题面 题解 考虑kruscal的过程 对于三个点\(x, y, x + 1\), 我们可以将\((x, y, z), (y, x + 1, z + 1)\)看做\((x, y, z), (x, x + ...
- mini dc(选做)
一.任务详情 提交测试截图和码云练习项目链接,实现Linux下dc的功能,计算后缀表达式的值 二.源代码 1.MyDC类 import java.util.StringTokenizer; impor ...
- FinalCutPro快捷键
FinalCutPro快捷键使用 FinalCutPro的快捷键使用十分有用,特对一些基本的快捷键进行了总结 1)i:截取片段开始Initial 2)o: 截取片段结束Over i和o可以在一个素材片 ...
- 黑马vue---8-10、v-cloak、v-text、v-html、v-bind、v-on的基本使用
黑马vue---8-10.v-cloak.v-text.v-html.v-bind.v-on的基本使用 一.总结 一句话总结: v-bind等这些东西都是用的vue.data里面的变量 1.使用 v- ...
- zeppelin 无法连接一个已有的standalone模式的spark集群
SparkInterpreter.java 这个文件里面读取master的属性有些问题: 原来代码中"master"属性的获取的地方应该是错了.设置和读取这个属性的对象不是同一个 ...
- AppiumLibrary移动APP测试
使用Genymotion模拟器结合RF执行 前提搭建环境参考<python_Appium测试环境搭建>文章详细介绍. 常用关键字 关 键 字 描 述 Click Button 点击 ...
- Selenium chromeDriver 下载地址
http://chromedriver.storage.googleapis.com/ http://npm.taobao.org/mirrors/chromedriver/
- Android应用开发编译框架流程与IDE及Gradle概要
1 背景 建议阅读本文之前先阅读<Android Studio入门到精通>和<Groovy脚本基础全攻略>及<Gradle脚本基础全攻略>三篇博客作为背景知识,这样 ...