[洛谷1231]教辅的组成

思路：

最大流跑三分图匹配，注意每本书只能用一次，所以把每本书拆成两个点，连一条边。
不能直接用EdmondsKarp算法，也不能直接用不加优化的Dinic算法，这样会TLE7个点。
本题正解是Dinic算法加上当前弧优化。

#include<queue>
#include<cstdio>
#include<cctype>
#include<vector>
#include<cstring>
#include<cstdlib>
inline int getint() {
    char ch;
    while(!isdigit(ch=getchar()));
    int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int inf=0x7fffffff;
int s,t;
struct Edge {
    int from,to;
    bool remain;
};
const int E=1400000,V=400002;
Edge e[E];
std::vector<int> g[V];
int sz=0;
inline void add_edge(const int u,const int v,const bool w) {
    e[sz]=(Edge){u,v,w};
    g[u].push_back(sz);
    sz++;
}
int d[V];
inline void SPFA() {
    d[s]=0;
    for(register int i=1;i<=t;i++) d[i]=inf;
    std::queue<int> q;
    q.push(s);
    while(!q.empty()) {
        int x=q.front();
        q.pop();
        if(x==t) return;
        for(register unsigned i=0;i<g[x].size();i++) {
            Edge &y=e[g[x][i]];
            if((d[x]+1<d[y.to])&&y.remain) {
                d[y.to]=d[x]+1;
                q.push(y.to);
            }
        }
    }
}
unsigned cur[V];
bool Augment(const int x,const bool remain) {
    if(x==t) return remain;
    for(unsigned &i=cur[x];i<g[x].size();i++) {
        Edge &y=e[g[x][i]];
        if((d[y.to]>d[x])&&y.remain) {
            if(Augment(y.to,std::min(remain,y.remain))) {
                e[g[x][i]].remain^=true;
                e[g[x][i]^1].remain^=true;
                return true;
            }
        }
    }
    return false;
}
inline int Dinic() {
    int maxflow=0;
    for(SPFA();d[t]!=inf;SPFA()) {
        memset(cur,0,sizeof cur);
        while(Augment(s,true)) maxflow++;
    }
    return maxflow;
}
int main() {
    int n2=getint(),n1=getint(),n3=getint();
    s=0,t=n1+n2*2+n3+1;
    for(register int i=1;i<=n1;i++) {
        add_edge(s,i,true);
        add_edge(i,s,false);
    }
    for(register int m=getint();m;m--) {
        int a=getint(),b=getint();
        add_edge(b,n1+a,true);
        add_edge(n1+a,b,false);
    }
    for(register int i=1;i<=n2;i++) {
        add_edge(n1+i,n1+n2+i,true);
        add_edge(n1+n2+i,n1+i,false);
    }
    for(register int m=getint();m;m--) {
        int a=getint(),b=getint();
        add_edge(n1+n2+a,n1+n2*2+b,true);
        add_edge(n1+n2*2+b,n1+n2+a,false);
    }
    for(register int i=1;i<=n3;i++) {
        add_edge(n1+n2*2+i,t,true);
        add_edge(t,n1+n2*2+i,false);
    }
    printf("%d
",Dinic());
    return 0;
}