POJ 2296 Map Labeler(2-sat)

题意：在一个地图上贴正方形标签，标签不能重合，地图上点的位置在标签的下边中点，或者上边中点，求标签的最大边长。

思路：一个标签可以在点的上方或者下方。设在上方的为i+n，下方为i。。二分边长，设边长为L  则：

当两点的x差值小于L时： 

1.两个y值一样时，一个标签在上方一个标签在下方，i -> -j     j -> -I   -i -> j   -j -> i

2.i点的y值大于j时，a.当两点y的差值小于L时：i的标签只能在上j的标签只能 在下： i -> -i    -j -> j

   b当两点y的差值小于2*L时，i的标签在下则j的标签只能在下： i -> j   -j  -> -i

3.j点的y值大于i时。同上：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
using namespace std;
const int INF = 0x3f3f3f3f;
struct cpoint{
    int x,y;
    void get(){scanf("%d%d",&x,&y);}
} re[109];
struct LT{
    int to,nex;
} L[100009];
int F[209],cnt;
void add(int f,int t)
{
    L[cnt].nex = F[f];
    L[cnt].to = t;
    F[f] = cnt++;
}
int n;
int ABS(int f){return f>=0?f:-f;}
int dfn[209],low[209],col[209],post[209],color,ind;
stack<int> S;
void tdfs(int k)
{
    dfn[k] = low[k] = ind++;
    S.push(k);post[k] = true;
    for(int i=F[k];i;i=L[i].nex)
    {
        int to = L[i].to;
        if(!dfn[to])
        {
            tdfs(to);
            low[k] = min(low[k],low[to]);
        }else if(post[to]&&dfn[to]<low[k])
        low[k] = dfn[to];
    }
    if(dfn[k]==low[k])
    {
        int i;color++;
        for(i=S.top(),S.pop();i!=k;i=S.top(),S.pop())
        {
            post[i] = 0,col[i] = color;
        }
        post[k] = 0,col[k] = color;
    }
}
bool solve()
{
    ind = 1;
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    for(int i=0;i<(n<<1);i++)
    if(!dfn[i]) tdfs(i);
    for(int i=0;i<n;i++)
    if(col[i]==col[i+n]) return false;
    return true;
}
int main()
{
    freopen("in.txt","r",stdin);
    int cas;
    scanf("%d",&cas);
    while(cas--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++) re[i].get();
        int l = 0,r = INF,mid;
        int ans = -1;
        while(l<=r)
        {
            mid = (l+r)>>1;
            cnt = 1;
            memset(F,0,sizeof(F));
            for(int i=0;i<n;i++)
            {
                for(int j=0;j<n;j++)
                {
                    if(i!=j&&ABS(re[i].x-re[j].x)<mid)
                    {
                        if(re[i].y==re[j].y)
                            add(i,j+n),add(j,i+n),add(j+n,i),add(i+n,j);
                        else if(re[i].y>re[j].y&&(re[i].y-re[j].y<mid))
                            add(i,i+n),add(j+n,j);
                        else if(re[j].y>re[i].y&&(re[j].y-re[i].y<mid))
                            add(j,j+n),add(i+n,i);
                        else if(re[i].y>re[j].y&&(re[i].y-re[j].y)<2*mid)
                            add(i,j),add(j+n,i+n);
                        else if(re[j].y>re[i].y&&(re[j].y-re[i].y)<2*mid)
                            add(j,i),add(i+n,j+n);
                    }
                }
            }
            if(solve())
                ans = mid,l = mid+1;
            else r = mid-1;
        }
        printf("%d\n",ans);
    }
    return 0;
}