HDU 4292 Food (拆点最大流)

题意：N个人，F种食物，D种饮料，给定每种食物和饮料的量。每个人有自己喜欢的食物和饮料，如果得到自己喜欢的食物和饮料才能得到满足。求最大满足的人数。

分析：如果只是简单地N个人选择F种食物的话可以用二分图匹配来完成，但是该题种一个人要选择两样东西。

采取以下操作建图：

　　增加源点和汇点，源点与F个食物相连，每条弧流量为第i种食物的数量；D个食物与汇点相连，每条弧流量为其数量。

　　所以将N个人拆成两个点，对应两点之间连一条边，流量为1。

然后跑出最大流就是答案。

#include<bits/stdc++.h>
using namespace std;
#define captype int
const int MAXN = ;
const int MAXM = ;
const int INF = <<;
struct EDGE{
    int to,next;
    captype cap;
} edg[MAXM];
int eid,head[MAXN];
int gap[MAXN];
int dis[MAXN];
int cur[MAXN];
int pre[MAXN];

void init(){
    eid=;
    memset(head,-,sizeof(head));
}
void AddEdge(int u,int v,captype c,captype rc=){
    edg[eid].to=v; edg[eid].next=head[u];
    edg[eid].cap=c;  head[u]=eid++;
    edg[eid].to=u; edg[eid].next=head[v];
    edg[eid].cap=rc; head[v]=eid++;
}
captype maxFlow_sap(int sNode,int eNode, int n){//n是包括源点和汇点的总点个数，这个一定要注意
    memset(gap,,sizeof(gap));
    memset(dis,,sizeof(dis));
    memcpy(cur,head,sizeof(head));
    pre[sNode] = -;
    gap[]=n;
    captype ans=;
    int u=sNode;
    while(dis[sNode]<n){
        if(u==eNode){
            captype Min=INF ;
            int inser;
            for(int i=pre[u]; i!=-; i=pre[edg[i^].to])
            if(Min>edg[i].cap){
                Min=edg[i].cap;
                inser=i;
            }
            for(int i=pre[u]; i!=-; i=pre[edg[i^].to]){
                edg[i].cap-=Min;
                edg[i^].cap+=Min;
            }
            ans+=Min;
            u=edg[inser^].to;
            continue;
        }
        bool flag = false;
        int v;
        for(int i=cur[u]; i!=-; i=edg[i].next){
            v=edg[i].to;
            if(edg[i].cap> && dis[u]==dis[v]+){
                flag=true;
                cur[u]=pre[v]=i;
                break;
            }
        }
        if(flag){
            u=v;
            continue;
        }
        int Mind= n;
        for(int i=head[u]; i!=-; i=edg[i].next)
        if(edg[i].cap> && Mind>dis[edg[i].to]){
            Mind=dis[edg[i].to];
            cur[u]=i;
        }
        gap[dis[u]]--;
        if(gap[dis[u]]==) return ans;
        dis[u]=Mind+;
        gap[dis[u]]++;
        if(u!=sNode) u=edg[pre[u]^].to;  //退一条边
    }
    return ans;
}

int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    int S,T,N,M,F,D,u,v,tmp;
    while(scanf("%d%d%d",&N,&F,&D)==){
        S=,T=*N+F+D+;            //源点，汇点
        init();
        for(int i=;i<=F;++i){
            scanf("%d",&tmp);
            AddEdge(S,i,tmp);
        }
        for(int i=;i<=D;++i){
            scanf("%d",&tmp);
            AddEdge(F+N*+i,T,tmp);
        }
        for(int i=;i<=N;++i){
            AddEdge(i+F,i+F+N,);
            char op[MAXN];
            scanf("%s",op);
            for(int j=;j<=F;++j){
                if(op[j-]=='Y') AddEdge(j,i+F,);
            }
        }
        for(int i=;i<=N;++i){
            char op[MAXN];
            scanf("%s",op);
            for(int j=;j<=D;++j){
                if(op[j-]=='Y') AddEdge(F+N+i,F+N*+j,);
            }
        }
        printf("%d\n",maxFlow_sap(S,T,T+));
    }
    return ;
}