题解报告：hdu 1392 Surround the Trees（凸包入门）

Problem Description

There are a lot of trees in an area. A peasant wants to buy a rope to surround all these trees. So at first he must know the minimal required length of the rope. However, he does not know how to calculate it. Can you help him?
The diameter and length of the trees are omitted, which means a tree can be seen as a point. The thickness of the rope is also omitted which means a rope can be seen as a line.

There are no more than 100 trees.

Input

The input contains one or more data sets. At first line of each input data set is number of trees in this data set, it is followed by series of coordinates of the trees. Each coordinate is a positive integer pair, and each integer is less than 32767. Each pair is separated by blank.
Zero at line for number of trees terminates the input for your program.

Output

The minimal length of the rope. The precision should be 10^-2.

Sample Input

12 7

24 9

30 5

41 9

80 7

50 87

22 9

45 1

50 7

Sample Output

243.06

解题思路：入门凸包！有一篇博文讲得非常好非常易懂，链接：凸包详解。再贴一张用叉积来判断栈顶的点是否为凸包上合法的顶点的依据。时间复杂度主要用在排序上为O(nlogn)。

以下是Graham算法的扫描过程：

AC代码：

 1 #include<iostream>
 2 #include<string.h>
 3 #include<algorithm>
 4 #include<cstdio>
 5 #include<cmath>
 6 using namespace std;
 7 const int maxn=105;
 8 const double PI=acos(-1.0);
 9 struct node{int x,y;};
10 node vex[maxn];//存入所有坐标点
11 node stackk[maxn];//凸包中所有的点
12 bool cmp1(node a,node b){//按点的坐标排序
13     if(a.y==b.y)return a.x<b.x;//如果纵坐标相同，则按横坐标升序排
14     else return a.y<b.y;//否则按纵坐标升序排
15 }
16 bool cmp2(node a,node b){//以基点为坐标原点，极角按升序排，这里可用atan2函数或者叉积来进行极角排序，但是用atan2函数来排序效率高时间快，不过精度比叉积低
17     double A=atan2(a.y-stackk[0].y,a.x-stackk[0].x);//返回的是原点至点(x,y)的方位角，即与x轴的夹角
18     double B=atan2(b.y-stackk[0].y,b.x-stackk[0].x);
19     if(A!=B)return A<B;//逆时针方向为正值，极角小的排在前面
20     else return a.x<b.x;//如果极角相同，则横坐标在前面的靠前排列
21 }
22 int cross(node p0,node p1,node p2){//计算两个向量a、b（a=(x1,y1),b=(x2,y2)）的叉积公式：a×b=x1y2-x2y1 ===> p0p1=(x1-x0,y1-y0),p0p2=(x2-x0,y2-y0)
23     return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
24 }
25 double dis(node a,node b){//计算两点之间的距离
26     return sqrt((a.x-b.x)*(a.x-b.x)*1.0+(a.y-b.y)*(a.y-b.y));
27 }
28 int main(){
29     int t;
30     while(~scanf("%d",&t)&&t){
31         for(int i=0;i<t;++i)//输入t个点
32             scanf("%d%d",&vex[i].x,&vex[i].y);
33         if(t==1)printf("%.2f\n",0.00);//如果只有一个点，则周长为0.00
34         else if(t==2)printf("%.2f\n",dis(vex[0],vex[1]));//如果只有两个点，则周长为两个点的距离
35         else{
36             memset(stackk,0,sizeof(stackk));//清0
37             sort(vex,vex+t,cmp1);//先按坐标点的位置进行排序
38             stackk[0]=vex[0];//取出基点
39             sort(vex+1,vex+t,cmp2);//将剩下的坐标点按极角进行排序，以基点为坐标原点
40             stackk[1]=vex[1];//将凸包中的第二个点存入凸集中
41             int top=1;//当前凸包中拥有点的个数为top+1
42             for(int i=2;i<t;++i){//不断地找外围的坐标点
43                 while(top>0&&cross(stackk[top-1],stackk[top],vex[i])<0)top--;//如果叉积为负数，则弹出栈顶元素
44                 //虽然第2个凸点显然是最外围的一点，但加上top>0保证了栈中至少有2个凸点
45                 stackk[++top]=vex[i];
46             }
47             double s=0;
48             for(int i=1;i<=top;++i)//计算凸包的周长
49                 s+=dis(stackk[i-1],stackk[i]);
50             s+=dis(stackk[top],vex[0]);//最后一个点和第一个点之间的距离
51             printf("%.2f\n",s);
52         }
53     }
54     return 0;
55 }