poj2079 Triangle （旋转卡壳之最大三角形）

题目链接：http://poj.org/problem?id=2079
题意：给出点集，求这个点集所组成的三角形的面积最大。
题解：求出凸包（不共线）后，遍历凸包的每个点做为三角形的第一个顶点，然后按凸包的顺序（我的代码是顺时针），依次为第二个点和第三个点。利用叉积的性质可以快速求得面积，按照面积的大小先更新第三个 点，直到下一个第三点与前两个点所形成的面积不大之前的面积，同理更新第二个点；

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<math.h>
#define eps 1e-8
#define inf 1e10
#define maxn 50020
struct point {double x,y;}p[maxn],save[maxn];
int n;
double xmult(point p1,point p2,point p0)
{
	return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
double max(double a,double b){return a>b? a:b;}
bool cmp(const point &a,const point &b )
{
	if(a.y==b.y)return a.x<b.x;
	return a.y<b.y;
}
int Graham()
{
	int i;
	std::sort(p,p+n,cmp);
	save[0]=p[0];
	save[1]=p[1];
	int top=1;
	for(i=0;i<n;i++)
	{
		while(top>0&& xmult(save[top],p[i],save[top-1])>=0)top--;
		save[++top]=p[i];
	}
	int mid=top;
	for(i=n-2;i>=0;i--)
	{
		while(top>mid&&xmult(save[top],p[i],save[top-1])>=0)top--;
		save[++top]=p[i];
	}
	return top;
}
double rc_tri()
{
	if(n<3)return 0;
	double maxx=0 ,temp;
	int i;
	int b=1,c=2;
	save[n]=save[0];
	save[n+1]=save[1];
	save[n+2]=save[2];
	for(i=0;i<n;i++)
	{
		while(xmult(save[c+1],save[b],save[i])
			-(temp=xmult(save[c],save[b],save[i]))>eps)
			c=(c+1)%n;
		maxx=max(maxx,temp);
		while(xmult(save[c],save[b+1],save[i])
			-(temp=xmult(save[c],save[b],save[i]))>eps)
			b=(b+1)%n;
		maxx=max(maxx,temp);
	}
	return maxx;
}
int main()
{
	while(scanf("%d",&n),n!=-1)
	{
		int i;
		for(i=0;i<n;i++)scanf("%lf%lf",&p[i].x,&p[i].y);
		n=Graham();
		printf("%.2lf\n",rc_tri()/2);
	}
	return 0;
}