洛谷 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

思路: 题目所给数据非常小,所以就算枚举出每一种方案也只需要枚举6!次,所以用dfs枚举出每一种情况.油滴扩展的最大值就是从该点到墙体或是到已经扩展的油滴的边界的最小值.

特殊点的处理: 如果有油滴被其他包含在内的话,需要特判,直接跳过该油滴,进行下一步搜索

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

const int maxn = 6+2, inf = 0x3f3f3f3f;
const double pi = acos(-1);
double px[maxn], py[maxn]; //点的坐标
double dis[maxn]; //点与墙/点之间的距离
double R[maxn], ans; //扩展后点的半径
double l, r, u, d; //框的左右上下
int n;
bool used[maxn];

void dfs( int s, double cnt ) {
    //printf("this is %d step, cnt = %lf\n", s, cnt);
    if( s == n ) {
        ans = max( ans, cnt );
        return;
    }
    bool isinclude = 0;
    for( int i = 1; i <= n; i++ ) {
        if( !used[i] ) {
            used[i] = 1;
            double Min = inf;
            for( int j = 0; j <= n; j++ ) {
                //如果是被扩展过的点
                if( ( used[j] && i != j ) ) {
                    if( sqrt( (px[i]-px[j])*(px[i]-px[j]) + (py[i]-py[j])*(py[i]-py[j]) ) - R[j] >= 0 )
                        Min = min( Min, sqrt( (px[i]-px[j])*(px[i]-px[j]) + (py[i]-py[j])*(py[i]-py[j]) ) - R[j] );
                    //如果被其他油滴包含住的话不使用这个点,直接进行下一层搜索
                    else {
                        isinclude = 1;
                        break;
                    }
                }
                //如果是墙体
                else if( j == 0 ) {
                    Min = min( Min, dis[i] );
                }
            }
            if( isinclude ) {
                dfs( s+1, cnt );
                used[i] = 0; //回溯
                continue;
            }
            else {
                R[i] = Min;
                dfs( s+1, cnt + R[i]*R[i]*pi );
            }
            used[i] = 0; //回溯
        }
    }
}

int main() {
    /***************************init****************************/
    scanf("%d", &n); //n代表几个点
    scanf("%lf %lf %lf %lf", &l, &u, &r, &d); //框的边界
    for( int i = 1; i <= n; i++ ) scanf("%lf %lf", &px[i], &py[i]);
    if( l > r ) swap( l, r );
    if( d > u ) swap( d, u );
    for( int i = 1; i <= n; i++ ) dis[i] = min( abs(px[i]-l), min( abs(r-px[i]), min( abs(u-py[i]), abs(py[i]-d )) ) );
    /***********************************************************/
    dfs(0, 0);
    cout << (int)round( (r-l)*(u-d)-ans ) << endl;
    return 0;
}