题目大意:
给你一个二分图,左边的点是1到n,右边的点是1到m,给你k条边,问你能否选出一些边,让每个点的度数在L,R之间
分析:
有上下界的网络流,很明显的一道模板题,也是我的第一道上下界网络流的题目,设置一个原点,一个汇点,在设置一个附件源y,和附加汇x,进行dinic算法。
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <queue>
#include <set>
#include <map>
#include <vector>
using namespace std;
const int maxn = 4000 +50;
const int INF=0x3f3f3f3f;
struct Edge
{
int from,to,cap,flow;
};
struct Dinic
{
int n,m,s,t;
vector<Edge>edges;
vector<int>G[maxn];
bool vis[maxn];
int d[maxn];
int cur[maxn];
void init(int n)
{
this->n=n;
for(int i=0;i<=n;i++) G[i].clear();
edges.clear();
}
void addEdge(int from,int to,int cap)
{
edges.push_back((Edge){from,to,cap,0});
edges.push_back((Edge){to,from,0,0});
m=edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
bool bfs()
{
memset(vis,0,sizeof(vis));
queue<int>Q;
Q.push(s);
d[s]=0;
vis[s]=1;
while(!Q.empty())
{
int x=Q.front(); Q.pop();
for(int i=0;i<G[x].size();i++)
{
Edge &e=edges[G[x][i]];
if(!vis[e.to]&&e.cap>e.flow)
{
vis[e.to]=1;
d[e.to]=d[x]+1;
Q.push(e.to);
}
}
}
return vis[t];
}
int dfs(int x,int a)
{
if(x==t||a==0) return a;
int flow=0,f;
for(int &i=cur[x];i<G[x].size();i++)
{
Edge &e=edges[G[x][i]];
if(d[x]+1==d[e.to]&&(f=dfs(e.to,min(a,e.cap-e.flow)))>0)
{
e.flow+=f;
edges[G[x][i]^1].flow-=f;
flow+=f;
a-=f;
if(a==0) break;
}
}
return flow;
}
int maxflow(int s,int t)
{
this->s=s; this->t=t;
int flow=0;
while(bfs())
{
memset(cur,0,sizeof(cur));
flow+=dfs(s,INF);
}
return flow;
}
void solve(int y,int kase)
{
bool flag=false;
for(int i=0;i<G[y].size();i++)
{
Edge &e=edges[G[y][i]];
// cout<<e.flow<<"..."<<e.cap<<"..."<<e.to<<endl;
if(e.cap!=e.flow) flag=true;
}
printf("Case %d: ",kase);
if(flag) printf("No\n");
else printf("Yes\n");
}
}D;
int main()
{
int n,m,k,kase=0;
while(~scanf("%d%d%d",&n,&m,&k))
{
int l,r;
scanf("%d%d",&l,&r);
D.init(4010);
int x=0,y=n+m+1;
int s=n+m+2,t=n+m+3;
D.addEdge(t,s,INF);
for(int i=1;i<=n;i++)
{
D.addEdge(s,i,r-l);
D.addEdge(s,x,l);
D.addEdge(y,i,l);
}
for(int i=1;i<=m;i++)
{
D.addEdge(i+n,t,r-l);
D.addEdge(i+n,x,l);
D.addEdge(y,t,l);
}
int u,v;
for(int i=1;i<=k;i++)
{
scanf("%d%d",&u,&v);
D.addEdge(u,v+n,1);
}
D.maxflow(y,x);
D.solve(y,++kase);
}
return 0;
}