poj1113 第一道凸包，一次a

#include <iostream>
#include<math.h>
#include<stdio.h>
#include <algorithm>
using namespace std;
const double pi=acos(-1.0);
const double eps=1E-9;
typedef struct
{
  double x,y;
} Point;
Point point[1100];
Point ch[1100];
int length;
int N,L;
double multiply(Point p1,Point p2,Point p0)
{
  return((p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y));
}
double dis(Point p1,Point p2)
{
  return(sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)));
}
void Graham(Point PointSet[],Point ch[],int n,int &length)
{
  int i,j,k=0,top=2;
  Point tmp;
  //找到最下且偏左的那个点
  for(i=1; i<n; i++)
    if ((PointSet[i].y<PointSet[k].y)||((PointSet[i].y==PointSet[k].y)&&(PointSet[i].x<PointSet[k].x)))
      k=i;
  //将这个点指定为PointSet[0]
  tmp=PointSet[0];
  PointSet[0]=PointSet[k];
  PointSet[k]=tmp;
  //按极角从小到大,距离偏短进行排序
  for (i=1; i<n-1; i++)
    {
      k=i;
      for (j=i+1; j<n; j++)
        if((multiply(PointSet[j],PointSet[k],PointSet[0])>0)||((multiply(PointSet[j],PointSet[k],PointSet[0])==0)
            &&(dis(PointSet[0],PointSet[j])<dis(PointSet[0],PointSet[k]))))
          k=j; //k保存极角最小的那个点,或者相同距离原点最近
      if(i!=k)
        {
          tmp=PointSet[i];
          PointSet[i]=PointSet[k];
          PointSet[k]=tmp;
        }
    }


  //前三个点先入栈
  ch[0]=PointSet[0];
  ch[1]=PointSet[1];
  ch[2]=PointSet[2];
  //判断与其余所有点的关系
  for (i=3; i<n; i++)
    {
      //不满足向左转的关系,栈顶元素出栈
      while(multiply(PointSet[i],ch[top],ch[top-1])>=0) top--;
      //当前点与栈内所有点满足向左关系,因此入栈.
      ch[++top]=PointSet[i];
    }
  length=top+1;
}






int main()
{


    double sum=0;
    scanf("%d%d",&N,&L);
    for (int i=0;i<N;i++ )
        scanf("%lf%lf",&point[i].x,&point[i].y);
    Graham(point,ch,N,length);
        for (int i=0;i<length-1;i++ )
        sum+=dis(ch[i],ch[i+1]);
        sum+=dis(ch[length-1],ch[0]);
        printf("%.0f",sum+2*pi*L);
return 0;


}