POJ 2296 Map Labeler

二分答案 + 2-SAT验证，判断正方形是否相交写起来有点烦，思路还是挺简单的。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<stack>
#include<queue>
#include<algorithm>
using namespace std;

const int maxn=;
struct Point
{
    double x,y;
}p[maxn];
int N;

stack<int>S;
vector<int>G[maxn];
vector<int>FG[maxn];
int Belong[maxn];
int flag[maxn];
int Block;

int left,right,mid,ans;

void init()
{
    for(int i=; i<maxn; i++) G[i].clear();
    for(int i=; i<maxn; i++) FG[i].clear();
    memset(Belong,,sizeof Belong);
    memset(flag,,sizeof flag);
    while(!S.empty()) S.pop();
    Block=;
}

void addEdge(int x,int y)
{
    G[x].push_back(y);
    FG[y].push_back(x);
}

void dfs1(int now)
{
    flag[now]=;
    for(int i=; i<G[now].size(); i++)
        if(!flag[G[now][i]])
            dfs1(G[now][i]);
    S.push(now);
}

void dfs2(int now)
{
    Belong[now]=Block;
    for(int i=; i<FG[now].size(); i++)
        if(!Belong[FG[now][i]])
            dfs2(FG[now][i]);
}

bool judge()
{
    for(int i=; i<*N; i++) if(!flag[i]) dfs1(i);
    while(!S.empty())
    {
        int Top=S.top();
        S.pop();
        if(!Belong[Top])
        {
            Block++;
            dfs2(Top);
        }
    }
    for(int i=; i<N; i++)
        if(Belong[*i]==Belong[*i+])
            return ;
    return ;
}

bool F(int p1,int p2,int f1,int f2)
{
    double leftTop_p1_x,leftTop_p1_y;
    double rightBottom_p1_x,rightBottom_p1_y;

    double leftTop_p2_x,leftTop_p2_y;
    double rightBottom_p2_x,rightBottom_p2_y;

    if(f1==)
    {
        leftTop_p1_x=p[p1].x-1.0*mid/;
        leftTop_p1_y=p[p1].y;

        rightBottom_p1_x=p[p1].x+1.0*mid/;
        rightBottom_p1_y=p[p1].y-mid;
    }
    else
    {
        leftTop_p1_x=p[p1].x-1.0*mid/;
        leftTop_p1_y=p[p1].y+mid;

        rightBottom_p1_x=p[p1].x+1.0*mid/;
        rightBottom_p1_y=p[p1].y;
    }

    if(f2==)
    {
        leftTop_p2_x=p[p2].x-1.0*mid/;
        leftTop_p2_y=p[p2].y;

        rightBottom_p2_x=p[p2].x+1.0*mid/;
        rightBottom_p2_y=p[p2].y-mid;
    }

    else
    {
        leftTop_p2_x=p[p2].x-1.0*mid/;
        leftTop_p2_y=p[p2].y+mid;

        rightBottom_p2_x=p[p2].x+1.0*mid/;
        rightBottom_p2_y=p[p2].y;
    }

    if(leftTop_p1_x<leftTop_p2_x&&leftTop_p2_x<rightBottom_p1_x
        && rightBottom_p1_y<leftTop_p2_y&&leftTop_p2_y<leftTop_p1_y) return ;

    if(leftTop_p1_x<leftTop_p2_x&&leftTop_p2_x<rightBottom_p1_x
        && rightBottom_p1_y<rightBottom_p2_y&&rightBottom_p2_y<leftTop_p1_y) return ;

    if(leftTop_p1_x<rightBottom_p2_x&&rightBottom_p2_x<rightBottom_p1_x
        && rightBottom_p1_y<leftTop_p2_y&&leftTop_p2_y<leftTop_p1_y) return ;

    if(leftTop_p1_x<rightBottom_p2_x&&rightBottom_p2_x<rightBottom_p1_x
        && rightBottom_p1_y<rightBottom_p2_y&&rightBottom_p2_y<leftTop_p1_y) return ;

    if(leftTop_p2_x<leftTop_p1_x&&leftTop_p1_x<rightBottom_p2_x
        && rightBottom_p2_y<leftTop_p1_y&&leftTop_p1_y<leftTop_p2_y) return ;

    if(leftTop_p2_x<leftTop_p1_x&&leftTop_p1_x<rightBottom_p2_x
        && rightBottom_p2_y<rightBottom_p1_y&&rightBottom_p1_y<leftTop_p2_y) return ;

    if(leftTop_p2_x<rightBottom_p1_x&&rightBottom_p1_x<rightBottom_p2_x
        && rightBottom_p2_y<leftTop_p1_y&&leftTop_p1_y<leftTop_p2_y) return ;

    if(leftTop_p2_x<rightBottom_p1_x&&rightBottom_p1_x<rightBottom_p2_x
        && rightBottom_p2_y<rightBottom_p1_y&&rightBottom_p1_y<leftTop_p2_y) return ;

    if(leftTop_p1_x==leftTop_p2_x&&rightBottom_p1_x==rightBottom_p2_x)
    {
        if(rightBottom_p1_y<leftTop_p2_y&&leftTop_p2_y<leftTop_p1_y) return ;
        if(rightBottom_p1_y<rightBottom_p2_y&&rightBottom_p2_y<leftTop_p1_y) return ;

        if(rightBottom_p2_y<leftTop_p1_y&&leftTop_p1_y<leftTop_p2_y) return ;
        if(rightBottom_p2_y<rightBottom_p1_y&&rightBottom_p1_y<leftTop_p2_y) return ;
    }

    if(leftTop_p1_y==leftTop_p2_y&&rightBottom_p1_y==rightBottom_p2_y)
    {
        if(leftTop_p1_x<leftTop_p2_x&&leftTop_p2_x<rightBottom_p1_x) return ;
        if(leftTop_p1_x<rightBottom_p2_x&&rightBottom_p2_x<rightBottom_p1_x) return ;

        if(leftTop_p2_x<leftTop_p1_x&&leftTop_p1_x<rightBottom_p2_x) return ;
        if(leftTop_p2_x<rightBottom_p1_x&&rightBottom_p1_x<rightBottom_p2_x) return ;
    }

    if(leftTop_p1_x==leftTop_p2_x&&rightBottom_p1_x==rightBottom_p2_x)
        if(leftTop_p1_y==leftTop_p2_y&&rightBottom_p1_y==rightBottom_p2_y)
            return ;

    return ;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&N);
        for(int i=;i<N;i++) scanf("%lf%lf",&p[i].x,&p[i].y);

        left=;right=;
        while(left<=right)
        {
            mid=(left+right)/;
            init();
            for(int i=;i<N;i++)
            {
                for(int j=i+;j<N;j++)
                {
                    if(F(i,j,,)) {addEdge(*i,*j+);addEdge(*j,*i+);}
                    if(F(i,j,,)) {addEdge(*i,*j);addEdge(*j+,*i+);}
                    if(F(i,j,,)) {addEdge(*i+,*j+);addEdge(*j,*i);}
                    if(F(i,j,,)) {addEdge(*i+,*j);addEdge(*j+,*i);}
                }
            }
            if(judge()) {ans=mid;left=mid+;}
            else right=mid-;
        }
        printf("%d\n",ans);
    }
    return ;
}