• <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
              呵呵,算法已經很清楚了吧,就是實現的問題了。  回復  更多評論   
             
            一本一本久久A久久综合精品| 91精品国产9l久久久久| 99久久精品免费看国产| Xx性欧美肥妇精品久久久久久| 久久93精品国产91久久综合| 少妇熟女久久综合网色欲| 欧美激情一区二区久久久| 91精品国产综合久久精品| 亚洲成色WWW久久网站| 久久久久久久久久久免费精品| 欧美成人免费观看久久| 香蕉久久夜色精品升级完成| 久久精品无码一区二区三区免费 | 一级做a爰片久久毛片看看| 久久精品综合网| 久久亚洲国产午夜精品理论片| 思思久久99热免费精品6| 99麻豆久久久国产精品免费| 伊人久久国产免费观看视频| 94久久国产乱子伦精品免费| 人妻精品久久无码区| 深夜久久AAAAA级毛片免费看| 国产精品一久久香蕉国产线看| 中文精品99久久国产| 精品久久国产一区二区三区香蕉| 国产三级久久久精品麻豆三级 | 久久久久亚洲AV无码专区首JN | 久久一区二区三区免费| 精品综合久久久久久97超人| 亚洲国产欧美国产综合久久| 一本大道久久东京热无码AV| 久久97久久97精品免视看秋霞 | 久久国产精品一国产精品金尊| 超级碰碰碰碰97久久久久| 久久精品国产精品亜洲毛片| 久久国产成人午夜AV影院| 久久国产福利免费| 人妻丰满?V无码久久不卡| 久久精品中文字幕大胸| 日产精品久久久久久久| 久久综合给合久久国产免费|