LG P1378 油滴扩展

题目描述

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

注：圆的面积公式V=pi*r*r，其中r为圆的半径。

输入输出格式

输入格式：

第1行一个整数N。

第2行为长方形边框一个顶点及其对角顶点的坐标，x,y,x’,y’。

接下去N行，每行两个整数xi,yi，表示盒子的N个点的坐标。

以上所有的数据都在[-1000，1000]内。

输出格式：

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

输入输出样例

输入样例#1： 复制

2
20 0 10 10
13 3
17 7

输出样例#1： 复制

50

一开始写的时候忘了考虑一个油滴位点在另一个油滴内，所以会得到r<0的情况- -

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define For(i,l,r) for(int i=l;i<=r;++i)
using namespace std;
int n,px[2],py[2],nx[7],ny[7],tx,ty,tx2,ty2,ans=2147483647,squ;//p[0]左下角的顶点,p[1]右上角的顶点 
bool flag[7];
double pi=acos(-1),tots=0,r[7],tmp;
void dfs(int num,int depth)
{
    double temp;
    r[num]=nx[num]-px[0];
    if(ny[num]-py[0]<r[num])r[num]=ny[num]-py[0];
    if(px[1]-nx[num]<r[num])r[num]=px[1]-nx[num];
    if(py[1]-ny[num]<r[num])r[num]=py[1]-ny[num];
    For(i,1,n)
    {
        if(!flag[i])
         continue;
        tx2=max(nx[i],nx[num]);
        tx=min(nx[i],nx[num]);
        ty2=max(ny[i],ny[num]);
        ty=min(ny[i],ny[num]);
        temp=sqrt((tx2-tx)*(tx2-tx)+(ty2-ty)*(ty2-ty));
        temp-=r[i];//如果两个油滴位点间距离小于已经扩展开的的r，小于0
        if(temp<r[num])
         r[num]=temp;
    }
    if(r[num]<0)//处理掉上一句注释的情况
     r[num]=0;
    temp=pi*r[num]*r[num];
    tots+=temp;
    if(depth==n)
    {
        tmp=squ-tots;
        if(tmp-(int)tmp/1>=0.5)
        tmp+=1;
        tmp=(int)tmp;
        if(tmp<ans)
         ans=tmp;
    }
    else
    {
        flag[num]=1;
        For(i,1,n)
         if(!flag[i])
          dfs(i,depth+1);
        flag[num]=0;
    }
    tots-=temp;
}
int main()
{
    scanf("%d",&n);
    scanf("%d %d %d %d",&px[0],&py[0],&px[1],&py[1]);
    if(px[0]>px[1])
     swap(px[0],px[1]);
    if(py[0]>py[1])
     swap(py[0],py[1]);
    squ=(py[1]-py[0])*(px[1]-px[0]);
    if(!n)
    {
        cout<<squ;
        return 0;
    }
    For(i,1,n)
     scanf("%d %d",&nx[i],&ny[i]);
    For(i,1,n)
     dfs(i,1);
    cout<<ans;
    return 0;
}