最大流:HDU-4292(Food)

最新推荐文章于 2021-05-10 20:27:32 发布

Teisei

最新推荐文章于 2021-05-10 20:27:32 发布

阅读量603

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/Teisei11/article/details/21610683

版权

5 篇文章 0 订阅

订阅专栏

#include <iostream>

#include <cstdio>

#include <cstring>

#define INF 0x7fffffff

#define MAXN 1000

#define MAXE 200000

using namespace std;

int N,F,D,Max;

struct Edge {
    int to,next,c;
}e[MAXE];
int size,head[MAXN];

void addedge(int from,int to,int c)

{
   e[size].to=to;
    e[size].c = c;
    e[size].next=head[from];
    head[from]=size++;

    e[size].to=from;
    e[size].c=0;
    e[size].next=head[to];
    head[to]=size++;

int dis[MAXN],gap[MAXN],cur[MAXN],pre[MAXN];

int SAP_GAP()

{
     int start=0,end=Max+1,NN=Max+2;

       int cur_flow,flow_ans=0,i,u,tmp,neck;

       memset(dis,0,sizeof(dis));

        memset(gap,0,sizeof(gap));

        memset(pre,-1,sizeof(pre));

      for(i=0;i<=Max+1;i++)
        cur[i]=head[i];

       u=start;
    gap[0]=NN;

      while(dis[start]<NN)

      {
   if(u==end)

           {
   cur_flow=INF;

                   for(i=start;i!=end;i=e[cur[i]].to)

                      {
   if(cur_flow>e[cur[i]].c)

                             {
                    neck=i;
                    cur_flow=e[cur[i]].c;
                }

                    for(i=start;i!=end;i=e[cur[i]].to)

                                 tmp=cur[i];

                                 e[tmp].c-=cur_flow;

                                      e[tmp^1].c+=cur_flow;

                           u=neck;
            flow_ans+=cur_flow;

           for(i=cur[u];i!=-1;i=e[i].next)

                         if(e[i].c&&dis[u]==dis[e[i].to]+1)
                break;

           if(i!=-1)
        {
            cur[u]=i;
            pre[e[i].to]=u;
            u=e[i].to;
        }

           else
        {

                    if(--gap[dis[u]]==0)    break;

                       cur[u]=head[u];

                    for(tmp=NN,i=head[u];i!=-1;i=e[i].next)

                                 if(e[i].c)  tmp=min(tmp,dis[e[i].to]);

                            dis[u]=tmp+1;

                           gap[dis[u]]++;

                    if(u!=start)    u=pre[u];

  return flow_ans;

void init()

{
   char str[MAXN];
    int k;
    size=0;

    memset(head,-1,sizeof(head));

     Max=N*2+F+D;

     for(int i=1;i<=F;i++)
    {
        scanf("%d",&k);
        addedge(0,i,k);
    }

     for(int i=1;i<=D;i++)
    {
        scanf("%d",&k);
        addedge(N*2+F+i,Max+1,k);
    }

     for(int i=1;i<=N;i++)
        addedge(F+i,F+N+i,1);

     for(int i=1;i<=N;i++)

        {
     scanf("%s",str);

              for(int j=1;j<=F;j++)
            if(str[j-1]=='Y')
                addedge(j,F+i,1);
    }

             for(int i=1;i<=N;i++)

                {   scanf("%s",str);

                  for(int j=1;j<=D;j++)

                              if(str[j-1]=='Y')
                addedge(F+N+i,F+N*2+j,1);

int main()

{
    //freopen("in.txt","r",stdin);

 while(scanf("%d%d%d",&N,&F,&D)!=EOF)

         {
        init();
        printf("%d\n",SAP_GAP());
    }

     return 0;

关注

专栏目录