• <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 閱讀(890) 評論(2)  編輯 收藏 引用 所屬分類: 計算幾何
            Comments
            • # re: [計算幾何]pku1031
              panjinww
              Posted @ 2007-09-23 15:03
              很漂亮的代碼,可是少了注釋!!  回復  更多評論   
            • # re: [計算幾何]pku1031
              Felicia
              Posted @ 2007-09-23 15:49
              @panjinww
              呵呵,算法已經很清楚了吧,就是實現的問題了。  回復  更多評論   
             
            久久久久av无码免费网| 久久久91精品国产一区二区三区| 成人国内精品久久久久影院| 97久久久精品综合88久久| 久久久久亚洲AV片无码下载蜜桃| 伊人丁香狠狠色综合久久| 人妻少妇精品久久| 久久AV高清无码| 91精品婷婷国产综合久久| 国产一区二区精品久久凹凸| 午夜精品久久久久久久| 精品国产热久久久福利| 欧美久久久久久| 精品人妻伦九区久久AAA片69| 天天做夜夜做久久做狠狠| …久久精品99久久香蕉国产| 精品无码久久久久久久动漫| 亚洲AV无码久久| 一本久久a久久精品综合香蕉| 久久播电影网| 久久午夜无码鲁丝片| 无码人妻精品一区二区三区久久久| 精品久久久久久久| 91精品国产综合久久四虎久久无码一级 | 久久久久久国产精品美女 | 国产真实乱对白精彩久久| 伊人久久精品无码二区麻豆| 久久狠狠爱亚洲综合影院| 久久国产成人| 国产69精品久久久久9999| 国产女人aaa级久久久级| 久久久久亚洲AV无码网站| 三级三级久久三级久久| 蜜臀av性久久久久蜜臀aⅴ麻豆 | 久久天天躁狠狠躁夜夜不卡| 2021国产精品久久精品| 欧美久久久久久精选9999| 国产香蕉97碰碰久久人人| 国内精品久久久久久中文字幕| 久久精品夜夜夜夜夜久久| 国产精品久久亚洲不卡动漫|