学习凸包（四）：Graham 扫描法

Graham扫描法

基本思想：通过设置一个关于候选点的堆栈来解决凸包问题。

操作：输入集合P中的每一个点都被压入栈一次，非凸包中的顶点的点最终将被弹出堆栈，
当算法终止时，堆栈中仅包含凸包中的顶点，其顺序为个各顶点在边界上出现的逆时针方向排列的顺序。


(1)设P0是P中Y坐标最小的点，如果有多个这样的点则取最左边的点作为P0；

(2) 设<P1,P2,……,Pn >是P中剩余的点，对其按逆时针方向相对P0 的极角进行排序，
如果有数个点有相同的极角，则去掉其余的点，只留下一个与P0 距离最远的那个点；

(3)
//前三个点先入栈
ch[0] = p[0];
ch[1] = p[1];
ch[2] = p[2];

//判断与其余所有点的关系
for (int i = 3; i < n; i++) {
//不满足向左转的关系,栈顶元素出栈
while (top > 0 && multiply(p[i], ch[top], ch[top - 1]) >= 0)
top--;
//当前点与栈内所有点满足向左关系,因此入栈.
ch[++top] = p[i];
}

原理：沿逆时针方向通过凸包时，在每个顶点处应该向左转。因此，while循环每次发现在一个顶点处没有向左转时，就把该顶点从堆栈中弹出。）当算法向点pi推进、在已经弹出所有非左转的顶点后，就把pi压入堆栈中。


下面是POJ1113的AC代码：关于POJ1113请参见 http://128kj.iteye.com/blog/1748635

Java代码

1. import java.util.Scanner;
3. class Point {
4. double x;
5. double y;
7. public Point(int x, int y) {
8. this.x = x;
9. this.y = y;
10. }
11. }
13. public class Main {
14. Point[] ch; //点集p的凸包
15. Point[] p ; //给出的点集
16. int n;
17. int l;
18. int len=0;
20. public Main(Point[] p,int n,int l){
21. this.p=p;
22. this.n=n;
23. this.l=l;
24. ch= new Point[n];
25. }
27. //小于0,说明向量p0p1的极角大于p0p2的极角
28. public double multiply(Point p1, Point p2, Point p0) {
29. return ((p1.x - p0.x) * (p2.y - p0.y) - (p2.x - p0.x) * (p1.y - p0.y));
30. }
32. //求距离
33. public double distance(Point p1, Point p2) {
34. return (Math.sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y)
35. * (p1.y - p2.y)));
36. }