点这里 |
---|
题意: 有n头牛、f种食物和d种饮料。一种食物或一种饮料只能喂给一头牛,一头牛也只能吃一种食物和一种饮料。问最多能投喂几头牛。
题解: 看上去很想二分图的匹配问题,但是牛既要对应食物也要对应饮料,所以我只能用最大流的模板来解决。手握模板,什么题目都能套,但是如何建图?
需要用s-t流将一头牛和它喜欢的食物和饮料连接,而食物和饮料之间没有直接关系,自然就想到把牛放在中间,两边是食物和饮料,即s_食物_牛_饮料_t,边的权值都为1,故如果形成一条这样的流,就构成一种分配方案。但是为了限制一头牛只能吃一种食物和饮料,需要把牛拆分,一头与食物连接,一头与饮料连接,再把这头牛的两部分连接起来,以达到限制流的目的。
图片来自这篇博客点这里。
过程中犯的错:
- 初始化: 因为设置了一个超级源点S,而S的编号为0,各种数组初始化的时候,千万不要忘了0.
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
const int inf = 0x3f3f3f3f;
const int maxn = 410;
const int maxm = 30000;
int N, F, D;
int s, t, cnt;
int cur[maxn];
int head[maxn];
int depth[maxn];
struct edge{ int v, w, next;} e[maxm];
void add(int u, int v, int w){ e[cnt].v = v, e[cnt].w = w, e[cnt].next = head[u], head[u] = cnt++;}
void addedge(int u, int v, int w){ add(u, v, w); add(v, u, 0);} //有向图,注意加的第二条边
void init(){ s = 0, t = 401, cnt = 0; for(int i = 0; i < maxn; i++) head[i] = -1;}
bool bfs(){
for(int i = 0; i < maxn; i++) depth[i] = 0;
queue<int> Q;
Q.push(s); depth[s] = 1;
while(!Q.empty()){
int u = Q.front(); Q.pop();
for(int i = head[u]; i != -1; i = e[i].next){
int v = e[i].v, w = e[i].w;
if(depth[v] == 0 && w > 0){ //点v没有去过并且还有残留容量
depth[v] = depth[u] + 1;
Q.push(v);
}
}
}
if(depth[t] > 0) return 1;
return 0;
}
int dfs(int u, int flow){
if(u == t) return flow;
for(int& i = cur[u]; ~i; i = e[i].next){ //&符号,改变i的同时也会改变cur[u]
int v = e[i].v, w = e[i].w;
if(depth[v] == depth[u] + 1 && w != 0){ //点v没有去过并且还有残留容量
int minflow = dfs(v, min(flow, w));
if(minflow > 0){
e[i].w -= minflow;
e[i ^ 1].w += minflow; //更新残留容量
return minflow;
}
}
}
return 0;
}
int Dinic(){
int ans = 0;
int k = 0;
while(bfs()){
for(int i = 0; i < maxn; i++) cur[i] = head[i];
while(int d = dfs(s, inf)){
ans += d;
}
}
return ans;
}
int main(){
while(~scanf("%d%d%d", &N, &F, &D)){
init();
for(int i = 1; i <= N; i++){
int f, d; scanf("%d%d", &f, &d);
while(f--){
int x; scanf("%d", &x);
addedge(x, i + 100, 1);
}
addedge(i + 100, i + 200, 1);
while(d--){
int x; scanf("%d", &x);
addedge(i + 200, x + 300, 1);
}
}
for(int i = 1; i <= F; i++) addedge(0, i, 1);
for(int i = 1; i <= D; i++) addedge(i + 300, 401, 1);
printf("%d\n", Dinic());
}
return 0;
}