POJ1113——WALL 凸包

调了两个小时（吐

poj1113

凸包有毒啊？？？

我闲的蛋疼用了一个函数叫atan2

然而这个东西用法有点神奇。。是用来求极角的

atan2(y,x） 注意是y,x 得到的是(x,y)点对对于0,0的弧度！注意是弧度！

弧度转化角度：

#define atan2(x,y) atan2(x,y)*180/3.1415926535

当然其实这题。。不转化角度也没关系，显然弧度与角度是增函数的关系没毛病啊

凸包过程：

1、找到点中左下角（先下，再左）的点，以它为原点

2、求出所有点的极角

3、按照极角排序

4、维护一个栈，保持凸多边形性质，即新加入的点不能在原先的边的右侧

5、完了。

这里我判断是否在右侧是用叉积（因为刚学）

其实也可以直接用解析式判断，具体不细讲

接下来说一下这个题的几个坑点….

1、角度要开double=_=

好吧这个是我傻逼了

2、共线情况要先按纵坐标从小到大排序，再按横坐标从大到小排序

画图理解一下就行，横坐标。。实际上是一样的。。。。

CODE：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<string>
#include<map>
#include<cstring>
#include<vector>
#define inf 1e9
#define ll long long
#define atan2(x,y) atan2(x,y)*180/3.1415926535
#define For(i,j,k) for(ll i=j;i<=k;i++)
#define Dow(i,j,k) for(ll i=k;i>=j;i--)
using namespace std;
struct point
{
int  x,y;
double aph;
}a[2001];
int n;
int d;
int sta[5001];
inline bool cmp(point t1,point t2)
{
if(t1.aph!=t2.aph)return t1.aph<t2.aph;
if(t1.y!=t2.y) t1.y<t2.y;
return t1.x<t2.x;
}
inline double sqr(double x)
{
return x*x;
}
int main()
{
scanf("%d%d",&n,&d);
int t=0;
a[0].y=1e9;
For(i,1,n)
{
scanf("%d%d",&a[i].x,&a[i].y);
if(a[i].y<a[t].y||(a[i].y==a[t].y&&a[i].x<a[t].x))
t=i;
}
swap(a[t],a[1]);
For(i,2,n)  a[i].x-=a[1].x,a[i].y-=a[1].y;
a[1].x=a[1].y=0;
For(i,1,n)
a[i].aph=atan2(a[i].y,a[i].x);
sort(a+2,a+n+1,cmp);
sta[1]=1;sta[2]=2;
int top=2;
For(i,3,n)
{
while(1)
{
int p1=sta[top],p2=sta[top-1];
int S=(a[p1].x-a[p2].x)*(a[i].y-a[p2].y)-(a[i].x-a[p2].x)*(a[p1].y-a[p2].y);
if(S<0)
top--;else break;
}
sta[++top]=i;
}
double  ans=0;
For(i,1,top)
ans+=sqrt(sqr((double)a[sta[i%top+1]].x-(double)a[sta[i]].x)+sqr((double)a[sta[i%top+1]].y-(double)a[sta[i]].y));

ans+=3.1415926*(double)d*2.00;
printf("%.0f\n",ans);
}