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 阅读(869) 评论(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
      回复  更多评论   


只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   博问   Chat2DB   管理


 
Copyright © Fucker Powered by: 博客园 模板提供:沪江博客