• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            先求凸包,答案是凸包周長(zhǎng) + 2πl。因?yàn)楹?jiǎn)單多邊形的轉(zhuǎn)角是360度,所以加上一個(gè)圓的周長(zhǎng)。

             

            /*************************************************************************
            Author: WHU_GCC
            Created Time: 2007-8-21 15:02:36
            File Name: pku1113.cpp
            Description: 
            ***********************************************************************
            */

            #include 
            <iostream>
            #include 
            <cmath>
            using namespace std;

            #define out(x) (cout << #x << ": " << x << endl)
            typedef 
            long long int64;
            const int maxint = 0x7FFFFFFF;
            const int64 maxint64 = 0x7FFFFFFFFFFFFFFFLL;
            template 
            <class T> void show(T a, int n) for (int i = 0; i < n; ++i) cout << a[i] << ' '; cout << endl; }
            template 
            <class T> void show(T a, int r, int l) for (int i = 0; i < r; ++i) show(a[i], l); cout << endl; }

            const int maxn = 1001;

            typedef 
            struct point_t
            {
                
            int x, y;
            }
            ;

            typedef 
            struct polygon_t
            {
                
            int n;
                point_t p[maxn];
            }
            ;

            int operator <(const point_t &a, const point_t &b)
            {
                
            return a.y < b.y || a.y == b.y && a.x < b.x;
            }


            point_t 
            operator -(const point_t &a, const point_t &b)
            {
                point_t ret;
                ret.x 
            = a.x - b.x;
                ret.y 
            = a.y - b.y;
                
            return ret;
            }


            double dist(const point_t &a, const point_t &b)
            {
                
            return sqrt(double((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y)));
            }


            int cross(const point_t &a, const point_t &b)
            {
                
            return a.x * b.y - a.y * b.x;
            }


            int turn_left(const point_t &a, const point_t &b, const point_t &c)
            {
                
            return cross(b - a, c - b) > 0;
            }


            class point_set_c
            {
            public:
                
            void init(int _n, point_t _p[]);
                polygon_t convex_hull();
            private:
                
            int n;
                point_t p[maxn];
            }
            ;

            void point_set_c::init(int _n, point_t _p[maxn])
            {
                n 
            = _n;
                
            for (int i = 0; i < n; i++)
                    p[i] 
            = _p[i];
            }


            polygon_t point_set_c::convex_hull()
            {
                
            int stack[maxn];
                
            int top = 1;
                stack[
            0= 0;

                sort(p, p 
            + n);

                
            for (int i = 1; i < n;)
                
            {
                    
            if (top == 1 || turn_left(p[stack[top - 2]], p[stack[top - 1]], p[i]))
                        stack[top
            ++= i++;
                    
            else top--;
                }

                
            int t_top = top;
                
            for (int i = n - 2; i >= 0;)
                
            {
                    
            if (top == t_top || turn_left(p[stack[top - 2]], p[stack[top - 1]], p[i]))
                        stack[top
            ++= i--;
                    
            else top--;
                }

                
                polygon_t ret;
                ret.n 
            = 0;
                
            for (int i = 0; i < top - 1; i++)
                    ret.p[ret.n
            ++= p[stack[i]];
                
            return ret;
            }


            point_t p[maxn];
            int n;
            point_set_c ps;
            polygon_t poly;

            int main()
            {
                
            int l;
                
            while (scanf("%d%d"&n, &l) != EOF)
                
            {
                    
            for (int i = 0; i < n; i++)
                        scanf(
            "%d%d"&p[i].x, &p[i].y);
                    ps.init(n, p);
                    poly 
            = ps.convex_hull();

                    
            double ans = 0.0;
                    
            for (int i = 0; i < poly.n; i++)
                        ans 
            += dist(poly.p[i], poly.p[(i + 1% poly.n]);
                    ans 
            += 2 * l * acos(-1.0);
                    printf(
            "%.0lf\n", ans);
                }

                
            return 0;
            }
            posted on 2007-08-21 15:43 Felicia 閱讀(495) 評(píng)論(0)  編輯 收藏 引用 所屬分類(lèi): 計(jì)算幾何
             
            久久99毛片免费观看不卡| 久久天天躁狠狠躁夜夜2020一| 无码人妻少妇久久中文字幕蜜桃| 色诱久久av| 无码国产69精品久久久久网站| 久久精品国产亚洲av影院| 国产亚州精品女人久久久久久 | 国产精品九九久久精品女同亚洲欧美日韩综合区 | 国产精品对白刺激久久久| www亚洲欲色成人久久精品| 国产精品久久久久久久人人看| 久久不射电影网| 少妇熟女久久综合网色欲| 亚洲欧美精品伊人久久| 午夜天堂av天堂久久久| 久久中文精品无码中文字幕| 99久久精品国内| 精品久久久久久国产潘金莲| 久久久久99这里有精品10| 国内精品久久久久国产盗摄| 精品国际久久久久999波多野| 无码8090精品久久一区| 久久久综合九色合综国产| 婷婷五月深深久久精品| 久久人人添人人爽添人人片牛牛| 久久久人妻精品无码一区| 亚洲国产成人久久精品影视| 国产精品18久久久久久vr| 久久久久久久亚洲Av无码| 国产精品久久国产精品99盘 | 久久99精品国产一区二区三区| 午夜精品久久久久久久无码| 久久AⅤ人妻少妇嫩草影院| 精品久久777| 色综合久久综精品| 大蕉久久伊人中文字幕| 91久久精品无码一区二区毛片| 97超级碰碰碰久久久久| 久久久精品午夜免费不卡| 97久久精品人人澡人人爽| 久久免费国产精品|