poj 1873 The Fortified Forest 搜索+凸包

原创 2016年08月31日 13:00:34

链接


题意:

给出一些树(2<=n<=15),每个树有一个二维坐标,现在要砍下一些树,去做成围栏把其它树围起来。每棵树都有价值,还有砍了后能做成多少米的围栏。
为了节约材料,尽量使围栏的周长小。
现在要求使砍树的价值和最小,在此情况下要求砍树的数目最小。
要求输出砍掉哪些树,并且做完围栏后多余的材料可以围多少米。


解法:

枚举所有状态(选中哪些树,不选哪些树),然后进行凸包。


代码:

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;

#define all(x) (x).begin(), (x).end()
#define for0(a, n) for (int (a) = 0; (a) < (n); (a)++)
#define for1(a, n) for (int (a) = 1; (a) <= (n); (a)++)
#define mes(a,x,s)  memset(a,x,(s)*sizeof a[0])
#define mem(a,x)  memset(a,x,sizeof a)
#define ysk(x)  (1<<(x))
#define sqr(x)  ((x)*(x))

typedef long long ll;
typedef pair<int, int> pii;
const int INF =0x3f3f3f3f;
const int maxn= 15    ;

const double PI=acos(-1.0);
const double eps=1e-10;
int dcmp(double x)
{
    if(fabs(x)<eps)  return 0;
    else return x<0?-1:1;
}
struct Point
{
    double x,y;
    int val,len,id;
    Point(double x=0,double y=0):x(x),y(y) {};
    bool operator ==(const Point B)const {return dcmp(x-B.x)==0&&dcmp(y-B.y)==0;}

    bool operator<(const Point& b)const
    {
        return dcmp(x-b.x)<0|| dcmp(x-b.x)==0 &&dcmp(y-b.y)<0;
    }
};
typedef Point Vector;
Vector operator -(Vector A,Vector B) {return Vector(A.x-B.x,A.y-B.y); }


double Cross(Vector A,Vector B){return A.x*B.y-A.y*B.x;}

double Dot(Vector A,Vector B){return A.x*B.x+A.y*B.y;}


Vector operator +(Vector A,Vector B) {return Vector(A.x+B.x,A.y+B.y); }
Vector operator *(Vector A,double p) {return Vector(A.x*p,A.y*p); }
Vector operator /(Vector A,double p) {return Vector(A.x/p,A.y/p); }
Vector operator -(Vector A)  {return  Vector(-A.x,-A.y);}

double dis(Point A,Point B)
{
    return sqrt(sqr(A.x-B.x)+sqr(A.y-B.y));
}

Point po[maxn+10],chp[maxn+10];
bool use[maxn+4];
int n;
double ConvexHell(Point *p,int n,Point *ch,int rea)
{

     int m=0;
    for(int i=0;i<n;i++) if(use[i])
    {
        while(m>1&&Cross(ch[m-1]-ch[m-2],p[i]-ch[m-1]   )<=0   )  m--;
        ch[m++]=p[i];
     }
     int k=m;
    for(int i=n-2;i>=0;i--) if(use[i])
    {
        while(m>k &&Cross(ch[m-1]-ch[m-2],p[i]-ch[m-1] )<=0)   m--;
        ch[m++]=p[i];
    }
    if(rea>1)  m--;
//    if(m<=1)  return 0;
//    if(m==2)  return dis(ch[0],ch[1]);
    double D=0;
    for0(i,m)
    {
        D+=dis(ch[i],ch[ (i+1)%m ]);
    }
//    cout<<"Convex "<<m<<endl;
    return D;

}


void solve()
{
    int ed=ysk(n)-1;
    int minVal=-1,num=-1,best=-1;
    double extr=-1;
    for0(s,ed+1)
    {
        int v=0,cnt=0,L=0;
        for0(i,n)
        {
            if(s&ysk(i)) use[i]=1;
            else  {use[i]=0;L+=po[i].len;v+=po[i].val,cnt++;}
        }
        double D=ConvexHell(po,n,chp,n-cnt);
//        cout<<s<<" "<<D<<endl;
        if(dcmp(L-D)<0)  continue;
        if(minVal<0||minVal>v||minVal==v&&num>cnt )  {minVal=v;num=cnt;best=s;extr=L-D;}

    }
    printf("Cut these trees:");
    vector<int >ve;
    for0(i,n) if( !(ysk(i)&best) ) ve.push_back(po[i].id+1);
    sort(all(ve));
    for0(i,ve.size())  printf(" %d",ve[i]);
    putchar('\n');
    printf("Extra wood: %.2f\n\n",extr);

}
int main()
{
   std::ios::sync_with_stdio(false);
   int kase=0;
   while(cin>>n&&n)
   {

       for0(i,n) {cin>>po[i].x>>po[i].y>>po[i].val>>po[i].len;po[i].id=i;}
       sort(po,po+n);
      printf("Forest %d\n",++kase);
       solve();
   }
   return 0;
}

