SSL2796 2017年10月27日提高组 油滴拓展（贪心）

2017年10月27日提高组 油滴拓展

Description

在一个长方型框子里，最多有N（0≤N≤6）个相异的点。在其中任何一个点上放一个很小的油滴，那么这个油滴会一直扩展，直到接触到其它油滴或者框子的边界。必须等一个油滴扩展完毕才能放置下一个油滴。那么应该按照怎样的顺序在这N个点上放置油滴，才能使放置完毕后所有油滴占据的总体积最大呢？（不同的油滴不会相互融合）
注：圆的面积公式V=pi*r*r，其中r为圆的半径。

Input

第一行一个整数N。
第二行为长方形边框一个顶点及其对角顶点的坐标，x，y，x’，y’。
接下去N行，每行两个整数xi，yi，表示盒子内N个点的坐标。
以上所有的整数都在[-1000, 1000]内。

分析：对于一个油滴来说能扩大的最大半径就是圆心到边界的最小距离，如果存在其他油滴，那能拓展的半径为两个圆圆心距离减去其中一个圆的半径，那么只需要统计一下每个圆的半径，先滴大的然后更新其他半径继续滴即可。

代码

#include <cstdio>
#include <cmath>
#include <algorithm>
#define N 10
using namespace std;

int xi,yi,xj,yj,n;
double ans,x[N],y[N];
struct arr
{
    double r;
    int n;
}a[N];

double fsqr(double x)
{
    return x*x;
}

int fabs(int x)
{
    return x>0?x:-x;
}

double minf(double x,double y)
{
    return x<y?x:y;
}

double distance(int i,int j,double r)
{
    if (sqrt(fsqr(x[i]-x[j])+fsqr(y[i]-y[j]))-r<0) return 0;
    return sqrt(fsqr(x[i]-x[j])+fsqr(y[i]-y[j]))-r;
}

int so(arr p,arr q)
{
    return p.r>q.r;
}

void sum(int xx)
{
    sort(a+1,a+xx+1,so);
    ans+=M_PI*a[1].r*a[1].r;
    a[0]=a[1];
    for (int i=1;i<xx;i++)
    {
        a[i]=a[i+1];
        a[i].r=min(distance(a[i].n,a[0].n,a[0].r),a[i].r);
    }
}

int main()
{
    scanf("%d",&n);
    scanf("%d%d%d%d",&xi,&yi,&xj,&yj);
    for (int i=1;i<=n;i++)
        scanf("%lf%lf",&x[i],&y[i]);
    for (int i=1;i<=n;i++)
    {
        a[i].r=min(min(fabs(x[i]-xi),fabs(x[i]-xj)),min(fabs(y[i]-yi),fabs(y[i]-yj)));
        a[i].n=i;
    }
    for (int i=n;i>0;i--) 
        sum(i);
    ans=fabs(xi-xj)*fabs(yi-yj)-ans;
    printf("%.0lf",ans);
}