洛谷 P1231 教辅的组成（最大流） 题解

题目来源：

https://www.luogu.org/problemnew/show/P1231

题目描述：

 

题目背景

滚粗了的HansBug在收拾旧语文书，然而他发现了什么奇妙的东西。

题目描述

蒟蒻HansBug在一本语文书里面发现了一本答案，然而他却明明记得这书应该还包含一份练习题。然而出现在他眼前的书多得数不胜数，其中有书，有答案，有练习册。已知一个完整的书册均应该包含且仅包含一本书、一本练习册和一份答案，然而现在全都乱做了一团。许多书上面的字迹都已经模糊了，然而HansBug还是可以大致判断这是一本书还是练习册或答案，并且能够大致知道一本书和答案以及一本书和练习册的对应关系（即仅仅知道某书和某答案、某书和某练习册有可能相对应，除此以外的均不可能对应）。既然如此，HansBug想知道在这样的情况下，最多可能同时组合成多少个完整的书册。

输入输出格式

输入格式：

 

第一行包含三个正整数N1、N2、N3，分别表示书的个数、练习册的个数和答案的个数。

第二行包含一个正整数M1，表示书和练习册可能的对应关系个数。

接下来M1行每行包含两个正整数x、y，表示第x本书和第y本练习册可能对应。（1<=x<=N1，1<=y<=N2）

第M1+3行包含一个正整数M2，表述书和答案可能的对应关系个数。

接下来M2行每行包含两个正整数x、y，表示第x本书和第y本答案可能对应。（1<=x<=N1，1<=y<=N3）

 

输出格式：

 

输出包含一个正整数，表示最多可能组成完整书册的数目。

 

输入输出样例

输入样例#1： 复制

5 3 4
5
4 3
2 2
5 2
5 1
5 3
5
1 3
3 1
2 2
3 3
4 3

输出样例#1： 复制

2

说明

样例说明：

如题，N1=5，N2=3，N3=4，表示书有5本、练习册有3本、答案有4本。

M1=5，表示书和练习册共有5个可能的对应关系，分别为：书4和练习册3、书2和练习册2、书5和练习册2、书5和练习册1以及书5和练习册3。

M2=5，表示数和答案共有5个可能的对应关系，分别为：书1和答案3、书3和答案1、书2和答案2、书3和答案3以及书4和答案3。

所以，以上情况的话最多可以同时配成两个书册，分别为：书2+练习册2+答案2、书4+练习册3+答案3。

数据规模：

对于数据点1, 2, 3，M1，M2<= 20

对于数据点4~10，M1，M2 <= 20000

解题思路：

     最大流的应用，dinic是比较好写，套个模板就行，关键就是建图，首先先建一个源点（0）和汇点（n1*2+n2+n3+1)，然后把每一本书拆成两本书中间连流量为1的边，然后练习册（1-n2）往书1（n2+1--n2+n1)连边，书2（n1+n2+1--n1*2+n2）在向答案(n1*2+n2+1--n1*2+n2+n3）连边,接下来就可以跑一遍最大流了。。。。。

代码：

#include <iostream>
#include <string>
#include <cstring>
#include <vector>
#include <queue>
#include <stack>
#include <algorithm>
#include <cmath>
#include <cctype>
#include <iomanip>
#define inf 1e9
const int maxn=20005;
using namespace std;
struct newt{
	int to,next,val;
}edge[15*maxn];
int n1,n2,n3,m1,m2,S,T,head[5*maxn],deep[5*maxn],cur[5*maxn],cnt;
void addedge(int u,int v,int w)
{
	edge[cnt].to=v;
	edge[cnt].val=w;
	edge[cnt].next=head[u];
	head[u]=cnt++;
}
inline bool bfs()
{
	memset(deep,-1,sizeof(deep));
	queue<int>q;
	deep[S]=0;q.push(S);
	//for(int i=0;i<=n1*2+n2+n3+1;i++)cur[i]=head[i];
	while(!q.empty())
	{
		int now=q.front();q.pop();
		for(int i=head[now];i!=-1;i=edge[i].next){
			int v=edge[i].to;
			if(deep[v]==-1&&edge[i].val)deep[v]=deep[now]+1,q.push(v);
		}
	}
	return deep[T]!=-1;
}
inline int dfs(int u,int op)
{
	if(u==T||op==0)return op;
	int flow=0,used=0;
	for(int i=head[u];i!=-1;i=edge[i].next)
	{
		//cur[u]=i;
		int v=edge[i].to;
		if(deep[v]==deep[u]+1&&edge[i].val){
			used=dfs(v,min(op,edge[i].val));
			if(!used)continue;
			flow+=used;op-=used;
			edge[i].val-=used;edge[i^1].val+=used;
			if(!op)break;
		}
	}
	if(!flow)deep[u]=-1;
	return flow;
}
int main()
{
	ios::sync_with_stdio(false);
	cnt=0;
	cin>>n1>>n2>>n3;
	cin>>m1;
	S=0;T=n1*2+n2+n3+1;
	int ans=0;
	memset(head,-1,sizeof(head));
	int x,y,z;
	for(int i=1;i<=n2;i++)
	addedge(0,i,1),addedge(i,0,0);
	for(int i=1;i<=m1;i++)
	{
		cin>>x>>y;
		addedge(y,x+n2,1);
		addedge(x+n2,y,0);
	}
	for(int i=1;i<=n1;i++)
	addedge(i+n2,i+n2+n1,1),addedge(i+n2+n1,i+n2,0);
	cin>>m2;
	for(int i=1;i<=m2;i++)
	{
		cin>>x>>z;
		addedge(x+n1+n2,z+n1*2+n2,1);
		addedge(z+n1*2+n2,x+n1+n2,0);
	}
	for(int i=1;i<=n3;i++)
	addedge(i+n1*2+n2,n1*2+n2+n3+1,1),addedge(n1*2+n2+n3+1,n1*2+n2+i,0);
	while(bfs()){
		ans+=dfs(S,2*inf);
	}
	cout<<ans<<endl;
	return 0;
}