POJ 2079/HDU 3934 旋转卡壳

题意：

给一些点，任选任选三个构成三角形，求这个三角形的最大面积。

题解：
旋转卡壳，水题了。。

n^2明显的过不了的，但是数据水，能过。

有没有nlogn的算法啊？？

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>

#define N 100010
#define EPS 1e-7

using namespace std;

struct PO
{
    double x,y;
}p[N],stk[N],o;

int top,n;
double ans;

inline PO operator -(PO a,PO b)
{
    PO c;
    c.x=a.x-b.x;
    c.y=a.y-b.y;
    return c;
}

inline double cross(PO a,PO b,PO c)
{
    return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}

inline int dc(double x)
{
    if(x>EPS) return 1;
    else if(x<-EPS) return -1;
    return 0;
}

inline bool cmp(const PO &a,const PO &b)
{
    if(dc(a.x-b.x)==0) return a.y<b.y;
    return a.x<b.x;
}

inline double getangle(PO &a,PO &b,PO &c,PO &d)
{
    return cross(o,b-a,d-c);
}

inline void read()
{
    for(int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y);
}

inline void graham()
{
    sort(p+1,p+1+n,cmp);
    top=-1;
    stk[++top]=p[1]; stk[++top]=p[2];
    for(int i=3;i<=n;i++)
    {
        while(top>=1&&dc(cross(stk[top-1],stk[top],p[i]))<=0) top--;
        stk[++top]=p[i];
    }
    int tmp=top;
    for(int i=n-1;i>=1;i--)
    {
        while(top>=tmp+1&&dc(cross(stk[top-1],stk[top],p[i]))<=0) top--;
        stk[++top]=p[i];
    }
}

inline void rotating_calipers()
{
    ans=0.0;
    for(int i=0;i<top;i++)
    {
        int s=i+1;
        for(int j=i+1,k=1;k<top;k++,j=(j+1)%top)
        {
            if(s==j) s=(s+1)%top;
            while(s!=i&&dc(getangle(stk[i],stk[j],stk[s],stk[(s+1)%top]))>0) s=(s+1)%top;
            if(s!=i) ans=max(ans,0.5*fabs(cross(stk[i],stk[j],stk[s])));
        }
    }
}

inline void go()
{
    graham();
    rotating_calipers();
    printf("%.2lf\n",ans);
}

int main()
{
    while(scanf("%d",&n)!=EOF) read(),go();//HDU
    //while(scanf("%d",&n)&&n!=-1) read(),go();//POJ
    return 0;
}

一个代码交两个OJ真是爽~