1069: [SCOI2007]最大土地面积
Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 2506 Solved: 955
[ Submit][ Status][ Discuss]
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
Source
对角线可以枚举,然后就是在两边各找一个最大面积三角形,显然三角形顶点具有单调性(同旋转卡壳)故旋转卡壳
写的时候。。一开始Cross函数左右傻傻分不清楚,然后,凸包里面的点是不能有重复的
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
const int maxn = 2E3 + 20;
typedef double DB;
const DB eps = 1E-10;
struct P{
DB x,y;
bool operator < (const P &b) const {
if (x < b.x) return 1;
if (x > b.x) return 0;
return y < b.y;
}
P operator - (const P &b) {
P ret; ret.x = x - b.x; ret.y = y - b.y;
return ret;
}
}p[maxn],q[maxn*2];
DB Cross(P v,P w) {return v.x*w.y - v.y*w.x;}
DB Dot(P v,P w) {return v.x*w.x + v.y*w.y;}
int n,cnt;
DB ans;
bool cmp(DB A)
{
if (A < 0) return 1;
if (abs(A) < eps) return 1;
return 0;
}
int main()
{
#ifdef YZY
freopen("land10.in","r",stdin);
#endif
cin >> n;
for (int i = 1; i <= n; i++) scanf("%lf%lf",&p[i].x,&p[i].y);
sort(p + 1,p + n + 1);
for (int i = 1; i <= n; i++) {
while (cnt > 1 && cmp(Cross(q[cnt]-q[cnt-1],p[i]-q[cnt-1]))) --cnt;
q[++cnt] = p[i];
}
int k = cnt;
for (int i = n - 1; i; i--) {
while (cnt > k && cmp(Cross(q[cnt]-q[cnt-1],p[i]-q[cnt-1]))) --cnt;
q[++cnt] = p[i];
}
for (int i = 0; i < cnt; i++) q[i] = q[i+1]; --cnt;
for (int x = 0; x < cnt; x++) {
int a,b; a = (x+1)%cnt; b = (x+3)%cnt;
for (int y = x + 2; y < cnt; y++) {
while ((a+1)%cnt != y && Cross(q[a]-q[x],q[y]-q[x]) < Cross(q[(a+1)%cnt]-q[x],q[y]-q[x]))
a = (a+1)%cnt;
while ((b+1)%cnt != x) {
DB A = Cross(q[y]-q[x],q[b]-q[x]);
DB B = Cross(q[y]-q[x],q[(b+1)%cnt]-q[x]);
if (A < B) b = (b+1)%cnt; else break;
}
ans = max(ans,Cross(q[a]-q[x],q[y]-q[x]) + Cross(q[y]-q[x],q[b]-q[x]));
}
}
printf("%.3f",ans/2);
return 0;
}