【洛谷P1378】油滴扩展

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

【问题分析】
简单地模拟一下顺序，然后贪心最大化当前点半径，重新判断其他点的最大可能半径即可，注意的是一旦当前点已被覆盖，则半径不可为负数，害得我WA了一次。

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const int N=7;
const double Pi=3.14159265;
int n,xxx,yyy,xx,yy;  int x[N],y[N],step[N];   long long s;
bool used[N];   double ans,zero=0.00;  double dis[N][N],closest[N],r[N],d[N][N],c[N];
double mymin(double a,double b)
{
    if (a>b)
        return b;
    return a;
}
double mymax(double a,double b)
{
    if (a>b)
        return a;
    return b;
}
void read()
{
    int i,j; double tempx,tempy;
    scanf("%d",&n);
    scanf("%d%d%d%d",&xxx,&yyy,&xx,&yy);
    s=abs(xx-xxx)*abs(yy-yyy);
    for (i=1;i<=n;i++)
        scanf("%d%d",&x[i],&y[i]);
    for (i=1;i<=n;i++)
        for (j=1;j<=n;j++)
        {
            tempx=abs(x[i]-x[j]);   tempx*=tempx;
            tempy=abs(y[i]-y[j]);   tempy*=tempy;
            dis[i][j]=sqrt(tempy+tempx);
            //俩点距离
        }
    for (i=1;i<=n;i++)
        closest[i]=mymin(mymin(abs(x[i]-xxx),abs(x[i]-xx)),mymin(abs(y[i]-yyy),abs(y[i]-yy)));//矩形边界
    return;
}
void mycopy()
{
    int i,j;
    for (i=1;i<=n;i++)
        for (j=1;j<=n;j++)
        d[i][j]=dis[i][j];
    for (i=1;i<=n;i++)
        c[i]=closest[i];
    return;
}
void work()
{
    int i,j;    double temp=0.00;
    mycopy();
    r[step[1]]=c[step[1]];
    for (i=2;i<=n;i++)
    {
        for (j=1;j<i;j++)
            d[step[i]][step[j]]-=r[step[j]];//修改距离
        for (j=1;j<i;j++)
            c[step[i]]=mymin(c[step[i]],d[step[i]][step[j]]);
        r[step[i]]=c[step[i]];
        r[step[i]]=mymax(r[step[i]],zero);//不为负
    }
    for (i=1;i<=n;i++)
        temp=temp+r[i]*r[i];
    ans=max(ans,temp);
    return;
}
void dfs(int k)
{
    int i;
    if (k==n+1)
    {
        work();
        return;
    }
    for (i=1;i<=n;i++)
        if (!used[i])
        {
            step[k]=i;
            used[i]=1;
            dfs(k+1);
            used[i]=0;
        }
    return;
}
int main()
{
    freopen("oilbox.in","r",stdin);
    freopen("oilbox.out","w",stdout);
    read();
    dfs(1);
    printf("%.0lf\n",s-ans*Pi);
    return 0;
}