HDU6214 Smallest Minimum Cut 【最大流求最小割边】

来源:互联网 发布:淘宝网白色运动鞋 编辑:程序博客网 时间:2024/05/17 01:02

传送门


因为最大流=最小割边的最大流量限制之和
因为m<=1000
将每条边流量*2000+1
不难发现 假如原图有2个最小割,大小分别为x,y(x < y)
改流量后只有x能流满
跑出来的最大流%2000=最小割边数量

#include<stdio.h>#include<bits/stdc++.h>#define ll long long#define pii pair<int,int>#define pll pair<ll,ll>#define MEM(a,x) memset(a,x,sizeof(a))#define lowbit(x) ((x)&-(x))using namespace std;const int INF = 1e9+7;const int inf=INF;const int N = 250 + 5;//点数const int M = 2e3 + 50;//边数*2(包括反向边)int level[N];//标号 level[i]:s到i的最短距离struct Edge{    int to,c,next;}edge[M];int head[N];inline void add_edge(int k,int u,int v,int c){    edge[k].to = v;    edge[k].c = c;    edge[k].next = head[u];    head[u] = k;}bool bfs(int s,int t,int n){//标号 计算level    deque<int>que;    fill(level,level+n+1,-1);    que.push_back(s);    level[s] = 0;    while(!que.empty()){        int u = que.front();        if(u == t){            return true;        }        que.pop_front();        for(int i = head[u];i!=-1;i = edge[i].next){            if(edge[i].c > 0 && level[edge[i].to] == -1){                level[edge[i].to] = level[u] + 1;                que.push_back(edge[i].to);            }        }    }    return false;}int dfs(int u,int t,int maxf){//u:所在的点 t:汇点 maxf:能流到u的流量    if(u == t){        return maxf;    }    int sum = 0;    for(int i = head[u];i!=-1;i = edge[i].next){        Edge&e = edge[i];        if(e.c>0 && level[e.to]>level[u]){            int f = dfs(e.to,t,min(maxf - sum,e.c));            sum += f;            edge[i].c -= f;            edge[i^1].c += f;            if(sum == maxf){//流量用完了                break;            }        }    }    level[u]=-1;    return sum;}int dinic(int s,int t,int n){//s:源点 t:汇点 n:点数    int ans = 0;    while(bfs(s,t,n)){        ans += dfs(s,t,INF);    }    return ans;}int slove(int n,int m){    fill(head,head+n+1,-1);    int S,T;    scanf("%d%d",&S,&T);    int nume=0;    for(int i=0;i<m;++i){        int u,v,w;        scanf("%d%d%d",&u,&v,&w);        add_edge(nume++,u,v,w*2000+1);        add_edge(nume++,v,u,0);    }    return dinic(S,T,n+1)%2000;}int main(){    //freopen("/home/lu/code/r.txt","r",stdin);    //freopen("/home/lu/code/w.txt","w",stdout);    int T;    scanf("%d",&T);    while(T--){        int n,m;        scanf("%d%d",&n,&m);        printf("%d\n",slove(n,m));    }    return 0;}
阅读全文
0 0
原创粉丝点击