您的位置:首页 > 理论基础 > 计算机网络

HihoCoder #1369 : 网络流一·Ford-Fulkerson算法

2018-02-13 21:23 453 查看

题目描述:

描述

小Hi和小Ho住在P市,P市是一个很大很大的城市,所以也面临着一个大城市都会遇到的问题:交通拥挤。

小Ho:每到周末回家感觉堵车都是一种煎熬啊。

小Hi:平时交通也还好,只是一到上下班的高峰期就会比较拥挤。

小Ho:要是能够限制一下车的数量就好了,不知道有没有办法可以知道交通系统的最大承受车流量,这样就可以限制到一个可以一直很顺畅的数量了。

小Hi:理论上是有算法的啦。早在1955年,T.E.哈里斯就提出在一个给定的网络上寻求两点间最大运输量的问题。并且由此产生了一个新的图论模型:网络流。

小Ho:那具体是啥?

小Hi:用数学的语言描述就是给定一个有向图G=(V,E),其中每一条边(u,v)均有一个非负数的容量值,记为c(u,v)≥0。同时在图中有两个特殊的顶点,源点S和汇点T。

举个例子:

其中节点1为源点S,节点6为汇点T。

我们要求从源点S到汇点T的最大可行流量,这个问题也被称为最大流问题。

在这个例子中最大流量为5,分别为:1→2→4→6,流量为1;1→3→4→6,流量为2;1→3→5→6,流量为2。

小Ho:看上去好像挺有意思的,你让我先想想。

提示:Ford-Fulkerson算法

输入

第1行:2个正整数N,M。2≤N≤500,1≤M≤20,000。

第2..M+1行:每行3个整数u,v,c(u,v),表示一条边(u,v)及其容量c(u,v)。1≤u,v≤N,0≤c(u,v)≤100。

给定的图中默认源点为1,汇点为N。可能有重复的边。

输出

第1行:1个整数,表示给定图G的最大流。

样例输入

6 7

1 2 3

1 3 5

2 4 1

3 4 2

3 5 3

4 6 4

5 6 2

样例输出

5

题解:

DINIC第二题,也是模板题,用来熟练DINIC写法。

代码如下:

#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
const int maxn=505,maxm=40005,inf=1<<30;
int n,m,tot,lnk[maxn],son[maxm],nxt[maxm],w[maxm],dep[maxn],que[maxn];
inline int read(){
int x=0; char ch=getchar();
while (ch<'0'||ch>'9') ch=getchar();
while (ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
return x;
}
void add(int x,int y,int z){
son[++tot]=y,w[tot]=z,nxt[tot]=lnk[x],lnk[x]=tot;
}
int dfs(int x,int sum){
if (x==n) return sum;
for (int j=lnk[x];j>-1;j=nxt[j])
if (dep[son[j]]==dep[x]+1&&w[j]>0) {
int sum1=dfs(son[j],min(sum,w[j]));
if (sum1>0) {w[j]-=sum1,w[j^1]+=sum1; return sum1;}
}
return 0;
}
int bfs(){
int head=0,tail=1; que[tail]=1;
memset(dep,0,sizeof(dep)); dep[1]=1;
while (head!=tail) {
head++;
for (int j=lnk[que[head]];j!=-1;j=nxt[j])
if (dep[son[j]]==0&&w[j]>0) que[++tail]=son[j],dep[son[j]]=dep[que[head]]+1;
}
return dep
;
}
int dinic(){
int ans=0;
while (bfs()) {
while (int sum=dfs(1,inf)) ans+=sum;
}
return ans;
}
int main(){
n=read(),m=read();
tot=-1; memset(lnk,255,sizeof(lnk)); memset(nxt,255,sizeof(nxt));
for (int i=1;i<=m;i++) {int x=read(),y=read(),z=read(); add(x,y,z); add(y,x,0);}
printf("%d\n",dinic());
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: