• <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年2月>
            303112345
            6789101112
            13141516171819
            20212223242526
            272812345
            6789101112
            統計
            • 隨筆 - 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 閱讀(876) 評論(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: 博客園 模板提供:滬江博客
            日韩欧美亚洲综合久久影院d3| 亚洲AV无码久久精品成人| 国产精品一区二区久久精品涩爱 | 亚洲AV伊人久久青青草原| 热久久最新网站获取| 99久久精品国产免看国产一区| 久久精品亚洲精品国产欧美| 99久久无色码中文字幕人妻| 亚洲国产成人久久精品影视| 综合久久一区二区三区 | 一级做a爰片久久毛片毛片| 久久精品国产亚洲av水果派| 一本一道久久精品综合| 91麻豆国产精品91久久久| 99久久免费国产精精品| 2021国内久久精品| 一本大道久久a久久精品综合| 久久精品国产亚洲AV香蕉| 热久久国产精品| 亚洲精品午夜国产VA久久成人| 精品国产乱码久久久久久浪潮| 亚洲国产另类久久久精品黑人| 国产精品日韩欧美久久综合| 久久亚洲精精品中文字幕| 久久久久久国产a免费观看不卡| 精品国产乱码久久久久久郑州公司| 久久无码AV中文出轨人妻| 精品国产91久久久久久久| 久久99热这里只有精品国产| 久久97久久97精品免视看秋霞 | 久久99热这里只有精品国产| 国产三级久久久精品麻豆三级| 欧美一区二区久久精品| 久久国产福利免费| 久久久青草青青亚洲国产免观| AV无码久久久久不卡蜜桃| 性做久久久久久免费观看| 国产三级精品久久| 狠狠色丁香婷综合久久| 久久狠狠高潮亚洲精品| 亚洲色大成网站www久久九|