CF 50C Happy Farm 5题解

最新推荐文章于 2019-01-25 14:05:33 发布

gh546okmy

最新推荐文章于 2019-01-25 14:05:33 发布

阅读量469

点赞数

分类专栏：计算几何

本文链接：https://blog.csdn.net/gh546okmy/article/details/42030855

版权

计算几何专栏收录该内容

14 篇文章 0 订阅

订阅专栏

求个简单的凸包。。看坐标都是整数，用整点表示最后求出的凸包居然还有错，改成浮点数就a了，最后答案一定是从凸包起始点到终点的，每两个相邻点的x差与y差取最大值，

然后求和+4，+4的原因是你要绕着走，相当于转一圈。

还有一种简单方法是直接每枚举一个点，求（x+y，x-y，y-x，-y-x），取n个点中算出的这四个的最大值，最后求和+4，具体为什么不清楚，大神都是这么做的=。=

#include<stdio.h>
#include<iostream>
#include<math.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define eps 1e-8
int sig(double x) {return (x>eps)-(x<-eps);}
struct P
{
    double x,y;
    P(double a=0,double b=0):x(a),y(b){}
    P operator +(const P &a)const{
        return P(x+a.x,y+a.y);
    }
    P operator -(const P &a)const{
        return P(x-a.x,y-a.y);
    }
    double operator *(const P &a)const{
        return x*a.x+y*a.y;
    }
    double operator ^(const P &a)const{
        return x*a.y-y*a.x;
    }
    int operator <(const P &a)const{
        if(a.x!=x) return x<a.x;
        else return y<a.y;
    }
};
int Graham(P *p,int n,P *q)
{
    sort(p,p+n);
    p[n]=p[0];
    int m=0;
    for(int i=0;i<n;i++)
    {
        while(m>1&&sig((q[m-1]-q[m-2])^(p[i]-q[m-2]))<=0) m--;
        q[m++]=p[i];
    }
    int k=m;
    for(int i=n-2;i>=0;i--)
    {
        while(m>k&&sig((q[m-1]-q[m-2])^(p[i]-q[m-2]))<=0) m--;
        q[m++]=p[i];
    }
    if(n>1) m--;
//    for(int i=0;i<m;i++)
//       printf("%d %d\n",q[i].x,q[i].y);
    return m;
}
P p[100100],q[100100];
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<n;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        int m=Graham(p,n,q);
        if(m==1)
        {
            printf("4\n");
            continue;
        }
        else if(m==2)
        {
            int sum=max(fabs(q[0].x-q[1].x),fabs(q[0].y-q[1].y))*2+4;
            printf("%d\n",sum);
            continue;
        }
        double ans=0;
        q[m]=q[0];
        for(int i=0;i<m;i++)
        {
            ans+=max(fabs(q[i].x-q[i+1].x),fabs(q[i].y-q[i+1].y));
            //cout<<ans<<endl;
        }
        printf("%.lf\n",ans+4);
    }
    return 0;
}