poj 2296(二分+2-SAT)

poj 2296
教训：
1、注意复杂情况的分类讨论
2、$a<x<b$这样的不要写错了

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
#define ms(i,j) memset(i, j, sizeof i);

const int MAXM = 1000 + 5; 

struct TwoSAT
{
    vector<int> G[MAXM*2];
    int S[MAXM*2];
    bool mark[MAXM*2];
    int n, c;
    void init(int n)
    {
        this->n = n;
        for (int i=0;i<=2*n;i++)
        {
            G[i].clear();
            mark[i] = false;
        }
    }
    void ins(int u, int v)
    {
        G[u].push_back(v);
    }
    int dfs(int x)
    {
        if (mark[x^1]) return false;
        if (mark[x]) return true;
        mark[x] = true;
        S[++c] = x;
        for (int i=0;i<G[x].size();i++)
        {
            if (!dfs(G[x][i])) return false;
        }
        return true;
    }
    int solve()
    {
        for (int i=0;i<2*n;i+=2)
        if (!mark[i]&&!mark[i+1])
        {
            c = 0;
            if (!dfs(i))
            {
                while (c>0) mark[S[c--]] = false;
                if (!dfs(i+1)) return false;
            }
        }
        return true;
    }
}ts;

int m, xi[MAXM], yi[MAXM];

int abss(int x) {return (x<0) ? (-x) : (x);}
void init()
{
    scanf("%d", &m);
    for (int i=0;i<m;i++)
    {
        scanf("%d%d", &xi[i], &yi[i]);
    }
}
void solve()
{
    int l = 0, r = 80000, ans = 0;
    while (l < r)
    {
        int mid = (l+r)/2;
        ts.init(m);
        for (int i=0;i<m;i++)
        for (int j=i+1;j<m;j++)
        {
            int ax = abss(xi[i]-xi[j]);
            int ay = abss(yi[i]-yi[j]);
            if (ax<mid)//2*i is up, 2*i+1 is down
            {
                if (ay<mid)
                {
                    if (yi[i]==yi[j])//特判差值为0的情况 
                    {
                        ts.ins(2*i+1, 2*j);
                        ts.ins(2*j+1, 2*i);
                        ts.ins(2*i, 2*j+1);
                        ts.ins(2*j, 2*i+1);
                    } else if (yi[i]>yi[j])//A是上面的点，B是下面的点 
                    {
                        ts.ins(2*i+1, 2*i);
                        ts.ins(2*j, 2*j+1);
                    } else//B是上面的点，A是下面的点 
                    {
                        ts.ins(2*j+1, 2*j);
                        ts.ins(2*i, 2*i+1);
                    }
                } else if (ay<2*mid)
                {
                    if (yi[i]>yi[j])//A是上面的点，B是下面的点 
                    {
                        ts.ins(2*i+1, 2*j+1);
                        ts.ins(2*j, 2*i);
                    } else//B是上面的点，A是下面的点 
                    {
                        ts.ins(2*j+1, 2*i+1);
                        ts.ins(2*i, 2*j);
                    }
                }
            }
        }
        if (ts.solve())
        {
            ans = mid;
            l = mid+1;
        }else r = mid;
    }
    printf("%d\n", ans);
}
int main()
{
    int t;
    scanf("%d", &t);
    while (t--)
    {
        init();
        solve();
    }
    return 0;
}