上下界网络流总结

最新推荐文章于 2020-04-05 18:28:53 发布
最新推荐文章于 2020-04-05 18:28:53 发布
阅读量112 收藏
点赞数
原文链接：http://www.cnblogs.com/mleautomaton/p/11240916.html
版权

前言

最近发现网络流很菜,所以就做点题目练练手.

无源汇上下界可行流

题面

我们不妨先把所有的边的容量设为\([0,upp-low]\),那么这个时候显然一条边的实际流量应该是\(flow_i+low_i\)对吧.

我们对于新建的图跑流,肯定是能够保证\(flow_i\)平衡,但是我们无法保证\(flow_i+low_i\)平衡,这时我们引入一个\(delta\)数组.

对于\(delta_i\),它的值是所有出边的下界-入边的下界.

那么这个时候我们对于\(delta_i>0\)的就只需要增加它的入边,\(delta<0\)的增加它的出边.

这时我们引入超级源点和超级汇点,然后直接按照上述方法连边即可.

/*
  mail: mleautomaton@foxmail.com
  author: MLEAutoMaton
  This Code is made by MLEAutoMaton
*/
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<iostream>
using namespace std;
#define ll long long
#define re register
#define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
inline int gi(){
    int f=1,sum=0;char ch=getchar();
    while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
    return f*sum;
}
const int N=50010,Inf=1e9+10;
int front[N],cnt,n,m,s,t,dep[N],in[N],out[N];
queue<int>Q;
struct node{int to,nxt,w;}e[N<<1];
void Add(int u,int v,int w){
    e[cnt]=(node){v,front[u],w};front[u]=cnt++;
    e[cnt]=(node){u,front[v],0};front[v]=cnt++;
}
bool bfs(){
    Q.push(s);memset(dep,0,sizeof(dep));dep[s]=1;
    while(!Q.empty()){
        int u=Q.front();Q.pop();
        for(int i=front[u];~i;i=e[i].nxt){
            int v=e[i].to;
            if(!dep[v] && e[i].w){
                dep[v]=dep[u]+1;Q.push(v);
            }
        }
    }
    return dep[t];
}
int dfs(int u,int flow){
    if(u==t || !flow)return flow;
    for(int i=front[u];~i;i=e[i].nxt){
        int v=e[i].to;
        if(dep[v]==dep[u]+1 && e[i].w){
            int di=dfs(v,min(flow,e[i].w));
            if(di){
                e[i].w-=di;e[i^1].w+=di;return di;
            }
            else dep[v]=0;
        }
    }
    return 0;
}
int Dinic(){
    int flow=0;
    while(bfs())while(int d=dfs(s,Inf))flow+=d;
    return flow;
}
int sum,Up[N];
int main(){
#ifndef ONLINE_JUDGE
    freopen("in.in","r",stdin);
#endif
    n=gi();m=gi();memset(front,-1,sizeof(front));
    s=0;t=n+1;
    for(int i=1;i<=m;i++){
        int u=gi(),v=gi(),low=gi(),upp=gi();
        in[v]+=low;out[u]+=low;Up[i]=upp;
        Add(u,v,upp-low);
    }
    for(int i=1;i<=n;i++)
        if(in[i]>out[i]){Add(s,i,in[i]-out[i]);sum+=in[i]-out[i];}
        else Add(i,t,out[i]-in[i]);
    int flow=Dinic();
    if(flow<sum)return puts("NO"),0;
    puts("YES");
    for(int i=1;i<=m;i++)
        printf("%d\n",Up[i]-e[(i-1)*2].w);
    return 0;
}

有源汇上下界最大流

题面

还是先按照上述的方法,将\(t->s\)连边\(\infty\)保证流量平衡,然后跑出来一个可行流.

接着在残量网络上面跑\(s->t\)的最大流就是答案.

