凸包模板

2021-04-09 23:02:00 阅读：128 来源： 互联网

标签：return int double cross 凸包 -- && 模板


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;}
};
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);
    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;
    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;
}

板子题：
P2742
给n 个点，求形成的凸包的周长。
样例1

输出
12.00

#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;}
};
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);
    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;
    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;
}

标签：return,int,double,cross,凸包,--,&&,模板
来源： https://www.cnblogs.com/yesuweiYYYY/p/14639193.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

凸包模板