• <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
              呵呵,算法已經很清楚了吧,就是實現的問題了。  回復  更多評論   
             
            久久精品国产精品青草app| 伊人伊成久久人综合网777| 久久91亚洲人成电影网站| 久久精品一区二区三区不卡| 精品久久久无码中文字幕天天| 久久笫一福利免费导航 | 亚洲欧美一区二区三区久久| 无码伊人66久久大杳蕉网站谷歌| 久久精品国产亚洲av水果派| 国内精品久久久久久麻豆| 亚洲色婷婷综合久久| 久久99精品久久久久久野外 | 精品国产乱码久久久久久浪潮| 亚洲精品视频久久久| 伊人丁香狠狠色综合久久| 伊人久久亚洲综合影院| 国产精品久久久99| 精品久久久噜噜噜久久久 | 精品永久久福利一区二区| 日本精品久久久久影院日本| 国产一区二区精品久久| 无码国内精品久久人妻蜜桃| 狠狠色丁香婷婷久久综合五月 | 要久久爱在线免费观看| 国产亚洲成人久久| 97超级碰碰碰碰久久久久| 久久福利青草精品资源站免费 | A狠狠久久蜜臀婷色中文网| 久久国产免费直播| 久久久久久综合网天天| 99久久这里只精品国产免费| 亚洲精品国产综合久久一线| 香蕉99久久国产综合精品宅男自 | 热re99久久6国产精品免费| 久久青青色综合| 久久妇女高潮几次MBA| 综合久久国产九一剧情麻豆| 久久久噜噜噜久久中文字幕色伊伊 | 久久精品国产99国产精品| 国内精品久久久久久中文字幕 | 久久99国产精品久久久|