凸包模板

struct pt{
    double x,y;
}k[MAXN];
bool cross(const pt&a,const pt&b,const pt&z){//za叉乘zb向量
    return (a.x-z.x)*(b.y-z.y)-(a.y-z.y)*(b.x-z.x)>=0;
}
bool cmp(const pt&a,const pt&b){return cross(a,b,k[0]);}
int ans[MAXN],num,n;
void graham(){
    int cnt=0;//记录y值最小的index
    for(int i=1;i<n;++i)if((k[cnt].y>k[i].y)||
    (k[cnt].y==k[i].y&&k[cnt].x>k[i].x))cnt=i;
    swap(k[0],k[cnt]);//放到第一个
    sort(k+1,k+n,cmp);//排序
    ans[0]=0;
    ans[1]=1;
    num=1;//先确定了两个点。
    for(int i=2;i<n;++i){
        while(num>1&&!cross(k[ans[num]],k[i],k[ans[num-1]]))
            num--;
            //如果vector k[ans[num-1]]->k[i]不在k[ans[num-1]]->k[ans[num]]左侧
            //出栈 
        ans[++num]=i;
    }//num+1就是凸包点的个数
}
上一篇:JSP页面显示乱码


下一篇:计算机图形学复习