poj - 3683 - Priest John's Busiest Day(2-SAT)

题意:有N场婚礼,每场婚礼的开始时间为Si,结束时间为Ti,每场婚礼有个仪式,历时Di,这个仪式要么在Si时刻开始,要么在Ti-Di时刻开始,问能否安排每场婚礼举行仪式的时间,使主持人John能参加所有的这些仪式的全过程。

题目链接:http://poj.org/problem?id=3683

——>>每场婚礼的仪式,要么在开始段举行,要么在结束段举行,且一定要举行,要求各场婚礼仪式没冲突——>>2-SAT。。。

2-SAT挺神,针对此类问题,可谓手到擒来。。。

LJ《训练指南》上的写法挺容易理解的。。。于是用上了。。。(相对于2003年伍昱论文中O(n)的算法,在时间上《训练指南》中的写法比不上)。。。

对于2-SAT,建图非常重要。。。

——>>设一场婚礼为i,mark[2*i] == 1表示在其开始段举行仪式,mark[2*i+1] == 1表示在其结束段举行仪式。。。

建图思路:对于一个仪式i和另一个仪式j,若i与j有冲突,则说明i不能举行或者j不能举行。。。即i == 0 || j == 0,所以,i' -> j,j' -> i。。。

 

#include <cstdio>

#include <cstring>

#include <vector>



using namespace std;



const int maxn = 1000 + 10;



int N;

int S[maxn], S1[maxn], S2[maxn], T[maxn], T1[maxn], T2[maxn], D[maxn];

int t[maxn][2];



struct Twoset {

    int n;

    vector<int> G[maxn*2];

    bool mark[maxn*2];

    int stk[maxn*2], c;



    void init(int n) {

        this->n = n;

        for(int i = 0; i < 2*n; 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^1].push_back(y);

        G[y^1].push_back(x);

    }



    bool dfs(int x) {

        if(mark[x^1]) return false;

        if(mark[x]) return true;

        mark[x] = true;

        stk[++c] = x;

        int sz = G[x].size();

        for(int i = 0; i < sz; i++) {

            int v = G[x][i];

            if(!dfs(v)) return false;

        }

        return true;

    }



    bool YES() {

        for(int i = 0; i < 2*n; i += 2) if(!mark[i] && !mark[i+1]) {

            c = 0;

            if(!dfs(i)) {

                while(c) mark[stk[c--]] = false;

                if(!dfs(i+1)) return false;

            }

        }

        return true;

    }



    void solve() {

        if(YES()) {

            puts("YES");

            for(int i = 0; i < 2*n; i++) if(mark[i]) {

                int id = i / 2;

                if(i&1) printf("%02d:%02d %02d:%02d\n", (T[id]-D[id])/60, (T[id]-D[id])%60, T[id]/60, T[id]%60);

                else printf("%02d:%02d %02d:%02d\n", S[id]/60, S[id]%60, (S[id]+D[id])/60, (S[id]+D[id])%60);

            }

        }

        else puts("NO");

    }

} solver;



bool isok(int L1, int R1, int L2, int R2) {     //判断两个仪式是否冲突

    return R1 <= L2 || R2 <= L1;

}



void read() {

    for(int i = 0; i < N; i++) {

        scanf("%d:%d %d:%d %d", &S1[i], &S2[i], &T1[i], &T2[i], &D[i]);

        S[i] = S1[i] * 60 + S2[i];

        T[i] = T1[i] * 60 + T2[i];

        t[i][0] = S[i];

        t[i][1] = T[i] - D[i];

    }

}



void build() {

    for(int i = 0; i < N; i++) for(int a = 0; a < 2; a++)

        for(int j = i+1; j < N; j++) for(int b = 0; b < 2; b++)

            if(!isok(t[i][a], t[i][a]+D[i], t[j][b], t[j][b]+D[j]))     //有冲突,不能同时赋a,b,所以a^1或者b^1成立

                solver.add_clause(i, a^1, j, b^1);

}



int main()

{

    while(scanf("%d", &N) == 1) {

        solver.init(N);

        read();

        build();

        solver.solve();

    }

    return 0;

}


 

 

你可能感兴趣的:(poj)