最大流SAP算法

       距离标记最短增广路算法，复杂度是O(n^2*m),比EK要快得多，看了书上的原理，然后理解了伪代码，一下子就写过了，1459跑了130ms， 还需要改进啊， 挺高兴的终于学了网络流快一点的算法了，半期后把dinic和预留推进算法学了。

       其实大致思路就是每次增广最短路径， 维护一个dis数组表示每个点到汇点的最小跳数，dis数组必须有这个性质： dis[i] <= dis[j] + 1 ( r[i][j] > 0); 定义一条允许弧：如果满足：dis[i] = dis[j] + 1( r[i][j] > 0);那么连接i，j的边就叫允许弧， 书上有这个定理：从源点到汇点的最短路一定是用允许弧构成。所以每次扩展路径都找允许弧，如果i没有允许弧就更新dis[i] = min{ dis[j] + 1 | r[i][j] > 0); 如此容易知道当dis[source] = n是就增广完毕，因为从源点到汇点最多有n-1步。最后dis数组可以用一个BFS初始化。

#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； //c1是c的反向网络，用于dis数组的初始化 int dis[Max]; void initialize()// 初始化dis数组 { int q[Max],head = 0, tail = 0;//BFS队列 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]; // top是当前增广路中最前面一个点。 memset(low,0,sizeof(low));//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];//更新low 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++)//找和top相邻dis最小的点 { if(r[top][j] > 0 && mindis > dis[j] + 1) mindis = dis[j] + 1; } dis[top] = mindis;//更新top的距离值 if(top != source) top = pre[top];// 回溯找另外的路 } } return(flow); }