bzoj1670: [Usaco2006 Oct]Building the Moat护城河的挖掘计算几何凸包

最新推荐文章于 2018-12-19 14:56:00 发布

原创最新推荐文章于 2018-12-19 14:56:00 发布 · 620 阅读

0 ·

CC 4.0 BY-SA版权

本文介绍了一种裸凸包算法的实现方式，并提供了一个简洁的C++代码模板。该算法首先寻找最低点作为起始点，然后通过极角排序其余点，并使用单调栈来去除那些不会出现在最终凸包上的点，最后计算出构成凸包的边的总长度。

裸凸包，上个模板。

#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
#define maxn 10000
struct node
{
    double x,y;
}p[maxn];
int s[maxn],top,n;
double ans;
double sqr(double x)
{
    return x*x;
}
double dis(node a,node b)
{
    return sqr(a.x-b.x)+sqr(a.y-b.y);
}
inline node operator-(node a,node b)
{
    node t;t.x=a.x-b.x;t.y=a.y-b.y;return t;
}
inline double operator*(node a,node b)
{
    return a.x*b.y-a.y*b.x;
}
bool cmp(node a,node b)
{
    double t=(a-p[1])*(b-p[1]);
    if(t==0)return dis(p[1],a)<dis(p[1],b);
    return t>0;
}
void graham()
{
    int t=1;
    for(int i=2;i<=n;i++)
        if(p[i].y<p[t].y||(p[i].y==p[t].y&&p[i].x<p[t].x))t=i;
    swap(p[1],p[t]);
    sort(p+2,p+n+1,cmp);
    top=2;
    s[1]=1;s[2]=2;
    for(int i=3;i<=n;i++)
    {
        for(;(p[s[top]]-p[s[top-1]])*(p[i]-p[s[top-1]])<=0;top--);
        top++;s[top]=i;
    }
    s[top+1]=1;
    for(int i=1;i<=top;i++)
        ans+=sqrt(dis(p[s[i]],p[s[i+1]]));
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%lf",&p[i].x,&p[i].y);
    }
    graham();
    printf("%.2lf\n",ans);
    return 0;
}