• <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>
            如果兩點的連線不和墻相交,那么在圖中為這兩點連一條邊,權值為這兩點的距離
            然后做 Dijkstra

            /*************************************************************************
            Author: WHU_GCC
            Created Time: 2007-8-12 19:53:33
            File Name: pku1556.cpp
            Description: 
            ***********************************************************************
            */

            #include 
            <iostream>
            #include 
            <cmath>
            #include 
            <vector>
            #include 
            <map>

            #define maxn 1010
            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-9;

            typedef 
            double weight;

            class graph_c
            {
            public:
                
            void init(int _n);
                
            void dijkstra(int S);
                
            void add_edge(int u, int v, weight w);
                weight dist[maxn];
            private:
                
            int n;
                vector 
            <int> r[maxn];
                vector 
            <weight> e[maxn];
                
            int pa[maxn];
                multimap 
            <weight, int> h;
            }
            ;

            void graph_c::init(int _n)
            {
                n 
            = _n;
                
            for (int i = 0; i < n; i++)
                
            {
                    r[i].clear();
                    e[i].clear();
                }

            }


            void graph_c::add_edge(int u, int v, weight w)
            {
                r[u].push_back(v);
                e[u].push_back(w);
            }


            void graph_c::dijkstra(int S)
            {
                weight d, tmp;
                
            int v;
                multimap
            <weight, int>::iterator it;
                h.clear();
                
            for (int i = 0; i < n; i++) dist[i] = -1;
                dist[S] 
            = 0;
                pa[S] 
            = -1;
                h.insert(multimap
            <weight, int>::value_type(0, S));
                
            while (!h.empty())
                
            {
                    it 
            = h.begin();
                    v 
            = it->second;
                    d 
            = it->first;
                    h.erase(it);
                    
            for (int i = 0; i < r[v].size(); i++)
                    
            {
                        tmp 
            = d + e[v][i];
                        
            int j = r[v][i];
                        
            if (dist[j] < 0 || tmp < dist[j])
                        
            {
                            dist[j] 
            = tmp;
                            pa[j] 
            = v;
                            h.insert(multimap
            <weight, int>::value_type(tmp, j));
                        }

                    }

                }

            }


            typedef 
            struct point_t
            {
                
            double x, y;
            }
            ;

            typedef 
            struct line_seg_t
            {
                point_t s, e;
            }
            ;

            double dist(const point_t &a, const point_t &b)
            {
                
            return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
            }


            int dblcmp(double d)
            {
                
            if (abs(d) < eps) return 0;
                
            return d > 0 ? 1 : -1;
            }


            double det(double x1, double y1, double x2, double y2)
            {
                
            return x1 * y2 - x2 * y1;
            }


            double cross(const point_t &a, const point_t &b, const point_t &c)
            {
                
            return det(b.x - a.x, b.y - a.y, c.x - a.x, c.y - a.y);
            }


            bool seg_intersect(const line_seg_t &a, const line_seg_t &b)
            {
                
            return (dblcmp(cross(a.s, b.s, b.e)) ^ dblcmp(cross(a.e, b.s, b.e))) == -2
                    
            && (dblcmp(cross(b.s, a.s, a.e)) ^ dblcmp(cross(b.e, a.s, a.e))) == -2;
            }


            line_seg_t wall[
            100];
            int cnt_wall;
            point_t p[
            100];
            int cnt_p;

            graph_c g;

            int main()
            {
                
            int n;
                
            while (scanf("%d"&n), n != -1)
                
            {
                    cnt_wall 
            = 0;
                    cnt_p 
            = 2;
                    p[
            0].x = 0.0;
                    p[
            0].y = 5.0;
                    p[
            1].x = 10.0;
                    p[
            1].y = 5.0;
                    
            for (int i = 0; i < n; i++)
                    
            {
                        
            double t1, t2, t3, t4, t5;
                        scanf(
            "%lf%lf%lf%lf%lf"&t1, &t2, &t3, &t4, &t5);
                        point_t pp;
                        pp.x 
            = t1;
                        pp.y 
            = t2;
                        p[cnt_p
            ++= pp;
                        pp.y 
            = t3;
                        p[cnt_p
            ++= pp;
                        pp.y 
            = t4;
                        p[cnt_p
            ++= pp;
                        pp.y 
            = t5;
                        p[cnt_p
            ++= pp;

                        line_seg_t t;
                        t.s.x 
            = t1;
                        t.s.y 
            = 0.0;
                        t.e.x 
            = t1;
                        t.e.y 
            = t2;
                        wall[cnt_wall
            ++= t;

                        t.s.x 
            = t1;
                        t.s.y 
            = t3;
                        t.e.x 
            = t1;
                        t.e.y 
            = t4;
                        wall[cnt_wall
            ++= t;

                        t.s.x 
            = t1;
                        t.s.y 
            = t5;
                        t.e.x 
            = t1;
                        t.e.y 
            = 10.0;
                        wall[cnt_wall
            ++= t;
                    }

                    g.init(cnt_p);
                    
            for (int i = 0; i < cnt_p; i++)
                        
            for (int j = i + 1; j < cnt_p; j++)
                        
            {
                            line_seg_t ls;
                            ls.s 
            = p[i];
                            ls.e 
            = p[j];
                            
            int flag = 1;
                            
            for (int k = 0; k < cnt_wall && flag; k++)
                                
            if (seg_intersect(ls, wall[k])) flag = 0;
                            
            if (flag)
                            
            {
                                g.add_edge(i, j, dist(p[i], p[j]));
                                g.add_edge(j, i, dist(p[i], p[j]));
                            }

                        }

                    g.dijkstra(
            0);
                    printf(
            "%.2lf\n", g.dist[1]);
                }

                
            return 0;
            }
            posted on 2007-08-13 10:34 Felicia 閱讀(479) 評論(0)  編輯 收藏 引用 所屬分類: 計算幾何
             
            新狼窝色AV性久久久久久| 一本色道久久综合狠狠躁篇| 久久久久亚洲av无码专区| 超级97碰碰碰碰久久久久最新| 久久99国产精品久久99小说| 99久久国产精品免费一区二区 | 久久精品国产半推半就| 国产精品日韩深夜福利久久 | 青青草原综合久久大伊人导航| 久久天天躁狠狠躁夜夜不卡| 久久99精品国产麻豆宅宅| 少妇久久久久久被弄到高潮| 久久天天躁狠狠躁夜夜躁2O2O| 精品水蜜桃久久久久久久| 久久精品国产99久久久古代| 久久伊人精品青青草原高清| 久久婷婷色综合一区二区| 91精品婷婷国产综合久久| 久久久久久久久无码精品亚洲日韩| 久久久久国产亚洲AV麻豆| 国产精品九九九久久九九| 亚洲va久久久噜噜噜久久狠狠| 久久久艹| 久久国产免费直播| 伊人久久综合热线大杳蕉下载| 久久精品国产亚洲av日韩| 久久99久久99精品免视看动漫| 久久www免费人成看国产片| 成人久久久观看免费毛片| 久久精品蜜芽亚洲国产AV| 久久精品国产亚洲AV香蕉| 久久亚洲电影| 久久只有这精品99| 人妻无码αv中文字幕久久琪琪布| 精品久久久久久久中文字幕| 99热都是精品久久久久久| 一本大道加勒比久久综合| 国产日产久久高清欧美一区| 久久精品国产91久久麻豆自制| 久久精品国产亚洲AV麻豆网站| 97r久久精品国产99国产精|