教辅的组成_洛谷1231_网络流

题目背景

---

滚粗了的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）

输出格式：

---

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

Analysis

---

第一眼是想用匹配做的，然后我就卡住了
因为是以每本书延伸的关系图，所以对于一本书来讲它只能跟一本答案和一本练习册配套，那么就拆点
连边st->练习册->书的入点->出点->答案->ed
跑最大流，dinic不会TLE

Code

---

#include <stdio.h>
#include <string.h>
#include <queue>
#define rep(i, a, b) for (int i = a; i <= b; i++)
#define fill(x, t) memset(x, t, sizeof(x))
#define min(x, y) x<y?x:y
#define INF 0x3f3f3f3f
#define N 50001
#define E 160001
using namespace std;
struct edge{int y, w, next;}e[E];
int dis[N], ls[N], edgeCnt;
inline int read(){
    int x = 0, v = 1; char ch = getchar();
    while (ch < '0' || ch > '9'){
        if (ch == '-'){
            v = -1;
        }
        ch = getchar();
    }
    while (ch >= '0' && ch <= '9'){
        x = x * 10 + ch - '0';
        ch = getchar();
    }
    return x * v;
}
int add(int x, int y, int w){
    e[++edgeCnt] = (edge){y, w, ls[x]}; ls[x] = edgeCnt;
    e[++edgeCnt] = (edge){x, 0, ls[y]}; ls[y] = edgeCnt;
    return 0;
}
int bfs(int st, int ed){
    queue<int>q;
    fill(dis, 0);
    q.push(st);
    dis[st] = 1;
    while (!q.empty()){
        int now = q.front(); q.pop();
        for (int i = ls[now]; i; i = e[i].next){
            if (!dis[e[i].y] && e[i].w > 0){
                dis[e[i].y] = dis[now] + 1;
                q.push(e[i].y);
                if (e[i].y == ed){
                    return 1;
                }
            }
        }
    }
    return 0;
}
int find(int now, int ed, int mn){
    if (now == ed || !mn){
        return mn;
    }
    int ret = 0;
    for (int i = ls[now]; i; i = e[i].next){
        if (e[i].w > 0 && dis[now] + 1 == dis[e[i].y]){
            int d = find(e[i].y, ed, min(mn - ret, e[i].w));
            e[i].w -= d;
            e[i ^ 1].w += d;
            ret += d;
            if (ret == mn){
                break;
            }
        }
    }
    return ret;
}
int dinic(int st, int ed){
    int mxFlow = 0;
    while (bfs(st, ed)){
        mxFlow += find(st, ed, INF);
    }
    return mxFlow;
}
int main(void){
    int n1 = read(), n2 = read(), n3 = read();
    int m1 = read();
    edgeCnt = 1;
    rep(i, 1, m1){
        int x = read(), y = read();
        add(y, n2 + x, 1);
    }
    int m2 = read();
    rep(i, 1, m2){
        int x = read(), y = read();
        add(n2 + n1 + x, n2 + n1 * 2 + y, 1);
    }
    rep(i, 1, n1){
        add(n2 + i, n2 + n1 + i, 1);
    }
    int st = 0, ed = n2 + n1 * 2 + n3 + 1;
    rep(i, 1, n2){
        add(st, i, 1);
    }
    rep(i, 1, n3){
        add(n2 + n1 * 2 + i, ed, 1);
    }
    int ans = dinic(st, ed);
    printf("%d\n", ans);
    return 0;
}