网络流-SAP算法模板

node 表示当前讨论到的节点； flow 表示该节点被传入/要传出的流量；s, t, ver 分别是起点，终点，节点总数； SAP 函数返回增广路成功增大的流量； dlt 是当前已传出的流量； c 是残量。

const int INF = ~0U >> 1;

int SAP(int node, int flow)
{
    if (node == t)
        return flow;
    int p, tmp, dlt = 0;
    for (p = fst[node]; p; p = nxt[p]) {
        if (w[p] <= 0 || dis[v[p]] + 1 != dis[node])
            continue;
        tmp = dfs(v[p], min(flow - dlt, c[p]));
        c[p] -= tmp;
        c[p + 1] += tmp;
        if ((dlt += tmp) == flow || dis[s] == ver)
            return dlt;
    }
    if (!--cnt[dis[node]])
        dis[s] = ver;
    ++cnt[++dis[node]];
    return dlt;
}

int main()
{
    /* input */
    while (dis[s] < ver)
        ans += SAP(s, INF);
    printf("%d\n", ans);
    return 0;
}