POJ 1873 The Fortified Forest 凸包+枚举组合

下面是我做这道题时的一些分析  大部分是错的  大家参考参考   我也是写下自己错的好好反思一下以后不犯这种错误
题意：这道题找的是花费最小的价值把树木围起来  花多少树木都没事  但是最后至少能剩下一颗  剩下两颗也不一定能用前面的所有树围起来  所以必须枚举到只剩一颗的情况
然后问题就来了：
1：怎么枚举所有情况?
2：怎么确定花费的最小价值?
3：枚举到只剩一棵树的情况   则找出价值最大的那一颗数就行了
4:枚举到只剩两颗棵树的时候也不用求凸包    计算剩余点的距离（两倍）和价值  找到一个价值最大的输出
5:其他情况 计算凸包的周长 和凸包包围树木的价值（标记数组看看树木是不是在凸包中  枚举的时候赋值  每次枚举都要初始化）   找到能满足篱笆条件的最大的一个组合且价值最大的情况  保存当前最佳状态 (思考一下能不能不用每次都计算凸包的周长和价值  找一找有什么剪枝优化没 )  
6:把价值从小到大排   价值相等高度从小到大排
这样排序后从头到尾枚举  每种情况只要找到满足条件就可以保证这是当前枚举情况的最优解（比如枚举的的是砍掉3颗树的情况 
只要有一个满足就可以输出了） 
枚举N个当中选一个的反向就是N个中选8个的选法   所以枚举操作可以大大简化


最后做的时候意识到既然我都枚举了所有情况那就都算一遍得了 反正数据量小也不会超时   而且不用判断特殊情况  因为所有情况在计算凸包的时候都计算了
刚开始做的时候自己写递归求组合弄了一上午没弄出来 主要是应该返回到什么位置没写出来要不都用自己写的了  后来就去看了别人的   代码用的是别人写的递归求的组合 
要是只剩两颗树就是两棵树的距离的2倍      一棵树的情况也一样    没有漏判的情况

下面看看代码   有注释

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define eps 1e-8
#define maxn 1000000
using namespace std;

struct Tree
{
    int x,y,v,l,n,ok;
}t[20],num[20],s[20];

struct H  //保存最优解
{
    int v,number[16],len;
    double l;
}h;

int T,e=1;

typedef struct Tree tree;

int cmp(tree a,tree b)//排序参数 求凸包用
{
    if(a.x<b.x)return 1;
    else if(a.x==b.x&&a.y<b.y)return 1;
    return 0;
}

double Dis(tree a, tree b)//计算距离
{
    return sqrt((double)((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)));
}

//极角比较， 返回-1: p0p1在p0p2的右侧，返回0:p0,p1,p2共线
int Cross(tree p1, tree p2, tree pb)
{
    int delta=(p1.x-pb.x)*(p2.y-pb.y)-(p2.x-pb.x)*(p1.y-pb.y);
    if (delta<0) return 1;
    else if (delta==0) return 0;
    else return -1;
}

int ConvexHull(tree *p,int n,tree *ch)//Andrew法寻找凸包  必须先排序
{
    int m = 0;
    for(int i = 0; i < n; i++)//寻找下凸包
    {
        if(p[i].ok)continue;
        while(m > 1 && Cross(ch[m-1],p[i],ch[m-2]) <= 0 )m--;
        ch[m++] = p[i];
    }
    int k = m;
    for(int i = n-2; i >= 0; i--)//寻找上凸包
    {
        if(p[i].ok)continue;
        while(m > k && Cross(ch[m-1],p[i],ch[m-2]) <= 0 )m--;
        ch[m++] = p[i];
    }
    if(n > 1)m--;
    return m;
}

void Solve(int m,int v,double l)
{
    double sum = 0;
    for(int i = 0; i < m; i++)//计算凸包的周长
        sum +=Dis(s[i],s[i+1]);
    if(v<h.v&&l-sum>0||fabs(l-sum)<eps)//价值不比原来的小就退出当前组合继续计算下一种组合
    {
        h.len = 0;
        h.v = v;//保存花费的最小价值
        h.l = l-sum;//保存剩余的树木长度
        for(int i = 0; i < T; i++)//更新砍掉的树木，最后排序从小到大输出
        if(t[i].ok)h.number[h.len++]=t[i].n;
    }
    for(int i = 0; i < T; i++)//这一定要初始化
        t[i].ok = 0;
}

void combine(int n, int m,const int M )
{
    for(int i=n; i>=m; i--)  // 注意这里的循环范围
    {
        num[m-1].n = i-1;
        if (m > 1)
            combine(i-1,m-1,M);
        else      // m == 1, 输出一个组合
        {
            int v = 0;
            double l = 0;
            for(int j=M-1; j>=0; j--)
            {
                t[num[j].n].ok=1;//标记此树已被砍掉
                v +=t[num[j].n].v;//计算当前组合砍掉的树的价值;
                l +=(double)t[num[j].n].l;//计算当前组合砍掉的树的长度
            }
	    //其实这里num这个定义成数组就行了   因为就只是递归记录下标而已
            int m = ConvexHull(t,T,s);//找出凸包
            Solve(m,v,l);
        }
    }
}

int main()
{
    #ifdef LOCAL
    freopen("in.txt","r",stdin);
    freopen("o.txt","w",stdin);
    #endif // LOCAL
    while((scanf("%d",&T),T)!=0)
    {
        for(int i = 0; i < T; i++)
        {
            scanf("%d%d%d%d",&t[i].x,&t[i].y,&t[i].v,&t[i].l);
            t[i].n = i+1;//标号
        }
        sort(t,t+T,cmp);//排序求凸包的必要步骤
        h.v = maxn;
        for(int i = 1; i < T; i++)
        combine(T,i,i);
        /*枚举砍掉的树的所有组合找出凸包*/
        sort(h.number,h.number+h.len);
        printf("Forest %d\nCut these trees:",e++);
        for(int i = 0; i < h.len; i++)
            printf(" %d",h.number[i]);
        printf("\nExtra wood: %.2f\n\n",h.l);
    }
    return 0;
}