hdu 4292 网络流

#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
const int inf=1<<29;
const int maxn=2000,maxm=500000;

#define min(a,b) a>b?b:a

struct edge
{
    int v,f,nxt;
};


bool vis[maxn+10];        //是否被访问
int dist[maxn+10];   //层数

int src,sink;
int g[maxn+10];
int nume;
edge e[maxm*2+10];
void add(int u,int v,int f)
{
    e[++nume].v=v;
    e[nume].f=f;
    e[nume].nxt=g[u];
    g[u]=nume;
    e[++nume].v=u;
    e[nume].f=0;
    e[nume].nxt=g[v];
    g[v]=nume;
}

void bfs()
{
    queue<int >que;
    memset(dist,0,sizeof(dist));
    while(!que.empty()) que.pop();
    que.push(src);
    vis[src]=true;
    while(!que.empty()){
        int u=que.front(),i;
        que.pop();
        for(i=g[u];    i;    i=e[i].nxt){
            if(e[i].f&&!vis[e[i].v]){
                que.push(e[i].v);
                dist[e[i].v]=dist[u]+1;
                vis[e[i].v]=true;
            }
        }
    }
}

int dfs(int u ,int delta)
{
    int i;
    if(u == sink) return delta;
    else{
        int ret=0;
        for(i=g[u];  i&δ  i=e[i].nxt){
            if(e[i].f&&dist[e[i].v]==dist[u]+1){
                int dd=dfs(e[i].v,min(e[i].f,delta));
                e[i].f-=dd;
                e[i^1].f+=dd;
                delta-=dd;
                ret+=dd;
            }
        }
        return ret;
    }
}
int dinic(){
    int ret=0;
    while(true){
        memset(vis,0,sizeof(vis));
        bfs();
        if(!vis[sink]) return ret;
        ret+=dfs(src,inf);
    }
}
int main()
{
    int food[300],drink[300];
    int i,j,n,f,d;
    while(scanf("%d%d%d",&n,&f,&d)!=EOF){
        nume=1;
        memset(g,0,sizeof(g));
        for(i=0;i<f;i++)   scanf("%d",&food[i]);
        for(i=0;i<d;i++)   scanf("%d",&drink[i]);
        getchar();    
        int scr=0;
        sink=f+2*n+d+1;
        char temp;
        for(i=0;i<f;i++) add(scr,i+1,food[i]);
        for(i=0;i<d;i++) add(2*n+f+1+i,f+2*n+d+1,drink[i]);
        for(i=0;i<n;i++) add(f+1+i,f+1+n+i,1);
        for(i=0;i<n;i++){
            for(j=0;j<f;j++){
                scanf("%c",&temp);
                if(temp=='Y')
                    add(j+1,f+1+i,1);
            }
            getchar();
        }
        for(i=0;i<n;i++){
            for(j=0;j<d;j++){
                scanf("%c",&temp);
                if(temp=='Y')
                    add(f+n+1+i,f+2*n+j+1,1);
            }
            getchar();
        }
        printf("%d\n",dinic());
    }
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值