pku 3207 Ikki's Story IV-Panda's Trick 2-sat実行可能解の有無判定

http://poj.org/problem?id=3207
タイトル:
1つの円盤の縁にn個の点があり、下付き文字は0から、m本の線が2 m個の互いに異なる点を接続し、線は円盤の内部にあってもよく、円盤の外部にあってもよく、任意の2本の線が交差できないことを要求する.m本の線(内外随意)を与え、各線が交差しないかどうかを聞く.
考え方:
i番目の線は、2*iと2*i+1のペアの頂点と見なすことができ、線分iとjが交差すると、2*iと2*j+1、および2*i+1と2*jの間に双方向のエッジが接続される.次に2−satに変換して実行可能解が存在するか否かを判断する.

//#pragma comment(linker,"/STACK:327680000,327680000")

#include <iostream>

#include <cstdio>

#include <cmath>

#include <vector>

#include <cstring>

#include <algorithm>

#include <string>

#include <set>

#include <functional>

#include <numeric>

#include <sstream>

#include <stack>

#include <map>

#include <queue>



#define CL(arr, val)    memset(arr, val, sizeof(arr))



#define ll long long

#define inf 0x7f7f7f7f

#define lc l,m,rt<<1

#define rc m + 1,r,rt<<1|1

#define pi acos(-1.0)

#define ll long long

#define L(x)    (x) << 1

#define R(x)    (x) << 1 | 1

#define MID(l, r)   (l + r) >> 1

#define Min(x, y)   (x) < (y) ? (x) : (y)

#define Max(x, y)   (x) < (y) ? (y) : (x)

#define E(x)        (1 << (x))

#define iabs(x)     (x) < 0 ? -(x) : (x)

#define OUT(x)  printf("%I64d
", x)

#define lowbit(x)   (x)&(-x)

#define Read()  freopen("din.txt", "r", stdin)

#define Write() freopen("dout.txt", "w", stdout);





#define N 1007

#define M 1007

using namespace std;



struct side

{

    int s,e;

}sid[M];



int n,m;

struct node

{

    int v;

    int next;

}g[M*M];

int head[M],ct;



int dfn[M],low[M];

int belong[M],stk[M];

bool isn[M];

int idx,cnt,top;







bool isok(int i,int j)

{

    //    ，     

    if ((sid[i].s < sid[j].s && sid[i].e > sid[j].s && sid[i].e < sid[j].e)

        || (sid[i].s > sid[j].s && sid[i].s < sid[j].e && sid[i].e > sid[j].e))

        return true;

    else return false;

}

void add(int u,int v)

{

    g[ct].v = v;

    g[ct].next = head[u];

    head[u] = ct++;



    g[ct].v = u;

    g[ct].next = head[v];

    head[v] = ct++;

}

void tarjan(int u)

{

    int i,j;

    dfn[u] = low[u] = ++idx;

    stk[++top] = u;

    isn[u] = true;

    for (i = head[u]; i != - 1; i = g[i].next)

    {

        int v = g[i].v;

        if (dfn[v] == -1)

        {

            tarjan(v);

            low[u] = min(low[u],low[v]);

        }

        else if (isn[v])

        {

            low[u] = min(low[u],dfn[v]);

        }

    }

    if (dfn[u] == low[u])

    {

        cnt++;

        do

        {

            j = stk[top--];

            belong[j] = cnt;

            isn[j] = false;

        }while (j != u);

    }

}



void solve()

{

    int i;

    for (i = 0; i < 2*m; ++i)

    {

        dfn[i] = low[i] = -1;

        isn[i] = false;

        belong[i] = 0;

    }

    idx = cnt = top = 0;



    for (i = 0; i < 2*m; ++i)

    {

        if (dfn[i] == -1) tarjan(i);

    }

    bool flag = false;

    for (i = 0; i < m; ++i)

    {

        if (belong[2*i] == belong[2*i + 1])

        {

            flag = true;

            break;

        }

    }

    if (flag) printf("the evil panda is lying again
");

    else printf("panda is telling the truth...
");



}

int main()

{

   // Read();

    int i,j;

    while (~scanf("%d%d",&n,&m))

    {

        for (i = 0; i < m; ++i)

        {

            scanf("%d%d",&sid[i].s,&sid[i].e);

            if (sid[i].s > sid[i].e)

            swap(sid[i].s,sid[i].e);

        }

        CL(head,-1); ct = 0;

        for (i = 0; i < m; ++i)

        {

            for (j = i + 1; j < m; ++j)

            {

                if (isok(i,j))

                {

                    add(2*i,2*j + 1);

                    add(2*j,2*i + 1);

                }

            }

        }

        solve();

    }

    return 0;

}