syhd142  
日历
<2011年7月>
262728293012
3456789
10111213141516
17181920212223
24252627282930
31123456
统计
  • 随笔 - 23
  • 文章 - 122
  • 评论 - 31
  • 引用 - 0

导航

常用链接

留言簿(2)

随笔档案(23)

文章分类(270)

文章档案(122)

我的豆瓣

搜索

  •  

最新评论

阅读排行榜

评论排行榜

 
题意:给定多个城市的网络,每个城市之间的通信有花费,要求使得首都和最大城市之间的通信断掉的最小花费。要求输出砸掉的边。
解法:给出了源点和汇点,果断最小割。也许是许久没有写图论的题了,最小割输出出错,忘记了可能有多个最小割存在。
#include <stdio.h>
#include 
<stdlib.h>
#include 
<string.h>

#define N 55
#define INF 1 << 28

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

int ISAP(int src, int sink, int n)
{
    
int maxflow = 0, flow, u = src, v, d;
    
for(int i = 0; i < n; i++)
        dis[i] 
= cnt[i] = 0;
    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];
            u 
= pre[u];
        }
    }
    
return maxflow;
}

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);
    }
}

int main()
{
    
int n, m;
    
while(scanf("%d %d"&n, &m), n + m)
    {
        memset(cap, 
0sizeof(cap));
        memset(g, 
0sizeof(g));
        memset(mk, 
0sizeof(mk));
        
for(int i = 0; i < m; i++)
        {
            
int a , b, c;
            scanf(
"%d %d %d"&a, &b, &c);
            a
--, b--;
            cap[a][b] 
= cap[b][a] = c;
            g[a][b] 
= g[b][a] = 1;
        }
        
int ans = ISAP(01, n);
        Dfs(
0, n);
        
for(int i = 0; i < n; i++)
        {
            
for(int j = 0; j < n; j++)
            {
                
if(g[i][j] && !cap[i][j] && mk[i] != mk[j])
                    printf(
"%d %d\n", i + 1, j + 1);
            }
        }
        printf(
"\n");
    }
    
return 0;
}
posted on 2010-07-07 16:39 Fucker 阅读(389) 评论(0)  编辑 收藏 引用 所属分类: ACM/ICPC图论

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


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