HDU 1392 Surround the Trees 构造凸包

又是一道模板题

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <cmath>
#include <algorithm>

using namespace std;

struct P
{
    int x,y;
    double angle;
} p[50010];

bool cmp(P p1,P p2)
{
    return p1.angle > p2.angle;
}

double calangle(P p,P t)
{
    return ((t.x-p.x)*(t.x-p.x) + 1 - (t.x-p.x-1)*(t.x-p.x-1))/(2*sqrt((t.x-p.x)*(t.x-p.x) + (t.y-p.y)*(t.y-p.y)));
}

P sp[50010];

bool judgesite(P p1,P p2, P p)
{
    P l1,l2;
    l1.x = p2.x - p1.x;
    l1.y = p2.y - p1.y;

    l2.x = p.x - p1.x;
    l2.y = p.y - p1.y;

    return (l1.x*l2.y - l1.y*l2.x) > 0 ? true : false ;
}

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

int main()
{
    int n,i;

    int top;

    double MinLen;

    while(scanf("%d",&n) && n)
    {

        top = 0;

        for(i = 1; i <= n; ++i)
        {
            scanf("%d %d",&p[i].x,&p[i].y);
        }

        P temp;

        for(temp = p[1],i = 2; i <= n; ++i)
        {
            if(temp.y > p[i].y)
            {
                temp = p[i];
            }
            else if(temp.y == p[i].y && temp.x > p[i].x)
            {
                temp = p[i];
            }
        }

        for(i = 1; i <= n; ++i)
        {
            if(temp.x != p[i].x || temp.y != p[i].y)
            {
                p[i].angle = calangle(temp,p[i]);
            }
            else p[i].angle = -2;
        }

        sort(p+1,p+n+1,cmp);

        sp[top++] = p[n];
        sp[top++] = p[1];

        for(i = 2; i <= n;)
        {
            if(judgesite(sp[top-2],sp[top-1],p[i]))
            {
                sp[top++] = p[i];
                ++i;
            }
            else top--;
        }

        top--;//最后一个点重复了

        MinLen = callen(sp[top-1],sp[0]);

        for(i = 1; i < top; ++i)
        {
            MinLen += callen(sp[i-1],sp[i]);
        }

        printf("%.2lf
",MinLen);
    }
    return 0;
}