http://poj.org/problem?id=2079
题意:求平面点集最大三角形面积。
思路和求平面最远点对类似,先固定一个点 p [ i ] p[i] p[i],然后根据单调性旋转 p [ j ] p[j] p[j], p [ k ] p[k] p[k]记录一个点为 p [ i ] p[i] p[i]的三角形的最大面积,再枚举下一个 p [ i ] p[i] p[i]直到找到最大三角形面积,而确定三点后的三角形面积可通过叉积方便求出。
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cmath>
using namespace std;
struct Point
{
int x,y;
Point(){}
Point(int _x,int _y)
{
x=_x;y=_y;
}
Point operator -(const Point &b)const
{
return Point(x-b.x,y-b.y);
}
int operator ^(const Point &b)const
{
return x*b.y-y*b.x;
}
int operator *(const Point &b) const
{
return x*b.x+y*b.y;
}
void input()
{
scanf("%d%d",&x,&y);
}
};
//距离的平方
int dist2(Point a,Point b)
{
return (a-b)*(a-b);
}
//二维凸包(int)
const int MAXN=50010;
Point list[MAXN];
int Stack[MAXN],top;
//相对于list[0]进行极角排序
bool _cmp(Point p1,Point p2)
{
double tmp=(p1-list[0])^(p2-list[0]);
if(tmp>0)
return true;
else if(tmp==0&&dist2(p1,list[0])<=dist2(p2,list[0]))
return true;
else return false;
}
void Graham(int n)
{
Point p0;
int k=0;
p0=list[0];
//找最左下的点
for(int i=1;i<n;i++)
{
if((p0.y>list[i].y)||(p0.y==list[i].y&&p0.x>list[i].x))
{
p0=list[i];
k=i;
}
}
swap(list[k],list[0]);
sort(list+1,list+n,_cmp);
if(n==1)
{
top=1;
Stack[0]=0;
return ;
}
if(n==2)
{
top=2;
Stack[0]=0;
Stack[1]=1;
return ;
}
Stack[0]=0;
Stack[1]=1;
top=2;
for(int i=2;i<n;i++)
{
while(top>1&&((list[Stack[top-1]]-list[Stack[top-2]])^(list[i]-list[Stack[top-2]]))<=0)
top--;
Stack[top++]=i;
}
}
//旋转卡壳,求两点间距离的平方的最大值
int rotating_calipers(Point p[],int n)
{
int ans=0;
for(int i=0;i<n;i++)
{
int j=(i+1)%n;
int k=(j+1)%n;
while(j!=i&&k!=i)
{
ans=max(ans,abs((p[j]-p[i])^(p[k]-p[i])));
while(((p[i]-p[j])^(p[(k+1)%n]-p[k]))<0)
k=(k+1)%n;
j=(j+1)%n;
}
}
return ans;
}
Point p[MAXN];
int main()
{
int n;
while(1)
{
scanf("%d",&n);
if(n==-1) break;
for(int i=0;i<n;i++)
list[i].input();
Graham(n);
for(int i=0;i<top;i++)
p[i]=list[Stack[i]];
printf("%.2f\n",(double)rotating_calipers(p,top)/2);
}
return 0;
}