[hdu 1348] Wall · 凸包

大概就是，给你n个点和一个值L，求这n个点的凸包，然后求与凸包相距l的外圈的周长。

结果=凸包的周长+半径的L的圆的周长

求凸包：

首先将n个按x坐标排序，然后选第一个点作为一个确定的点（最左边的点肯定在凸包上）

然后分两步，先求上凸包，再求下凸包。

求上凸包时，我们将目前已经在凸包里的点压入栈中，当做到第i号点的时候，如果发现

向量<top-1,i>在向量<top-1,top>的逆时针方向，那么说明top是凸包内的点而不是凸包上的点（现在求的是上凸包），所以将top从栈中弹出。

求下凸包的时候可以当成翻过来求上凸包，从后往前再扫一遍就行啦。

判断两个向量的相对方向可以用叉积，设有两个向量A(x1,y1),B(x2,y2)，则A×B=x1y2-x2y1，

若B在A的逆时针方向，则A×B>0；

若B在A的顺时针方向，则A×B<0。

#include <cstdio>#include <algorithm>#include <cstring>#include <cmath>using namespace std;#define pi 3.1415926535#define sqr(x) ((x)*(x))const int N=1005;struct arr{int x,y;}a[N];int stack[N],top;int T,n,l;double ans;bool cmp(const arr A,const arr B){return A.x<B.x || A.x==B.x && A.y<B.y;} double cha(int i,int j,int k){int x1,x2,y1,y2;x1=a[j].x-a[i].x;y1=a[j].y-a[i].y;x2=a[k].x-a[i].x;y2=a[k].y-a[i].y;return x1*y2-x2*y1;}void CC(){top=0;stack[++top]=1;//求上凸包 for (int i=2;i<=n;i++){while (top>1 && cha(stack[top-1],stack[top],i)>=0) top--;stack[++top]=i;}int tmp=top;//求下凸包  for (int i=n-1;i>=1;i--){while (top>tmp && cha(stack[top-1],stack[top],i)>=0) top--;stack[++top]=i;}ans=2*pi*l;for (int i=1;i<top;i++){int p,q;p=stack[i];q=stack[i+1];ans+=sqrt(sqr(a[p].x-a[q].x)+sqr(a[p].y-a[q].y));}ans+=sqrt(sqr(a[stack[top]].x-a[1].x)+sqr(a[stack[top]].y-a[1].y));}int main(){scanf("%d",&T);while (T--){scanf("%d%d",&n,&l);for (int i=1;i<=n;i++)scanf("%d%d",&a[i].x,&a[i].y);sort(a+1,a+n+1,cmp);CC();printf("%.0lf\n",ans);if (T) puts("");//注意格式。。。。。囧 }return 0;}