【bzoj 2502】清理雪道（有源汇上下界最小流）

传送门biu~
题目相当于给一个拓扑图，每个边的流量下界为1，求最小流。
首先跑出可行流。建超级源点SS和超级汇点TT。原网络中的T->S连一条流量为inf的弧。把每个点必须流入的流量从SS流入，必须流出的流量流进TT。从SS到TT跑最大流即为原网络满足所有下界的可行流，我们就可以知道原网络中的流量是多少。
对于原网络在求完可行流后的残量网络，无法保证是最小流。这时把原网络以外的弧全删掉，从T至S跑最大流，相当于把尽量多的流流回S，把这部分流删掉即为最小流。

#include<bits/stdc++.h>
using namespace std;
const int inf=1e8,N=150,M=30050;
int n,head[N],degree[N],fir[N],dep[N],nex[M],to[M],cap[M],tp=1;
inline void add(int x,int y,int c){
    nex[++tp]=head[x];
    head[x]=tp;
    to[tp]=y;
    cap[tp]=c;
}
inline void Insert(int x,int y,int c){add(x,y,c);add(y,x,0);}
queue<int>q;
inline int bfs(int s,int t){
    memset(dep,0,sizeof dep);
    dep[s]=1;q.push(s);
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=nex[i]){
            if(!dep[to[i]] && cap[i]){
                dep[to[i]]=dep[x]+1;
                q.push(to[i]);
            }
        }
    }
    return dep[t];
}
int dfs(int x,int now,int t){
    if(x==t || !now)    return now;
    int c=0;
    for(int &i=fir[x];i;i=nex[i]){
        if(dep[to[i]]==dep[x]+1 && cap[i]){
            int f=dfs(to[i],min(now,cap[i]),t);
            now-=f;
            c+=f;
            cap[i]-=f;
            cap[i^1]+=f;
            if(!now)    break;
        }
    }
    return c;
}
inline int Dinic(int s,int t){
    int c=0;
    while(bfs(s,t)){
        for(int i=1;i<=n+4;++i) fir[i]=head[i];
        c+=dfs(s,inf,t);
    }
    return c;
}
int main(){
    scanf("%d",&n);int S=n+1,T=n+2,SS=n+3,TT=n+4;
    for(int i=1;i<=n;++i){
        int num;
        scanf("%d",&num);
        while(num--){
            int x;
            scanf("%d",&x);
            Insert(i,x,inf);
            ++degree[x],--degree[i];
        }
    }
    for(int i=1;i<=n;++i)   Insert(S,i,inf),Insert(i,T,inf);
    for(int i=1;i<=n;++i){
        if(degree[i]>0)     Insert(SS,i,degree[i]);
        if(degree[i]<0)     Insert(i,TT,-degree[i]);
    }
    Insert(T,S,inf);
    Dinic(SS,TT);
    int sum=cap[tp];cap[tp]=cap[tp^1]=0;
    for(int i=head[SS];i;i=nex[i])  cap[i]=cap[i^1]=0;
    for(int i=head[TT];i;i=nex[i])  cap[i]=cap[i^1]=0;
    printf("%d\n",sum-Dinic(T,S));
    return 0;
}