今天刚接触网络流,感觉还不错,刚刚学完用EK和KK求网络最大流,前面的道依然艰辛,不过我会好好努力,最为网络流初学者我总结今日所学。至于网络流的一些基础认识我就没时间多述了。先了解网络流的基本概念之后,网络最大流的求法就有如下的步骤了:
1.首先根据原网络来构造其残流网络:
残流网络的构造呢就很简单满足一下条件就可以了:当通过边(u,v)的流f(u,v)>0则在残流网络上加上边(u,v),其上面的流量为c-f,如果c-f=0则这条边不添加,并且在残流网络上添加(v,u)边其上面的流量大小f(v,u)=f(u,v)。如图:
图2就是图1网络的残流网络。(图1中(4,4)代表(f,c));
2.根据残留网络寻找增广路径,在增广路径中找到最小的f,然后将路径上所有的边(u,v)上流量f(u,v)-=min(f),它相反的边f(v,u)+=min(f),然后循环的找和更新知道找不到增广路径。还能加载在原网络上的最大流就是之前所找到的所有的增光路经上的min(f).s。
POJ 1273是个练手的号题目,对于初学网络流者来说的基础题!
#include<iostream>
#include<string.h>
using namespace std;
#define MAXN 205
#define INF 0x7fffffff
bool mark[MAXN];
int que[MAXN*10],map[MAXN][MAXN],pre[MAXN];
int n,m,mim;
__int64 total;
void Init()
{
memset(mark,0,sizeof(mark));
memset(que,0,sizeof(que));
for( int i=0; i<=m; i++)
pre[i]=-1;
}
bool Edmonds_Karp();
void Update();
int main()
{
int a,b,c;
while( scanf("%d%d",&n,&m)!=EOF){
memset(map,0,sizeof(map));
while( n--){
scanf("%d%d%d",&a,&b,&c);
map[a][b]+=c;
}
total=0;
while( Edmonds_Karp())
Update();
printf("%I64d\n",total);
}
return 0;
}
bool Edmonds_Karp()
{
int first,rear,u,v;
first=rear=0;
Init();
mark[1]=true;
que[rear++]=1;
while( first!=rear){
int u=que[first++];
for( v=1; v<=m; v++){
if( !mark[v]&&map[u][v]>0){
que[rear++]=v;
mark[v]=true;
pre[v]=u;
if( v==m) return true;
}
}
}
return false;
}
void Update()
{
mim=INF;
int u,v;
v=m;
while( v!=1&&v!=-1){
u=pre[v];
if( map[u][v]<mim)
mim=map[u][v];
v=pre[v];
}
total+=mim;
v=m;
while( v!=-1&&v!=1){
u=pre[v];
map[u][v]-=mim;
map[v][u]+=mim;
v=pre[v];
}
}