POJ 1149

纪念一下我的第一个网络流;最基本的标号法FORD算法,基本上能把图建立出来就是模板了;具体解释见代码:

#include<stdio.h>
#include<algorithm>
#include<queue>
#include<string.h>
#define MAXN 105
#define MAXM 1005
#define INF 100000000
using namespace std;
int last[MAXN],new_flow[MAXN];
int point_c[MAXN][MAXN],point_f[MAXN][MAXN];
int house[MAXM],judge[MAXM];
int n,m,s,t;
int ford()//基本上就是模板,只不过我的书上使用数组我仅仅只是改成了队列;
{
    int res=0;
    memset(point_f,0,sizeof(point_f));
    queue<int>Q;
    while(true)
    {
        while(!Q.empty()) Q.pop();//感觉数组用memset处理应该要快一些,不过应该查不到哪里去;
        memset(last,-1,sizeof(last));
        Q.push(0);
        new_flow[0]=INF;
        while(!Q.empty()&&last[t]==-1)
        {
            int p,v=Q.front();
            Q.pop();
            for(int i=1;i<=t;i++)
            {
                if(last[i]==-1&&(p=point_c[v][i]-point_f[v][i]))
                {
                    last[i]=v;
                    Q.push(i);
                    new_flow[i]=min(new_flow[v],p);
                }
            }
        }
        if(last[t]==-1) break;
        int k1=t,k2=last[k1];
        while(k2!=-1&&k1!=-1)
        {
            point_f[k2][k1]+=new_flow[t];
            point_f[k1][k2]=-point_f[k2][k1];//主要是用于“反悔”,这一步建议与我一样的初学者仔细研究一下,这里困扰了我N久,最后学长给我举了一个例子就秒懂了;在博客末尾我会详细说明一下;
            k1=k2;
            k2=last[k2];
        }
    }
    for(int i=0;i<t;i++)
        res+=point_f[i][t];
    return res;
}
int main()
{
    int num,x,ans;
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        s=0;
        t=n+1;
        memset(judge,0,sizeof(judge));
        memset(point_c,0,sizeof(point_c));
        for(int i=1;i<=m;i++) scanf("%d",&house[i]);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&num);
            for(int j=0;j<num;j++)
            {
                scanf("%d",&x);
                if(judge[x]==0)
                {
                    point_c[s][i]+=house[x];//如果i是某个猪圈的第一个购买者,那么把他与源点连一条边;
                }
                else
                    point_c[judge[x]][i]=INF;//如果不是,那么与上一个购买者连一条边,这里权值应该是无限大
                judge[x]=i;
            }
            scanf("%d",&point_c[i][t]);//购买的数量就是每个购买者与汇点之间那条边的流量上限;
        }
        ans=ford();
        printf("%d\n",ans);
    }
    return 0;
}

首先有一个图:

源点是 0 汇点是3;

边为:

点 0 到 1 权值: 2

点 0 到 2 权值:1

点 1 到 2 权值: 1

点 1 到 3 权值 :1

点 2 到 3 权值:2

自己手动跑一遍就知道了。。。。。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值