const int INF=0xfffffff ;
struct Point{
int x,y ;
} ;
Point p[50005],s[50005] ; // 1->top为凸包中的点。
int top ;
int direction(Point p1,Point p2,Point p3)
{
return (p3.x-p1.x)*(p2.y-p1.y)-(p2.x-p1.x)*(p3.y-p1.y) ;
}
int dis(Point p1,Point p2)
{
return (p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y) ;
}
int cmp(Point p1,Point p2)//极角排序
{
int temp=direction(p[0],p1,p2) ;
if(temp<0)return 1 ;
if(temp==0 && dis(p[0],p1)<dis(p[0],p2))return 1 ;
return 0 ;
}
void Graham(int n)
{
int pos,minx,miny ;
minx=miny=INF ;
for(int i=0 ;i<n ;i++)
if(p[i].x<minx || (p[i].x==minx && p[i].y<miny))
{
minx=p[i].x ;
miny=p[i].y ;
pos=i ;
}
swap(p[0],p[pos]) ;
sort(p+1,p+n,cmp) ;
p[n]=p[0] ;
s[0]=p[0] ;s[1]=p[1] ;s[2]=p[2] ;
top=2 ;
for(int i=3 ;i<=n ;i++)
{
while(direction(s[top-1],s[top],p[i])>=0 && top>=2)top-- ;
s[++top]=p[i] ;
}
}
扫描线凸包模板
猜你喜欢
转载自blog.csdn.net/qq_41713256/article/details/81347230
今日推荐
周排行