HDU 1392-Surround the Trees(求凸包周长)

---

求凸包比较容易是实现的一种方法：Graham扫描算法(基于平面扫描法)；

                                                 

基本过程：

1、寻找极点（x最小的点如果x相同则y最小的点）

2、极角排序（通过叉乘来判断其余点与极点连线的位置，按照逆时针顺序来排序）

3、用数组模拟栈来（因为此过程中需要栈顶和栈顶下一元素，so用数组比较方便获取）将p0,p1,p2三点放入栈中，之后依次遍历p3~pn各点，若栈顶的两个点和当前的点pi这三点连线的方向向顺时针方向偏转，表明栈顶的点是一个凹陷处，应删除，则栈顶元素出栈。并将pi压入栈中。如果向逆时针方向偏转，则直接压入栈中。 

4.重复上述过程，直到遍历完成，在栈中的点就是凸包上的顶点。（结果图是最上边的图）

---

 代码实现

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e3+10;
struct point{
	double x,y;
}p[maxn],stk[maxn];//用数组模拟栈 
double dis(point a,point b){//两点间距离
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int cross(point a,point b){//叉乘
	return a.x*b.y-a.y*b.x;
}
int crosst(point a,point b,point c){//ab x ac;
	return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
	
}
bool cmp(point a,point b){
	int flag = (a.x-p[1].x) * (b.y-p[1].y) - (b.x-p[1].x) * (a.y-p[1].y);
	if(flag > 0) return true; //逆时针排序
	if(flag == 0 && dis(a,p[1]) < dis(b,p[1])) //如果极角度数相同，按其到极点距离从大到小排序 
		return true; 
	return false;
}
double perimeter(int n){
	double ans = 0;
	for(int i=1;i<=n;i++){
		ans+=dis(stk[i],stk[i+1]);
	}
	return ans;
}
int main(){
	int n;
	while(cin>>n&&n){
		for(int i=1;i<=n;i++)
			cin>>p[i].x>>p[i].y;
		if(n==1) { printf("0\n");continue;}
		if(n==2) { printf("%.2lf\n",dis(p[1],p[2]));continue;}
        //1 2 特判
		int pole=1; 
		for(int i=2;i<=n;i++){//find pole
			if(p[i].x < p[pole].x || (p[i].x == p[pole].x) && p[i].y < p[pole].y) 
				pole=i;
		}
		swap(p[1],p[pole]);//将p[1]作为极点
		sort(p+2,p+n+1,cmp);//极角排序
		
		stk[1] = p[1]; stk[2] = p[2];
		int tot = 2;  
		for(int i = 3;i<=n;i++){//这个while有点迷
			while(tot >= 2 && crosst(stk[tot-1],stk[tot],p[i]) <= 0) tot--;
			stk[++tot]=p[i];					
		}//stk剩余元素则是凸包各顶点
		stk[tot+1]=p[1];//方便计算周长
		double ans = perimeter(tot);	
		printf("%.2lf\n",ans);
	}
	
	return 0;
} 

上述代码还可以继续完善，如重载运算符

#include<bits/stdc++.h>
using namespace std;

struct Point {
	double x, y;
	Point friend operator - (Point a, Point b) {
		return {a.x - b.x, a.y - b.y};
	}
} p[105], Stack[105];

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

double multi(Point p1, Point p2, Point p3) {
	return X(p2 - p1, p3 - p1);
}

double dis(Point a, Point b) {
	Point c = a - b;
	return sqrt(c.x * c.x + c.y * c.y);
}
int cmp(Point a, Point b) {
	double x = X(a - p[1], b - p[1]);

	if(x > 0) return 1;
	if(x == 0 && dis(a, p[1]) < dis(b, p[1])) return 1;

	return 0;
}
int main() {
	int N;
	while(scanf("%d", &N) && N) {
		for(int i = 1; i <= N; i++)
			scanf("%lf %lf", &p[i].x, &p[i].y);
		// N ==1or2情况特判
		if(N == 1) {
			printf("0.00\n");
			continue;
		} else if(N == 2) {
			printf("%.2f\n", dis(p[1], p[2]));
			continue;
		}
		int k = 1;
		for(int i = 2; i <= N; i++)
			if(p[i].y < p[k].y || (p[i].y == p[k].y && p[i].x < p[k].x))
				k = i;
		swap(p[1], p[k]);
		
		sort(p + 2, p + N + 1, cmp);
		Stack[1] = p[1];
		Stack[2] = p[2];
		int tot = 2;
		for(int i = 3; i <= N; i++) {
		    
			while(tot >= 2 && multi(Stack[tot - 1], Stack[tot], p[i]) <= 0) tot--;
			Stack[++tot] = p[i];
		}
		//便于计算
		Stack[tot + 1] = Stack[1];
		double sum = 0;
		for(int i = 1; i <= tot; i++) {
			sum += dis(Stack[i], Stack[i + 1]);
		}
		printf("%.2f\n", sum);
	}
	return 0;
}