Control HDU - 4289

没什么难的,不难看出这是求最小割。
关键是解决点权问题。
我们可以拆点,拆成in和out
然后连边权值为cost

之后的边正常连,只不过权值都为inf代表不可以割这里的边。

然后跑个最大流就好了。

#include<iostream>
#include<algorithm>
using namespace std;
const int max_n = 500;
const int max_m = 20500;
const int inf = 1e9;
struct edge {
    int to,cap,next;
}E[max_m<<2];
int n;
int head[max_n],cnt=1,gap[max_n],last[max_n],d[max_n],que[max_n],ql,qr;
void add(int u,int to,int cap) {
    E[++cnt].to=to;E[cnt].cap=cap;
    E[cnt].next=head[u];head[u]=cnt;
    E[++cnt].to=u;E[cnt].cap=0;
    E[cnt].next=head[to];head[to]=cnt;
}int aug(int x,int s,int t,int mi) {
    if (x==t) return mi;
    int flow=0;
    for (int &i=last[x],v=E[i].to;i;i=E[i].next,v=E[i].to) if (d[x]==d[v]+1) {
        int tmp=aug(v,s,t,min(mi,E[i].cap));
        flow+=tmp,mi-=tmp,E[i].cap-=tmp,E[i^1].cap+=tmp;
        if (!mi) return flow;
    }
    if (!(--gap[d[x]])) d[s]=n+1;
    ++gap[++d[x]],last[x]=head[x];
    return flow;
}
int maxflow(int s,int t) {
    fill(gap,gap+n+10,0);fill(d,d+n+10,0);
    ++gap[d[t]=1];
    for (int i=1;i<=n;++i) last[i]=head[i];
    que[ql=qr=1]=t;
    while (ql<=qr) {
        int x=que[ql++];
        for (int i=head[x],v=E[i].to;i;i=E[i].next,v=E[i].to) if (!d[v]) ++gap[d[v]=d[x]+1],que[++qr]=v;
    }

    int ret=aug(s,s,t,inf);
    while (d[s]<=n)    ret+=aug(s,s,t,inf);
    return ret;
}
int N,M;
int main(){
    while(~scanf("%d %d",&N,&M)){
        fill(head,head+max_n,0);cnt=1;
        int s,t;scanf("%d %d",&s,&t);
        t+=N;n = N<<1;
        for (int i=1;i<=N;++i){
            int cost;scanf("%d",&cost);
            add(i,i+N,cost);
        }
        for (int i=1;i<=M;++i){
            int u,v;scanf("%d %d",&u,&v);
            add(u+N,v,inf);add(v+N,u,inf);
        }
        printf("%d\n",maxflow(s,t));
    }
}
kuangbin题单刷题详解(网络流) 文章被收录于专栏

题单:https://vjudge.net/article/371

全部评论

相关推荐

Hello_WordN:咱就是说,除了生命其他都是小事,希望面试官平安,希望各位平时也多注意安全
点赞 评论 收藏
分享
贺兰星辰:不要漏个人信息,除了简历模板不太好以外你这个个人简介是不是太夸大了...
点赞 评论 收藏
分享
1 收藏 评论
分享
牛客网
牛客企业服务