poj 2296 Map Labeler（2-SAT+二分,构图）

http://poj.org/problem?id=2296
题意：

给你n个点，要你在这n个点上放一个正方形，点
只能在正方形的上边或下边的中点上，所有正方形大小一样，
不能重叠，求最大的正方形。。。

分析：

如果abs(s[i].x-s[j].x)>=r则可以随便放
如果 abs[s[i].x-s[j].y)< r;
如果abs(s[i].y-s[j].y)< r,如果s[i].y==s[i].y则要求一个放上面一个放下面。
否则只能是上面的点放上面，下面的点放下面。
如果r<=abs(s[i].y-s[j].y)< 2*r,则除了上面的点放下方、下面的点放上方的情况都是可以的。

using namespace std;const int N = 1000 + 10;vector<int> g[N];int pre[N], lowlink[N], sccno[N], dfs_clock, scc_cnt;stack<int> S;void dfs(int u) {    pre[u] = lowlink[u] = ++dfs_clock;    S.push(u);    for(int i = 0; i < g[u].size(); i++) {        int v = g[u][i];        if(!pre[v]) {            dfs(v);            lowlink[u] = min(lowlink[u], lowlink[v]);        } else if(!sccno[v]) {            lowlink[u] = min(lowlink[u], pre[v]);        }    }    if(lowlink[u] == pre[u]) {        ++scc_cnt;        int x;        do {            x = S.top();            S.pop();            sccno[x] = scc_cnt;        } while(x!=u);    }}void find_scc(int n) {    dfs_clock = scc_cnt = 0;    memset(sccno, 0, sizeof(sccno));    memset(pre, 0, sizeof(pre));    for(int i = 0; i < n; i++)        if(!pre[i]) dfs(i);}//2*i上面//2*i+1 下面int a[N],b[N],n;bool ok(int d) {    //cout<<d<<endl;    for(int i=0; i<=2*n; i++)g[i].clear();    for(int i=0; i<n; i++)        for(int j=i+1; j<n; j++) {            if(abs(a[i]-a[j])>=d)continue;            if(abs(b[i]-b[j])<d) {                if(b[i]==b[j]) { //一个在上面一个在下面                    g[2*i].push_back(2*j+1);                    g[2*j].push_back(2*i+1);                    g[2*i+1].push_back(2*j);                    g[2*j+1].push_back(2*i);                } else if(b[i]>b[j]) { //i在上面 j在下面                    g[2*i+1].push_back(2*i);                    g[2*j].push_back(2*j+1);                } else {                    g[2*i].push_back(2*i+1);                    g[2*j+1].push_back(2*j);                }            } else if(abs(b[i]-b[j])<2*d) {                if(b[i]>b[j]) {  //i在下j在上不允许                    g[2*i+1].push_back(2*j+1);                    g[2*j].push_back(2*i);                } else { //                     g[2*j+1].push_back(2*i+1);                    g[2*i].push_back(2*j);                }            }        }    find_scc(n*2);    for(int i=0; i<2*n; i++)        if(sccno[i]==sccno[i^1]) {            return 0;        }    return 1;}int main() {    int T;   // freopen("f.txt","r",stdin);    scanf("%d",&T);    while(T--) {        scanf("%d",&n);        for(int i=0; i<n; i++)            scanf("%d%d",&a[i],&b[i]);        int l=0,r=100000;        while(l<r) {            int mid=l+(r-l+1)/2;            if(ok(mid))l=mid;            else r=mid-1;        }        printf("%d\n",l);    }    return 0;}