网络最大流ISAP算法

模板

int ISAP(int s,int t,int n){
    for(int i=0;i<=n;i++) cur[i]=h[i],d[i]=INF;
    queue<int>q;
    q.push(t);
    d[t]=0;
    int flow=0;
    while(!q.emtpy()){
        int u=q.top();q.pop();
        for(int i=h[u];~i;i=e[i].xt)
            if(d[e[i].t]==INF) num[d[e[i].t]=d[u]+1]++,q.push(e[i].t);
    }
    for(int x=s,ok=0;d[s]<n;ok=0){
        if(x==t){
            int a=INF;
            for(;x!=s;x=e[p[x]].f) 
                a=min(a,e[p[x]].capa-e[p[x]].flow);
            for(x=t,flow+=a;x!=s;x=e[p[x]].f) 
                e[p[x]].flow+=a,e[p[x]^1].flow-=a;
        }
        for(int i=cur[x];~i;i=e[i].xt)
            if(e[i].capa>e[i].flow&&d[e[i].t]==d[u]-1){
                ok=1;
                cur[x]=i;
                p[x=e[i].t]=i;
                break;
            }
        if(!ok){
            int m=n-1;
            for(int i=cur[x]=h[x];~i;i=e[i].xt)
                if(e[i].capa>e[i].flow) m=min(m,d[e[i].t]);
            if(--num[d[x]]==0) break;
            num[d[x]=m+1]++
            if(x!=s) x=e[p[x]].f;
        }
    }
    return flow;
}

转载于:https://www.cnblogs.com/lsq647vsejgfb/p/9480417.html