2-SAT总结

最新推荐文章于 2022-05-14 10:31:14 发布
最新推荐文章于 2022-05-14 10:31:14 发布
阅读量498 收藏
点赞数
分类专栏： ACM 文章标签： 二分 2-SAT SCC
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/randygx/article/details/25561491
版权

看到HDOJ上有一个2-SAT专题，刷完之后顺便总结一下。

最简单的2-SAT问题就是判断是否有解，回答YES或NO就行了，不需要输出具体的解。这种问题建图完毕后一个SCC（强连通分量）就好了。如果要输出解，就多一个拓扑排序。代码中涉及多次建图，注意细节就好了。

1.HDOJ 3062 Party

很裸的2-SAT模板，而且不用输出路径。 
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
        
        
#include 
        
        
          #include 
          #include 
          
            #include 
           
             #include 
            
              #include 
             
               #include 
              
                #include 
               
                 #include 
                
                  #include 
                 
                   #include 
                  
                    #include 
                   
                     #define mp make_pair #define MEMSET(a, b) memset(a, b, sizeof(a)) using namespace std; typedef unsigned int ui; typedef long long ll; typedef unsigned long long ull; typedef pair 
                    
                      pii; typedef vector 
                     
                       vi; typedef vi::iterator vi_it; typedef map 
                      
                        mii; typedef priority_queue 
                       
                         pqi; typedef priority_queue 
                        
                          , greater 
                         
                           > rpqi; const int MAX_N = 2000 + 10; int link[MAX_N][MAX_N], rlink[MAX_N][MAX_N]; int stk[MAX_N]; bool vis[MAX_N]; int belong[MAX_N]; int top, cnt; void dfs_visit(int (*lk)[MAX_N], int u, bool ft) { vis[u] = true; for (int i = 1; i <= lk[u][0]; ++i) { if (!vis[lk[u][i]]) { dfs_visit(lk, lk[u][i], ft); } } if (ft) { stk[top++] = u; } else { belong[u] = cnt; } } void scc(int n) { top = 0; MEMSET(vis, 0); for (int i = 0; i < n; ++i) { if (!vis[i]) { dfs_visit(link, i, true); } } cnt = 0; MEMSET(vis, 0); while (top) { if (!vis[stk[--top]]) { dfs_visit(rlink, stk[top], false); ++cnt; } } } bool two_sat(int n) { scc(n << 1); for (int i = 0; i < n; ++i) { if (belong[i] == belong[i + n]) { return false; } } return true; } int main(int argc, char *argv[]) { int n, m, i, a1, a2, c1, c2; while (cin >> n >> m) { for (i = 0; i < (n << 1); ++i) { rlink[i][0] = link[i][0] = 0; } while (m--) { scanf("%d%d%d%d", &a1, &a2, &c1, &c2); int b1 = a1 + n * (c1 ^ 1); int b2 = a2 + n * (c2 ^ 1); a1 += n * c1; a2 += n * c2; link[a1][++link[a1][0]] = b2; link[a2][++link[a2][0]] = b1; rlink[b2][++rlink[b2][0]] = a1; rlink[b1][++rlink[b1][0]] = a2; } if (two_sat(n)) { puts("YES"); } else { puts("NO"); } } return 0; }

