poj1113 Wall

题目链接：http://poj.org/problem?id=1113

题意：平面上有一些点，求包围这些点且距离最外层的点的距离为L的图形的周长，所谓最外层即这些点的凸包，如下图：

         

拐点处为包围拐点的半径为L的扇形，所有的拐点处包围图形的总长度正好为一个半径为L的圆形的周长，所以所求答案为凸包长度+2*Pi*L

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;

const double inf = 999999;
const double eps = 1e-8;
const double pi = 3.141592653;

struct Point {
    double x, y;
} p[1005];

int stk[1005], top;

int dblcmp(double k) {
    if (fabs(k) < eps) return 0;
    return k > 0 ? 1 : -1;
}

double multi(Point p0, Point p1, Point p2) {
    return (p1.x - p0.x) * (p2.y - p0.y) - (p1.y - p0.y) * (p2.x - p0.x);
}

double getDis(Point a, Point b) {
    return  sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

bool cmp(const Point& a, const Point& b) {
    double tmp = multi(p[0], a, b);
     if (!dblcmp(tmp)) return getDis(p[0], a) < getDis(p[0], b);
     return dblcmp(tmp) > 0;
}

int main()
{
    int n, l, k, i;
    double tx, ty, tmp, ans;

    scanf ("%d%d", &n, &l);
    tx = ty = inf;
    for (i = 0; i < n; i++) {
        scanf ("%lf%lf", &p[i].x, &p[i].y);
        tmp = p[i].y-ty;
        if (dblcmp(tmp) < 0) {
            ty = p[i].y;
            tx = p[i].x;
            k = i;
        }
        else if (!dblcmp(tmp) && dblcmp(p[i].x-tx) < 0) {
            tx = p[i].x;
            k = i;
        }
    }
    Point tp;
    tp.x = p[0].x; tp.y = p[0].y; //将最右下角的点交换到p[0]来
    p[0].x = p[k].x; p[0].y = p[k].y;
    p[k].x = tp.x; p[k].y = tp.y;
    sort(p+1, p+n, cmp); //极角排序
    
/*此处与黑书上有所不同，根据算法的思想，应当是让堆栈中存储一个临时凸包，然后根据新的点不断更新凸包，
  所以黑书上初始化时加入p[]的前三个点，堆栈中初始化为一个凸包，但是以下在初始化时只加入了两个点，那是因为
  黑书中的排序采用的是合并排序，排序过程中当极角相同时，只保留了最远点，那样就能保证p[]中的前三个点必定为
  一个凸包，但是此处我采用的排序只是普通快排，并未合并，所以p[]中的前三个点并不一定能构成一个凸包(三点共线时)，
  当然由于此题是求凸包周长，三点共线并不会影响最终结果，但是如果求的是凸包的顶点，则会有影响了
*/
    top = 0;
    stk[0] = 0;
    stk[++top] = 1;
    for (i = 2; i < n; i++) {
        while (top >= 1 && dblcmp(multi(p[stk[top-1]], p[i], p[stk[top]])) >= 0) top--;
        stk[++top] = i;
    }
    ans = 0;
    for (i = 1; i <= top; i++)
        ans += getDis(p[stk[i-1]], p[stk[i]]);
    ans += getDis(p[stk[0]], p[stk[top]]) + 2*pi*l;
    printf ("%.0lf\n", ans);
    return 0;
}