/*
  mail: mleautomaton@foxmail.com
  author: MLEAutoMaton
  This Code is made by MLEAutoMaton
*/
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<iostream>
using namespace std;
#define ll long long
#define re register
#define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
inline int gi(){
    int f=1,sum=0;char ch=getchar();
    while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
    return f*sum;
}
const int N=100010,Inf=1e9+10;
int n,m,s,t,front[N],cnt,dep[N],ss,tt,delta[N],sum;queue<int>Q;
struct node{int to,nxt,w;}e[N<<1];
void Add(int u,int v,int w){
    e[cnt]=(node){v,front[u],w};front[u]=cnt++;
    e[cnt]=(node){u,front[v],0};front[v]=cnt++;
}
bool bfs(){
    Q.push(ss);memset(dep,0,sizeof(dep));dep[ss]=1;
    while(!Q.empty()){
        int u=Q.front();Q.pop();
        for(int i=front[u];~i;i=e[i].nxt){
            int v=e[i].to;
            if(!dep[v] && e[i].w){
                dep[v]=dep[u]+1;Q.push(v);
            }
        }
    }
    return dep[tt];
}
int dfs(int u,int flow){
    if(u==tt || !flow)return flow;
    for(int i=front[u];~i;i=e[i].nxt){
        int v=e[i].to;
        if(dep[v]==dep[u]+1 && e[i].w){
            int di=dfs(v,min(e[i].w,flow));
            if(di){
                e[i].w-=di;e[i^1].w+=di;return di;
            }
            else dep[v]=0;
        }
    }
    return 0;
}
int Dinic(){
    int flow=0;
    while(bfs())while(int d=dfs(ss,Inf))flow+=d;
    return flow;
}
int main(){
    n=gi();m=gi();s=gi();t=gi();memset(front,-1,sizeof(front));
    ss=n+1;tt=ss+1;
    for(int i=1;i<=m;i++){
        int u=gi(),v=gi(),low=gi(),upp=gi();
        Add(u,v,upp-low);
        delta[u]-=low;delta[v]+=low;
    }
    for(int i=1;i<=n;i++)
        if(delta[i]>0)Add(ss,i,delta[i]),sum+=delta[i];
        else Add(i,tt,-delta[i]);
    Add(t,s,Inf);
    int flow=Dinic();
    if(flow<sum)puts("please go home to sleep");
    else{
        ss=s;tt=t;
        printf("%d\n",Dinic());
    }
    return 0;
}

有源汇上下界最小流

题面

还是先按照上述的方法,将\(t->s\)连边\(\infty\)保证流量平衡,然后跑出来一个可行流.

接着在残量网络上面跑出\(t->s\)的最大流.最后那可行流减去最大流就是答案.

/*
  mail: mleautomaton@foxmail.com
  author: MLEAutoMaton
  This Code is made by MLEAutoMaton
*/
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<iostream>
using namespace std;
#define ll long long
#define re register
#define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
inline int gi(){
    int f=1,sum=0;char ch=getchar();
    while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
    return f*sum;
}
const int N=1000010,Inf=1e9+10;
int n,m,s,t,cur[N],front[N],cnt,dep[N],ss,tt,delta[N],sum;queue<int>Q;
struct node{int to,nxt,w;}e[N<<1];
void Add(int u,int v,int w){
    e[cnt]=(node){v,front[u],w};front[u]=cnt++;
    e[cnt]=(node){u,front[v],0};front[v]=cnt++;
}
bool bfs(){
    Q.push(ss);memset(dep,0,sizeof(dep));dep[ss]=1;
    while(!Q.empty()){
        int u=Q.front();Q.pop();
        for(int i=front[u];~i;i=e[i].nxt){
            int v=e[i].to;
            if(!dep[v] && e[i].w){
                dep[v]=dep[u]+1;Q.push(v);
            }
        }
    }
    return dep[tt];
}
int dfs(int u,int flow){
    if(u==tt || !flow)return flow;
    for(int &i=cur[u];~i;i=e[i].nxt){
        int v=e[i].to;
        if(dep[v]==dep[u]+1 && e[i].w){
            int di=dfs(v,min(e[i].w,flow));
            if(di){
                e[i].w-=di;e[i^1].w+=di;return di;
            }
            else dep[v]=0;
        }
    }
    return 0;
}
int Dinic(){
    int flow=0;
    while(bfs()){
        for(int i=1;i<=max(tt,n);i++)cur[i]=front[i];
        while(int d=dfs(ss,Inf))flow+=d;
    }
    return flow;
}
void del(int u){
    for(int i=front[u];~i;i=e[i].nxt)
        e[i].w=e[i^1].w=0;
}
int main(){
#ifndef ONLINE_JUDGE
    freopen("in.in","r",stdin);
#endif
    n=gi();m=gi();s=gi();t=gi();memset(front,-1,sizeof(front));
    ss=n+1;tt=ss+1;
    for(int i=1;i<=m;i++){
        int u=gi(),v=gi(),low=gi(),upp=gi();
        Add(u,v,upp-low);
        delta[u]-=low;delta[v]+=low;
    }
    for(int i=1;i<=n;i++)
        if(delta[i]>0)Add(ss,i,delta[i]),sum+=delta[i];
        else Add(i,tt,-delta[i]);
    Add(t,s,Inf);
    int flow=Dinic();
    if(flow<sum)puts("please go home to sleep");
    else{
        flow=e[cnt-1].w;
        del(ss);del(tt);ss=t;tt=s;e[cnt-1].w=e[cnt-2].w=0;
        printf("%d\n",flow-Dinic());
    }
    return 0;
}

