POJ 3207 Ikki's Story IV - Panda's Trick（2 - sat啊）

题目链接：http://poj.org/problem?id=3207

Description

liympanda, one of Ikki’s friend, likes playing games with Ikki. Today after minesweeping with Ikki and winning so many times, he is tired of such easy games and wants to play another game with Ikki.

liympanda has a magic circle and he puts it on a plane, there are n points on its boundary in circular border: 0, 1, 2, …, n − 1. Evil panda claims that he is connecting m pairs of points. To connect two points, liympanda either places the link entirely inside the circle or entirely outside the circle. Now liympanda tells Ikki no two links touch inside/outside the circle, except on the boundary. He wants Ikki to figure out whether this is possible…

Despaired at the minesweeping game just played, Ikki is totally at a loss, so he decides to write a program to help him.

Input

The input contains exactly one test case.

In the test case there will be a line consisting of of two integers: n and m (n ≤ 1,000, m ≤ 500). The following m lines each contain two integers a_i and b_i, which denote the endpoints of the ith wire. Every point will have at most one link.

Output

Output a line, either “panda is telling the truth...” or “the evil panda is lying again”.

Sample Input

4 20 13 2

Sample Output

panda is telling the truth...

Source

POJ Monthly--2007.03.04, Ikki

题意：

平面上有一个圆，圆的边上按顺时针放着0..n-1共n个点。现在要连m条边，比如a，b，那么a到b可以从圆的内部连接，也可以从圆的外部连接。给你的信息中，每个点最多只能连一条边。问是否可以连接这m条边，使这些边都不相交。

PS:

对于每条Link，要么在圆外，要么在圆内，且不可同时满足，
只能两者取一，判断这M条Link是否合法，也就是M条Link不冲突，
这就是典型的2-sat问题了。 将每条Link i 看做一个点，如果Link在圆内，
则选做i ，如果在圆外， 则选做i'。对于两条线(i,j) ，如果i，j不能同时
在圆内，也就可以推出两者不能同时在圆外，这个证明很容易，读者可
以自行证明。i, j不能同时在圆内，则有边(i, j') 、(j ,i')、(i'，j)、(j' ,i)
(这是由2-sat的构图方式决定的，具体可以看《由对称性解2-SAT问题》
 这篇论文）。建图完了之后，本题就是判断2-sat问题是否有解，
 先求原图的强连通分量，并缩点，(这里我们称：(i,i')属于同一组)，
 判断是否存在(i,i')属于同一组，若存在，则不可能，若不存在则可能。

代码如下：

#include <cstdio>#include <cstring>#include <queue>#include <vector>#include <algorithm>#include <iostream>using namespace std;//2-SAT 强连通缩点const int MAXN = 2010;//点数的两倍const int MAXM = 2000010;struct Edge{    int to;    int next;} edge[MAXM];int head[MAXN],tot;void init(){    tot = 0;    memset(head,-1,sizeof(head));}void addedge(int u,int v){    edge[tot].to = v;    edge[tot].next = head[u];    head[u] = tot++;}int Low[MAXN],DFN[MAXN],Stack[MAXN],Belong[MAXN];//Belong数组的值1~sccint Index,top;int scc;bool Instack[MAXN];int num[MAXN];void Tarjan(int u)//tarjan求强连通分量{    int v;    Low[u] = DFN[u] = ++Index;    Stack[top++] = u;    Instack[u] = true;    for(int i = head[u]; i != -1; i = edge[i].next)    {        v = edge[i].to;        if( !DFN[v] )        {            Tarjan(v);            if(Low[u] > Low[v])Low[u] = Low[v];        }        else if(Instack[v] && Low[u] > DFN[v])            Low[u] = DFN[v];    }    if(Low[u] == DFN[u])    {        scc++;        do        {            v = Stack[--top];            Instack[v] = false;            Belong[v] = scc;        }        while(v != u);    }}//判断是否有解bool solvable(int n)//n是总个数,需要选择一半{    memset(DFN,0,sizeof(DFN));    memset(Instack,false,sizeof(Instack));    Index = scc = top = 0;    for(int i = 1; i <= 2*n; i++)        if(!DFN[i])            Tarjan(i);    for(int i = 1; i <= n; i++)        if(Belong[2*i] == Belong[2*i+1]) return false;    return true;}int a[MAXN], b[MAXN];int n, m;void build_grap(){    for(int i = 1; i <= m; i++)    {        for(int j = i+1; j <= m; j++)        {            if((a[i]<=a[j]&&b[i]>=a[j]&&b[i]<=b[j]) || (a[i]>=a[j]&&a[i]<=b[j]&&b[i]>=b[j]))            {                addedge(2*i,2*j+1);//内 外                addedge(2*j,2*i+1);//内 外                addedge(2*i+1,2*j);//外 内                addedge(2*j+1,2*i);//外 内            }        }    }}int main(){    while(~scanf("%d%d",&n,&m))    {        if(n == 0 && m == 0)            break;        init();        for(int i = 1; i <= m; i++)        {            scanf("%d%d",&a[i],&b[i]);            if(a[i] > b[i])            {                int tt = a[i];                a[i] = b[i];                b[i] = tt;            }        }                build_grap();        if(solvable(n))//判断是否有解        {            printf("panda is telling the truth...\n");        }        else        {            printf("the evil panda is lying again\n");        }                }    return 0;}