洛谷P1378 油滴扩展

题目描述

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

输入格式：

第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

思路

预处理出点与边和其余的的距离，dfs，注意只有已被放的油滴会有影响。

#include<cstdio>
#include<cmath>
using namespace std;
int n,x,y,xx,yy,a[100],b[100];
double f[100][100],max;
bool t[100];
double min(double v,double e)
{
    if (v<e) return v;
    else return e;
}
int abs(int v)
{
    if (v<0) return -v;
    else return v;
}
int sqr(int c)
{
    return c*c;
}
int ss(int c,double p)
{
    if (c>n) 
    {
        if (p>max) max=p;
        return 0;
    }
    for (int i=1;i<=n;i++)
    {
        if (t[i])
        {
            t[i]=false;
            double o=0;
            o=min(abs(a[i]-x),abs(xx-a[i]));
            o=min(o,abs(b[i]-y));
            o=min(o,abs(yy-b[i]));
            for (int j=1;j<=n;j++)
            {
                if (f[i][j]<=0)
                    o=min(f[i][j],o);
                if (not t[j]) 
                    o=min(f[i][j],o);
            }
            if (o<=0) 
            {
                t[i]=true;
                ss(c+1,p);
            }
            else
            {
                for (int j=1;j<=n;j++)
                {
                    f[i][j]=f[i][j]-o;
                    f[j][i]=f[i][j];
                }
                ss(c+1,p+3.1415926*o*o);
                t[i]=true;
                for (int j=1;j<=n;j++)
                {
                    f[i][j]=f[i][j]+o;
                    f[j][i]=f[i][j];
                }
            }
        }
    }
}
int main()
{
    scanf("%d",&n);
    scanf("%d%d%d%d",&x,&y,&xx,&yy);
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            f[i][j]=0xfffffff;
    for (int i=1;i<=n;i++)
    {
        scanf("%d%d",&a[i],&b[i]);
        for (int k=1;k<i;k++)
        {
            f[i][k]=sqrt(sqr(abs(a[i]-a[k]))+sqr(abs(b[i]-b[k])));
            f[k][i]=f[i][k];
        }
    }
    for (int i=1;i<=n;i++)
        t[i]=true;              
    ss(1,0);
    max=abs(xx-x)*abs(yy-y)-max;
    printf("%d",int(max+0.5));  
}