• <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中加堆優(yōu)化,復雜度會降為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 閱讀(542) 評論(0)  編輯 收藏 引用 所屬分類: graph theory

            <2012年2月>
            2930311234
            567891011
            12131415161718
            19202122232425
            26272829123
            45678910

            導航

            統(tǒng)計

            公告

            統(tǒng)計系統(tǒng)

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            国产成人综合久久精品红| 久久人妻少妇嫩草AV无码蜜桃| 99久久这里只精品国产免费| 人妻无码αv中文字幕久久琪琪布| 2021国产精品久久精品| 久久久久综合网久久| 日本亚洲色大成网站WWW久久| 亚洲欧美日韩中文久久| 国产精品美女久久久久AV福利 | 99精品国产在热久久无毒不卡| 国内精品久久久久影院一蜜桃| 国产亚洲成人久久| 99久久精品国产一区二区| 久久99精品久久久久久噜噜| 久久精品无码一区二区无码| 亚洲国产成人久久综合碰| 1000部精品久久久久久久久| 久久只这里是精品66| 久久AAAA片一区二区| 2021少妇久久久久久久久久| 综合人妻久久一区二区精品| 久久97久久97精品免视看| 久久精品男人影院| 久久青青草原亚洲av无码app| 国产美女亚洲精品久久久综合 | 中文精品久久久久人妻不卡| 国产精品成人无码久久久久久| AV无码久久久久不卡网站下载 | 久久亚洲日韩看片无码| 热久久国产欧美一区二区精品| 国产AⅤ精品一区二区三区久久| 日韩精品久久久肉伦网站| 亚洲国产精品18久久久久久| 亚洲国产精品无码久久久久久曰| 色婷婷狠狠久久综合五月| 久久久国产精品| 亚洲精品久久久www| 免费久久人人爽人人爽av| 国产精品美女久久福利网站| 亚洲va久久久久| 久久精品国产亚洲AV高清热|