【CSP模拟赛】starway(玄学建边 最小生成树)

问題描述

   小w伤心的走上了 Star way to heaven。
   到天堂的道路是一个笛卡尔坐标系上一个n×m的长方形通道(顶点在(0,0))和(n,m)),小w从最左边任意一点进入,从右边任意一点走到天堂。
   最左最右的距离为n,上下边界距离为m。
   其中长方形内有k个Star,每个Star都有一个整点坐标,Star的大小可以忽略不计,
   每个Star以及长方形上下两个边缘(宇宙的边界)都有引力,所以为了成功到达heaven小w离他们越远越好。
   请问小V走到终点的路径上,距离所有星星以及边界的最小距离最大值可以为多少?

输入格式

   一行三个整数n,m,k。
   接下来k行,每行两个整数xi,yi表示一个点的坐标。

输出格式

    一行一个整数表示答案,绝对误差不能超过10-6
输入样例
    10 5 2
    1 1
    2 3
输出样例
    1.11803399
    对于100%的数据 $kleq 6000 ,n,mleq 10^6$

分析

  没什么好分析的,noi.ac上的原题都不会做了,可见这个Oier有多菜

  根据套路,题目中说到最小距离最大,我们考虑二分这个最小距离然后检验

  设最小距离为len

  对于一个点,如果到它的最小距离为len,那么以它为圆心,len为半径形成的圆就是雷区,是不能进入的

  对于一条线,雷区就是一个贼长的矩形

  下图是len为1时样例的雷区

  显然,如果小w只走白色区域就可以从左边走到右边,那么这个len就是可以在增加一些的

  那么如何判断这些雷区可以封死小w走的路呢?

  我们来看一下封死的情况,让len为$sqrt 2$

 

  如果把上下边界也看成点,距离小于len*2的点连边

 

  我们发现,若是上下两个边界联通了,就从左边到右边的路就被封死了。

  于是我们可以二分len然后O(n^2)连边,O(n)判定,复杂度为O(nlogn)

  不过6000的数据还是过不了。

  由联通,所以想到最小生成树,在最小生成树中,从下边界到上边界路径上最长的边应该就是限制上下边界联通的边,拿它除以2就是答案

  然而身为年度憨憨人物的我忘了prim只记得kruskal,暴力最小生成树

  prim的时间复杂度为O(n^2),kruskal的时间复杂度为O(mlogm),此题m=n^2

  所以选择prim算法。

  如果怕精度问题可以用longlong存下距离的平方

#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
int n,m,k,vis[6005],pre[6005],p[6005][2];long long d[6005],ans;
long long dis(int a,int b){return 1ll*(p[a][0]-p[b][0])*(p[a][0]-p[b][0])+1ll*(p[a][1]-p[b][1])*(p[a][1]-p[b][1]);}
int main()
{
    scanf("%d%d%d",&n,&m,&k);d[k+2]=1ll*m*m;
    for(int i=1;i<=k;i++)scanf("%d%d",&p[i][0],&p[i][1]),d[i]=1ll*p[i][1]*p[i][1];
    for(int i=1;i<=k+2;i++)
    {
        int nx=0;long long mn=1000000000000ll;
        for(int j=1;j<=k+2;j++)if(d[j]<=mn&&!vis[j])mn=d[nx=j];
        if(nx==k+2){d[nx]=ans=mn;break;}vis[nx]=1;
        for(int j=1;j<=k+2;j++)
        if(j==k+2)
            pre[j]=d[j]>1ll*(m-p[nx][1])*(m-p[nx][1])?nx:pre[j],
            d[j]=min(d[j],1ll*(m-p[nx][1])*(m-p[nx][1]));
        else if(!vis[j]&&d[j]>dis(j,nx))d[j]=dis(j,nx),pre[j]=nx;
    }
    int x=k+2;while(x)ans=max(ans,d[x]),x=pre[x];
    printf("%.9lf
",sqrt(ans)/2.0);
}
原文地址:https://www.cnblogs.com/firecrazy/p/11772126.html