洛谷P1378 油滴扩展

题目描述

在一个长方形框子里，最多有 𝑁 个相异的点，在其中任何一个点上放一个很小的油滴，那么这个油滴会一直扩展，直到接触到其他油滴或者框子的边界。必须等一个油滴扩展完毕才能放置下一个油滴。那么应该按照怎样的顺序在这 𝑁 个点上放置油滴，才能使放置完毕后所有油滴占据的总面积最大呢？（不同的油滴不会相互融合）

注：圆的面积公式 𝑆=𝜋𝑟^2，其中 𝑟 为圆的半径。

输入格式

第一行，一个整数 𝑁。

第二行，四个整数 𝑥,𝑦,𝑥′,𝑦′，表示长方形边框一个顶点及其对角顶点的坐标。

接下来 𝑁 行，第 𝑖 行两个整数 𝑥𝑖,𝑦𝑖，表示盒子内第 𝑖 个点的坐标。

输出格式

一行，一个整数，长方形盒子剩余的最小空间（结果四舍五入输出）。

输入输出样例

输入 #1

2
20 0 10 10
13 3
17 7

输出 #1

50

说明/提示

对于 100% 的数据，1≤𝑁≤6，坐标范围在 [−1000,1000] 内。

大致思路：dfs所有未使用过的油滴，然后检查这个油滴是否在任意一个油滴形成的圆内，如果通过验证，则遍历下一个油滴，直到所有油滴用完，如果这个方案大于原先占据的面积，更新占据的面积。

错误代码（60pts）：

#include<bits/stdc++.h>
using namespace std;
const int N=10;
const double PI=3.1415926;
int n;
double a,b,a1,b1;
double x[N],y[N];
double maxr[N];
bool used[N];
double occupied;
double r[N];
template<typename T>
inline T min3(T a,T b,T c)
{
	return min(min(a,b),c);
}
template<typename T>
inline T min4(T u,T v,T u1,T v1)
{
	return min(min(u,v),min(u1,v1));
}
inline double dist(double x,double y,double xx,double yy)
{
	return sqrt((xx-x)*(xx-x)+(yy-y)*(yy-y));
}
inline double max_radius(double x,double y)//保证a<=x<=a1,b<=y<=b1 
{
	return min4(y-b,b1-y,x-a,a1-x);
}
void dfs(int u,double sum)
{
	if(u>=n)
	{
		occupied=max(occupied,sum);
		return;
	}
	for(int i=0;i<n;i++)
		if(!used[i])
		{
			r[i]=min(r[i],maxr[i]);
			bool is_legal=true;
			for(int j=0;j<n;j++)
				if(used[j])
				{
					if(dist(x[i],y[i],x[j],y[j])<=r[j])
					{
						is_legal=false;
						break;
					}
					r[i]=min(r[i],dist(x[i],y[i],x[j],y[j])-r[j]);
				}
			if(!is_legal) continue;
			used[i]=true;
			dfs(u+1,sum+PI*r[i]*r[i]);
			used[i]=false;
			r[i]=maxr[i];
		}
}
int main()
{
	cin>>n;
	cin>>a>>b>>a1>>b1;
	if(a>a1) swap(a,a1);
	if(b>b1) swap(b,b1);
	for(int i=0;i<n;i++) cin>>x[i]>>y[i];
	double size=fabs(a1-a)*fabs(b1-b);
	//cout<<"size = "<<size<<endl;
	for(int i=0;i<n;i++) r[i]=maxr[i]=max_radius(x[i],y[i]);
	dfs(0,0.0);
	cout<<int(size-occupied+0.5)<<endl;
	return 0;
}

正确代码：

#include<bits/stdc++.h>
using namespace std;
const int N=10;
const double PI=acos(-1);
int n;
double a,b,a1,b1;
double x[N],y[N];
bool used[N];
double occupied;
double r[N];
template<typename T>
inline T min4(T u,T v,T u1,T v1)
{
	return min(min(u,v),min(u1,v1));
}
inline double dist(double x,double y,double xx,double yy)
{
	return sqrt((xx-x)*(xx-x)+(yy-y)*(yy-y));
}
double eval(int j)
{
	double ans=min4(fabs(x[j]-a),fabs(x[j]-a1),fabs(y[j]-b),fabs(y[j]-b1));
	for(int i=0;i<n;i++)
		if(i!=j&&used[i])
		{
			double d=dist(x[i],y[i],x[j],y[j]);
			ans=min(ans,max(d-r[i],0.0));
		}
	return ans;
}
void dfs(int u,double sum)
{
	if(u>=n)
	{
		occupied=max(occupied,sum);
		return;
	}
	for(int i=0;i<n;i++)
		if(!used[i])
		{
			r[i]=eval(i);
			used[i]=true;
			dfs(u+1,sum+PI*r[i]*r[i]);
			used[i]=false;
		}
}
int main()
{
	cin>>n;
	cin>>a>>b>>a1>>b1;
	if(a>a1) swap(a,a1);
	if(b>b1) swap(b,b1);
	for(int i=0;i<n;i++) cin>>x[i]>>y[i];
	double size=fabs(a1-a)*fabs(b1-b);
	dfs(0,0.0);
	cout<<int(size-occupied+0.5)<<endl;
	return 0;
}

麻烦看到的大佬不要直接划走，帮我看一下那个错误的代码错在哪里了，谢谢！