Graham算法的思路,大概如下:对平面上的点的集合,从中找到有最小的y坐标值的点p,然后根据其它点和p的连线与正x轴所成的角度将平面上的点进行排序,排序后,扫描从p开始的有序列表,如果所有的这些点都在凸包上,那么每三个相继的点,会组成一个左旋,从另一方面说,如果相继的三个点,p1,p2,p3,组成了一个右旋,则可以立即去除p2,因为它不可能在凸包上!如此扫描到p3等于p时,扫描就结束,剩下的点,就全都是凸包上的点了!
这里面有几个难点:
1、如何根据他们的夹角,来对双链表进行排序
2、如何判断左右旋
3、如何扫描
针对第一个问题,在下面的程序中,采用的是简单的选择排序,即每次在无序区内找到最大的元素,将其放到有序区内,因为这是对双链表进行排序,就有双链表的排序特点,但我认为采用插入,即断开结点之间的联系始终不是好的办法,所以采用了这种简单的交换的方法,虽然效率慢了点,但是重在说明问题!
针对第二个问题,判断左右旋,其实可以用一个行列式的值来判定
|x1 x2 x3|
|y2 y2 y3|
|1 1 1 |
若该行列式为正,则为左旋,反之,则为右旋。其实这个行列式是根据点是在直线的上方还是下方而推断出来的,用直线方程也可以判断左右旋。
针对第三个问题看下面的代码中的scan()函数就很清楚了!
下面是用c++写的具体代码:
#include <iostream>
#include <cmath>
using namespace std;
struct point
{
float x,y;//坐标
float cos;//该点和x轴的正向所形成的角的余弦
struct point *prev,*next;
};
typedef struct point Type;
/*创建一个双向链表*/
void creatLinklist_Du(Type *L,int n)
{
Type *p;
Type *s;
p=L;
cin>>p->x>>p->y;
p->prev=NULL;
p->next=NULL;//初始化第一个结点
int i;
for(i=0;i<n-1;i++)
{
s=new Type;//这里忘记给S分配空间了,所以出错了。
cin>>s->x>>s->y;
s->next=NULL;
s->prev=p;
p->next=s;
p=s;
}
}
void calculateCos(Type *L)//找到y坐标最小的元素,让它是第一个点,并计算每个点的cos值
{
Type *p=L,*pMin=L;
float min=L->y;
while(p)
{//找到y坐标最小的结点
if(min>p->y)
{
pMin=p;
min=p->y;
}
p=p->next;
}
float x0=pMin->x;
float y0=pMin->y;
p=L;
while(p)
{
float x=p->x;
float y=p->y;
if(x==x0 && y==y0)
p->cos=1;
else
p->cos=(x-x0)/sqrt((x-x0)*(x-x0)+(y-y0)*(y-y0));
p=p->next;
}
}
/*简单选择排序*/
void simpleSlectSort(Type *L)//按cos值进行降序排序
{
Type *p,*pMax,*q,*temp=new Type;
for(p=L;p->next;p=p->next)
{
pMax=p;
for(q=p->next;q;q=q->next)//找到无序区的最大的cos
{
if(q->cos > pMax->cos)
pMax=q;
}
if(p!=pMax)//将此最大值放进有序区,这里只需要交换两个结点的内容就可以了
{
temp->cos=pMax->cos;
temp->x=pMax->x;
temp->y=pMax->y;
pMax->cos=p->cos;
pMax->x=p->x;
pMax->y=p->y;
p->cos=temp->cos;
p->x=temp->x;
p->y=temp->y;
}
}
}
float area(float x1,float y1,float x2,float y2,float x,float y)//计算三个点的左旋还是右旋
{
return (x1*y2+x2*y+x*y1)-(x*y2+x1*y+x2*y1);
}
void scan(Type *L)//顺序扫描双链表,判断相继的三个点是否组成左旋或是右旋
{
Type *p=L,*p1=L,*p2,*p3;
float temp;
do{
p2=p1->next;
if(p2->next)
p3=p2->next;
else
p3=p;
temp=area(p1->x,p1->y,p2->x,p2->y,p3->x,p3->y);
if(temp>=0.0)//如果temp为正,则为左旋
p1=p1->next;
else
{
p1->next=p3;
p3->prev=p1;
delete p2;
p1=p1->prev;
}
}while(!(p3==p && temp>=0.0));
}
/*输出双向链表中的元素*/
void outputLinklist_Du(Type *L)
{
Type *p;
p=L;
while(p)
{
cout<<p->x<<" "<<p->y<<"/t"<<p->cos<<endl;
p=p->next;
}
}
int main()
{
//freopen("in.txt","r",stdin);
Type *L=new Type;
creatLinklist_Du(L,4);
calculateCos(L);
simpleSlectSort(L);
scan(L);
outputLinklist_Du(L);
return 0;
}