• <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>

            O(1) 的小樂

            Job Hunting

            公告

            記錄我的生活和工作。。。
            <2010年10月>
            262728293012
            3456789
            10111213141516
            17181920212223
            24252627282930
            31123456

            統計

            • 隨筆 - 182
            • 文章 - 1
            • 評論 - 41
            • 引用 - 0

            留言簿(10)

            隨筆分類(70)

            隨筆檔案(182)

            文章檔案(1)

            如影隨形

            搜索

            •  

            最新隨筆

            最新評論

            閱讀排行榜

            評論排行榜

            最短路Dijkstra算法

            Dijkstra算法

            普通實現復雜度O(n^2)  BinaryHeaps優化的是O(mlogn) Fib堆優化的是O(m+nlogn)   m是圖的邊數,n是圖的點數

            特點: 基本思想是貪心思想,對于每個節點v[i],維護估計最短路長度最大值 每次取出一個使得該估計值最小的t,并采用與t相連的邊對其余點的估計值進行更新,更新后不再考慮t。在此過程中,估計值單調遞減,所以可以得到確切的最短路。

            使用松弛技術,

            優點: 簡單,實現快速!  具有貪心算法的特性:最短路徑的子路徑是最短路徑!   相比于和他同類型的Bellman ford,它更快,缺點是邊權值非負。

            缺點: 邊的權值必須為非負,Dijkstra算法假定輸入圖中的所有邊的權值都是非負的

             

            樸素Dijkstra算法實現:

            #include <iostream>
            #include <vector>
            #include <algorithm>
            using namespace std;

            #define INF 329999   //定義此圖的最大距離
            #define MaxN 100           //一個圖最大的點數

            int cost[MaxN][MaxN];                      //此為Dijkstra的數組實現,復雜度為O(n^2)
            int pathCon[MaxN];                         //記錄最短路徑獲得的時候的前向節點 路徑重建
            int pathDis[MaxN];                         //路徑距離
            int vis[MaxN];                             //標記是否結束訪問

            void Dijkstra(int n, int root)

                memset(vis, 0, sizeof(vis));     //標記vis為全部未訪問
                vis[root] = 1;                   //標記root已訪問
                for (int i=0; i<n; i++)             
                    pathDis[i] = cost[root][i], pathCon[i] = root;
                pathDis[root] = 0; pathCon[root] = -1; // 樹根的標記
                int pre = root;

                for (int i=1; i<n; i++)        //遍歷剩下的n-1個點,每執行一次能夠確定一個點的最短距離
                {
                    int Min = INF;
                    for (int j=0; j<n; j++)
                    {
                        if (vis[j]==0 && cost[pre][j]!=INF && pathDis[pre]+cost[pre][j]<pathDis[j])
                            pathDis[j] = pathDis[pre] + cost[pre][j],pathCon[j] = pre; 
                    }
                    for (int j=0; j<n; j++)
                    {
                        if (vis[j] == 0 && pathDis[j] < Min)
                            Min = pathDis[j], pre = j;
                    }
                    vis[pre] = 1;
                }
            }
            void init(int n)                      //初始化構建矩陣的時候,不相鄰元素距離為INF,自身到自身距離為0
            {
                //n=4
                cost[0][0]=0,  cost[0][1]=5,cost[0][2]=1,  cost[0][3]=4,  cost[0][4]=INF;
                cost[1][0]=5,  cost[1][1]=0,cost[1][2]=3,  cost[1][3]=1,  cost[1][4]=2;
                cost[2][0]=1,  cost[2][1]=3,cost[2][2]=0,  cost[2][3]=INF,cost[2][4]=6;
                cost[3][0]=4,  cost[3][1]=1,cost[3][2]=INF,cost[3][3]=0,  cost[3][4]=3;
                cost[4][0]=INF,cost[4][1]=2,cost[4][2]=6,  cost[4][3]=3,  cost[4][4]=0;

            }
            void reshape(int n,int root)
            {
                for(int i=0;i<n;i++)
                {
                    if(pathDis[i]==INF) cout<<"node "<<root<<" to node "<<i<<"is INF"<<endl;
                    else cout<<"node "<<root<<" to node "<<i<<"is "<<pathDis[i]<<endl;
                }
                vector<int> Path[MaxN];            //構建每條最短路徑
                for(int i=n-1;i>=0;i--)
                {
                    vector<int> temp;
                    int j=i;temp.push_back(j);
                    while(pathCon[j]!=-1)
                    {
                        temp.push_back(pathCon[j]);
                        j=pathCon[j];
                    }
                    reverse(temp.begin(),temp.end());
                    Path[i]=temp;
                }

                for(int i=0;i<n;i++)
                {
                    cout<<"Path "<<root<<"    "<<"to "<<i<<"  is ";
                    for(int j=0;j<Path[i].size();j++)
                        cout<<Path[i][j]<<"   ";
                    cout<<endl;
                }
            }
            int main()     //n個點,所以n是5
            {
                init(5);
                Dijkstra(5,0);
                reshape(5,0);
            }

             

             

             

            加入優先隊列優化的Dijkstra算法:

            #include <iostream>
            #include <vector>
            #include <algorithm>
            #include <queue>
            using namespace std;

            // Dijkstra O(m*log n)

            #define MaxN 100           //一個圖最大的點數
            #define MaxM 1000          //一個圖最大的邊數
            #define V   MaxN
            #define E   MaxM
            #define INF 329999         //定義此圖的最大距離
            #define typec int                          // type of cost

            typec cost[E], pathDis[V];                 // cost[E] 每條邊的代價   pathDis[V]最短路徑的值
            int e, pnt[E], cluster[E], head[V];
            int pathCon[V], vis[V];
            /*
            pnt[]  保存每條邊的終點
            cost[] 保存每條邊的權值
            cluster[]  把邊按照起始點聚類
            head[u] u號節點的邊結束邊序號
            */
            struct qnode {
                int v; typec c;
                qnode (int vv = 0, typec cc = 0) : v(vv), c(cc) {}
                bool  operator < (const qnode& r)  const {  return c>r.c; }
            };
            void dijkstra(int n, const int root)
            {
                qnode mv;
                priority_queue<qnode> que;
                vis[root] = 1; pathDis[root] = 0;
                que.push(qnode(root, 0));
                for (int pre = root, i=1; i<n; i++)                  //循環n次
                {
                    for (int j = head[pre]; j != -1; j = cluster[j])    //展開節點
                    {
                        int k = pnt[j];                              //j為邊索引號
                        if (vis[k] == 0 &&
                            pathDis[pre] + cost[j] < pathDis[k])
                        {
                                pathDis[k] = pathDis[pre] + cost[j];
                                que.push(qnode(pnt[j], pathDis[k]));
                                pathCon[k] = pre;
                        }
                    }
                    while (!que.empty() && vis[que.top().v] == 1)
                        que.pop();
                    if (que.empty()) break;
                    mv = que.top(); que.pop();
                    vis[pre = mv.v] = 1;
                }
            }

            inline void addedge(int u, int v, typec c){
                pnt[e] = v; cost[e] = c; cluster[e] = head[u]; head[u] = e++;
            }
            void init(int n, int m)           //n為點數, m 為邊數
            {
                int i, u, v; typec c;
                e = 0;
                memset(head, -1, sizeof(head));
                memset(vis, 0, sizeof(vis));
                memset(pathCon, -1, sizeof(pathCon));
                for (i = 0; i < n; i++) pathDis[i] = INF;
                //for (i = 0; i < n; ++i)
                //{
                    //scanf("%d%d%d", &u, &v, &c); // %d: type of cost
                    //addedge(u, v, c);            // 這里添加的是單向邊,如果題目描述不同,需要加入addedge(v,u,c);
                //}
                    addedge(0,2,1);
                    addedge(2,0,1);
                    addedge(0,1,4);
                    addedge(1,0,4);
                    addedge(1,3,1);
                    addedge(3,1,1);
                    addedge(2,3,4);
                    addedge(3,2,4);
                    addedge(2,4,2);
                    addedge(4,2,2);
                    addedge(3,5,2);
                    addedge(6,3,2);
                    addedge(4,5,3);
                    addedge(5,4,3);
            }

            void reshape(int n,int root)
            {
                for(int i=0;i<n;i++)
                {
                    if(pathDis[i]==INF) cout<<"node "<<root<<" to node "<<i<<"is INF"<<endl;
                    else cout<<"node "<<root<<" to node "<<i<<"is "<<pathDis[i]<<endl;
                }
                vector<int> Path[MaxN];            //構建每條最短路徑
                for(int i=n-1;i>=0;i--)
                {
                    vector<int> temp;
                    int j=i;temp.push_back(j);
                    while(pathCon[j]!=-1)
                    {
                        temp.push_back(pathCon[j]);
                        j=pathCon[j];
                    }
                    reverse(temp.begin(),temp.end());
                    Path[i]=temp;
                }

                for(int i=0;i<n;i++)
                {
                    cout<<"Path "<<root<<"    "<<"to "<<i<<"  is ";
                    for(int j=0;j<Path[i].size();j++)
                        cout<<Path[i][j]<<"   ";
                    cout<<endl;
                }
            }

            int main()
            {
                init(6,7); //n = 6 m = 7
                dijkstra(6, 3);
                reshape(6,3);
            }

             

            上述實現的圖是:

            Dijkstra

            用LatexDraw 畫的圖,效果還可以。。。

            posted on 2010-10-14 13:47 Sosi 閱讀(1503) 評論(0)  編輯 收藏 引用

            統計系統
            精品久久久久久久中文字幕| 久久精品国产亚洲AV忘忧草18| 亚洲日韩中文无码久久| 99久久99久久精品免费看蜜桃| 色婷婷综合久久久久中文字幕| a级成人毛片久久| 国产精品亚洲美女久久久| 精品久久久噜噜噜久久久 | 伊人丁香狠狠色综合久久| 亚洲精品无码久久毛片| 情人伊人久久综合亚洲| 久久精品国产亚洲av高清漫画| 欧美伊人久久大香线蕉综合| 午夜精品久久久内射近拍高清| 久久国产免费观看精品3| 亚洲伊人久久成综合人影院 | 精品久久久久久中文字幕| 亚洲va久久久噜噜噜久久男同| 久久久久亚洲av成人网人人软件 | 久久精品一区二区三区中文字幕| 国产精品久久国产精品99盘| 无码国内精品久久人妻麻豆按摩| 国产精品一久久香蕉产线看| 亚洲国产日韩欧美综合久久| 国产精品免费久久久久影院| 久久久噜噜噜久久熟女AA片| 亚洲伊人久久成综合人影院 | 久久久国产99久久国产一| 婷婷久久五月天| 久久精品中文字幕第23页| 精品久久久久久久久中文字幕| 国内精品九九久久久精品| 亚洲∧v久久久无码精品| 中文精品久久久久人妻| 日韩va亚洲va欧美va久久| 青青草原综合久久大伊人导航 | 99国产精品久久| 国产欧美久久一区二区| 99久久精品九九亚洲精品| 秋霞久久国产精品电影院| 国产精品成人久久久久久久|