poj 1113 凸包问题
我的第一道凸包题，终于弄懂了Graham算法，对我来说弄懂一个算法是一个痛苦的过程，我痛苦了我好几天……
Grahma算法模板：
1、 寻找一个基准点，这个点通常是最左边的最下面那个点（其实只要这个点一定会在凸包边界上就行），然后将之个点置于第一个位置。
2、 将所有点按相对于基准点的极角排序。
3、 因为排序之后前两个点一定会在凸包上，假设第三个点也在凸包上，因此我们可以从第四个点开始依次扫描每个点。对于每个当前点他在当前时刻是极角最大的点，所以在当前时刻它必定在凸包上，那么此时就可能改变了它前一个点是不是在凸包上。怎么判断呢？我们可以咱当前点到推三个点，假设是p1,p2,p3，判断这三个点所形成的角是不是大于180度，小于180度继续扫描下一个点，若大于180度说明第二个点一定在凸包内部，然后我们暂时停止扫描下一个点，将p1,p2各自后退一个点继续判断p1,p2,p3所组成的角是不是大于180度，（这也就是所谓的“退化”问题，看了《算法导论》后才知道这儿需要有个循环）若仍然大于180度循环这种操作直到三个点组成的角小于等于180度，最后继续扫描下一个点，直到结束。在这个过程中我们可以用栈来保存凸包边上的点。
该题需要最后算出各边长的和再加上圆的周长。
#include<iostream>
#include<cmath>
using namespace std;
#define MAX_INT 123456789
struct point
{
int x,y;
};
point vertex[1001],stack[1000];
int angle(point p1,point p2,point p0)//判断极角大小
{
return (p1.y-p0.y)*(p2.x-p0.x)-(p2.y-p0.y)*(p1.x-p0.x);
}
int dist(point p1,point p2)
{
return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
int cmp(const void* pp,const void* qq)//排序判断条件
{
point p1,p2;
p1=*(point*) pp; p2=*(point*) qq;
if(angle(p1,p2,vertex[0])>0)
return 1;
else if(angle(p1,p2,vertex[0])==0 && dist(p1,vertex[0])>dist(p2,vertex[0]))
return 1;
return -1;
}
int judge(point p1,point p2,point p3)//判断三点组成的角度
{
return (p3.y-p2.y)*(p2.x-p1.x)-(p2.y-p1.y)*(p3.x-p2.x);
}
double Graham(int n,int m)
{
int sum_max=0,i=0;
point p1,p2;
double sum=0;
p1=vertex[1],p2=vertex[2];
stack[0]=vertex[0];
stack[1]=vertex[1];
stack[2]=vertex[2];
sum_max=2;
for(i=3;i<n;i++)
{
while(1)//解决退化问题
{
if(judge(p1,p2,vertex[i])>=0)
{
p1=p2,p2=vertex[i];
stack[++sum_max]=vertex[i];
break;
}
else
{
--sum_max;
p2=stack[sum_max];
p1=stack[sum_max-1];
}
}
}
for(i=0;i<sum_max;i++)
sum+=sqrt((double)dist(stack[i],stack[i+1]));
sum+=sqrt((double)dist(stack[sum_max],vertex[0]))+3.1415926*2*m;
return sum;
}
int main()
{
int i,m,n,min;
point minPoint;
while(scanf("%d%d",&n,&m)!=EOF)
{
memset(vertex,0,sizeof(vertex));
minPoint.x=minPoint.y=MAX_INT;
for(min=i=0;i<n;i++)
{
scanf("%d%d",&vertex[i].x,&vertex[i].y);
if(minPoint.x>vertex[i].x)
{
minPoint=vertex[i];
min=i;
}
else if(minPoint.x==vertex[i].x && minPoint.y>vertex[i].y)
{
minPoint=vertex[i];
min=i;
}
}
vertex[min]=vertex[0]; vertex[0]=minPoint;
qsort(vertex+1,n-1,sizeof(vertex[0]),cmp);
printf("%.lf\n",Graham(n,m));
}
return 0;
}