给出一些树的坐标,这些树能围成一片区域,每50单位的区域可以养一头牛,求最多能养多少头牛。
就是最大区域除以50。
根据凸包的定义,这些点形成的最大面积当然就是凸包的面积咯。
先求出凸包,然后用叉乘求出面积。
AC代码:
#include<iostream>
#include<algorithm>
using namespace std;
struct Point{
int x,y;
Point(){} //构造函数
Point(int xx,int yy){
x = xx;
y = yy;
}
Point operator - (const Point a)const{
return Point(x-a.x,y-a.y);
}
}p[11000],ch[11000]; //元序列和凸包序列
bool cmp(Point a,Point b){ //点排序
if(a.x == b.x) return a.y < b.y;
return a.x < b.x;
}
double det (Point a, Point b) { //返回叉乘的摸 亦可以用重载*的方法
return a.x * b.y - a.y * b.x;
}
int ConvexHull ( Point * p,int n, Point * ch ){ //构造凸包
sort(p,p+n,cmp);
int m=0;
for (int i=0;i<n;i++){
while (m >1&& det(ch[m -1] - ch[m -2] , //回溯
p[i]-ch[m -2]) <=0) m --;
ch[m++]= p[i];
}
int k=m;
for (int i=n-2;i >=0;i--){
while (m>k && det(ch[m -1] - ch[m -2] ,
p[i]-ch[m -2]) <=0) m --;
ch[m ++]= p[i];
}
if(n >1) m--;
return m;
}
int main(void){
int n,m,ans;
double area = 0;
cin>>n;
for(int i=0;i<n;i++) cin>>p[i].x>>p[i].y;
m = ConvexHull(p,n,ch);
for(int i=0;i<m;i++){
area += ch[i].x*ch[i+1%m].y - ch[i+1%m].x*ch[i].y;
}
area /= 2;
if(area <= 0 ) area = -area;
ans = area/50;
cout<<ans<<endl;
return 0;
}