HDU 2022(凸包中的最大三角形面积)

HDOJ 2022 原题
题意：给定N个二维平面的x，y坐标，求出任意三个点组成三角形的最大面积。
题解: （分治法）用凸包算法求出其中若干组成的凸包,然后枚举栈中的元素组成三角形面积，得到最大值。

凸包（Convex Hull）是一个计算几何（图形学）中的概念。
在一个实数向量空间V中，对于给定集合X，所有包含X的凸集的交集S被称为X的凸包。X的凸包可以用X内所有点(X1，…Xn)的凸组合来构造.
在二维欧几里得空间中，凸包可想象为一条刚好包著所有点的橡皮圈。
用不严谨的话来讲，给定二维平面上的点集，凸包就是将最外层的点连接起来构成的凸多边形，它能包含点集中所有的点。

#include <bits/stdc++.h>

using namespace std;
int sta[50005];//使用手写栈存储排序后点的序号.
int T,cnt;
struct node
{
    double x,y;
    node friend operator -(node &p,node &q)//重载'-'号.
    {
        return {p.x-q.x,p.y-q.y};
    }
} p[50005];
double cross(node a,node b)//向量积(叉积)，用来判断两条向量线段的位置关系与求其距离.(叉积小于0,向量a在b的左边)
{
    return a.x*b.y-b.x*a.y;
}
double dis(node a,node b)//两点间的距离.
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cmp1(node &a,node &b)//以两向量叉积的正负为依据排序，使得结构数组中的元素与p[0]的连线从小到大逆时针旋转。
{
    int k=cross(a-p[0],b-p[0]);
    if(k>0)return 1;
    if(k==0)return dis(a,p[0])<dis(b,p[0]);
    return 0;
}
/*bool cmp2(node &a,node &b)//使用排序找y坐标最小点会超时QAQ....
{
    if(a.y==b.y)
        return a.x<b.x;
    return a.y<b.y;
}*/
double multi(node &a,node &b,node &c)//三个点的叉积(以a为向量起始点).
{
    return cross(b-a,c-a);
}
void Graham()//凸包
{
    cnt=1;
    sta[0]=0;
    sta[1]=1;
    for(int i=2; i<T; i++)
    {
        while(cnt>=1&&multi(p[sta[cnt-1]],p[sta[cnt]],p[i])<0)
            cnt--;
        sta[++cnt]=i;
    }
    return ;
}
double rotating_calipers()
{
    int i,j,k;
    double max1=-1;
    for(i=0; i<=cnt; i++)
        for(j=i+1; j<=cnt; j++)
            for(k=j+1; k<=cnt; k++)
                if(max1<abs(multi(p[sta[i]],p[sta[j]],p[sta[k]])))
                {
                    max1=abs(multi(p[sta[i]],p[sta[j]],p[sta[k]]));
                }
    return max1*0.5;
}
int main()
{
    int i;
    while(~scanf("%d",&T))
    {
        for(i=0; i<T; i++)
        {
            scanf("%lf %lf",&p[i].x,&p[i].y);
        }
        int k=0;
        for(int i=0; i<T; i++)//找y坐标最小的点.
        {
            if(p[i].x<p[k].x||(p[i].x==p[k].x&&p[i].y<p[k].y)) k=i;
        }
        swap(p[0],p[k]);
        sort(p+1,p+T,cmp1);
        Graham();
        printf("%.2lf\n",rotating_calipers());
    }
    return 0;
}

一篇比较容易理解凸包算法的文章.

OJBK