• P1378 油滴扩展


    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
    暴力+枚举
     1 #include<iostream>
     2 #include<cstdio>
     3 #include<algorithm>
     4 #include<cmath>
     5 #include<cstring>
     6 
     7 using namespace std;
     8 const int N = 10;
     9 const double pi = acos(-1);
    10 
    11 double ans,r[N];
    12 int n,sx,sy,ex,ey;
    13 int x[N],y[N],s[N];
    14 bool v[N];
    15 
    16 double calc(int a,int b,int c,int d)
    17 {
    18     return sqrt((double)(a-c)*(double)(a-c)+(double)(b-d)*(double)(b-d));
    19 }
    20 
    21 void renew()
    22 {
    23     double sum = 0;
    24     for(int i=1;i<=n;++i)
    25         sum += pi*r[i]*r[i];
    26     ans = max(ans,sum);
    27 } 
    28 
    29 void dfs(int dep)
    30 {
    31     if(dep>n) renew();
    32     for(int i=1;i<=n;++i)
    33         if(!v[i])
    34         {
    35             r[i] = s[i];
    36             v[i] = true;
    37             for(int j=1;j<=n;++j)
    38                 if(v[j] && i!=j)
    39                 {
    40                     double dis = calc(x[j],y[j],x[i],y[i]);
    41                     dis = dis-r[j];
    42                     r[i] = min(r[i],dis);
    43                     if(r[i]<0) r[i]=0;
    44                 }
    45             dfs(dep+1);
    46             v[i] = false;
    47         }
    48 }
    49 int main()
    50 {
    51     scanf("%d%d%d%d%d",&n,&sx,&sy,&ex,&ey);
    52     memset(s,0x3f,sizeof s);
    53     for(int i=1;i<=n;++i)
    54     {
    55         scanf("%d%d",&x[i],&y[i]);
    56         s[i] = min(min(abs(sx-x[i]),abs(ex-x[i])),min(abs(sy-y[i]),abs(ey-y[i])));
    57     }
    58     dfs(1);
    59     double tmp = (double)(sx-ex)*(double)(sy-ey);
    60     if(tmp<0) tmp = -tmp;
    61     ans = tmp-ans;
    62     int anss = floor(ans);
    63     if(double(ans-anss)>=0.5) anss++;
    64     printf("%d
    ",anss);
    65     return 0;
    66 }
  • 相关阅读:
    smarty对网页性能的影响--开启opcache
    1stopt8.0 代码示例
    1stopt、matlab和python用morris、sobol方法实现参数敏感性分析
    MATLAB 实现sobol参数敏感性分析
    matlab中自带的sobol的函数提供的sobol序列
    matlab和fortran混合编程
    mathematic语法基础
    fortran常用语句--读写带注释文档、动态数组等语法
    fortran语言调用fortran写的dll
    C语言函数指针与 c#委托和事件对比
  • 原文地址:https://www.cnblogs.com/mjtcn/p/6868750.html
Copyright © 2020-2023  润新知