POJ 2296 - Map Labeler 构图2-sat..注意细节...

              题意:

                       有一个地图...地图上有N个city...每个city是都是一个边长相同的平行于x/y轴的正方形..现在给了每个city的一个点.这个点可能是正方形上边的中点,也可能是下边的中点..每个city的区域不能有交集(边可以重合)..问正方形的边最长可以是多少..

              题解:

                       很明显的二分查找+2-sat判断....关键就是如何建图...先将所有city按照x从小到大排序..判断冲突的前提两者横坐标之差小于所二分的长度M..那么就只要考虑纵坐标了..分为两个大方向...两个点之间的距离小于M...两个点之间的距离小于2*M...而这里必须要注意一个特殊的情况..两个点的纵坐标相同.拿出来单独讨论..

Program:

#include<iostream>
#include<stdio.h>
#include<cmath>
#include<queue>
#include<stack>
#include<string.h>
#include<map>
#include<set>
#include<algorithm>
#define oo 1000000007
#define MAXN 20005<<1
#define MAXM 100000<<2
#define ll long long
using namespace std;  
struct city
{
       int x,y;
       bool operator <(city a) const
       {
               return x<a.x;
       }
}C[MAXN];
struct node
{
       int y,next;
}line[MAXM];
int Lnum,_next[MAXN],dfn[MAXN],low[MAXN],tp[MAXN],tpnum,DfsIndex;
bool instack[MAXN];
stack<int> mystack;
void addline(int x,int y)
{
       line[++Lnum].next=_next[x],_next[x]=Lnum,line[Lnum].y=y;
}
void tarjan(int x)
{
       instack[x]=true,mystack.push(x);
       dfn[x]=low[x]=++DfsIndex;
       for (int k=_next[x];k;k=line[k].next)
       {
               int y=line[k].y;
               if (!dfn[y])
               {
                     tarjan(y);
                     low[x]=min(low[x],low[y]);
               }else
               if (instack[y])
                     low[x]=min(low[x],dfn[y]);
       }
       if (low[x]==dfn[x])
       {
               tpnum++;
               do
               {
                       x=mystack.top();
                       mystack.pop();
                       instack[x]=false;
                       tp[x]=tpnum;
               }while (low[x]!=dfn[x]);
       }
}
bool _2sat(int N,int M)
{
       int i,j;
       Lnum=0,memset(_next,0,sizeof(_next));
       for (i=0;i<N;i++)  // i<<1上方  i<<1|1 下方 
          for (j=i+1;j<N;j++)
             if (C[j].x-C[i].x<M && abs(C[j].y-C[i].y)<(M<<1))
             {
                    if (abs(C[j].y-C[i].y)<M)
                    {
                           if (C[j].y==C[i].y)
                           {
                                 addline(i<<1,j<<1|1),addline(i<<1|1,j<<1);
                                 addline(j<<1,i<<1|1),addline(j<<1|1,i<<1);
                           }else
                           if (C[j].y>C[i].y) addline(j<<1|1,j<<1),addline(i<<1,i<<1|1);
                                        else  addline(i<<1|1,i<<1),addline(j<<1,j<<1|1);
                    }else
                    if (abs(C[j].y-C[i].y)<(M<<1))
                    {
                           if (C[j].y>C[i].y) addline(j<<1|1,i<<1|1),addline(i<<1,j<<1);
                                        else  addline(i<<1|1,j<<1|1),addline(j<<1,i<<1);
                    }
             }
       memset(dfn,0,sizeof(dfn));
       memset(instack,false,sizeof(instack));
       while (!mystack.empty()) mystack.pop();
       DfsIndex=tpnum=0;
       for (i=0;i<(N<<1);i++)
           if (!dfn[i]) tarjan(i);
       for (i=0;i<N;i++)
          if (tp[i<<1]==tp[i<<1|1]) return false;       
       return true;
}
int main()
{          
       int cases,N,i; 
       scanf("%d",&cases);
       while (cases--)
       {
               scanf("%d",&N);
               for (i=0;i<N;i++) scanf("%d%d",&C[i].x,&C[i].y);
               sort(C,C+N);
               int l=-1,r=oo,mid;
               while (r-l>1)
               {
                      mid=l+r>>1;
                      if (_2sat(N,mid)) l=mid;
                           else r=mid;
               }
               printf("%d\n",l);
       }
       return 0;
}