转载于:https://www.cnblogs.com/mleautomaton/p/11240916.html

araw94333
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
上下界网络流
努力奋斗才能梦想成真
03-13 2972
犹记得上一年的时候做过上下界网络流的题目，但是那个时候只会套网上的建图方法，对内部的原理一知半解，因此前几天碰到这个题的时候就有些凌乱了，于是，重新学习了下，深刻理解了，估计不会再忘记了，呵呵。 好了，言归正传。 上下界网络流的话这篇博客写的很不错的，介绍了怎么建图，怎么求解。然后这篇就是更加的数学化，理论化一点，不过也很好理解的额 嗯，然后就站在巨人的肩膀上，总结一下。 上下界网络流的问
网络流总结·小编
08-31
网络流算法分析，谢谢~ 里面有各种实现的算法和代码，另外 SAP 为主要算法； 另外网络流的上下界、费用流的讲解都有哦。
有上下界的网络流学习笔记
weixin_30824479的博客
01-08 129
//其实主要还是自己复习用 //假定读者能够熟练打dinic的板子 有上下界的网络流的核心是”调整”,我们通过一个初始的未必可行的流调整出一个可行流,还可以从可行的未必最大/最小的流调整出最大/最小流. 另一个常用技巧是有源汇的流和无源汇的流(循环流)的转换.除了无源汇可行流的求解,其他有源汇的上下界网络流都要用到这个技巧. 无源汇有上下界可行流(也就是循环流) 模型:一个网络,...
网络流教程(max_flow)
09-18
### 网络流教程(max_flow) #### 一、网络流概述 本教程旨在系统地介绍网络流（Network Flow）的基本概念、理论基础及其在计算机科学领域的应用。网络流问题是一类重要的图论问题，在很多实际场景中都有广泛的应用...
有流量上下界的最大流最小流算法实现
03-04
总结来说，有流量上下界的最大流最小流算法是在网络流问题中的一个重要工具，通过Ford-Fulkerson方法及其变种，如Edmonds-Karp算法，我们可以有效地解决这类问题。在实现过程中，关键在于正确处理流量的增加和路径的...
第六节课-分支定界法（）
最新发布
03-16
由于这类规划不能有效地被连续模型表示，因此大部分情况下无法直接使用传统的线性规划或网络流模型来进行处理。尽管如此，它们的重要性并没有因此减少，相反，这些模型在工程和管理领域中往往代表了关键的决策问题。...
周源-有上下界的网络流周源-有上下界的网络流
10-01
周源-有上下界的网络流周源-有上下界的网络流周源-有上下界的网络流周源-有上下界的网络流
上下界网络流 cf704D
rgndao的博客
04-05 166
cf704D 给棋盘上n个棋子染色，红黑两种，代价不同，有m个约束，要求某一行或某一列的红黑棋子数差不能超过某一数值。求最小代价的染色方案。 上下界网络流的思路就是分别减去最低流量，剩余的放在图中跑最大流，对于原有大于零下界的边就通过新建超级源点ss和汇点tt解决，即出点多出的连向汇点，入点缺少的从源点补充。合法要求：源点s发出的必须能跑满，否则原有的平衡条件就打破了，无可行流。 对于原本就...
网络流模型应用实例与总结
"这篇资料主要汇总了网络流模型在解决各类问题中的应用，包括最大流、最小割、有上下界网络流和最小费用流等，并列举了一些经典的编程竞赛题目作为实例进行讲解。" 1. 最大流是网络流问题的基础，用于寻找网络中从...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值