【NOIP模拟考三】水资源——二分答案+Floyd验证 day1 first 信号连接

本文链接：https://blog.csdn.net/c20190101zjx/article/details/77163656

题目描述

输入

输出

样例输入

(如果复制到控制台无换行，可以先粘贴到文本编辑器，再复制)

样例输出

5.00000000

提示

n≤100

精度要求1e-6,建议保留8位输出

我们传说中的4次模拟赛终于已经结束了三场，接受了三位大神的的神题洗礼，感觉整个人都萌比了。。。不过题目也是状况频出，数据有问题、标程有问题、题目描述出现错误都有。。。好了不说了，我们先来看看这道水资源吧

写这道题只是顾及其他题要写，不得已而为之。显然，现在的半径合法，半径再加大一点也依然合法；反之，现在的半径不合法，那么小一点也一定不会合法，满足函数单调性，所以可以进行二分答案！由于最大只需要让最远的两个点能够触碰到即可，所以最大半径约为14142（sqrt(20000^2+20000^2)/2），由于要保留8位小数，还要乘上10^8再log，这样算出来是40左右。

然后考虑验证，我们可以将小于半径两倍的边加入图，用Floyd算出两点间最短距离，只要最长距离不超过T即可，这样做耗费是100^3，也就是100W，两个乘起来就是4000W，可以卡过。

后来看了题解就萌比了，为什么要枚举半径？显然，最优时一定会有两点距离恰好是两倍半径（这样可以保证最优），那么我们就只要枚举是哪一条边，边数最多有9900，log(9900)就是13左右，要快很多（不过上面的方法也可以A就是了）

好了，放代码了。。。

//%>_<% ~(>_<)~
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define inf 0x3f3f3f3f
#define N 105
#define eps 1e-9
#define sqr(a) ((a)*(a))
int n,T,x[N],y[N];
double d(int a,int b){return sqrt(sqr(x[a]-x[b])+sqr(y[a]-y[b]));}
double dis[N][N];
double A[N][N];
bool equal(double a,double b)
{
	if(a-b<eps&&b-a<eps)
		return 1;
	return 0;
}
bool check(double r)
{
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if(dis[i][j]<=2*r)
				A[i][j]=dis[i][j];
	for(int k=1;k<=n;k++)
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				A[i][j]=min(A[i][j],A[i][k]+A[k][j]);
	double maxx=0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			maxx=max(maxx,A[i][j]);
			A[i][j]=inf;
		}
	return maxx<=T;
}
int main()
{
	scanf("%d%d",&n,&T);
	for(int i=1;i<=n;i++)
		scanf("%d%d",&x[i],&y[i]);
	double l=0,r=0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			dis[i][j]=d(i,j);
			A[i][j]=inf;
			if(dis[i][j]-eps>T)
			{
				printf("-1");
				return 0;
			}
			r=max(r,dis[i][j]/2);
		}
	while(!equal(l,r))
	{
		double mid=(l+r)/2;
		if(check(mid))
			r=mid;
		else
			l=mid;
	}
	printf("%.8lf",l);
}