BZOJ 3007 拯救小云公主 二分答案+对偶图

标签：bzoj   bzoj3007   二分答案   对偶图   计算几何   

题目大意：给定一个矩形和矩形内的一些点，求一条左下角到右上角的路径，使所有点到这条路径的最小距离最大

最小距离最大，果断二分答案

现在问题转化成了给定矩形中的一些圆形障碍物求左下角和右上角是否连通

然后就是对偶图的问题了

左下角和右上角连通等价于对偶图中左上两条边和右下两条边不连通

因此将所有相交的圆之间连边，从左上两条边广搜即可

时间复杂度O(n^2log(min(r,l)/EPS))

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 3030
#define S 0
#define T (n+1)
using namespace std;
struct Point{
	double x,y;
	friend istream& operator >> (istream &_,Point &p)
	{
		scanf("%lf%lf",&p.x,&p.y);
		return _;
	}
	friend double Distance(const Point &p1,const Point &p2)
	{
		return sqrt( (p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y) );
	}
}pos,points[M];
struct abcd{
	int to,next;
}table[M*M];
int head[M],tot;
int n;
double dis[M][M];
void Initialize()
{
	memset(head,0,sizeof head);
	tot=0;
}
void Add(int x,int y)
{
	table[++tot].to=y;
	table[tot].next=head[x];
	head[x]=tot;
}
bool BFS()
{
	static int q[M];
	static bool v[M];
	int i,r=0,h=0;
	memset(v,0,sizeof v);
	v[S]=true;q[++r]=S;
	while(r!=h)
	{
		int x=q[++h];
		for(i=head[x];i;i=table[i].next)
			if(!v[table[i].to])
			{
				v[table[i].to]=true;
				q[++r]=table[i].to;
				if(table[i].to==T)
					return true;
			}
	}
	return false;
}
bool Judge(double r)
{
	int i,j;
	Initialize();
	for(i=1;i<=n;i++)
	{
		if( points[i].x-1<r || pos.y-points[i].y<r )
			Add(S,i);
		if( points[i].y-1<r || pos.x-points[i].x<r )
			Add(i,T);
	}
	for(i=1;i<=n;i++)
		for(j=1;j<i;j++)
			if( dis[i][j]<2*r )
				Add(i,j),Add(j,i);
	return !BFS();
}
double Bisection()
{
	double l=0,r=min(pos.x-1,pos.y-1);
	while(r-l>1e-4)
	{
		double mid=(l+r)/2.0;
		if( Judge(mid) )
			l=mid;
		else
			r=mid;
	}
	return (l+r)/2.0;
}
int main()
{
	int i,j;
	cin>>n>>pos;
	for(i=1;i<=n;i++)
		cin>>points[i];
	for(i=1;i<=n;i++)
		for(j=1;j<i;j++)
			dis[i][j]=Distance(points[i],points[j]);
	printf("%.2lf\n", Bisection() );
}

BZOJ 3007 拯救小云公主 二分答案+对偶图

标签：bzoj   bzoj3007   二分答案   对偶图   计算几何   

原文地址：http://blog.csdn.net/popoqqq/article/details/44224815