凸包算法(Graham扫描法)详解

先说下基础知识,不然不好理解后面的东西

两向量的X乘p1(x1,y1),p2(x2,y2)

 

 p1Xp2如果小于零则说明  p1在p2的逆时针方向

如果大于零则说明 p1在p2的顺时针方向

struct node{
    double x,y;
    node friend operator -(node a,node b)//对减法符号进行重载 
    {
        return {a.x-b.x,a.y-b.y};
    }
}p[10000],s[10000];
double X(node a,node b){
    return a.x*b.y-a.y*b.x;
}

 

 这个方法很有用处。比如判断一个点是否在一条线段的左边还是右边,可以用X乘来判断,或者判断两条线段是否相交

 

 

接着说说凸包    Graham扫描法

1.在平面上一些散乱的点,首先  找找到这些点中处于最左下方的点

        for(int i=1;i<=N;i++)
        cin>>p[i].x>>p[i].y;
        int k=1;
        for(int i=2;i<=N;i++)
        {
            if(p[i].y<p[k].y||(p[k].y==p[i].y&&p[i].x<p[k].x))
            k=i;
        } 
        swap(p[1],p[k]);        

2.对这些点进行排序。把按照极角(polar angle)从小到大排序(以 p1为极点),极角相同的点按照到的距离从小到大排序。

int cmp(node a,node b)
{
    double x=X(a-p[1],b-p[1]);//以p[1]为极点,通过X乘来判断 
    
    if(x>0) return 1;//让a处于b的顺时针 
    if(x==0&&dis(a,p[1])<dis(b,p[1]))return 1;//角度相同看距离 
    return 0;
}
    

sort(p+2,p+N+1,cmp);

3.再开一个结构体数组s 来储存凸包最外围的点,也就是结果,这个有点容易让人搞迷。

遍历剩下的点,while循环把发现不是凸包顶点的点移除出去,因为当逆时针遍历凸包时,我们应该在每个顶点向左转。因此当while循环发现在一个顶点处没有向左转时,就把该顶点移除出去。

至于如何判断向左向右则是根据叉积来判断,前面我们已经解决过这个问题了

double multi(node a,node b,node c)
{
    return X(b-a,c-a);
}




      s[1]=p[1];
        s[2]=p[2];
        int t=2;
        for(int i=3;i<=N;i++)
        {
            // 发现在栈里边一个顶点处没有向左转时,就把该顶点移除出去
            while(t>=2&&multi(s[t-1],s[t],p[i])<=0) t--;
            s[++t]=p[i];
        }

这个是求凸包的周长的

hdu1392    http://acm.hdu.edu.cn/showproblem.php?pid=1392

算是模板题吧

#include<bits/stdc++.h>
using namespace std;

struct point{
    double x,y;
    point friend operator -(point a,point b)
    {return {a.x-b.x,a.y-b.y};}
}p[105],s[105];
double dis(point a,point b)
{
    point c=a-b;
    return sqrt(c.x*c.x+c.y*c.y);
}
double X(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}
int cmp(point a,point b)
{
    double x=X(a-p[1],b-p[1]);
    
    if(x>0) return 1;
    if(x==0&&dis(a,p[1])<dis(b,p[1])) return 1;
    return 0;
}
double multi(point p1,point p2,point p3)
{
    return X(p2-p1,p3-p1);
}
int main()
{
    int N;
    while(scanf("%d",&N),N)
    {
        for(int i=1;i<=N;i++) cin>>p[i].x>>p[i].y;
        
        if(N==1)
        {
            printf("0.00\n");
            continue;
        }
        else if(N==2)
        {
            printf("%.2lf\n",dis(p[1],p[2]));
            continue;
        }
        
        int k=1;
        for(int i=2;i<=N;i++)
        if(p[i].y<p[k].y||(p[i].y==p[k].y&&p[i].x<p[k].x))k=i;
        swap(p[1],p[k]);
        
        sort(p+2,p+1+N,cmp);
        
        s[1]=p[1];
        s[2]=p[2];
        int t=2;
        for(int i=3;i<=N;i++)
        {
            while(t>=2&&multi(s[t-1],s[t],p[i])<=0) t--;
            s[++t]=p[i];
        }
        double sum=0;
        for(int i=1;i<t;i++)
        {
            sum+=dis(s[i],s[i+1]);
        }
        printf("%.2lf\n",sum+dis(s[1],s[t]));
    }
    return 0;
}

 

emmm  再来个求任意多边形的面积

struct Point {
    double x, y;
};

//计算任意多边形的面积,顶点按照顺时针或者逆时针方向排列
double polygon_area(Point *p, int n)
{
    if(n < 3) return 0;

    double sum = 0;
    p[n + 1] = p[1];
    for(int i = 1; i <= n; i++)
        sum += p[i].x * p[i + 1].y - p[i].y * p[i + 1].x;//可以理解为不管这个多边形在哪,都以原点为分割点,就算原点在外面也可以算出,因为有正负可以抵消掉多余的
    sum = fabs(sum / 2.0);
    return sum;
}

 再来个求面积均匀的多边形重心

需要把多边形以p[0]为分界点  分成n-2个三角形,求出这些三角形的重心(i,j),乘以该三角形的面积,如上图公式

#include<bits/stdc++.h>
using namespace std;
struct node{
    double x,y;
    node friend operator -(node a,node b)
    {
        return {a.x-b.x,a.y-b.y};
    }
    double friend operator *(node a,node b)//对*进行重载  node*node 相当于X乘 
    {
        return a.x*b.y-a.y*b.x;
    }
}a[1000010];
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n;
        cin>>n;
        for(int i=1;i<=n;i++) cin>>a[i].x>>a[i].y;
        
        double S=0,X=0,Y=0;
        for(int i=2;i<n;i++)
        {
            double x=(a[i]-a[1])*(a[i+1]-a[1]);//这个乘和下面的不一样,这时X乘,求出三角形面积 
            X+=(a[1].x+a[i].x+a[i+1].x)*x;//重心(没除以3)乘以面积 
            Y+=(a[1].y+a[i].y+a[i+1].y)*x;
            S+=x;
        }
        printf("%.2lf %.2lf\n",X/S/(double)3,Y/S/(double)3);//除以3为重心 
     } 
    return 0;
}

 

posted @ 2018-08-10 17:07  踩在浪花上  阅读(18407)  评论(0编辑  收藏  举报