[USACO5.1]Fencing the Cows

最新推荐文章于 2025-01-21 17:30:14 发布

原创最新推荐文章于 2025-01-21 17:30:14 发布 · 309 阅读

0 ·

CC 4.0 BY-SA版权

usaco 同时被 2 个专栏收录

5 篇文章

订阅专栏

凸包

0 篇文章

订阅专栏

本文介绍了一种计算二维平面上点集凸包的算法实现。通过定义特定的数据结构和比较函数，采用排序和扫描策略来构造凸包。该实现包括点的定义、比较函数、插入和删除操作等关键步骤，并最终计算并输出凸包周长。

裸的二维凸包

namespace point
{
    struct point
    {
        double x,y;
        point()
        {
            x=y=0;
        }
        point makepoint(double x,double y)
        {
            point k;
            k.x=x;
            k.y=y;
            return k;
        }
        point operator - (point a)
        {
            return makepoint(x-a.x,y-a.y);
        }
    }g[10010];
    int len,s[10010];
    bool cmp1(point a,point b)
    {
        return a.y<b.y;
    }
    bool cmp2(point a,point b)
    {
        return a.x*b.y<a.y*b.x;
    }
    bool cmp3(point c,point b,point a)
    {
        return !cmp2(a-c,b-c);
    }
    void insert(int x)
    {
        s[++len]=x;
    }
    void popfrom(int x)
    {
        while(len>2&&cmp3(g[x],g[s[len]],g[s[len-1]]))
            len--;
    }
    double sqr(double x)
    {
        return x*x;
    }
    double dis(point a,point b)
    {
        return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
    }
    double calc()
    {
        double ans=dis(g[s[1]],g[s[len]]);
        fr(i,2,len)
            ans+=dis(g[s[i]],g[s[i-1]]);
        return ans;
    }
}
int n;
int main()
{
    n=read();
    point::point *a=point::g;
    fr(i,1,n)
        scanf("%lf%lf",&(a+i)->x,&(a+i)->y);
    sort(a+1,a+n+1,point::cmp1);
    fd(i,n,1)
        a[i]=a[i]-a[1];
    sort(a+2,a+n+1,point::cmp2);
    fr(i,1,2)
        point::insert(i);
    fr(i,3,n)
    {
        point::popfrom(i);
        point::insert(i);
    }
    printf("%.2lf\n",point::calc());
    return 0;
}