最大流(邻接阵形式)-- 求网络最大流,邻接阵形式

最新推荐文章于 2021-07-16 15:23:07 发布

weixin_30757793

最新推荐文章于 2021-07-16 15:23:07 发布

阅读量73

点赞数

原文链接：http://www.cnblogs.com/agpro/archive/2010/07/05/1771529.html

版权

//最大流(邻接阵形式)-- 求网络最大流,邻接阵形式
//求网络最大流,邻接阵形式
//返回最大流量,flow返回每条边的流量
//传入网络节点数n,容量mat,源点source,汇点sink

const int MAXN = 302;
const int inf = 1000000000;

int max_flow(int n,int mat[][MAXN],int source,int sink,int flow[][MAXN])
{
	int pre[MAXN],que[MAXN],d[MAXN],p,q,t,i,j;
	if (source==sink) 
		return inf;
	for (i=0;i<n;i++)
		for (j=0;j<n;flow[i][j++]=0);
	for (;;)
	{
		for (i=0;i<n;pre[i++]=0);
		pre[t=source]=source+1,d[t]=inf;
		for (p=q=0;p<=q&&!pre[sink];t=que[p++])
			for (i=0;i<n;i++)
				if (!pre[i]&&(j=mat[t][i]-flow[t][i]))
					pre[que[q++]=i]=t+1,d[i]=d[t]<j?d[t]:j;
				else if (!pre[i]&&(j=flow[i][t]))
					pre[que[q++]=i]=-t-1,d[i]=d[t]<j?d[t]:j;
		if (!pre[sink]) 
			break;
		for (i=sink;i!=source;)
			if (pre[i]>0)
				flow[pre[i]-1][i]+=d[sink],i=pre[i]-1;
			else
				flow[i][-pre[i]-1]-=d[sink],i=-pre[i]-1;
	}
	for (j=i=0;i<n;j+=flow[source][i++]);
	return j;
}