POJ2296.Map Labeler——2-sat二分最大值+可行性判断

http://poj.org/problem?id=2296

题目描述：
在二维坐标在有n个点，在这n个点上放正方形，每个点只能放在正方形上/下边界的中点，要使得这n个点放置的正方形互不相交，求正方形的最大边长

分析：
每个正方形要么放上面，要么放下面，2种限制。。。
对第i个点，2*i表示放上面，2*i+1表示放下面.设边长为r
1.当abs(x[i]-x[j])>=r，则这两个正方形可以任意放
2.当abs(y[i]-y[j]) < r
2.1当y[i]==y[j]时，i和j的放法相反 i xor j=1
2.2当y[i]>y[j]时，i必须取上，j必须取下
2.3当y[i] < y[j]时，i必须取下，j必须取上
3.当r< =abs(y[i]-y[j])<2*r时
3.1当y[i]>y[j]时，如果i放下，j一定也得放下
3.2当y[i]< y[j]时，如果i放上，j一定也得放上

构图贴模板

//320K  16MS    C++
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
const int MAXN=210;
const int MAXM=40010;
using namespace std;
int low[MAXN],dfn[MAXN],Stack[MAXN],sccno[MAXN];
int dfs_clock,scc,top;
int x[MAXN],y[MAXN];
bool instack[MAXN];
struct Edge{
    int to,next;
}edge[MAXM<<2];int head[MAXN],tot;
void addedge(int u,int v){
    edge[tot].to=v;edge[tot].next=head[u];head[u]=tot++;
}
void init(){
    tot=0;
    memset(head,0xff,sizeof(head));
}
void dfs(int u){
    int v;
    low[u]=dfn[u]=++dfs_clock;
    Stack[top++]=u;
    instack[u]=1;
    for(int i=head[u];i!=-1;i=edge[i].next){
        v=edge[i].to;
        if(!dfn[v]){
            dfs(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;
        for(;;){
            v=Stack[--top];
            instack[v]=0;
            sccno[v]=scc;
            if(v==u) break;
        }
    }

}
bool solveable(int n){
    memset(dfn,0,sizeof(dfn));
    memset(instack,false,sizeof(instack));
    dfs_clock=top=scc=0;
    for(int i=0;i<n;++i){
        if(!dfn[i]) dfs(i);
    }
    for(int i=0;i<n;i+=2){
        if(sccno[i]==sccno[i+1])
            return false;
    }
    return true;
}
void build_graph(int n,int mid){
    init();
    for(int i=0;i<n;++i){
        for(int j=i+1;j<n;++j){
            if(abs(x[i]-x[j])>=mid) continue;
            if(abs(y[i]-y[j])<mid){
                if(y[i]==y[j]){
                    addedge(i<<1,j<<1|1);
                    addedge(j<<1,i<<1|1);
                    addedge(i<<1|1,j<<1);
                    addedge(j<<1|1,i<<1);
                }
                if(y[i]>y[j]){
                    addedge(i<<1|1,i<<1);
                    addedge(j<<1,j<<1|1);
                }
                if(y[i]<y[j]){
                    addedge(i<<1,i<<1|1);
                    addedge(j<<1|1,j<<1);
                }
            }
            else if(abs(y[i]-y[j])<2*mid){
                if(y[i]>y[j]){
                    addedge(i<<1|1,j<<1|1);
                    addedge(j<<1,i<<1);
                }
                else{
                    addedge(i<<1,j<<1);
                    addedge(j<<1|1,i<<1|1);
                }
            }
        }
    }
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("in.cpp","r",stdin);
#endif // ONLINE_JUDGE
    int T,n,u,v;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i=0;i<n;++i){
            scanf("%d%d",&x[i],&y[i]);
        }
        int l=0,r=20000,mid,ans;
        while(l<=r){
            mid=(l+r)>>1;
            build_graph(n,mid);
            if(solveable(2*n)) {ans=mid,l=mid+1;}
            else r=mid-1;
        }
        printf("%d\n",ans);
    }
    return 0;
}