【POJ】3348 - Cows(经典凸包)

给出一些树的坐标,这些树能围成一片区域,每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;
}

猜你喜欢

转载自blog.csdn.net/qq_41009682/article/details/81268501
今日推荐