• <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 閱讀(895) 評論(2)  編輯 收藏 引用 所屬分類: 計算幾何
            Comments
            • # re: [計算幾何]pku1031
              panjinww
              Posted @ 2007-09-23 15:03
              很漂亮的代碼,可是少了注釋!!  回復  更多評論   
            • # re: [計算幾何]pku1031
              Felicia
              Posted @ 2007-09-23 15:49
              @panjinww
              呵呵,算法已經很清楚了吧,就是實現的問題了。  回復  更多評論   
             
            午夜欧美精品久久久久久久| 亚洲一区中文字幕久久| 久久天天躁狠狠躁夜夜avapp| 久久精品欧美日韩精品| 蜜臀久久99精品久久久久久小说 | 久久人人爽人人爽人人片AV麻豆| 久久99国产精品久久99| 国产成人无码精品久久久久免费 | 久久午夜伦鲁片免费无码| 久久九九亚洲精品| 亚洲性久久久影院| 久久精品草草草| 99久久这里只精品国产免费| 久久精品免费一区二区三区| 久久久午夜精品| 无码人妻久久一区二区三区| 国产成人久久激情91| 2019久久久高清456| 国产精品xxxx国产喷水亚洲国产精品无码久久一区 | 久久99国产精一区二区三区| 久久影视综合亚洲| 精品久久无码中文字幕| 欧美一区二区久久精品| 久久综合综合久久97色| 久久99热这里只有精品国产| 大蕉久久伊人中文字幕| 久久棈精品久久久久久噜噜| 久久笫一福利免费导航| 精品久久久久久无码国产| 99久久久精品| 99久久99久久久精品齐齐 | 麻豆精品久久精品色综合| 一个色综合久久| 亚洲欧洲久久av| 中文字幕久久亚洲一区| 久久无码国产| 久久综合偷偷噜噜噜色| 久久无码国产专区精品| 亚洲欧美伊人久久综合一区二区| 亚洲精品乱码久久久久久蜜桃不卡| 亚洲欧洲中文日韩久久AV乱码|