这个题要求输出一组解，另外需要根据区间的交来建图，代码比较长。判断区间是否相交时我先排序，这样内层循环不用遍历所有j，不过效率其实差不多，只快了一点点。。。

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
         #include 
         #include 
         
           #include 
          
            #include 
           
             #include 
            
              #include 
             
               #include 
              
                #include 
               
                 #include 
                
                  #include 
                 
                   #include 
                  
                    #include 
                   
                     #define mp make_pair #define X first #define Y second #define MEMSET(a, b) memset(a, b, sizeof(a)) using namespace std; typedef unsigned int ui; typedef long long ll; typedef unsigned long long ull; typedef pair 
                    
                      pii; typedef vector 
                     
                       vi; typedef vi::iterator vi_it; typedef map 
                      
                        mii; typedef priority_queue 
                       
                         pqi; typedef priority_queue 
                        
                          , greater 
                         
                           > rpqi; const int MAX_N = 2000 + 10; int link[MAX_N][MAX_N], rlink[MAX_N][MAX_N]; int stk[MAX_N], ans[MAX_N], op[MAX_N]; bool vis[MAX_N], flag[MAX_N][MAX_N]; int belong[MAX_N]; int top, cnt; pair 
                          
                            cere[MAX_N]; void dfs_visit(int (*lk)[MAX_N], int u, bool ft) { vis[u] = true; for (int i = 1; i <= lk[u][0]; ++i) { if (!vis[lk[u][i]]) { dfs_visit(lk, lk[u][i], ft); } } if (ft) { stk[top++] = u; } else { belong[u] = cnt; } } void scc(int n) { top = 0; MEMSET(vis, 0); for (int i = 0; i < n; ++i) { if (!vis[i]) { dfs_visit(link, i, true); } } cnt = 0; MEMSET(vis, 0); while (top) { if (!vis[stk[--top]]) { dfs_visit(rlink, stk[top], false); ++cnt; } } } bool two_sat(int n) { scc(n << 1); for (int i = 0; i < n + n; i += 2) { if (belong[i] == belong[i + 1]) { return false; } } return true; } void color(int u) { ans[u] = 2; for (int i = 1; i <= rlink[u][0]; ++i) { if (!ans[rlink[u][i]]) { color(rlink[u][i]); } } } bool cmp(const pair 
                           
                             &p1, const pair 
                            
                              &p2) { return p1.Y < p2.Y; } int main(int argc, char *argv[]) { // freopen("D:\\in.txt", "r", stdin); int n, i, j; cin >> n; for (i = 0; i < n; ++i) { int s1, s2, t1, t2, d; scanf("%d:%d %d:%d %d", &s1, &s2, &t1, &t2, &d); cere[i + i] = mp(mp(s1 * 60 + s2, s1 * 60 + s2 + d), i + i); cere[i + i + 1] = mp(mp(t1 * 60 + t2 - d, t1 * 60 + t2), i + i + 1); } sort(cere, cere + n + n); for (i = 0; i < n + n; ++i) { for (j = i + 1; j < n + n && cere[j].X.X < cere[i].X.Y; ++j) { int u = cere[i].Y, v = cere[j].Y ^ 1; if (!flag[u][v]) { flag[u][v] = true; link[u][++link[u][0]] = v; rlink[v][++rlink[v][0]] = u; } u = cere[j].Y, v = cere[i].Y ^ 1; if (!flag[u][v]) { flag[u][v] = true; link[u][++link[u][0]] = v; rlink[v][++rlink[v][0]] = u; } } } if (!two_sat(n)) { puts("NO"); return 0; } else { puts("YES"); } MEMSET(flag, 0); for (i = 0; i < n + n; ++i) { rlink[i][0] = 0; op[belong[i]] = belong[i ^ 1]; } for (i = 0; i < n + n; ++i) { int u = belong[i]; for (j = 1; j <= link[i][0]; ++j) { int v = belong[link[i][j]]; if (u != v && !flag[v][u]) { flag[v][u] = true; rlink[v][++rlink[v][0]] = u; } } } MEMSET(vis, 0); top = 0; for (i = 0; i < cnt; ++i) { if (!vis[i]) { dfs_visit(rlink, i, true); } } while (top) { if (!ans[stk[--top]]) { ans[stk[top]] = 1; color(op[stk[top]]); } } sort(cere, cere + n + n, cmp); for (i = 0; i < n + n; ++i) { if (ans[belong[i]] == 1) { printf("%02d:%02d %02d:%02d\n", cere[i].X.X / 60, cere[i].X.X % 60, cere[i].X.Y / 60, cere[i].X.Y % 60); } } return 0; }

和上题差不多的。 
#include 
   
   
    
    
#include 
    
    
     
     
 
const int maxn=10010;
int n, m, nxt[maxn], head[maxn], pnt[maxn], ne, e, a, b, a0, b0;
char c1, c2;
int nnxt[maxn], nhead[maxn], npnt[maxn];
int order[maxn], norder, id[maxn], v[maxn];
int ans[maxn], op[maxn];
 
void dfs(int d){
    v[d] = 1;
    for(int i=head[d]; i!=-1; i=nxt[i])
        if(!v[pnt[i]])
            dfs(pnt[i]);
    order[norder++] = d;
}
void ndfs(int d, int k){
    v[d] = 1;
    id[d] = k;
    for(int i=nhead[d]; i!=-1; i=nnxt[i])
        if(!v[npnt[i]])
            ndfs(npnt[i], k);
}
void addedge(int s, int t){
    pnt[e] = t; nxt[e] = head[s]; head[s] = e++;
}
void addnedge(int t, int s){
    npnt[ne] = s; nnxt[ne] = nhead[t]; nhead[t] = ne++;
}
void color(int d){
    ans[d] = 2;
    for(int i=head[d]; i!=-1; i=nxt[i])
        if(!ans[pnt[i]])
            color(pnt[i]);
}
int main(){
    while(1){
        norder = e = ne = 0;
        memset(head, -1, sizeof head);
        memset(nhead, -1, sizeof nhead);
        scanf("%d%d", &n, &m);
        if(!n&&!m)
            break;
        for(int i=0; i
     
     
      
      =0; --i)
            if(!v[order[i]])
                ndfs(order[i], k++);
        int mark = 1;
        for(int i=0; i
      
      
       
       =0; --i)
            if(!ans[order[i]]){
                ans[order[i]] = 1;
                color(op[order[i]]);
            }
        for(int i=1; i

4.HDOJ 3715 Go Deeper

二分答案之后判断2-SAT问题是否有解。SCC用tarjan写的。

#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
        
        
#include 
        
        
          #include 
          #include 
          
            #include 
           
             #include 
            
              #include 
             
               #include 
              
                #include 
               
                 #include 
                
                  #include 
                 
                   #include 
                  
                    #include 
                   
                     #include 
                    
                      #define mp make_pair #define X first #define Y second #define MEMSET(a, b) memset(a, b, sizeof(a)) using namespace std; typedef unsigned int ui; typedef long long ll; typedef unsigned long long ull; typedef pair 
                     
                       pii; typedef vector 
                      
                        vi; typedef vi::iterator vi_it; typedef map 
                       
                         mii; typedef priority_queue 
                        
                          pqi; typedef priority_queue 
                         
                           , greater 
                          
                            > rpqi; const int MAX_N = 20000 + 10; int head[MAX_N]; int edge[MAX_N * 100][2]; int dfn[MAX_N], low[MAX_N]; int stk[MAX_N]; int belong[MAX_N]; int top, cnt, idx, edgecnt; int a[MAX_N], b[MAX_N], c[MAX_N]; int N; void dfs(int u) { dfn[u] = low[u] = ++idx; stk[top++] = u; int v; for (int i = head[u]; i != -1; i = edge[i][1]) { if (!dfn[v = edge[i][0]]) { dfs(v); low[u] = min(low[u], low[v]); } else if (!belong[v]) { low[u] = min(low[u], dfn[v]); } } if (dfn[u] == low[u]) { belong[u] = ++cnt; while (stk[--top] != u) { belong[stk[top]] = cnt; } } } void tarjan(int n) { MEMSET(belong, 0); MEMSET(dfn, 0); cnt = idx = top = 0; for (int i = 1; i <= n; ++i) { if (!dfn[i]) { dfs(i); } } } bool two_sat(int n) { tarjan(n + n); for (int i = 1; i <= n; ++i) { if (belong[i] == belong[i + n]) { return false; } } return true; } inline void add_edge(int u, int v) { edge[edgecnt][0] = v; edge[edgecnt][1] = head[u]; head[u] = edgecnt++; } void build(int n) { int i; MEMSET(head, -1); edgecnt = 0; for (i = 1; i <= n; ++i) { if (c[i] == 0) { add_edge(a[i], b[i] + N); add_edge(b[i], a[i] + N); } else if (c[i] == 1) { add_edge(a[i], b[i]); add_edge(b[i], a[i]); add_edge(a[i] + N, b[i] + N); add_edge(b[i] + N, a[i] + N); } else { add_edge(a[i] + N, b[i]); add_edge(b[i] + N, a[i]); } } } int main(int argc, char *argv[]) { // freopen("D:\\in.txt", "r", stdin); int t, n, m, i, j; cin >> t; while (t--) { scanf("%d%d", &n, &m); N = n; for (i = 1; i <= m; ++i) { scanf("%d%d%d", a + i, b + i, c + i); } int low = 0, high = m, ans = 0; while (low <= high) { int mid = (low + high) >> 1; build(mid); if (two_sat(N)) { low = mid + 1; ans = mid; } else { high = mid - 1; } } printf("%d\n", ans); } return 0; }

5.POJ 3207 Ikki's Story IV - Panda's Trick

也是只需要判断，不过注意弦相交的判断方法。

#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
        
        
#include 
        
        
          #include 
          #include 
          
            #include 
           
             #include 
            
              #include 
             
               #include 
              
                #include 
               
                 #include 
                
                  #include 
                 
                   #include 
                  
                    #include 
                   
                     #define mp make_pair #define MEMSET(a, b) memset(a, b, sizeof(a)) using namespace std; typedef unsigned int ui; typedef long long ll; typedef unsigned long long ull; typedef pair 
                    
                      pii; typedef vector 
                     
                       vi; typedef vi::iterator vi_it; typedef map 
                      
                        mii; typedef priority_queue 
                       
                         pqi; typedef priority_queue 
                        
                          , greater 
                         
                           > rpqi; const int MAX_M = 1000 + 2; int line[MAX_M][2]; vector 
                          
                            link[MAX_M]; int dfn[MAX_M], low[MAX_M], belong[MAX_M]; int stk[MAX_M]; int cnt = 0, top = 0, idx = 0; bool intersect(int *a, int *b, int n) { if (a[0] < b[0] && b[0] < a[1] && a[1] < b[1]) { return true; } if (b[0] < a[0] && a[0] < b[1] && b[1] < a[1]) { return true; } return false; } void dfs(int u) { dfn[u] = low[u] = ++idx; stk[top++] = u; int v; for (int i = 0; i < link[u].size(); ++i) { if (!dfn[v = link[u][i]]) { dfs(v); low[u] = min(low[u], low[v]); } else if (!belong[v]) { low[u] = min(low[u], dfn[v]); } } if (dfn[u] == low[u]) { belong[u] = ++cnt; while (stk[--top] != u) { belong[stk[top]] = cnt; } } } void tarjan(int n) { for (int i = 0; i < n; ++i) { if (!dfn[i]) { dfs(i); } } } bool two_sat(int n) { tarjan(n + n); for (int i = 0; i < n; ++i) { if (belong[i] == belong[i + n]) { return false; } } return true; } int main(int argc, char *argv[]) { // freopen("D:\\in.txt", "r", stdin); int n, m; cin >> n >> m; for (int i = 0; i < m; ++i) { cin >> line[i][0] >> line[i][1]; if (line[i][0] > line[i][1]) { swap(line[i][0], line[i][1]); } for (int j = 0; j < i; ++j) { if (intersect(line[i], line[j], n)) { link[i].push_back(j + m); link[j].push_back(i + m); link[i + m].push_back(j); link[j + m].push_back(i); } } } if (two_sat(m)) { puts("panda is telling the truth..."); } else { puts("the evil panda is lying again"); } return 0; }

下面几题还没做，先给出链接




randygx
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
2-sat思想及入门
qq_43041641的博客
07-14 971
2-sat思想及入门什么是2-sat例题分析建边判断有解输出方案完整代码 什么是2-sat 要知道2-sat是什么，我们先要知道什么是适定性（Satisfiability）问题，适，就是合适，定，就是确定，适定性问题通俗的说就是确定是否可以满足所有的条件？或者说就是确定一个满足所有条件的方案。取英文的前三个字母，简称sat问题。 通俗的sat问题表述一般是这样的： 有很多个集合，每个集合里面有若...
2-sat算法总结
m725kk的博客
08-18 416
2-sat算法总结 定义：2−sat2-sat2−sat类问题，简单来说就是有nnn个集合，每个集合有两个元素。已知若干个<a,b><a,b><a,b>，表示aaa与bbb矛盾（a、ba、ba、b不同集合），然后从每个集合中取出一个元素，判断是否能够取出nnn个元素。 思路：假设有集合a、ba、ba、b，分别有a1、a2、b1、b2a_1、a_2、b_1、b_2a1​、a2​、b1​、b2​。若a1a_1a1​与b2b_2b2​存在矛盾，那么我们拉两条有向边(a1,b1)
2-SAT问题详解
qq_42376883的博客
05-14 1682
2-sat问题详解
2-SAT问题合集-ybtoj
Mint_hexagram的博客
02-24 2123
2-SAT的合集
2-sat总结
夏天的伞
05-06 536
sat问题 关于定义，详见赵爽的《2-SAT解法浅析》，大概意思是有n个bool变量，一个全集为所有的他们以及他们的非，取其中的若干个变量或其非组成或的表达式，将这些表达式相与，求使与出的值为真的解集，为n-sat问题。 可证明，n>2时为NP完全问题，所以只讨论n=2的情况。 对于一个表达式a+b， a+b=!(!a*!b)  ，即，如果在解集中选择了a，那么必须选择!b，选择了b
2-SAT学习总结
scarlyw的博客
01-15 471
2-SAT 图论算法：   理解： 给定一个布尔方程，判断是否存在一组布尔变量的取值方案，使得整个方程值为真的问题，被称为布尔方程的可满足性问题 (SAT) 。 SAT 问题是 NP 完全的，但对于一些特殊形式的 SAT 问题我们可以有效求解。 我们将下面这种布尔方程称为合取范式：   其中称为文字，它是一个布尔变量或其否定。像这样用连接的部分称为子句。如果合取范式的每个子句中
ACM- 2-SAT问题
胖亚亚的小酒馆
03-03 631
前言：这篇文章是参考着饶齐的总结写出来的，但只有一些文字性的描述类似。 现在有一个由N个布尔值组成的序列A，给储户一些限制关系 比如 A[x] AND A[y] =0 、A[x] OR A[y] OR A[z] =1等， 要确定A[0...N-1]的值，使其满足所有限制关系。这个问题称为2-SAT问题 特别的，若每种限制关系中最多只对两个元素进行限制，则称为2-SAT问题。 由于在...
3.3（差分约束+2-SAT总结）
又菜又爱玩
03-03 274
终于把题基本补完了，头大，，， 发现几个以前理解的不对的地方，比如2-SAT tarjan解法反向拓扑求可行解的过程，再有之前 关于2-SAT中如果p-&gt;q和非q-&gt;非p一定同时存在的时候可以用扩展域并查集代替tarjan在线判断合法性的理解现在看来是不全面的。 （写到后面发现并查集有个地方写错了，改后就A了，有点戏剧，不过目前来看之前的思考还是...
2-sat情况总结
Endeavor_G的博客
11-10 194
2-SAT简介
Rorelynn的博客
01-30 807
前方是2-SAT蒟蒻版（大雾） 蒟蒻不会打代码啊QAQ,代码有参考 配合洛谷食用风味更佳 2-SAT, 就是2-Satisfiability。 Satisfiability，可满足性。 背景故事 我要请大家开party！–ObliaObliaOblia ObliaObliaOblia这样想。 因为ObliaObliaOblia是萌妹子同时是OIer和MOer，所以她可以同时邀请到OI和MO的神...
2-SAT学习资料
09-07
"【2-sat】专题～ - Amb@HDU - 博客园_files"和"2-SAT总结 - kuangbin - 博客园_files"这两个文件夹可能包含了与博客文章相关的图片、代码示例或其他辅助材料，有助于读者更好地理解2-SAT的实现和应用。 总的来说，...
2-sat.rar_2sat
09-21
《2-SAT问题详解与解决方案》 在计算机科学领域，特别是在图论和逻辑满足问题中，2-SAT（2-CNF-SAT）是一个基础且重要的问题。2-SAT问题源于布尔逻辑中的二元约束满足问题（CNF-SAT），但其简化了原问题的复杂性，...
2-sat.zip_2sat算法
09-21
总结，2-SAT算法是解决布尔逻辑约束问题的一个实用工具，它的高效性使得它在实际应用中具有广泛价值。理解2-SAT问题的结构和算法原理，以及如何通过编程实现，对于理解和解决更复杂的布尔逻辑问题具有基础性的指导...
2-sat.rar_2sat_SAT求解
最新发布
09-24
总结来说，2-SAT问题是一个重要的布尔逻辑问题，它的解决方案对于理解和处理布尔逻辑和图论问题具有深远意义。通过深入研究"2-sat.cpp"代码，我们可以学习到如何运用算法高效地解决这类问题，这对于计算机科学和相关...
ACM-2-sat资料
05-24
总结，2-SAT问题是一个在算法竞赛和理论计算机科学中常见的问题，它的理解与解决对于提升逻辑思维和算法能力有着重要作用。通过学习2-SAT，不仅可以掌握基本的图论知识，还能为解决更复杂的SAT问题打下坚实的基础。...
线段树求LIS并统计最长子序列个数
randygx的专栏
05-04 2217
求最长单增子序列(LIS)长度并统计取得最大长度的子序列数。扩展：把LIS换成最长先增后减子序列
2014嘉杰信息杯ACM/ICPC湖南程序设计邀请赛暨第六届湘潭市程序设计竞赛 解题报告
randygx的专栏
06-08 1658
A simple problem w
ZOJ Monthly, June 2014 解题报告
randygx的专栏
06-01 1368
A.Another Recurrence Sequence B.Gears 题目大意：有n个齿轮，一开始各自为一组，之后进行m次操作，包括以下4种类型： 1.合并两组齿轮，合并的两个应该反向旋转 2.把某个齿轮从所在组删除，自为一组，但不影响同组其它齿轮的状态与关系 3.询问两个齿轮是同向、反向或无关系(即不在同一组) 4.询问某个齿轮所在组的齿轮总数 分析：
TOJ 3883: Roadblock -- 单源最短路
randygx的专栏
05-02 1073
3883: Roadblock Time Limit(Common/Java):1000MS/3000MS     Memory Limit:65536KByte Total Submit: 17            Accepted:7 Description Every morning, FJ wakes up and walks across the farm fro
2-SAT问题解析与对称性解法
总结来说，2-SAT问题可以通过对称性解法进行有效解决，特别是在和平委员会这样的问题中，这种解法能够帮助我们理解问题的本质，并提供了一种实用的求解策略。在实际应用中，2-SAT问题的解法可以广泛应用于逻辑推理、...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值