ACM-ICPC 2018 沈阳赛区网络预赛 F题 Fantastic Graph

       题目大意：给定一个二分图，左边有n个点，右边有m个点，图中有k条边，问能否删去其中若干条边或不变，使得每个点所连的边数量都在给定的区间[L,R]内。

       输入：多组数据，最多三十组。对于每组数据，第一行三个数字n(1<=n<=2000),m(1<=m<=2000)和k(0<=k<=6000)；第二行两个数字L和R(0<=L<=R<=300)；接下来k行每行两个数字，表示从左边的某个点到右边的某个点有一条边。

       输出：每组数据输出一行，第i组数据格式为Case i: Yes/No。

       题解：

       一道有上下界网络流裸题。

       教程可以参照这个网址：https://www.cnblogs.com/kane0526/archive/2013/04/05/3001108.html

       具体到这道题的话，建模就是：二分图的每条边流量为1，建一个源点和一个汇点，由源点到左边每个点连一条流量为R-L的边，再由右边每个点到汇点连一条流量为R-L的边，由汇点向源点连一条流量无限大的边。之后再建一个超级源点和一个超级汇点，由超级源点向左边每个点连一条流量为L的边，由右边每个点向超级汇点连一条流量为L的边，之后从超级源点到超级汇点跑一遍最大流，如果超级源点和超级汇点所连的边全部跑满，则存在可行解，否则不存在。

       （比赛时再次把锅甩给了队友……以下代码是队友写的）

        代码：

/*To The End Of The Galaxy*/
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<iomanip>
#include<stack>
#include<map>
#include<set>
#include<cmath>
#define debug(x) cerr<<#x<<"="<<x<<endl
#define INF 0x7f7f7f7f
#define llINF 0x7fffffffffffll
#define P(x,y) (((x-1)*m)+y)
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
inline int init()
{
    int now=0,ju=1;char c;bool flag=false;
    while(1)
    {
        c=getchar();
        if(c=='-')ju=-1;
        else if(c>='0'&&c<='9')
        {
            now=now*10+c-'0';
            flag=true;
        }
        else if(flag)return now*ju;
    }
}
inline long long llinit()
{
    long long now=0,ju=1;char c;bool flag=false;
    while(1)
    {
        c=getchar();
        if(c=='-')ju=-1;
        else if(c>='0'&&c<='9')
        {
            now=now*10+c-'0';
            flag=true;
        }
        else if(flag)return now*ju;
    }
}
struct edge
{
    int from,to,cap,flow,pre;
}Edge[100005];
int head[4010],dis[4010],cur[4010];
bool vis[4010];
int n,m,k,cnt,L,R;
inline void addedge(int from,int to,int cap)
{
    ++cnt;
    Edge[cnt]=((edge){from,to,cap,0,head[from]});
    head[from]=cnt;
    ++cnt;
    Edge[cnt]=((edge){to,from,0,0,head[to]});
    head[to]=cnt;
}
int S,T,s,t;
queue<int> q;
bool bfs()
{
    int now;
    while(!q.empty())q.pop();
    memset(vis,0,sizeof(vis));
    dis[S]=0;
    q.push(S);
    while(!q.empty())
    {
        now=q.front();q.pop();
        if(now==T)return true;
        if(vis[now])continue;
        vis[now]=1;
        for(int j=head[now];j;j=Edge[j].pre)
        {
            if(!vis[Edge[j].to]&&Edge[j].cap>Edge[j].flow)
            {
                dis[Edge[j].to]=dis[now]+1;
                q.push(Edge[j].to);
            }
        }
    }
    return false;
}
inline int dfs(int now,int maxflow)
{
    if(now==T||maxflow==0)return maxflow;
    int &j=cur[now];
    int flow=0,f;
    for(;j;j=Edge[j].pre)
    {
        if(dis[Edge[j].to]==dis[now]+1&&(f=dfs(Edge[j].to,min(maxflow,Edge[j].cap-Edge[j].flow)))>0)
        {
            flow+=f;maxflow-=f;
            Edge[j].flow+=f;Edge[((j-1)^1)+1].flow-=f;
            if(maxflow==0)break;
        }
    }
    return flow;
}
int dinic()
{
    int ans=0;
    while(bfs())
    {
        for(register int i=1;i<=T;i++)cur[i]=head[i];
        ans+=dfs(S,INF);
    }
    return ans;
}
int main()
{
	int fr,to,ans,d=0;
	bool flag;
	while(scanf("%d%d%d",&n,&m,&k)!=EOF)
	{
		d++;
		memset(head,0,sizeof(head));
		memset(vis,0,sizeof(vis));
		memset(cur,0,sizeof(cur));
		memset(dis,0,sizeof(dis));
		flag=true;
		s=n+m+1;t=s+1;S=t+1;T=S+1;
		cnt=0;
		L=init();R=init();
		for(int i=1;i<=k;i++)
		{
			fr=init();to=init();
			addedge(fr,to+n,1);
		}
		for(int i=1;i<=n;i++)
		{
			addedge(s,i,R-L);
			addedge(S,i,L);
		}
		for(int i=1;i<=m;i++)
		{
			addedge(i+n,T,L);
			addedge(i+n,t,R-L);
		}
		ans=dinic();if(ans!=L*m)flag=false;
		for(int j=head[S];j;j=Edge[j].pre)
		{
			if(Edge[j].cap-Edge[j].flow!=0)
			{
				flag=false;
			}
		}
		if(flag)printf("Case %d: Yes\n",d);
		else printf("Case %d: No\n",d);
	}
	return 0;
}
/*
srO xudyh davidlee1999WTK linkct1999 zlser Orz
*/