HDU 1392 求凸包

最新推荐文章于 2022-01-29 11:58:16 发布

lentty

最新推荐文章于 2022-01-29 11:58:16 发布

阅读量377

点赞数

分类专栏：大学ACM记忆文章标签： ini c

本文链接：https://blog.csdn.net/lentty1452/article/details/7782298

版权

大学ACM记忆专栏收录该内容

85 篇文章 0 订阅

订阅专栏

凸包在计算几何中的位置相当重要。一定要掌握，关于凸包模板，网上的看着太。。。，还是写了个自己能理解的SORT，还是不习惯QSORT

这里求凸包用的是graham扫描法

#include<stdio.h>
#include<iostream>
#include<math.h>
#include<algorithm>
#define maxn 100010
using namespace std;

struct point
{
  double x,y;
}p[maxn],stack[maxn];

double dis(point p1,point p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}

double area(point a,point b,point c)
{ 
   return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}

int cmp(point a,point b)
{
    if(area(p[0],a,b)<0)
    return 1;
    if(area(p[0],a,b)==0&&dis(a,p[0])>dis(b,p[0]))
    return 1;
    return 0;
}

int graham(point p[],point stack[],int n)
{
    int i,k=0;
    int top=1;
    point tmp;
    for(i=1;i<n;i++)
    {
       if(p[i].y<p[k].y||(p[i].y==p[k].y&&p[i].x<p[k].x))
       k=i;
    }
    tmp=p[0];
    p[0]=p[k];
    p[k]=tmp;
    sort(p+1,p+n,cmp);
    stack[0]=p[0];
    stack[1]=p[1];
    stack[2]=p[2];
    for(i=2;i<n;i++)
    {
      while(top&&area(stack[top-1],stack[top],p[i])>=0)
      top--;
      stack[++top]=p[i];
    }
    return top;
}

int main()
{
    int n,i,k;
    double sum;
    while(scanf("%d",&n)==1&&n)
    {
       sum=0.0;
       for(i=0;i<n;i++)
       scanf("%lf%lf",&p[i].x,&p[i].y);
       if(n==1)
       printf("0.00\n");
       else if(n==2)
       printf("%.2lf\n",dis(p[0],p[1]));
       else
       {
           k=graham(p,stack,n);
           stack[k+1]=stack[0];
           for(i=0;i<=k;i++)
             sum+=dis(stack[i],stack[i+1]);
           printf("%.2lf\n",sum);
       }
    }
    return 0;
}