COGS896圈奶牛

1. 圈奶牛
   ★★☆ 输入文件：fc.in 输出文件：fc.out 简单对比
   时间限制：1 s 内存限制：128 MB
   USACO/fc(译by Felicia Crazy)
   描述
   农夫约翰想要建造一个围栏用来围住他的奶牛，可是他资金匮乏。他建造的围栏必须包括他的奶牛喜欢吃草的所有地点。对于给出的这些地点的坐标，计算最短的能够围住这些点的围栏的长度。
   PROGRAM NAME: fc
   INPUT FORMAT(file fc.in)
   输入数据的第一行包括一个整数 N。N（0 <= N <= 10,000）表示农夫约翰想要围住的放牧点的数目。接下来 N 行，每行由两个实数组成，Xi 和 Yi,对应平面上的放牧点坐标（-1,000,000 <= Xi,Yi <= 1,000,000）。数字用小数表示。
   OUTPUT FORMAT(file fc.out)
   输出必须包括一个实数，表示必须的围栏的长度。答案保留两位小数。
   SAMPLE INPUT (file fc.in)
   4
   4 8
   4 12
   5 9.3
   7 8
   SAMPLE OUTPUT (file fc.out)
   12.00
   裸凸包。。
   要进攻计算几何了233
   是一个值得庆祝的日子。。
   附上本蒟蒻的代码：

#include<cstdio>#include<iostream>#include<cmath>#include<algorithm>using namespace std;int n,top;double ans;struct node{    double x,y;};node p[10001],s[10001];double dis(const node a,const node b){    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}double mul(const node a,const node b,const node p0){    return (a.x-p0.x)*(b.y-p0.y)-(b.x-p0.x)*(a.y-p0.y);}bool cmp(const node a,const node b){    if (mul(a,b,p[0])==0)      return dis(a,p[0])<dis(b,p[0]);    return mul(a,b,p[0])>0;}void graham(){    int k=0,i;    node t;    top=2;    for (i=1;i<n;i++)      if ((p[k].y>p[i].y) || (p[k].y==p[i].y && p[k].x>p[i].x))        k=i;    t=p[0];    p[0]=p[k];    p[k]=t;    sort(p+1,p+n,cmp);    s[0]=p[0];    s[1]=p[1];    s[2]=p[2];    for (i=3;i<n;i++)      {        while (top && mul(p[i],s[top],s[top-1])>=0)          top--;        top++;        s[top]=p[i];      }    top++;    s[top]=p[0];    for (i=0;i<top;i++)      ans+=dis(s[i],s[i+1]);}int main(){    freopen("fc.in","r",stdin);    freopen("fc.out","w",stdout);    int i;    scanf("%d",&n);    for (i=0;i<n;i++)      scanf("%lf%lf",&p[i].x,&p[i].y);    graham();    printf("%.2lf",ans);    fclose(stdin);    fclose(stdout);    return 0;}