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
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;
}