bzoj 3007: 拯救小云公主 最小生成树

题意

英雄又即将踏上拯救公主的道路……
这次的拯救目标是——爱和正义的小云公主。
英雄来到boss的洞穴门口，他一下子就懵了，因为面前不只是一只boss，而是上千只boss。当英雄意识到自己还是等级1的时候，他明白这就是一个不可能完成的任务。
但他不死心，他在想，能不能避开boss去拯救公主呢，嘻嘻。
Boss的洞穴可以看成一个矩形，英雄在左下角（1,1），公主在右上角（row，line）。英雄为了避开boss，当然是离boss距离越远越好了，所以英雄决定找一条路径使到距离boss的最短距离最远。
Ps:英雄走的方向是任意的。
你可以帮帮他吗？
当英雄找到了美丽漂亮的小云公主，立刻就被boss包围了！！！英雄缓闭双眼，举手轻挥，白光一闪后使用了回城卷轴，回到了城堡，但只有小云公主回去了……因为英雄忘了进入回城的法阵了。
n<=3000；

分析

看完题的第一个想法就是二分答案，但其实这题不用。
题解

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;

typedef long long LL;

const int N=3015;

int n,c,r,cnt,f[N],x[N],y[N];
struct edge{int x,y;LL len;}e[N*N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

LL get_dis(int i,int j)
{
    return (LL)(x[i]-x[j])*(x[i]-x[j])+(LL)(y[i]-y[j])*(y[i]-y[j]);
}

LL sqr(int x)
{
    return (LL)x*x;
}

bool operator < (edge a,edge b)
{
    return a.len<b.len;
}

int find(int x)
{
    if (f[x]==x) return x;
    f[x]=find(f[x]);
    return f[x];
}

int main()
{
    n=read();r=read();c=read();
    for (int i=1;i<=n;i++) x[i]=read(),y[i]=read();
    for (int i=1;i<n;i++)
        for (int j=i+1;j<=n;j++)
        {
            e[++cnt].x=i;e[cnt].y=j;e[cnt].len=get_dis(i,j);
        }
    for (int i=1;i<=n;i++)
    {
        e[++cnt].x=i;e[cnt].y=n+1;e[cnt].len=sqr(min(y[i]-1,r-x[i])*2);
        e[++cnt].x=i;e[cnt].y=n+2;e[cnt].len=sqr(min(x[i]-1,c-y[i])*2);
    }
    sort(e+1,e+cnt+1);
    LL mx=0;
    for (int i=1;i<=n+2;i++) f[i]=i;
    for (int i=1;i<=cnt;i++)
    {
        int p=find(e[i].x),q=find(e[i].y);
        if (p!=q)
        {
            f[p]=q;
            mx=max(mx,e[i].len);
            if (find(n+1)==find(n+2)) break;
        }
    }
    printf("%.2lf",(double)sqrt(mx)*0.5);
    return 0;
}