[WIKIOI1298]凸包周长[裸凸包]

http://www.wikioi.com/problem/1298/ 这是题目。

是一道裸的凸包算法题，但是我一开始还是写WA了。先贴代码

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#include <queue>
#include <iostream>
#define mod 9999973
#define N 100005
using namespace std;
int n , top;
double ans;
struct point
{
    int x , y;
    bool operator < (const point & a) const
    {
        if(x != a.x) return x < a.x;
        return y < a.y;
    }
}e[N] , sta[N];
double get(point p1 , point p2 , point p3)
{
    return (p1.x - p3.x) * (p2.y - p3.y) - (p2.x - p3.x) * (p1.y - p3.y);
}
double dist(point p1 , point p2)
{
    return sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y) * (p1.y - p2.y));
}
void work()
{
    scanf("%d", &n);
    for(int i = 1 ; i <= n ; ++i) scanf("%d%d", &e[i].x, &e[i].y);
    sort(e + 1 , e + 1 + n);
    sta[++top] = e[1] , sta[++top] = e[2];
    for(int i = 3 ; i <= n ; ++i)
    {
        while(top >= 2 && get(sta[top] , e[i] , sta[top - 1]) <= 0) --top;
        sta[++top] = e[i];
    }
    int len = top;
    sta[++top] = e[n - 1];
    for(int i = n - 2 ; i >= 1 ; --i)
    {
        while(len != top && get(sta[top] , e[i] , sta[top - 1]) <= 0) --top;
        sta[++top] = e[i];
    }
    ans += dist(sta[1] , sta[top]);
    for(int i = 1 ; i < top ; ++i) ans += dist(sta[i] , sta[i + 1]);
    printf("%.1f", ans);
}
int main()
{
    work();
    return 0;
}

错误原因：
1、一开始把求叉积的代码写挫了，return (p1.x - p3.x) * (p2.y - p3.y) - (p2.x - p3.x) * (p1.y - p3.y);
把(p2.y - p3.y) 写成了(p1.x - p3.x)

2、在while循环退栈的过程中，叉积>=0都可以退，不仅仅只是>0，这是关键错误，改了这个程序就过了。

求凸包的方法：

先将所有的点按照x大小排序，x相等就用y来排序。

然后从第三个点开始往后扫，扫到某个点i，求叉积，判断左右。。

第一遍从前往后扫出下凸壳，然后从后往前弄出上凸壳，再求周长就可以了

然后先将第一个点和第二个点入栈(这个很好理解)