Scan法求凸包

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1348

给一个半径和n个点

求圆的周长 + n个点的凸包的周长

#include<bits/stdc++.h>
using namespace std;
const int maxn = ;
const double pi = acos(-1.0);
struct node {
    double x,y;
}p[maxn],P[maxn];
int n,tot;
double l,ans;

//向量AB 和 AC求x积 如果X(A,B,C)>0 则AC在AB的左边
double X(node A,node B,node C) {
    return (B.x-A.x)*(C.y-A.y)-(B.y-A.y)*(C.x-A.x);
}
double len(node A,node B) {
    return sqrt((B.x-A.x)*(B.x-A.x)+(B.y-A.y)*(B.y-A.y));
}
bool cmp(node A,node B) {
    double pp = X(p[],A,B);
    if(pp>) return true;
    if(pp<) return false;
    return len(p[],A) < len(p[],B);
}

void solve() {
    for(int i=;i<n;i++) {
        if(p[i].y < p[].y) swap(p[],p[i]);
        else if(p[i].y==p[].y && p[i].x < p[].x)
            swap(p[],p[]);
    }
    sort(p+,p+n,cmp);
    P[]=p[];
    P[]=p[];
    tot=;
    for(int i=;i<n;i++) {
        while (tot> && X(P[tot-],P[tot],p[i])<=) tot--;
        tot++;
        P[tot]=p[i];
    }
}

int main () {
    int T; cin >> T;
    for(int cas=;cas<=T;cas++) {
        if(cas!=) puts("");
        scanf("%d%lf",&n,&l);
        ans = *pi*l;
        for(int i=;i<n;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        if(n==)
            printf("%.0f\n",ans);
        if(n==)
            printf("%.0f\n",ans+len(p[],p[]));
        else {
            solve();
            for(int i=;i<tot;i++) {
                ans += len(P[i],P[i+]);
            }
            ans += len(P[],P[tot]);
            printf("%.0f\n",ans);
        }

    }
    return ;
}