BZOJ 1069: [SCOI2007]最大土地面积

Description

在某块平面土地上有N个点，你可以选择其中的任意四个点，将这片土地围起来，当然，你希望这四个点围成的多边形面积最大。

Input

第1行一个正整数N，接下来N行，每行2个数x,y，表示该点的横坐标和纵坐标。

Output

最大的多边形面积，答案精确到小数点后3位。

Sample Input

5
0 0
1 0
1 1
0 1
0.5 0.5

Sample Output

1.000

HINT

数据范围 n<=2000, |x|,|y|<=100000

题解

枚举四边形对角线，两边旋转卡壳找出面积最大的三角形。

注意：在凸包上共线的点最好删到只剩两个，否则有很多问题，比如精度和一些奇奇怪怪的问题。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
int n,zz;
struct dian {double x,y;} d[2005],q[2005];
void init()
{
	scanf("%d",&n);
	int i;
	for(i=1;i<=n;i++) scanf("%lf%lf",&d[i].x,&d[i].y);
}
double dis(dian a,dian b)
{
	return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
double operator*(dian a,dian b)
{
	return a.x*b.y-a.y*b.x;
}
dian operator-(dian a,dian b)
{
	dian t;
	t.x=a.x-b.x; t.y=a.y-b.y;
	return t;
}
bool operator<(dian a,dian b)
{
	double t=(a-d[1])*(b-d[1]);
	if(t==0) return dis(a,d[1])<dis(b,d[1]);
	if(t>0) return true;
	else return false;
}
void gradam()
{
	int i,t=1;
	for(i=2;i<=n;i++)
	   {if(d[i].y<d[t].y||(d[i].y==d[t].y&&d[i].x<d[t].x)) t=i;}
	swap(d[1],d[t]); sort(d+2,d+n+1);
	//for(i=1;i<=n;i++) printf("%.3lf %.3lf\n",d[i].x,d[i].y);
	q[++zz]=d[1]; q[++zz]=d[2];
	for(i=3;i<=n;i++)
	   {while(zz>1&&(d[i]-q[zz-1])*(q[zz]-q[zz-1])>=0)
	      zz--;
	    q[++zz]=d[i];
	   }
	//for(i=1;i<=zz;i++) printf("%.3lf %.3lf\n",q[i].x,q[i].y);
}
void rotating_c()
{
	q[zz+1]=d[1];
	double ans=0;
	int xx,yy,i,j;
	for(i=1;i<=zz;i++)
	   {xx=i%zz+1; yy=(i+2)%zz+1;
	    for(j=i+2;j<=zz;j++)
	       {while(xx%zz+1!=j&&(q[xx]-q[i])*(q[j]-q[i])<(q[xx+1]-q[i])*(q[j]-q[i]))
	           xx=xx%zz+1;
	        while(yy%zz+1!=i&&(q[j]-q[i])*(q[yy]-q[i])<(q[j]-q[i])*(q[yy+1]-q[i]))
	           yy=yy%zz+1;
	        ans=max(ans,(q[xx]-q[i])*(q[j]-q[i])+(q[j]-q[i])*(q[yy]-q[i]));
		   }
	   }
	printf("%.3lf\n",ans/2);
}
int main()
{
	freopen("land.in","r",stdin);
	freopen("land.out","w",stdout);
	init(); gradam(); rotating_c();
	return 0;
}