2-SAT+二分——POJ2723

题目描述：

题目描述很简单，有m扇门，每扇门有两个锁，两个锁对应两把钥匙，任意打开一个锁就可以打开一扇门，其中钥匙也是俩俩配对的，也就是说一对钥匙中只要使用了一把，另一把就不能使用了，但是使用过的钥匙可以反复使用。求最多能开多少门。

大致思路：

看过描述之后还是可以看出来，这个题目还是比较明显的2-SAT的题目，主要问题还是在建边。钥匙不用说，所有的钥匙对都要加到图中，而且两个变量的关系应该是如果a是1，就能推出b是0，或者是如果b是1就能推出a是0的关系。

之后就要看门的问题。因为需要求最多能开多少门，所以多多少少会想到将多少门关系边加入到图中决定了答案的多少，但是m的范围是1-2048，如果枚举m的话n^3算法还是不能接受的。之后注意到了一个条件，门必须按顺序开，这样的话就好办了，可以二分m，每次对于二分的m，将前m个门关系对加入到图中，即可。

最后就是如何建立门关系边了。有很多人会想，是不是应该跟要是一样呢？开了a的锁，能推出不用开b，其实不然，题目中说两把锁只要开一个就能打开一扇门，言外之意应该是如果不开a，必须要开b。所以关系确定了，建边也变得简单了，就是如果a是0，那么b就是1，如果b是0，那么a是1，这样的关系。建好边之后求解就好了。

代码：

#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>

using namespace std;

const int maxn = 1024 + 10;
const int maxm = 2048 + 10;
struct TwoSAT {
    int n;
    vector<int> G[maxn*4];
    bool mark[maxn*4];
    int S[maxn*4],c;
    bool dfs(int x) {
        if (mark[x^1]) return false;
        if (mark[x]) return true;
        mark[x] = true;
        S[c++] = x;
        for (int i = 0; i < G[x].size(); i++)
            if (!dfs(G[x][i])) return false;
        return true;
    }
    void init(int n) {
        this->n = n;
        for (int i = 0; i < n*4; i++) G[i].clear();
        memset(mark,0,sizeof(mark));
    }
    void add_clause(int x,int xval,int y,int yval) {
        x = x * 2 + xval;
        y = y * 2 + yval;
        G[x].push_back(y);
        //G[y^1].push_back(x);
    }
    bool solve() {
        for (int i = 0; i < n*2; i += 2)
        if (!mark[i] && !mark[i+1]) {
            c = 0;
            if (!dfs(i)) {
                while (c > 0) mark[S[--c]] = false;
                if (!dfs(i+1)) return false;
            }
        }
        return true;
    }
};

int ka[maxn],kb[maxn],da[maxm],db[maxm];
int n,m;

int main() {
    while (cin>>n>>m,n) {
        for (int i = 0; i < n; i++) scanf("%d%d",&ka[i],&kb[i]);
        for (int i = 0; i < m; i++) scanf("%d%d",&da[i],&db[i]);
        int l = 0,r = m+1;
        while (r - l > 1) {
            int mid = (l + r) / 2;
            TwoSAT ts;
            ts.init(n);
            for (int i = 0; i < n; i++) {
                ts.add_clause(ka[i],1,kb[i],0);
                ts.add_clause(kb[i],1,ka[i],0);
            }
            for (int i = 0; i < mid; i++) {
                ts.add_clause(da[i],0,db[i],1);
                ts.add_clause(db[i],0,da[i],1);
            }
            if (ts.solve()) l = mid;
            else r = mid;
        }
        cout<<l<<endl;
    }
}