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;
}