hdu6214-最小割边数&最大流最小割-Smallest Minimum Cut

来源:互联网 发布:网络信息扩展平台 编辑:程序博客网 时间:2024/05/16 19:39

http://acm.hdu.edu.cn/showproblem.php?pid=6214
建边的时候每条边权 w = w * (E + 1) + 1;
这样得到最大流 maxflow / (E + 1) ,最少割边数 maxflow % (E + 1)
道理很简单,如果原先两类割边都是最小割,那么求出的最大流相等
但边权变换后只有边数小的才是最小割了
大佬的解释:
乘(E+1)是为了保证边数叠加后依然是余数,不至于影响求最小割的结果
因为假设最小割=k,那么现在新图的最小割为k*(E+1)+p,p为割的边数,本质上是,原来你割一条边,需要代价,
由于你要求边数最小 所以你多割一条边,就多一的代价,但是这个代价不足以影响到原来的代价。
原来割一条边,代价xi,现在割一条边,代价xi*A+1,只要让A>m+1,m为边数,即使割了所有的边,自己加上去的代价也就m。
用两遍dinic 一直错,不过好像先跑完,再用残量网络再跑一遍对了。。

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#define INF 1<<30#define maxn 22100#define maxm 500000using namespace std;int v[maxm],Next[maxm],w[maxm];int ww[maxn];int first[maxn],d[maxn],q[maxn];int e,S,T;int n;void init(){    e = 0;    memset(first,-1,sizeof(first));}void add_edge(int a,int b,int c){    v[e] = b;Next[e] = first[a];w[e] = c*(n+1)+1;first[a] = e++;    v[e] = a;Next[e] = first[b];w[e] = 0;first[b] = e++;    ww[e]=w[e];}int bfs(){    int rear = 0;    memset(d,-1,sizeof(d));    d[S] = 0;q[rear++] = S;    for(int i = 0;i < rear;i++){        for(int j = first[q[i]];j != -1;j = Next[j])            if(w[j] && d[v[j]] == -1){                d[v[j]] = d[q[i]] + 1;                q[rear++] = v[j];                if(v[j] == T)   return 1;            }    }    return 0;}int dfs(int cur,int maxflow){    if(cur == T)    return maxflow;    for(int i = first[cur];i != -1;i = Next[i]){        if(w[i] && d[v[i]] == d[cur] + 1)            if(int t = dfs(v[i],min(maxflow,w[i]))){                w[i] -= t;w[i^1] += t;                return t;            }    }    return 0;}int dinic(){    int ans = 0;    while(bfs()){       ans += dfs(S,INF);    }    return ans;}int main(){    int m,ncase,cas = 1;    scanf("%d",&ncase);    while(ncase--)    {        init();         //scanf("%d%d",&S,&T);        scanf("%d%d",&n,&m);                 scanf("%d%d",&S,&T);        //S = 0,T = n-1;        for(int i = 0;i < m;i++){            int a,b,c,d;            scanf("%d%d%d",&a,&b,&c);            add_edge(a,b,c);             //add_edge(b,a,c);        }        //dinic();        printf("%d\n",dinic()%(n+1));    }    return 0;}
阅读全文
0 0
原创粉丝点击