poj 3281 Dining 最大流

题目链接


一开始的想法就是food-cow-water,但是wa了,看了别人的思路知道每一只牛牛用一条边隔开,然后就达到了限流的效果。

所以 食物——>牛牛---->牛牛---->water建完图就成了裸的最大流

代码:

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define inff 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define PI 3.14159265358979323846
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define pii make_pair
#define pr pair<int,int>
const int dir[4][2]={0,1,0,-1,1,0,-1,0};
typedef long long ll;
const ll inFF=9223372036854775807;
typedef unsigned long long ull;
using namespace std;
const int maxm=1e5+5;
const int maxn=1e5+5;
int cur[maxn],d[maxn];
int head[maxn],sign;
int n,f,w;
struct node
{
    int to,p,val;
}edge[maxm];
void init()
{
    sign=0;
    for(int i=0;i<=f+w+2*n+10;i++)
        head[i]=-1;
}
void add(int u,int v,int val)
{
    edge[sign]=node{v,head[u],val};
    head[u]=sign++;
}
int bfs(int st,int ed)
{
    for(int i=0;i<=f+w+2*n+10;i++) d[i]=0;
    d[st]=1;
    queue<int> q;
    q.push(st);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u];~i;i=edge[i].p)
        {
            int v=edge[i].to;
            if(!d[v]&&edge[i].val)
            {
                d[v]=d[u]+1;
                q.push(v);
            }
        }
    }
    if(d[ed]==0) return false;
    return true;
}
int dfs(int u,int ed,int low)
{
    if(u==ed) return low;
    int temp=low,a;
    for(int &i=cur[u];~i;i=edge[i].p)
    {
        int v=edge[i].to;
        if(d[v]==d[u]+1&&edge[i].val)
        {
            if(a=dfs(v,ed,min(low,edge[i].val)))
            {
                edge[i].val-=a;
                edge[i^1].val+=a;
                temp-=a;
                if(temp==0) break;
            }
        }
    }
    return low-temp;
}
int main()
{
    int num1,num2,x;
    while(cin>>n>>f>>w)
    {
        init();
        int s=f+w;
        for(int p=s+1;p<=s+n;p++)
        {
            scanf("%d %d",&num1,&num2);
            while(num1--)
            {
                scanf("%d",&x);
                add(x,p,1),add(p,x,0);
            }
            while(num2--)
            {
                scanf("%d",&x);
                add(p+n,x+f,1),add(x+f,p+n,0);
            }
        }
        for(int i=s+1;i<=s+n;i++) add(i,i+n,1),add(i+n,i,0);
        for(int i=1;i<=f;i++) add(0,i,1),add(i,0,0);
        for(int i=f+1;i<=s;i++) add(i,s+2*n+1,1),add(s+2*n+1,i,0);
        int ans=0;
        while(bfs(0,s+2*n+1))
        {
            for(int i=0;i<=s+2*n+1;i++)
                cur[i]=head[i];
            ans+=dfs(0,s+2*n+1,inff);
        }
        printf("%d\n",ans);
    }
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值