网络流最大流ISAP算法[草地排水][运输问题1]

1、距离标号：

所谓距离标号 ，就是某个点到汇点的最少的弧的数量（即边权值为1时某个点到汇点的最短路径长度）。

设点i的标号为level[i]，那么如果将满足level[i]=level[j]+1的弧(i,j)叫做允许弧 ，且增广时只走允许弧。

2、断层（本算法的Gap优化思想）：

gap[i]数组表示距离标号为i的点有多少个，如果到某一点没有符合距离标号的允许弧，那么需要修改距离标号来找到增广路；

如果重标号使得gap数组中原标号数目变为0，则算法结束。

SAP算法框架：

1、初始化；

2、不断沿着可行弧找增广路。可行弧的定义为{( i , j ) , level[i]==level[j]+1};

3、当前节点遍历完以后，为了保证下次再来的时候有路可走，重新标号当前距离，level[i]=min(level[j]+1);

该算法最重要的就是gap常数优化了。

转自：http://www.cnblogs.com/wally/archive/2013/05/03/3054778.html

模板：

描述

在农夫约翰的农场上，每逢下雨，贝茜最喜欢的三叶草地就积聚了一潭水。这意味着草地被水淹没了，并且小草要继续生长还要花相当长一段时间。因此，农夫约翰修建了一套排水系统来使贝茜的草地免除被大水淹没的烦恼（不用担心，雨水会流向附近的一条小溪）。作为一名一流的技师，农夫约翰已经在每条排水沟的一端安上了控制器，这样他可以控制流入排水沟的水流量。 

农夫约翰知道每一条排水沟每分钟可以流过的水量，和排水系统的准确布局（起点为水潭而终点为小溪的一张网）。需要注意的是，有些时候从一处到另一处不只有一条排水沟。

根据这些信息，计算从水潭排水到小溪的最大流量。对于给出的每条排水沟，雨水只能沿着一个方向流动，注意可能会出现雨水环形流动的情形。

格式

PROGRAM NAME: ditch

INPUT FORMAT:

(file ditch.in)

第1行:	两个用空格分开的整数N (0 <= N <= 200)和M (2 <= M <= 200)。N是农夫约翰已经挖好的排水沟的数量，M是排水沟交叉点的数量。交点1是水潭，交点M是小溪。
第二行到第N+1行:	每行有三个整数，Si, Ei,和Ci。Si和Ei (1 <= Si, Ei <= M)指明排水沟两端的交点，雨水从Si流向Ei。Ci (0 <= Ci <= 10,000,000)是这条排水沟的最大容量。

 

OUTPUT FORMAT:

(file ditch.out)

输出一个整数，即排水的最大流量。

 

SAMPLE INPUT

5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10

SAMPLE OUTPUT

50 

 

#include<iostream>
#include<cstdio>
#include<cstring>
#define MM 1000000000
int n,m,map[201][201],level[201]={0},gap[201],pre[201],l[100000];
bool f[201];
using namespace std;
int ISAP(int x,int y){	
    int v,minn,maxn,u,ans=0,i;
	memset(gap,0,sizeof(gap));
	memset(pre,-1,sizeof(pre));
	pre[x]=x;u=x;gap[0]=m;
	while(level[x]<m){
		for(v=1;v<=m;++v){
			if(map[u][v]>0&&level[v]+1==level[u]){
				break;
			}
		}
		if(v<=m){
			pre[v]=u;u=v;
			if(v==m){
				maxn=MM;
				for(i=v;i!=x;i=pre[i]){
					maxn=min(maxn,map[pre[i]][i]);
				}
				ans+=maxn;
				for(i=v;i!=x;i=pre[i]){
					map[pre[i]][i]-=maxn;
					map[i][pre[i]]+=maxn;
				}
				u=x;
			}
		}
		else{
		    minn=m;
			for(v=1;v<=m;++v){
				if(map[u][v]>0){
					minn=min(minn,level[v]);
				}
			}	
			minn+=1;
			gap[level[u]]-=1;
			if(!gap[level[u]]) break;
			level[u]=minn;
			gap[level[u]]+=1;
			u=pre[u];
		}
	}
	return ans;
}
int main()
{
	freopen("ditch.in","r",stdin);
	freopen("ditch.out","w",stdout);
	int i,j,x,y,z,h,t;
	memset(map,0,sizeof(map));
	memset(f,1,sizeof(f));
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;++i){
		scanf("%d%d%d",&x,&y,&z);
		map[x][y]+=z;
	}
