The Fortified Forest POJ - 1873 （二进制枚举+凸包）

题目地址：http://poj.org/problem?id=1873
二进制枚举所有被删的点，求出能构成围墙的树木的最大价值组合。
另外，G++不要用 %.2lf 今天又在这W了 换成C++ A了

#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<iostream>
#include<string.h>
using namespace std;
const int MAXN=1005;
const double PI=acos(-1.0);

struct point
{
    int x,y,l,v;
};
point list[MAXN];
point list1[MAXN];
int stack[MAXN],top;

int cross(point p0,point p1,point p2) //计算叉积  p0p1 X p0p2
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
double dis(point p1,point p2)  //计算 p1p2的 距离
{
    return sqrt((double)((p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y)));
}
bool cmp(point p1,point p2) //极角排序函数 ， 角度相同则距离小的在前面
{
    int tmp=cross(list[0],p1,p2);
    if(tmp>0) return true;
    else if(tmp==0&&dis(list[0],p1)<dis(list[0],p2)) return true;
    else return false;
}
void init(int n) //输入，并把  最左下方的点放在 list[0]  。并且进行极角排序
{
    int i,k;
    point p0;
   // scanf("%d%d",&list[0].x,&list[0].y);
    p0.x=list[0].x;
    p0.y=list[0].y;
    k=0;
    for(i=1;i<n;i++)
    {
       // scanf("%d%d",&list[i].x,&list[i].y);
        if( (p0.y>list[i].y) || ((p0.y==list[i].y)&&(p0.x>list[i].x)) )
        {
            p0.x=list[i].x;
            p0.y=list[i].y;
            k=i;
        }
    }
    list[k]=list[0];
    list[0]=p0;
    sort(list+1,list+n,cmp);
}
void graham(int n)
{
    int i;
    if(n==1) {top=0;stack[0]=0;}
    if(n==2)
    {
        top=1;
        stack[0]=0;
        stack[1]=1;
    }
    if(n>2)
    {
        for(i=0;i<=1;i++) stack[i]=i;
        top=1;

        for(i=2;i<n;i++)
        {
            while(top>0&&cross(list[stack[top-1]],list[stack[top]],list[i])<0) top--;
            top++;
            stack[top]=i;
        }
    }
}
int main()
{

  int n,ans,sum,sumv,xx,t;
  double suml,ans1,res;
  t=1;
  while(scanf("%d",&n)!=EOF)
  {
      if(n==0) break;
      for(int i=0;i<n;i++)
        scanf("%d%d%d%d",&list1[i].x,&list1[i].y,&list1[i].v,&list1[i].l);
    sum=990000;
      for(int i=1;i<(1<<n)-1;i++)
      {
          ans=0;
          suml=0;sumv=0;
          for(int j=0;j<n;j++)
          {
              if(i&(1<<j))
              {
                  sumv+=list1[j].v;
                  suml+=(double)list1[j].l;
              }
              else
              {
                  list[ans].x=list1[j].x;
                  list[ans++].y=list1[j].y;
              }
          }
          init(ans);
          graham(ans);
          top++;
          res=0;
            for(int k=0;k<top;k++)
            {
                 int j=(k+1)%top;
                  res+=dis(list[stack[k]],list[stack[j]]);
            }

            if(suml>res)
            {
                if(sum>sumv)
                {
                    ans1=suml-res;
                    sum=sumv;
                    xx=i;
                }
            }
      }
      printf("Forest %d\n",t);
      printf("Cut these trees: ");
      for(int j=0;j<n;j++)
      {
          if(xx&(1<<j))
            printf("%d ",j+1);
      }
      printf("\n");
      printf("Extra wood: %.2lf\n",ans1);
      printf("\n");
      t++;
  }
    return 0;
}
/*
Forest 1
Cut these trees: 2 4 5
Extra wood: 3.16
*/