P1378 油滴擴充套件

自為風月馬前卒發表於2017-05-18
套件

題目描述

在一個長方形框子裡，最多有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


首先DFS出一個完整的排列
然後在進行判斷，取能取的值的最小值
注意當一個伸展距離是負數的時候要改成0

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<cmath>
  5 #include<algorithm>
  6 using namespace std;
  7 double dis[7][7];
  8 int map[2001][2001];
  9 struct node
 10 {
 11     int bh;//   編號 
 12     int wx,wy;//位置 
 13     double kzjl;// 擴充套件距離 
 14 }a[2001];
 15 int x11,y11,x22,y22,n,ans=0x7ffffff,minx,miny,maxx,maxy;
 16 int how[7];
 17 int vis[7];
 18 int fact[10]={0,1,2,6,24,120,720};
 19 void deal_dis()
 20 {
 21     for(int i=1;i<=n;i++)
 22         for(int j=1;j<=n;j++)
 23             dis[i][j]=sqrt((a[i].wx-a[j].wx)*(a[i].wx-a[j].wx)+(a[i].wy-a[j].wy)*(a[i].wy-a[j].wy));
 24 }
 25 
 26 void pd()
 27 {
 28     for(int i=1;i<=n;i++)
 29     a[i].kzjl=0;
 30     for(int i=1;i<=n;i++)
 31     {
 32         a[how[i]].kzjl=min(min(abs(maxy-a[how[i]].wy),abs(a[how[i]].wy-miny)),min(abs(maxx-a[how[i]].wx),abs(a[how[i]].wx-minx)));
 33         for(int j=1;j<i;j++)
 34             a[how[i]].kzjl=min((dis[how[i]][how[j]]-a[how[j]].kzjl),a[how[i]].kzjl);
 35         if(a[how[i]].kzjl<0)
 36         a[how[i]].kzjl=0;
 37     }
 38     double now=0;
 39     for(int i=1;i<=n;i++)
 40     {
 41         now+=a[i].kzjl*a[i].kzjl*acos(-1);
 42     }
 43     now=(abs(maxx-minx))*(abs(maxy-miny))-now;
 44     now=(int)(now+0.5);
 45     ans=min(ans,(int)now);
 46 }
 47 
 48 void dfs(int p,int num)// num：第i滴，p：第i滴的編號 
 49 {
 50     how[num]=p;
 51     if(num==n)
 52     {
 53         pd();
 54         return ;
 55     }
 56     for(int i=1;i<=n;i++)
 57     {
 58         if(vis[i]==0)
 59         {
 60             vis[i]=1;
 61             dfs(i,num+1);
 62             vis[i]=0;
 63         }
 64     }
 65     
 66 }
 67 
 68 int main()
 69 {
 70     //freopen("oilbox.in","r",stdin);
 71     //freopen("oilbox.out","w",stdout);
 72     
 73     
 74     scanf("%d%d%d%d%d",&n,&x11,&y11,&x22,&y22);
 75     minx=min(x11,x22);miny=min(y11,y22);maxx=max(x11,x22);maxy=max(y11,y22);
 76     
 77     if(n==0)
 78     {
 79         printf("%d",(maxx-minx)*(maxy-miny));
 80         return 0;
 81     }
 82     
 83     for(int i=1;i<=n;i++)
 84     {
 85         scanf("%d%d",&a[i].wx,&a[i].wy);
 86         how[i]=i;
 87         a[i].bh=i;
 88         a[i].kzjl=0;
 89     }
 90     
 91     deal_dis();
 92     
 93     for(int i=1;i<=fact[n];i++)
 94     {
 95         vis[i]=1;
 96         dfs(i,1);
 97         vis[i]=0;
 98     }
 99     
100     printf("%d",ans);
101     return 0;
102 }

 











相關文章