POJ1873

题意：给你n个树，每个树都有自己的坐标和价值、长度。现在要砍掉一些树，来把剩下来的树给围起来。

优先砍掉价值便宜的树，如果有几种砍掉树方案切价值一样则选择砍掉数目最少的方案。然后价值和数量一样随便输出一种就行了。

思路：由于n比较小，就可以用dfs枚举出每一颗树的状态。在用剩下树 使用凸包求出需要围城的周长。

这道花了一个早上，烦死了，才学习了凸包。好多东西使用的时候出现了问题。比如求凸包中的while写成了if，莫名其妙就

原数组给替换了，而且题目给的数据还非常的水，根本检测不出来我的问题，盯着自己的代码，不断找错，不过最后终于A了。

#include<stdio.h>#include<math.h>#include<string.h>#include<algorithm>using namespace std;const int maxn=1<<16;const int maxm=1<<29;int n,e,k;double sum1=0,sum2=0,sum3=0;double maxn1=0,maxn2=0,maxn3=0;int ans[maxn],kj[maxn];struct node{    double x,y;    double val,l;} st[maxn],num[maxn],num2[maxn];int cnt=0;double dis(node p1,node p2){    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);}double multi(node p1,node p2,node p3){    return (p1.x-p3.x)*(p2.y-p3.y)-(p2.x-p3.x)*(p1.y-p3.y);}bool cmp(node p1,node p2){    if(multi(p1,p2,num2[0])>0) return true;    if(multi(p1,p2,num2[0])==0&&dis(p1,num2[0])<dis(p2,num2[0]))return true;    return false;}void dfs(int cur){    if(cur==n)    {        sum1=0,sum2=0,sum3=0;        cnt=0;        k=0;        for(int i=0; i<n; i++)        {            if(ans[i]==0)//0代表没有砍掉，保存在另一个数组中            {                num2[cnt++]=num[i];            }            else            {                sum1+=num[i].val;                sum2+=num[i].l;            }        }        sum3=n-cnt;        if(cnt==0||cnt==n) return;        if(cnt==1)//当只剩下来一个树的时候，不用花费长度来围它。下面可以和下面的合并，懒得合并了        {            if(maxn1>sum1)            {                maxn1=sum1;                maxn2=sum2;                maxn3=sum3;                for(int i=0; i<n; i++)                {                    kj[i]=ans[i];                }            }            else if(maxn1==sum1)            {                if(sum3<maxn3)                {                    maxn1=sum1;                    maxn2=sum2;                    maxn3=sum3;                    for(int i=0; i<n; i++)                        kj[i]=ans[i];                }            }            return ;        }        double sum=0;        if(cnt==2)//当剩下来两颗树的时候，要花费2倍的他们的距离长度。这要分开求，不然会有错误。        {            sum=dis(num2[0],num2[1]);            sum=2*sqrt(sum);        }        else        {            for(int i=0; i<cnt; i++)            {                if(num2[i].x<num2[k].x||(num2[i].x==num2[k].x&&num2[i].y<num2[k].y)) k=i;            }            swap(num2[0],num2[k]);            sort(num2+1,num2+cnt,cmp);            st[0]=num2[0];            st[1]=num2[1];            st[2]=num2[2];            e=2;            for(int i=3; i<cnt; i++)            {                while(e>1&&multi(num2[i],st[e],st[e-1])>=0) e--;                st[++e]=num2[i];            }            for(int i=0; i<=e; i++)                sum+=sqrt((double)dis(st[i],st[i==e?0:i+1]));        }        if(sum<=sum2)        {            if(maxn1>sum1)            {                /*                for(int i=0; i<n; i++)                    printf("%d ",ans[i]);                puts("");                printf("%lf %lf %lf\n",maxn1,sum,sum2);                */                maxn1=sum1;                maxn2=sum2-sum;                maxn3=sum3;                for(int i=0; i<n; i++)                    kj[i]=ans[i];            }            else if(maxn1==sum1)            {                if(sum3<maxn3)                {                    maxn1=sum1;                    maxn2=sum2-sum;                    maxn3=sum3;                    for(int i=0; i<n; i++)                        kj[i]=ans[i];                }            }        }        return ;    }    ans[cur]=0;    dfs(cur+1);    ans[cur]=1;    dfs(cur+1);}int main(){    int tt=1;    while(scanf("%d",&n)!=EOF)    {        if(n==0) break;        maxn1=maxn,maxn2=maxm;        for(int i=0; i<n; i++)            scanf("%lf %lf %lf %lf",&num[i].x,&num[i].y,&num[i].val,&num[i].l);        dfs(0);        printf("Forest %d\n",tt++);        printf("Cut these trees:");        for(int i=0; i<n; i++)            if(kj[i])                printf(" %d",i+1);        puts("");        printf("Extra wood: %.2lf\n\n",maxn2);    }    return 0;}