凸包问题

#include<cstdio>
#include<algorithm>
#include<fstream>
#include<cstring>
#include<iostream>
using namespace std;
const int MAX=550;
struct point
{
    int x;
    int y;
    int v;
    int i;
    int tmp;
} Ini[MAX],res[MAX],p[MAX];
bool flag[MAX];
int ans[MAX];
bool cmp(point A,point B)  //同速度的排序，仍然是找x小的和y小的；
{
    if(A.y==B.y)return A.x<B.x;
    return A.y<B.y;
}
bool cmp1(point A,point B)//排序从大到小，以及同速度的x小点，同速度同x的y小点；这里是对原数组的排序；
{
    if(A.v>B.v)return true;
    if(A.v==B.v)
    {
        if(A.x<B.x)return true;
        if(A.x==B.x&&A.y<B.y)return true;
    }
    return false;
}
int cross(point A,point B,point C)//点积；
{
    return (B.x-A.x)*(C.y-A.y)-(C.x-A.x)*(B.y-A.y);
}
int Graham(point *p,int n)
{
    //if (n<3)return n;
    sort(p,p+n,cmp);
    memset(flag,false,sizeof(flag));
    int i;
    int top=0;
    for(i=0; i<n; i++)
    {
        while(top>=2&&cross(res[top-2],res[top-1],p[i])<0)/*构造下边界，若不在凸包上则将点出栈；(这里是判断如果点积小于0继续循环，知道找到凸包的点，而=0即共线情况也是凸包的点，也要跳出，所以这里不能取等)知道满足凸包条件跳出循环，然后记录；*/
        {
            top--;
        }
        res[top]=p[i];   //res数组存下在凸包上的点；
        res[top].tmp=i; //记录入栈点的编号；

        top++;
    }
    for(i=0;i<top;i++)
        flag[res[i].tmp]=true; 


    int t=top+1;
    for(i=n-2; i>=0; i--) //i>=0
    {
        while(top>=t&&cross(res[top-2],res[top-1],p[i])<0)  //构造上边界；

            top--;
        if(!flag[i]) res[top++]=p[i];  //满足条件且不在已构造的凸包内；
        //if(i==0)top--;
    }
    /*for(i=0; i<top; i++)
        printf("$%d %d\n",res[i].x,res[i].y);*/
    return top;//返回凸包上点的个数；
}
int main()
{
    int n;
    int i,j;
    //ifstream cin("A.txt");
    int __case=0;
    while(cin>>n&&n)
    {
        for(i=0; i<n; i++)
        {
            cin>>Ini[i].x>>Ini[i].y>>Ini[i].v;
            Ini[i].i=i;//标记位置id
        }
        sort(Ini,Ini+n,cmp1);//排序找出速度最大点（前面的好几位都是速度最大点），并按坐标排序；
        //cout<<Ini[i-1].x<<Ini[i-1].y;
        int tmp=0;
        for(i=0; i<n; i++)
        {
            if(Ini[0].v>0&&Ini[0].v==Ini[i].v) tmp++;
            else break;
        }//筛出最大速度的点；

        memset(ans,0,sizeof(ans));
        point po;
        for(i=0,j=0;i<tmp;i++)
        {
            po=Ini[i];
            //flag=1;
            if(i<tmp-1&&Ini[i+1].x==po.x&&Ini[i+1].y==po.y)//找重点 ，因为已经按排序过，所以重点一点是出现在相邻点之间；
            {
                while(i<tmp-1&&Ini[i+1].x==po.x&&Ini[i+1].y==po.y)
                {
                    ans[Ini[i].i]=-1;
                    i++;
                }
                ans[Ini[i].i]=-1;   //重点则不可能无限，先标记掉；
            }
            p[j++]=po;//保留所有的符合的点，注意此时重点也要留下一个点。
        }
        int m=Graham(p,j);//传入数组时不仅要有地址，还要传入长度；此时传入的数组是已经去完重点的，在函数中就不用判断了；

        printf("Case #%d: ",++__case);
        for(i=0; i<m; i++)
        {
            if(!ans[res[i].i])ans[res[i].i]=1;//如果不是重点且在凸包上，则使其为1；
        }
        for(i=0; i<n; i++)
        {
            if(ans[i]==1)printf("1");
            else printf("0");
        }
        printf("\n");
    }
    return 0;
}