[USACO5.1]二维凸包模板

传送门

Description

求\(n\)个点凸包的周长

Solution 

计算几何打暴力必备

Code 

#include<bits/stdc++.h>
#define reg register
#define ll long long
#define db double
using namespace std;

const int MN=1e5+5;
const db eps=1e-8;
int n;
struct Point{
    db x,y;
    Point(db x=0,db y=0):x(x),y(y){}
    db mo(){return sqrt(1.*x*x+1.*y*y);}
    Point operator-(const Point&o)const{return Point(x-o.x,y-o.y);}
    db operator^(const Point&o)const{return x*o.y-y*o.x;}
    bool operator <(const Point&o)const{return x<o.x||(x==o.x&&y<o.y);}
}p[MN];
int u[MN],d[MN],tu,td;

int main()
{
    scanf("%d",&n);
    reg int i;
    for(i=1;i<=n;++i) scanf("%lf%lf",&p[i].x,&p[i].y);
    std::sort(p+1,p+n+1);
    for(i=1;i<=n;++i)
    {
        while(tu>1&&((p[u[tu]]-p[u[tu-1]])^(p[i]-p[u[tu-1]]))>-eps) --tu;
        u[++tu]=i;
    }
    for(i=n;i;--i)
    {
        while(td>1&&((p[d[td]]-p[d[td-1]])^(p[i]-p[d[td-1]]))>-eps) --td;
        d[++td]=i;
    }
    db ans=0.;
    for(i=1;i<tu;++i) ans+=(p[u[i+1]]-p[u[i]]).mo();
    for(i=1;i<td;++i) ans+=(p[d[i+1]]-p[d[i]]).mo();
    return 0*printf("%.2lf\n",ans); 
}

---

Blog来自PaperCloud，未经允许，请勿转载，TKS！