教辅的组成——dinic多路增广

　　这题就是裸的dinic，随便模板打一个就是了。

　　因为每本书只能选一次，我们就考虑在书中连一条边，就是：书的入点——边——书的出点。

　　然后整个图建成这样：超级源点s->练习册->书的入点->出点->答案->超级汇点t。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector> 
#include<queue>
using namespace std;
inline int read(){
    int t=1,num=0;
    char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')t=-1;c=getchar();}
    while(c>='0'&&c<='9'){num=num*10+c-'0';c=getchar();}
    return num*t;
}//快读，不说了 
const int INF=10000010;
struct edge{
    int to;int c;int rev;
};//to表示边的终点，c表示边的流量，rev表示该边的反向边 
vector<edge> g[40010];
int n,m,iter[40010],level[40010];
void add(int f,int t,int c){//读边 
    g[f].push_back((edge){t,c,g[t].size()});
    g[t].push_back((edge){f,0,g[f].size()-1});
}
void bfs(int s){//形成一个分层图 
    memset(level,0,sizeof(level));
    //level是点的层数，从1开始保存 
    queue<int>q;
    level[s]=1;
    q.push(s);
    while(!q.empty()){
        int v=q.front();q.pop();
        for(int i=0;i<g[v].size();i++){
            edge &e=g[v][i];
            if(e.c>0&&level[e.to]==0){
                level[e.to]=level[v]+1;
                q.push(e.to);
            }
        }
    }
}
int dfs(int v,int t,int f){//dinic 多路增广 
    if(v==t) return f;
    int used=0,d;
    for(int &i=iter[v];i<g[v].size();i++){
        edge &e=g[v][i];
        if(e.c>0&&level[v]+1==level[e.to]){
            d=f-used;
            d=dfs(e.to,t,min(d,e.c));
            if(d>0){
                e.c-=d;
                g[e.to][e.rev].c+=d;
                used+=d;
                if(used==f)return f;
            }
        }
    }
    return used;
}
int flow(int s,int t){//求解 
    int flow=0;
    while(1){
        bfs(s);
        if(level[t]==0) return flow;
        memset(iter,0,sizeof(iter));
        flow+=dfs(s,t,INF);
    }
}
int main()
{
    int n1,n2,n3,m;
    n1=read();n2=read();n3=read();
    m=read();
    for(int i=1;i<=m;i++){
        int x,y;
        x=read();y=read();
        add(y,x+n2,1);
    }
    //写得有些繁琐，请不要介意 
    for(int i=1;i<=n1;i++) add(i+n2,i+n2+n1,1);
    m=read();
    for(int i=1;i<=m;i++){
        int x,y;
        x=read();y=read();
        add(x+n1+n2,y+n1+n1+n2,1);
    }
    for(int i=1;i<=n2;i++) add(0,i,1);
    for(int i=1;i<=n3;i++) add(i+n1+n1+n2,n1+n1+n2+n3+1,1);
    printf("%d\n",flow(0,n1+n1+n2+n3+1));
    return 0;
}

  本文由Yzyet编写，网址为www.cnblogs.com/Yzyet。非Yzyet同意，禁止转载，侵权者必究。