POJ 3281 Dining (最大流)

思路:

类似与hihocoder的二分图多重匹配。
虚拟源点,源点连食物,食物连牛,牛连对应着的的自己,对应着的牛连水,水连汇点。
边权全部为1.
拆点是为了保证每头牛只通过一条流,即每头牛只吃一份套餐~
我是这么编号的:
源点:0
食物:1~f
牛:1+f ~n+f
新牛:1+n+f~n+n+f
水:1+n*2+f~d+n*2+f
汇点:2*n+f+d+1

ps:

接下来就是套板子了,我这个是dinic+当前弧优化。跑的比xx记者都快,不慌。(实在不行就换isap就好了啊2333。数组大小我乱开的,读者不要在意这些细节。。)

#include <iostream>
#include <cstdio>
#include <queue>
#include <string.h>
typedef long long int lli;
using namespace std;

const int maxn = 5200;
struct edge{
    lli to,next,v;
}ed[600000+20000];
lli ans;
lli head[maxn],dis[maxn],cnte,cur[maxn];
void ini(){
    memset(head,-1,sizeof(head));
    ans = 0;cnte = 0;
}
inline void ae(int x,int y,int v){
    ed[cnte].to = y;
    ed[cnte].v = v;
    ed[cnte].next=head[x];
    head[x]=cnte++;

    ed[cnte].to = x;
    ed[cnte].v = 0;
    ed[cnte].next=head[y];
    head[y]=cnte++;
}
int bfs(int s,int t){
    memset(dis,-1,sizeof(dis));
    queue<int> q;q.push(s);
    dis[s] = 0;
    int now;
    while(!q.empty()){
        now = q.front();q.pop();
        for(int i = head[now];i!=-1;i=ed[i].next){
            int v = ed[i].to;
            if(dis[v] == -1 && ed[i].v > 0){
                dis[v] = dis[now]+1;q.push(v);
                if(v == t) return 1;
            }
        }
    }
    if(dis[t] < 0) return 0;
    return 1;
}
lli dfs(int s,int t,lli low){
    if(s==t || !low) return low;
    lli a,ret = 0;
    for(lli &i = cur[s];i!=-1;i=ed[i].next){
        int v= ed[i].to;
        if(dis[v]==dis[s]+1 && ed[i].v>0 && (a=dfs(v,t,min(low,ed[i].v)))){
            ed[i].v -= a; ed[i^1].v += a;
            low -= a; ret += a;
            if(!low) break;
        }
    }
    return ret;
}

void dinic(int s,int t){
    while(bfs(s,t)){
        memcpy(cur,head,sizeof(head));
        ans += dfs(s,t,0x3f3f3f3f);
    }
}

int main(){
    int n,a,b,v,f,d,fi,di;
    ini();
    scanf("%d%d%d",&n,&f,&d);
    int s = 0,t = 1+n*2+f+d;
    for(int i = 1;i <= f;i++){//源点向食物建边
        ae(s,i,1);
    }
    for(int i = 1;i <= d;i++){//水向汇点建边
        ae(i+f+n*2,t,1);
    }
    for(int i = 1;i <= n;i++){//牛向牛建边 以此来实现对流量的控制
        ae(i+f,i+f+n,1);
    }
    for(int i = 1;i <= n;i++){
        scanf("%d%d",&fi,&di);
        for(int j = 1;j <= fi;j++){
            scanf("%d",&v);
            ae(v,f+i,1);
        }
        for(int j = 1;j <= di;j++){
            scanf("%d",&v);
            ae(i+f+n,v+f+n*2,1);
        }
    }
    dinic(s,t);
    printf("%lld\n",ans);

    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值