相关文章推荐

POJ 1873 The Fortified Forest(凸包+DFS枚举)

题目链接: POJ 1873 The Fortified Forest 题意: 有n棵树,知道了每棵树的坐标和高度以及价值,需要砍掉一些树能用这些树的长度对没砍掉的树做成一个围栏, 要保证剩下...
  • Ramay7
  • Ramay7
  • 2016年04月19日 01:18
  • 153

POJ 1873 The Fortified Forest (计算几何,枚举+凸包)

看了很久了的题了,集训的时候就看过,每棵树有自己的价值和长度,现在要砍掉一些书做成栅栏,来围起剩余的树,并要使砍掉的树价值总和最小。 当时觉得有价值,树也有砍与不砍2种状态,就觉得这是DP,然后在死...

POJ 1873 The Fortified Forest 凸包+枚举

题意:一个人有好多树,这些树有坐标,有价值,以及砍下这颗树能造多长的篱笆。 现在他要保护自己的树不被别人砍,就需要用篱笆把树围起来,但是需要从自己的树中砍下来几棵树造成篱笆,所以问题就来了,如何砍最...

POJ 1873 The Fortified Forest (状态压缩+凸包 WF 推荐)

POJ 1873 The Fortified Forest (状态压缩+凸包 WF 推荐)

POJ 1873 The Fortified Forest 暴力凸包

题意:给出n棵树的坐标,树的高度和树的价值,从这些树中砍掉一些(整棵整棵的)做围栏把剩余的树围起来,使得消耗的树的价值最小。输出应砍掉哪里些树以及剩余的材料的长度。(如果砍掉的价值相同,则取砍掉数目少...
  • hcbbt
  • hcbbt
  • 2013年11月15日 11:18
  • 2160

POJ 1873 The Fortified Forest(凸包+枚举)

POJ 1873 The Fortified Forest(凸包+枚举) http://poj.org/problem?id=1873 题意:        有n颗树,每颗树有x,y坐标和它的价...

POJ 1873 The Fortified Forest (枚举+凸包)

看了很久了的题了,集训的时候就看过,每棵树有自己的价值和长度,现在要砍掉一些书做成栅栏,来围起剩余的树,并要使砍掉的树价值总和最小。 当时觉得有价值,树也有砍与不砍2种状态,就觉得这是DP,然后在死...

POJ 1873 The Fortified Forest 计算凸包

题目描述:http://poj.org/problem?id=1873 题目大意: 一片森林中有n棵树,现在需要砍掉其中的一些树做成围栏把剩下的树围起来。给定每棵树的(x,y)坐标、价值和砍掉之后...
  • f_cpp
  • f_cpp
  • 2014年11月09日 08:06
  • 158

poj 1873 The Fortified Forest(凸包+枚举)

【题目大意】:给出n棵树,每棵树有坐标,权值,和长度,现在要砍下其中的某些树,来做围栏来把剩下的树围起来,问损失的数最小价值是多少。 【解题思路】:明显求围栏就是一个凸包问题,由于树最多只有...
  • new_wu
  • new_wu
  • 2012年04月12日 01:26
  • 376

poj 1873 The Fortified Forest (位运算枚举 + 凸包周长)

题目链接:http://poj.org/problem?id=1873 大意:有一片N棵树的森林,要从中砍掉几棵树做成篱笆,把剩下的树围起来 输入:给N课树,每棵树的坐标是x,y,每棵树有一个vi和l...
  • TheHide
  • TheHide
  • 2016年07月23日 16:41
  • 157
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:poj 1873 The Fortified Forest 搜索+凸包
举报原因:
原因补充:

(最多只允许输入30个字)