最大网络流Edmond-Karp算法

#include<iostream>

using namespace std;

const int MAXV=2000;   //最大顶点数量
const int MAXE=2000;   //最大边数量
const int INFINITY = 0x7fffffff;   //无穷大
int capacity[MAXV][MAXV]; //记录残流网络的容量
int flow[MAXV];  //标记从源点到当前节点实际还剩多少流量可用
int pre[MAXV];  //标记在这条路径上当前节点的前驱,同时标记该节点是否在队列中
int Queue[MAXV];   //求最短增广路算法需要用到的队列

int MaxFlow(int src, int des, int n);
int Edmonds_Karp(int src, int des, int n);

int main()
{
int  m, n, u, v;

printf("请输入顶点数量：");
scanf("%d", &n);
printf("\n请输入边数量：");
scanf("%d", &m);

for(int i=0; i<m; ++i)
{
cin >> u >> v;
cin >> capacity[u][v];
}

cout << MaxFlow(0, n-1, n) << endl;

system("pause");
return 0;
}

int MaxFlow(int src, int des, int n)
{
int sumFlow = 0; //存储 src到des的最大流
int inc = 0; //存储最短增广路的流量增量
int u, v;

while (inc=Edmonds_Karp(src, des, n) != -1)
{
//逆向调整增广路（即修改对应残流网络的容量）
v = des;
while (v != src)
{
u = pre[v];
capacity[u][v] -= inc;
capacity[v][u] += inc;
v = u;
}

sumFlow += inc;
}

return sumFlow;
}

int Edmonds_Karp(int src, int des, int n)//广度优先搜索寻找最短增广路
{
int u, front = 0, rear = 0; //清空队列

for(int i=0; i<n; ++i) //初始化列表
{
pre[i] = -1;
}
//源点加入队列
flow[src] = INFINITY;
pre[src] = src;
Queue[rear++] = src;

while (front < rear) //队列非空
{
u = Queue[front++];
for(int i=0; i<n; ++i) //寻找未访问过的邻接点，并计算最小增加流量
{
if (pre[i] == -1 && capacity[u][i])
{
pre[i] = u;
flow[i] = (flow[u] < capacity[u][i]) ? flow[u] : capacity[u][i];

if (i == des) //找到汇点，直接返回最小增加流量
{
return flow[des];
}

Queue[rear++] = i;
}
}
}

//判断残流网络中是否存在增广路径
return (pre[des] == -1) ? -1 : flow[des];
}