poj 3207 Ikki's Story IV - Panda's Trick

類型：2-sat

題目：http://poj.org/problem?id=3207

來源：POJ Monthly--2007.03.04, Ikki

思路：對每一條link有兩種選擇，放圓內或者圓外。求滿足條件的情況是否存在，屬於2-sat問題。

（1）將每條link看做兩個點，分別表示在圓內鏈接和在圓外鏈接，只能選一個，用i和i`表示。

（2）如果兩條link必須一個在圓內，一個在圓外。那麼

若：選擇i，j`必選

        選擇j，i`必選

        選擇i`，j必選

        選擇j`，i必選

據此構造新邊，然後通過Tarjan算法求出所有強連通分量，如果某i和i`屬於同一強連通分量，結果爲假，反之爲真
參考：http://blog.sina.com.cn/s/blog_64675f540100k13v.html

// poj 3207 Ikki's Story IV - Panda's Trick
// ole 1656K 32MS
#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 MAXN = 1100;
const int MAXM = 500 * 500 * 4 + 2;
bool in_stack[MAXN];
int num, cnt, m, n, cnt_num;
int step[MAXN];
int low[MAXN];
int head[MAXN];
int belong[MAXN];
int g[MAXN][2];
stack<int> st;
struct edge{
    int v, nxt;
}e[MAXM];

void addedge(int u, int v) {
    e[cnt].v = v;
    e[cnt].nxt = head[u];
    head[u] = cnt++;
}

void Tarjan(int u) {
    int v;
    int i,j;

    step[u] = low[u] = ++num;
    st.push(u);
    in_stack[u] = true;
    for(i = head[u]; i != -1; i = e[i].nxt){
        v = e[i].v;
        if(!step[v]){
            Tarjan(v);
            low[u] = min(low[u],low[v]);
        }
        else if(in_stack[v])
            low[u] = min(low[u],step[v]);
    }
    if(step[u] == low[u]){
        cnt_num++;
        do{
            j = st.top();
            st.pop();
            in_stack[j] = false;
            belong[j] = cnt_num;
        }while(j != u);
    }
}

void init() {
    int i, j;

    CLR(head, -1);
    CLR(belong, 0);
    CLR(g, 0);
    CLR(step, 0);
    CLR(low, 0);
    CLR(in_stack, false);
    cnt = num = cnt_num = 0;
    FORE(i, 1, m) {
        scanf("%d %d", &g[i][0], &g[i][1]);
        if(g[i][0] > g[i][1])
            swap(g[i][0], g[i][1]);
    }
    FORE(i, 1, m - 1) FORE(j, i + 1, m) {
        if(g[i][0] < g[j][0] && g[j][0] < g[i][1] && g[i][1] < g[j][1] ||
           g[j][0] < g[i][0] && g[i][0] < g[j][1] && g[j][1] < g[i][1]) {
                addedge(i, j + m);
                addedge(i + m, j);
                addedge(j, i + m);
                addedge(j + m, i);
           }
    }
}

void solve() {
    int i;

    FORE(i, 1, 2 * m) if(step[i] == 0)
        Tarjan(i);
    FORE(i, 1, m)
        if(belong[i] == belong[i + m])
            break;
    (i <= m) ? printf("the evil panda is lying again\n") : printf("panda is telling the truth...\n");
}

int main() {
    while(scanf("%d %d", &n, &m) != EOF) {
        init();
        solve();
    }
    return 0;
}