题目描述
在一个长方形框子里,最多有 𝑁 个相异的点,在其中任何一个点上放一个很小的油滴,那么这个油滴会一直扩展,直到接触到其他油滴或者框子的边界。必须等一个油滴扩展完毕才能放置下一个油滴。那么应该按照怎样的顺序在这 𝑁 个点上放置油滴,才能使放置完毕后所有油滴占据的总面积最大呢?(不同的油滴不会相互融合)
注:圆的面积公式 𝑆=𝜋𝑟^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;
}
麻烦看到的大佬不要直接划走,帮我看一下那个错误的代码错在哪里了,谢谢!