——————byMBA智库 http://wiki.mbalib.com/wiki/%E5%87%B8%E5%8C%85
算法
这里就介绍一个Graham扫描算法
Graham扫描法
由最底的一点A1开始,计算它跟其他各点的连线和x轴的角度,按小至大将这些角度排序,称它们的对应点为A_2,A_3,\ldots,A_n。这里的时间复杂度可达O(nlogn)。
考虑最小的角度对应的点A3。若由A2到A3的路径相对A1到A2的路径是向右转的(可以想象一个人沿A1走到A2,他站在A2时,是向哪边改变方向),表示A3不可能是凸包上的一点,考虑下一点由A2到A4的路径;否则就考虑A3到A4的路径是否向右转……直到回到A1。
这个算法的整体时间复杂度是O(nlogn),注意每点只会被考虑一次,而不像Jarvis步进法中会考虑多次。
这个算法由葛立恒在1972年发明。它的缺点是不能推广到二维以上的情况。
代码如下:
#include<cstdio>
#include<stack>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn = 1000;
struct Node{
int x,y;
}u,node[maxn];
int N;
int cmp1(const Node &a, const Node &b)
{
return a.y < b.y;
}
/*极角排序*/
int cmp2(const Node &a, const Node &b)
{
return atan2(1.0*(a.y-u.y), 1.0*(a.x-u.x)) < atan2(1.0*(b.y-u.y), 1.0*(b.x-u.x));
}
/*计算叉积*/
int Direction(Node a, Node b, Node c)
{
return (c.x-a.x) * (b.y-a.y) - (b.x-a.x) * (c.y - a.y);
}
/*判断是否向右*/
bool right(Node a, Node b, Node c)
{
if(Direction(a, b, c) < 0)
{
return true;
}
else
return false;
}
int main()
{
scanf("%d",&N);
for(int i=0; i<N; i++)
{
scanf("%d%d",&node[i].x, &node[i].y);
}
sort(node, node+N, cmp1);
stack<Node> s;
s.push(node[0]);
u = node[0];
sort(node+1, node+N, cmp2);
s.push(node[1]);
s.push(node[2]);
for(int i=3; i<N; i++)
{
Node top = s.top();
s.pop();
Node n_top = s.top();
if(right(node[i], top, n_top))
{
s.push(node[i]);
}
else{
s.push(top);
s.push(node[i]);
}
}
while(!s.empty())
{
printf("%d %d\n",s.top().x, s.top().y);
s.pop();
}
return 0;
}