★★☆ 输入文件:fc.in
输出文件:fc.out
简单对比
时间限制:1 s 内存限制:128 MB
描述
农夫约翰想要建造一个围栏用来围住他的奶牛,可是他资金匮乏。他建造的围栏必须包括他的奶牛喜欢吃草的所有地点。对于给出的这些地点的坐标,计算最短的能够围住这些点的围栏的长度。
输入数据的第一行包括一个整数 N。N(0 <= N <= 10,000)表示农夫约翰想要围住的放牧点的数目。接下来 N 行,每行由两个实数组成,Xi 和 Yi,对应平面上的放牧点坐标(-1,000,000 <= Xi,Yi <= 1,000,000)。数字用小数表示。
输出必须包括一个实数,表示必须的围栏的长度。答案保留两位小数。
4
4 8
4 12
5 9.3
7 8
12.00
题解:裸凸包求长度,上模板。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstdlib> 4 #include<cstring> 5 #include<algorithm> 6 #include<cmath> 7 #include<queue> 8 #include<vector> 9 using namespace std; 10 const int eps=1e-7; 11 int N,top; 12 double ans; 13 struct P{ 14 double x,y; 15 friend P operator-(P a,P b){ 16 P t; t.x=a.x-b.x; t.y=a.y-b.y; 17 return t; 18 } 19 friend double operator*(P a,P b){ 20 return ((a.x*b.y)-(b.x*a.y)); 21 } 22 friend double dis(P a,P b){ 23 return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); 24 } 25 }p[20000],stac[20000]; 26 inline bool operator<(P a,P b){ 27 double t=(a-p[1])*(b-p[1]); 28 if(fabs(t)<=eps) return dis(a,p[1])<dis(b,p[1]); 29 return t>eps; 30 } 31 inline void graham(){ 32 int tmp=1; 33 for(int i=2;i<=N;i++){ 34 if(p[i].y<p[tmp].y||(fabs(p[i].y-p[tmp].y)<=eps&&p[i].x<p[tmp].x)) tmp=i; 35 } 36 swap(p[1],p[tmp]); 37 sort(p+2,p+N+1); 38 stac[++top]=p[1]; stac[++top]=p[2]; 39 for(int i=3;i<=N;i++){ 40 while((stac[top]-stac[top-1])*(p[i]-stac[top-1])<=eps) top--; 41 stac[++top]=p[i]; 42 } 43 stac[top+1]=p[1]; 44 for(int i=1;i<=top;i++){ 45 ans+=dis(stac[i],stac[i+1]); 46 } 47 printf("%.2lf",ans); 48 } 49 int main(){ 50 freopen("fc.in","r",stdin); 51 freopen("fc.out","w",stdout); 52 scanf("%d",&N); 53 for(int i=1;i<=N;i++) scanf("%lf%lf",&p[i].x,&p[i].y); 54 graham(); 55 return 0; 56 }
原文:http://www.cnblogs.com/CXCXCXC/p/5244262.html