网络流之SAP算法

在SAP算法中，我们定义每个顶点的距离标号(Distance Labels)，即残留网络中这个点到汇点的距离。然后，我们只在距离标号相邻的点间寻找增广路径。如果从一个点出发没有容许边，就对重标记并回溯（类似预流推进）。首先介绍一些定义：

残留容量：容量－流量。即，残留容量r[i,j]=c[i,j]-f[i,j]。如果只求流值，则在网络流算法中只需记录残量，而不必记录容量和流量。

残留边、残留网络：有残留容量的边成为残留边，由网络的点集、残留边集构成残留网络。

距离标号 ：对于每个点i，定义其在当前残留网络中到汇点的距离为d[i]。如d[t]=0。

容许边：当且仅当r[i,j]>0且d[i]=d[j]+1时，边(i,j)是容许边。我们的算法只考虑容许边。

算法的一开始，我们从汇点进行一次BFS求出距离标号。然后从源点开始递归操作。对于一个点i的操作，如果存在容许边(i,j)，当j是汇点时增广并从源点开始递归，否则就令i是j的前驱，并递归j。否则对i重标号，使得d[i]=min{d[j]}+1 ( r[i,j] )>0，并回溯。当d[s]>=n时算法停止，此时不存在任何一条增广路。

模板一

#include <iostream>

#include <cstring>

#include <cstdlib>

using namespace std;

const int Max = 225;

const int oo = 210000000;

int n,m,c[Max][Max],r[Max][Max],c1[Max][Max],source,sink；

int dis[Max];

void initialize()

{

     int q[Max],head = 0, tail = 0;

     memset(dis,0,sizeof(dis));

     q[++head] = sink;

     while(tail < head)

     {

         int u = q[++tail], v;

         for(v = 0; v <= sink; v++)

         {

               if(!dis[v] && c1[u][v] > 0)

               {

                   dis[v] = dis[u] + 1;

                   q[++head] = v;

               }

         }

     }

}

int maxflow_sap()

{

    initialize();

    int top = source, pre[Max], flow = 0, i, j, k, low[Max];

 

    memset(low,0,sizeof(low));

    while(dis[source] < n)

    {

        bool flag = false;

        low[source] = oo;

        for(i = 0; i <= sink; i++)

        {

            if(r[top][i] > 0 && dis[top] == dis[i] + 1)

            {

                flag = true;

                break;

            }

        }

        if(flag)

        {

            low[i] = r[top][i];

            if(low[i] > low[top]) low[i] = low[top];

            pre[i] = top; top = i;

            if(top == sink)

            {

                flow += low[sink];

                j = top;

                while(j != source)

                {

                    k = pre[j];

                    r[k][j] -= low[sink];

                    r[j][k] += low[sink];

                    j = k;

                }

                top = source;

                memset(low,0,sizeof(low));

            }

        }

        else

        {

            int mindis = 10000000;

            for(j = 0; j <= sink; j++)

            {

                if(r[top][j] > 0 && mindis > dis[j] + 1)

                    mindis = dis[j] + 1;

            }

            dis[top] = mindis;

            if(top != source) top = pre[top];

        }

    }

    return(flow);

}

模板二

#include <iostream>
#include <queue>
#define msize 1024      //最大顶点数目
using namespace std;
 
int d[msize];           //标号
int r[msize][msize];    //残留网络，初始为原图
int num[msize];         //num[i]表示标号为i的顶点数有多少
int pre[msize];         //记录前驱 
int n,m,s,t;            //m个顶点，n条边，从源点s到汇点t
 //此处别忘了在使用最大流之前先调用此函数初始化标号
void ini_d() //BFS计算标号，汇点t标号为0
{
    int k;
    queue<int>Q;
 
    memset(d,1,sizeof(d));          //将距离设置成为无穷大,此处亦可以使用for循环实现 
    memset(num,0,sizeof(num));
 
    Q.push(t);
    d[t]=0;                   //汇点的标号为0 
    num[0]=1;
    while (!Q.empty())
    {
        k=Q.front(),Q.pop();
        for (int i=0;i<m;i++)       //遍历所有的结点 
        {
            if (d[i]>=m&&r[i][k]>0) //此处要特别注意，通过frontint的值改变其他的距离标号
            {
                d[i]=d[k]+1;
                Q.push(i);
                num[d[i]]++;
            }
        }
    }
}
 
int findAlowArc(int i)       //从i出发寻找允许弧
{
    int j;
    for (j=0;j<m;j++) if (r[i][j]>0&&d[i]==d[j]+1) return j;
 
    return -1;
}
 
int reLable(int i)         //重新标号
{
    int mm=INT_MAX;
    for (int j=0;j<m;j++)
        if (r[i][j]>0) mm=min(mm,d[j]+1);
 
    return mm==INT_MAX?m:mm;
}
 
int maxFlow(int s,int t)      //从源点s出发的最大流
{
    int flow=0,i=s,j;
    int delta;              //增量
 
    memset(pre,-1,sizeof(pre));
    while (d[s]<m)
    {
        j=findAlowArc(i);
        if (j>=0)
        {
            pre[j]=i;                     
            i=j;                //从前往后找 
            if (i==t)           //更新残留网络
            {
                delta=INT_MAX;
                for (i=t;i!=s;i=pre[i]) delta=min(delta,r[pre[i]][i]); //找到增广路径的增量 
                for (i=t;i!=s;i=pre[i]) r[pre[i]][i] -= delta, r[i][pre[i]] += delta; //更改流量 
                flow += delta;
            }
        }
        else
        {
            int x=reLable(i);       //重新标号
            num[x]++;
            num[d[i]]--;
            if (num[d[i]]==0) return flow;      //间隙优化
            d[i]=x;
            if (i!=s) i=pre[i];
        }
    } 
    return flow;
}