• <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>

            對多邊形的每個邊計算和原點的夾角范圍,求并得到總夾角范圍 alpha
            答案就是 min(alpha, 2 * pi) * k * h

            /*************************************************************************
            Author: WHU_GCC
            Created Time: 2007-8-11 20:01:23
            File Name: pku1031.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 double eps = 1e-10;
            const double pi = acos(-1.0);
            const int maxn = 110;

            typedef 
            struct point_t
            {
                
            double x, y;
            }
            ;

            typedef 
            struct interval_t
            {
                
            double l, r;
            }
            ;

            double get_angle(double x, double y)
            {
                
            if (abs(x) < eps)
                
            {
                    
            if (y > 0return pi / 2.0;
                    
            else return 3.0 * pi / 2.0;
                }

                
            else
                
            {
                    
            double t = atan(y / x);
                    
            if (x < 0)
                        t 
            += pi;
                    
            if (t < 0) t += 2.0 * pi;
                    
            return t;
                }

            }


            interval_t get_range(point_t a, point_t b)
            {
                
            if (a.x * b.y - a.y * b.x < eps)
                
            {
                    swap(a, b);
                }

                
            double t1 = get_angle(a.x, a.y);
                
            double t2 = get_angle(b.x, b.y);
                
            if (t2 < t1 - eps) t2 += 2.0 * pi;
                interval_t ret;
                ret.l 
            = t1;
                ret.r 
            = t2;
                
            return ret;
            }


            bool operator <(const interval_t &a, const interval_t &b)
            {
                
            return a.l < b.l || a.l == b.l && a.r < b.r;
            }


            bool between(const double &a, const double &x, const double &b)
            {
                
            return x > a - eps && x < b + eps;
            }


            point_t p[maxn];
            interval_t interval[maxn];

            int main()
            {
                
            double k, h;
                
            int n;
                scanf(
            "%lf%lf%d"&k, &h, &n);
                
            {
                    
            for (int i = 0; i < n; i++)
                        scanf(
            "%lf%lf"&p[i].x, &p[i].y);
                    
            for (int i = 0; i < n; i++)
                        interval[i] 
            = get_range(p[i], p[(i + 1% n]);
                    
                    
            double l = interval[0].l, r = interval[0].r;
                    
            for (int i = 1; i < n; i++)
                    
            {
                        
            for (int j = -1; j <= 1; j++)
                        
            {
                            
            if (between(l, interval[i].l + j * 2 * pi, r) || between(l, interval[i].r + j * 2 * pi, r))
                            
            {
                                l 
            = min(l, interval[i].l + j * 2 * pi);
                                r 
            = max(interval[i].r + j * 2 * pi, r);
                                
            break;
                            }

                        }

                    }

                    printf(
            "%.2lf\n", min(r - l, 2.0 * pi) * k * h);
                }

                
            return 0;
            }
            posted on 2007-08-12 12:06 Felicia 閱讀(900) 評論(2)  編輯 收藏 引用 所屬分類: 計算幾何
            Comments
            • # re: [計算幾何]pku1031
              panjinww
              Posted @ 2007-09-23 15:03
              很漂亮的代碼,可是少了注釋!!  回復  更多評論   
            • # re: [計算幾何]pku1031
              Felicia
              Posted @ 2007-09-23 15:49
              @panjinww
              呵呵,算法已經很清楚了吧,就是實現的問題了。  回復  更多評論   
             
            99久久精品毛片免费播放| 97久久综合精品久久久综合| 精品一区二区久久| 精品久久久久久无码中文野结衣| 人人狠狠综合久久亚洲| 久久夜色精品国产噜噜亚洲AV| 91精品久久久久久无码| 精品国产乱码久久久久久呢| 一本久久a久久精品综合夜夜| 久久久久亚洲AV片无码下载蜜桃 | 久久久久久毛片免费看| 亚洲色大成网站WWW久久九九| 亚洲午夜久久久精品影院| 久久一日本道色综合久久| 性做久久久久久久久浪潮| 婷婷综合久久中文字幕| 亚洲伊人久久成综合人影院| 久久亚洲AV成人出白浆无码国产| 欧美久久天天综合香蕉伊| 99国产精品久久久久久久成人热| 色老头网站久久网| 久久人人爽人爽人人爽av| 狠狠人妻久久久久久综合蜜桃 | 久久久久无码精品国产| 久久国产欧美日韩精品免费| 久久久91人妻无码精品蜜桃HD| 欧美日韩中文字幕久久伊人| 国产精品99久久免费观看| 日产精品久久久一区二区| 亚洲国产精品无码久久一线| 久久精品国产亚洲av麻豆图片 | 精品一区二区久久| 欧美一区二区精品久久| 国产精品免费久久久久电影网| 国产精品久久久久无码av| 精品久久777| 伊人久久精品线影院| 亚洲国产成人乱码精品女人久久久不卡 | 久久精品免费观看| 久久久久99精品成人片三人毛片 | 国产精品久久久天天影视|