凸包模板

来源:互联网 发布:建筑图纸设计软件 编辑:程序博客网 时间:2024/06/11 06:01

抄了一下kuangbin大神的模板ORZ

#include<stdio.h>#include<math.h>#include<algorithm>#include<iostream>using namespace std;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{    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--;            stack[++top]=i;        }    }}