凸包面积--Graham扫描法

凸包面积–Graham扫描法

1.选取p0作为y坐标最小的点，如果y坐标相等，选取x坐标最小的点（即最左下的点）

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

2.对剩余的点相对与p0点的极角进行排序，除极角相等的点，即距离最远的点

bool cmp(point a,point b)
{
  float xx=(a.x-p[0].x)*(b.y-p[0].y)-(a.y-p[0].y)*(b.x-p[0].x);
  if(xx)
  {
    return (a.x-p[0].x)*(b.y-p[0].y)<(a.y-p[0].y)*(b.x-p[0].x)
  }
  else return a.x<b.x;
}

3.折线段拐向判断，若没有左转的时候出栈点，否则入栈

差乘

float ss(point a,point b,point c) 
{     
     float x1=b.x-a.x;
     float y1=b.y-a.y;
     float x2=c.x-a.x;
     float y2=c.y-a.y;
     return x1*y2-x2*y1; 
}

判断是否入栈

  pp[0]=p[0];
        pp[1]=p[1];
        pp[2]=p[2];
        int top=2;
        for(int i=3;i<n;i++)
        {
            while(ss(pp[top-1],pp[top],p[i])>=0)top--;
            pp[++top]=p[i];
        }

4.计算面积

       float ans=0;
        while(top>=2)
        {
            float m=ss(pp[0],pp[top-1],pp[top-2])/2.0;
            if(m<0)m=0-m;
            ans+=m;
            top--;
        }

就想用一个白色的凸多边形把墙上的墨点盖住。你能告诉麦兜最小需要面积多大的凸多边形才能把这些墨点盖住吗？
现在，给出了这些墨点的坐标，请帮助麦兜计算出覆盖这些墨点的最小凸多边形的面积。
输入：

多组测试数据。第一行是一个整数T，表明一共有T组测试数据。
每组测试数据的第一行是一个正整数N(0< N < = 105)，表明了墨点的数量。接下来的N行每行包含了两个整数Xi和Yi（0<=Xi,Yi<=2000），表示每个墨点的坐标。每行的坐标间可能包含多个空格。

输出：

每行输出一组测试数据的结果，只需输出最小凸多边形的面积。面积是个实数，小数点后面保留一位即可，不需要多余的空格。

样例输入：

2
4
0 0
1 0
0 1
1 1
2
0 0
0 1

样例输出：

1.0
0.0

 #include<bits/stdc++.h> 
 using namespace std; 
 struct point
  {     
      float x;
      float y; 
  }p[110],pp[110]; 
  
//根据角排序
bool cmp(point a,point b)
{
    float xx=(a.x-p[0].x)*(b.y-p[0].y)-(a.y-p[0].y)*(b.x-p[0].x);
    if(xx>0||xx==0&&a.x<b.x)
    {
        return true;
    }
    else return false;
}//差积
float ss(point a,point b,point c)
{
    float x1=b.x-a.x;
    float y1=b.y-a.y;
    float x2=c.x-a.x;
    float y2=c.y-a.y;
    return x2*y1-x1*y2;}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n,k=0;
        cin>>n;
        for(int i=0;i<n;i++)
        {
            cin>>p[i].x>>p[i].y;
            if(p[i].y<p[k].y||(p[i].y==p[k].y&&p[i].x<p[k].x))
            {
                k=i;
            }
        }
        swap(p[0],p[k]);
        sort(p+1,p+n,cmp);
        pp[0]=p[0];
        pp[1]=p[1];
        pp[2]=p[2];
        int top=2;
        for(int i=3;i<n;i++)
        {
            while(ss(pp[top-1],pp[top],p[i])>=0)top--;
            pp[++top]=p[i];
        }
        float ans=0;
        if(n<=2)
        {
            printf("0.0\n");
            continue;
        }
        while(top>=2)
        {
            float m=ss(pp[0],pp[top-1],pp[top])/2.0;
            if(m<0)m=0-m;
            ans+=m;
            top--;
        }
        printf("%.1lf\n",ans);
    }
}