HDU 4292 Food（网络流） #by zh

我今天深深的给这题跪了，模版不行是硬伤啊TAT，这题一开始就感觉是网络流题，也想了一个初步的建图策略，不过WA了。后来和人讨论后发现建图有bug，然后果断把人拆点，这样能保证一个人只能计算一次，但是这样建图还是WA了啊……然后我就觉得特别不科学，自己写的数据也能过，但是交上去就是WA，怎么交怎么WA，还以为是自己的建图又错了，但是感觉没有错误的地方。后来实在无奈去问过了的人，发现建图方法和我是一样的有木有！我只能大呼模版坑爹了，不过这不科学啊，之前找这个模版的时候A过一个题，今天的02题也是用这个模版A的，但是05就一直WA，一直WA……模版不行是硬伤啊，无力吐槽了，还好发现了模版的问题也不算太亏，以后就不会再出现这样的问题了，伤不起要不今天能有一个不错的名次的，伤不起，贴一下用了好甜的模版就AC的代码……

#include<stdio.h>
#include<string.h>
const int N=1000;
const int M=400010;
const int inf=1<<30;
struct Edge
{    int v,next,w;    }e[M];
int head[N],cnt;    int n,f,d,ans;    int s,t;    int pre[N],cur[N],dis[N],gap[N];
int q[N],open,tail;
void addedge(int u,int v,int w)
{
    e[cnt].v=v;    e[cnt].w=w;    e[cnt].next=head[u];    head[u]=cnt++;
    e[cnt].v=u;    e[cnt].w=0;    e[cnt].next=head[v];    head[v]=cnt++;
}
void BFS()
{
    int i,u,v;
    memset(gap,0,sizeof(gap));
    memset(dis,-1,sizeof(dis));
    open=tail=0;
    q[open]=t;
    dis[t]=0;
    while (open<=tail)
    {
        u=q[open++];
        for(i=head[u];i!=-1;i=e[i].next)
        {
            v=e[i].v;
            if(e[i].w!=0 || dis[v]!=-1) continue;
            q[++tail]=v;
            ++gap[dis[v]=dis[u]+1];
        }
    }
}
int sap(int n)
{
    int i,v,u,flow=0,aug=inf;
    int flag;
    BFS();
    gap[0]=1;
    for (i=1;i<=n;i++) cur[i]=head[i];
    u=pre[s]=s;
    while (dis[s]<n)
    {   flag=0;
        for (int j=cur[u];j!=-1;j=e[j].next)
     {   v=e[j].v;
            if (e[j].w>0&&dis[u]==dis[v]+1)
           {
                flag=1;    if (e[j].w<aug) aug=e[j].w;     pre[v]=u;      u=v;
                if (u==t)
               {
                    flow+=aug;
                    while (u!=s)
                   {  u=pre[u];   e[cur[u]].w-=aug;    e[cur[u]^1].w+=aug;    }
                    aug=inf;
                }
                break;
            }
            cur[u]=e[j].next;
        }
        if (flag) continue;
        int mindis=n;
        for (int j=head[u];j!=-1;j=e[j].next)
        {
            v=e[j].v;
            if (e[j].w>0&&mindis>dis[v])
           {   mindis=dis[v];    cur[u]=j;    }
        }
        if (--gap[dis[u]]==0) break;
        ++gap[dis[u]=mindis+1];     u=pre[u];
    }
   return flow;
}
int main()
{
    //freopen("input.txt","r",stdin);
    while(scanf("%d%d%d",&n,&f,&d)!=EOF)
    {
        s=f+d+n+n+1; t=s+1;
        cnt = 0;
        memset(head, -1, sizeof(head));
        //memset(vis,0,sizeof(vis));
        //memset(map,0,sizeof(map));
        for(int i=1; i<=f; i++)
        {
            int temp;
            scanf("%d",&temp);
            addedge(s,i,temp);
        }
        for(int i=1; i<=d; i++)
        {
            int temp;
            scanf("%d",&temp);
            addedge(f+i,t,temp);
        }
        for(int i=1; i<=n; i++)
        {
            char c;
            getchar();
            //scanf("%s",map[i]);
            for(int j=0; j<f; j++)
            {
                scanf("%c",&c);
                if(c=='Y')
                {
                    addedge(j+1,f+d+i,1);
                }
            }
        }
        for(int i=1; i<=n; i++)
        {
            char c;
            getchar();
            for(int j=0; j<d; j++)
            {
                scanf("%c",&c);
                if(c=='Y')
                {
                    addedge(f+d+n+i,f+j+1,1);
                }
            }
        }
        for(int i=1; i<=n; i++)
        {
            addedge(f+d+i,f+d+n+i,1);
        }
        n=f+d+2*n+2;
        printf("%d\n", sap(t));
    }
}