码迷,mamicode.com
首页 > 其他好文 > 详细

poj3281(最大流)

时间:2015-02-16 12:56:37      阅读:132      评论:0      收藏:0      [点我收藏+]

标签:

 

传送门:Dining

题意:一些牛,一些食物,一些饮料,每头牛都有其喜欢的几种食物和几种饮料,求最多能给多少头牛即找到食物又找到饮料~也就是有多少个 牛---食物---饮料 的匹配,而且满足一一匹配,每个牛,食物,或饮料都只能使用一次。

分析:如果直接源点-食物-牛-汇点这样建图,那么可能会很多份食物通过某头牛传送过去,不能保证每头牛吃一份食物,因此拆点,把每头牛拆成牛-牛,边权为1,这样就能保证每头牛只能吃一份了。

技术分享
#pragma comment(linker,"/STACK:1024000000,1024000000")
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <limits.h>
#include <iostream>
#include <algorithm>
#include <queue>
#include <cstdlib>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define LL long long
#define mod 100000000
#define inf 0x3f3f3f3f
#define eps 1e-6
#define N 500
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define PII pair<int,int>
using namespace std;
inline int read()
{
    char ch=getchar();
    int x=0,f=1;
    while(ch>9||ch<0){if(ch==-)f=-1;ch=getchar();}
    while(ch<=9&&ch>=0){x=x*10+ch-0;ch=getchar();}
    return x*f;
}
int n,d,f,vs,vt,tot;
int pre[N],q[N],cur[N],h[N];
struct edge
{
    int v,w,next;
    edge(){}
    edge(int v,int w,int next):v(v),w(w),next(next){}
}e[N*N];
void addedge(int u,int v,int w)
{
    e[tot]=edge(v,w,pre[u]);
    pre[u]=tot++;
}
void init()
{
    memset(pre,-1,sizeof(pre));
    tot=0;
}
/*******************dinic************************/
int bfs()
{
    int head=0,tail=1;
    memset(h,-1,sizeof(h));
    q[0]=vs;h[vs]=0;
    while(head!=tail)
    {
        int u=q[head++];
        for(int i=pre[u];~i;i=e[i].next)
        {
            int v=e[i].v,w=e[i].w;
            if(w&&h[v]==-1)
            {
                h[v]=h[u]+1;
                q[tail++]=v;
            }
        }
    }
    return h[vt]!=-1;
}
int dfs(int u,int flow)
{
    if(u==vt)return flow;
    int used=0;
    for(int i=cur[u];~i;i=e[i].next)
    {
        int v=e[i].v,w=e[i].w;
        if(h[v]==h[u]+1)
        {
            w=dfs(v,min(flow-used,w));
            e[i].w-=w;e[i^1].w+=w;
            if(e[i].w)cur[u]=i;
            used+=w;
            if(used==flow)return flow;
        }
    }
    if(!used)h[u]=-1;
    return used;
}
int dinic()
{
    int res=0;
    while(bfs())
    {
        for(int i=vs;i<=vt;i++)cur[i]=pre[i];
        res+=dfs(vs,inf);
    }
    return res;
}
/********************dinic***********************/
void build()
{
    vs=0;vt=n*2+f+d+1;
    for(int i=1;i<=f;i++)addedge(vs,i,1),addedge(i,vs,0);
    for(int i=1;i<=d;i++)addedge(i+2*n+f,vt,1),addedge(vt,i+2*n+f,0);
    for(int i=1;i<=n;i++)addedge(f+i,f+n+i,1),addedge(f+n+i,f+i,0);
    for(int i=1;i<=n;i++)
    {
        int m1,m2,x;
        m1=read();m2=read();
        while(m1--)
        {
            x=read();
            addedge(x,f+i,1);
            addedge(f+i,x,0);
        }
        while(m2--)
        {
            x=read();
            addedge(f+n+i,f+2*n+x,1);
            addedge(f+2*n+x,f+n+i,0);
        }
    }
}
int main()
{
    while(scanf("%d%d%d",&n,&f,&d)>0)
    {
        init();
        build();
        printf("%d\n",dinic());
    }
}
View Code

 

poj3281(最大流)

标签:

原文地址:http://www.cnblogs.com/lienus/p/4293982.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!