【有源汇上下界可行流】ACM-ICPC 2018 沈阳赛区网络预赛 - F - Fantastic Graph

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/monochrome00/article/details/82533495

题目链接<https://nanti.jisuanke.com/t/31447>


题意:

给出一张二分图,初始每个节点的度数都为零。选择若干条边,使得每个节点的度数范围再[L,R]范围内。每选一条边,边上两端的节点度数+1。


题解:

关于有上下界的网络流:https://blog.csdn.net/u011008379/article/details/38306477

源点与左边每个节点连[L,R]的边。右边每个节点与汇点连[L,R]的边。

左右两边按照题意连权值为1的边。

最后判断有无可行流即可。


#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<iostream>
using namespace std;
const int N=1e5+7;
const int inf=1<<26;
struct Edge{
    int u,v,w,nxt;
    Edge(int u=0,int v=0,int w=0,int nxt=0):u(u),v(v),w(w),nxt(nxt){}
}edge[30*N];
int n,m,edn,sp,tp;
int p[N],d[N],c[N];
void add(int u,int v,int w){
    edge[edn]=Edge(u,v,w,p[u]);p[u]=edn++;
    edge[edn]=Edge(v,u,0,p[v]);p[v]=edn++;
}
bool bfs(){
    memset(d,-1,sizeof(d));d[sp]=0;
    queue<int>q;q.push(sp);
    while(!q.empty()){
        int u=q.front();
        q.pop();
        for(int i=p[u];~i;i=edge[i].nxt){
            int v=edge[i].v;
            if(d[v]==-1&&edge[i].w){
                d[v]=d[u]+1;
                q.push(v);
                if(v==tp) return true;
            }
        }
    }
    return ~d[tp];
}
int dfs(int u,int b){
    if(u==tp) return b;
    int r=0;
    for(int i=c[u];~i;i=edge[i].nxt){
        int v=edge[i].v;
        if(edge[i].w&&d[v]==d[u]+1){
            int x=min(edge[i].w,b-r);
            c[u]=i;
            x=dfs(v,x);
            r+=x;
            edge[i].w-=x;
            edge[i^1].w+=x;
            if(r==b) break;
        }
    }
    if(!r)d[u]=-2;
    return r;
}
int dinic(){
    int total=0,t;
    while(bfs()){
        memcpy(c,p,sizeof(p));
        while(t=dfs(sp,inf))
        total+=t;
    }
    return total;
}
int x[N];
int main()
{
    int k,lo,hi,cs=0;
    while(scanf("%d%d%d",&n,&m,&k)!=EOF){
        memset(p,-1,sizeof(p));edn=0;
        sp=n+m+1,tp=n+m+2;
        int ss=n+m+3,tt=n+m+4;
        scanf("%d%d",&lo,&hi);
        int u,v;
        for(int i=1;i<=k;i++){
            scanf("%d%d",&u,&v);
            add(u,v+n,1);
        }
        for(int i=1;i<=n;i++){
            add(ss,i,hi-lo);
            add(sp,i,lo);
            add(ss,tp,lo);
        }
        for(int i=1;i<=m;i++){
            add(i+n,tt,hi-lo);
            add(sp,tt,lo);
            add(i+n,tp,lo);
        }
        add(tt,ss,inf);
        printf("Case %d: ",++cs);
        if(dinic()==(n+m)*lo) puts("Yes");
        else puts("No");
    }
    return 0;
}

猜你喜欢

转载自blog.csdn.net/monochrome00/article/details/82533495