Section 5.1 Fencing the Cows

/*
ID: niepeng1
PROG: fc
LANG: C++
*/
/*
对那个text的讲解没有什么印象，就又参考了别人的程序。
典型的凸包题目，就是选定一个定点（左下），然后从这个定点对所有的定点排序，
逆时针添加定点，每次添加后如果最后三个点成的角度大于180°，就循环取出最后
那个定点。
最后凸包已经构成，直接计算就行了。
此时方法要求没有任何限制，如果有定点个数或凸包形状，就......
*/
#include <iostream>
#include <math.h>
#include <algorithm>
using namespace std;
const int NODENUM=10001;
const double PZERO=1e-6;
struct Node{
 double x;
 double y;
 double angle;
 double dis;
}ver[NODENUM],tem;
int stack[NODENUM+1],n;
double ans;
bool cmp(Node a,Node b)
{
 return a.angle<b.angle;
}
double dis(int s,int t)
{
 return sqrt((ver[s].x-ver[t].x)*(ver[s].x-ver[t].x)+(ver[s].y-ver[t].y)*(ver[s].y-ver[t].y));
}
void Calc(int s)
{
 ver[s].dis=dis(s,0);
 ver[s].angle=acos((ver[s].x-ver[0].x)/ver[s].dis);
}
int cross(int o,int f,int t)
{//判断夹角大小，以180°为参考
 double ret;
 ret=(ver[f].x-ver[o].x)*(ver[t].y-ver[o].y)-(ver[t].x-ver[o].x)*(ver[f].y-ver[o].y);
 if (ret>PZERO) return 1;
 else if (ret<-PZERO) return -1;
 else return 0;
}
void CutVer()
{
 int now,i;
 for (now=1,i=1;i<n;i++)
 {//从出发点引出线，点在同一条直线上的取距离最远的点
  while (i<n&&fabs(ver[i].angle-ver[now].angle)<=PZERO)
  {
   if (ver[i].dis>ver[now].dis) ver[now]=ver[i];
   i++;
  }
  if (i==n) break;
  now++;ver[now]=ver[i];
 }
 n=now+1;
}
void Graham()
{//每次添加一个点后，从最后一个点向前扫描，去除那些使角度大于180度的角
 int i;
 stack[0]=3;stack[1]=0;stack[2]=1;stack[3]=2;
 for (i=3;i<n;i++)
 {
  while (cross(stack[stack[0]-1],stack[stack[0]],i)<0) stack[0]--;
  stack[++stack[0]]=i;
 }
}

int main()
{
 int i,j;
 freopen("fc.in","r",stdin);
 freopen("fc.out","w",stdout);
 scanf("%d",&n);
 j=0;
 for (i=0;i<n;i++)
 {
  scanf("%lf%lf",&ver[i].x,&ver[i].y);
  if (ver[i].y<ver[j].y||(ver[i].y==ver[j].y&&ver[i].x<ver[j].x))
   j=i;
 }//取出最左下角的定点，作为出发点
 tem=ver[0];ver[0]=ver[j];ver[j]=tem;
 for (i=1;i<n;i++) Calc(i);
 ver[0].angle=0;
 sort(ver+1,ver+n,cmp);//按照角度排序
 CutVer();
 Graham();
 ans=dis(stack[stack[0]],stack[1]);//计算凸包的周长
 for (i=1;i<stack[0];i++)
  ans+=dis(stack[i],stack[i+1]);
 printf("%.2lf/n",ans);
 return 0;

}