首页 > 代码库 > 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 }
View Code

 

POJ 1113 Wall(Graham求凸包周长)