• poj3281(最大流)


    传送门: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
    ",dinic());
        }
    }
    View Code
  • 相关阅读:
    Java线上应用故障排查之一:高CPU占用【转】
    JAVA 之 Tomcat知识框架【转】
    FileZilla 配置备份与还原【转】
    linux删除第几天日志【原创】
    sar命令使用【转】
    Mysql Group Replication 简介及单主模式组复制配置【转】
    grep 中的正则表达式【转】
    java程序out of memory【转】
    Java笔记18:JUnit单元测试
    Java笔记17:导出可执行jar包
  • 原文地址:https://www.cnblogs.com/lienus/p/4293982.html
Copyright © 2020-2023  润新知