最大流模板

最新推荐文章于 2019-10-07 09:42:23 发布

weixin_30314813

最新推荐文章于 2019-10-07 09:42:23 发布

阅读量119

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/ACAC/archive/2010/05/18/1738717.html

博客给出了最大流算法的代码实现。定义了相关常量和二维数组表示网络，通过BFS搜索增广路径，不断更新网络中的流，直到不存在增广路径，最终返回最大流结果。代码包含队列操作、路径记录等内容。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <iostream>
using namespace std;
#define MAXN 205
#define INF 2110000000
#define MIN(x,y) (x<y?x:y)
int map[MAXN][MAXN];

int max_flow(int num,int map[][MAXN],int source,int sink)//参数含义：结点数量网络源点汇点
{
    int my_queue[MAXN],queue_first,queue_end;//数组做队列实现BFS搜索路径
    int pre[MAXN],min_flow[MAXN];//记录结点的父节点当前路径中最小的一段的值，也即限制值
    int flow[MAXN][MAXN];//记录当前网络中的流
    int ans=0;//最终结果
    memset(flow,0,sizeof(flow));
    while(1)//一直循环，直到不存在增广路径
    {
        queue_first=0;//初始化队列
        queue_end=0;
        my_queue[queue_end++]=source;
        memset(pre,-1,sizeof(pre));
        min_flow[source]=INF;
        pre[source]=-2;//源点的父节点需特殊标示
        while(queue_first<queue_end)//BFS寻找增广路径
        {
            int temp=my_queue[queue_first++];//出队列
            for(int i=0;i<num;i++)//由结点temp往外扩展
            {
                if(pre[i]==-1&&flow[temp][i]<map[temp][i])//当结点i还未被探索到，并且还有可用流量
                {
                    my_queue[queue_end++]=i;//加入队列
                    pre[i]=temp;//标示父节点
                    min_flow[i]=MIN(min_flow[temp],(map[temp][i]-flow[temp][i]));//求得min_flow
                }
            }
            if(pre[sink]!=-1)//sink的父节点不为初始值，说明BFS已经找到了一条路径
            {
                int k=sink;
                while(pre[k]>=0)
                {
                    flow[pre[k]][k]+=min_flow[sink];//将新的流量加入flow
                    flow[k][pre[k]]=-flow[pre[k]][k];
                    k=pre[k];
                }
                break;
            }
        }
        if(pre[sink]==-1) return ans;//不存在增广路径，返回
        else ans+=min_flow[sink];
    }
}
本文来自优快云博客，转载请标明出处：http://blog.youkuaiyun.com/logic_nut/archive/2009/08/03/4403751.aspx