Food HDU - 4292 网络流

题意:
有n个人,f种食物和d种饮料,每种饮料和食物的个数已知,假设一个人拿到一个自己喜欢的食物或饮料,那个人就会感到满足。
最多能有几个人满足?

思路:拆点建图。

一个人必须又要饮料又要食物,假设我们直接用食物和饮料把一个人连接起来,是无法保证每个人只会给一组食物饮料的,有可能会出现很多组食物和饮料通过着一个人经过的情况。

所以我们把点转化成边进行限制,让边的流量为1,这样就可以每个人只吃一组食物饮料了。

图如下所建。

在这里插入图片描述

/*Keep on going Never give up*/
#pragma GCC optimize(3,"Ofast","inline")
#include<bits/stdc++.h>
#define endl '\n'
//#define int long long
using namespace std;

const int maxn = 2000;
const int maxm = 1e6+10;

struct edge{
    long long to, next, cap, flow;  // 容量 和 流量
}e[maxm << 1]; int head[maxn], tot;
int n, m, s, t, d[maxn], cur[maxn]; // d[]分层,cur[]弧优化

inline void addedge(int u, int v, int w) {
    e[++tot] = {v, head[u], w, 0}, head[u] = tot;
    e[++tot] = {u, head[v], 0, 0}, head[v] = tot;
}

inline bool bfs() { // 分层 t-s路径优化
    memset(d, 0x3f, sizeof d);
    queue<int> q; q.push(t);
    d[t] = 0;
    while(!q.empty()) {
        int x = q.front(); q.pop();
        for(int i = head[x]; ~i; i = e[i].next) {
            int j = e[i].to;
            if (d[j] == d[0] && e[i^1].cap > e[i^1].flow) {
                d[j] = d[x] + 1;
                if (j == s) return true;
                q.push(j);
            }
        }
    }
    return false;
}

long long dfs(int x, long long flow, long long used = 0) {   //used多路增广
    if (x == t) return flow;
    //for(int i = cur[x]; ~i && flow != used; i = e[ cur[x] = i ].next) {  // 当前弧优化
    for(int i = head[x]; ~i && flow != used; i = e[i].next) {
        long long j = e[i].to, f;
        if (d[x] - 1 == d[j] && e[i].cap > e[i].flow) {
            f = dfs(j, min(e[i].cap-e[i].flow, flow-used));
            if (f) e[i].flow += f, e[i^1].flow -= f, used += f;
        }
    }
    if (!used) d[x] = 0;   // 剪枝优化
    return used;
}

int dinic(){
    long long maxflow = 0;
    while(bfs()) {   // 每次dfs分层都能找到多条增广路
        //memcpy(cur, head, sizeof head);
        maxflow += dfs(s, INT_MAX);
    }
    return maxflow;
}
int food[210],drink[210];
string s1[210],s2[210];
int f,dd;
void init(){
    memset(head,-1,sizeof head);
    tot=-1;
}


int sol(){

    for(int i=1;i<=f;i++) addedge(s,i,food[i]);   //边连吃的
    for(int i=1;i<=n;i++){    //吃的连人
        for(int j=0;j<s1[i].size();j++){
            if(s1[i][j]=='Y'){
                addedge(j+1,i+f,1);
            }
        }
    }
    for(int i=1;i<=n;i++){   //人连人
        addedge(f+i,f+n+i,1);
    }
    for(int i=1;i<=n;i++){
        for(int j=0;j<s2[i].size();j++){
            if(s2[i][j]=='Y'){
                addedge(f+n+i,f+2*n+j+1,1);
            }
        }
    }
    for(int i=1;i<=dd;i++){
        addedge(f+2*n+i,t,drink[i]);
    }
    int res=dinic();
    return res;
}



int main(){

    while(cin>>n>>f>>dd){
        init();
        s=0,t=f+dd+2*n+1;
        for(int i=1;i<=f;i++)  cin>>food[i];
        for(int i=1;i<=dd;i++) cin>>drink[i];
        for(int i=1;i<=n;i++) cin>>s1[i];
        for(int i=1;i<=n;i++) cin>>s2[i];

        cout<<sol()<<endl;
    }


}


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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值