二维凸包

#include<bits/stdc++.h>
using namespace std;
#define N 1000010
const int mod = 1000000007;
const int inf = 1000000007;
const double eps = 1e-9;
void in(int &x){
    x=0;char c=getchar();
    int y=1;
    while(c<'0'||c>'9'){if(c=='-')y=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+c-'0';c=getchar();}
    x*=y;
}
void o(int x){
    if(x<0){x=-x;putchar('-');}
    if(x>9)o(x/10);
    putchar(x%10+'0');
}
int sgn(double x){
    if (x>eps)return 1;
    if (x<-eps)return -1;
    return 0;
}
// 向量
struct vec{
    double x,y;
    vec(){x=y=0;}
    vec(double _x,double _y){x=_x,y=_y;}
//    向量四则运算
    vec operator + (vec v){return vec(x+v.x,y+v.y);}
    vec operator - (vec v){return vec(x-v.x,y-v.y);}
    vec operator * (double v){return vec(x*v,y*v);}
    vec operator / (double v){return vec(x/v,y/v);}
    double operator * (vec v){return x*v.x+y*v.y;}
};
// 按照 x 为第一关键字, y 为 第二关键字排序
bool cmpXY(vec a,vec b){
    if (sgn(a.x-b.x))return a.x<b.x;
    return a.y<b.y;
}
// 叉乘
// 叉乘 大于 相同 表示 对于一个点来说 时针 顺逆相同
double cross(vec a ,vec b){
    return a.x*b.y-a.y*b.x;
}
// 计算凸包
int convec_hull(vec*v,int n,vec* z){
    // 先安装关键字排个序
    sort(v,v+n,cmpXY);
    // z 为手动模拟的一个栈
    int m = 0;
//
    for(int i=0;i<n;i++){
//  位于逆时针方向 且 栈非空 出栈
        while (m>1 && cross(z[m-1]-z[m-2],v[i]-z[m-2])<=0)--m;
//        入栈
        z[m++] = v[i];
    }
    int k = m;
//    从 n - 2 开始 因为 n-1 即队尾必在凸包中
    for(int i=n - 2;i>=0;--i) {
//        同理,判栈
        while( m > k && cross(z[m-1]-z[m-2],v[i]-z[m-2]) <= 0 ) --m;
//        入栈
        z[m++] = v[i];
    }
// 初始点重复入栈,退出
    if(n>1)--m;
    return m;
}
double Length(vec x){
    return sqrt(x.x*x.x+x.y*x.y);
}
int n;
vec a[N],ans[N];
signed main(){
    in(n);
    for(int i=0;i<n;i++){
        cin>>a[i].x>>a[i].y;
    }
    int m=convec_hull(a,n,ans);
    double ANS=0.0;
//    输出长度
    for(int i=0;i<m;i++){
        ANS += Length(ans[i+1]-ans[i]);
    }
    printf("%.2lf\n",ANS);
    return 0;
}
posted @ 2021-11-11 17:03  yesuweiYYYY  阅读(26)  评论(0编辑  收藏  举报