首页 > 代码库 > POJ 1113 Wall(Graham求凸包周长)
POJ 1113 Wall(Graham求凸包周长)
题目链接
题意 : 求凸包周长+一个完整的圆周长。
因为走一圈,经过拐点时,所形成的扇形的内角和是360度,故一个完整的圆。
思路 : 求出凸包来,然后加上圆的周长
1 #include <stdio.h> 2 #include <string.h> 3 #include <iostream> 4 #include <cmath> 5 #include <algorithm> 6 7 const double PI = acos(-1.0) ; 8 using namespace std ; 9 10 struct point11 {12 double x,y ;13 }p[110],p1[110];14 int n ,m;15 16 double dis(point a,point b)17 {18 return sqrt(double((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y))) ;19 }20 double cross(point a,point b,point c)21 {22 return (a.x-c.x)*(b.y-c.y) - (a.y-c.y)*(b.x-c.x) ;23 }24 bool cmp(point a,point b)25 {26 if(cross(a,b,p[0]) > 0 || cross(a,b,p[0]) == 0 && dis(a,p[0]) < dis(b,p[0]))27 return true ;28 return false ;29 }30 void Graham()31 {32 int cnt;33 sort(p+1,p+n,cmp) ;34 p[n] = p[0] ;35 p1[0] = p[0] ;36 p1[1] = p[1] ;37 cnt = 1 ;38 for(int i = 2 ; i <= n ; i++)39 {40 while(cnt >= 1 && cross(p1[cnt-1],p1[cnt],p[i]) <= 0) cnt -- ;41 p1[ ++cnt] = p[i] ;42 }43 double sum = 0.0 ;44 for(int i = 0 ; i < cnt ; i++)45 {46 sum += dis(p1[i],p1[i+1]) ;47 }48 printf("%.0lf\n",sum+2*PI*m) ;49 }50 int main()51 {52 while(~scanf("%d %d",&n,&m))53 {54 int pos = 0 ;55 for(int i = 0 ; i < n ; i++)56 {57 scanf("%lf %lf",&p[i].x,&p[i].y) ;58 if(p[pos].y > p[i].y || (p[pos].y == p[i].y && p[i].x < p[pos].x))59 pos = i ;60 }61 // if(n == 1) {puts("0.00");continue ;}62 // if(n == 2)63 // {64 // printf("%.2lf\n",dis(p[1],p[0])) ;65 // continue ;66 // }67 swap(p[pos],p[0]) ;68 //printf("%d %d\n",p[0].x,p[0].y) ;69 Graham() ;70 }71 return 0 ;72 }
POJ 1113 Wall(Graham求凸包周长)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。