poj2079 triangle 凸包 旋转卡壳

传送门

题意就是从给定的点中找出面积最大的三角形。

很明显三个点都是凸包上的点，要求一下凸包，但是三角形的边长不一定是凸包的边！！！

我们一般做旋转卡壳的时候是固定两个点，来找第三个点，

但此题由于边长不一定是凸包的边长，所以底边还要判断一遍……

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAXN=1000005;
const double PI=acos(-1.0);
typedef long long ll;
struct point
{
    double x,y;
};
point list[MAXN];
point stack[MAXN];
int top;

/**叉积**/
double cross(point p0,point p1,point p2)
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
/**点积**/
double dot(point p0,point p1,point p2)
{
    return (p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y);
}
/**距离的平方**/
double dis(point p1,point p2)
{
    return (double)(p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y);
}
/**极角排序**/
bool cmp(point p1,point p2)
{
    double 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 graham(int n)
{
    int i;
    if(n==1) {top=0;stack[0]=list[0];}
    if(n==2)
    {
        top=1;
        stack[0]=list[0];
        stack[1]=list[1];
    }
    if(n>2)
    {
        for(i=0;i<=1;i++) stack[i]=list[i];
        top=1;

        for(i=2;i<n;i++)
        {
            while(top>0&&cross(stack[top-1],stack[top],list[i])<=0) top--;
            top++;
            stack[top]=list[i];
        }
    }
}
/**旋转卡壳**/
double rotating_calipers(point p[],int n)
{
    int k1=1,k2=2;
    p[n]=p[0];
    double ans=0;
    for(int i=0;i<n;i++)
    {
        while(fabs(cross(p[i],p[k1],p[k2]))<fabs(cross(p[i],p[k1+1],p[k2])))
            k1=(k1+1)%n;
        ans=max(ans,0.5*fabs(cross(p[i],p[k1],p[k2])));
        while(fabs(cross(p[i],p[k1],p[k2]))<fabs(cross(p[i],p[k1],p[k2+1])))
            k2=(k2+1)%n;
        ans=max(ans,0.5*fabs(cross(p[i],p[k1],p[k2])));
    }
    return ans;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=-1)
    {
        if(n==-1) break;
        scanf("%lf%lf",&list[0].x,&list[0].y);
        int k=0;
        for(int i=1;i<n;i++)
        {
            scanf("%lf%lf",&list[i].x,&list[i].y);
            if(list[k].y<list[i].y||(list[k].y==list[i].y&&list[k].x<list[i].x))
                k=i;
        }
        swap(list[k],list[0]);
        sort(list+1,list+n,cmp);
        graham(n);
        printf("%.2f\n",rotating_calipers(stack,top+1));
    }
    return 0;
}