P2498 [SDOI2012]拯救小云公主
题目描述
英雄又即将踏上拯救公主的道路……
这次的拯救目标是——爱和正义的小云公主。
英雄来到\(boss\)的洞穴门口,他一下子就懵了,因为面前不只是一只\(boss\),而是上千只\(boss\)。当英雄意识到自己还是等级1的时候,他明白这就是一个不可能完成的任务。
但他不死心,他在想,能不能避开\(boss\)去拯救公主呢,嘻嘻。
\(Boss\)的洞穴可以看成一个矩形,英雄在左下角\((1,1)\),公主在右上角(row,line)\(。英雄为了避开\)boss\(,当然是离\)boss\(距离越远越好了,所以英雄决定找一条路径使到距离\)boss$的最短距离最远。
Ps:英雄走的方向是任意的。
你可以帮帮他吗?
当英雄找到了美丽漂亮的小云公主,立刻就被\(boss\)包围了!!!英雄缓闭双眼,举手轻挥,白光一闪后使用了回城卷轴,回到了城堡,但只有小云公主回去了……因为英雄忘了进入回城的法阵了。
输入输出格式
输入格式:
第一行,输入三个整数,\(n\)表示\(boss\)的数目,\(row\),\(line\)表示矩形的大小;
接下来\(n\)行,每行分别两个整数表示\(boss\)的位置坐标。
输出格式:
输出一个小数,表示英雄的路径离\(boss\)的最远距离,精确到小数点后两位。
感觉这题好神,没见过模型的话。
70pts:
二分答案最远距离,并查集维护boss所在圆是否连通。
复杂度:\(O(log?N^2)\)
期望得分:70 实际得分:70~100
100pts:
(1) 最小生成树
这样去看,英雄通过这个平面图,实际上只有两种类型。
一是穿过个两个boss之间的地方,这时我们走两个boss之间最中间的地方是最划算的。
二是穿过边界和boss之间的地方,这时我们贴着边界走是最划算的。
我们有很多个这样的地方可以走,走哪些我们可以确定既可以穿过去所有的boss又能使它距离我们最远呢?
如果我们把这些点(包括边界抽象出来的点)连成一条链
我们可以选择一个最优的方案,从中间穿过去
然而仔细想想,如果通过了另外的两点比最大距离还小的点,那不就废了吗?
换一种角度理解好了
按照\(kruskal\)的方法,我们一个一个从小到大加边
当上面的边和下面的边所代表的虚点(实际上有四种情况)连通时,我们就不得不对这条链(实际形态应该是树的连接两个虚点的链)上的某一边做出选择了,这时候我们选择之中最大的走,也不会遇到更小的两点距离了
注意如果是boss和Boss之间的边权要除2(走中间)
然而这个题只可以用prim算法(稠密图)
复杂度:\(O(N^2)\)
期望得分:100 实际得分:100
(2) SPFA
实质上SPFA也只是求得了这条在最小生成树上的最大边的权值。
我们跑0点到n+1点(两个虚点)之间的每条路径上最长的边其中的最小值
注意松弛目的是最小值,但松弛时是从前面的边的最小值和当前边选出一个最大的
复杂度:\(O(N^2)\) (建图)
期望得分:100 实际得分:100
贴一下SPFA的代码吧
#include
#include
#include
double min(double x,double y){return xy?x:y;}
const int N=3050;
double g[N][N],dis[N],x[N],y[N],n,m;
int num,used[N];
double get(double x1,double y1,double x2,double y2)
{
return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2))/2;
}
void init()
{
scanf("%d%lf%lf",&num,&n,&m);
for(int i=1;i<=num;i++)
scanf("%lf%lf",x+i,y+i);
for(int i=1;i<=num;i++)
for(int j=i+1;j<=num;j++)
g[i][j]=g[j][i]=get(x[i],y[i],x[j],y[j]);
for(int i=1;i<=num;i++)
{
g[0][i]=g[i][0]=min(x[i]-1,m-y[i]);
g[i][num+1]=g[num+1][i]=min((n-x[i]),y[i]-1);
}
g[0][num+1]=g[num+1][0]=1e233;
dis[0]=0;
for(int i=1;i<=num+1;i++)
dis[i]=1e233;
}
void spfa()
{
std::queue q;
q.push(0);
while(!q.empty())
{
int u=q.front();
q.pop();used[u]=0;
for(int i=1;i<=num+1;i++)
{
if(i==u) continue;
if(dis[i]>max(dis[u],g[u][i]))
{
dis[i]=max(dis[u],g[u][i]);
if(!used[i])
{
used[i]=1;
q.push(i);
}
}
}
}
}
int main()
{
init();
spfa();
printf("%.2lf",dis[num+1]);
return 0;
}
2018.7.12