最大流(Dinic模版)——HDOJ 1532

HDOJ 1532 Drainage Ditches

这道题目是用来做Dinic算法的模版题目，而上一篇博客 最大流——HDOJ 3549 是用来作为EK算法的模版题目。

关于Dinic算法的讲解，网上有很多，我这里也不多讲什么。推荐一篇讲得比较容易懂的 http://comzyh.tk/blog/archives/568/#Dinic-Code

/*
HDOJ 1532
最大流经典入门题目 下面是Dinic算法，矩阵存储图，易于理解
*/

#include <iostream>
using namespace std;

#define min(a,b) (a)<(b)?a:b
#define MAX 0x5fffffff

int map[250][250];	//邻接矩阵 
int dis[250];	//距源点距离,分层图 
int que[2000],front,rear;	//BFS队列 ,首,尾 
int N,M,ANS;	//N:点数;M,边数

int BFS()
{
	int i,j;
	memset(dis,0xff,sizeof(dis));//以-1填充 
	dis[1]=0;
	front=rear=0;
	que[0]=1;
	while (front <= rear)
	{
		j=que[front++];
		for (i=1;i<=N;i++)
			if (dis[i]<0 && map[j][i]>0)
			{
				dis[i]=dis[j]+1; 
				que[++rear]=i;
			}
	}
	if (dis[N]>0)
        return 1;
	else
        return 0;//汇点的DIS小于零,表明BFS不到汇点 
}

//Find代表一次增广,函数返回本次增广的流量,返回0表示无法增广 
int find(int x,int low)//Low是源点到现在最窄的(剩余流量最小)的边的剩余流量
{
    int i,a=0;
    if (x==N)
		return low;//是汇点 
    for(i=1;i<=N;i++)
		if (map[x][i] >0 //联通 
			&& dis[i]==dis[x]+1 //是分层图的下一层 
			&&(a=find(i,min(low,map[x][i]))))//能到汇点(a <> 0) 
		{
			map[x][i] -= a;
			map[i][x] += a;
			return a;
		}
		return 0;
		
}

int main()
{
    int i,j,f,t,flow,tans;
    while (scanf("%d%d",&M,&N)!=EOF)
	{
		memset(map,0,sizeof(map));
		for(i=1;i<=M;i++)
		{
			scanf("%d%d%d",&f,&t,&flow);
			map[f][t] += flow;
		}		
		ANS=0;
		while(BFS())//要不停地建立分层图,如果BFS不到汇点才结束 
		{
			while(tans=find(1,0x7fffffff))
				ANS+=tans;//一次BFS要不停地找增广路,直到找不到为止 
		}
		printf("%d\n",ANS);
    }
    return 0;
}