類型:2-sat【模闆】
題目:http://poj.org/problem?id=3683
來源:POJ Founder Monthly Contest – 2008.08.31, Dagger and Facer
思路:本題為需要輸出結果的2-sat問題。令i和i + n表示節點i,對于有向邊的構造,如果i和j有沖突,那麼選擇i,j + n必選,選擇j, i + n必選,構造兩條有向邊。其他類似。
然後就是Tarjan求強連通分量,縮點,判斷i,i + n是否屬于同一強連通分量,如果屬于,那麼結果為“NO”。
對于存在解的情況,縮點後,構造新圖【反向建圖】,同時記錄互相沖突的節點,對新圖topo排序,記錄topo序列,然後染色,對于沒有沖突的點染為紅色,其他為藍色。
記錄染色情況,最後列印輸出。
// poj 3683 Priest John's Busiest Day
// re re ac 12196K 125MS
#include <iostream>
#include <string>
#include <queue>
#include <stack>
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;
#define FOR(i,a,b) for(i = (a); i < (b); ++i)
#define FORE(i,a,b) for(i = (a); i <= (b); ++i)
#define FORD(i,a,b) for(i = (a); i > (b); --i)
#define FORDE(i,a,b) for(i = (a); i >= (b); --i)
#define CLR(a,b) memset(a,b,sizeof(a))
#define PB(x) push_back(x)
const int N = 5000;
const int M = 4000010;
int low[N], dfn[N];
int belong[N];
bool inStack[N], vis[N];
stack<int> st;
int n, m;
int step, t;
int conflict[N];//縮點後有沖突的點
int du[N];//入度
int color[N];//1為紅色,-1為藍色,紅色輸出
int cnt1, cnt2;
int head1[N], head2[N];
int top[N], cnt;
struct node {
int s, e;
int len;
}c[N];
struct edge1 {
int v, nxt;
}E1[M];
struct edge2 {
int v, nxt;
}E2[M];
void addedge1(int u, int v) {
E1[cnt1].v = v;
E1[cnt1].nxt = head1[u];
head1[u] = cnt1++;
}
void addedge2(int u, int v) {
E2[cnt2].v = v;
E2[cnt2].nxt = head2[u];
head2[u] = cnt2++;
}
void tarjan(int u) {
int i;
step++;
st.push(u);
low[u] = dfn[u] = step;
vis[u] = 1;
inStack[u] = 1;
for(i = head1[u]; i != -1; i = E1[i].nxt) {
int x = E1[i].v;
if(!vis[x]) {
tarjan(x);
low[u] = min(low[u], low[x]);
}
else if(inStack[x])
low[u]=min(low[u], dfn[x]);
}
if(low[u] == dfn[u]) {
t++;
while(1) {
int x = st.top();
st.pop();
belong[x] = t;
inStack[x] = 0;
if(x == u)
break;
}
}
}
void init() {
cnt1 = cnt2 = step = t = 0;
CLR(head1, -1);
CLR(head2, -1);
CLR(du, 0);
CLR(color, 0);
CLR(vis, 0);
CLR(inStack, 0);
}
bool isConflict(int a, int b, int c, int d) {
if(a >= d || b <= c)
return false;
return true;
}
void Build() {
int i,j;
init();
FORE(i, 1, n - 1)
FORE(j, i + 1, n) {
if(isConflict(c[i].s, c[i].s + c[i].len, c[j].s, c[j].s + c[j].len)) {
addedge1(i, j + n);
addedge1(j, i + n);
}
if(isConflict(c[i].s, c[i].s + c[i].len, c[j].e - c[j].len, c[j].e)) {
addedge1(i, j);
addedge1(j + n, i + n);
}
if(isConflict(c[i].e - c[i].len, c[i].e, c[j].e - c[j].len, c[j].e)) {
addedge1(i + n, j);
addedge1(j + n, i);
}
if(isConflict(c[i].e - c[i].len, c[i].e, c[j].s, c[j].s + c[j].len)) {
addedge1(i + n, j + n);
addedge1(j, i);
}
}
}
void Rebuild(){//逆圖
int i, j;
FORE(i, 1, 2 * n) {
for(j = head1[i]; j != -1; j = E1[j].nxt) {
int a = belong[i], b = belong[E1[j].v];
if(a != b) {
addedge2(b, a);
du[a]++;
}
}
}
FORE(i, 1, n) {
int a = belong[i], b = belong[i + n];
conflict[a] = b;//縮點後有沖突的點
conflict[b] = a;
}
}
void topsort() {
int i,j;
queue<int> q;
FORE(i, 1, t)
if(du[i] == 0)
q.push(i);
cnt = 0;
while(!q.empty()) {
int x = q.front();
top[++cnt] = x;
q.pop();
for(i = head2[x]; i != -1; i = E2[i].nxt) {
int tmp = E2[i].v;
du[tmp]--;
if(du[tmp] == 0)
q.push(tmp);
}
}
}
void dfs_Blue(int u) {
int i;
color[u] = -1;
for(i = head2[u]; i != -1; i = E2[i].nxt) {
int x = E2[i].v;
if(color[x] == 0)
dfs_Blue(x);
}
}
void dfs_Red() {
int i, j;
FORE(i, 1, cnt) {
int x = top[i];
if(color[x] == 0) {
color[x] = 1;//Red
dfs_Blue(conflict[x]);//把所有與x有沖突的點及其子孫染藍色
}
}
}
void output() {
int i, j, t1, t2;
printf("YES\n");
FORE(i, 1, n) {
int x = belong[i], y = belong[i + n];
if(color[x] == 1) {
t1 = c[i].s;
t2 = c[i].s + c[i].len;
}
if(color[y] == 1) {
t1 = c[i].e - c[i].len;
t2 = c[i].e;
}
printf("%02d:%02d %02d:%02d\n", t1 / 60, t1 % 60, t2 / 60, t2 % 60);
}
}
void solve() {
int i, j;
FORE(i, 1, 2 * n)
if(!vis[i])
tarjan(i);
FORE(i, 1, n)
if(belong[i] == belong[i + n]) {
printf("NO\n");
return ;
}
Rebuild();
topsort();
dfs_Red();
output();
}
int main() {
while(scanf("%d", &n) != EOF) {
int i, j;
FORE(i, 1, n) {
int h1, m1, h2, m2;
scanf("%d:%d %d:%d %d", &h1, &m1, &h2, &m2, &c[i].len);
c[i].s = h1 * 60 + m1;
c[i].e = h2 * 60 + m2;
}
Build();
solve();
}
return 0;
}