网络流的三个性质:
1、容量限制: f[u,v]<=c[u,v]
2、反对称性:f[u,v] = - f[v,u]
3、流量平衡: 对于不是源点也不是汇点的任意结点,流入该结点的流量和等于流出该结点的流量和。
只要满足这三个性质,就是一个合法的网络流.
最大流问题,就是求在满足网络流性质的情况下,源点 s 到汇点 t 的最大流量。
1、容量限制: f[u,v]<=c[u,v]
2、反对称性:f[u,v] = - f[v,u]
3、流量平衡: 对于不是源点也不是汇点的任意结点,流入该结点的流量和等于流出该结点的流量和。
只要满足这三个性质,就是一个合法的网络流.
最大流问题,就是求在满足网络流性质的情况下,源点 s 到汇点 t 的最大流量。
算法的关键在于
1)如何找出增广路径。
2)如何更新流量
1)如何找出增广路径。
2)如何更新流量
何为增广路径,说的直白些,所谓增广路径,就是找到这样一条路径,其流量不满,未达到容量上限。所有的可能的增广路径在一起便构成了残量网络。
残量网络:为了更方便算法的实现,一般根据原网络定义一个残量网络。其中r(u,v)为残量网络的容量。r(u,v) = c(u,v) – f(u,v),通俗地讲:就是对于某一条边(也称弧),还能再有多少流量经过。
增广路算法:每次用BFS找一条最短的增广路径,然后沿着这条路径修改流量值(实际修改的是残量网络的边权)。当没有增广路时,算法停止,此时的流就是最大流。
图解最大流:
可以证明,可行流为最大流,当且仅当不存在新的增广路径。
下面是代码实现:
#include<cstdio>
#include<vector>
#include<queue>
#include<cstring>
#include<iostream>
using namespace std;
#define INF 0x3f3f3f
#define maxn 100000
struct Edge
{
int from,to,cap,flow;
Edge(int u,int v,int c,int f):from(u),to(v),cap(c),flow(f) {}
};
struct EdmondsKarp
{
int m;
vector<Edge> edges;///边数的两倍
vector<int> G[maxn];///邻接表,G[i][j]表示结点i的第j条边在e数组中的序号
int a[maxn];///起点到i的可改进量
int p[maxn];///最短路树上p的入弧编号
void init(int n)
{
for(int i=0; i<n; i++)
G[i].clear();
edges.clear();
}
void AddEgde(int from,int to,int cap)
{
edges.push_back(Edge(from,to,cap,0));
edges.push_back(Edge(to,from,0,0));///反向弧
m=edges.size();
G[from].push_back(m-2);///记录每条弧的编号
G[to].push_back(m-1);///记录对应反向弧的编号
}
int Maxflow(int s,int t)
{
int flow=0;
for(;;)
{
memset(a,0,sizeof(a));
queue<int>Q;
a[s]=INF;
Q.push(s);
while(!Q.empty())///寻找増广路
{
int x=Q.front();
Q.pop();
for(int i=0; i<G[x].size(); i++)
{
Edge& e=edges[G[x][i]];
if(!a[e.to]&&e.cap>e.flow)
{
p[e.to]=G[x][i];
a[e.to]=min(a[x],e.cap-e.flow);
Q.push(e.to);
}
}
if(a[t])break;
}
if(!a[t])break;
for(int u=t; u!=s; u=edges[p[u]].from)///沿着増广路更新流
{
edges[p[u]].flow+=a[t];
edges[p[u]^1].flow-=a[t];
}
flow+=a[t];
}
return flow;
}
} w;
int main()
{
int n,i,from,to,cap,s,t;
while(~scanf("%d%d",&n,&t))///n为边数t为结点数
{
w.init(n);
for(i=0; i<n; i++)
{
scanf("%d%d%d",&from,&to,&cap);
w.AddEgde(from,to,cap);
}
cout<<w.Maxflow(1,t)<<endl;
}
}
/*
输入
8 6
1 2 6
2 3 4
1 4 7
4 5 2
4 3 3
2 5 4
5 6 8
3 6 4
输出
10
*/