关键概念与性质:
距离函数(distance function),我们说一个距离函数是有效的当且仅当满足有效条件(valid function)
(1)d(t)= 0;
(2)d(i)<= d(j)+1(如果弧rij在残余网络G(x)中);
性质1:
如果距离标号是有效的,那么d(i)便是残余网络中从点i到汇点的距离的下限;
证明:
性质2:
允许弧(边):如果对于边rij>0,且d(i)= d(j)+1,那么称为允许边
允许路:一条从源点到汇点的且有允许弧组成的路径
允许路是从源点到汇点的最短增广路径。
证明:
(1)因为rij>0所以必然是一条增广路
(2)假设路径p是一条允许路包含k条弧,那么由d(i) = d(j)+1 可知,d(s)= k;
又因为d(s)是点s到汇点的距离下限,所以距离下限为k,所以p便是一条最短路。
性质3:
在sap算法中距离标号始终是正确,有效的。并且每次的冲标号都会是距离标号严格递增
证明:略;
伪代码:
代码模板:
运用gap优化:
即当标号中出现了不连续标号的情况时,即可以证明已经不存在新的增广流,此时的流量即为最大流。
简单证明下:
假设不存在标号为k的结点,那么这时候可以将所有的结点分成两部分,一部分为d(i)>k,另一部分为d(i)<k
如此分成两份,因为性质2可知,允许路为最短的增广路,又因为不存在从>k部分到<k部分的增广流,那么有最
大流最小割定理可知此时便是最大流。
优化代码:要注意在标号的时候不能直接把所有的初始为0,而应该为-1,否则会在gap优化的时候出现问题,不满足递增的性质,切记!
注意:在运用sap的时候必须要时刻的保证标号的两个性质,因此,不能再初始标号的时候将全部初始为0层,因为有些点是不具有层数的,或者说是层数是无穷大的,不可达的。
网上找的比较清晰地模板sap,有各种优化