题意:有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;
}