最大流 Dinic + Sap 模板

不说别的,直接上模板。

 

Dinic+当前弧优化:

struct Edge{
    int x,y,c,ne;
}e[M*2];
int be[N],all;
int d[N],q[N];
int stack[N],top;//栈存的是边
int cur[N];//当前弧优化
void add(int x, int y, int z)//需保证相反边第一个为偶数 { e[all].x=x; e[all].y=y; e[all].c=z; e[all].ne=be[x]; be[x]=all++; e[all].x=y; e[all].y=x; e[all].c=0; e[all].ne=be[y]; be[y]=all++; }
bool BFS(int s, int t)//化为层次图 使得边数从m降低为n 复杂度随之下降 { memset(d,-1,sizeof(d)); int head=0,tail=0; q[++tail]=s; d[s]=0; while(head!=tail) { int u=q[++head]; for(int i=be[u]; i!=-1; i=e[i].ne) if(e[i].c>0 && d[e[i].y]==-1){ d[e[i].y]=d[u]+1; q[++tail]=e[i].y; if(tail==N-1) tail=0; if(e[i].y==t) return 1; } } return 0; } int Dinic(int s, int t)//防止爆栈 用stack模拟递归 { int ans=0; while(BFS(s,t)) { memcpy(cur,be,sizeof(be)); int u=s; top=0;//dfs开始 清空栈 while(1) { if(u==t) { int minc=1000000000,mini; for(int i=0; i<top; i++) if(minc>e[stack[i]].c) { minc=e[stack[i]].c; mini=i;//以便之后回到这继续增广 } for(int i=0; i<top; i++) { e[stack[i]].c-=minc; e[stack[i]^1].c+=minc;//第一个二进制取反 即取相反边 } ans+=minc; top=mini; u=e[stack[mini]].x; } for(int i=cur[u]; i!=-1; cur[u]=i=e[cur[u]].ne) if(e[i].c>0 && d[e[i].y]==d[e[i].x]+1) break; if(cur[u]!=-1) { stack[top++]=cur[u]; u=e[cur[u]].y; }else { if(top==0) break; //循环结束标志 d[u]=-1;//当前节点不在增广路中 删除 u=e[stack[--top]].x;//回溯 } } } return ans; }

ISAP+GAP+当前弧优化:

struct Edge{
    int x,y,c,ne;
}e[M*2];
int x,y,z,n,m,s,t;
int be[N],all;
int d[N],q[N];
int stack[N];//模拟递归
int gap[N],cur[N];//gap优化+当前弧优化

void add(int x, int y, int z)//保证第一个为偶数
{
    e[all].x=x; e[all].y=y; e[all].c=z;
    e[all].ne=be[x];
    be[x]=all++;
    e[all].x=y; e[all].y=x; e[all].c=0;
    e[all].ne=be[y];
    be[y]=all++;
}

void BFS(int s, int t)
{
    memset(d,-1,sizeof(d));
    memset(gap,0,sizeof(gap));
    gap[0]=1;
    int head=0,tail=0;
    q[++tail]=t;
    d[t]=0;
    while(head!=tail)
    {
        int u=q[++head];
        for(int i=be[u]; i!=-1; i=e[i].ne)
            if(d[e[i].y]==-1)
            {
                d[e[i].y]=d[u]+1;
                q[++tail]=e[i].y;
                gap[d[e[i].y]]++;
            }
    }
}

int sap(int s, int t, int n)
{
    int ans=0;
    BFS(s,t);
    memcpy(cur,be,sizeof(be));
    int top=0;
    int u=s;
    while(d[s]<n)
    {
        if(u==t)
        {
            int minc=1000000000,mini;
            for(int i=0; i<top; i++)
                if(minc>e[stack[i]].c)
                {
                    minc=e[stack[i]].c;
                    mini=i;
                }
            for(int i=0; i<top; i++)
            {
                e[stack[i]].c-=minc;
                e[stack[i]^1].c+=minc;
            }
            ans+=minc;
            top=mini;
            u=e[stack[mini]].x;
            continue;
        }
        for(int i=cur[u]; i!=-1; cur[u]=i=e[i].ne)//当前弧优化
            if(e[i].c>0 && d[e[i].y]+1==d[u]) break;
        if(cur[u]!=-1)
        {
            stack[top++]=cur[u];
            u=e[cur[u]].y;
        }else
        {
            int mind=n;
            for(int i=be[u]; i!=-1; i=e[i].ne)//更新距离标号
                if(e[i].c>0 && mind>d[e[i].y])
                {
                    mind=d[e[i].y];
                    cur[u]=i;
                }
            gap[d[u]]--;
            if(!gap[d[u]]) return ans;//gap表示当前距离的点有多少个 一旦==0 说明断层直接退出循环
            d[u]=mind+1;
            gap[d[u]]++;
            if(u!=s) u=e[stack[--top]].x;
        }
    }
    return ans;
}

void init()
{
    all=0;
    memset(be,-1,sizeof(be));
}

 

posted @ 2014-08-06 00:02  Estimator  阅读(383)  评论(0编辑  收藏  举报