• <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>
            syhd142  
            日歷
            <2011年11月>
            303112345
            6789101112
            13141516171819
            20212223242526
            27282930123
            45678910
            統計
            • 隨筆 - 23
            • 文章 - 122
            • 評論 - 31
            • 引用 - 0

            導航

            常用鏈接

            留言簿(2)

            隨筆檔案(23)

            文章分類(270)

            文章檔案(122)

            我的豆瓣

            搜索

            •  

            最新評論

            閱讀排行榜

            評論排行榜

             
            很裸的題目,即求所有點對之間的最小割。有Gomory-Hu和Gusfield兩種算法,貌似Gomory-Hu的難實現一點,效率搞一點,Gusfield的容易實現一點,但是效率沒那么高?
            看著網上的偽代碼寫出來的。。具體論文還沒有仔細看。
            貼出論文Dan Gusfield: Very Simple Methods for All Pairs Network Flow Analysis. SIAM J. Comput. 19(1): 143-155 (1990)
            #include <iostream>
            using namespace std;

            #define N 205
            #define INF 1 << 28

            int g[N][N], mk[N];
            int cap[N][N], dis[N], pre[N], cnt[N];

            int ISAP(int src, int sink, int n)
            {
                
            int maxflow = 0, flow, u = src, v, d;
                memset(dis, 
            0sizeof(dis));
                memset(cnt, 
            0sizeof(cnt));
                pre[src] 
            = src;
                cnt[src] 
            = n;
                
            while(dis[src] < n)
                {
                    
            if(u == sink)
                    {
                        flow 
            = INF;
                        
            for(int i = sink; i != src; i = pre[i])
                            
            if(cap[pre[i]][i] < flow)
                                flow 
            = cap[pre[i]][i];
                        maxflow 
            += flow;
                        u 
            = src;
                        
            for(int i = sink; i != src; i = pre[i])
                        {
                            cap[pre[i]][i] 
            -= flow;
                            cap[i][pre[i]] 
            += flow;
                        }
                    }
                    
            for(d = n, v = 0; v < n; v++)
                    {
                        
            if(!cap[u][v]) continue;
                        d 
            = d < dis[v] ? d : dis[v];
                        
            if(dis[v] + 1 == dis[u]) break;
                    }
                    
            if(v < n) pre[v] = u, u = v;
                    
            else
                    {
                        
            if(!(--cnt[dis[u]])) break;
                        
            ++cnt[dis[u] = d + 1];
                        
            if(u != src) u = pre[u];
                    }
                }
                
            return maxflow;
            }

            int parent[N], cut[N][N];

            void dfs(int u, int n)
            {
                mk[u] 
            = 1;
                
            for(int i = 0; i < n; i++)
                {
                    
            if(!mk[i] && cap[u][i])
                        dfs(i, n);
                }
            }

            void Solve(int n)
            {
                memset(parent, 
            0sizeof(parent));
                memset(cut, 
            127sizeof(cut));
                
            for(int i = 1; i < n; i++)
                {
                    memset(mk, 
            0sizeof(mk));
                    memcpy(cap, g, 
            sizeof(g));
                    
            int mincut = ISAP(i, parent[i], n);
                    dfs(i, n);
                    
            for(int j = i + 1; j < n; j++)
                        
            if(mk[j] && parent[j] == parent[i])
                        {
                            parent[j] 
            = i;
                        }
                    cut[i][parent[i]] 
            = cut[parent[i]][i] = mincut;
                    
            for(int j = 0; j < i; j++)
                        cut[i][j] 
            = cut[j][i] = min(mincut, cut[parent[i]][j]);
                }
                
                
            for(int i = 0; i < n; i++)
                {
                    
            for(int j = 0; j < n; j++)
                    {
                        
            if(i == j) cut[i][j] = 0;
                        printf(
            "%d", cut[i][j]);
                        
            if(j != n - 1) printf(" ");
                    }
                    printf(
            "\n");
                }
            }

            int main()
            {
                
            int t, n, cas = 1;
                scanf(
            "%d"&t);
                
            while(t--)
                {
                    scanf(
            "%d"&n);
                    
            for(int i = 0; i < n; i++)
                    
            for(int j = 0; j < n; j++)
                        scanf(
            "%d"&g[i][j]);
                    printf(
            "Case #%d:\n", cas++);
                    Solve(n);
                }
                
            return 0;
            }
            posted on 2010-07-09 17:38 Fucker 閱讀(886) 評論(3)  編輯 收藏 引用 所屬分類: ACM/ICPC圖論
            評論:
            • # re: UVA 11594 All Pairs Maximum Flow  Ziwen Posted @ 2011-02-12 11:29
              能不能把論文給我發一份
              THX
              Very Simple Methods for All Pairs Network Flow Analysis
              caocao9926@163.com  回復  更多評論   

            • # re: UVA 11594 All Pairs Maximum Flow  xiaohuo Posted @ 2011-11-20 13:01
              xiaohuo_2008@yahoo.cn 論文能發給我一下嗎?跪求啊!
              Very Simple Methods for All Pairs Network Flow Analysis
              謝謝!  回復  更多評論   

            • # re: UVA 11594 All Pairs Maximum Flow  Fucker Posted @ 2011-11-20 21:33
              論文找不著了,自己去教育網下吧。。@xiaohuo
                回復  更多評論   

             
            Copyright © Fucker Powered by: 博客園 模板提供:滬江博客
            色欲综合久久躁天天躁蜜桃| 久久香蕉国产线看观看99| 一级做a爰片久久毛片看看| 久久精品亚洲欧美日韩久久| 久久人人爽人人爽人人片AV高清| AV狠狠色丁香婷婷综合久久| 狠狠色伊人久久精品综合网| 99蜜桃臀久久久欧美精品网站| 久久噜噜电影你懂的| 中文字幕无码免费久久| 久久久久人妻一区精品| 久久久久人妻一区精品色 | 91精品国产91久久久久久| 亚洲欧美精品一区久久中文字幕| 91精品国产高清91久久久久久| 国内精品久久久久影院老司| 国产L精品国产亚洲区久久| 麻豆成人久久精品二区三区免费| 中文精品99久久国产| 久久99精品久久久久久野外| 女人香蕉久久**毛片精品| 亚洲欧美伊人久久综合一区二区| 波多野结衣久久一区二区| 国产69精品久久久久99尤物| 国产麻豆精品久久一二三| 久久亚洲精精品中文字幕| 亚洲国产精品无码成人片久久| 久久精品国产亚洲AV香蕉| 九九精品久久久久久噜噜| 青青草原综合久久大伊人导航| 久久久久这里只有精品| 久久精品免费网站网| 久久天天躁狠狠躁夜夜av浪潮| 久久精品国产精品亚洲下载| 国产精品美女久久久久av爽| 国产午夜精品久久久久九九| 久久激情亚洲精品无码?V| 欧美大战日韩91综合一区婷婷久久青草 | 亚洲欧洲日产国码无码久久99| 中文国产成人精品久久不卡| 国产婷婷成人久久Av免费高清|