hdu 4292(拆点的网络流)

题意：有F种食物 D种饮料 它们都有一定的数量 有N个人 每个人都有自己喜欢吃的食物和饮料 （每个人至少要一种食物和饮料） 只有能满足他的要求时他才会接服务 求最大能满足多少人？

思路：网络流 建一超级源点 汇点 源点与食物相连 边权为其数量，汇点与饮料相连 边权也为其数量 把人分成两个点 之间的边权为1 每个人与之需要的食物和饮料相连 边权为1 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<set>
#include<cstdio>
#include<string>
using namespace std;
const int maxn=1005;
const int inf=1<<30;
struct node
{
    int v,next;
    int val;
} s[maxn*maxn*2];
int level[maxn];//顶点的层次
int p[maxn];
int que[maxn*10];//BFS中用于遍历的顶点，DFS求增广中记录边
int out[10*maxn];//DFS用于几乎定点的分支
int ind;
inline void insert(int x,int y,int z)
{
    s[ind].v=y;
    s[ind].val=z;
    s[ind].next=p[x];
    p[x]=ind++;
    s[ind].v=x;
    s[ind].val=0;
    s[ind].next=p[y];
    p[y]=ind++;
}
int n,f,d;

void init()
{
    ind=0;
    memset(p,-1,sizeof(p));
}
int max_flow(int n,int source,int sink)
{
    int ret=0;
    int h=0,r=0;
    while(1)//DFS
    {
        int i;
        for(i=0; i<=n; ++i)
            level[i]=0;
        h=0,r=0;

        level[source]=1;
        que[0]=source;
        while(h<=r)//BFS
        {
            int  t=que[h++];
            for(i=p[t]; i!=-1; i=s[i].next)
            {
                if(s[i].val&&level[s[i].v]==0)
                {
                    level[s[i].v]=level[t]+1;
                    que[++r]=s[i].v;
                }
            }
        }
        //topset=r;//记录原点的集合个数
        if(level[sink]==0)break;//找不到汇点
        for(i=0; i<=n; ++i)
            out[i]=p[i];
        int  q=-1;
        while(1)
        {
            if(q<0)
            {
                int  cur=out[source];
                for(; cur!=-1; cur=s[cur].next)
                {
                    if(s[cur].val&&out[s[cur].v]!=-1&&level[s[cur].v]==2)
                    {
                        break;
                    }
                }
                if(cur>=0)
                {
                    que[++q]=cur;
                    out[source]=s[cur].next;
                }
                else
                {
                    break;
                }
            }


            int  u=s[que[q]].v;


            if(u==sink)//一条增广路
            {
                int  dd=inf;
                int  index=-1;
                for(i=0; i<=q; i++)
                {
                    if(dd>s[que[i]].val)
                    {
                        dd=s[que[i]].val;
                        index=i;
                    }
                }
                ret+=dd;
                //cout<<ret<<endl;
                for(i=0; i<=q; i++)
                {
                    s[que[i]].val-=dd;
                    s[que[i]^1].val+=dd;
                }
                for(i=0; i<=q; i++)
                {
                    if(s[que[i]].val==0)
                    {
                        q=index-1;
                        break;
                    }
                }
            }
            else
            {
                long cur=out[u];
                for(; cur!=-1; cur=s[cur].next)
                {
                    if(s[cur].val&&out[s[cur].v]!=-1&&level[u]+1==level[s[cur].v])
                    {
                        break;
                    }
                }
                if(cur!=-1)
                {
                    que[++q]=cur;
                    out[u]=s[cur].next;
                }
                else
                {
                    out[u]=-1;
                    q--;
                }
            }
        }
    }
    return ret;
}

int main()
{
    while( scanf("%d%d%d",&n,&f,&d)!=EOF )
    {
        ind=0;
        for(int i=0; i<=2+f+2*n+d; i++)
            p[i]=-1;
        int start=0;
        int end=1+f+2*n+d;
        int aa;
        for (int i = 1; i <= f; i++)
        {
            scanf("%d",&aa);
            insert(start, i,aa);
            // maxflow.addedge(i, s, 0);
        }
        for (int i = 1; i <= d; i++)
        {
            scanf("%d",&aa);
            insert(f + n * 2 + i, end, aa);
        }
        //for(int i=1;i<=n;i++)
        //{
        /*int fn,dn;
        cin>>fn>>dn;
        for(int j=0;j<fn;j++)
        {
            int ff;
            cin>>ff;
            //insert(0,ff,1);
            insert(ff,f + i * 2 - 1,1);
        }
        for(int j=0;j<dn;j++)
        {
            int dd;
            cin>>dd;
            //insert(dd+2*n+f,end,1);
            insert(f + i * 2, f + n * 2 + dd, 1);
        }*/
        //}
        char ch[300];
        for(int i=1; i<=n; i++)
        {
            cin>>ch;
            for(int j=0; j<f; j++)
            {
                if(ch[j]=='Y')
                  insert(j+1,f + i * 2 - 1,1);
            }
        }
        
        for(int i=1; i<=n; i++)
        {
            char ch[300];
            cin>>ch;
            for(int j=0; j<d; j++)
            {
                if(ch[j]=='Y')
                  insert(f + i * 2, f + n * 2 + j + 1, 1);
            }
        }

        for(int i=1; i<=n; i++)
            insert(f + i * 2 - 1, f + 2 * i,1);
            
        int ans=max_flow(2+2*n+f+d,start,end);
        cout<<ans<<endl;
    }
    return 0;
}