• <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 閱讀(533) 評論(0)  編輯 收藏 引用 所屬分類: graph theory

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

            導航

            統計

            公告

            統計系統

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            久久97久久97精品免视看| 99久久精品国内| 国产香蕉久久精品综合网| 国内精品久久国产| 99久久精品国产麻豆| 国产精品xxxx国产喷水亚洲国产精品无码久久一区 | 久久午夜无码鲁丝片午夜精品| 久久99精品九九九久久婷婷| 久久99精品久久久大学生| 久久免费视频观看| 久久伊人精品一区二区三区| 国产欧美一区二区久久| 中文字幕久久精品| 亚洲国产精品婷婷久久| 国产精品久久久久久久久软件 | 久久91精品国产91久久小草| 久久青青草原精品国产不卡| 国产精品欧美久久久天天影视 | 久久久精品久久久久影院| 7777久久亚洲中文字幕| 国产精品中文久久久久久久| 四虎国产精品免费久久5151| 久久精品人人槡人妻人人玩AV| 久久久精品人妻无码专区不卡 | 无码伊人66久久大杳蕉网站谷歌| 色综合久久久久网| 久久精品aⅴ无码中文字字幕不卡 久久精品aⅴ无码中文字字幕重口 | 久久99精品久久久久久水蜜桃| 国产午夜免费高清久久影院| 亚洲精品乱码久久久久久蜜桃图片 | 久久国产精品偷99| 亚洲乱亚洲乱淫久久| 久久久久久久99精品免费观看| 久久99精品久久久久久久不卡| 麻豆av久久av盛宴av| 噜噜噜色噜噜噜久久| 日本WV一本一道久久香蕉| 久久久久久免费视频| 久久综合狠狠综合久久| 久久99精品国产自在现线小黄鸭| 中文国产成人精品久久不卡|