/*	l[1]=m;f[m]=false;level[m]=0;h=t=1;
	while(h<=t){
		for(i=1;i<=m;++i){
			if(map[i][l[h]]>0&&f[i]){
				t+=1;
				l[t]=i;
				level[i]=level[l[h]]+1;
				f[i]=false;
			}
		}
		h+=1;
	}*/
	printf("%d\n",ISAP(1,m));
}

提示:

        保存距离的level数组的反向广搜不需要在初始化的时候做，在做ISAP的时候就已经动态的维护了，我开始很逗比的初始化了结果WA了两遍！！

11. 运输问题1

【问题描述】

    一个工厂每天生产若干商品，需运输到销售部门进行销售。从产地到销地要经过某些城镇，有不同的路线可以行走，每条两城镇间的公路都有一定的流量限制。请你计算，在不考虑其它车辆使用公路的前提下，如何充分利用所有的公路，使产地运输到销地的商品最多，最多能运输多少商品。

【输入格式】

输入文件有若干行
第一行，一个整数n，表示共有n个城市(2<=n<=100),产地是1号城市，销地是n号城市。
下面有n行,每行有n个数字。第p行第q列的数字表示城镇p与城镇q之间有无公路连接。数字为0表示无，大于0表示有公路，且该数字表示该公路流量。

【输出格式】

输出文件有一行
第一行，1个整数max，表示最大流量为max。

【输入输出样例】

输入文件名： maxflowa.in

6
0 4 8 0 0 0
0 0 4 4 1 0
0 0 0 2 2 0
0 0 0 0 0 7
0 0 0 6 0 9
0 0 0 0 0 0

输出文件名：maxflowa.out

8

 

这次用next数组敲了一遍模板：

#include<iostream>
#include<cstdio>
using namespace std;
struct use{
	int st,en,va;
}edge[20001];
int point[101]={0},next[20001]={0},pre[101]={0},dis[101]={0},gap[101]={0};
int ISAP(int stt,int enn)
{
	int i,j,x,y,u,ans=0;
	bool f=false;
	gap[0]=enn;u=stt;
	while(dis[stt]<enn){
		f=false;y=point[u];     
		while (y!=0){
			if (dis[u]==dis[edge[y].en]+1&&edge[y].va>0){
				f=true;break;
			}
			y=next[y];
		}
		if (f){
			pre[edge[y].en]=y;u=edge[y].en;y=2100000000;
			if (u==enn){
				for (i=enn;i!=stt;i=edge[pre[i]].st)
					y=min(y,edge[pre[i]].va);
				ans+=y;
				for (i=enn;i!=stt;i=edge[pre[i]].st){
					  x=(pre[i]-1)^1+1;
					  edge[pre[i]].va-=y;
					  edge[x].va+=y;
				}
				u=stt;
			}
		}
		else{
			--gap[dis[u]];
			if (gap[dis[u]]==0) return ans;
			y=enn;
			for (i=point[u];i!=0;i=next[i])
			  if (edge[i].va>0)
				y=min(y,dis[edge[i].en]);
			dis[u]=y+1;
			++gap[dis[u]];
			if (u!=stt) u=edge[pre[u]].st;
	    }
	}
	return ans;
}
int main()
{
	freopen("maxflowa.in","r",stdin);
	freopen("maxflowa.out","w",stdout);
	int n,i,j,tot=0,x;
	scanf("%d",&n);
	for (i=1;i<=n;++i)
	  for (j=1;j<=n;++j){
	  	 scanf("%d",&x);
	  	 if (x){
	  	   ++tot;next[tot]=point[i];point[i]=tot;
	  	   edge[tot].st=i;edge[tot].en=j;edge[tot].va=x;
	  	   ++tot;next[tot]=point[j];point[j]=tot;
	  	   edge[tot].st=j;edge[tot].en=i;edge[tot].va=0;
	     }
	  }
	printf("%d\n",ISAP(1,n));
}