POJ 1696 Space Ant（卷包裹法）

 http://poj.org/problem?id=1696

这道题我花了将近两个小时才把他整出来，还是不错的，顺便学习了下卷包裹法，卷包裹法的讲解参考下面的链接

http://www.cnblogs.com/Booble/archive/2011/02/28/1967179.html

主要是我写的时候犯了一个错误，参见程序

#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

const int maxn=55;

struct node
{
   int idx;
   int x,y;
   bool flag;
}a[maxn],b[maxn];

//计算叉积
int det(int i,int num,int k)
{
   return (a[num].x-b[k].x)*(a[i].y-b[k].y)-(a[i].x-b[k].x)*(a[num].y-b[k].y);
}

//当三点共线的时候，选择距离较远的一点
bool dist(int i,int num,int k)
{
   return (a[i].x-b[k].x)*(a[i].x-b[k].x)+(a[i].y-b[k].y)*(a[i].y-b[k].y)<(a[num].x-b[k].x)*(a[num].x-b[k].x)+(a[num].y-b[k].y)*(a[num].y-b[k].y);
}

int main()
{
   int m,n,k,num;
   scanf("%d",&m);
   while(m--)
   {
      int minx=1005,miny=1005;
      scanf("%d",&n);
      for(int i=0;i<n;i++)
      {
         scanf("%d%d%d",&a[i].idx,&a[i].x,&a[i].y);
         a[i].flag=1;
         if(miny>a[i].y||(miny==a[i].y&&minx>a[i].x))
         {
            minx=a[i].x,miny=a[i].y;
            num=i;
         }     
      }
      //cout<<num<<endl;
      a[num].flag=0;
      k=0;
      b[k++]=a[num];//像这样写的时候要注意下面的k
      while(k!=n)
      {
         for(int i=0;i<n;i++)
            if(a[i].flag)
            {
               num=i;break;
            }
         for(int i=0;i<n;i++)
         {
            //下面这句话里面的k-1我刚开始写的是k，一直得不到正确数据，找了好久才找到，纠结~
            if(i!=num&&a[i].flag&&(det(i,num,k-1)<0||det(i,num,k-1)==0&&dist(i,num,k-1)))
               num=i;
         }
         a[num].flag=0;
         b[k++]=a[num];
      }
      printf("%d",k);
      for(int i=0;i<k;i++)
         printf(" %d",b[i].idx);
      printf("\n");
   }
  
   return 0;
}