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

            Stoer_Wagner 算法,求無向圖的最小割

              prim算法不僅僅可以求最小生成樹,也可以求“最大生成樹”。最小割集Stoer-Wagner算法就是典型的應用實例。

                求解最小割集普遍采用Stoer-Wagner算法,不提供此算法證明和代碼,只提供算法思路:

            1.min=MAXINT,固定一個頂點P

            2.從點P用類似prim的s算法擴展出“最大生成樹”,記錄最后擴展的頂點和最后擴展的邊

            3.計算最后擴展到的頂點的切割值(即與此頂點相連的所有邊權和),若比min小更新min

            4.合并最后擴展的那條邊的兩個端點為一個頂點(當然他們的邊也要合并,這個好理解吧?)

            5.轉到2,合并N-1次后結束

            6.min即為所求,輸出min

            prim本身復雜度是O(n^2),合并n-1次,算法復雜度即為O(n^3)

            如果在prim中加堆優化,復雜度會降為O((n^2)logn)


            #include <cmath>

            #include 
            <cstdio>

            #include 
            <memory.h>

            #include 
            <algorithm>

            #include 
            <iomanip>

            #include 
            <iostream>

            #include 
            <vector>

            #include 
            <string>

            #include 
            <queue>

             

            using namespace std;

             

            const int N = 500 + 3;

             

            int n, m;

            int mat[N][N];

            int dist[N];

            int visited[N];

            int del[N];  // true表示該點已經被刪掉

             

            // 結點~n

            int Stoer_Wagner()

            {

                 
            int minCut = INT_MAX;  // 無向圖最小割

                 
            int tmp;

                 
            int i, t, j, k, pre;

                 
            int s = 1;   // 源點

                 memset(del, 
            0sizeof(del));

             

                 
            for (t = 1; t < n; t++)  // n - 1次Maximum Adjacency Search

                 {

                     
            for (i = 1; i <= n; i++)

                          
            if (!del[i])

                               dist[i] 
            = mat[s][i];

             

                     memset(visited, 
            0sizeof(visited));

                     visited[s] 
            = 1;

                     k 
            = s;

                     
            for (i = 1; i <= n - t; i++)  // 每次剩下n - t + 1個結點

                     {

                          tmp 
            = -1e9;

                          pre 
            = k;

                          k 
            = 0;

                          
            for (j = 1; j <= n; j++)

                          {

                               
            if (!del[j] && !visited[j] && dist[j] > tmp)

                               {

                                   k 
            = j;

                                   tmp 
            = dist[j];

                               }

                          }

                          
            if (!k) return 0;  // 不連通

             

                          visited[k] 
            = 1;

                          
            for (j = 1; j <= n; j++)

                               
            if (!del[j] && !visited[j])

                                   dist[j] 
            += mat[k][j];

                     }

             

                     minCut 
            = min(minCut, dist[k]);

                     del[k] 
            = 1;  // 刪除k點

             

                     
            // 合并k點和源點

                     

                     
            for (i = 1; i <= n; i++)

                          
            if (!del[i] && i != pre)

                          {

                               mat[pre][i] 
            += mat[k][i];

                               mat[i][pre] 
            = mat[pre][i];

                          }

                 }

             

                 
            return minCut;

            }

             

            int main ()

            {

                 
            int u, v, w, i;

                 
            while (scanf("%d%d"&n, &m) != EOF)

                 {

                     memset(mat, 
            0sizeof(mat));

                     
            while (m--)

                     {

                          scanf(
            "%d%d%d"&u, &v, &w);

                          
            if (u == v) continue;  

                          mat[u 
            + 1][v + 1+= w;

                          mat[v 
            + 1][u + 1+= w;

                     }

                     printf(
            "%d\n", Stoer_Wagner());

                 }

            }


            posted on 2010-10-25 23:46 yzhw 閱讀(548) 評論(0)  編輯 收藏 引用 所屬分類: graph theory

            <2010年12月>
            2829301234
            567891011
            12131415161718
            19202122232425
            2627282930311
            2345678

            導航

            統計

            公告

            統計系統

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            日韩精品久久久肉伦网站| 伊人久久综在合线亚洲2019| 久久久精品久久久久影院| 午夜精品久久久久久久无码| 亚洲va久久久久| 久久97精品久久久久久久不卡| 99久久精品国产综合一区| 久久久久久伊人高潮影院| 大伊人青草狠狠久久| 欧美成a人片免费看久久| 久久久久久毛片免费播放| 成人午夜精品久久久久久久小说| 热久久最新网站获取| 嫩草影院久久99| 午夜精品久久久久久久久| 久久久久国产| 91精品国产91久久久久久青草| 国内精品人妻无码久久久影院导航| 色综合合久久天天综合绕视看| 无码人妻少妇久久中文字幕蜜桃| 久久婷婷五月综合成人D啪| 国产精品久久影院| 日韩久久久久久中文人妻 | 久久人人爽人人人人爽AV| 国产精品久久国产精品99盘| 久久婷婷五月综合国产尤物app| 精品视频久久久久| 久久婷婷久久一区二区三区| 久久精品国产99久久无毒不卡 | 日本三级久久网| 狠狠色丁香久久综合五月| 久久精品无码一区二区WWW| 色综合合久久天天给综看| 久久久久久A亚洲欧洲AV冫| 久久综合狠狠色综合伊人| 久久久久夜夜夜精品国产| 69久久精品无码一区二区| 久久99精品国产99久久6男男| WWW婷婷AV久久久影片| 久久久久免费精品国产| 久久国产影院|