SAP求最大流

学最大流,在网上和书上看了别人的代码,感觉不太好理解,有些乱。整理了一下,形成了自己的模板。


#define INF 0xffffff
#define MAX_NODE 205
#define MAX_EDGE 205


struct node{
    int next,c,to;
    node(){}
    node(int c,int to,int next):c(c),to(to),next(next){}
}edge[MAX_EDGE];

int head[MAX_NODE],dis[MAX_NODE],gap[MAX_NODE],cur[MAX_NODE],pre[MAX_NODE];
int top;
int n,m,START,END;

void add(int a,int b,int c){
    edge[top]=node(c,b,head[a]);
    head[a]=top++;
    edge[top]=node(0,a,head[b]);
    head[b]=top++;
}


int SAP(){
    int max_flow=0,flow=INF,u=START,v;
    for(int i=1;i<=n;++i)dis[i]=gap[i]=0,cur[i]=head[i];
    gap[0]=n;pre[START]=START;
    while(dis[START]<n){
        for(int &i=cur[u];i!=-1;i=edge[i].next){
                v=edge[i].to;
            if(edge[i].c&&dis[u]==dis[v]+1){
                break;
            }
        }


        if(cur[u]!=-1){
            flow=min(flow,edge[cur[u]].c);
            pre[v]=u;u=v;
            if(v==END){
                max_flow+=flow;
                int i=u;
                while(i!=START){
                    i=pre[i];
                    edge[cur[i]].c-=flow;
                    edge[cur[i]^1].c+=flow;
                }
                flow=INF;
                u=START;
            }
        }

        else{
            int minds=n;
            for(int i=head[u];i!=-1;i=edge[i].next){
                    v=edge[i].to;
                if(edge[i].c&&minds>dis[v]){
                    minds=dis[v];
                    cur[u]=i;
                }
            }
             if(--gap[dis[u]]==0)break;
            ++gap[dis[u]=minds+1];
            u=pre[u];
        }
    }
    return max_flow;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值