BZOJ 1069: [SCOI2007]最大土地面积 旋转卡壳

1069: [SCOI2007]最大土地面积

Time Limit: 1 Sec  Memory Limit: 128 MB
Submit: 3680  Solved: 1460
[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

Sample Output

1.000

HINT

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

直接枚举对踵点 然后枚举左右两个点

苦苦调试好久 发现叉积竟然写错了。。
#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;

typedef double db;
typedef long long ll;

inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=2010,inf=0X3f3f3f3f;
const db eps=1e-8;

struct node
{
db x,y;

friend db cross(const node &x,const node &y,const node &bas)
{
db x1=x.x-bas.x,x2=y.x-bas.x;
db y1=x.y-bas.y,y2=y.y-bas.y;
return x1*y2-x2*y1;
}
friend db dis(const node &x,const node &y)
{return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));}

}p
;

int n;

inline bool cmp(const node &x,const node &y)
{
db tmp=cross(x,y,p[1]);
return tmp==0 ? dis(x,p[1])<dis(y,p[1]) : tmp<0;
}

int st
,top;

void graham()
{
register int i,tmp=1;
for(i=1;i<=n;++i)
if(p[i].x<p[tmp].x || (p[i].x==p[tmp].x&&p[i].y<p[tmp].y))
tmp=i;
swap(p[1],p[tmp]);
sort(p+2,p+n+1,cmp);
st[1]=1;st[2]=2;top=2;
for(i=3;i<=n;++i)
{
while(top>1&&cross(p[st[top]],p[i],p[st[top-1]])>-eps)
top--;
st[++top]=i;
}
st[++top]=1;
}

db ans;

void RC()
{
register int i,now=2,j;
for(i=1;i<top;++i)
{
while( cross(p[st[now+1]],p[st[i+1]],p[st[i]])>cross(p[st[now]],p[st[i+1]],p[st[i]]) )
{
now++;
if(now==top)now=1;
}
db tmp[2];
tmp[0]=tmp[1]=0;
j=i;
while(j!=now)
{
j++;
if(j==top)j=1;
tmp[0]=max( tmp[0],abs(cross(p[st[i]],p[st[now]],p[st[j]])) );
}
while(j!=i)
{
j++;
if(j==top)j=1;
tmp[1]=max( tmp[1],abs(cross(p[st[i]],p[st[now]],p[st[j]])) );
}
ans=max(ans,(tmp[0]+tmp[1])*0.5);
}
}

int main()
{
n=read();
register int i;
for(i=1;i<=n;++i)
scanf("%lf%lf",&p[i].x,&p[i].y);
graham();
RC();
printf("%.3lf\n",ans);
}