凸包,模板---hdu1348

kuangbin大神的模板

const int MAXN=1000;

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

int cross(point p0,point p1,point p2){ //计算向量叉积 p0p1 X p0p2 (sin)
    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);   //list[0] 不变
}

void graham(int n){            //算法核心部分
    int i;
    if(n==1) {top=0;stack[0]=0;}
    if(n==2){       //只有0,1两个点
        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循环 是核心
            while(top>0 && cross(list[stack[top-1]],list[stack[top]],list[i])<=0)
                top--;        //若是新加入的这个点导致与前一个点不满足条件,说明这个店可以代替前一点
            top++;
            stack[top]=i;
        }

    }
}

你可能感兴趣的:(凸包,模板---hdu1348)