hdu-1392(二维凸包模板)

#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
typedef struct
{
    double x;
    double y;
} POINT;
POINT result[102];//保存凸包上的点,相当于所说的栈S
POINT a[102];
int n,top;
double Distance(POINT p1,POINT p2)//两点间的距离
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double Multiply(POINT p1,POINT p2,POINT p3)//叉积
{
    return((p2.x-p1.x)*(p3.y-p1.y)-(p2.y-p1.y)*(p3.x-p1.x));
}
bool cmp(POINT aa,POINT bb)//按点乘排序,cos=a*b/|a||b|
{
    int ax=aa.x-a[0].x;
    int ay=aa.y-a[0].y;
    int bx=bb.x-a[0].x;
    int by=bb.y-a[0].y;
    float cosa=ax/sqrt(ax*ax+ay*ay);
    float cosb=bx/sqrt(bx*bx+by*by);
    return cosa>cosb;
}
//寻找凸包的过程,p0,p1,p2..的寻找过程在下面main中进行了
void Tubao()
{
    int i;
    result[0].x=a[0].x;
    result[0].y=a[0].y;
    result[1].x=a[1].x;
    result[1].y=a[1].y;
    result[2].x=a[2].x;
    result[2].y=a[2].y;
    top=2;
    for(i=3; i<=n; i++)
    {
        while(Multiply(result[top-1],result[top],a[i])<=0&&top>=2)
            top--;
        result[top+1].x=a[i].x;
        result[top+1].y=a[i].y;
        top++;
    }
}
int main()
{
    int i,p;
    double px,py,len,temp;
    while(scanf("%d",&n)!=EOF,n)//
    {
        for(i=0; i<n; i++)
            scanf("%lf%lf",&a[i].x,&a[i].y);
        if(n==1)
        {
            printf("0.00\n");
            continue;
        }
        else if(n==2)
        {
            printf("%.2lf\n",Distance(a[0],a[1]));
            continue;
        }
//这里的目的好像是找出y坐标最小的点,之后把他定义为P0
        py=-1;
        for(i=0; i<n; i++)
        {
            if(py==-1||a[i].y<py)
            {
                px=a[i].x;
                py=a[i].y;
                p=i;
            }
            else if(a[i].y==py&&a[i].x<px)
            {
                px=a[i].x;
                py=a[i].y;
                p=i;
            }
        }
//swap(a[0],a[p])
        temp=a[0].x;
        a[0].x=a[p].x;
        a[p].x=temp;
        temp=a[0].y;
        a[0].y=a[p].y;
        a[p].y=temp;
        sort(a+1,a+n,cmp);
        a[n].x=a[0].x;
        a[n].y=a[0].y;
//调用tubao()
        Tubao();
        len=0.0;
        for(i=0; i<top; i++)
            len=len+Distance(result[i],result[i+1]);
        printf("%.2lf\n",len);
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值