凸包（convex hull）

凸包（convex hull）

凸包是指对于平面上给定的一些点，选取一个最小的凸多边形使得这些点或者在其内部，或者在其边上。

求凸包的Graham扫描法 

 

针对一个有三个或以上点的点集Q 

 

基本算法：
1. 已知各点坐标，将他们放入一个二维坐标系
2. 算出各点极角，按降序排序 ， 取最小点做起点
3. 三点扫描一次：
    if        对于不符合的点删除，然后继续扫描当前三点
    else    符合的点就找下一个点继续，直到扫描了所有的点

 

代码分析：
1. 建立一个点结构体point， 包括x，y坐标，和angle表示极角
2. 考虑到要排序，将初始的点用multiset 保存
3. 保存处理的点用deque（有的方法是用stack，因为操作中有类似出栈入栈的操作。但是判断点的过程是要访问到两个点，就是说不仅仅是要访问栈顶元素，而且还要访问栈顶的下一个元素，个人觉得不符合stack 的定义，所以我就用了deque）
4. 最后结果是保存在deque中 ， deque中的点元素就是组成所求凸包的点 

 

代码：

#include <iostream> #include <deque> #include <cmath> #include <set> using namespace std; #define PI 3.141592654 const double P = PI/180.0; struct point { double x; double y; double angle; point(double x1, double y1):x(x1) , y(y1) {angle = atan (x/y);} bool operator <(const point& other) const //插入SET 这个是关键 , 指定排序规则 { return angle < other.angle; } }; //输入点求极角并且排序 void Input (multiset<point> &p , int &number) { cin>>number; double x, y; for (int i = 0 ; i < number; i++) { cin>>x>>y; point tmp(x, y); p.insert(tmp); } } //计算凸包函数 void Caculate (deque<point> &deq , multiset<point> st , int number) { //插入三个 multiset<point>::iterator i = st.begin(); for (int j = 0; j < 3; j++ , i++) deq.push_back(*i); for (; i != st.end(); i++) { //核心代码：判断合适的点 while (1) { point top = deq.back(); point second = *(deq.end() - 2); point Pi = *i; //计算向量P1（x1 ， y1） double x1 = second.x - top.x; double y1 = second.y - top.y; //计算向量P2（x2 ， y2） double x2 = Pi.x - top.x; double y2 = Pi.y - top.y; //p1 和 p2 的叉乘< 0 if ((x1 * y2 - x2 * y1) < 0) deq.pop_back(); else break; } deq.push_back(*i); } } int main () { int point_number; //点的个数 multiset<point> st; //初始时的点保存在st里头 deque<point> deq; //保存最后的点 Input(st,point_number); //输入点并且排序 Caculate (deq , st , point_number); //计算凸包 //输出组成凸包的点 for (int i = 0; i < deq.size(); i++) { cout<<deq[i].x<<ends<<deq[i].y<<ends<<deq[i].angle<<endl; } return 0; } //测试数据 /* 10 2 2 3 3 4 4 4 1 5 2 5 3 6 5 5 6 3 6 1 4 */